?? fsxdnewton.m
字號:
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;
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -