千百次分離,一個反覆出現的形狀
上一篇導讀裡,你已看見整級的寓意:分離變數法從不只交給你一條空間常微分方程——它交給你的是一個空間特徵值問題,一條會問「對哪些特殊的數,才存在一個能符合邊界條件的非零形狀?」的微分方程。在均勻的棒上,那問題是 X'' + λ X = 0 配上 X(0) = X(L) = 0,答案是熟悉的正弦。但只要你的幾何一彎曲、或材質一變化,空間常微分方程就會長出額外的項,從一個問題到下一個問題,看起來開始真的各不相同。這篇導讀要講的,是一個安靜的奇蹟:所有這些長相互異的問題,骨子裡都是同一類問題,而有一個標準形式能把它揭露出來。
在極座標下分離拉普拉斯方程,會吐出一條徑向方程,帶著一個 1/r 項與一個 -m^2/r^2 項——這就是貝索方程。在球面上分離,會吐出勒讓德方程,二階導數上騎著一個 (1 - x^2) 因子。分離量子諧振子,會吐出厄米方程。這些看起來毫不相像。然而 史特姆-劉維問題 所揭示的觀察是:經過誠實的代數運算後,它們每一個都能寫成同一個範本——自伴形式——而那個範本,正是它們共有的良好行為所棲身之處。
自伴範本
這就是那個通用形狀。自伴史特姆-劉維形式 把空間算子寫成「導數的導數」,外面裹著三個係數函數:(p(x) y')' + q(x) y + λ w(x) y = 0,定義在區間 a < x < b 上。慢慢讀它。整個二階導數的部分,都被收進一個外層的導數 (p y')' 之內,其中 p(x) > 0 是一個光滑的係數。接著 q(x) 是一個不帶導數的單純乘子。而 λ 就是我們要獵捕的特徵值,它附著的不只是 y,而是 y 乘上一個為正的 權函數 w(x) > 0。那個被收攏的形式 (p y')',就是整個訣竅所在;一切美好都從它流淌而出。
Self-adjoint Sturm-Liouville form: d / dy \ -- | p(x) -- | + q(x) y + lambda w(x) y = 0, a < x < b dx \ dx / p(x) > 0 (smooth) ... lives inside ONE outer derivative q(x) ... plain multiplier, no derivative w(x) > 0 (weight) ... rides with the eigenvalue lambda Examples reduced to this form: plain string: p=1, q=0, w=1 -> y'' + lambda y = 0 Bessel: p=x, q=-m^2/x, w=x Legendre: p=1-x^2, q=0, w=1
何必這樣費勁把導數收攏起來?因為幾乎任何二階線性常微分方程——a(x) y'' + b(x) y' + c(x) y + λ d(x) y = 0——都能被逼進這個形式。你去找一個積分因子 μ(x),使得乘過去之後,y'' 與 y' 兩項恰好坍縮成 (p y')'。配方是 p(x) = exp(∫ b/a dx),其餘一切隨之而定。所以自伴形式並非少數幸運方程的特權;它是幾乎每一個正則二階算子都能被趕往的終點。那條樸素的弦方程 y'' + λ y = 0,不過是 p = 1、q = 0、w = 1 的平凡情形,不費吹灰之力就已端坐在這個形式裡。
為什麼自伴才是真正的重點
現在來看那個收攏形式裡藏著的對稱性。定義算子 L[y] = (p y')' + q y,於是問題寫成 L[y] = -λ w y。那個關鍵詞是 自伴(也叫對稱):對於兩個夠好的函數 u 與 v,u 乘 L[v] 的積分,等於 v 乘 L[u] 的積分。用符號寫,便是在 [a, b] 上 ∫ (u L[v] - v L[u]) dx = 0。這正是定義對稱矩陣的那條矩陣恆等式 u·(Av) = (Au)·v 在微分算子世界裡的孿生兄弟。下一篇導讀所要證明的、實特徵值與正交特徵函數的整座大廈,全都立足於這一道等式之上。
那道等式從哪裡來?來自分部積分,做兩次,對著那收攏的 (p y')'。當你把 u (p v')' 分部積分,會剝下導數,並在 x = a 與 x = b 兩端冒出一個邊界項 p(u v' - v u')。對 v (p u')' 做同樣的事再相減:內部的部分恰好因為這形式的構造方式而相消,只留下那個邊界項。所以自伴性並非自動成立——它是兩件事協同運作:收攏形式讓內部相消,而恰當的邊界條件必須讓剩下的邊界項歸零。兩半缺一不可。
讓它閉合的那些邊界條件
於是自伴性要的是一件很具體的事:讓邊界項 p(u v' - v u') 在 b 處的值減去 a 處的值,對每一對可容許的函數都化為零。最友善的消滅方式是分離型邊界條件——在每一端各自獨立的一個齊次條件。一個狄利克雷條件 y = 0、一個諾伊曼條件 y' = 0,或一個羅賓混合 α y + β y' = 0,置於每個端點上,都能辦到,因為每一個都讓那一端的括號各自歸零。這些,恰恰就是分離變數法早已要求的齊次條件,所以並沒有對你提出任何新要求——理論只是解釋了,為什麼那些一直以來都是對的條件。
第二種、較為微妙的閉合邊界項的方式,是用週期邊界條件把兩端黏合起來:y(a) = y(b) 且 y'(a) = y'(b),並要求 p(a) = p(b)。如今 b 處的括號不是各自歸零,而是靠匹配與 a 處的括號相互抵消。這是一個圓環或圓形區域的情形,也說明了為何在一整個週期上的傅立葉級數同時帶有正弦與餘弦——週期問題同樣是自伴的,只是靠另一種抵消。無論走哪條路,分離型或週期型,都會把你帶進同一個對稱性成立的安全港灣。
正則、奇異,以及隨之而來的權重
誠實地說,並非每個史特姆-劉維問題都同樣溫馴。一個 正則史特姆-劉維問題 是行為良好的情形:一個有限區間 [a, b],其上 p、q、w 全都連續,而且至關緊要的是,p(x) > 0 與 w(x) > 0 處處成立,連端點也不例外。樸素的弦方程 y'' + λ y = 0 在 [0, L] 上就是正則的,而對正則問題而言,特徵值會排成一道乾淨的無窮階梯,邁向 +∞,上方不會擁擠堆積。這是黃金標準的設定,中心定理在此以它最強、最令人安心的形式適用。
那些著名的特殊函數,住在更艱難的世界裡——奇異史特姆-劉維問題,其中正則性在某個端點崩壞。貝索方程有 p(x) = x,它在中心 r = 0 處消失;勒讓德方程有 p(x) = 1 - x^2,它在兩極 x = ±1 處都消失;有些問題則坐落在無窮區間上。在這樣的點上,通常的邊界條件被一個更溫和的要求所取代——只要求解在那裡保持有界。值得驚嘆的是,許多良好行為在這道弱化之下依然存活,這恰恰是為何貝索函數與勒讓德多項式之類的正交族仍能繼承史特姆-劉維的回報。接下來的導讀,正是建立在這之上。
這個形式還有最後一塊值得一提:權函數 w(x)。它不是裝飾。當中心定理說特徵函數是正交的,意思是 對權函數正交——內積是 ∫ y_m(x) y_n(x) w(x) dx,那個 w(x) 因子就嵌在裡頭。在彎曲或不均勻的區域上,自然的正交性是帶權的,而權重,就是你把方程梳理成自伴形式時隨 λ 一同帶來的那個東西。對貝索函數 w(x) = x;對樸素的弦 w = 1,權重隱形了,這也是為何初學者常常從未察覺它的存在。只要把形式擺對,正確的權重便自動交到你手上。
這個形式接下來所許下的承諾
退一步,看看你都鋪設了些什麼。你把分離總會咳出來的那條雜亂空間常微分方程,梳理成單一的自伴史特姆-劉維形式 (p y')' + q y + λ w y = 0,並釘下了能消滅邊界項的邊界條件。握著這一個結構性事實,史特姆-劉維定理 便蓄勢待發,要交付三項保證,而這一級接下來的篇幅都在拆解它們:特徵值是實的,並排成一道遞增的階梯;特徵函數對 w 正交;以及——最深刻的部分——它們是完備的,於是任何合理的函數都能用它們展開,化為一個廣義傅立葉級數。
這正是傅立葉級數在均勻棒上給過你的回報,只是如今從筆直的區間與常係數中解放了出來。誠實的提醒,與上一級的相同:這整套機器需要一個定義在可分離幾何上的線性、齊次問題——它不是萬用的偏微分方程求解器,而你無法分離的彎曲區域、或疊加失效的非線性項,根本就在它的射程之外。但凡它射程所及,自伴形式便是那把萬能鑰匙。認出那一個形狀,驗明那兩處分部積分的相消,理論的每一份回報便都名正言順地屬於你。