优化设计-鲍威尔法程序(c语言).doc

#include #include #define m 10 /*数组长度m = 维数n */float f(float x);void mjtf(int n,float x0,float h,float s,float a,float b);void mhjfgf(int n,float a,float b,float flag,float x);void mbwef(int n,float x0,float h,float flag,float a,float b,float x);float f(float x) float result; result=60-10*x0-4*x1+x0*x0+x1*x1-x0*x1; return result;/*多维进退法子程序*/void mjtf(int n,float x0,float h,float s,float a,float b) int i; float x1m,x2m,x3m,f1,f2,f3; for(i=0;i=f1) /*判断搜索方向*/ /*搜索方向为反向，转身*/ h=(-1)*h; for(i=0;in;i+) x3i=x1i; f3=f1; for(i=0;in;i+) x1i=x2i; f1=f2; for(i=0;in;i+) x2i=x3i; f2=f3; /*搜索方向为正向*/ for(i=0;in;i+) /*计算第三试点*/ x3i=x2i+h*si; f3=f(x3); while(f3f2) /*判断是否未完成搜索*/ /*未完成，继续搜索*/ h=2*h; for(i=0;in;i+) x1i=x2i; f1=f2; for(i=0;in;i+) x2i=x3i; f2=f3; for(i=0;in;i+) x3i=x2i+h*si; f3=f(x3); /*已完成*/ for(i=0;in;i+) /*输出初始搜索区间*/ if(x1ix3i) ai=x1i; bi=x3i; else ai=x3i; bi=x1i; /*多维黄金分割法子程序*/void mhjfgf(int n,float a,float b,float flag,float x) int i; float x1m,x2m,f1,f2,sum; for(i=0;in;i+) /*计算初始两试点*/ x1i=bi-(float)0.618*(bi-ai); f1=f(x1); for(i=0;in;i+) x2i=ai+(float)0.618*(bi-ai); f2=f(x2); do if(f1=f2) /*判断消去区间*/ /*消去右*/ for(i=0;in;i+) bi=x2i; for(i=0;in;i+) x2i=x1i; f2=f1; for(i=0;in;i+) x1i=bi-(float)0.618*(bi-ai); f1=f(x1); else /*消去左*/ for(i=0;in;i+) ai=x1i; for(i=0;in;i+) x1i=x2i; f1=f2; for(i=0;in;i+) x2i=ai+(float)0.618*(bi-ai); f2=f(x2); sum=0; for(i=0;iflag*0.1); for(i=0;in;i+) xi=(float)0.5*(bi+ai); /*鲍威尔法子程序*/void mbwef(int n,float x0,float h,float flag,float a,float b,float x) int i,j,k,r; float x1m,x2m,f0,f1,f2,fnm,smm,sum; for(i=0;in;i+) for(k=0;kn;k+) if(i=k) sik=1; else sik=0; k=1; while(1) for(i=0;in;i+) x1i=x0i; for(i=0;in;i+) mjtf(n,x1,h,si,a,b); mhjfgf(n,a,b,flag,x1); fni=f(x0)-f(x1); for(i=0;in;i+) x2i=2*x1i-x0i; for(i=1;in;i+) if(fn0=f0|(f0-2*f1+f2)*(f0-f1-fn0)*(f0-f1-fn0)=0.5*fn0*(f0-f2)*(f0-f2) sum=0; for(i=0;in;i+) sum+=(x1i-x0i)*(x1i-x0i); if(f1=f2) for(i=1;in;i+) x0i=x1i; else for(i=1;in;i+) x0i=x2i; else for(i=r;in;i+) for(j=0;jn;j+) sij=si+1j; for(i=0;in;i+) sni=x1i-x0i; mjtf(n,x1,h,sn,a,b); mhjfgf(n,a,b,flag,x1); sum=0; for(i=0;in;i+) sum+=(x1i-x0i)*(x1i-x0i); for(i=0;in;i+) x0i=x1i; if(sqrt(sum)=flag) break; else k+=1; for(i=0;in;i+) xi=x1i;/*鲍威尔法主程序*/void main() int i,n; float h,flag,x0m,am,bm,xm; printf(nn); printf(请输入维数:n); scanf(%d,&n); printf(请输入初始点:); for(i=0;in;i+) printf(nx0%d=,i); scanf(%f,&x0i); printf(n请输入初始步长:n); scanf(%f,&h); printf(n请输入精度: