跳轉到

7.1 阻尼振盪與共振#

本章導言#

\[ m\,\ddot{x} + c\,\dot{x} + k\,x = F_0 \cos\omega t. \]

這一行方程是工程的《聖經》：汽車避震器、建築抗震、收音機選台、原子鐘、雷射腔，全部活在它的參數空間裡。6.1 節已解其齊次部分（三種阻尼），4.3 節已備好特解技術；本節把兩者合併，聚焦物理最深刻的現象——共振：當強迫頻率踩中系統的固有頻率，微小的推動能累積成巨大的響應。1940 年 Tacoma 海峽大橋的倒塌使這一數學現象成為大眾記憶；本節還原其真相（那是自激，不是純共振——但共振與自激的共同源頭正是本節的方程）。

齊次解的回顧：三種阻尼的相簿#

\(m\ddot{x} + c\dot{x} + kx = 0\) ，記 \(\omega_0 = \sqrt{k/m}\) （固有頻率）、 \(2\delta = c/m\) （阻尼率）：

情形 條件 解 \(x(t)\) 相平面（6.1 節）
欠阻尼 \(\delta < \omega_0\) \(A e^{-\delta t}\cos(\omega_d t - \phi)\) ， \(\omega_d = \sqrt{\omega_0^2 - \delta^2}\) 穩定焦點（螺線）
臨界阻尼 \(\delta = \omega_0\) \((C_1 + C_2 t) e^{-\omega_0 t}\) 重根結點
過阻尼 \(\delta > \omega_0\) \(C_1 e^{r_1 t} + C_2 e^{r_2 t}\) （兩負實根） 結點

臨界阻尼的工程崇拜：門閉器、汽車避震、磅秤指針——「最快回到靜止而不過衝」的設計目標正是 \(\delta = \omega_0\) （判別式零點，4.2 節）。工程師在特徵方程的判別式上調校世界。

品質因子： \(Q = \frac{\omega_0}{2\delta}\) ——欠阻尼系統振盪多少個週期才顯著衰減的量度。鐘擺 \(Q \sim 10^3\) 、音叉更高、汽車避震 \(Q < 1\) 。一個純量概括振盪的「壽命」。

受迫振動與頻率響應#

對 \(m\ddot{x} + c\dot{x} + kx = F_0 \cos\omega t\) ，用 4.3 節待定係數法（不撞形，因 \(\omega \neq \omega_0\) 一般成立）：特解 \(x_p = A\cos(\omega t - \psi)\) ，代入解得振幅響應

\[ A(\omega) = \frac{F_0/m}{\sqrt{(\omega_0^2 - \omega^2)^2 + (2\delta\omega)^2}}, \]

相位滯後 \(\tan\psi = \frac{2\delta\omega}{\omega_0^2 - \omega^2}\) 。

讀圖（歸一化 \(A \cdot \omega_0^2 m / F_0\) 對 \(\omega/\omega_0\) ）：

  • 低頻（ \(\omega \ll \omega_0\) ）： \(A \to F_0/k\) ——跟隨準靜態變形；
  • 共振峰： \(\omega \approx \omega_r = \sqrt{\omega_0^2 - 2\delta^2}\) （阻尼小時 \(\approx \omega_0\) ），峰值 \(A_{\max} \approx \frac{F_0}{2m\delta\omega_0} = Q \cdot \frac{F_0}{k}\) ——峰值放大 \(Q\) 倍！
  • 高頻（ \(\omega \gg \omega_0\) ）： \(A \approx \frac{F_0}{m\omega^2}\) ——質量跟不上，響應按 \(\omega^{-2}\) 消失。

共振的物理：強迫每週期向系統輸入能量 \(\Delta E \propto A F_0 \sin\psi\) ；在 \(\psi = 90°\) （共振相位）時輸入功率最大、恰好與耗散平衡——振幅停在 \(Q \times\) 靜態值。無阻尼時 \(\psi\) 恆向 \(90°\) 收斂而耗散為零，能量單調累積、振幅線性增長（4.3 節的 \(t\cos\omega_0 t\) ）——阻尼不是共振的敵人，是共振的終結者。

