1、精選優質文檔-傾情為你奉上/無約束最優化坐標輪換法/C語言#include <stdio.h> #include <math.h> #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,flo

2、at 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<n;i+) x1i=x0i; x2i=x0i+h*si; f1=f(x1); f2=f(x2); if(f2>=f1) h=(-1)*h; for(i=

3、0;i<n;i+) x3i=x1i; for(i=0;i<n;i+) x1i=x2i; for(i=0;i<n;i+) x2i=x3i;for(i=0;i<n;i+) x3i=x1i+h*si; f3=f1; f1=f2;f2=f3; f3=f(x3); h=2*h;for(i=0;i<n;i+) x3i=x2i+h*si; f3=f(x3); while(f3<f2) for(i=0;i<n;i+) x1i=x2i; f1=f2; for(i=0;i<n;i+) x2i=x3i; f2=f3; for(i=0;i<n;i+) x3i=x2

4、i+h*si; f3=f(x3); for(i=0;i<n;i+) if(h>0) 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;i<n;i+) x1i=ai+(float)0.382*(bi-ai); f1=f(x1); for(i=0;i<n;i+) x2i=ai+(float)0.618*(bi-ai); f2=f(

5、x2); do if(f1<f2) for(i=0;i<n;i+) bi=x2i; for(i=0;i<n;i+) x2i=x1i; f2=f1; for(i=0;i<n;i+) x1i=ai+(float)0.382*(bi-ai); f1=f(x1); else for(i=0;i<n;i+) ai=x1i; for(i=0;i<n;i+) x1i=x2i; f1=f2; for(i=0;i<n;i+) x2i=ai+(float)0.618*(bi-ai); f2=f(x2); sum=0; for(i=0;i<n;i+) sum+=(bi

6、-ai)*(bi-ai); while(sqrt(sum)>flag); for(i=0;i<n;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;i<n

7、;i+) for(k=0;k<n;k+) if(k=0)sik=c; else sik=d; k=1; do for(i=0;i<n;i+) x1i=x0i; for(i=0;i<n;i+) mjtf(n,x1, h,si,a,b); mhjfgf(n, a, b, flag1, x1); sum=0; for(i=0;i<n;i+) sum= sum+(x1i-x0i)*(x1i-x0i); if(sqrt(sum) <= flag2) break; else k=k+1; for(i=0;i<n;i+) x0i=x1i; while(1); for(i=

8、0;i<n;i+) xi=x1i; /*坐標輪換法主程序*/ void main() int i,n; float h,flag1,flag2,x0m,am,bm,xm; printf("n<坐標輪換法>n"); printf("請輸入維數:n"); scanf("%d",&n); printf("請輸入初始點:"); for(i=0;i<n;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,flag