混合法求极值程序C语言.docVIP

优化设计作业 混合惩罚函数法求极值 轮机研 李凯 2011129026 混合惩罚函数法求极值程序  #includestdio.h #includemath.h double dv[2]; double x1[2]; //目标函数 double func(double x1,double x2) { return ((x1-2)*(x1-2)+(x2-1)*(x2-1)); } //约束条件 double hj(double x1,double x2) { return (x1-2*x2+1); } double gi(double x1,double x2) { double g=0.25*x1*x1+x2*x2-1; if (g0) return 0; else return g; } //惩罚函数 double fun(double x1,double x2,double r) { double f=func(x1,x2); double h=hj(x1,x2); double g=gi(x1,x2); return (f+r/g+(g*g+h*h)/r); } //单纯形法求最优解 void F(double x1,double x2,double dv[2],double r) { double x[7][2]; double max,min,t,xh[2],xe[2],f[6],fe,fh,fl; int i,l,m,n,k=0; x[0][0]=x1; x[0][1]=x2; x[1][0]=2.5; x[1][1]=3; x[2][0]=3; x[2][1]=2.5; do { for(i=0;i=2;i++) f[i]=fun(x[i][0],x[i][1],r); /*****判断好点、坏点*****/ max=min=f[0]; m=n=0; for(i=0;i=2;i++) { if(f[i]max) {max=f[i]; m=i;} if(f[i]min) {min=f[i]; n=i;} } xh[0]=x[m][0]; xh[1]=x[m][1]; fh=fun(xh[0],xh[1],r); xe[0]=x[n][0]; xe[1]=x[n][1]; fe=fun(xe[0],xe[1],r); /*****求中心点、反射点*****/ for(i=0;i2;i++) { x[3][i]=0.5*(x[0][i]+x[1][i]+x[2][i]-xh[i]); x[4][i]=x[3][i]+1*(x[3][i]-xh[i]); } f[3]=fun(x[3][0],x[3][1],r); f[4]=fun(x[4][0],x[4][1],r); /*****判断除好点、坏点以外的点l*****/ for(l=0;l=2;l++) if(l!=ml!=n)break; fl=fun(x[l][0],x[l][1],r); /*****分类判断*****/ if(f[4]fe) /*****反射点函数值比好点差*****/ if(f[4]fl) /*****反射点函数值比除坏点以外的点都差*****/ { if(f[4]=fh) /*****反射点函数值比坏点还差，反射失败*****/ { x[6][0]=x[3][0]+0.5*(xh[0]-x[3][0]); x[6][1]=x[3][1]+0.5*(xh[1]-x[3][1]); f[6]=fun(x[6