當初等技巧用盡時
至今你所信賴的每一種方法,都悄悄假設了係數的某種特殊性質。對於常係數方程 a y'' + b y' + c y = 0,你猜 y = e^(rx),微積分便坍縮成一條特徵多項式。對於柯西—歐拉方程,你猜 y = x^r,再一次走運。但這些都是*特殊*係數所賜的禮物。一旦你面對一條老老實實的變係數方程,如 y'' - 2x y' + 2 y = 0,其中係數 2x 確實隨 x 而動,每一個巧妙的猜測都失效了:沒有單一的指數、沒有單一的冪、也沒有任何初等函數的有限組合能解它。
這並非罕見的窘境;它是事物的*常態*。誠實地說:大多數微分方程都沒有可用封閉形式表達的解,而變係數正是那堵牆出現之處。穿牆而過的辦法,是改變我們的雄心。我們不再強求一個完成的公式,而是接受解以無窮冪級數的形式呈現——一個和式 y = a0 + a1 (x - x0) + a2 (x - x0)^2 + ...——於是我們的工作變成逐一找出係數 a0、a1、a2、……。畢竟,一條冪級數不過是一個永不停止的多項式;只要我們能釘住它的每一個係數,我們就釘住了那個函數。
什麼使一個點成為常點
這個方法並非處處可用,所以第一件事是找一個安全的地方來展開。把二階線性方程化為標準形 y'' + p(x) y' + q(x) y = 0,首項係數除掉成 1。當 p(x) 與 q(x) 在 x0 處都*解析*時——意即各自在 x0 附近都有自己收斂的冪級數,無爆裂、無除以零、無撕裂——點 x0 便是一個常點。實務上,若 p 與 q 是多項式,或是分母在 x0 處不為零的多項式之比,則 x0 為常點。其餘各處——但凡某個係數爆炸之處——皆為奇點,本篇那套乾淨的理論觸及不到那裡。
繫於常點的那個承諾,既精確又令人安心。在常點處,冪級數解總是存在,而且你能找到兩個獨立的解——憑你早已信賴的疊加,這足以構成二階方程的完整通解。這個方法不會悄然失敗,也不會吐出一條發散的胡言級數;常點條件,正是「級數會乖乖收斂」的那紙保證書。正是這份保證,使此處成為合適的起點;至於在奇點*處*展開那更棘手的情形,則是弗羅貝尼烏斯法,本節的下一個大步。
這部機器:代入、平移、對齊
這就是整套冪級數法的引擎。寫 y = 對 n 求和的 a_n (x - x0)^n。冪級數的求導是逐項進行的,因此 y' 與 y'' 本身也是冪級數,只是指標經過平移、前面多了額外的因子。接著你把這三條級數全部代入微分方程。把這團亂麻化為方法的訣竅,是重新編號指標:你在每條級數中平移求和指標,使每一項都帶有相同冪次的 (x - x0)。一旦對齊,整條方程便成為單獨一條等於零的冪級數。
而現在是決定性的事實:一條冪級數在某區間中對所有 x 都等於零,唯有當它的每一個係數都等於零時方能成立。(這就是冪級數的恆等定理——正是這份唯一性,讓你能比對方程兩邊的係數。)於是把每個冪次 (x - x0)^k 的係數設為零,便交給你一條關聯諸 a 的方程。將所有 k 蒐羅在一起,這些方程構成一條遞迴關係:一條由先前的係數算出每個新係數的規則。微分方程,就此被點化成一道你能徒手一項一項搖出來的代數題。
Equation: y' = y (use it as a tiny demo; ordinary point x0 = 0) Try: y = a0 + a1 x + a2 x^2 + a3 x^3 + ... Then: y' = a1 + 2 a2 x + 3 a3 x^2 + ... Match powers of x (coefficient of x^k on each side must agree): x^0 : a1 = a0 x^1 : 2 a2 = a1 -> a2 = a1 / 2 = a0 / 2! x^2 : 3 a3 = a2 -> a3 = a2 / 3 = a0 / 3! general recurrence : (n+1) a_(n+1) = a_n So a_n = a0 / n! and y = a0 ( 1 + x + x^2/2! + ... ) = a0 e^x. The series rebuilt e^x from scratch -- no exponential was ever assumed.
由兩個自由常數得出兩個解
在那個小演示裡,唯有 a0 是自由的,因為 y' = y 是一階的。對一條真正的二階方程而言,遞迴關係會留下兩個未定的係數——幾乎總是 a0 與 a1——而其後的每一個係數都由它們算出。這絕非偶然:a0 恰是 y(x0),a1 恰是 y'(x0),正是二階問題一向需要的那兩份初值資料。這兩個自由常數,就是通解中那兩個任意常數的冪級數面貌。
要萃取出兩塊獨立的構件,就用乾淨的取值把遞迴跑兩遍。設 a0 = 1、a1 = 0,磨出一條級數——稱它 y1;再設 a0 = 0、a1 = 1,磨出另一條——y2。由於它們的起始資料不同,y1 與 y2 自動獨立,於是憑疊加原理,通解為 y = a0 y1 + a1 y2,此處 a0、a1 如今扮演自由常數的角色。你不曾以封閉形式寫下它們,卻已構築出一組基本解組。
這條級數可信的範圍
一條級數解的好壞,僅止於它收斂的那塊區域,所以最後一個誠實的問題是:離 x0 多遠之內你還能信它?答案優美而幾何。每個解的級數,至少在一個延伸到方程*最近奇點*的圓內收斂——而收斂半徑至少是那段距離。你甚至不必把級數算出來,就能知道它能載你走多遠:只消找出係數 p(x) 與 q(x) 在何處失常,最近的那個麻煩點便為你那個有保證的半徑封頂。
一個微妙之處,使這件事不只是有用,而是優美:最近的奇點可能是一個*複數*,完全落在實軸之外,它卻照樣決定半徑。經典的例子是 y'' + y = 0,其解為 cos x 與 sin x。它們在 0 處的冪級數對每一個實數 x 都收斂——半徑無窮——而確實,p = 0、q = 1 在整個複平面上都沒有奇點。相對地,一條唯一的麻煩坐落在 x = i 與 x = -i 的方程,其以實數為中心的級數只在 |x| 小於 1 時有保證,因為那兩個虛數點離原點的距離為 1,縱然它們從不出現在你正操作的實軸上。
用一行話握住本篇的全程。在常點處,你用「對公式的搜尋」換來「對係數的耐心培育」:代入一條冪級數,對齊冪次得出一條遞迴關係,由兩個自由常數把它跑出兩個獨立解,並信任這結果一直延伸到最近的奇點。下一篇將磨利這門手藝的後半段——如何讀出一條遞迴的規律、並精確釘住收斂半徑——然後我們才敢在奇點*處*展開,去會一會那些初等方法永遠觸及不到的貝索函數與特殊函數。