9.2 波方程與 d'Alembert 解#

本章導言#

\[ \frac{\partial^2 u}{\partial t^2} = c^2\, \frac{\partial^2 u}{\partial x^2}. \]

1747 年，Jean le Rond d'Alembert（達朗貝爾，1717–1783） 為振動弦寫下這個方程並解出：

\[ u(x, t) = \frac{1}{2}\big[f(x - ct) + f(x + ct)\big] + \frac{1}{2c}\int_{x - ct}^{x + ct} g(s)\, ds. \]

波分兩半、各自奔逃——這個閉式解（d'Alembert 解）是偏微分方程史上第一個通解。但它引爆了十八世紀最激烈的函數概念之爭（與 Euler、Daniel Bernoulli 的「三角級數」路線對撞，1.2 節、9.1 節的伏筆在此引爆），並留下至今活躍的問題：波與熱，兩種時間性格。本節講述 d'Alembert 解的推導、波的傳播幾何與 Huygens 原理。

推導：座標的魔法#

引入特徵座標 \(\xi = x - ct\) 、 \(\eta = x + ct\) 。鏈鎖律給

\[ u_{xx} - \frac{1}{c^2} u_{tt} = 4\, \frac{\partial^2 u}{\partial \xi\, \partial \eta}\cdot\frac{1}{?} \]

（精確計算： \(u_{xx} = u_{\xi\xi} + 2u_{\xi\eta} + u_{\eta\eta}\) ， \(u_{tt} = c^2(u_{\xi\xi} - 2u_{\xi\eta} + u_{\eta\eta})\) ，故 \(u_{tt} - c^2 u_{xx} = -4c^2 u_{\xi\eta}\) 。）方程化為

\[ \frac{\partial^2 u}{\partial \xi\, \partial \eta} = 0. \]

先對 \(\eta\) 積分： \(u_\xi = F(\xi)\) （與 \(\eta\) 無關）；再對 \(\xi\) 積分： \(u = \Phi(\xi) + \Psi(\eta)\) ——一般解：

\[ u(x, t) = \Phi(x - ct) + \Psi(x + ct). \]

兩個任意函數（二階方程的自由度，1.1 節「階數 = 常數個數」的函數版！）：

\(\Phi(x - ct)\) ：形狀不變、以速度 \(c\) 右行的行波；
\(\Psi(x + ct)\) ：左行波。

波的疊加：兩波相遇時相加、錯開後原形重現——無阻尼線性波的「長記性」（9.1 節對照）。

初值匹配： \(u(x, 0) = \Phi + \Psi = f(x)\) 、 \(u_t(x, 0) = -c\Phi' + c\Psi' = g(x)\) 。積分第二式聯立第一式解出（d'Alembert 公式）：

\[ \Phi = \frac{1}{2}\Big(f - \frac{1}{c} \int^x g\Big), \quad \Psi = \frac{1}{2}\Big(f + \frac{1}{c}\int^x g\Big) \quad\Longrightarrow\quad u = \frac{f(x-ct) + f(x+ct)}{2} + \frac{1}{2c}\int_{x-ct}^{x+ct} g. \]

解讀：初位移 \(f\) 對半分裂為左右兩波；初速度 \(g\) 貢獻區間 \([x - ct, x + ct]\) 上的積分項（全域初速度在時間 \(t\) 後的「累積援軍」）。初位移給形狀，初速度給動量。

依賴域與影響域：波的因果幾何#

依賴域：點 \((x, t)\) 的解只取決於初始線段 \([x - ct, x + ct]\) ——光錐（特徵線 \(x \pm ct =\) 常數）內的過去影響現在。
影響域：初始點 \(x_0\) 的擾動只影響 \(|x - x_0| \le ct\) 的錐形區域——擾動以速度 \(c\) 傳播、不越雷池。

對照熱方程（9.1 節熱核 \(K \sim e^{-x^2/4\alpha t}\) 處處非零）：波有有限傳播速度，熱瞬時全域滲透。愛因斯坦的因果律（無信息超光速）與波方程天然相容、與熱方程根本衝突（Cattaneo 修正方程 1948 補救）——「有限傳播速度」是相對論時代篩選物理方程的第一道關卡。

有限傳播速度 = d'Alembert 解的直接推論，也是雙曲型方程（本節）與拋物型（9.1）、橢圓型（9.3）三分天下的本質差異：

型	範例	傳播	時間箭頭
雙曲	波方程	有限速度（光錐）	可逆（記憶）
拋物	熱方程	瞬時全域	不可逆（遺忘）
橢圓	Laplace 方程	無時間（靜態）	—

