特別的從來不是矩形——而是座標
在這一整級裡,你在區間與矩形上分離了熱方程與拉普拉斯方程,而每一次空間因子回來時都是正弦或餘弦。人很容易就此斷定:偏微分方程根本偏愛三角函數。並非如此。正弦之所以出現,是因為邊界是一段筆直的線段、算子是樸素的 d^2/dx^2——在那樣的幾何上,正弦正是既在夾住的端點處消失、又仍是特徵函數的那些形狀。換掉房間的形狀,房間最愛的形狀也隨之改變。
大多數真實的區域並不是矩形。一張鼓面是一個圓盤;管子裡的空氣裹著一個圓柱;一顆恆星或一團電荷周圍的重力與電位住在一個球面之外。要在那裡分離變數,第一件差事就是把拉普拉斯算子改寫成貼合邊界的座標,好讓邊界成為被釘住的單一座標。在圓盤上,這意味著極座標 (r, theta),邊緣不過就是 r = a。代價是算子不再長得像 u_xx + u_yy,而長出新的、依賴 r 的項——而正是那些項,召喚出新的特殊函數。
在圓盤上:角度仍給正弦,半徑卻給出新東西
在極座標下,拉普拉斯算子變成 Laplacian u = u_rr + (1/r) u_r + (1/r^2) u_(theta theta)。代入乘積 u = R(r) T(theta) 再相除,第一篇導讀那個拱頂石論證再次點火:一個只含 r 的函數等於一個只含 theta 的函數,唯有兩者皆為常數方可。角向那塊給出 T'' + n^2 T = 0。但請留意 theta 的邊界條件不是夾子——而是你轉滿一整圈後必須回到出發點,即週期邊界條件 T(theta + 2*pi) = T(theta)。這份週期性逼使 n 為整數,且 T(theta) = cos(n theta) 或 sin(n theta)。所以角度仍說著傅立葉語——一個繞圓的傅立葉級數。
徑向方程,正是幾何咬人之處。對拉普拉斯方程,它寫作 r^2 R'' + r R' - n^2 R = 0,一個等維方程,由調和的基本積木 R = r^n 解出(以及 r^(-n),因在中心爆掉而被丟棄)。這就是為何圓盤上的穩態溫度朝中心如 r^n 般衰減。但對一面振動的鼓——也就是波動方程,時間分離後留下亥姆霍茲方程 Laplacian u + k^2 u = 0——徑向方程多出一項 k^2 r^2,變成 r^2 R'' + r R' + (k^2 r^2 - n^2) R = 0。那就是貝索方程,其乖巧的解是貝索函數 J_n(k r),根本不是正弦或餘弦。
Polar Laplacian: u_rr + (1/r) u_r + (1/r^2) u_(theta theta)
Guess u = R(r) T(theta) :
angle -> T'' + n^2 T = 0 periodic -> cos(n*theta), sin(n*theta), n = 0,1,2,...
radius -> Laplace: r^2 R'' + r R' - n^2 R = 0 -> R = r^n
Helmholtz: r^2 R'' + r R' + (k^2 r^2 - n^2) R = 0 -> R = J_n(k r) (Bessel)貝索函數是上過健身房的正弦
別讓這個陌生的名字嚇著你:貝索函數 J_n(r) 在精神上就是一個適應了向外擴散的二維世界的正弦。它一次又一次地振盪、穿越零點,正像正弦,但每一個接續的隆起都略矮一些、零點也不再等距,因為向外鋪過越來越長周長的能量必然變稀薄。把鼓的邊緣夾住,即 R(a) = 0,就要求 k a 落在 J_n 的某個零點上——而 J_n 有一道無窮的零點階梯,於是齊次邊界條件再一次挑出一串離散的容許 k 值。那些正是鼓的自然頻率,也是鼓聽起來不同於弦的緣由。
而這個妙語把整級綁在一起。正如你曾把任意的初始輪廓展開成正弦,你也能把任意的徑向輪廓展開成貝索函數——一個貝索函數展開——而係數靠的正是第三篇導讀那同一個把戲:正交性。固定 n 的那些貝索函數在圓盤上彼此正交,但有一個值得標明的轉折:這份正交性帶著一個 r 的權重,所以內積積的是 R_1(r) R_2(r) 乘以 r dr,而非僅僅 dr。那個 r 是面積元的徑向部分,也是一個喬裝的史特姆-劉維問題的標誌——下一級會把這套統一的框架攤開來講。
在球面上:緯度、經度與球諧函數
踏進三維、配上一個球形邊界,同樣的故事會多富一個維度地上演。在球座標 (r, theta, phi)——半徑、極角(緯度)、方位角(經度)——中,乘積假設變成 u = R(r) Y(theta, phi)。分離把徑向因子 R 乾淨地剝下,留下角向部分 Y 在球面上滿足它自己的特徵值問題。那個曲面問題的自然特徵函數就是球諧函數 Y_l^m(theta, phi),它們是 sin(n x) 真正的三維對應物:一個球面不可再約的振動模式。
在角向部分內部,再一次分離把經度與緯度拆開。方位角 phi 是週期的,所以它給出熟悉的 e^(i m phi),又是傅立葉。極角 theta 才是深的那個:經過代換 x = cos(theta),它的方程變成勒壤得方程,而在兩極保持有限的解是勒壤得多項式 P_l(x)(以及它們的伴隨表親)。正如正弦與貝索函數,勒壤得多項式也是正交的——這裡是在區間 [-1, 1] 上——所以球面上的一個函數展開成一個球諧函數展開,係數同樣靠正交性讀出。本篇導讀裡的每一個特殊函數,都是同一個念頭披著不同的外衣。
一個方法,多套字母表——以及它依然止步之處
退一步,看見其中的統一。配方在這一整級從未改變——變的只是房間,以及隨之而來的自然字母表。每一次的步驟都一模一樣,而每一個特殊函數,不過是所選幾何恰好偏愛的那個空間特徵函數。
- 把拉普拉斯算子改寫成貼合邊界的座標,好讓邊界成為被釘住的單一座標。
- 代入乘積假設並相除,讓方程分離成每個座標一條常微分方程,由一個分離常數綁在一起。
- 施加齊次邊界與有限性條件;它們把常數量子化,並挑出一道離散的特徵函數階梯(正弦、貝索函數、勒壤得多項式)。
- 把模態疊加,並藉正交性——配上正確的權重——定出係數,以匹配初始或邊界資料。
誠實面對這台美麗機器的邊界,正如第一篇導讀所言。它依然需要一條線性方程,因為疊加是承重的樑柱,而非線性項會讓模態相互干涉而非相加。它依然需要沿分離方向齊次的條件——一個非零的邊界值或一個源,必須先藉移至穩態減掉。它也依然需要一個可分離的幾何:沒有任何座標系能把一般團塊上的拉普拉斯算子梳理整齊,所以一張腎形的薄膜或一塊 L 形的板,產不出乾淨的乘積解。對那些區域,我們交棒給更一般的機器——用史特姆-劉維理論去認證特徵函數的完備與正交、用數值算出的特徵函數展開,以及把傅立葉級數的念頭推廣成一個以算子所提供的任意基底寫出的廣義傅立葉級數。