ceq=[];
在命令窗口键入：
x(0)=[13,21,0.8,42,12,16]; A=; b=;
LBnd=[];
[x,f]=fmincon(‘myfun’,x(0),A,b,[],[],LBnd,[],’constraint’) 优化结果为：
2、使用复合型法进行优化： 在命令窗口键入：
f=0.78539815*(4.75*x(1)*x(2)^2*x(3)^2+85*x(1)*x(2)*x(3)^2-85*x(1)*x(3)^2+0.92*x(1)*x(6)^2-x(1)*x(5)^2+0.8*x(1)*x(2)*x(3)*x(6)-1.6*x(1)*x(3)*x(6)+x(4)*x(5)^2+x(4)*x(6)^2+28*x(5)^2+32*x(6)^2);
g=[x(2)*x(3)-33,16-x(1)/x(3),x(1)/x(3)-35,44163/(x(2)*x(3)*sqrt(x(1)))-855,-261+7098/(x(1)*x(2)*x(3)^2*(0.169+0.6666/100*x(2)-0.854/10000*x(2)^2)),-213+7098/(x(1)*x(2)*x(3)^2*(0.2824+0.177/100*x(2)-0.314/10000*x(2)^2)),-0.03*x(4)+0.01233*x(4)^3/x(2)/x(3)/(x(5)^4),-55+1/(x(5)^3)*sqrt((29050*x(4)/x(3)/x(2))^2+(0.58*27300)^2),-55+1/(x(6)^3)*sqrt((29050*x(4)/x(3)/x(2))^2+(0.58*5*27300)^2),0.2-x(3),10-x(5),13-x(6),x(5)-15,x(6)-20,4-x(4)+x(1)+0.5*x(6)]; X=;
[x,minf] = minconSimpSearch(f,g,x,1.2,0.5,8,0.3); 所得结果为：
3、使用内点罚函数法进行优化：