目前，我确实使用assimulos solver suite来求解形式为0 = F（t，y（t），y'（t））的隐式微分方程 . 我想使用一个带有原生scipy安装的解算器，并查看scipy.integrate.solve_ivp的文档 . 文档说解算器解决了以下形式的系统：
解决ODE系统的初始值问题 . 给定初始值dy / dt = f（t，y）y（t0）= y0时，该函数对常微分方程组进行数值积分
我确实错过了解决 dy / dt = f(t, y, y') 形式的系统的能力 . 这可能是 scipy.integrate.solve_ivp 吗？
1 回答
这样的等式称为differential algebraic equation（DAE） . 据我所知，SciPy没有内置的DAE求解器 .
This answer至this similar question提供了使用
fsolve结合odeint解决问题的代数部分的示例 . 这可能适用于您的问题 .您可能也对ODES package感兴趣，它提供了额外的求解器作为SciPy的扩展，包括DAE求解器 .