机械优化设计程序(精选).docVIP

应用变尺度法计算 计算 f(x1,x2)=x1^2+2*x2^2-4*x1-2*x1*x2 的无约束极值，初始点x0=[1,1]。 /* tt ---- 一维搜索初始步长 ff ---- 差分法求梯度时的步长 ac ---- 终止迭代收敛精度 ad ---- 一维搜索收敛精度 n ----- 设计变量的维数 xk[n] -- 迭代初始点 */ #includestdio.h #includestdlib.h #includemath.h #includeconio.h #define tt 0.01 #define ff 1.0e-6 #define ac 1.0e-6 #define ad 1.0e-6 #define n 2 double ia; double fny(double *x) { double x1=x[0],x2=x[1]; double f; f=x1*x1+2*x2*x2-4*x1-2*x1*x2; return f; } double * iterate(double *x,double a,double *s) { double *x1; int i; x1=(double *)malloc(n*sizeof(double)); for(i=0;in;i++) x1[i]=x[i]+a*s[i]; return x1; } double func(double *x,double a,double *s) { double *x1; double f; x1=iterate(x,a,s); f=fny(x1); return f; } void finding(double a[3],double f[3],double *xk,double *s) { double t=tt; int i; double a1,f1; a[0]=0;f[0]=func(xk,a[0],s); for(i=0;;i++) { a[1]=a[0]+t; f[1]=func(xk,a[1],s); if(f[1]f[0]) break; if(fabs(f[1]-f[0])=ad) { t=-t; a[0]=a[1];f[0]=f[1]; } else { if(ia==1) return; //break t=t/2;ia=1; } } for(i=0;;i++) { a[2]=a[1]+t; f[2]=func(xk,a[2],s); if(f[2]f[1]) break; t=2*t; a[0]=a[1];f[0]=f[1]; a[1]=a[2];f[1]=f[2]; } if(a[0]a[2]) { a1=a[0]; f1=f[0]; a[0]=a[2]; f[0]=f[2]; a[2]=a1; f[2]=f1; } return; } double lagrange(double *xk,double *ft,double *s) { int i; double a[3],f[3]; double b,c,d,aa; finding(a,f,xk,s); for(i=0;;i++) { if(ia==1) { aa=a[1]; *ft=f[1]; break; } d=(pow(a[0],2)-pow(a[2],2))*(a[0]-a[1])-(pow(a[0],2)-pow(a[1],2))*(a[0]-a[2]); if(fabs(d)==0) break; c=((f[0]-f[2])*(a[0]-a[1])-(f[0]-f[1])*(a[0]-a[2]))/d; if(fabs(c)==0) break; b=((f[0]-f[1])-c*(pow(a[0],2)-pow(a[1],2)))/(a[0]-a[1]); aa=-b/(2*c); *ft=func(xk,aa,s); if(fabs(aa-a[1])=ad) {if(*ftf[1]) aa=a[1];break;} if(aaa[1]) { if(*ftf[1]) {a[2]=aa;f[2]=*ft;} else if(*ftf[1]) {a[0]=a[1];a[1]=aa;f[0]=f[1];f[1]=*ft;} else if(*ft==f[1]) { a[2]=aa;a[0]=a[1]; f[2]=*ft;f[0]=f[1];