fsxdnewton.m

来自「非线性规划中的最速下降法,黄金分割法,阻尼牛顿法,牛顿切线法.」· M 代码 · 共 21 行

function x=fsxdnewton(f,e,a,b)
% fsxdnewton函数 阻尼newton法
% x=fsxdnewton(f,e,a,b)为输入函数 f为函数 e为允许误差 (a,b)为初始点

% fsx TJPU 2008.6.16

syms x1 x2; 
v=[x1,x2]; 
df=jacobian(f,v); 
df=df.'; 
Q=jacobian(df,v);
x0=[a,b]';
g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)}); %求梯度
Q1=subs(Q,{x1,x2},{x0(1,1),x0(2,1)});  %求hesee矩阵 
	while(norm(g1)>e) 
	p=-Q1\g1; 
	x0=x0+p; 
	g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)}); 
  Q1=subs(Q,{x1,x2},{x0(1,1),x0(2,1)});
  end; 
x=x0;