把函數當向量,把積分當點積
在本階梯的第一篇裡,我們用一個看起來像魔法的動作求出了傅立葉係數:乘上某一個波、在一個週期上積分,然後眼看著其餘每一項都消失。那個動作並不是魔法——它是喬裝打扮的幾何,而本篇就把這層幾何挑明,好讓你看清為什麼這一切都立得住。要抓住的那一個想法是:定義在一段區間上的函數,行為恰恰像一個向量,只不過它有無窮多個分量。要把這句口號變成能拿來計算的東西,你需要一種給兩個函數取點積的辦法。
對空間中的普通向量 u 和 v,點積是 u.v = u_1 v_1 + u_2 v_2 + u_3 v_3——你把對應的分量兩兩配對再求和。對一段區間上的兩個函數 f 和 g,把「對應分量兩兩配對再求和」換成「把對齊的取值相乘、再在整段區間上加起來」,這正好就是一個定積分:定義內積為 <f, g> = f(x) g(x) 在一個週期上的積分。離散的求和變成連續的積分,是因為函數在每一點都有取值,而不只是在第 1、2、3 號槽位上有。有了這一個定義,幾何的全部詞彙——長度、夾角、垂直、投影——就整套地搬到了函數身上。
正交性:三角波是一組互相垂直的坐標軸
現在就為我們目錄裡的兩個波算一算這個內積。那個核心事實——三角函數系的正交性——是說,任意兩個不同的波,其內積恰好為零。在從 -L 到 L 的區間上:只要 m 不等於 n,sin(m pi x / L) sin(n pi x / L) dx 的積分就是 0;一個餘弦配上另一個不同的餘弦也是 0;而每一個正弦配上每一個餘弦,全是 0。原因是具體的,並不神秘:兩個不同頻率的波之乘積本身就是一段振盪,它待在零線之上的時間恰與待在零線之下的時間一樣多,於是在一整個週期裡,正面積與負面積相消為烏有。
內積唯一不為零的時候,是一個波遇上它自己的孿生兄弟。cos^2(n pi x / L) 在 -L 到 L 上的積分是 L,sin^2 也一樣——因為一個被平方的波從不取負值,它的面積無從相消。用向量的話來說:每個波的平方長度都是 L(它不是單位向量,只是有一個固定長度),而任意兩個不同的波彼此垂直。於是函數族 { 1, cos(pi x / L), sin(pi x / L), cos(2 pi x / L), sin(2 pi x / L), …… } 就是函數空間裡一組互相垂直的坐標軸,正像標準基向量分別指向 x、y、z 那樣——只不過這裡有無窮多根。
無窮維裡的畢達哥拉斯:帕塞瓦爾定理
正交性解鎖的回報就在這裡。對一個以垂直分量寫出的普通向量,它的平方長度等於各分量平方之和:|v|^2 = v_1^2 + v_2^2 + v_3^2。這無非就是畢達哥拉斯定理。因為三角波是一組互相垂直的坐標軸,完全相同的定理對函數也成立——而且它有個名字,叫帕塞瓦爾定理。f 的「平方長度」是 <f, f> = [f(x)]^2 在一個週期上的積分,帕塞瓦爾說它等於每個諧波各自平方貢獻之和,沒有任何交叉項倖存。正是這份垂直殺死了那些交叉項,正如 |v|^2 裡沒有 x.y 這一項才讓畢達哥拉斯定理這麼乾淨。
PARSEVAL'S THEOREM (trig Fourier series, period 2L)
(1/L) * integral_{-L}^{L} [f(x)]^2 dx
= a_0^2 / 2 + sum_{n>=1} ( a_n^2 + b_n^2 )
WHY it falls out of orthogonality:
square the series f = a_0/2 + sum [ a_n cos_n + b_n sin_n ]
integrate term by term over one period.
- cross terms cos_m cos_n (m != n), sin_m sin_n, cos sin -> integrate to 0
- only the SELF terms survive: integral cos_n^2 = L, integral sin_n^2 = L
divide by L -> each a_n, b_n contributes its own square, nothing else.
ENERGY reading: left side = average power of the signal over a period;
right side = power summed harmonic by harmonic.為什麼把左邊叫做「能量」?在物理裡,極多的量都按振幅的平方變化:振動弦的能量、電流耗散的功率、聲音的強度。[f(x)]^2 在一個週期上的時間平均,恰恰就是那個平均功率。所以帕塞瓦爾說:一個訊號的總功率等於它各個諧波所攜帶功率之和。從彎彎曲曲的時間圖景過渡到柱狀圖般的頻率圖景——也就是離散頻譜——其間沒有任何損耗,也沒有任何重複計數。正是這一句話,讓工程師可以談論「60 赫茲處坐著多少能量」,彷彿各個頻率是彼此分隔、互不相干的格子。是正交性讓這些格子彼此獨立。
一筆做出來的帳,以及白送給級數的一份禮物
我們把這筆帳落實到第一篇裡的那個方波上——那個奇波,在 (0, L) 上等於 +1、在 (-L, 0) 上等於 -1,它的係數是:n 為奇數時 b_n = 4/(n pi),其餘為零。帕塞瓦爾的左邊平凡得很:[f(x)]^2 處處就是 1(因為 f 不是 +1 就是 -1),於是它的平均就是 1。右邊是 b_n^2 對奇數 n 求和,也就是 (4/pi)^2 乘以 [ 1 + 1/9 + 1/25 + 1/49 + …… ]。令兩邊相等就逼出 1 + 1/9 + 1/25 + …… = pi^2 / 8。從一個方波的能量平衡裡,竟掉出一個著名無窮和的精確值——一個本來沒道理這麼容易求到的數。
這是帕塞瓦爾貨真價實的一項超能力,值得停下來體會:它把一個你也許永遠湊不出初等技巧去求和的難無窮級數,轉化成一個你能算的常規積分。挑一個函數、求出它的傅立葉係數、寫下能量平衡,你就把一個和式以閉形式求出來了。三角波給你 1/n^4 之和;別的選擇能夠到 1/n^6 乃至更高。這一整窩偶數次冪的和(也就是偶數 p 處那個收斂的 p 級數之值,歷史上的巴塞爾問題及其表親)都是這樣掉出來的。你並沒有發明什麼新的求和技巧;你只是換個角度重新用了那條能量恆等式。
- 挑一個其平方容易積分的週期函數 f(常常是一個分段常數的波,或像 x、|x| 這樣的簡單多項式)。
- 照普通辦法把它的傅立葉係數 a_n、b_n 算出來一次。
- 寫下帕塞瓦爾:(1/L) 乘以 f^2 的積分 = a_0^2/2 + sum (a_n^2 + b_n^2)。左邊是一個容易的單積分。
- 解出右邊那個未知的級數——它現在等於一個你透過積分算出的、以閉形式給出的數。
完備性:誠實的細則
現在輪到課本常常一筆帶過的那條警示。單憑正交性只能保證一個不等式——貝塞爾不等式——即各諧波能量之和至多等於總能量:a_0^2/2 + sum(a_n^2 + b_n^2) 小於或等於 (1/L) 乘以 f^2 的積分。從幾何上看這是顯然的:把一個向量投影到若干根垂直軸上,所捕獲的長度絕不可能多過這個向量本身的長度。帕塞瓦爾把這個「小於或等於」升級成精確的等式,而這次升級是一個獨立的、更強的論斷。它說三角函數系不只是正交的,而且是完備的:沒有任何一塊剩餘的 f 藏在你的坐標軸指不到的某個方向裡。你在所有諧波上收集到的能量,一點都沒漏。
能量和裡還埋著一條實用的教訓。因為各諧波能量加起來是一個有限的總數,所以項 a_n^2 + b_n^2 必須隨 n 增大趨於零——它們不可能永遠持續貢獻下去。f 越光滑,它們衰減得越快(連續函數的係數比帶跳躍的函數衰減得更快),於是絕大部分能量都住在頭幾個諧波裡。這就是有損壓縮背後的全部原理:留下響亮的低次諧波、丟掉微弱的高次諧波,而帕塞瓦爾預先精確地告訴你,你正在丟棄多少能量——丟棄多少保真度。同一套帳目會在非週期訊號那裡以帕塞瓦爾-普朗歇爾定理的面貌重現,那裡諧波之和變成了對一段連續頻譜的積分。
為什麼這是本階梯的承重思想
退一步,看看你手裡握著什麼。整套傅立葉方法立在一個幾何事實——正交性——之上,它讓你獨立地讀出每一個係數;以及一個帳目事實——帕塞瓦爾——它讓你在從不離開頻域的情況下追蹤能量。兩者合起來說:一個訊號和它的頻譜,是同一個對象的兩份忠實描述;換一組基,什麼也不丟。這與普通幾何裡挑選一組方便的垂直坐標軸完全是同一招,只不過「坐標軸」如今是一些波,而「空間」是無窮維的。
而且這個想法並不止步於正弦和餘弦。正交性加完備性這一格局,是本卷往後一切內容的模板:當你求解熱方程或一張振動的膜時,天然的搭建砌塊未必是三角波,而是別的正交族——勒讓德多項式、貝塞爾函數、某個邊值問題的本徵函數。它們每一個都帶著自己的帕塞瓦爾式能量恆等式,而其抽象版本——本徵函數的完備性——正是保證那些展開能夠捕獲整個函數的東西。在這裡把三角情形當作純幾何來吃透,往後每一個「廣義傅立葉級數」都不過是同一個故事,換了一套字母表來講。