一架命名曲線的參考書架
你早就懂得如何在理解一個函數之前先認識它。在任何人給你看證明之前很久,你就認識正弦了:一道升到 1、落到 -1、在 pi 的整數倍處穿過零點、永遠重複的波。你能想象它、勾畫它、用它推理。直到後來,泰勒級數和微分方程 y'' = -y 才解釋了它為何如此。本篇要為一架著名的曲線書架做同樣的事——正交多項式與貝塞爾函數——它們出現在高等物理與工程所到之處。我們現在把它們當作有已知圖像和幾條鮮明事實的命名對象來認識;以後的台階才會從它們的微分方程把它們推導出來。
下面這個統一的想法讓整架書架連成一體,在任何名字出現之前就值得直白地說清楚。你在線性代數裡學過用互相垂直的坐標軸方向來搭一個向量,因為互相垂直的分量彼此不干擾:每個坐標都能單獨讀出。同樣的把戲對函數也奏效,只要我們決定「垂直」對兩條曲線意味著什麼。答案是:把它們相乘、積分,當這個積分為零時就稱它們垂直——正交。這架書架上的每個家族,不過是一組多項式(對貝塞爾來說是一組振盪曲線),它們被精心挑選得彼此正交。它們就是函數空間的垂直坐標軸。
四個多項式家族,按它們的棲息地
區分這四個家族最乾淨的方式,不是看它們的公式,而是看每一個住在哪裡——它的區間和它的權。把它們想成適應了四種棲息地的四個物種。勒讓德多項式住在有限區間 [-1, 1] 上,配最樸素的權 w(x) = 1;它們是球面的天然函數,所以才在球座標下主宰物理。厄米多項式住在整條直線 (-infinity, infinity) 上,配鐘形權 e^{-x^2}——正是前幾篇裡那個高斯形狀;它們是量子諧振子的函數。參見勒讓德多項式與厄米多項式。
拉蓋爾多項式住在半直線 [0, infinity) 上,配衰減權 e^{-x};它們是氫原子的徑向函數,那裡電子的概率向外指數式地消退。切比雪夫多項式回到 [-1, 1] 上,但配權 1 over square root of (1 - x^2),這個權在兩個端點附近堆起極大的強調。這種端點強調使它們成為數值逼近的冠軍:它們把誤差均勻地鋪開,馴服了樸素高次擬合中那種狂野的振盪。參見拉蓋爾多項式與切比雪夫多項式。四種棲息地、四種權、四個家族——而棲息地預示了用途。
Family Interval Weight w(x) Famous home
---------- --------------- ------------------ -----------------------
Legendre [-1, 1] 1 sphere / spherical coords
Chebyshev [-1, 1] 1/sqrt(1 - x^2) numerical approximation
Hermite (-inf, +inf) e^{-x^2} quantum harmonic oscillator
Laguerre [0, +inf) e^{-x} hydrogen atom (radial part)
First few Legendre: P0 = 1, P1 = x, P2 = (3x^2 - 1)/2, P3 = (5x^3 - 3x)/2
First few Chebyshev: T0 = 1, T1 = x, T2 = 2x^2 - 1, T3 = 4x^3 - 3x
Beautiful Chebyshev fact: T_n(cos theta) = cos(n theta) (a disguised cosine!)把圖像想出來,它們就不再抽象。在 [-1, 1] 上,P_n 是一條起伏的曲線,歸一化得使 P_n(1) = 1,它在區間內部恰好穿過零點 n 次——P_2 在軸下方下沉一次,P_3 擺過零點三次,依此類推;下標越大,起伏越多,正像弦上更高的諧波。切比雪夫的 T_n 更加生動:因為 T_n(cos theta) = cos(n theta),它的圖像簡直就是一道被彎到 [-1, 1] 上的餘弦波,在 +1 與 -1 之間以等高的漣漪振盪。一旦你能勾畫出它們,這些多項式就和它們所推廣的三角函數一樣具體了。
兩台能生成任何家族的機器
你不必背這些多項式,也不必靠試錯去找它們。有兩台標準機器能把它們搖出來,現在把兩台都見識一遍,神秘感就蕩然無存了。第一台是一行求導公式。對勒讓德它寫作 P_n(x) = (1 over (2^n n!)) times (x^2 - 1)^n 的 n 階導數。把簡單的多項式 (x^2 - 1)^n 恰好求導 n 次,再縮放,P_n 就掉了出來——完美正交,無需運氣。這就是一條羅德里格斯公式,書架上的每個家族都有一條(各有自己的內函數與權)。參見羅德里格斯公式。
第二台機器更漂亮:一個單獨的封閉形式函數,它關於一個啞變量 t 的冪級數,把整個家族都裝在係數裡。對勒讓德,把 1 over square root of (1 - 2 x t + t^2) 展開成關於 t 的冪級數,t^n 的係數恰好就是 P_n(x)。一個整潔的公式,悄悄存下了無窮多個多項式。這就是一個生成函數,它是一舉證明家族性質的主力——把生成函數求導或相乘,整條整條的恆等式立刻掉出來。
為何正交性才是全部的要點
下面就是為整架書架正名的回報。因為家族是正交的,你可以把區間上幾乎任何函數展開成家族成員之和——一個廣義傅立葉級數——並通過單獨一個積分讀出每個係數,完全不必解代數方程組。這正是前面某個台階上普通傅立葉級數的精確回響:那裡正弦餘弦彼此正交,每個係數都來自一個積分。它之所以奏效,靠的是那幅垂直坐標軸圖景:當你向 P_k 投影時,其餘每個成員都貢獻為零,因為它與 P_k 的重疊積分等於零。
- 選出區間與權和你的問題相匹配的家族(球面 -> 勒讓德,振子 -> 厄米,半直線衰減 -> 拉蓋爾,逼近 -> 切比雪夫)。
- 要找到展開式中 P_k 的係數 c_k,就把你的函數與 P_k 帶權相乘並積分,再除以 P_k 平方的已知積分(該家族的歸一化)。一個積分,一個係數。
- 把 c_0 P_0 + c_1 P_1 + c_2 P_2 + ... 加起來,在項變小處截斷。因為各部分正交,一個截斷和就是該長度下最好的逼近——當你添加更多項時,沒有哪個係數需要重新調整。
在你太過信任它之前,有兩條誠實的提醒。第一,展開只在對函數加上條件時才好好收斂——粗略地說,函數應相對該權平方可積,並且和普通傅立葉級數一樣,一處跳躍間斷會在跳躍附近產生一個頑固的過衝,它不會隨著你添項而消失(吉布斯現象)。第二,單憑正交性並不能保證家族豐富到足以表示每一個函數;那個額外的性質叫完備性,是一條你應當知道正在被假定的獨立定理,而不是正交性免費奉送給你的東西。
貝塞爾函數:一道衰減的鼓面正弦
現在來認識這架書架上最著名的、卻不是多項式的居民。往圓形池塘裡丟一顆石子,或敲一面圓鼓面,漣漪會以一圈圈的環向外擴散。這種圓幾何的天然形狀,就是第一類貝塞爾函數,記作 J_n(x)。這幅心象精確而易記:J_n(x) 看上去像一道正弦波,永遠振盪下去,但振幅緩慢縮小,大致按 1 over square root of x 衰減。它是一道在向外行進時耗盡能量的正弦——正是水面上一圈漣漪在它那一環的高度鋪到越來越長的圓周上時所做的事。參見第一類貝塞爾函數。
幾條鮮明的事實讓 J_n 變得具體。J_0 在 x = 0 處從值 1 出發(像餘弦),而每個下標 n 至少為 1 的 J_n 都從 0 出發(像正弦);之後它們全都振盪並衰減。它們的零點——使 J_n(x) = 0 的那些值——並不像正弦的 pi 整數倍那樣等間距,但在遠處會趨於近乎等間距,而這些零點恰好就是被允許的鼓的頻率,正是圓鼓聽起來不和諧、不像弦那樣有確定音高的原因。同一個問題還有第二個解,即第二類貝塞爾函數 Y_n(x),它是那位在原點處爆向 -infinity 的搭檔;參見第二類貝塞爾函數。只有當你區域的中心被排除時你才保留 Y_n——對一整面鼓面,你把它丟掉,好讓答案在正中保持有限。
別把這層家族相似讀過頭。正交多項式是貨真價實的多項式——有限的、精確的、你能寫全的初等表達式。貝塞爾函數不是多項式,也沒有初等封閉形式;J_n(x) 由一個無窮冪級數(或等價地由一個積分)定義,是一個非初等函數,正如前幾篇裡的伽馬函數與誤差函數。統一整架書架的,不是它們的代數類型,而是它們的角色:每一個都是為某種特定幾何量身定制的一組正交積木。多項式用於平面與球面問題,貝塞爾函數用於圓形與柱形問題。
把它們全都繫在一起的東西
退後一步,這架書架便顯出一條統一的組織律,以後的台階會把它說精確。這些家族中的每一個——勒讓德、厄米、拉蓋爾、切比雪夫、貝塞爾——都是某個特定二階微分方程的解集,而這些方程並不是五個互不相干的偶然。它們全是同一個母模式(斯圖姆-劉維爾形式)的實例,這個母模式自動保證:解會帶著恰當的權正交地呈現出來。正交性不是我們手工強加的;它是從微分方程的結構裡掉出來的定理,正如振動模態的相互垂直從一個對稱矩陣裡掉出來一樣。
而且這些特定家族不斷出現、而非某個別的隨機清單,是有一個乾淨的理由的。當你在一個對稱區域裡取一個波動、熱傳導或量子問題並分離變量時,每個坐標裡剩下的常微分方程,恰恰就是這些命名方程之一:球面角給你勒讓德方程,柱面半徑給你貝塞爾方程,振子給你厄米方程。問題的幾何替你選定了家族。這就是為何同樣這五個名字會在靜電學、聲學、熱流與量子力學中反覆出現——它們是那些真實物理可解的對稱區域的天然詞彙。
所以,帶著這幅可用的圖景而非那些公式離開本篇。你現在能一眼認出五條命名曲線:四個在各自區間上像彎折諧波般起伏的多項式家族,以及像漸弱正弦般向外蕩漾的貝塞爾函數。你知道每一個都住在一處棲息地——一個區間和一個權——它使家族成員彼此正交,成為垂直坐標軸,你只需為每個係數算一個積分,就能沿著它們展開任何合理的函數。這恰好足夠你跟上後面的台階,在那裡這些同樣的對象將從它們的微分方程被推導出來,而那個母級的斯圖姆-劉維爾模式終將被徹底揭開。