weixin_56907349的博客

原创 matlab——Newton迭代极值出错

syms x1 x2 y=100*(x1-x2^2)^2+(1-x1)^2; pretty(y); %y1=diff(y,x1) %y2=diff(y,x2) y1=jacobian(y,[x1,x2]); y2=jacobian(y1,[x1,x2]); h=y1*inv(y2) %% x1=1; x2=0; y1=jacobian(y,[x1,x2]); y2=jacobian(y1,[x1,x2]); h=y1*inv(y2) 前半段没问题，我算出了Newton的迭代表达式，但是后半段无论是用