python - 求解一阶微分方程组和二阶微分方程组（非线性）

Question

问题

我目前有一个由四个方程组成的系统。两个是二阶微分方程，两个是一阶微分方程：

四个方程

初始条件为：

x = 0 |
y = 0.3 | 
f(x) = 2.05 |
f(y) = 0.55 |

所有角度都以度为单位。

我试过的

我曾尝试使用 Google Colabs 并使用 SciPy 和 NumPy。不幸的是，我无法弄清楚如何对其进行编程，因为这些方程是非线性的。有人可以给我一些关于使用其他模块的提示吗？

标签： pythonnumpygoogle-colaboratorydifferential-equations

Answer 1

Scipy 有一个Runge Kutta 求解器。首先，您必须将 ODE 转换为一阶系统（您始终可以通过设置 z=y' 来做到这一点），然后尝试 RK 求解器。