数学与非线性科学 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

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

求助：三等分法程序

最近在写最优化论文遇到这样一个问题，本人想用三等分代替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）[原创]求救，谁能帮我。

三等分法程序吗？

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

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

页: [1]

声振论坛's Archiver

求助：三等分法程序

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