python计算hessian矩阵

要计算函数的Hessian矩阵，可以使用SymPy库中的hessian函数。以下是一个示例代码：

from sympy import symbols, Matrix, hessian
from sympy.abc import x, y

# 定义函数
f = x**2 + y**3

# 定义变量
vars = [x, y]

# 计算Hessian矩阵
H = hessian(f, vars)

# 打印结果
print("Hessian矩阵:")
print(Matrix(H))

输出结果为：

Hessian矩阵:
Matrix([[2, 0], [0, 6*y]])

其中，第一行第一列的元素为$f$对$x$的二阶偏导数，第一行第二列的元素为$f$对$x$和$y$的混合偏导数，第二行第一列的元素为$f$对$y$和$x$的混合偏导数，第二行第二列的元素为$f$对$y$的二阶偏导数。