Hessian在Python中使用有限差分方法

Question

可以说，我有一个函数z = x^2 + y^2。现在，我想在python中实现此功能z的粗麻布。到目前为止，我已经使用如下所示的有限差分方法计算了导数-

def derivative(x, y, f, h):
    return [(f(x + h, y) - f(x - h, y )) / (2*h), (f(x , y + h) - f(x , y - h)) / (2*h)]

我不确定如何重用它来计算Hessian。

Answer 1

有一个更简单的方法：使用numdifftools包。它具有Hessian函数，该函数可以执行您想要的操作。 Documentation