用Mathematica模拟双摆

1、双摆系统在任意时刻的受力分析图。

3、把上面的四个微分方程，视为一个微分方程组：fcz = 辘腋粪梯{Subscript[m, 1] x1''[t] == (\[Lambda]1[t柯计瓤绘]/Subscript[l, 1]) x1[t] - (\[Lambda]2[t]/Subscript[l, 2]) (x2[t] - x1[t]), Subscript[m, 1] y1''[t] == (\[Lambda]1[t]/Subscript[l, 1]) y1[t] - (\[Lambda]2[t]/Subscript[l, 2]) (y2[t] - y1[t]) - Subscript[m, 1] g, Subscript[m, 2] x2''[t] == (\[Lambda]2[t]/Subscript[l, 2]) (x2[t] - x1[t]), Subscript[m, 2] y2''[t] == (\[Lambda]2[t]/Subscript[l, 2]) (y2[t] - y1[t]) - Subscript[m, 2] g};

5、作为微分系统，需要有明确的初始状态：begin={x1[0] == 1, y1[0] == 0, y1'[0] == 0, x2[0] == 1, y2[0] == -1, y2'[0] == 0};

7、这样就可以模拟双摆的摆动过程。

9、实时的画出第二个摆锤的移动轨迹。