端點反射與駐波#

半無限弦（ \(x \ge 0\) 、端點固定 \(u(0, t) = 0\) ）：入射波 \(\Phi(x - ct)\) 到達端點後反相反射——解法（奇延拓）：把 \(f, g\) 延拓為奇函數，d'Alembert 公式照用，端點條件自動滿足。「鏡像法」：固定端 = 掛一面「符號鏡」。

兩端固定的弦（長 \(L\) ）：無限次反射使波往復疊加，產生駐波—— \(\sin\frac{n\pi x}{L}\cos\frac{n\pi c t}{L}\) 型（9.1 節分離變數的三角親戚；事實上分離變數解與行波解等價：任意行波可展開為駐波之和——兩種解法的統一是 Fourier 1811 年工作的核心洞察）。吉他弦的每一個音，都是一次反射鏈的長期駐留。

高維：Huygens 原理#

三維波方程 \(\Box u = 0\) （ \(\Box = \partial_{tt} - c^2\Delta\) ，d'Alembert 算子——相對論的幾何語言）的基本解：球面殼 \(u = \frac{\delta(r - ct)}{4\pi c^2 t}\) ——擾動是嚴格在半徑 \(ct\) 的球面上的薄殼，殼外殼內皆平靜。Huygens 原理（1678 光的波動說、1890 年代 Volterra–Kirchhoff 嚴格化）：三維波的擾動「不留尾跡」——我們聽到的聲音是乾淨的、無殘響的（若有尾跡，聲音將糊成一片）。

二維不行：擾動產生持續的內部波（尾跡）——池塘的漣漪久久不平。Hadamard（1923）指出 Huygens 原理在偶數維失效、奇數維（ \(\ge 3\) ）成立（n 維的精細結構： \(n\) 奇時「銳利」）。我們住在三維空間，所以能聽清語言——聲學的清晰性是空間維數的禮物（8.3 節場景的延續）。

一維：擾動分裂兩半後永久駐留（ \(f(x - ct)\) 永不消失）——一維波的「完美記憶」。

一個具體例題#

例題 9.2.1：無限長弦，初速度為零、初位移為三角形

\[ f(x) = \max\big(0,\ 1 - |x|\big). \]

用 d'Alembert 公式求 \(u(x, t)\) ，描述 \(t = \frac12\) 與 \(t = 1\) （ \(c = 1\) ）時的波形，並討論長期行為。

解： \(g = 0\) ，故

\[ u(x, t) = \frac{f(x - t) + f(x + t)}{2}. \]

** \(t = 0\) **：三角形原樣（兩半重疊，各承一半，合為 \(f\) ）。

** \(t = 1\) ：兩個半高三角中心在 \(\pm 1\) ，剛好不再重疊（支撐域 \([-2, 0]\) 與 \([0, 2]\) 僅觸於原點）——兩個獨立的半高三角波**，背道而馳，原點處 \(u = 0\) 。

長期： \(t \to \infty\) 兩波各自以速度 1 奔向 \(\pm\infty\) ，形狀永不改變（無色散：所有頻率同速 \(c\) ）——非色散介質（9.1 節熱方程的 \(n^2\) 差異在此為零：所有模式同速）。∎

對照色散：深水波、光在介質中、Schrodinger 波——不同頻率不同速度，波包散開（dispersion）；波方程無此現象。色散是否存在，決定「信號能傳多遠不失真」——光纖通訊的色散補償（DSF 光纖、色散補償光纖）是 d'Alembert 理想的工程救贖。

歷史啟示#

第一個 PDE 通解：d'Alembert 1747 的解展示了「換座標使方程平凡化」的至高技巧——特徵線方法（1784 Charpit、Lagrange 的一階 PDE 理論）由此生長。
函數之爭的引爆點：d'Alembert 的「任意函數」 \(\Phi, \Psi\) （可微嗎？可表為公式嗎？）與 Bernoulli 的「三角級數」展開，把「什麼是函數」逼上梁山——爭論的最終裁決是「兩者都對」：行波與駐波是同一解的兩種表示（Fourier 級數 = 行波的頻譜分解）。
波熱之別 = 可逆與不可逆：同一根弦、同一批數學工具，卻分出「記憶」與「遺忘」兩種命運——PDE 的類型（雙曲/拋物）是物理時間之箭的數學形式。

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

