JOVANA
Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

遞迴關係與收斂半徑

遞迴是真正逐個係數建造級數解的引擎,而僅僅一條定理,就能在你動手算任何東西之前,告訴你答案能觸及多遠。學會轉動這台引擎,並學會從奇點讀出收斂半徑。

從一整面方程到一條規則

在上一篇導引中,你看過了冪級數法的整套儀式:把未知函數寫成 y = a_0 + a_1 x + a_2 x^2 + ...,逐項微分,代入方程,再依 x 的冪次收攏。回給你的不是一個整潔的答案,而是一整面條件——x^0 的係數一條方程、x^1 一條、x^2 一條,如此無止境。面對無窮多條方程,你或許會絕望。出路在於它們並非彼此獨立:它們全是同一條規則的重複。那條不斷重複的規則,就是遞迴關係

遞迴關係是一條用較早的係數表示每個新係數的公式,例如用 a_n 表示 a_(n+2)。把它想成一列骨牌:你親手推倒頭一兩張(那是你的自由起始值),其後每個係數便自動接連倒下,各由它前面的決定。方程把一張無窮長的要求清單,壓縮成了一條你可以像搖柄一樣轉動的可攜規則。在這整個主題裡,解微分方程,其實就是找到並轉動這支搖柄。

轉動遞迴

讓我們真的在一個乾淨的例子上轉動搖柄:y'' + y = 0——一個你早已知道答案(正弦與餘弦)的方程,正好適合用來觀看機器運作。代入級數、令 x^n 的係數為零,便產生遞迴 a_(n+2) = -a_n / ((n + 2)(n + 1)),對每個 n >= 0 都成立。這是一個兩項遞迴:它只把一個係數連到落後它兩步的那個,從不連到緊鄰的那個。這個相隔兩格的間隙,馬上就要做出一件漂亮的事。

  1. 選定兩個自由常數 a_0 與 a_1;其餘一切都由它們決定。
  2. 餵入偶數指標:a_2 = -a_0/(2·1) = -a_0/2!,再 a_4 = -a_2/(4·3) = a_0/4!,再 a_6 = -a_0/6!,...
  3. 另外餵入奇數指標:a_3 = -a_1/3!,再 a_5 = a_1/5!,再 a_7 = -a_1/7!,...
  4. 讀出這兩條規律:a_0 那條鏈加總成 cos(x),a_1 那條鏈加總成 sin(x)。通解即 a_0 cos(x) + a_1 sin(x)。

留意兩項遞迴免費奉送的這個美妙結構事實:由於每一步都跳兩格,偶指標的係數自成一條鏈(全由 a_0 建起),奇指標的係數則另成一條完全分離的鏈(全由 a_1 建起)。級數俐落地分裂成一個偶冪解與一個奇冪解——而這兩塊正是你的兩個獨立解。你不必去搜尋第二個解;遞迴與第一個解一口氣就把它一併交給了你。

當鏈條停下:多項式現身

正是在這裡,遞迴的價值超越了重建正弦與餘弦。在許多物理方程裡,遞迴帶著一個參數——一個烤進方程中的常數——而對該參數的某些特殊值,其中一條鏈會撞上一個被迫為零的係數,自此往下的一切也都為零。無窮級數終止成一個有限的多項式。這並非令人愉快的偶然,而是物理學中最重要的那些特殊函數背後的機制。

勒讓德方程 (1 - x^2) y'' - 2x y' + l(l + 1) y = 0,其遞迴算出來是 a_(n+2) = a_n · (n(n+1) - l(l+1)) / ((n+2)(n+1))。盯著分子看:當行進指標 n 等於參數 l 時,因子 n(n+1) - l(l+1) 恰好變為零,於是 a_(l+2) = 0,自該處起整條鏈便死去。若 l 是非負整數,兩條鏈之一便終止,那個解就是一個 l 次多項式——一個勒讓德多項式。另一條鏈則永遠跑下去,且恰好在 x = +-1 處爆掉。

答案能觸及多遠?

一個級數解只有在那塊無窮和真正收斂到有限值之處才算誠實。在那區域之外,項無界增長,「解」不過是一堆毫無意義的符號。級數的觸及範圍——從你的展開點向外、級數忠實代表一個真正解的那段距離——就是它的收斂半徑。一個常微分方程的級數並不是通行全部 x 的免費通行證;它是一個局部物件,只在一個圓盤上有效,而你得知道那個圓盤有多大。

美妙之處——也是本階段的頭條結果——在於你能在算出任何一個係數之前就預測半徑。定理說:對一個在常點 x0 附近展開的方程,其收斂半徑至少等於 x0 到方程最近奇點的距離。所以這份食譜純粹是記帳:化為標準形式,找出係數爆掉之處,量出 x0 到最近那一處的距離,那段距離就是你保證的最小觸及範圍。你早先學的奇點分類並非抽象的歸類學——它字面上把你的答案圍了起來。

  equation                          x0     nearest singular pt      guaranteed R
  --------                          --     -------------------      ------------
  y'' + y = 0                       0      none (entire)            infinity
  (1 - x^2) y'' - 2x y' + 6y = 0    0      x = 1, x = -1            1
  (1 + x^2) y'' + ... = 0           0      x = i, x = -i  (!)       1
  (x - 3) y'' + y = 0               0      x = 3                    3
直接從奇點讀出 R:量出 x0 到最近一個的距離。第三列的奇點是複數——在實線上看不見,卻仍把半徑封頂在 1。

你看不見的那個奇點

再看一次在 x0 = 0 附近展開的 (1 + x^2) y'' + ... = 0。沿實軸,係數 1 + x^2 從不為零——它是一道平滑友善的隆起,處處看不出任何麻煩。然而半徑卻只有 1,而非無窮。原因是 1 + x^2 = 0 發生在 x = i 與 x = -i,兩個高踞於複平面上的點,各距原點 1。定理是在複平面裡量距離的,所以一個躲在實線之外的奇點,仍掌控著一個看起來完全是實數的級數能觸及多遠。

這解釋了一個確實令新手困惑的現象。你可以有一個在整條實線上看起來無比平滑的實函數——沒有尖峰、沒有角點、沒有無窮——然而它的泰勒級數卻在某個有限半徑處突然停止收斂,從實數圖形上看不出任何理由。理由總是存在的,只是藏在複平面裡。它誠實地提醒我們:實線只是一個更豐富幾何的一片薄切,而收斂聽命於那個幾何,並非聽命於你的眼睛沿著 x 所能看見的。

兩點誠實的提醒為此作結。第一,定理給的是保證的最小值,而非精確值——真正的半徑可能更大(級數或許越過最近的奇點仍繼續收斂),但絕不會更小。第二,這整套乾淨的故事只活在常點處;在奇點處,平凡級數可能失效,你得動用下幾篇導引中更重型的機器。在有效區間之內,級數解是一個真實可靠的答案——只要記得它是局部的,且其觸及範圍可以事先算出。