使用PyDSTool

Question

我正在尝试使用PyDSTool包解决具有多个自由度（DOF）的动态系统。 ODE定义为：

  

M ddx + D dx + K x = f，

其中dx和ddx是x的一阶和二阶导数。为了简单起见，我首先要解决一个2 DOF系统（以及后来的大型系统），其中包含以下矩阵：

M = [[10,  0],
     [ 0,  8]]

K = [[  4000000.,  -3000000.],
     [ -3000000.,  10000000.]]

D = [[  400.,  -300.],
     [ -300.,  1000.]]

f = [[300 * sin(100*t)],
     [0]]

并且位移矢量定义为：

x = [[x1],
     [x2]]

我尝试了文档中的1 DOF示例，可以找到here，但是我对如何将其扩展到更大的系统感到有点迷茫。有人可以帮我解决多自由度案件吗？