9.1 節：駐波解 = 分離變數；兩法的統一。
8.3 節：基本解 \(K = \frac{\delta(r - ct)}{\text{殼}}\) 是 \(\delta\) 的 PDE 顯形。
7.3 節：d'Alembert 算子 \(\Box\) 在相對論與 Lorenz 觀念中的角色。
5.2 節：球諧函數展開 = 三維波的角向分解（球面波輻射、原子軌道）。
相對論：光錐 = 本節依賴域；Minkowski 幾何以 \(\Box\) 為度量算子（1908）；「無信息超光速」是雙曲傳播的物理化。
醫學成像：CT 與地震成像的反散射問題 = 從邊界測量反演波方程係數——波方程是「看見內部」的數學載體。
數值方法：波方程的 CFL 條件（Courant–Friedrichs–Lewy 1928）：數值波速必須 \(\ge\) 物理波速（步長比 \(\frac{c\Delta t}{\Delta x} \le 1\) ）——顯式格式被光錐管轄。

練習#

用 d'Alembert 公式解 \(u_{tt} = u_{xx}\) 、 \(u(x, 0) = e^{-x^2}\) 、 \(u_t(x, 0) = 0\) ，描述 \(t = 1\) 時的波形。
初位移為零、初速度 \(g(x) = \max(0, 1 - |x|)\) ：寫出解並討論 \(t\) 大時波形（積分項主導，波形為「坡」）。
驗證 \(\Phi(x - ct) + \Psi(x + ct)\) 滿足波方程（直接求偏導）。
半無限弦固定端反射：初位移 \(f(x) = \max(0, 1 - x)\) （支撐在 \(x > 0\) ）。用奇延拓 + d'Alembert 描述波到達端點與反射後的行為（畫示意圖）。
思考題：三維波無尾跡（Huygens）、二維有尾跡。請討論：若我們活在二維空間，聲音通訊會有什麼本質困難？這說明物理定律的哪些性質對「生命可能性」是敏感的？（可自由發揮，但需以波傳播的數學為據。）

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

場景一：色散關係與波包。平面波 \(e^{i(kx - \omega t)}\) 代入得色散關係 \(\omega = ck\) （波方程：線性、無色散）；Schrodinger 方程給 \(\omega = \frac{\hbar k^2}{2m}\) （二次、強色散——波包立即擴散）；深水波 \(\omega = \sqrt{gk}\) 。群速度 \(v_g = \frac{d\omega}{dk}\) （波包行進速度）vs 相速度 \(\frac{\omega}{k}\) ——信息以群速度傳播（Lighthill 1960s 理論化）。光纖、電漿、量子力學共用這套語言。

場景二：非線性波與孤立子。KdV 方程 \(u_t + 6u u_x + u_{xxx} = 0\) （1895 Korteweg–de Vries 為解釋 Russell 1834 年觀察的「平移不散的孤立水峰」）：非線性（陡化）與色散（展寬）精確平衡 ⇒ 孤立子（soliton）——碰撞後恢復原形（Zabusky–Kruskal 1965 計算發現，GGKM 1967 反散射法精確求解）。非線性波不是亂流，可以是精確可積的——孤立子理論是二十世紀應用數學最漂亮的篇章（2.3 節場景五的 Bäcklund 變換是其工具）。

場景三：激波與熵條件。非線性守恆律 \(u_t + f(u)_x = 0\) （如 Burgers 方程 \(u_t + u u_x = 0\) ）：特徵線交叉 ⇒ 解的梯度有限時間爆破（3.3 節場景五）⇒ 激波（shock）誕生。弱解不唯一，熵條件（Lax 1954、Oleinik 1957；Hopf–Lax 公式 1950）挑選物理解。車流的激波（堵車的移動邊界）、超音速音爆、河口的潮湧——都是這套數學。「爆破後如何續命」的問題，波理論比 ODE 更早面對。

場景四：數值雙曲——CFL 與迎風。顯式格式穩定性受光錐管轄： \(\frac{c\Delta t}{\Delta x} \le 1\) （CFL 1928）；非線性激波需要迎風格式（upwind）與人工黏性（von Neumann–Richtmyer 1950）、通量限制器（1980s TVD）——曼哈頓計畫的流體模擬催生了整個數值雙曲學科。今天的星系模擬、碰撞安全氣囊 CFD，都是 CFL 與熵條件的孫子。

場景五：廣義相對論的波。Einstein 1916 預言引力波（ \(\Box\) 算子在時空彎曲幾何上的本徵波）；百年後 LIGO 2015 年 9 月 14 日直接探測 GW150914（兩黑洞並合，13 億光年外）——時空的漣漪以光速傳播、無尾跡（四維時空的 Huygens 性質）。波方程從一根振動的弦，走過 268 年，抵達了宇宙的耳膜。d'Alembert 1747 年如果知道他的特徵座標有一天被用於測量時空本身……

