Python中的微分方程

Question

我正在使用一个微分方程系统，想知道有没有哪个Python库是最常用来解决微分方程的。

我的方程是非线性的一阶方程。

Answer 1

如果你想进行符号计算，可以使用sympy这个工具。

Answer 2

你可以使用SciPy中的integrate，它可以和标准的LAPACK程序进行配合，来处理类似的事情。

可以看看这个教程，这是我在谷歌上找到的一个例子。这里是相关的文档。

Answer 3

如果你需要解决大型的非线性系统（特别是那些比较复杂的），那么使用scipy工具可能会比较慢，而且操作起来也不太方便。在这种情况下，PyDSTool这个软件包现在已经被很多人使用了。它可以把你的方程自动转换成C语言代码，并且能和优秀的求解器结合使用。如果你想定义一些状态事件，比如阈值交叉，或者从数组中添加外部输入信号，甚至进行其他分析（比如分岔分析，因为这个软件包还提供了与AUTO的接口），它都特别适合。