我发现即使我在代码中使用了 F.Jacobian('x')，在 Sympy 中计算符号矩阵的雅可比的速度也很慢。但是我发现当简单地用函数 diff 计算时，Symengine 比 Sympy 快得多。可悲的是，当我尝试使用 Symengine 进行 Jacobian 矩阵时，计算 Jacobian 矩阵的速度是相同的，因为我只能使用 for 循环来计算矩阵。
Sympy 代码：a = Fx.jacobian(change) （Fx 是目标矩阵，change 是变量列表）
符号引擎代码：
J = symengine.zeros( , ) # 初始化雅可比矩阵
对于范围内的 j (***)：
a = symengine.diff(Fx, change[j])
我认为如果 symengine 有雅可比函数，速度会快得多。