Search code examples
pythonvectorsimulinkodeodeint

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


I'm trying to "translate" a Simulink model (built to solve a second order ODE and shown in the picture below) into python.

enter image description here

inW and inT are vectors of dimensions [1601x2], gains MFW and MFT have dimension [2x1] and gains MC and MK have dimension [2x2]. All initial conditions have been set to 0.

out_acc, out_rate and out_disp all have dimension of [1601x2].

Any help in setting up this problem in python?

Thank you!

\Simone.


Solution

  • Solved with a for loop that updates the initial condition at each time step!