python - python中具有耦合导数的微分方程

Question

我正在尝试使用 sympy 和 scipy 解决一组微分方程，但无法弄清楚如何以适当的形式将它们带入。简而言之，我有一组两个耦合的二阶微分方程，我可以将它们重新写成一个由四个一阶微分方程组成的系统，形式如下：

点(x1) = x2

点（x2）= f（x1，x2，x3，x4，点（x4））

点 (x3) = x4

点（x4）= f（x1，x2，x3，x4，点（x2））

我的问题是一阶导数的耦合。在 Matlab 中，我可以使用时间和状态变量相关的质量矩阵来解决这个问题，但在 python 中，我只看到了具有时间不变质量矩阵的方法。任何建议将不胜感激。

标签： pythonnumpyscipysympyenthought