通解的全貌：

\[ x(t) = \underbrace{e^{-\delta t}(\cdots)}_{\text{暫態：由初值決定，指數死亡}} + \underbrace{A\cos(\omega t - \psi)}_{\text{穩態：由強迫決定，永生}}. \]

「暫態死、穩態留」（4.3 節）的振動學版本：收音機開機後的雜音是暫態，選台後的音樂是穩態； \(Q\) 越高，暫態死得越慢（選台的「餘音」）。

純共振與參數共振#

純共振（無阻尼、 \(\omega = \omega_0\) ）： \(x_p = \frac{F_0}{2m\omega_0}\, t\sin\omega_0 t\) ——振幅線性增長，4.3 節撞形修正的物理顯影。

參數共振：強迫的不是力，而是參數（鞦韆擺長的週期性變化）：

\[ \ddot{x} + \omega_0^2 (1 + h\cos\gamma t)\, x = 0 \qquad \text{（Mathieu 方程，1868）}. \]

Floquet 理論（6.1 節場景三）顯示：在 \(\gamma \approx \frac{2\omega_0}{n}\) （半頻的整分數，尤其 \(n = 1\) ：強迫頻率 = 兩倍固有頻率）處發生指數放大。鞦韆的正確盪法（每個頂點站起一次）正是參數共振——能量的注入不靠外力方向，靠系統剛度的調變。橋梁風振中的某些模式、離子阱質譜儀、參量放大器（光通訊）都是 Mathieu 的用戶。

Tacoma 的真相：共振 vs 自激#

1940 年 11 月 7 日，Tacoma 海峽大橋在 42 mph 風中扭轉發散、崩塌（影片至今是工程學第一教材）。大眾敘事說「共振」；1980–1990 年代的研究（Billah–Scanlan 1991）釐清：橋的崩壞是單自由度扭轉的自激振盪（aeroelastic flutter）——氣流對橋面做功的速率隨振幅增長，等效負阻尼：

\[ \ddot{x} + \underbrace{(c_{\rm mech} - c_{\rm aero})}_{< 0}\,\dot{x} + kx = 0. \]

特徵根實部變正（6.1 節：焦點反轉）——振盪指數放大。這與共振（頻率匹配）不同：自激在任何頻率匹配下都可能，關鍵是能量收支轉正。二者統一於本節的能量視角：

\[ \frac{dE}{dt} = P_{\rm in} - P_{\rm diss} \begin{cases} < 0 & \text{衰減} \\ = 0 & \text{穩態振動} \\ > 0 & \text{發散（共振暫態 / 顫振）} \end{cases} \]

「能量收支」是振盪理論的總帳本——7.2 節的極限環（Van der Pol）就是「收支在小振幅轉正、大振幅轉負」的自限系統。

一個具體例題#

例題 7.1.1：系統 \(\ddot{x} + 0.4\dot{x} + 25.0 x = 2\cos\omega t\) （ \(m = 1\) ， \(c = 0.4\) ， \(k = 25\) ）。(a) 求 \(Q\) 與共振頻率；(b) 求共振時穩態振幅與靜態偏移之比；(c) 說明若去掉阻尼（ \(c = 0\) ）會發生什麼。

解：

(a) \(\omega_0 = \sqrt{25} = 5\) rad/s， \(\delta = \frac{0.4}{2} = 0.2\) 。 \(Q = \frac{\omega_0}{2\delta} = \frac{5}{0.4} = 12.5\) 。共振頻率 \(\omega_r = \sqrt{\omega_0^2 - 2\delta^2} = \sqrt{25 - 0.08} \approx 4.992\) ——幾乎等於 \(\omega_0\) （小阻尼）。

