相关疑难解决方法(0)

这是情况:我有一个符号函数lamb,它是变量z的元素和变量h的函数元素的函数.这是羊羔象征功能的图像

.

现在我想根据变量eta和xi计算此函数的Gradient和Hessian.当然我用谷歌搜索它,但我找不到直接的方法来做到这一点.我发现的是这里,但正如我所说,它并不是这种情况的最佳方法.任何的想法？贝娄,源代码.谢谢.

from sympy import Symbol, Matrix, Function, simplify

eta = Symbol('eta')
xi = Symbol('xi')

x = Matrix([[xi],[eta]])

h = [Function('h_'+str(i+1))(x[0],x[1]) for i in range(3)]
z = [Symbol('z_'+str(i+1)) for i in range(3)]

lamb = 0
for i in range(3):
    lamb += 1/(2*sigma**2)*(z[i]-h[i])**2
simplify(lamb)