補充：d'Alembert 的身世與「三角函數之父」的陰影#

d'Alembert 是私生子，被遺棄於巴黎聖讓勒龍教堂台階，取名自教堂；養父為玻璃工匠。他成為數學家、百科全書派副主編（與 Diderot 並肩），寫下「極限」條目（1.1 節）。1747 年弦振動論文的爭論（與 Euler、D. Bernoulli、Lagrange）實質是函數概念之爭——d'Alembert 要求解的初值是「可微的」，Bernoulli 的級數可以是「不連續的」，Euler 中間搖擺。這場爭論的和解需要一個新的函數概念（Dirichlet 1837），比爭論本身晚了九十年。d'Alembert 晚年說：「我的工作在於尋找可以縮短的長路。」——特徵座標正是這句話的注腳。

補充：為什麼「初速度的積分」出現在解中#

\(\frac{1}{2c}\int_{x-ct}^{x+ct} g\, ds\) 項的直觀：初速度 \(g\) 是「連續不斷的推動」在時間 \([0, t]\) 內的累積——每一時刻 \(s\) 的推動以波的形式在時刻 \(t\) 到達 \(x \pm c(t - s)\) 的位置，疊加恰為該區間積分。「因果的時空堆疊」——這個圖像經 Duhamel 原理（4.3 節場景二）推廣為：任何時空源 \(F(x, t)\) 的響應 = 源在光錐上的積分。波方程的 Green 函數是「光錐上的 \(\delta\) 」——波只記住它過去光錐上的事。

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

d'Alembert 解： \(u = \frac{f(x-ct)+f(x+ct)}{2} + \frac{1}{2c}\int g\) 。
行波： \(\Phi(x - ct)\) 、 \(\Psi(x + ct)\) ；形狀不變的平移。
特徵線： \(x \pm ct =\) 常數；波的因果軌道。
依賴域 / 影響域：光錐內的過去與未來。
有限傳播速度：擾動以 \(c\) 為上限；雙曲型的本質。
駐波：兩端反射形成的固定形狀振盪；行波的疊加等價物。
Huygens 原理：3D（奇維）波無尾跡；偶維失效。
色散關係： \(\omega(k)\) ；群速度 \(d\omega/dk\) 。
孤立子：非線性平衡色散的平移不變波（KdV）。
激波與熵條件：特徵交叉後的物理弱解選擇。
CFL 條件：數值雙曲的光錐約束。

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

1678：Huygens 波動說。
1747：d'Alembert 弦振動方程與解。
1753：D. Bernoulli 駐波疊加主張。
1784：Charpit–Lagrange 特徵線法。
1834：Russell 觀察孤立水波。
1895：KdV 方程。
1905–1908：Einstein–Minkowski 時空幾何。
1928：CFL 條件。
1950：Hopf–Lax 激波理論。
1965：Zabusky–Kruskal 孤立子；1967 GGKM。
2015：LIGO 探測引力波。

深化問答三則#

問：為什麼二階方程的「通解」是兩個任意函數，而非兩個任意常數？ 答：因為對 \(t\) 而言的初值是函數（整條初始線的位移與速度），不是兩個數——PDE 的初值是無窮維的（每一點 \(x\) 都有自己的初值）。ODE 的 \(n\) 個常數對應「單點的 \(n\) 個數」；PDE 的兩個函數對應「一條線上的兩個函數」。維數升級把常數升級為函數——這是 1.1 節「階數 = 常數個數」在無窮維的變奏。

問：行波解與駐波解（9.1 分離變數）矛盾嗎？ 答：不矛盾，是同一解空間的兩組基。任意行波 \(f(x - ct)\) 可展開為 \(\sum b_n \sin\frac{n\pi(x - ct)}{L}\cdot\frac{?}{}\) 型駐波組合（有界區域上）；任意駐波可寫成反向行波之和（ \(\sin kx\cos\omega t = \frac12[\sin(kx - \omega t) + \sin(kx + \omega t)]\) ）。公式 vs 級數之爭（1747–1822）的最終答案是「兩種基、同一空間」——這正是傅立葉分析的元命題。

問：有限傳播速度為什麼是「好」的性質？ 答：三重好處。（1）因果可分：區域外的事件不影響你（可局部建模、可並行計算）；（2）數值可解：CFL 允許顯式格式（每步只更新鄰域）；（3）物理可信：符合相對論因果律。熱方程的全域即時傳播雖然數學上方便（熱核全域光滑），物理上需要修正（Cattaneo 方程、相對論熱傳導）。好的 PDE 應該知道自己的光錐——這是現代場論的入門戒律。

（本節完）