优化设计--复合形法.docx

优化设计作业 复合形法 题目：课本例1-1。 程序如下： #include stdlib.h #include math.h #include stdio.h void com(int n,int k,int kg,float ep,float x[],float bl[],float bu[], float xc[][100],float *f) ; void constraint(float x[],float g[]); void xcent(int n,int k,int ll,int lh,float x0[],float xcom[][100]); void fxse(int n,int k,float x[],float xcom[][100],float fxk[]); float objfx(float x[]); int gau(float x[],float g[],int kg); float objfx(float x[]) { float f1; f1=120*x[0]+x[1]; return f1; } void constraint(float x[],float g[]) { g[0]=-x[0]+0.1; g[1]=-x[1]+1; g[2]=x[0]-3; g[3]=x[1]-8; g[4]=1-0.25*x[1]; g[5]=1-7*x[0]*x[1]/45; g[6]=1-7*x[0]*x[0]*x[0]*x[1]/45; } int gau(float x[],float g[],int kg) { int i; constraint(x,g); for(i=0;ikg;i++) { if(g[i]0) return 0;; } return 1; } void xcent(int n,int k,int ll,int lh,float x0[],float xcom[][100]) { int i,l; float xs; for(i=0;in;i++) { xs=0; for(l=0;lll;l++) { if(l!=lh) xs=xs+xcom[i][l]; } if(lh-1) x0[i]=xs/(ll-1); else x0[i]=xs/ll; } } void fxse(int n,int k,float x[],float xcom[][100],float fxk[]) { int l,lp,lp1,i; float te; for(l=0;lk-1;l++) for(lp=0;lpk-l;lp++) { lp1=lp+1; if(fxk[lp]=fxk[lp1]) { te=fxk[lp]; fxk[lp]=fxk[lp1]; fxk[lp1]=te; for(i=0;in;i++) { x[i]=xcom[i][lp]; xcom[i][lp]=xcom[i][lp1]; xcom[i][lp1]=x[i]; } } } } void com(int n,int k,int kg,float ep,float x[],float bl[],float bu[], float xcom[][100],float *f) { int i,iw,l,ll,lh,it; float fx,fx0,sdx,fxh,fxr,alp; float *x0=(float*)calloc(n,sizeof(float)); float *xh=(float*)calloc(n,sizeof(float)); float *xr=(float*)calloc(