混合法求極值程序C語言

1、優(yōu)化設計作業(yè) 混合懲罰函數(shù)法求極值 輪機研 李凱2011129026混合懲罰函數(shù)法求極值程序dv2,double r)#include#iiicludedouble dv2;double xl2;/目標函數(shù)double func(double xl.double x2)return (xl-2)+(x2-1 )*(x2-!);約束條件double lij(double xl,double x2)return (xl-2*x2+l);double gi(double xl,double x2)double g=0.25*xl*xl+x2*x2-l;if (g0)return 0;else ret

2、uin g;懲罰函數(shù)double fun(double xl,double x2,double r)double f=func(xl,x2);double h=hj(xhx2);double g=gi(xl,x2);return (frr/g+(g*g+h*h)/i);double x72;double inax.niiii,t,xli2,xe2.f6,fe.fli,fl;mt i,l,m,n.k=0;x00=xl;x0l=x2;xl0=2.5;xll=3；x 0=3；x2l=2.5;doffbr(i=0;i=2;i+)fi=fun(xi0,xilj);/*判斷好點、壞點*/max=niui=

3、flO;m=n=0:for(i=0;i=2;i+)max=fi; m=i; if(finun)mui=fi; n=i;xli0=xm0;xlil=xml;fli=fiin(xli0 ,xli 1 ,r);xe0=xn0;xel=xnl;fe=fun(xe0,xe 1 .r);/*求中心點、反射點*/foi(i=0;i2;i+)x3i=0.5*(x0i+xli+x 2i.xhi)；單純形法求最優(yōu)解void F(double xljoublex2,doublex4i=x3i+l*(x3國.xhi);)f3=fuii(x3O,x3l)；f!4=fim(x40,x4l,r);/*判斷除好點、壞點以外的

4、點*/fbr(l=0;lfe) /*反射點函數(shù)值 比好點差*/if(f!4fl) /*反射點函數(shù)值 比除壞點以外的點都差*/( if(f!4=fh) /*反射點函 數(shù)值比壞點還差，反射失敗*/x60=x30+0.5*(xh0-x30);x6l=x3l+0.5*(xhl-x3l);f!6=fbn(x60,x6l,r);if(f6fli) /* 收 縮失敗，縮小邊長*/for(i=0;i=2;i+)xi0=0.5*(xi0+xe0);xil=0.5*(xil+xel);else /* 收縮成功，形成新的單純矩陣*/xm0=x60;xml=x6l;else /*反射點僅比最壞點好*/xli0=x40

5、;xlil=x4l;x60=x30+0.5*(xli0-x30);x6l=x3l+0.5*(xlil-x3l);fl6=ftin(x60,x6l,r);fbr(i=0;ife)xm0=x40;xml=x4l; else xm0=x50;xml=x5l; fbr(i=0;ile-4);dvO=xnO;dvl=xnl;printf(MdvO=%f %firdvO,dvl);main()int n=l,c=10.i;double x2=( 1.5,2);double e=0.001 ,fl.f2,r=0.00000001 ,wl;dofbi(i=0;i=le-3)(i-=c*r;pnntf(Hi-=

6、%f ”,r);else break;wlule(l);w 1=func(x0 ,x 1 );pnnt 町 x 1 =%f.x2=%f,w 1 =%f ”,x0,xl,wl);運行結果:. 128000 xX =X .805968 . x2 . 000730 . 02 4000 dv0=X . 3381660 - 9*78X46t=8.192000d0=0 . 94923&0 933035dv =0-9492360.933035ht*8 - X92000 xX 949236 x2 =0 . 933035 bt*=G5.536000dv0=0.8404*760.9X4443klv =0 - 8 404760 - 9 A 4443.536000 x1 -0 . 8 4047G - x2 0 - .288000 du0=0.8255330-9AX96Xdv =0.825S330-9AX96Xht*S24.288000 xX -0 - 825533 x2 0 - 9 AX96X *=41.9 4 . 30400010=0 . 824X260 - 91X262du =0-82