兩個日常物件,各藏著一個方程
把一根吉他弦在兩端固定,拉到一邊,再鬆手——它便彈成一片歌唱的顫動。另外,給一根金屬棒的中段加熱,然後走開——暖意四散,熱點冷卻,給足時間,整根棒就漂向一個平坦的溫度。這看起來是兩個截然不同的故事:一個在鳴響,一個在歸於沉寂。然而它們受幾乎相同的數學支配,而傅立葉級數正是把二者一併撬開的那一件工具。本篇就是這一階梯那個抽象斷言——「任何訊號都是諧波之和」——終於在真正的物理上兌現價值的地方。
每個物件都藏著一個偏微分方程——它聯繫的不是普通導數,而是某個含兩個變量(這裡是位置 x 與時間 t)的未知函數的好幾個偏導數。振動弦服從波動方程,寫作 d^2u/dt^2 = c^2 d^2u/dx^2(常縮寫為 u_tt = c^2 u_xx),其中 u(x,t) 是弦在位置 x、時刻 t 處的橫向位移,c 是由張力與質量定下的波速。受熱細棒服從熱方程,du/dt = k d^2u/dx^2(u_t = k u_xx),此處 u(x,t) 是溫度,k 是熱擴散率。兩個物理世界,兩個偏微分方程,差別只在一個微妙處:一階時間導數對上二階時間導數。
另有兩個事實把每個問題釘死為唯一的答案。其一,邊界條件:兩端都被夾住。弦在 x = 0 與 x = L 處被繫住,故那裡不能動,u(0,t) = u(L,t) = 0;棒的兩端被維持在零溫,方程形式相同。其二,初始條件:你在 t = 0 時刻起步的形狀——弦被撥動後的輪廓,或棒的初始溫度分佈。偏微分方程加上邊界條件再加上初始狀態,構成一個邊值問題,而這整套恰好只有一個解。我們全部的活兒,就是把它找出來。
分離變量法:讓傅立葉變得必不可少的引擎
核心動作是分離變量法:賭答案能分解成「只依賴空間的部分」乘以「只依賴時間的部分」,即 u(x,t) = X(x) T(t)。這是一個猜測,而非定律——它無法表示兩個變量的每一個函數——但對這些乾淨的、線性的、帶著如此簡單邊界的問題,它恰恰是對的那個猜測。把它代入偏微分方程,一樁小小的奇蹟發生了:所有含 x 的東西堆到一邊,所有含 t 的東西堆到另一邊。
- 把 u = X(x) T(t) 代入偏微分方程。對棒而言:X T' = k X'' T。兩邊同除以 k X T,得到 T'/(k T) = X''/X。
- 論證分離常數。左邊只依賴 t,右邊只依賴 x,可它們對一切 x 與 t 都相等——所以二者必同等於某一個常數,記作 -lambda。一個含兩個變量的偏微分方程,就此裂成兩個常微分方程。
- 在夾緊端條件 X(0) = X(L) = 0 下求解空間方程 X'' + lambda X = 0。只有特殊值 lambda = (n pi / L)^2 才給出非零解,而那個解就是 X_n(x) = sin(n pi x / L)。是邊界親手挑出了這些諧波。
- 對每個 n 求解相應的時間方程,再把所有積木 u_n(x,t) = X_n(x) T_n(t) 疊加成一個無窮和,並用初始形狀釘住各係數——而這恰恰就是一個傅立葉級數。
請細看第 3 步剛發生了什麼,因為這就是整條階梯的秘密所在。我們並沒有「選擇」使用正弦——是夾緊的邊界把正弦強加給了我們。一個在 x = 0 與 x = L 處都為零的函數,必由 sin(n pi x / L) 搭成;餘弦不肯,因為 cos(0) 不為零。邊界條件替我們挑出了那個正交族,正是你早先學過其歐拉-傅立葉公式的那套正弦基。這就是傅立葉級數與這些偏微分方程密不可分的緣由:被夾緊物件的物理,字面意義上挑選了諧波,而傅立葉的配方是求出每種各需多少的唯一辦法。
弦:鳴響著的駐波
現在把弦做完。它的空間部分如前推得 X_n(x) = sin(n pi x / L)。它的時間部分來自波動方程的兩個時間導數,給出 T_n'' + (c n pi / L)^2 T_n = 0——一個簡諧振子的方程。於是 T_n(t) 是一個在時間上以頻率 omega_n = c n pi / L 搖擺的正弦餘弦組合。每個積木 u_n = sin(n pi x / L) [A_n cos(omega_n t) + B_n sin(omega_n t)] 都是一個駐波:一個在空間上固定的正弦拱形,其高度在時間裡上下脈動,從不行進,只在原地呼吸。
完整的運動把它們全部疊加:u(x,t) = sum over n of sin(n pi x / L)[A_n cos(omega_n t) + B_n sin(omega_n t)]。最慢的模式 n = 1 定下你聽到的音高——弦的基頻 omega_1 = c pi / L——而更快的模式 n = 2, 3, …… 是它的泛音,正是這些諧波讓小提琴和長笛在同一個音上有了不同的音色。要定出各個量 A_n,令 t = 0 並要求形狀匹配被撥動的輪廓 f(x):於是 f(x) = sum A_n sin(n pi x / L),這正是初始形狀的傅立葉正弦級數。這些 A_n 就是它的傅立葉係數;B_n 以同樣方式來自初始速度。
VIBRATING STRING -- length L, ends fixed, plucked into shape f(x), released from rest
PDE: u_tt = c^2 u_xx
Boundary: u(0,t) = u(L,t) = 0 (tied at both ends)
Initial: u(x,0) = f(x), u_t(x,0) = 0 (plucked, then let go)
Separate u = X(x) T(t) -> X'' + lambda X = 0 , T'' + c^2 lambda T = 0
Boundaries force X_n = sin(n pi x / L), lambda_n = (n pi / L)^2
Time part (released from rest) -> T_n = A_n cos(omega_n t), omega_n = c n pi / L
Solution (superpose the standing waves):
u(x,t) = sum_{n>=1} A_n sin(n pi x / L) cos(c n pi t / L)
Coefficients = Fourier sine series of the pluck:
A_n = (2/L) integral_0^L f(x) sin(n pi x / L) dx細棒:同樣的模式,卻是消逝而非鳴響
現在輪到棒,留意它在何處與弦分道揚鑣。空間部分一模一樣——X_n(x) = sin(n pi x / L),正是同樣的正弦模式,因為夾緊的邊界相同。唯一變的是時間方程。熱方程只有一個時間導數,而非兩個,所以分離給出 T_n' = -k (n pi / L)^2 T_n——一個一階方程,其解是一個衰減的指數,T_n(t) = exp(-k (n pi / L)^2 t)。一個導數對上兩個導數,就是音樂與沉寂之間的全部差別:弦的模式永遠振盪,棒的模式只是消融而去。
疊加起來便得到 u(x,t) = sum over n of B_n sin(n pi x / L) exp(-k (n pi / L)^2 t)。讀一讀這個指數:衰減率是 k(n pi / L)^2,它隨 n 的平方增長。所以高次諧波——那些尖銳、扭動、含精細細節的模式——死得最快,而緩慢的 n = 1 鼓包逗留最久。這正是熱為何把東西抹平:每一處帶稜帶角的特徵都是高諧波的,而高諧波幾乎瞬間蒸發,早在棒變得通體冰冷之前,就已留下一根光滑而平淡的棒。要求出 B_n,再次令 t = 0:初始溫度 f(x) 必等於 sum B_n sin(n pi x / L),所以 B_n 又一次是起始形狀的傅立葉正弦係數。
它為何管用、在哪裡失效、又通向何方
讓我們能夠把無窮多個駐波相加、卻仍然求解偏微分方程的那條邏輯,就是疊加原理:因為兩個方程都是線性的,所以任何解之和仍是一個解。我們孤立地求解單個諧波——這很容易,因為單一模式無非是一個振子或一個衰減指數——然後相信把它們全部疊加就給出真實的運動,以初始形狀的傅立葉係數作為混合權重。「先解一個模式再疊加」,正是你在線性常微分方程那裡見過的分而治之,如今被抬升到了空間與時間的函數上。傅立葉級數恰恰是負責「分」的那台機器。
退後一步,好好品味你如今能做的事。物理學兩個偉大的方程——波動與熱——都向同一套有紀律的流程屈服:分離變量,讓邊界把正弦模式遞到你手上,求解每個模式那個小小的時間方程,再讓初始形狀的傅立葉係數為求和加權。弦之所以鳴響,是因為它的模式振盪;棒之所以抹平,是因為它的模式衰減;而區分二者的那一個數,就是時間究竟出現一次還是兩次。從這裡起,道路朝兩個方向攀升。當幾何不再是一段平直的區間——一面振動的鼓膜、一根圓柱中的熱——正弦模式便讓位給別的正交族,比如你見過的、基於特殊函數的那些展開,它們全都受同一個「投影到一組基上」的思想支配。而當物體是無限長而非被夾緊時,離散的諧波模糊成一片連續區,傅立葉級數便化為傅立葉變換。你方才親眼看著一個抽象的級數,成為了物理世界的語言。