哪位帮忙看下这个程序（牛顿法求极值）

nvwu

这个程序能运行，但结果不对，不知道问题在哪

1. x=[-1.9,2]'
   
2. f=100* (x(2)-x(1)^2)^2+(1-x(1))^2;
   
3. grade=[ -400*(x(2)-x(1)^2)*x(1)-2+2*x(1);  200*x(2)-200*x(1)^2];
   
4. hessian=[1200*x(1)^2-400*x(2)+2,-400*x(1);-400*x(1),200];
   
5. invhessian=inv(hessian);
   
6. dk=-invhessian*grade;
   
7. x=x+dk
   
8. n=1;
   
9. while(abs(norm(grade))>=1e-6)&(n<10000)
   
10.    xnew=x
    
11.    grade=[ -400*(x(2)-x(1)^2)*x(1)-2+2*x(1);  200*x(2)-200*x(1)^2]
    
12.    dk=-inv([1200*x(2)^2-400*x(2),-400*x(1);-400*x(1),200])*grade;
    
13.    x=xnew+dk;
    
14.     n=n+1
    
15. end
    
16. f=100*(x(2)-x(1)^2)^2+(1-x(1))^2

复制代码

nvwu

编得不怎样，希望高手能指点一下　

xiejin17

你把原题目也发上来啊～，里面这么多数字～会不会是数字错了啊～

xuebx

能将原型说明白么，你想得到的结果又是什么呢？

xuebx

要不然无法对比分析阿

wei343

最好是把算法和结果传上来。

nvwu

原题是求香蕉函数极小值f=100*(x(2)-x(1)^2)^2+(1-x(1))^2
初始点x0= [-1.2 1]'（上面写成了x=[-1.9,2]'）.   其最小值为x*=[1 1]',f(x*)=0
算法是利用x(k+1)=x(k)-grad(f(x))/hessian(f(x))直到梯度 gradef(x)<误差 停止迭代

branch

几条建议:
注意使用分号,不要给我刷屏;
   我想你可能想看每一步修正的结果(刚好结果是2*1矩阵,可以用2*n矩阵保存结果);
语句尽量简单如第一部分(循环之前有很多重复),xnew也是重复,这样给阅读带来不便;
对于norm他好象就有求模的最大值的意思;

nvwu

谢谢