(b) 靜態偏移 \(\frac{F_0}{k} = \frac{2}{25} = 0.08\) 。共振振幅（近似 \(\omega = \omega_0\) ）：

\[ A_{\max} \approx \frac{F_0}{2m\delta\omega_0} = \frac{2}{2 \times 1 \times 0.2 \times 5} = 1.0. \]

比值 \(\frac{A_{\max}}{F_0/k} = \frac{1}{0.08} = 12.5 = Q\) ✓——共振放大恰為品質因子。

(c) 無阻尼時在 \(\omega = 5\) 處純共振： \(x_p = \frac{F_0}{2 m \omega_0}\, t\sin\omega_0 t = \frac{2}{10} t\sin 5t = 0.2\, t\sin 5t\) 。振幅 \(0.2t\) 線性增長——在 \(t = 50\) 秒時振幅達 10（125 倍靜態偏移）；線性模型本身早已失效（材料屈服、大變形）。數學上解存在（4.2 節：線性方程不爆破），物理上系統早已毀滅——「線性理論的無界」是「非線性災難」的預警器。

驗算： \(A(\omega)\) 公式在 \(\omega = 5\) ： \(A = \frac{2}{\sqrt{0 + (2 \times 0.2 \times 5)^2}} = \frac{2}{\sqrt{4}} = 1\) ✓。

歷史啟示#

  • ** \(Q\) 因子的文明史：教堂鐘、小提琴共鳴箱、雷射腔、原子鐘（ \(Q \sim 10^9\) ，GPS 之父）——人類用「提高 \(Q\) 」儲存振盪，用「降低 \(Q\) 」馴服振盪**；共振是把雙刃劍， \(Q\) 是握柄。
  • 模型的誠實時刻：線性理論預言無界共振，真實世界以非線性（材料、幾何）收場——線性模型的「無界」恰恰是它的價值：它告訴你系統會在什麼尺度上失效。
  • 事故改寫理論：Tacoma（1940）之後，氣動彈性力學（aeroelasticity）成為橋梁規範的必備章節；工程災難是應用數學最貴也最有效的教科書。

與其他章節及領域的關聯#

  • 4.2、6.1 節：特徵根與阻尼三態（本節的數學底座）。
  • 4.3 節：待定係數與撞形修正（共振的代數面向）。
  • 7.2 節：極限環——「負阻尼」與正阻尼並存時的穩態振盪（Van der Pol）。
  • 7.3 節：多自由度與非線性下的共振交互作用（內共振、混沌的通路）。
  • 8.3 節：脈衝響應與頻率響應是一對 Laplace/Fourier 變換； \(Q\) 對應極點靠近虛軸的程度。
  • 量子物理：共振吸收（原子躍遷）、雷射（受激輻射的 \(Q\) 儲能）；NMR（核磁共振成像）把本節方程搬進人體。
  • 音樂與聲學：共鳴腔設計、樂器的模式頻率（Bessel 模式，5.3 節的鼓）。
  • 電子學：RLC 諧振電路（1.2 節）是本節的電學雙胞胎： \(Q = \frac{\omega_0 L}{R}\) ，選台濾波的核心。

練習#

  1. 系統 \(\ddot{x} + 4\dot{x} + 16x = 0\) ：判斷阻尼類型，求 \(Q\) （欠阻尼則用 \(\omega_d\) 修正定義），並求衰減到 5% 振幅所需時間。
  2. 求 \(\ddot{x} + 2\dot{x} + 5x = 10\cos 2t\) 的穩態響應（振幅與相位），並計算靜態偏移比較之。
  3. 證明共振峰值頻率 \(\omega_r = \sqrt{\omega_0^2 - 2\delta^2}\) （對 \(A(\omega)\) 求極值），並說明 \(\delta > \frac{\omega_0}{\sqrt2}\) 時頻率響應無峰值（單調遞減）——「重阻尼系統没有共振峰」。
  4. 設計題：你要為品質 200 kg 的汽車車身設計避震器（每輪 \(k = 20000\) N/m，四輪），使車身臨界阻尼。求所需每輪阻尼係數 \(c\) ，並討論實務上為何選 \(0.3\) – \(0.5\) 倍臨界（舒適 vs 操控的折衷）。
  5. 思考題：Mathieu 方程 \(\ddot{x} + \omega_0^2(1 + h\cos\gamma t)x = 0\) 的最強參數共振在 \(\gamma = 2\omega_0\) 。請從「鞦韆」的角度解釋為什麼是兩倍頻率（提示：盪鞦韆者每個半週期起身一次；蹲下在何時做功最大？），並連結 Floquet 乘子。

