优化设计-单纯形法

《优化设计-单纯形法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、£=0.005o1)用单纯形法求法min(x/+2x22・4xi・2xiX2),已知a=l,卩=0・5,y=2,程序如下：#include"stdio.h"#include"math.h"floatabs(float);floatmax(float,float);floatf(float*);floatremain();voidfs();voidxs();voidss();voidsx();voideh();intk=0;float*xh;float*xe;floatx0[2]={0,0},xl[2]={0,l},x2[2]={l,0},xnl[2],xn2[2],xn3[2],xn4[2]

2、,p[3];floatfh,fe,e=0.005;floatfO,fl,f2,fn1,fn2,fn3,fn4;voideh(){f()=f(x());fl=f(xl);f2=f(x2);if(fO>fl){if(fO>f2){if(fl>⑵{fh=fO;fc=f2;xh=xO;xc=x2;p[0]=0;p[1]=1；p[2]=1;}else{fh=fO;fe=fl;xh=xO;xe=x1;p[0]=0;p[l]=l;p[2]=1;J}else{fh=f2;fe=fl;xh=x2;xe=xl;p[0]=l;p[l]=l;p[2]=0;J}else{if(f0

3、f2;fe=f();xh=x2;xe=x();p[0]=l;p[ll=l;p[2]=0;}else{fh=fl;fc=fO;xh=x1;xc=xO;p[0]=1;p[l]=0;p[2]=1;}}else{fh=fl;fe=f2;xh=x1;xe=x2;p[0]=l;p[l]=0;p[2]=1;J}xnHOJ=(xOLOJ*p[OJ+xlL0J*pLl]+x2L0J*仇2[)/2;xnl[lJ=(xOll]*p[O]+x1[1]*p[1]+x2[1]*p[2])/2;fnl=f(xnl);}voidfs(){xn2[0J=xnl[0]+(xn1[0]-xh[0]);xn2[1]=xn1[1]+

4、(xn1[1]-xh[1]);fn2=f(xn2);}voidxs(){xn3[0J=xn1[0J+2*(xn2[0J-xn1[OJ);xn3[l]=xn1[l]+2*(xn2[l]-xn1[1]);fn3=f(xn3);}voidss(){xn4[0J=xn1[0]+0.5*(xh[0卜xn1[0]);xn4[l]=xn1[l]+0.5*(xh[l]-xn1[1]);fn4=f(xn4);}voidsx(){x0[0]=0.5*(x0[0]+xe[0]);x0[l]=0.5*(x0[l]+xe[l]);xl[0]=0.5*(xl[0]+xe[0]);xl[l]=0.5*(xl[1]+xe[

5、1]);x2l0J=0.5*(x2[0]+xe[0]);x2[l]=0.5*(x2[l]+xe[l]);)floatmax(floatx,floaty){if(x>y)retumx;elsereturny;}floatremain(){if(p[0]==0)rcturnmax(f1,f2);elseif(p[l]==O)returnmax(f0,f2);elseif(p[2J==0)returnmax(f0,fl);elsereturn0;}floatf(float*x){floatr;r=x[0]*x[0]+2*x[l]*x[l]-4*x[0]-2*x[0]*x[l];returnr;}v

6、oidmain(){floati,j,s;printf("kxOxlx2f(x(n+l))");eh();printf(n%dt[%.4f%.4f][%.4f%.4f][%.4f%・4f]%.4f,,,k,x0[0],x0[l],xl[0],xl[l],x2[0],x2[l],fnl);k++;do{step1:fs();if((fn2>fe)&&(fn2

7、3[1];fh=f(xh);gotostep4;}else{xh[01=xn2[0];xh[1]=xn2[1];fh=f(xh);gotostep4;}}else{if(fn2>remain());gotostep2;}stcp2:{if(fn2>=fh){ss();gotostep3;}else{xh[0]=xn2[0];xh[l]=xn2[l];fh=f(xh);ss();gotostep3;}}step