机械优化设计复合型法求解

1、复合 型 法 实 验 报 生 告.复合形法求约束思想在可行区内构造一个具有k个顶点的初始复合形.对该复合形的 顶点的日标函数值进行比较，找到日标函数值最大的顶点（称最坏点）, 然后按一定的法则求出日标函数值的所有下降的可行的新点,并用此 点代替最坏点，构成新的复合形，复合形的形状每改变一次，就向最优 点移动一步，直至逼近最优点.二.程序框图三.程序#include double objfx(double x)double ff;ff=(x0-5)*(x0-5)+4*(x1-6)*(x1-6);return ff;#include void main()void comple(int n,int

2、 k,int kg,double ep,double x,double bl,double bu, double xcom100,double *f);double a=0,0,b=10,20,f,x2,xcom2100;comple(2,3,3,0.00001,x,a,b,xcom,&f);printf( nnn输出最优解及目标函数值：n);printf(n x1=%.5fx2=%.5ff(x1,x2)=%.5fnn ,x0,x1,f);#include #include #include double objfx(double x);void constraint(double x,dou

3、ble g);void constraint(double x,double g)g0=64-x0*x0-x1*x1;g1=x1-x0-10;g2=x0-10;int gau(double x,double g,int kg)int i;constraint(x,g);for(i=0;i0)goto s33;return 1;s33:return 0;void xcent(int n,int k,int ll,int lh,double x0,double xcom100)int i,l;double xs;for(i=0;in;i+)xs=0;for(l=0;l-1)x0i=xs/(ll-1

4、);elsex0i=xs/ll;void fxse(int n,int k,double x,double xcom100,double fxk)int l,lp,lp1,i;double temp;for(l=0;lk-1;l+)for(lp=0;lpk-l;lp+)lp1=lp+1;if(fxklp=fxklp1)temp=fxklp;fxklp=fxklp1;fxklp1=temp;for(i=0;in;i+)xi=xcomilp;xcomilp=xcomilp1;xcomilp1=xi;void comple(int n,int k,int kg,double ep,double x,

5、double bl,double bu, double xcom100,double *f)int i,iw,l,ll,lh,it;double fx,fx0,sdx,fxh,fxr,alp;double *x0=(double*)calloc(n,sizeof(double);double *xh=(double*)calloc(n,sizeof(double);double *xr=(double*)calloc(n,sizeof(double);double *fxk=(double*)calloc(k,sizeof(double);double *g=(double*)calloc(k

6、g,sizeof(double);s5:for(i=0;in;i+)xi=bli+rand()/40000.0*(bui-bli);iw=gau(x,g,kg);if(iw=0)goto s5;for(i=0;in;i+)xcomi0=xi;for(l=1;lk;l+)for(i=0;in;i+)xcomil=bli+rand()/50000.0*(bui-bli);lh=-1;for(ll=1;llk;ll+)xcent(n,k,ll,lh,xO,xcom);iw=gau(xO,g,kg)；if(iw=O)goto s5;for(i=0;in;i+)xi=xcomill+l;s24:iw=g

7、au(x,g,kg);if(iw=0)(for(i=0;in;i+)xi=x0i+0.5*(xi-x0i);goto s24;else(for(i=0;in;i+)xcomill+l=xi;for(l=0;lk;l+)(for(i=0;in;i+)xi=xcomil;fx=objfx(x);fxkl=fx;it=0;s14: it=it+1;printf(nnn N= %2dn,it);lh=-1;xcent(n,k,k,lh,x0,xcom);fx0=objfx(x0);iw=gau(x0,g,kg);for(i=0;in;i+)printf( X(%d)min=%.5f ,i+1,x0i)

8、;printf( Fmin=%.5fn,fx0);for(i=0;ikg;i+)printf( G(%d)min=%.5f,i,gi);sdx=0;for(l=0;lk;l+)sdx=sdx+(fx0-fxkl)*(fx0-fxkl);sdx=sqrt(sdx/(double)k);if(sdxep) goto s38;fxse(n,k,x,xcom,fxk);lh=0;s22: fxh=fxklh;for(i=0;in;i+)xhi=xcomilh;xcent(n,k,k,lh,x0,xcom);iw=gau(x0,g,kg);if(iw=0)goto s36;alp=1.3;s12: fo

9、r(i=0;i=fxh)if(alp1.0e-4)alp=alp*0.5;goto s12;lh=lh+1;if(lh3)goto s22;for(i=0;in;i+)xcomilh=xri;fxklh=fxr;goto s14;s36: for(i=0;in;i+)bli=xcomik;bui=x0i;goto s5;s38: for(i=0;in;i+)xi=x0i;*f=objfx(x);free(x0);free(xh);free(xr);free(g);free(fxk);四.程序结果N =2Fnin=49.80302G2)min=-9.59616Xnin=0.40384Gnin=-

10、11.46400Xnin=8.6?61GCDnin =-1.72624N =3Xnin=0.66460G(0)inin=-5.48670Xnin=8.30933Gl)nin=-2.35527Fnin=40.12775G2)min=-9.33540N =4Xnin=1.15847G0)nin=-6.53102Xnin=8.31?99Gl)nin=-2.84048Fnin=36.24963G2)nin=-8.84153N =5Xnin=l.37966Gnin=-3.70029Xnin=8.11152Gmin=-3.26814Fnin=30.94100Gnin=-8.62034N =6Xnin=2.05487Gnin=-5.69989Xnin=8.09181Gnin=-3.96306Fnin=26.17644Gnin=-7.94513N =7Xnin=2.44370Gnin=-3.13106Xnin=7.82045Gnin=-4.62325Fnin=19.79079Gnin=-7.5563