优化设计-单纯形法.docx

1) 用单纯形法求法min（x12+2x22-4x1-2x1x2），已知α=1，β=0.5，γ=2，ε=0.005。 程序如下： #include stdio.h #include math.h float abs(float); float max(float,float); float f(float*); float remain(); void fs( ); void xs( ); void ss( ); void sx( ); void eh( ); int k=0; float* xh; float* xe; float x0[2]={0,0},x1[2]={0,1},x2[2]={1,0},xn1[2],xn2[2],xn3[2],xn4[2],p[3]; float fh,fe,e=0.005; float f0,f1,f2,fn1,fn2,fn3,fn4; void eh( ) { f0=f(x0); f1=f(x1); f2=f(x2); if(f0f1) { if(f0f2) { if(f1f2) {fh=f0; fe=f2; xh=x0; xe=x2; p[0]=0;p[1]=1;p[2]=1;} else {fh=f0; fe=f1; xh=x0; xe=x1; p[0]=0;p[1]=1;p[2]=1;} } else { fh=f2; fe=f1; xh=x2; xe=x1; p[0]=1;p[1]=1;p[2]=0;} } else { if(f0f2) { if(f1f2){fh=f2; fe=f0; xh=x2; xe=x0; p[0]=1;p[1]=1;p[2]=0;} else {fh=f1; fe=f0; xh=x1; xe=x0; p[0]=1;p[1]=0;p[2]=1;} } else { fh=f1; fe=f2; xh=x1; xe=x2; p[0]=1;p[1]=0;p[2]=1;} } xn1[0]=(x0[0]*p[0]+x1[0]*p[1]+x2[0]*p[2])/2; xn1[1]=(x0[1]*p[0]+x1[1]*p[1]+x2[1]*p[2])/2; fn1=f(xn1); } void fs( ) { xn2[0]=xn1[0]+(xn1[0]-xh[0]); xn2[1]=xn1[1]+(xn1[1]-xh[1]); fn2=f(xn2); } void xs( ) { xn3[0]=xn1[0]+2*(xn2[0]-xn1[0]); xn3[1]=xn1[1]+2*(xn2[1]-xn1[1]); fn3=f(xn3); } void ss( ) { xn4[0]=xn1[0]+0.5*(xh[0]-xn1[0]); xn4[1]=xn1[1]+0.5*(xh[1]-xn1[1]); fn4=f(xn4); } void sx( ) { x0[0]=0.5*(x0[0]+xe[0]); x0[1]=0.5*(x0[1]+xe[1]); x1[0]=0.5*(x1[0]+xe[0]); x1[1]=0.5*(x1[1]+xe[1]); x2[0]=0.5*(x2[0]+xe[0]); x2[1]=0.5*(x2[1]+xe[1]); } float max(float x,float y) {if (xy)return x; else return y; } float remain() {if(p[0]==0)return max(f1,f2); else if(p[1]==0)return max(f0,f2); else if(p[2]==0)return max(f0,f1); else return 0; } float f(float* x) { float r; r=x[0]*x[0]+2*x[1]*x[1]-4*x[0]-2*x[0]*x[1]; return r; } void main() { float i,j,s;printf(\n k x0 x1 x2 f(x(n+1))\n\n); eh(); printf(%d\t[%.4f%.4f][%.4f%.4f] [%.4f %.4f] %.4f\n,k,x0[0],x0[1],x1[0],x1[1],x2[0],x2[1],fn1); k++; do {step1: fs(); if((fn2fe)(fn2remain())){x