fsxhesse.m

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

function x=fsxhesse(f,a,b)
% fsxhesse函数 求函数的hesse矩阵；
% 本程序仅是简单的求二次函数的hesse矩阵！；
% x=fsxhesse(f)为输入函数 f为二次函数 x1,x2为自变量；

% fsx TJPU 2008.6.15

x1=a;x2=b;
fx=diff(f,'x1');     %求f对x1偏导数
fy=diff(f,'x2');     %求f对x2偏导数
fxx=diff(fx,'x1');   %求二阶偏导数 对x1再对x1
fxy=diff(fx,'x2');   %求二阶偏导数 对x1再对x2
fyx=diff(fy,'x1');   %求二阶偏导数 对x2再对x1
fyy=diff(fy,'x2');   %求二阶偏导数 对x2再对x2
fxx=subs(fxx);       %将符号变量转化为数值
fxy=subs(fxy);
fyx=subs(fyx);
fyy=subs(fyy);
x=[fxx,fxy;fyx,fyy]; %求hesse矩阵