在Matlab中求解具有六个初始条件的二阶紧耦合非线性ODE系统

Question

我正在解决流体动力学方面的问题；特别是紧密耦合的非线性常微分方程。以下是我实际问题的按比例缩小版本。 过去，我已经多次解决了耦合极点系统，但是这种情况有所不同，因为一个变量的双导数取决于另一个变量的双导数。如何在ode45中实现它？我需要3 x 2 = 6个x，x点和x-ddot相对于t的时间（从0到2）的图。所有所需的初始条件的值为零。

任何帮助将不胜感激。如果ODE没有最后的双导数项，那么我本人也会这样做。