无约束坐标轮换法编程.doc

。/无约束最优化坐标轮换法/C语言#include #include #define m 10 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 zblhf(int n,float x0,float h,float flag1,float flag2,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; for(i=0;in;i+) x1i=x2i; for(i=0;in;i+) x2i=x3i;for(i=0;in;i+) x3i=x1i+h*si; f3=f1; f1=f2;f2=f3; f3=f(x3); h=2*h;for(i=0;in;i+) x3i=x2i+h*si; f3=f(x3); while(f3f2) 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;i0) 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,f0; for(i=0;in;i+) x1i=ai+(float)0.382*(bi-ai); f1=f(x1); for(i=0;in;i+) x2i=ai+(float)0.618*(bi-ai); f2=f(x2); do if(f1f2) for(i=0;in;i+) bi=x2i; for(i=0;in;i+) x2i=x1i; f2=f1; for(i=0;in;i+) x1i=ai+(float)0.382*(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); for(i=0;in;i+) xi=(float)0.5*(bi+ai); f0=f(x); /*坐标轮换法子程序*/ void zblhf(int n,float x0,float h,float flag1,float flag2,float a,float b,float x) int i,k; float x1m,smm,sum,c,d; printf(输入搜索方向：n);scanf(%f%f,&c,&d);for(i=0;in;i+) for(k=0;kn;k+) if(k=0)sik=c; else sik=d; k=1; do for(i=0;in;i+) x1i=x0i; for(i=0;in;i+) mjtf(n,x1, h,si,a,b); mhjfgf(n, a, b, flag1, x1); sum=0; for(i=0;in;i+) sum= sum+(x1i-x0i)*(x1i-x0i); if(sqrt(sum) = flag2) break; else k=k+1; for(i=0;in;i+) x0i=x1i; while(1); for(i=0;in;i+) xi=x1i; /*坐标轮换法主程序*/ void main() int i,n; float h,flag1,flag2,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请输入黄金分割法迭代精度:n); scanf(%f,&flag1); printf(n请输入坐标轮换法迭代精度:n); scanf(%f,&flag2); zblhf(n,x0,h,flag1,flag2,a,b,x); printf(n极小点坐标为:n); for(i=