求助：三等分法程序

x0416 · 发表于 2006-6-10 00:53

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

最近在写最优化论文遇到这样一个问题，本人想用三等分代替0.618法去解方程的最值点，然后作比较（虽然最终肯定还是0.618法要好点）但三等分法程序无法写出，那位大侠有想过这样的问题，有没有写过这方面　的程序，能否发给小弟，<a href="mailtx0416@163.com" target="_blank" >x0416@163.com</A>.

[ 本帖最后由 xinyuxf 于 2006-12-26 11:50 编辑 ]

tammy · 发表于 2006-6-11 07:06

三等分法程序吗？

x0416 · 发表于 2006-6-13 15:19

是呀，我自已写了一个，但有一个错误不知道怎么改，各位大侠请赐教。 #include<math.h> #include<stdio.h> //全局变量 double a[100],b[100]; double L,X; double R[100],U[100]; int k=1; //要计算的函数 double function(double x){ double y=(x+1)*(x-2)*(x-2); return y; } //初始化区间 void Initalize(){ printf("Input the a[1] and b[1]:\n"); scanf("%lf %lf",&a[1],&b[1]); if(a[1]>b[1]){ printf("error,the a[1] must be small than the b[1],please input again.\n"); exit(0); } printf("Input the L\n"); scanf("%lf",&L); if(L<0){ printf("error,the L must be biger than zero.\n"); exit(0); } } //step 3 void Step3(){ a[k+1]=R[k]; b[k+1]=b[k]; R[k+1]=U[k]; U[k+1]=a[k+1]+2/3*(b[k+1]-a[k+1]); } //step4 void Step4(){ a[k+1]=a[k]; b[k+1]=U[k]; U[k+1]=R[k]; R[k+1]=a[k+1]+1/3*(b[k+1]-a[k+1]); } //处理过程 void process(){ R[1]=a[1]+1/3*(b[1]-a[1]); U[1]=a[1]+2/3*(b[1]-a[1]); while((b[k]-a[k])>=L){ if(function(R[k])>function(U[k])){ Step3(); } else { Step4(); } k=k+1; } printf("the min area is :[%lf %lf]\n",a[k],b[k]); double X=(a[k]+b[k])/2; printf("the min point is:%lf\n",X); } //主函数 void main(){ Initalize(); process(); } getch();

账号		自动登录	找回密码
密码			我要加入

[应用数学] 求助：三等分法程序

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（x0416）[原创]求救，谁能帮我。