约束坐标轮换法.docxVIP

WORD格式 专业资料整理 #includestdio.h double minfun(double x1,double x2) {double f; f=x1*x1+2*x2*x2-4*x1-8*x2+15; return f; } double g1(double x1,double x2) {double f; f=9-x1*x1-x2*x2; return f; } double g2(double x1,double x2) {double f; f=x1+0*x2; return f; } double g3(double x1,double x2) {double f; f=0*x1+x2; return f; } int main() { double x0[2],h0,e,x1[2],x2[2],a0; double e1[2]={1,0},e2[2]={0,1}; double fx0,fx1,fx2,fx11[100],fx22[100]; int i,j,k; h0=0.5;e=0.0001; printf( 请输入初始点 x0： ); scanf(%lf%lf,x0[0],x0[1]); printf(%lf %lf\n,x0[0],x0[1]); if(g1(x0[0],x0[1])=0 g2(x0[0],x0[1])=0 g3(x0[0],x0[1])=0) {fx0=minfun(x0[0],x0[1]); printf(fx0=%lf\n,fx0); } else {printf( 请重新输入初始点 x0： ); scanf(%lf%lf,x0[0],x0[1]); } // 实现初始点的判断 ; for(k=0;;k++) { for(i=0;i2;i++) x1[i]=x0[i]+h0*e1[i]; printf(%lf %lf\n,x1[0],x1[1]); if(g1(x1[0],x1[1])=0 g2(x1[0],x1[1])=0 g3(x1[0],x1[1])=0) {fx1=minfun(x1[0],x1[1]); printf(fx1=%lf\n,fx1); if(fx1fx0) a0=h0; else a0=-h0; } else a0=-h0;// 首次加大步长 ; for(j=0;;) {for(i=0;i2;i++) x1[i]=x0[i]+a0*e1[i]; printf(x1[0]=%f x1[1]=%f\n,x1[0],x1[1]); if(g1(x1[0],x1[1])=0 g2(x1[0],x1[1])=0 g3(x1[0],x1[1])=0) {fx11[j]=minfun(x1[0],x1[1]); printf(fx1=%lf\n,fx11[j]); if(fx11[j]fx0) {fx0=fx11[j]; x0[0]=x1[0]; x0[1]=x1[1]; a0=2*a0; j++; }// 实现 X 方向的多次搜索； else break; } else break; } for(i=0;i2;i++) x2[i]=x0[i]+a0*e2[i]; printf(x2[0]=%f x2[1]=%f\n,x2[0],x2[1]); if(g1(x2[0],x2[1])=0 g2(x2[0],x2[1])=0 g3(x2[0],x2[1])=0) {fx2=minfun(x2[0],x2[1]); printf(fx2=%lf\n,fx2); if(fx2fx0) a0=a0; else a0=-a0; } else a0=-a0;//Y 正反方向 ; for(j=0; ;) {for(i=0;i2;i++) x2[i]=x0[i]+a0*e2[i]; printf(x2[0]=%f x2[1]=%f\n,x2[0],x2[1]); if(g1(x2[0],x2[1])=0 g2(x2[0],x2[1])=0 g3(x2[0],x2[1])=0) {fx22[j]=minfun(x2[0],x2[1]); printf(fx2=%lf\n,fx22[j]); if(fx22[j]fx0) { fx0=fx22[j]; x0[0]=x2[0]; x0[1]=x2[1]; a0=2*a0; j++; }// 实现 Y 方向的多次搜索； else break; } else break; } if(h0e) h0=0.5*h0; else break; } printf(%lf %lf\n,x0[0],x0[1]); printf(fx[min]=%lf\n,fx0); }// 输出结果 ,开始判断出去 .