python - scipy.integrate.solve_ivp 不清楚如何求解形式为 0=F(t, y(t), y'(t)) 的隐式 ODE

Question

目前，我确实使用 assimulos solver suite求解形式为 0 = F(t, y(t), y'(t)) 的隐式微分方程。我想使用 native scipy 安装附带的求解器，并查看了 scipy.integrate.solve_ivp 的文档.文档说解算器解决了以下形式的系统:

Solve an initial value problem for a system of ODEs. This function numerically integrates a system of ordinary differential equations given an initial value dy / dt = f(t, y) y(t0) = y0

我确实怀念解决 dy/dt = f(t, y, y') 形式系统的能力。这可能与 scipy.integrate.solve_ivp 吗？

Answer 1

这样的方程称为 differential algebraic equation (DAE)。据我所知，SciPy 没有内置的 DAE 求解器。

This answer至 this similar question提供了使用 fsolve 结合 odeint 解决问题的代数部分的示例。这可能适用于您的问题。

您可能还对 ODES package 感兴趣，它提供了额外的求解器作为 SciPy 的扩展，包括 DAE 求解器。