fun084.m

约束最优化

%min  f=100(x2-x1^2)^2+(1-x1)^2
%s.t. x1^2+x2^2<=1.5
%     x1+x2>=0

function  [f,df]=fun084(x)
f=100*(x(2)-x(1)^2)^2+(1-x(1))^2;
if nargout >1
df=[-400*x(1)*(x(2)-x(1)^2)-2*(1-x(1)),200*(x(2)-x(1)^2)]';
end