python - 基于sympy中向量导数的向量函数的导数

Question

我想使用 sympy 计算基于 Python 中向量的导数的向量函数的导数。例如，如果函数是：

f(g) = g(t) / (gT(t) * g(t)) ^ 0.5

然后我想要：

df/dt = [dg(t) / dt * (gT(t) * g(t)) ^ 0.5 - 0.5 * g(t) * (dgT(t) / dt * g(t) + gT(t) * dg(t) / dt) * (gT(t) * g(t)) ^ (-0.5)] / [gT(t) * g(t)]

其中 gT(t) 是向量函数 g(t) 的转置。我怎样才能做到这一点？

Answer 1

导数可以用

>>> from sympy.abc import t
>>> g = Function('g', commutative=False)
>>> gT = Function('gT', commutative=False)
>>> f = g(t) / sqrt(gT(t) * g(t))
>>> f.diff(t).simplify()
(-g(t)*(gT(t)*g(t))**(-1/2)*(gT(t)*Derivative(g(t), t) + Derivative(gT(t), t)*g(t))*(gT(t)*g(t))**(-1) + 2*Derivative(g(t), t)*(gT(t)*g(t))**(-1/2))/2