用MATLAB计算hessian矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

谁有代码呀~~或者思路也好~~~{:2_27:}

小女子拜谢~~

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：Hessian MATLAB matla atlab HESS MATLAB 矩阵 Hessian

The Hessian of a scalar valued function f:Rn®R is the n×n matrix of second order partial derivatives of f. In MATLAB we can obtain the Hessian of f by computing the Jacobian of the Jacobian of f. Consider once again the function f(x,y)=(4x2-1)e-x2-y2.

>> syms x y real
>> Hf=jacobian(jacobian(f));
>> Hf=simple(Hf)
Hf =
[2*exp(-x^2-y^2)*(2*x+1)*(2*x-1)*(2*x^2-5),  4*x*y*exp(-x^2-y^2)*(-5+4*x^2)]
[4*x*y*exp(-x^2-y^2)*(-5+4*x^2), 2*exp(-x^2-y^2)*(-1+2*y^2)*(2*x+1)*(2*x-1)]

2# pertain

谢谢。我试试之后回来反馈

我没运行出来。。==

直接在MATLAB的面板上输入

jacobian(jacobian((4*x^2-1)*exp(-x^2-y^2)))

得到结果

ans =

[ 8*exp(-x^2-y^2)-32*x^2*exp(-x^2-y^2)-2*(4*x^2-1)*exp(-x^2-y^2)+4*(4*x^2-1)*x^2*exp(-x^2-y^2),                                          -16*x*y*exp(-x^2-y^2)+4*(4*x^2-1)*x*y*exp(-x^2-y^2)]
[                                          -16*x*y*exp(-x^2-y^2)+4*(4*x^2-1)*x*y*exp(-x^2-y^2),                                     -2*(4*x^2-1)*exp(-x^2-y^2)+4*(4*x^2-1)*y^2*exp(-x^2-y^2)]

似乎不是矩阵。不过数目应该是对的了·~~~

下面一个问题是能apply subs 命令吗？

后面应该有变量吧。
jacobian(jacobian((4*x^2-1)*exp(-x^2-y^2), [x y ]), [x y])