在MATLAB ode45中，如何解决ODE系统，其中参数依赖于其他参数，而不依赖于时间？

Question

考虑运载火箭的方程式

dv / dt = T / m - D / m - g * sin（gam）

v * d（gam）/ dt = -g（-v ^ 2 / W）* cos（gam）其中 W是常数

dh / dt = v * sin（gam）。

其中v表示速度，T是推力。

我的推力值是离散的，因此我必须使用 interp1 ，这很简单。

然而，我的拖D **取决于**速度本身，我该如何合并？

如您所见，参数是另一个参数的函数。