rosen084.m

来自「matlab讲义」· M 代码 · 共 8 行

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

function  [df,dg]=rosen084(x)
df=[-400*x(1)*(x(2)-x(1)^2)-2*(1-x(1)),200*(x(2)-x(1)^2)];
dg=[2*x(1),-1;2*x(2),-1];