数值积分方程
\(\begin{cases} J_1\dot{\omega}_1-(J_1-J_3)\omega_2\omega_3=L_1\\ J_1\dot{\omega}_2+(J_1-J_3)\omega_3\omega_1=L_2\\ J_3\dot{\omega}_3=L_3 \end{cases}\)
\(L_1=mgl\sin\theta\cos\psi,\quad L_2=-mgl\sin\theta\sin\psi,\quad L_3=0\)
\(J_1=J_2\) 为横向主惯量,\(J_3\) 为对称轴惯量,\(mgl\) 为重力力矩系数。