应这样来求。
x0=[0;0];
A=[];b=[];
Aeq=[];beq=[];
lb=[0.01;0.01];ub=[];
[x,fval]=fmincon(@myobj,x0,[],[],[],[],lb,[],@mycon)
%目标函数
function f=myobj(x)
f=-pi*x(1)^2*x(2);
end
%约束函数
function [g,ceq]=mycon(x)
g=[];
ceq=300-(2*pi*x(1)*x(2)+2*pi*x(1)^2);
end
运行结果
x1=3.9894,x2= 7.9788,最大值f=398.9423