Matlab：ODE45 +拍摄方法似乎不起作用

Question

我试图在数值上同时解决三个ODE函数：

dw/ds = sin(theta)

dx1/ds = cos(theta)

dtheta/ds = - a * cos(theta) * (w-w0) + a * sin(theta) * x1

我正在解决一个被分成两半的列，我想对w1绘制w并且有一个正w（现在看来是不可能的）。主要有三个功能：

• first_sec_w：使用已经找到的theta初始值的RK4方法（使用shoot_first和shoot_second）。
• shoot_first：一个函数，使用RK4解决函数，假设θ的初始值，直到L0 / 2处的边界条件为0，并返回此初始值。
• shoot_second：几乎和上面一样，但它从L0开始直到L0 / 2。

任何人都可以帮助我吗？

我会上传代码，所以这里看起来不会太乱。

Indexing Strategies