双摆

**双摆**是将一根单摆连接在另一个单摆的尾部所构成的系统。双摆同时拥有着简单的构造和复杂的行为。高能量双摆的摆动轨迹表现出对于初始状态的极端敏感。两个初始状态差异极小的双摆在一段时间的运行后表现非常不同，是一种具有混沌性质的简单动力系统。

## 分析以及詮釋
可以考慮許多不同種類的双摆：二個擺的長度及重量可能相同，也可能不同。二個擺可能都是單擺，也有可能是複擺（compound pendulum），其運動可能限制在二維空間，也可以在三維空間內進行。在以下的分析中，二個擺的擺長及質量都相同，運動限制在二維空間內。

複擺的質量假設是延著其長度均勻分佈，則其複擺的質心是在中點，複擺的臂對中點的轉動慣量為。

比較方便定義系統位形空间的方式是用複擺臂和垂直線之間的夾角為廣義座標。角度名稱為及。二桿質心的位置可以用二個座標表示。若笛卡尔坐标系的原點是在第一個擺（最上方擺）的固定點，則其第一個擺的質心在：
:\begin{align}
x_1 &= \frac{l}{2} \sin \theta_1 \\
y_1 &= -\frac{l}{2} \cos \theta_1
\end{align}
第二個擺的質心在：
:\begin{align}
x_2 &= l \left ( \sin \theta_1 + \tfrac{1}{2} \sin \theta_2 \right ) \\
y_2 &= -l \left ( \cos \theta_1 + \tfrac{1}{2} \cos \theta_2 \right )
\end{align}
上述資訊已經可以建立拉格朗日量（Lagrangian）。
### 拉格朗日力學
双摆系統的拉格朗日量為
:
\begin{align}L & = \text{kinetic energy} - \text{potential energy} \\
 & = \tfrac{1}{2} m \left ( v_1^2 + v_2^2 \right ) + \tfrac{1}{2} I \left ( {\dot \theta_1}^2 + {\dot \theta_2}^2 \right ) - m g \left ( y_1 + y_2 \right ) \\
 & = \tfrac{1}{2} m \left ( {\dot x_1}^2 + {\dot y_1}^2 + {\dot x_2}^2 + {\dot y_2}^2 \right ) + \tfrac{1}{2} I \left ( {\dot \theta_1}^2 + {\dot \theta_2}^2 \right ) - m g \left ( y_1 + y_2 \right ) \end{align}

第一項是質心的平移动能，第二項是擺延著質心旋轉的轉動動能，最後一項是双摆在均勻重心場下的势能。其點標示表示變數的时间导数。

將以上的座標代入，重組後可得
:
L = \tfrac{1}{6} m l^2 \left ( {\dot \theta_2}^2 + 4 {\dot \theta_1}^2 + 3 {\dot \theta_1} {\dot \theta_2} \cos (\theta_1-\theta_2) \right ) + \tfrac{1}{2} m g l \left ( 3 \cos \theta_1 + \cos \theta_2 \right ).

這裡只有一個守恆量（能量），沒有守恆的動量，二個廣義的動量可以表示為

:\begin{align}
p_{\theta_1} &= \frac{\partial L}{\partial {\dot \theta_1}} = \tfrac{1}{6} m l^2 \left ( 8 {\dot \theta_1} + 3 {\dot \theta_2} \cos (\theta_1-\theta_2) \right ) \\
p_{\theta_2} &= \frac{\partial L}{\partial {\dot \theta_2}} = \tfrac{1}{6} m l^2 \left ( 2 {\dot \theta_2} + 3 {\dot \theta_1} \cos (\theta_1-\theta_2) \right ).
\end{align}

上式可以求得

:\begin{align}
{\dot \theta_1} &= \frac{6}{ml^2} \frac{ 2 p_{\theta_1} - 3 \cos(\theta_1-\theta_2) p_{\theta_2}}{16 - 9 \cos^2(\theta_1-\theta_2)} \\
{\dot \theta_2} &= \frac{6}{ml^2} \frac{ 8 p_{\theta_2} - 3 \cos(\theta_1-\theta_2) p_{\theta_1}}{16 - 9 \cos^2(\theta_1-\theta_2)}.
\end{align}

運動方程式為

:\begin{align}
{\dot p_{\theta_1}} &= \frac{\partial L}{\partial \theta_1} = -\tfrac{1}{2} m l^2 \left ( {\dot \theta_1} {\dot \theta_2} \sin (\theta_1-\theta_2) + 3 \frac{g}{l} \sin \theta_1 \right ) \\
{\dot p_{\theta_2}} &= \frac{\partial L}{\partial \theta_2} = -\tfrac{1}{2} m l^2 \left ( -{\dot \theta_1} {\dot \theta_2} \sin (\theta_1-\theta_2) + \frac{g}{l} \sin \theta_2 \right ).
\end{align}

最後四個方程是有系統目前狀態時，系統隨時間演進的顯式方程。不太可能再進一步求得方程的積分解析解，得到和時間顯函數的解。不過利用龙格－库塔法或其他數值方式，可以進行數值積分來求解。

## 混沌運動

双摆的運動是混沌運動，且對初始條件非常敏感。右圖是雙擺在不同初始條件下，是否會翻倒（成為倒擺）的圖。其初始值的範圍是在方向的−3到3，而初始值的範圍是在方向的−3到3。點的顏色說明擺在以下時間內會翻倒：
* （綠色）
* （紅色）
* （紫色）
* （蓝色）

若在時間後，仍然不會翻倒，其顏色為白色。

中心白色區域的邊界可以依能量守恆推得，為以下的曲線：
:3 \cos \theta_1 + \cos \theta_2 = 2. 

因此若

:3 \cos \theta_1 + \cos \theta_2 > 2, 

以能量的關係，双摆不可能翻倒。在此區域外，以能量來說，双摆有可能翻倒，但是否會翻倒本身是很複雜的問題。若雙擺的末端是點質量，不是質量均勻分佈的桿子，情形類似。

雙擺沒有自然共振頻率，因此可用在大樓中，大樓本身是主要的倒擺，而上面又有一個質量，形成倒雙擺。

## 相关條目
*单摆
*倒雙擺
*

## 参考資料

*
* Eric W. Weisstein, *[http://scienceworld.wolfram.com/physics/DoublePendulum.html Double pendulum] * (2005), ScienceWorld *(contains details of the complicated equations involved)* and "Double Pendulum " by Rob Morris, Wolfram Demonstrations Project, 2007 (animations of those equations).
* Peter Lynch, *[https://web.archive.org/web/20030608233118/http://www.maths.tcd.ie/~plynch/SwingingSpring/doublependulum.html Double Pendulum]*, (2001). *(Java applet simulation.)*
* Northwestern University, *[http://www.physics.northwestern.edu/vpl/mechanics/pendulum.html Double Pendulum] *, *(Java applet simulation.)*
* Theoretical High-Energy Astrophysics Group at UBC, *[https://web.archive.org/web/20070310213326/http://tabitha.phas.ubc.ca/wiki/index.php/Double_pendulum Double pendulum]*, (2005).
## 外部链接
*麻省理工学院TSG小组进行的双摆行为演示（实物）
*Youtube上的双摆轨迹演示（电脑模拟）