5.1 冪級數解法#

本章導言#

Liouville 1841 年的否定定理（2.3 節）宣告：積分法的武器庫有底。但數學沒有投降，而是換了彈藥——冪級數（power series）。變係數方程 \(y'' + p(x) y' + q(x) y = 0\) 的解若不是初等函數，就定義它為級數：Bessel 函數、Legendre 多項式、超幾何函數由此誕生。本節講述冪級數解法（Newton 1665 年首創、Frobenius 1873 年臻於完備）的機械程序、收斂理論與歷史。

為何級數能解方程：解析性的封閉性#

定理（存在性，解析版）：若 \(p(x), q(x)\) 在 \(x_0\) 的鄰域內解析（可展開為收斂冪級數），則初值問題的解在該鄰域內解析，且由 \(y(x_0), y'(x_0)\) 唯一確定。

這是 3.2 節 Picard–Lindelöf 的加強版（Cauchy 1820s–1840s 與 Sofia Kovalevskaya（柯瓦列夫斯卡婭，1850–1891） 1875 年的 PDE 版定理是理論高點；Kovalevskaya 因此工作與論文《論剛體繞定點旋轉》在 1884 年成為史上第一位女數學教授——斯德哥爾摩大學，此前還需特別立法讓她能授課）。

為什麼解析性帶來好運：解析函數由其在一點的全部導數決定（Taylor 級數），而方程恰好把所有高階導數遞推地算出：

\[ y'' = -p y' - q y \Rightarrow \text{求導得 } y''' = \text{(由 } p, q, y, y' \text{ 的更低階)}, \ldots \]

方程 = 導數的遞推機。冪級數解法只是把這個遞推機顯式運轉。

機械程序：三步走#

求 \(y'' - x y = 0\) （Airy 方程，1838 年 George Airy 研究光學焦散時引入）在 \(x_0 = 0\) 附近的解：

第一步：設級數。 \(y = \sum_{n=0}^{\infty} a_n x^n\) ，則 \(y' = \sum_{n \ge 1} n a_n x^{n-1}\) ， \(y'' = \sum_{n \ge 2} n(n-1) a_n x^{n-2}\) 。

第二步：代入對齊。

\[ \sum_{n \ge 2} n(n-1) a_n x^{n-2} - \sum_{n \ge 0} a_n x^{n+1} = 0. \]

換指數使兩個級數同冪：左式第一項改寫 \(\sum_{k \ge 0} (k+2)(k+1) a_{k+2} x^{k}\) ，第二項 \(\sum_{k \ge 1} a_{k-1} x^{k}\) ：

\[ 2 a_2 + \sum_{k \ge 1} \big[(k+2)(k+1) a_{k+2} + a_{k-1}\big] x^k = 0. \]

第三步：逐項清零（遞推）。 \(2 a_2 = 0 \Rightarrow a_2 = 0\) ；對 \(k \ge 1\) ：

\[ a_{k+2} = -\frac{a_{k-1}}{(k+2)(k+1)}. \]

自由參數 \(a_0, a_1\) （正是初值 \(y(0), y'(0)\) ！）。算出：

\[ a_3 = -\frac{a_0}{6}, \quad a_4 = -\frac{a_1}{12}, \quad a_6 = -\frac{a_3}{30} = \frac{a_0}{180}, \quad \ldots \]

兩個線性無關解：

\[ \begin{aligned} y_1 &= a_0\Big(1 - \frac{x^3}{6} + \frac{x^6}{180} - \cdots\Big), \\ y_2 &= a_1\Big(x - \frac{x^4}{12} + \frac{x^7}{504} - \cdots\Big). \end{aligned} \]

收斂半徑：無窮（遞推係數 \(\sim \frac{1}{n^2}\) 衰減快）。注意解的結構： \(y_1\) 含 \(x^{0,3,6,\ldots}\) 、 \(y_2\) 含 \(x^{1,4,7,\ldots}\) ——三步一跳。這是遞推「跨三步」的痕跡，也是 Airy 方程解（ \(\operatorname{Ai}(x)\) 等）的指紋。

收斂半徑：奇點決定命運#

定理（Cauchy–Hadamard）：級數解的收斂半徑 \(R\) = 從展開點到最近奇點（在複平面）的距離。

這把 3.3 節的奇點分類變成實用預言：

\(y'' + y = 0\) ：無有限奇點 \(\Rightarrow\) \(R = \infty\) （ \(\sin, \cos\) 全域解析）✓；
\(\frac{y'}{1 + x^2}\) 型： \(x = \pm i\) 是複奇點，雖在實軸上看不見，仍把 \(x_0 = 0\) 的收斂半徑限為 \(R = 1\) ！實軸上完美的係數，解級數卻在 \(|x| = 1\) 斷掉——十九世紀數學家對此困惑許久，複分析（Weierstrass、Cauchy 的複積分理論）給出最終解釋。看實軸不足以下判斷，必須望向複平面。

Frobenius 法：在奇點處繼續工作#

正則奇點 \(x_0\) 處（3.3 節： \((x - x_0) p\) 、 \((x - x_0)^2 q\) 解析），普通 Taylor 級數失效，Ferdinand Georg Frobenius（弗羅貝尼烏斯，1849–1917） 1873 年的方法：設

\[ y = x^r \sum_{n=0}^{\infty} a_n x^n \qquad (a_0 \neq 0), \]

代入後最低次冪給出指標方程（indicial equation）定 \(r\) ；其後逐階遞推。以 Bessel 方程（5.3 節主題） \(x^2 y'' + x y' + x^2 y = 0\) （ \(\nu = 0\) ）為例：指標方程 \(r^2 = 0\) 重根 \(r = 0\) ，一解為級數，另一解需對數修正 \(y_2 = y_1 \ln x + (\cdots)\) ——重根與共軛複指標、整數差指標三種情況各有劇本，5.3 節詳演。

Frobenius 法的威力：它不僅給出解，還從指標方程讀出解在奇點附近的行為（冪律 \(x^r\) 或對數）——奇點不再是障礙，而是資訊源（呼應 3.3 節）。

超幾何方程：級數解的統一帝國#

\[ x(1 - x) y'' + \big[c - (a + b + 1) x\big] y' - ab\, y = 0 \]

的解 \({}_2F_1(a, b; c; x) = \sum \frac{(a)_n (b)_n}{(c)_n n!} x^n\) （ \((a)_n = a(a+1)\cdots(a+n-1)\) ）——超幾何級數（Euler 1778、Gauss 1812 給出收斂的完整理論，Gauss 的成名作）。它的三個正則奇點 \(0, 1, \infty\) （Riemann P-方程，3.3 節）經參數與變數代換，收編了幾乎全部經典特殊函數：

\[ e^x = {}_2F_1(\cdot\,;\,x), \quad \arcsin x, \quad P_n(x)\ (\text{Legendre}), \quad J_\nu(x)\ (\text{Bessel 經代換}), \quad \ldots \]

Gauss 1812 年論文的副標題幾乎是宣言：「級數是函數的更一般形式」——函數概念從「公式」走向「級數定義」，現代特殊函數論由此奠基。

一個具體例題#

例題 5.1.1：用冪級數解 \(y'' - y = 0\) （已知真解 \(e^{\pm x}\) ），驗證級數法重現指數函數；並指出遞推結構。

解：設 \(y = \sum_{n \ge 0} a_n x^n\) 。代入： \(\sum_{n \ge 2} n(n-1) a_n x^{n-2} - \sum_{n \ge 0} a_n x^n = 0\) 。對齊（令 \(k = n - 2\) ）：

\[ \sum_{k \ge 0} \big[(k+2)(k+1) a_{k+2} - a_k\big] x^k = 0 \quad\Longrightarrow\quad a_{k+2} = \frac{a_k}{(k+2)(k+1)}. \]

偶鏈： \(a_2 = \frac{a_0}{2}, a_4 = \frac{a_0}{24}, a_6 = \frac{a_0}{720}, \ldots\) ，通式 \(a_{2m} = \frac{a_0}{(2m)!}\) 。奇鏈： \(a_3 = \frac{a_1}{6}, a_5 = \frac{a_1}{120}, \ldots\) ，通式 \(a_{2m+1} = \frac{a_1}{(2m+1)!}\) 。

\[ y = a_0 \sum_{m \ge 0} \frac{x^{2m}}{(2m)!} + a_1 \sum_{m \ge 0} \frac{x^{2m+1}}{(2m+1)!} = a_0 \cosh x + a_1 \sinh x. \]

（ \(\cosh x = \frac{e^x + e^{-x}}{2}\) ——級數法自動重現 4.2 節的指數解，但以雙曲形式。）

觀察：遞推跨兩步，解分成偶、奇兩條獨立鏈—— \(y'' - y = 0\) 在 \(x \to -x\) 下不變，偶解與奇解正是這個對稱性的兩個化身。遞推的跨度 = 方程對稱性的翻譯（呼應 1.3 節 Airy 的三步跳：Airy 方程無此對稱，故跳三步而非對稱分鏈）。

收斂覆核： \(\left|\frac{a_{k+2}}{a_k}\right| = \frac{1}{(k+2)(k+1)} \to 0\) ，收斂半徑無窮——與奇點理論一致（係數解析全域）。∎

歷史啟示#

定義即解：當積分法用盡，「新函數 = 某方程的級數解」成為正當的定義方式。Bessel、Airy、Painlevé 函數都是這樣出生的——數學對象的戶籍制度被級數改寫。
奇點是路標：Cauchy–Hadamard 理論把收斂半徑與複奇點綁定；「看不見的奇點也能砍斷實級數」是複分析強制性的最佳廣告。
統一之美：超幾何函數把幾十種特殊函數收進一個參數家族——Gauss 的 1812 年論文展示了「用最一般的參數族統一特例」的現代數學品味。

與其他章節及領域的關聯#

3.3 節：常點/正則奇點分類是本節方法的導航圖。
5.2 節：Legendre 方程是「常點級數 + 正則奇點多項式化」的完整案例。
5.3 節：Bessel 方程是 Frobenius 法的招牌演出（含對數修正）。
4.1 節：級數解的 Wronskian 檢驗線性無關。
2.3 節：Liouville 不可解定理是級數法登場的歷史理由。
9.1 節：分離變數產生的正交級數（傅立葉級數）是本節「冪級數」的三角親戚。
10.1 節：級數是數值方法的競爭者與互補者：近展開點用級數、遠處用步進。
數論： \(\cosh\) 級數與指數生成函數；Gauss 超幾何級數是生成函數方法（組合數學）的起點之一。

練習#

用冪級數解 \(y'' + y = 0\) ，驗證遞推給出 \(\sin, \cos\) 級數（連接 4.2 節）。
用冪級數解 \(y'' + x y' + y = 0\) ，求兩個線性無關解（提示：其中一條鏈會截斷為多項式——這是 Hermite 方程的親戚現象）。
對 \((1 - x) y'' = y\) ，求級數解並用比值判別法求收斂半徑；用 Cauchy–Hadamard 理論解釋（最近奇點在哪？）。
證明 Airy 方程的兩個級數解 \(y_1, y_2\) 線性無關（取 \(x = 0\) 處的 Wronskian = \(a_0 b_1 - a_1 b_0 \neq 0\) ）。
思考題：係數 \(p(x) = \frac{1}{1 + x^2}\) 在實軸上無窮光滑，但以 \(x_0 = 0\) 展開的解級數收斂半徑只有 1。請從複奇點 \(\pm i\) 的角度完整解釋，並討論：這對「只用實數做分析」的方法論提出什麼挑戰？

進階研讀：級數解的五個深化場景#

場景一：指標方程的三種情況。Frobenius 指標方程 \(r(r-1) + p_0 r + q_0 = 0\) （ \(p_0, q_0\) 為 \(xp, x^2 q\) 在奇點的值）： - 根差非整數：兩個純 Frobenius 級數解； - 根差為整數：大指標解是級數，小指標解可能含對數項； - 重根：必含對數（Bessel \(\nu = 0\) ）。對數項的出現機制：Wronskian \(\sim e^{-\int p} \sim x^{-p_0}\) ，若兩個冪律解的 Wronskian 無法匹配，必須 \(\ln x\) 出場——Abel 恆等式（4.1）預言了對數的必然。

場景二：Kovalevskaya 的迴響。Kovalevskaya 1888 年以剛體旋轉方程的「解皆解析」性質分類了可積情況（第三個可積情況以她命名）；「奇點的性質判別可積性」的傳統延續到 Painlevé 分析（3.3 節）與現代可積系統。一位數學家對「解是否優雅」的研究，成為可積系統理論的奠基石。

場景三：漸近級數——發散但有用。大 \(x\) 時，Airy、Bessel 函數的展開式一般發散（Poincaré 1886 年定義「漸近級數」給出嚴格意義：截斷誤差小於被丟棄首項）。Stokes 現象（1847 年 George Stokes 發現）：漸近展開跨過 Stokes 線時，被忽略的指數小項突然「出現」——這是 WKB（2.3 節）與數學物理的深水區。發散級數不是失敗，是另一種精確。

場景四：級數解的數值使用。Taylor 方法（高階數值積分器，10.1 節）：每步用 ODE 的自動微分算級數係數，步長可放大千倍。現代軟體（如天體力學積分器、Jet Propulsion Laboratory 的星曆計算）用冪級數逐步法達到極高精度——Newton 1665 年的原始方法，仍是最高精度的數值引擎之一。

場景五：生成函數的世界。把離散序列 \(a_n\) 編碼為 \(A(x) = \sum a_n x^n\) ，微分方程變成序列的遞推——組合數學的標準武器。例：Catalan 數滿足 \(C(x) = 1 + x C(x)^2\) （代數方程式級數）；整數分拆的生成函數滿足變係數 ODE。ODE 的級數觀點向上（分析）與離散（組合）兩個方向都開著門。

補充：Sofia Kovalevskaya 速寫#

Sofia Korvin-Krukovskaya（1850–1891）生於莫斯科貴族家庭，少女時代臥室以微積分講義糊牆（傳說由此愛上數學）。因女性不得入大學（俄國），她與 Vladimir Kovalevsky 行「名義婚姻」出國，1874 年哥廷根大學以 Weierstrass 背書授予博士（in absentia）——論文含 PDE 存在性定理（今稱 Cauchy–Kowalevskaya 定理）。1884 年任斯德哥爾摩大學教授（歐洲首位女性數學教授），1888 年以剛體旋轉論文獲巴黎科學院 Bordin 獎（獎金加倍，因評審認為工作超乎尋常）。她的定理名稱今天每天都在 PDE 課堂出現。

補充：Gauss 1812 年論文的收斂判準#

Gauss 在《無窮級數 \(\frac{1 + \frac{ab}{1c} x + \cdots}{\text{的一般研究}}\) 》中給出 \({}_2F_1\) 收斂的完備判準： \(|x| < 1\) 絕對收斂； \(|x| > 1\) 發散； \(|x| = 1\) 時視 \(\operatorname{Re}(a + b - c)\) 正負。這是收斂理論史上第一個「邊界情況精確分類」的成果，開創了「判準要打到邊界」的數學品味。同期 Cauchy（1821）建立一般級數理論，兩人共同把「收斂」從修辭變成定理。

名詞速查表（本節出現的核心術語）#

冪級數解法：設 \(y = \sum a_n (x - x_0)^n\) ，代入遞推。
解析函數：可展開為收斂冪級數的函數。
遞推關係： \(a_{k+2} = f(a_k, a_{k-1}, \ldots)\) ；方程的離散心跳。
Cauchy–Hadamard 定理：收斂半徑 = 到最近複奇點的距離。
Frobenius 法：正則奇點處設 \(y = x^r \sum a_n x^n\) 。
指標方程：定 \(r\) 的最低次方程；讀出奇點附近解的冪律。
對數修正：重根/整數差根時第二解含 \(\ln x\) 。
超幾何函數： \({}_2F_1(a, b; c; x)\) ；經典特殊函數的統一母體。
漸近級數：發散但截斷誤差可控的展開（Poincaré 1886）。
Cauchy–Kowalevskaya 定理：解析 PDE 初值問題的存在唯一（1875）。

年表速查（本節涉及的關鍵年份）#

1665–1676：Newton 以級數解微分方程（首創）。
1812：Gauss 超幾何級數與收斂理論。
1821：Cauchy《分析教程》級數理論。
1838：Airy 方程。
1841：Liouville 不可解定理。
1849–1917：Frobenius（1873 級數法）。
1857：Riemann P-方程。
1875：Cauchy–Kowalevskaya 定理。
1886：Poincaré 漸近級數理論。
1888：Kovalevskaya 剛體旋轉可積性。

深化問答三則#

問：級數解「解出」了方程嗎？還是只是把問題換了形式？ 答：換了形式，但質變：級數可計算（截斷即得數值）、可微分（逐項求導繼續滿足方程）、可分類（收斂半徑、奇點結構、對稱性）。十九世紀數學界接受了「級數就是函數」的新身份（Gauss 1812、Weierstrass 的解析延拓理論是其法典化）。「解」的標準從「寫成已知函數」降級為「寫成可計算的確定對象」——這一步讓微分方程理論覆蓋了全部實際問題。

問：遞推關係裡的自由參數為何恰好是 \(a_0, a_1\) ？ 答：因為二階方程需要兩個初值 \(y(0) = a_0\) 、 \(y'(0) = a_1\) （3.2 節）。級數法與初值理論的這個對應不是巧合：級數展開就是初值問題的「座標系」，每一個自由參數是一個初始導數。更高階方程的級數法有 \(n\) 個自由鏈，一一對應。

問：Frobenius 對數項「醜」嗎？ 答：醜而必要。它反映真實的解析結構：奇點處兩個解的 Wronskian（ \(\sim x^{-p_0}\) ）無法由兩個純冪律生成， \(\ln x\) 是唯一的調和者。而且對數項攜帶資訊：它告訴你解在繞奇點後會變號/混合（單值性破壞，3.3 節場景三）——Bessel 函數的 \(\ln\) 項正是其在原點單值性破壞的數學簽名。最「醜」的項，往往是結構最深的項。

（本節完）