傅立葉級數翻不過的那堵牆
前幾篇指南交給你一台強大的機器:任何表現良好、以 2L 為週期重複的函數,都是純波的加權和,是它那一組諧波的離散頻譜。但請注意那行小字——重複。傅立葉級數只由能在週期裡整數次塞滿的波搭成,所以它造出來的東西永遠只能是週期的。給它一個單脈衝、一個孤零零的高斯鼓包、一口敲響一次便歸於沉寂的鐘,級數會乖乖地在一個寬度為 2L 的窗口裡復現那個形狀——然後把它向左向右永遠複製貼上下去,變成你從未要求過的一串幽靈回聲。對一個一次性的訊號來說,除了原來那個窗口,它處處都是錯的。
然而我們真正關心的大多數訊號並不是週期的。一句說出口的話、一記砸在梁上的錘擊、一個光子的波包、一束短暫雷射脈衝帶來的溫度尖峰——每一個都只發生一次,便了無蹤跡。我們需要一個辦法,去求一個活在整條數軸上、永不重複的訊號的頻譜。訣竅不是發明一台新機器,而是對我們手上這台取極限:把週期調大、再大、還要更大,直到 2L 變成無窮,那一個窗口吞下整條數軸。看著諧波在這個過程中遭遇了什麼,就是本篇的全部故事。
從最乾淨的形式出發:複指數級數
如果我們從複指數傅立葉級數出發,而不是從那個分兩條公式的正弦餘弦版本出發,取極限就遠沒那麼痛苦。藉助歐拉恆等式 e^{i theta} = cos(theta) + i sin(theta),一個單獨的複波 e^{i n pi x / L} 就把一個餘弦和一個正弦打包在了一起,於是整個級數坍縮成一個乾淨的求和:f(x) = 對所有整數 n 求和的 c_n e^{i n pi x / L},其中那條單一的係數公式 c_n = (1/2L) 乘以 從 -L 到 L 的 f(x) e^{-i n pi x / L} dx 之積分,取代了兩條歐拉-傅立葉公式。指數上那個負號,是投影在「讀出」f 中含有多少波 e^{+i n pi x / L}——還是從前那個正交性動作,只是現在一筆到位。
在我們拉伸 L 之前,有一個關鍵的記帳動作。每個諧波 n 對應一個角頻率 omega_n = n pi / L。當 n 走過整數時,這些頻率坐落在一架間距均勻的梯子的橫檔上,相鄰兩檔之間的間隔是 delta-omega = pi / L。這個間距就是離散頻譜那把「梳子的齒寬」。盯住它:當我們把 L 調大,橫檔就彼此滑近。那道由許可頻率組成的柵欄,眼看就要密到變成一堵實心牆。
讓週期趨於無窮:求和化為積分
現在動手。把係數公式代回級數,並把各項重新分組,讓間距 delta-omega = pi / L 顯式地露出來。當 L 趨向無窮,兩件事同時發生。其一,頻率間距 delta-omega 收縮趨零,於是許可頻率不再是離散的一串 n pi / L,而填滿成一個連續變量 omega,它可以取任何值。其二——這才是要害——一個對 n 求和、每一項又乘以一個小間距 delta-omega 的式子,正是黎曼和的標準排場:在極限裡,「對 n 求和再乘以 delta-omega」化成了「對 d omega 積分」。離散的求和融化成了一個遍歷所有頻率的積分。
Complex series, period 2L, with omega_n = n pi / L and delta-omega = pi / L:
c_n = (1/2L) integral_{-L}^{L} f(x) e^{-i omega_n x} dx
f(x) = sum_n c_n e^{i omega_n x}
= sum_n [ (1/2L) integral_{-L}^{L} f(t) e^{-i omega_n t} dt ] e^{i omega_n x}
= (1/2pi) sum_n [ integral_{-L}^{L} f(t) e^{-i omega_n t} dt ] e^{i omega_n x} * delta-omega
\______ a Riemann sum in omega ______/
L -> infinity (so delta-omega -> 0, omega_n -> continuous omega):
define F(omega) = integral_{-infinity}^{infinity} f(x) e^{-i omega x} dx <-- FORWARD transform
then f(x) = (1/2pi) integral_{-infinity}^{infinity} F(omega) e^{i omega x} d omega <-- INVERSE變換對:一個函數與它的頻譜
在極限中倖存下來的,是一對配套的積分——傅立葉變換對。正變換 F(omega) = 從 -infinity 到 infinity 的 f(x) e^{-i omega x} dx 之積分,取活在空間或時間上的訊號 f,返回活在頻率上的 F:它回答的是「頻率 omega 在 f 裡藏了多少?」逆變換 f(x) = (1/2pi) 乘以 從 -infinity 到 infinity 的 F(omega) e^{i omega x} d omega 之積分,則把所有這些頻率碎片疊回去,重建出原函數。傅立葉變換與它的逆是同一個對象的兩種視角,就像一個和弦既能寫成時間上的波形,也能寫成一串音符——在兩者之間往返不丟失任何資訊。
請精確地看清從級數跨到變換時究竟變了什麼。對一個週期訊號,頻譜是那串離散的係數 c_n——能量只許出現在特定的諧波頻率上,是一排孤立的尖峰。對一個一次性訊號,頻譜則是光滑函數 F(omega)——一道在每個頻率上都有定義的連續頻譜,齒與齒之間不再有空隙。一旦看出緣由便很直觀:一個永不重複的波,無法由一份有限的精確諧波菜單做成;它需要一整段連綿不斷的頻率連續體,每個頻率貢獻一個無窮小的薄片,才能搭出一個只升起一次、再不回頭的形狀。正是週期性把一道頻譜量子化成離散的譜線;失去週期性,譜線就模糊成了一條頻帶。
讀懂頻譜:F(omega) 到底告訴了你什麼
F(omega) 一般是複數,而這是一種優點,不是麻煩——它在每個頻率上都攜帶兩份實資訊。它的模 |F(omega)|,即幅度譜,說明頻率 omega 出現得有多強;它的輻角,即相位譜,說明那個波的波峰落在哪裡。單看幅度,告訴你的是音色——一個聲音是由什麼構成的,為什麼長笛和小提琴奏同一個音聽起來不同。相位告訴你的是對齊方式——把相位打亂,你保住了每個頻率的強度,卻會把一記清脆的咔嗒聲抹成一片噪聲的混響。重建訊號需要這兩半都在,這正是逆積分要帶上完整複數 F 的原因。
有一個定量的錨點能讓頻譜變得具體:能量在變換對的兩端守恆。帕塞瓦爾-普朗謝雷爾定理說,|f(x)|^2 dx 的積分等於 (1/2pi) 乘以 |F(omega)|^2 d omega 的積分——你在空間裡算出的總能量,等於你在頻率裡算出的總能量。於是 |F(omega)|^2 簡直就是單位頻率上的能量密度:它在一段頻帶 [a, b] 上覆蓋的面積,就是那條頻帶所攜帶的能量。這才是「一個訊號的頻譜」的誠實含義——不是一個含糊的比喻,而是一份可測量的、能量在頻率上的分布,是你在傅立葉級數的帕塞瓦爾定理裡逐諧波求和的那份能量的、精確的連續繼承者。
你現在握著什麼,以及它開啟的那一家族
- 要求一個訊號的頻譜,就算正變換 F(omega) = f(x) e^{-i omega x} dx 在整條數軸上的積分——一個把 f 投影到每個頻率 omega 的波上的反常積分。
- 讀出結果:|F(omega)| 是每個頻率出現了多少(幅度),它的輻角是相位,而 |F(omega)|^2 是單位頻率上的能量密度——這就是頻譜的真正含義。
- 要重建訊號,就跑逆變換:f(x) = (1/2pi) 乘以 F(omega) e^{i omega x} d omega 在整條數軸上的積分——把那一連續體的頻率碎片疊回成原來的形狀。
你已經跨過了從週期世界通往整條實數軸的那座橋,而那個曾經讀出單個諧波的投影思想,如今讀出的是一整段連續體。這一階梯往後的一切都建立在這對積分上。變換把求導變成乘以 i omega,於是傅立葉變換撬開熱方程與波方程的方式,正如拉普拉斯變換撬開初值常微分方程的方式——用代數代替微積分。電腦把這個積分取樣成一份有限的清單,當作離散與快速傅立葉變換來跑。而把核 e^{-i omega x} 換成另一個,你就得到這一家族的其餘成員——梅林變換、漢克爾變換、希爾伯特變換——每一個都是為它自己那種幾何配的恰當鏡片。你已經學會把一個一次性訊號當作一個連續的和弦來聽——並叫得出其中的每一個頻率。