import sympy
import numpy
from sympy import ordered, Matrix
x1,x2=sympy.symbols('x1 x2')
f=x1**2-x1*x2-4*x1 x2**2-x2
X0=numpy.array([[1],[1]])
v = list(ordered(f.free_symbols))
gradient = lambda f, v: Matrix([f]).jacobian(v)
gradf=sympy.transpose(gradient(f, v))
gradfx0=(gradf.subs([(x1, X0[0]), (x2, X0[1])]))
print(gradfx0)
我想計算兩個變數函式在python. 我定義了函式并找到函式 ( grad) 的梯度向量。現在，當我嘗試替換X0為grad時，結果是
Matrix([[2*x1 - x2 - 4], [-x1 2*x2 - 1]]).
我想要的結果應該是
Matrix([[-3], [0]]).
如何將點替換為 sympy 陣列？
uj5u.com熱心網友回復：
使用.subs和derive_by_array：
from sympy import symbols
from sympy.tensor.array import derive_by_array
x1, x2 = symbols('x1 x2')
f = x1 ** 2 - x1 * x2 x2 ** 2 - 4 * x1 - x2
grad = derive_by_array(f, (x1, x2))
# = [2*x1 - x2 - 4, -x1 2*x2 - 1]
gradx0 = grad.subs({x1: 1, x2: 1})
# = [-3, 0]
如果您想先說明您的觀點x0，請使用變數x0：
x0 = (1, 1)
gradx0 = grad.subs(zip((x1,x2), x0))
# = [-3, 0]
轉載請註明出處，本文鏈接：https://www.uj5u.com/qiye/419285.html
標籤：
上一篇：當運算元之一為NaN時，哪個算術運算子不會產生NaN結果？
下一篇：CRC不回傳預期值