examp084.m

来自「约束最优化」· M 代码 · 共 14 行

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

x0=[-1.9,2];
A=[-1 -1];b=0;
[x1,f1,exitflag1,output1,lambda1] = fmincon('fun084',x0,A,b,[],[],[],[],'con084')     
pause

opt=optimset('GradObj','on','GradConstr','on');
[x2,f2,exitflag2,output2,lambda2]=fmincon('fun084',x0,A,b,[],[],[],[],'con084',opt)
pause