一個思想,許多戲服
到現在為止,你已經在三四套戲服裡見過同一個把戲了。普通的傅立葉級數把函數寫成正弦餘弦之和;勒讓德級數把它寫成勒讓德多項式之和;傅立葉-貝塞爾級數則用貝塞爾函數。每一次,配方都感覺一模一樣:把函數投影到一族構件上,每個構件一個積分。廣義傅立葉級數不過是把這個配方只陳述一遍,而把構件留作一個待填的空白。無論你的幾何遞給你哪一組正交本徵基,記賬方式都相同。
為什麼這件事會一再發生?因為所有這些基都源自同一處。前幾篇指南已說明:斯圖姆-劉維爾問題——帶合適邊界條件的本徵值方程 (d/dx)(p y') + q y + lambda w y = 0——是自伴的,而自伴性逼著它的本徵函數彼此正交。正弦、勒讓德多項式、貝塞爾函數並不是三個互不相干的奇蹟;它們是同一類方程在三個不同區間、帶三種不同權函數下的三個解。正交性使展開變得廉價,而正交性是有保證的。所以「在一組正交本徵基裡展開任意函數」並非僥倖的特例——它正是整套理論的要旨所在。
通用的係數公式
下面就是統御一切的那一個公式。設 {y_n} 是你的斯圖姆-劉維爾問題在 [a, b] 上的正交本徵函數,帶權函數 w(x)。你想把 f(x) 寫成 sum over n of c_n y_n(x)。要釘住某個特定係數 c_m,把兩邊同乘 y_m(x) w(x) 再從 a 到 b 積分。在右邊,正交性讓除 n = m 那一項外的每一項都消失——恰如把向量與 x 軸作點積會抹掉 y、z 分量。存活下來的是 c_m 乘以 y_m^2 w dx 的積分。解出便得通用法則:c_m =(integral from a to b of f y_m w dx)除以(integral from a to b of y_m^2 w dx)。
看清兩個部件。分子 integral of f y_m w dx 是投影——「f 中含有多少模式 m」。分母 integral of y_m^2 w dx 是 y_m 長度的平方,只起歸一化作用;若你預先把每個 y_m 標度成單位帶權範數,分母就成 1,c_m 直接是投影。決定性的角色是權 w,而它恰是初學者最愛漏掉的那一個。它不是可以無視的裝飾:它正是自伴形式中坐在 lambda y 前面的那個係數,而本徵函數是關於 w 正交的,從來不是關於裸的 dx。一旦漏掉 w,你的係數就會出錯,級數也重建不出 f。
Generalized Fourier series: f(x) = sum_n c_n y_n(x)
Project onto y_m, integrate with weight w over [a,b]:
integral_a^b f y_m w dx = sum_n c_n * integral_a^b y_n y_m w dx
|
orthogonality kills every n != m
v
integral_a^b f y_m w dx = c_m * integral_a^b y_m^2 w dx
==> c_m = ( integral_a^b f y_m w dx ) / ( integral_a^b y_m^2 w dx )
Special cases (just change a,b,w):
sines on [0,L], w=1 -> ordinary Fourier coefficient
Legendre on [-1,1], w=1 -> c_n = (2n+1)/2 * integral f P_n dx
Bessel on [0,a], w=r -> Fourier-Bessel coefficient一幅做出來的圖:球面上的一個台階
想像一個球,上半球被塗成熱的、下半球被塗成冷的——赤道以上溫度為 +1,以下為 -1。在變量 x = cos(theta) 下,這就是 [-1, 1] 上的階躍函數 f(x) = sign(x),而天然的基是勒讓德多項式(權 w = 1)。由對稱性,f 是奇函數,所以只有奇數次的 P_1, P_3, P_5, ... 存活;每個偶數係數在你動手算之前就恰為零,這正是「奇向量在偶座標軸上沒有分量」在函數空間裡的回聲。對第一個非零項套用那個通用公式,你得到 f 約等於 (3/2) P_1(x) 加上更小的修正,而其後每個係數都由各自獨立的積分算出。
- 辨認幾何的基與權:球面給出 [-1, 1] 上的勒讓德多項式 P_n,權 w = 1。
- 先利用對稱性:f(x) = sign(x) 是奇函數,所以所有偶數次係數無需任何積分就為零。
- 套用通用公式 c_n = ((2n+1)/2) 乘以 從 -1 到 1 對 f(x) P_n(x) dx 的積分,每個存活的 n 算一個積分。
- 把投影求和:f(x) = c_1 P_1 + c_3 P_3 + c_5 P_5 + ...,即冷熱球的勒讓德級數。
現在留意它與冪級數有何不同、又有何相同。泰勒級數也把 f 寫成無窮和,但它的係數通過某一點處的各階導數彼此耦合,而像 sign(x) 這樣的單個跳躍會瞬間把它摧毀——你在 x = 0 處連導數都求不出。廣義傅立葉級數沒有這種脆弱:每個係數都是 f 對 y_n 的一個全局平均,由一個貨真價實的定積分算出,它根本不在乎 f 有折角還是有跳躍。這種穩健——欣然吞下不連續的數據——正是這些展開(而非泰勒級數)成為處理帶粗糙初值的邊值問題之正確工具的原因。
級數在何種意義下等於 f?
下面這個問題,把仔細的讀者和馬虎的讀者分開。當我們寫 f(x) = sum of c_n y_n(x) 時,那個等號許諾了什麼?誠實的回答是均方意義下收斂(也叫均方收斂或 L^2 收斂):隨著你保留更多項,誤差平方的帶權積分 integral of (f - 部分和)^2 w dx 一路奔向零。這是一個關於累計總誤差的陳述,而不是關於任何單獨一點的陳述。級數可以在零散的幾個點上有一丁點偏差卻依然在均方意義下完美收斂,因為孤立的幾處誤差對積分毫無貢獻。
級數在均方意義下收斂到 f、不留任何殘餘,恰恰就是本徵基的完備性——對正則斯圖姆-劉維爾問題而言,那條深刻的保證:本徵函數觸及函數空間中的每一個方向、一個都不漏。完備性是比正交性更強、更難的定理,但對正則問題它成立。要把分工看清楚:正交性使係數容易算出;完備性使所得級數真正等於 f。兩者你都需要,而它們是不同的事實。
帕塞瓦爾:能量的賬本
每個正交展開裡都藏著一條守恆律,而它不過是被拉伸到無窮多維的勾股定理。在普通空間裡,向量長度的平方等於其各分量平方之和。帕塞瓦爾恆等式對函數說的是同一句話:f 帶權範數的平方等於其係數平方之和(每個再按其本徵函數的範數平方加權)。用符號寫,integral from a to b of f^2 w dx 等於 sum over n of c_n^2 乘以 integral of y_n^2 w dx——而若你已把 y_n 歸一化為單位範數,它就乾淨地寫作「f 的能量等於 c_n^2 之和」。函數的整個大小被逐模式精確核算,無所遺失,亦無所杜撰。
這不只是優雅——它是你的誤差表。只保留前 N 項,帕塞瓦爾告訴你殘餘的均方誤差恰是從 N+1 起的 c_n^2 尾和。於是你能預先決定為達到某容差需要多少模式:盯著尾和縮小,足夠小時就停。帕塞瓦爾還是完備性的精確憑證。其單邊版本貝塞爾不等式(c_n^2 之和至多等於範數)對任何正交集都成立;當 f 沒有任何部分逃出基底時,它恰好升級為帕塞瓦爾等式。所以若帕塞瓦爾成為等式,你的基底就是完備的;若它有所虧欠,你便漏了一些方向。
它還附贈一份令人愉快的紅利。因為帕塞瓦爾是一個精確的等式,你可以把它倒過來用,去求那些表面上與本徵函數毫不相干的和。把 f(x) = x 在 [0, pi] 上展為正弦級數,把係數餵進帕塞瓦爾的傅立葉版本,那條著名的恆等式 sum of 1/n^2 = pi^2/6 便掉了出來——一個關於整數的事實,竟由一個關於振動弦的事實變出。這與本徵函數情形下的同一恆等式,即本徵函數的帕塞瓦爾恆等式,應用於所有斯圖姆-劉維爾基中最簡單的那一組,是同一回事。
為什麼這是整卷書的引擎
退一步,看清下游的每個方法為何總要去抓正交的碎片。當你用分離變量法求解偏微分方程時,空間部分坍縮成一個斯圖姆-劉維爾問題,它的本徵基便白白送到你手上;於是本徵函數展開法說:把你的初始數據展為廣義傅立葉級數,又因為每個本徵函數都是一個按自己簡單時間表演化的純模式,難解的偏微分方程便裂成一摞各自易解、再加回去的單模式問題。熱量在棒中擴散、帶電球周圍的勢、鼓面振響——全是同一招。
兩條誠實的提醒,好讓你懷著正確的期待。其一,這一切都立在完備性之上,而完備性對有限區間上的正則斯圖姆-劉維爾問題有保證,在別處卻可能失效;在無界域上,離散譜可能鋪展成連續譜,求和變為積分,你就跨進了變換方法(下一階梯)。其二,均方收斂確實弱於逐點收斂:級數在能量意義下是正確答案,跳躍處的吉布斯過衝永不消失,而對級數逐項求導後的收斂可能遠比級數本身慢得多——所以對廣義傅立葉級數求導務必小心。這些都不是缺陷;它是這套方法所給出的、精確而誠實的契約。