应用变尺度法计算计算f(x1,x2)=x1^2+2*x2^2-4*x1-2*x1*x2的无约束极值，初始点x0=[1,1]。/*tt----一维搜索初始步长ff----差分法求梯度时的步长ac----终止迭代收敛精度ad----一维搜索收敛精度n-----设计变量的维数xk[n]--迭代初始点*/#include<stdio.h>#include<stdlib.h>#include<math.h>#include<conio.h>#definett0.01#defineff1.0e-6#defineac1.0e-6#definead1.0e-6#definen2doubleia;doublefny(double*x){doublex1=x[0],x2=x[1];doublef;f=x1*x1+2*x2*x2-4*x1-2*x1*x2;returnf;}double*iterate(double*x,doublea,double*s){double*x1;inti;x1=(double*)malloc(n*sizeof(double));for(i=0;i<n;i++)x1[i]=x[i]+a*s[i];returnx1;}doublefunc(double*x,doublea,double*s){double*x1;doublef;x1=iterate(x,a,s);f=fny(x1);returnf;}voidfinding(doublea[3],doublef[3],double*xk,double*s){doublet=tt;inti;doublea1,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;//breakt=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;}doublelagrange(double*xk,double*ft,double*s){inti;doublea[3],f[3];doubleb,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(*ft>f[1])aa=a[1];break;}if(aa>a[1]){if(*ft>f[1]){a[2]=aa;f[2]=*ft;}elseif(*ft<f[1]){a[0]=a[1];a[1]=aa;f[0]=f[1];f[1]=*ft;}elseif(*ft==f[1]){a[2]=aa;a[0]=a[1];f[2]=*ft;f[0]=f[1];a[1]=(a[0]+a[2])/2;f[1]=func(xk,a[1],s);}}else{if(*ft>f[1]){a[0]=aa;f[0]=*ft;}elseif(*ft<f[1]){a[2]=a[1];a[1]=aa;f[2]=f[1];f[1]=*ft;}elseif(*ft==f[1]){a[0]=aa;a[2]=a[1];f[0]=*ft;f[2]=f[1];a[1]=(a[0]+