使用 Runge-Kutta 求解数值四耦合微分方程答案

【问题标题】：Solving Numerical four coupled differential equations using Runge-Kutta使用 Runge-Kutta 求解数值四耦合微分方程
【发布时间】：2020-09-25 20:02:26
【问题描述】：

如何在 Python 中数值求解以下四个耦合方程？

我需要源代码！！！！！！

y'(1) = (-wa+wp)*y(2)+g*y(4)
y'(2) = (-wa+2wp)*y(1)-g*y(3)
y'(3) = (wa-wp)*y(4)+g*y(2)
y'(4) =-(wa+wp)*y(3)-g*y(1)

【问题讨论】：

您对数值求解器有何经验？是整个概念，还是只是 Python 语言？使索引移动，以便数组从索引0? 开始
我没有数值解的经验
您是刚刚了解了微分方程的整个主题，还是在遥远的过去，您最近发现自己有必要找到一个数值解？您是否认识到系统是线性的，并且知道您可以期待指数增长或下降的解决方案，可能是振荡的？
我对微分方程有足够的了解，现在我只想用数值求解上述方程。
根据您找到的关于 scipy.integrate 的答案，还有 scikit 和 pyDStool 等（jitcode 作为 scipy 的扩展，...）。从文档中，您是否了解如何对一阶系统进行编码以便将其传递给求解器例程？

标签： python numeric differential-equations runge-kutta

【解决方案1】：

查看Integration and ODEs (scipy.integrate)。您使用的具体求解器将取决于您拥有的边界条件：特别是对于初始值问题的 solve_ivp 和对于边界值问题的 solve_bvp。

【讨论】：