鲍威尔法编程-powell法编程-c语言编程-c++6.0.doc

《鲍威尔法编程-powell法编程-c语言编程-c++6.0.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、#include#defineN2floatgs(floatz[N]){floatf;//f=10*(z[0]+z[1]-5)*(z[0]+z[1]-5)+(z[0]-z[1])*(z[0]-z[1]);f=4+4.5*z[0]-4*z[1]+z[0]*z[0]+2*z[1]*z[1]-2*z[0]*z[1]+z[0]*z[0]*z[0]*z[0]-2*z[0]*z[0]*z[1];//f=1.5*z[0]*z[0]+0.5*z[1]*z[1]-z[0]*z[1]-2*z[0];return(f);}floatywyh(f

2、loatx[N],floatt,floats[N],floatz[N]){floatq=0.618,e,A[N],f[2],a,b,c=0.1,d=0.1;inti,j=0;a=0;b=t;f[0]=gs(x);for(i=0;if[1])a=b-t;elsebreak;}while(f[0]

3、>f[1]);}else{t=0-t;do{a=a+t;f[1]=f[0];for(i=0;if[0]){b=a-t;t=2*t;}elsebreak;}while(f[1]>f[0]);}do{j=j+1;A[0]=b-q*(b-a);for(i=0;i

4、{if(f[0]>f[1]){a=A[0];A[0]=A[1];f[0]=f[1];A[1]=a+q*(b-a);for(i=0;i50)break;}else{b=A[1];A[1]=A[0];f[1]=f[0];A[0]=b-q*(b-a);for(i=0;i50)break;}if((f[1]-f[0])/f[1]>0)e=(f[1]-f[0

5、])/f[1];elsee=0-(f[1]-f[0])/f[1];}while(e>c);if((A[1]-A[0])/A[0]>0)e=(A[1]-A[0])/A[0];elsee=0-(A[1]-A[0])/A[0];if(e<=d)break;else{a=A[0];b=A[1];}if((A[1]-A[0])/A[0]>0)e=(A[1]-A[0])/A[0];elsee=0-(A[1]-A[0])/A[0];}while(e>d);if(f[0]

6、return(f[0]);}else{for(i=0;i

7、,Ex=");scanf("%f,%f",&c,&d);printf("请输入步长T=");scanf("%f",&t);k=0;ks=0;for(i=0;i

8、,s,z);for(j=0;j