純冪級數在哪裡崩潰
在本層的前兩篇指南中,你把一個二階線性方程的解,建構成一個以某點為中心、乾淨俐落的冪級數 y = sum a_n x^n,再去追逐係數的遞迴關係。整套做法靠著一個悄悄的假設:你所展開的那個點是一個常點,一個方程表現良好的地方。本篇指南要問的是,在「另一種」點上會發生什麼——以及當那恰恰是你無法迴避的有趣之處時,如何拯救級數法。
把方程寫成標準形 P(x) y'' + Q(x) y' + R(x) y = 0,並把某個點 x0 放到顯微鏡下。麻煩在最高階係數 P(x0) = 0 的那一刻開始。為了化成正規形 y'' + p(x) y' + q(x) y = 0 而通除,會逼出一個除以零,於是 p(x) = Q/P 與 q(x) = R/P 可能在 x0 處爆掉。這樣的點叫做奇異點:方程正規化後的係數在那裡不是解析的。在奇異點上,那條保證冪級數解存在的舒服定理,根本就不適用。
兩種奇異:正則與非正則
並非所有奇異點都同樣兇暴。關鍵的問題是,當你趨近 x0 時,p 與 q 爆得有「多嚴重」。若它們只有「輕微」的奇異——p 至多是單極點,q 至多是雙極點——那場爆炸便溫和到可以馴服。精確的判準,是把作怪的係數乘上恰好足夠的 (x - x0) 來抵消爆破:構造 (x - x0) p(x) 與 (x - x0)^2 q(x),再問這「兩者」是否都在 x0 處解析(表現良好、可展成泰勒級數)。
若兩個乘積都通過該檢驗,x0 便是一個正則奇異點——的確奇異,卻是溫馴的,本篇指南的方法在那裡行得通。若任一乘積仍然失控,x0 就是一個非正則奇異點,一個真正狂野的所在,連弗羅貝尼烏斯都束手投降,你需要更重型的機具(例如漸近展開)。整條分類規則濃縮在一張短表裡,而學會在十秒內套用它,正是本篇指南的實用技藝。
classify x0 in P(x) y'' + Q(x) y' + R(x) y = 0
p(x) = Q/P, q(x) = R/P
P(x0) != 0 ........................ ORDINARY point -> plain power series
P(x0) = 0 and both
(x - x0)*p(x) analytic at x0
(x - x0)^2*q(x) analytic at x0 .. REGULAR singular -> Frobenius works
otherwise ......................... IRREGULAR singular -> Frobenius fails
example: x^2 y'' + x y' + (x^2 - v^2) y = 0 (Bessel, at x0 = 0)
p = 1/x -> x*p = 1 analytic OK
q = (x^2-v^2)/x^2 -> x^2*q = x^2 - v^2 analytic OK
=> x = 0 is a REGULAR singular point為何純級數行不通——以及那唯一的修補
要體會普通配方為何失靈,試試你能想到最簡單的奇異方程:柯西-歐拉方程 x^2 y'' - 2 y = 0。它的解是 y = x^2 與 y = x^(-1)。看看第二個。純冪級數 y = a_0 + a_1 x + a_2 x^2 + ... 只由非負整數次方搭成,所以它「永遠」無法重現 x^(-1)——答案根本就住在這個方法所搜尋的空間之外。硬把遞迴推下去只會塌成胡言,因為最高階係數的記帳在第一步就除以了零。
弗羅貝尼烏斯的想法簡單到令人卸下心防:若答案想要像 x^(-1) 或 x^(1/2) 這樣的次方,那就「讓」級數從一個我們不預先固定的次方開始。把一個普通冪級數乘上一個未知因子 x^r,其中指數 r 是個待定的數。試解便成為 y = x^r * sum a_n x^n = sum a_n x^(n+r),其中 a_0 不為零。這個唯一的新未知量 r 就是全部的修補——它讓級數在方程要求時,恰好從一個分數、負數或其他非整數次方開始。
操作此法:最低次方決定 r
弗羅貝尼烏斯法的操作,與你已熟悉的純冪級數法如出一轍,只多了一個新的開局動作。你把 y = sum a_n x^(n+r) 及其導數代入方程,再按 x 的次方收集各項。由於現在每一項都帶著偏移 r,出現的 x 最低次方就是 x^r 本身,而它的係數必須像其他每一項一樣歸零。這個唯一的最低次方條件很特別:它只牽涉 a_0 與 r,而既然 a_0 不為零,它便化成一個只關於 r 的方程。
那個關於 r 的方程就是指標方程——一個 r 的二次式,它的兩根 r1 與 r2 是容許的起始指數,是解可以起飛的「指標」。其餘各次方的 x 給出尋常的遞迴關係,但現在遞迴帶著 r 作為參數,於是每一個根的選擇都生成自己那條係數鏈。下方的走查顯示了這項計算的骨幹;至於兩根究竟告訴了你「什麼」,那番細究正是緊接著下一篇指南的主題。
- 用 (x - x0) p 與 (x - x0)^2 q 的檢驗確認 x0 是正則奇異點;唯有如此,弗羅貝尼烏斯才獲得授權。
- 把弗羅貝尼烏斯級數 y = sum a_n x^(n+r) 及其逐項導數代入方程。
- 按 x 的次方收集;令「最低」次方(a_0 項)的係數為零,從中讀出 r 的指標方程。
- 解指標方程得其兩根 r1 與 r2(當它們為實數時取 r1 >= r2)。
- 把一個根代回一般次方的條件以得到遞迴,再逐一輾出 a_1、a_2、……來建構該根的級數解。
此法承諾了什麼,又在何處沉默
對這份保證要誠實。弗羅貝尼烏斯保證「較大」的根 r1 永遠給出一個形如 x^r1 乘上一個解析級數的真正級數解,在與先前相同那類圓盤上有效——它的收斂半徑至少延伸到方程最鄰近的「另一個」奇異點。然而較小的根 r2 是微妙的:它只在幸運的情況下才產生一個乾淨的第二解。當 r1 - r2 是非負整數時,兩根可能碰撞或互相干擾,一個對數項 log(x) 便可能被逼進第二解中。這三種結局——相異根、重根、整數差的根——正是接下來編目的三種情況。
值得停下來談談最常見的誤解。人們以為,既然我們是「以級數」展開,解就必定是個在 x0 處毫無意外的漂亮解析函數。並非如此:因子 x^r 往往使解恰在奇異點處變得不解析——x^(1/2) 在那裡有垂直切線,x^(-1) 爆掉,而對數項發散。級數捕捉的是 x0 「附近」的解,但方程的奇異性在解本身上留下了指紋。那枚指紋就是資訊:它正是鼓膜的中心、或原子的原點,在數學中宣告自己特殊地位的方式。
最後,記住整套做法的邊界。弗羅貝尼烏斯只在「正則」奇異點獲得授權;在非正則奇異點處,指標機器可能無法給出兩個可用的根,而那裡自然的展開根本不是冪式的,而是指數式與漸近式的。所以這個方法強大卻不萬能——它恰恰是數學物理那些方程的正確工具,那些方程簡直是為了擁有正則奇異點而設計的,而它正是載你前往本層最後一篇指南中貝索方程與勒讓德方程的橋樑。