exact.m

来自「共轭梯度法」· M 代码 · 共 22 行

x0=[0, -1];
counts=0;
sigma=0.001;
F=sym('100*(x2-x1^2)^2+(1-x1)^2');
x=x0';
[f g]=myfun(x0)
p0=-g;
while norm(g)>sigma 
    norm(g)
    x0=x;
    f0=f;
    g0=g;
    %tmin=linesearch(x0,p0)
    s= newton(F,x0,p0)
    x=x0+s*(p0)
    [f g]=myfun(x')
    alpha=(g'*g)/(g0'*g0)
    p=(-g)+alpha*p0
    p0=p;
    counts=counts+1
end