p2.m

在matlab环境下

function fun=p2(x)
g(1)=85.334407+0.0056858*x(2)*x(5)+0.00026*x(1)*x(4)-0.0022053*x(3)*x(5);
g(2)=92-(85.334407+0.0056858*x(2)*x(5)+0.00026*x(1)*x(4)-0.0022053*x(3)*x(5));
g(3)=80.51249+0.0071317*x(2)*x(5)+0.0029955*x(1)*x(2)+0.0021813*x(3)^2-90;
g(4)=110-(80.51249+0.0071317*x(2)*x(5)+0.0029955*x(1)*x(2)+0.0021813*x(3)^2);
g(5)=9.300961+0.0047026*x(3)*x(5)+0.0012547*x(1)*x(3)+0.0019085*x(3)*x(4)-20;
g(6)=25-(9.300961+0.0047026*x(3)*x(5)+0.0012547*x(1)*x(3)+0.0019085*x(3)*x(4));
fun=5.3578547*x(3)^2+0.8356891*x(1)*x(5)+37.293239*x(1)-40792.141;
fun=-fun;
g_m=min(g);
g_min=min(0,g_m);
fun=fun+fun*g_min;