matlab - 如何在 Matlab 中的二阶 ODE 求解器中实现三个初始条件？

Question

免责声明：我是 matlab 新手。对于两个变量 y1 和 y2，我有一组两个耦合的 ODE。我有以下 IC： y1(0), y2(0), dy1(0)/dt, dy2(0)/dt, d^2y2(0)/dt^2 , d^2y1(0)/dt ^2

d^2y2/dt^2 = f(y2, dy2/dt, d^2y1/dt^2)

d^2y1/dt^2 = f(y1, dy1/dt, d^2y2/dt^2)

在 Matlab 中，二阶 ODE 的 ode 求解器的工作方式是它需要两个 IC：y1(0)、dy1(0)/dt 和 y2(0)、dy2(0)/dt，用于求解 y1 和 y2。

但我的问题需要三个 IC 才能实现二阶 ODE。如何实现这一点？

标签： matlabnumerical-methodsodenumerical-analysis