鲍威尔法C++源程序.doc

《鲍威尔法C++源程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、#include#include#include#include#include#include#definen12#definett0.005#definead0.//定义常量//tt　　初始迭代步长//ad　　收敛精度floatia;floatfny(float*x){floatf;f=10*pow((x[0]+x[1]-5),2)+pow((x[0]-x[1]),2);//目标函数return(f);}float*iterate(float*x,float

2、a,float*s){float*x1;x1=(float*)malloc(n1*sizeof(float));for(inti=0;i

3、;for(inti=0;;i++){a[1]=a[0]+t;f[1]=func(xk,a[1],s);if(f[1]=ad){t=-t;a[0]=a[1];f[0]=f[1];}else{if(ia==1)return;t=t/2;ia=1;}}for(i=0;;i++){a[2]=a[1]+t;f[2]=func(xk,a[2],s);if(f[2]>f[1])break;t=2*t;a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];}if(a[0]>a[2]){a1=a[0];f1

4、=f[0];a[0]=a[2];f[0]=f[2];a[2]=a1;f[2]=f1;}return;}//secondinsertfloatlagrange(float*xk,float*ft,float*s){floata[3],f[3];floatb,c,d,aa;finding(a,f,xk,s);for(inti=0;;i++){if(ia==1){aa=a[1];*ft=f[1];break;}d=(pow(a[0],2)-pow(a[2],2))*(a[0]-a[1])-(pow(a[0],2)-pow(a[1],2))*(a[0]-a[2]);if(fabs(d

5、)==0)break;c=((f[0]-f[2])*(a[0]-a[1])-(f[0]-f[1])*(a[0]-a[2]))/d;if(fabs(c)==0)break;b=((f[0]-f[1])-c*(pow(a[0],2)-pow(a[1],2)))/(a[0]-a[1]);aa=-b/(2*c);*ft=func(xk,aa,s);if(fabs(aa-a[1])<=ad){if(*ft>f[1])aa=a[1];break;}if(aa>a[1]){if(*ft>f[1]){a[2]=aa;f[2]=*ft;}elseif(*ft

6、[1]=aa;f[0]=f[1];f[1]=*ft;}elseif(*ft==f[1]){a[2]=aa;a[0]=a[1];f[2]=*ft;f[0]=f[1];a[1]=(a[0]+a[2])/2;f[1]=func(xk,a[1],s);}}else{if(*ft>f[1]){a[0]=aa;f[0]=*ft;}elseif(*ft

7、1]=func(xk,a[1],s);}}}if(*ft>f[1]){aa=a[1];*ft=f[1];}return(aa);}float*powell(float*xk){floath[n1][n1],s[n1]={0,0},ff[n1+1]={0,0,0};floatf1,f3,aa;floatdk[n1],*x0,xk1[n1];intm=0,i,j;for(i=0;i