gusucode.com > 学会不同搜索法(最速下降法,共轭梯度法,牛顿法,拟牛顿法)的计算步骤 > code/newton.m
%牛顿法极小化函数的通用子函数zuisu.m %输入变量为初始的迭代点,输出变量为极小值点 function x0=newton(x) %判断梯度范数是否满足计算精度的要求.是,标志变量设为1,输出结果;否,标志变量设为0,继续迭代 if sum(abs(tidu(x)).^2)<1e-4 flag=1; x0=x; else flag=0; end %循环求解函数的极小点 while flag==0 d=-inv(hesse(x))*tidu(x)'; a=gold(x,d); x=x+a*d'; %判断梯度范数是否满足计算精度的要求.是,标志变量设为1,输出结果;否,标志变量设为0,继续迭代 if sum(abs(tidu(x)).^2)<1e-4 flag=1; x0=x; else flag=0; end end