使用 SciPy 求解常微分方程

Question

我有一个追随者 ordinary differential equation 和数值参数 Sigma=0.4，x(0) = 4 和 dx(0)/dt = 0
我的任务是使用 ode 函数
得到微分方程的柯西问题解决方案（初始值问题解决方案） 有人能帮我吗？我什至不知道如何以正确的方式为 SciPy 编写方程式，尤其是数字参数。
P.S。不好意思发不了图，我刚注册。

Answer 1

正如沃伦所说，scipy.integrate.odeint 是 'SciPy' 解决这个问题的方法。

但在将问题提交给 SciPy（或您最终使用的任何求解器）之前，您需要使用类似以下方法将二阶 ODE 转换为一阶 ODE：http://tutorial.math.lamar.edu/Classes/DE/SystemsDE.aspx

为了让事情变成 SciPy 你需要让你的等式看起来像：

y' = f(y)

但是现在你的等式是这样写的：

y'' = f(y, y')

解决办法是给你的系统增加更多的变量，但是link会解释得更透彻。