🪀 耦合振子與正常模態

兩顆質量各自的牛頓方程長這樣：

麻煩在於 $x_1$ 的方程裡有 $x_2$，反之亦然。 它們綁在一起，沒辦法一個一個解。寫成矩陣：

綁在一起的原因，就是 $K$ 的非對角元素不是 0。而線性代數早就給了答案： 對稱矩陣可以正交對角化。找到 $K$ 的特徵向量， 在那組座標下 $K$ 是對角的，非對角項消失 — 兩個方程就此分家。

🌟 這個矩陣的特徵向量

因為 $K$ 對兩顆質量完全對稱，特徵向量不用算就能猜出來：

• 模態 1（同相） — 兩顆一起往同一邊。中間那根彈簧 從頭到尾沒有被拉伸過，所以 $k_c$ 完全不出現在 ${\omega }_1$ 裡。把 k_c 滑桿拉到底，ω₁ 紋風不動 — 你可以現在就試。
• 模態 2（反相） — 兩顆往相反方向。中間彈簧被拉到最兇， 所以它額外貢獻了 $2k_c$，頻率比較高。

頻率就是特徵值開根號除以質量：${\omega }_i=\sqrt{{\lambda }_i/m}$。 這不是巧合 — 特徵值是「這個方向有多硬」，特徵向量是「哪個方向」。

⚖️ 對角化的觀測後果

把位移投影到兩個特徵向量上，得到模態座標 $q_1=(x_1+x_2)/\sqrt{2}$、$q_2=(x_1-x_2)/\sqrt{2}$， 總能量就乾淨地拆成兩份：

兩份各自守恆。右邊面板上那兩條「模態能量」不會動， 而「各質量的局部能量」兩條一直在互相搬 — 同一個系統、同一個瞬間，換個座標看， 「能量在流動」這件事就消失了。這就是對角化真正買到的東西。

🥁 拍是怎麼來的

選「只推第一顆」，初始條件同時含有兩個模態。兩個餘弦疊加，用和差化積：

快的那個 $\bar{\omega }$ 是你看到的震盪，慢的那個 $\Delta \omega /2$ 是包絡線（畫面上的虛線）。當 $x_1$ 的包絡線壓到 0， 能量整包跑到了 $x_2$。耦合越弱，$\Delta \omega$ 越小，搬運越慢 — 把 k_c 調到 0.05 看看，那會慢到像在偷渡。

🎮 動手試試

1. 先按「同相」，然後把 k_c 從 0 拉到 2。ω₁ 完全沒變，圖形也沒變 — 耦合彈簧沒參與。
2. 再按「反相」，同樣拉 k_c。這次 ω₂ 一路往上，震得越來越快。同一根彈簧，兩種待遇。
3. 「只推第一顆」 — 看兩條虛線包絡線交錯呼吸，一條鼓起來另一條就癟下去。這就是能量在兩顆之間搬。
4. 把 k_c 調到 0.05，拍週期會拉到很長；調到 1.5，快到幾乎看不出「拍」，只剩雜亂。耦合強度決定資訊傳遞的速度。
5. 看「不對稱混合」：兩個模態的份量不相等，所以能量搬不乾淨 — 包絡線最低點不會碰到 0。

🚀 這個結構在哪裡還會出現

• 分子振動 — CO₂ 的對稱伸縮／反對稱伸縮／彎曲就是三個正常模態，紅外線光譜量到的是特徵值
• 建築與橋梁 — 結構的「模態分析」找的就是這些特徵向量，怕的是地震頻率剛好撞上某一個
• 電路 — 兩個耦合 LC 迴路是完全一樣的數學，能量在兩個迴路間來回就是拍
• 量子力學 — 兩能階耦合系統的 avoided crossing、以及分子軌域的成鍵／反鍵，是同一個 2×2 對角化
• PCA 與圖論 — 一樣是「找對稱矩陣的特徵向量」，只是那裡叫主成分或譜分群

📚 接下來

這頁的系統是線性的，所以有精確解，畫面上每個像素都是公式直接算出來的、沒有累積誤差。 下一頁換到另一個極端：一個沒有解析解、只能靠取樣去逼近的系統 — 而且它會在某個溫度突然變樣。