Numpy Scipy-求解形式为：A*dy/dt=f（t，y）的常微分方程组

在scipy ODE系统中，我需要解决哪些问题：

显然，如果A是可逆的，那么我可以等价地求解：

---

使用scipy.integrate中的标准方法。但仍然存在以下问题：效率如何？如果A是单数且不能反转怎么办？

是

fun

线性的吗？一般来说不是，但在我遇到的大多数问题中，它实际上是K*y+fIn的形式。一般情况下，你基本上有一个受代数约束的ODE集。在您的特定形式中，这个问题在控制文献中称为描述符系统。A和K的特征值起作用，例如，参见微分代数方程组

scipy

没有DAE解算器，您必须检查其他软件包，如DSTools和Sundals。谢谢您的提示！目前，我最终使用了assimulo，它为Sundails IDA解算器提供了pythonic接口。我也要试试petsc4py。也存在scikits.odes，但它与日晷的接口似乎不足以满足我（大规模）的需要。