進階研讀：振盪理論的五個深化場景#

場景一：多自由度與模態分析。 \(M\ddot{\mathbf{x}} + C\dot{\mathbf{x}} + K\mathbf{x} = \mathbf{F}\) 的無阻尼特徵問題 \(K\mathbf{v} = \omega^2 M\mathbf{v}\) 給固有模態（正交，5.2 節 Sturm–Liouville 的離散版）；強迫響應 = 模態分解 + 各模態單自由度響應。建築抗震（模態疊加譜分析）、橋梁監測（頻率指纹辨損傷）、手機馬達（LRA 諧振）都靠模態表。

場景二：內共振與非線性模態。非線性系統中，若兩個模態頻率成簡單比（ \(\omega_2 \approx 3\omega_1\) 等），能量在模態間週期性交換（內共振）——弦與擺的耦合實驗（1850s Melde）到太空索狀結構。非線性模態（Rosenberg 1960s、Vakakis 1990s）打破「模態獨立」的線性直覺，是奈米機械振子（NEMS）研究的日常。

場景三：同步現象。弱耦合振子系統的頻率鎖定（entrainment）：Huygens 1665 年發現掛在同一橫樑上的兩隻擺鐘反相同步（他病臥在床的觀察——科學史上最著名的臥床實驗）；Pikovsky–Rosenblum–Kuramoto 理論（Kuramoto 模型 1975）解釋螢火蟲同步閃爍、心臟起搏細胞、電網相鎖。共振的一對一擴充（振子對振子）成為複雜系統的第一通用語。

場景四：隨機振動與譜密度。實境中的「強迫」不是純弦而是隨機過程；線性系統的響應譜 = 輸入譜 × 傳遞函數模方（Wiener–Khinchin 1930s）。地震工程的反應譜、車輛平順性分析、電子噪聲（1/f 噪聲的共振放大）都用它。確定性振動學 + 隨機過程 = 隨機振動學（Crandall–Mark 1963 奠基）。

場景五：能量收集與負共振設計。振動能量收集器（piezo harvester）把環境振動轉電能——設計原則恰是「造一個小 \(Q\) 、寬頻的共振器」（與鐘錶的「高 \(Q\) 窄頻」相反）；非線性雙穩態收集器（2010s）用 6.2 節的分岔拓寬頻寬。同一份數學，既有「保持振盪」的傳統（鐘錶、原子鐘），也有「收割振盪」的新行業。

補充：Huygens 的鐘擺同步（1665）#

1656 年 Huygens 發明擺鐘（以單擺調節擒縱機構， \(T = 2\pi\sqrt{L/g}\) 由他 1659 年導出）。1665 年 2 月，他臥病觀察到並排兩鐘在半小時內達成精確反相同步，且受擾後自動回復——他把原因歸於橫樑的微幅共享振動。此「奇異的同情」（odd sympathy）現象直到 2002 年才被 Baltanás、Pikovsky 等以現代耦合振子理論完全解釋（耦合通過樑的位移而非衝量時，反相同步是唯一穩定態）。從觀察到理論，用了 337 年——同步科學的耐心。

補充：為什麼收音機能「選台」#

