如何用Python解决僵硬的ode？

Question

我是Python初学者。我正在尝试切换一些我在matlab中的程序。 我需要解一个方程，它的所有输入都是刚性矩阵。在matlab中我使用

[ttT,uT] = ode23s('SST',t,fT);

Answer 1

对于在Matlab中所做的大多数事情，都可以使用Python中的NumPy模块来完成。可以找到here。在

您可能还发现相关模块SciPy也很有用。在

PyDSTool也可能与您相关。它是雷达解算器的包装。在

然后您可以尝试使用matplotlib进行绘图。它的工作原理很像Matlab的绘图工具。在

以下链接也可能有帮助：

• http://www.ews.uiuc.edu/~mrgates2/ode/

• http://wiki.python.org/moin/NumericAndScientific?action=show&redirect=SciPy

• Integrate stiff ODEs with Python