gusucode.com > 学会不同搜索法(最速下降法,共轭梯度法,牛顿法,拟牛顿法)的计算步骤 > code/hesse.m
%hesse矩阵子函数,hesse.m,输入变量为四维向量,输出变量为hesse矩阵在x的数值形式 function h=hesse(x) %用符号表达形式来表示待求得函数 syms x1 x2 x3 x4 f='(x1+10*x2)^2+5*(x3-x4)^2+(x2-2*x3)^4+10*(x1-x4)^4'; g=[diff(f,'x1') diff(f,'x2') diff(f,'x3') diff(f,'x4')]; %函数f的hesse矩阵的符号表达形式 h=jacobian(g,[x1 x2 x3 x4]); x1=x(1);x2=x(2);x3=x(3);x4=x(4); h=eval(h);