本層樓把你帶到了哪裡
到此為止,本層樓講的是一個連貫的長故事。分離變數法總會吐出一個空間特徵值問題;這個問題總能寫成自伴的史特姆-劉維形式 -(p u')' + q u = lambda w u;而一旦寫成這樣,史特姆-劉維定理就保證特徵值為實數、特徵函數關於權函數 w 正交,以及最深刻的一點——完備性,於是任何合理的函數都能展開成廣義傅立葉級數 f = sum c_n u_n。這就是整部機器。最後這一篇要問一個很自然的後續問題:當定義域不是一段整齊的區間時,特徵函數到底長什麼樣子?
在區間 [0, L] 上,配上空間問題 X'' + lambda X = 0 與兩端歸零,特徵函數就是 sin(n pi x / L)——單純的正弦,於是你得到一個普通的傅立葉正弦級數。那看起來很特別,幾乎像魔法。但誠實的真相恰恰相反:正弦與餘弦一點也不特別。它們不過是最簡單的算子(p = 1、q = 0、w = 1)在一段直線上的史特姆-劉維特徵函數。換一種幾何,就換一個算子,於是換來另一族模態——但同樣正交、同樣完備。
彎曲的幾何讓 p、q、w 現身
為什麼圓盤不肯給出正弦?因為在極座標 (r, theta) 中分離拉普拉斯算子,跟在 (x, y) 中不是同一套代數。角向部分會繞一圈,所以 theta 仍然給出普通的正弦與餘弦——角度是週期的,它的特徵函數就是那些熟面孔。真正被幾何咬住的是徑向部分。把 Laplacian u 寫成極座標形式會冒出 1/r 與 1/r^2,當你整理徑向方程時,它落成 -(r R')' + (m^2/r) R = lambda r R。盯著它看:這正是史特姆-劉維形式,其中 p(r) = r、q(r) = m^2/r,而最關鍵的是有一個權函數 w(r) = r——這正是圓盤的面積元素 r dr 喬裝而成。
這裡有個轉折,讓它們成為更狂野、更耐人尋味的近親。在中心 r = 0 處,係數 p(r) = r 為零,於是算子的最高階項退化,你無法在那裡施加普通的邊界條件。這正是奇異史特姆-劉維問題的標記,而它在你於任何彎曲座標系中分離變數時都會冒出來。補救之道很優雅:在那個壞端點處,你不要求一個值或一個導數,而只要求解保持有界(有限)。光是這一個要求,就丟掉了會發散的解、留下了物理的那個——而且了不起的是,定理的那些承諾依然成立。
貝索函數:鼓面的模態
那個徑向方程的有界解就是 貝索函數 J_m(k r),由此產生的貝索函數展開正是傅立葉級數在圓形幾何中的近親。一個貝索函數有點像衰減的正弦:它會振盪、一次又一次穿過零點,但零點之間的間距會漂移,振幅也隨著 r 增大而緩緩淡去。想像漣漪在圓形池塘上擴散的剖面——那個一座座隆起、逐漸安頓下來的形狀,本質上就是 J_0。
做量子化的是邊界,跟弦的兩端所做的一模一樣。把半徑為 a 的鼓面邊緣夾住(即 r = a 處 u = 0),你就逼出 J_m(k a) = 0——波數 k 必須落在貝索函數的某個零點上。每個貝索函數都有自己一道無窮的零點階梯,這些零點挑出了容許的徑向波數 k_j,正如 n pi / L 在弦上挑出容許的正弦頻率。把你的初始輪廓展開成一串 J_m(k_j r),就是傅立葉-貝索級數,而係數靠正交性算出——只是現在的正交是關於權函數 r 成立的,這正是幾何交給我們的面積元素。
勒讓德與球面的諧波
上升到球面,同樣的故事就在球座標 (r, theta, phi) 中上演。經度 phi 會繞一圈,所以它給出普通的 e^(i m phi)——沿著赤道的正弦與餘弦。餘緯 theta 給出一個新的奇異史特姆-劉維問題,它的有界解就是 勒讓德多項式(以及它們的親戚,連帶勒讓德函數)。勒讓德方程是 -((1 - x^2) P')' = lambda P,定義在區間 [-1, 1] 上,其中 x = cos(theta);這裡 p(x) = 1 - x^2 在兩端 x = +/-1(球的兩極)都為零——又是奇異、又是有界性,而特徵函數出來就是勒讓德多項式 P_0、P_1、P_2、...,在 [-1, 1] 上以權函數 1 正交。
把經度因子與餘緯因子縫在一起,你就得到 球諧函數 Y_l^m(theta, phi):球面天然的駐波,一層薄薄的球殼會以這些形狀來鳴響。任何定義在球面上的函數——地球表面的溫度、宇宙微波背景天空的亮度、原子軌域的形狀——都能展開成它們的疊加,就跟區間上的函數展開成傅立葉級數的方式一模一樣。這就是球諧函數展開,而勒讓德多項式不過是它 m = 0 的特例,即只隨緯度變化的那些諧波。
interval [0,L] : eigenfunctions sin(n pi x/L) weight w = 1 (flat geometry) disk : eigenfunctions J_m(k_j r) weight w = r (Bessel, k a is a zero of J_m) sphere : eigenfunctions Y_l^m(theta,phi) weight = area element (Legendre is the m=0 slice) same recipe everywhere: f = sum c_n u_n , c_n = (integral f u_n w) / (integral u_n^2 w)
統一的圖像——以及它誠實的邊緣
- 在幾何所建議的座標系中分離偏微分方程;那些繞圈的角向因子給出普通的正弦與餘弦。
- 把剩下的(徑向或餘緯)因子讀成一個史特姆-劉維問題 -(p u')' + q u = lambda w u,並從座標變換中辨認出 p、q 與權函數 w。
- 凡是 p 在某端點為零之處(圓盤的中心、球的兩極),就把邊界條件換成有界性——這就是奇異情形,它只留下物理的、有限的解。
- 讓真正的邊界把特徵值量子化(一個貝索零點、一個整數階 l),然後用來自加權正交的係數展開你的資料,c_n = (integral f u_n w) / (integral u_n^2 w)。
一旦看穿這點,整個領域就坍縮成一個想法:各種特殊函數並不是一座收藏著彼此無關的奇珍的博物館,它們是拉普拉斯算子在不同定義域上的特徵函數,而加權正交正是把它們任何一個變成可用展開的那唯一一件工具。線上的厄米函數(量子諧振子)與半線上的拉蓋爾函數(氫原子)也加入同一個家族——全都是奇異史特姆-劉維特徵函數,全都在各自的加權空間中完備。同一套機制也驅動著用於振動與聲學的亥姆霍茲方程,那裡由球貝索函數擔任徑向的角色。
有兩個誠實的告誡,要帶著它們繼續往上爬。第一,這一切都需要一個可分離的幾何,配上齊次、常係數的邊界條件——圓盤、球、矩形。把邊界扭成花生的形狀,分離就直接失敗;你只能退回變換、格林函數或數值方法,而後面幾層樓正是在建這些。第二,特徵值那道離散的階梯,本身就是有界定義域的特徵:在無窮區間上,譜可能變成部分連續,求和變成積分,廣義傅立葉級數就交棒給傅立葉變換。各種特殊函數是分離變數法輝煌的終點——也正是你能感覺到它極限、看清這門學科其餘部分為何存在的地方。