tc_downhill_newton.m

来自「％ Nevill等插值算法地实现 不会花费大家的时间的 不希望其他人没有帐号」· M 代码 · 共 22 行

function x=tc_DownHill_newton(x0,f)
%newton法的变形，newton下山法，含有下山因子w初始值为1
%f=inline('[x^3-x-1,3*x^2-1]');
ff=f(x0);
w=1;
while 1
    ff=f(x0);
    x1=x0-w*ff(1)/ff(2);
    fff=f(x1);
    x=x1;
    if abs(fff(1))<abs(ff(1))
        w=2*w;
    else if abs(fff(1))>abs(ff(1))
        w=0.5*w;
        continue;   
        end
    end
     if abs(fff(1)-ff(1))<1e-12
        return;
     end
    x0=x1;
end  %while