当状态向量由向量组成时,如何求解 python 中的二阶 ODE?

How can I solve a second order ODE in python when state vector is made of vectors?

我正在尝试 "translate" 将一个 Simulink 模型(用于求解二阶 ODE 并如下图所示)转换为 python。

inW 和 inT 是维度为 [1601x2] 的向量,增益 MFW 和 MFT 的维度为 [2x1],增益 MC 和 MK 的维度为 [2x2]。所有初始条件都设置为0.

out_acc、out_rate 和 out_disp 的维度均为 [1601x2]。

在 python 中设置这个问题有什么帮助吗?

谢谢!

\西蒙。

解决了在每个时间步更新初始条件的 for 循环!