天線接收全部頻率的微弱電壓；LC 諧振電路的阻抗 \(|Z(\omega)|\) 在 \(\omega_0 = \frac{1}{\sqrt{LC}}\) 處峰起（本節振幅公式的電學版），帶寬 \(\Delta\omega = \frac{\omega_0}{Q}\) 。** \(Q\) 就是「選擇性」： \(Q\) 高 ⇒ 帶窄 ⇒ 相鄰電台互不干擾；但 \(Q\) 過高 ⇒ 信號的調變邊帶被削（失真）。收音機設計是 \(Q\) 的折衷藝術**——與避震器「軟硬折衷」同一數學、相反方向。

名詞速查表（本節出現的核心術語）#

  • **固有頻率 \(\omega_0\) **： \(\sqrt{k/m}\) ；系統的私語。
  • **阻尼率 \(\delta\) **： \(c/2m\) ；指數衰減率。
  • 欠阻尼 / 臨界 / 過阻尼：判別式三態。
  • **品質因子 \(Q\) **： \(\frac{\omega_0}{2\delta}\) ；振盪壽命與共振放大倍率。
  • **頻率響應 \(A(\omega)\) **：穩態振幅對強迫頻率。
  • **共振峰 / 共振頻率 \(\omega_r\) **： \(A\) 的極大值。
  • 暫態 / 穩態：齊次解（死亡）與特解（永生）。
  • 參數共振 / Mathieu 方程：剛度週期調變引起的放大；半頻規則。
  • 顫振（flutter）：氣動負阻尼自激；Tacoma 事故的真因。
  • 模態分析：多自由度系統的固有模態分解。

年表速查（本節涉及的關鍵年份）#

  • 1656：Huygens 擺鐘。
  • 1659：Huygens 導出 \(T = 2\pi\sqrt{L/g}\) 。
  • 1665：Huygens 觀察鐘擺同步。
  • 1739：Euler 非齊次解法（共振的代數誕生）。
  • 1868：Mathieu 方程。
  • 1877：Rayleigh《聲學理論》，受迫振動系統化。
  • 1932：Krylov–Bogoliubov 非線性振動方法。
  • 1940：Tacoma 海峽大橋倒塌。
  • 1948：頻率響應法（Bode）標準化於控制工程。
  • 1991：Billah–Scanlan 釐清 Tacoma 為顫振非共振。

深化問答三則#

問：共振一定是壞事嗎？ 答：共振是放大器，好壞看用途。壞：結構疲勞、噪音、Tokyo 的行人橋鎖步（2000 年倫敦千禧橋，側向鎖步同步——場景三）；好：收音機選台、MRI 成像、雷射、NMR 化學分析、石英錶。工程的第一課不是「避免共振」，而是「知道你的 \(Q\) 在跟誰共振」。

問：為什麼「共振振幅 \(\sim Q \times\) 靜態」而「無阻尼時 \(\sim t\) 」？ 答：無阻尼時能量只進不出， \(E \propto\) 功率 × \(t\) ，振幅 \(\propto \sqrt{E} \propto t\) （線性增長）；有阻尼時到達穩態「進出平衡」，振幅停在 \(P_{\rm in} = P_{\rm diss}\) 的水平，而該平衡位置由 \(Q\) 決定。穩態振幅 = 能量收支的交點——這個視角（7.2 節）正是極限環理論的開端：Van der Pol 振盪器把「收支交點」內生化，振盪不靠外力自我維持。

問：線性模型預言共振無界，是否意味線性理論「失敗」？ 答：相反，是線性理論的誠實。它精確回答了「線性世界會如何」（無界），於是精確指出非線性接手的時刻（振幅與尺寸可比時）。物理上，大振幅的彈簧硬化（Duffing）、材料屈服、幾何非線性接管，把無界化為有限響應或破壞。好模型的價值包括精確地報告自己的失效點——這是 1.2 節「線性化給趨勢，非線性定生死」的最佳註腳。

（本節完）