将下面语句在matlab编辑器中保存为mycon.m
function [c, ceq]=mycon (x)
c=(300/210000*(1-x(1)/300))^3/x(2)^2-(300/210000*(1-x(1)/300))^3/(0.2+x(2))^2-0.000075;
ceq=[ ]; %无等式约束
在matlab窗口中输入:
fun=’-x(1)*(0.02+x(2))/(2*(300-x(1)))+300/210000’; x0=[1 1];A=[];b=[]; Aeq=[ ]; beq=[ ];
lb=[0 0];ub=[300 1]; [x,fval,exitflag,output,lambda,grad,hessian]=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,@mycon)