机械优化设计坐标轮换发c语言.docVIP

#include stdio.h #include math.h #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 zblhf(int n,float x0[],float h,float flag1,float flag2,float a[],float b[],floatx[]); /*目标函数(n维)*/ /*入口参数： x ：n维数组，自变量 */ /*返回值 ：函数值 */ float f(float x[]) { float result; result=60-10*(x[0])-4*(x[1])+(x[0]*x[0])+(x[1]*x[1])-(x[0]*x[1]); return result; } /*多维进退法子程序*/ /*入口参数： n ：优化模型维数 x0 ：n维数组，初始点坐标 h ：初始搜索步长 s ：n维数组，搜索方向 */ /*出口参数： a ：n维数组，搜索区间下限 b ：n维数组，搜索区间上限*/ void mjtf(int n,float x0[],float h,float s[],float a[],float b[]) { int i; float x1[m],x2[m],x3[m],f1,f2,f3; for(i=0;in;i++) /*计算初始两试点*/ { x1[i]=x0[i]; x2[i]=x0[i]+h*s[i]; } f1=f(x1); f2=f(x2); if(f2=f1) /*判断搜索方向*/ { /*搜索方向为反向，转身*/ h=(-1)*h; for(i=0;in;i++) x3[i]=x1[i]; f3=f1; for(i=0;in;i++) x1[i]=x2[i]; f1=f2; for(i=0;in;i++) x2[i]=x3[i]; f2=f3; } /*搜索方向为正向*/ for(i=0;in;i++) /*计算第三试点*/ x3[i]=x2[i]+h*s[i]; f3=f(x3); while(f3f2) /*判断是否未完成搜索*/ { /*未完成，继续搜索*/ h=2*h; for(i=0;in;i++) x1[i]=x2[i]; f1=f2; for(i=0;in;i++) x2[i]=x3[i]; f2=f3; for(i=0;in;i++) x3[i]=x2[i]+h*s[i]; f3=f(x3); } /*已完成*/ for(i=0;in;i++) /*输出初始搜索区间*/ { if(h0) //if(x1[i]x3[i]) { a[i]=x1[i]; b[i]=x3[i]; } else { a[i]=x3[i]; b[i]=x1[i]; } } } /*多维黄金分割法子程序*/ void mhjfgf(int n,float a[],float b[],float flag,float x[]) { int i; float x1[m],x2[m],f1,f2,sum; for(i=0;in;i++) /*计算初始两试点*/ x1[i]=b[i]-(float)0.618*(b[i]-a[i]); f1=f(x1); for(i=0;in;i++) x2[i]=a[i]+(float)0.618*(b[i]-a[i]); f2=f(x2); do { if(f1=f2) /*判断消去区间*/ { /*消去右*/ for(i=0;in;i++) b[i]=x2[i]; for(i=0;in;i++) x2[i]=x1[i]; f2=f1; for(i=0;in;i++) x1[i]=b[i]-(float)0.618*(b[i]-a[i]); f1=f(x1); } else { /*消去左*/ for(i=0;in;i++) a[i]=x1[i]; for(i=0;in;i++) x1[i]=x2[i]; f1=f2; for(i=0;in;i++) x2[i]=a[i]+(float)0.618*(b[i]-a[i]); f2=f(x2); } sum=0; for(i=