有界面的程序.docVIP

有界面的程序 #includemath.h #define pi 3.1416 /*pp密度,pp1为p;pp2许用应力;pp3为焊缝成形系数;double ;*/ double pp1,pp2,V,pp,pp3; double objfx(double x[]) { double ff; ff=pp*x[1]*(x[0]+x[1])*4*V/x[0]/x[0]-0.67*pi*pp*x[0]*x[1]*(x[0]+x[1])+1.5*pp*pi*x[1]*(x[0]+x[2])*(x[0]+x[2]); return ff; } #includestdio.h void main() {printf(输入密度值(kg/m^3)\n); scanf(%lf,pp); printf(输入体积值(m^3)\n); scanf(%lf,V); printf(输入压力值(kpa)\n); scanf(%lf,pp1); printf(输入许用应力值(kpa)\n); scanf(%lf,pp2); printf(输入焊缝系数\n); scanf(%lf,pp3); void comple(int n,int k,int kg,double ep,double x[],double bl[],double bu[], double xcom[][100],double *f); double a[]={0,0,0},b[]={10,9,9},f,x[3],xcom[3][100]; comple(3,4,8,0.00001,x,a,b,xcom,f); printf( \n\n\n 输出最优解及目标函数值：\n); printf(\n x1=%f x2=%f x3=%f f(x1,x2,x3)=%f\n\n ,x[0],x[1],x[2],f); getchar(); getchar(); } /*wsd说明 n 优化模型维数 k复合形顶点数 kg 约束函数个数 ep 收敛精度 x bl初始点下限 bu初始点上限 xcom进程中的优化结果 *f */ #include stdlib.h #include math.h #include stdio.h double objfx(double x[]); void constraint(double x[],double g[]); void constraint(double x[],double g[]) { g[0]=pp1*x[0]-(2*pp2*pp3-pp1)*x[1]; g[1]=pp1*x[0]-(4*pp2-pp1)*x[2]; g[7]=12*x[0]*x[0]*x[0]+24*x[1]*x[0]*x[0]-1.27*V+2*x[0]*x[0]*x[0]; g[6]=1.27*V-0.67*x[0]*x[0]*x[0]-28*x[0]*x[0]*x[0]-56*x[1]*x[0]*x[0]; g[5]=2*x[1]-0.5*x[0]; g[2]=-x[0]; g[3]=-x[1]; g[4]=-x[2]; } int gau(double x[],double g[],int kg) { int i; constraint(x,g); for(i=0;ikg;i++) { if(g[i]0) goto s33; } return 1; s33:return 0; } void xcent(int n,int k,int ll,int lh,double x0[],double xcom[][100]) { int i,l; double 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,double x[],double xcom[][100],double fx