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

指標方程及其三種情形

當你在正則奇點處發動一個弗羅貝尼烏斯級數時,方程最先回給你的不是係數,而是領頭指數 r 的容許值。本篇就來讀那個小小的二次式——指標方程——並說明它兩個根之間的差距如何決定一切,包括那唯一一個必定逼出對數的情形。

指數先於係數出場

在上一篇你已在正則奇點處架好了[[method-of-frobenius|弗羅貝尼烏斯方法]]:不再用樸素的冪級數,而是試一個 y = x^r (a_0 + a_1 x + a_2 x^2 + ...) 形式的解,其中未知的領頭指數 r 容許是任何數——負的、分數的、什麼都行。正是那個多出來的因子 x^r,讓級數得以在奇點旁邊存活下來,而從常點出發的一個赤裸的冪級數在這裡毫無機會。現在我們問那個顯而易見的下一個問題:*r 究竟容許是什麼?*

這裡有一個小小的奇蹟。當你把 y = x^r (a_0 + ...) 代入方程並蒐集 x 的各次冪時,最先出現的最低次冪根本不混入任何更高階的係數——它只是 a_0 乘以一個關於 r 的多項式。由於 a_0 依定義不為零(它是第一個不為零的係數),那個多項式本身就必須為零。所以在你算出任何一個 a_1、a_2、…… 之前,方程就強加了一個只關於 r 的乾淨代數條件。那個條件就是[[ode-indicial-equation|指標方程]],它是守門人:唯有它所容許的 r 值,才可能起頭一個有效的弗羅貝尼烏斯解。

把指標方程讀出來

有一個快速的訣竅,省得你每次都重新推導它。把方程寫成標準的正則奇異形式 x^2 y'' + x p(x) y' + q(x) y = 0,其中 p(x) 與 q(x) 是在 x = 0 附近正常的解析函數。令 p_0 = p(0)、q_0 = q(0)——就只是那兩個函數的*常數項*。那麼指標方程恰好是 r(r - 1) + p_0 r + q_0 = 0。就這樣:一個簡單的二次式。誠實而略令人意外的事實是,只有 p_0 與 q_0 會進來——p 與 q 泰勒展開的其餘部分完全碰不到指數。它們稍後才要緊,是為了係數;而*領頭行為*由兩個數就決定了。

Put the ODE in standard form:

    x^2 y'' + x p(x) y' + q(x) y = 0      ( p, q analytic at x = 0 )

Read the constant terms:

    p_0 = p(0)        q_0 = q(0)

Indicial equation (a quadratic in r):

    r(r - 1) + p_0 r + q_0 = 0

Example - Bessel's equation of order nu:
    x^2 y'' + x y' + (x^2 - nu^2) y = 0
    p_0 = 1 ,  q_0 = -nu^2
    r(r-1) + r - nu^2 = r^2 - nu^2 = 0  ->  r = +nu , r = -nu
只有縮放後係數的常數項 p_0、q_0 抵達指標方程;級數的其餘部分決定係數,而非指數。

花點時間體會那兩個根 r_1 與 r_2 究竟*意味*著什麼。它們是兩個解在奇點附近的領頭指數,所以它們預告了各自的局部形狀:正的 r 表示解一開始很小且乖巧,負的 r 表示它像 1/x^|r| 般爆掉,而分數的 r 表示它帶上像 sqrt(x) 那樣的分支。對貝索方程,兩根算出來是 +nu 與 -nu,正是你將在最後一篇遇到的兩個貝索函數的階。先讀指標根,就像在寫完整句之前先讀出每個解的第一個字。

整齣戲只繫於一個數:r_1 - r_2

把兩根排成 r_1 >= r_2(在它們為實數時)。決定接下來整個故事的單一量,是那個差 r_1 - r_2。這正是[[frobenius-three-cases|弗羅貝尼烏斯根的三種情形]]的核心。較大的根 r_1 *永遠*給出一個乾淨的弗羅貝尼烏斯級數——無論如何你都能倚賴它。整個問題,唯一的問題,是那*第二個*獨立解長什麼樣:較小的根 r_2 會遞給你第二個乾淨的級數嗎,還是機器會卡住、逼出一個對數?那道差 r_1 - r_2 恰恰回答了這一點。

為什麼決定性的數是那個*差*,而不是兩根本身?這個直覺很犀利,值得記牢。當你從較小的根 r_2 建造第二個解時,遞迴會不斷地除以形如 (r_2 + n)(r_2 + n - 1) + p_0(r_2 + n) + q_0 的量——那正是把指標多項式在 r_2 + n 處求值。當 r_2 + n 恰好等於*另一個*根 r_1 時,也就是 n = r_1 - r_2 為正整數時,那個式子就為零。分母出現零,就是機器卡住。所以麻煩不在於單獨的 r_1 或 r_2;而在於從 r_2 一步步加整數,會不會某次正好落回 r_1。那個差,就是全部的答案。

三種情形,逐一走過

現在這套分類幾乎自己就寫出來了。情形一——差 r_1 - r_2 不是整數(這包含像 1/2 這樣的分數或無理數)。這是幸運、輕鬆的情形:從 r_2 一步步加整數永遠碰不到 r_1,所以沒有任何分母會消失。每個根都送上各自的獨立級數,y_1 = x^(r_1)(級數)與 y_2 = x^(r_2)(級數),你只要把它們組合起來即可。順帶一提,複共軛根也落在這裡——它們的差是純虛數,絕不會是正整數。

情形二——兩根相等,r_1 = r_2。一個二重根只能遞給你*一個*級數,所以第二個獨立解必須來自別處——而它總是帶著一個對數現身:y_2 = y_1 ln(x) + x^(r_1)(一個新級數)。這裡的對數並非可有可無,而是被逼出來的,正如常係數與柯西-尤拉方程的重根情形逼出了一個多餘的 t 或 ln x 因子。若你固執地去找兩個平凡級數,你只會找到一個,然後卡在半個通解上。

情形三——兩根相差一個正整數,r_1 - r_2 = N。這是真正微妙的那一個,唯一誠實的說法是:*要看情況,你得自己檢查。*較大的根 r_1 仍給出它乾淨的級數。較小的根 r_2 *可能*給出第二個獨立級數——也可能遞迴在第 N 步撞上那個致命的零,逼出形式 y_2 = c y_1 ln(x) + x^(r_2)(級數)。常數 c 可能恰為零(那就沒有對數,你走運了)或不為零(那就有一個真正的對數)。光憑那個差你無法判斷是哪一種;你必須真的把 c 算出來。提防那個誘人的半真理——整數差*總是*意味著對數,它並不是。

對數從何而來,又為何要緊

當某個情形確實逼出了對數,你不是偶然撞見它——而是刻意把它建進去。有兩條標準路線帶你到達。一是降階法:取已知解 y_1,去找一個 y_2 = u(x) y_1 形式的第二解;代入後問題化簡為關於 u' 的一階方程,而一個 1/y_1^2 的積分悄悄地產生出 ln(x)。另一條路線是把 r 當作連續變數,對 r 微分級數解,再在重根處求值——對 r 微分 x^r 給出 x^r ln x,那正是對數誕生之處。下面是這整套情形分析在實務中如何展開的逐步視圖。

  1. 確認 x = 0 是正則奇點(上一篇的分類),然後把方程寫成標準形式 x^2 y'' + x p(x) y' + q(x) y = 0。
  2. 讀出 p_0 = p(0)、q_0 = q(0),組出指標方程 r(r-1) + p_0 r + q_0 = 0。解出兩根 r_1 >= r_2。
  3. 依 r_1 - r_2 分類:非整數(情形一)、零(情形二),或正整數(情形三)。這唯一一個數,在任何係數之前就告訴你第二解的形式。
  4. 用本階段稍早的遞迴關係,從較大的根 r_1 建造 y_1——這個級數總是管用,且總是收斂到下一個奇點。
  5. 按該情形的規則建造 y_2:第二個乾淨級數(情形一)、被逼出的 y_1 ln(x) + 級數(情形二),或算出常數 c 以判定對數是否存在(情形三)。

為什麼要費這麼多心?因為指標根與其三種情形不是繁瑣記帳——它們是下一篇等著你的特殊函數的結構骨架。情形二的對數型第二解,正是整數階貝索函數取得其「第二類」夥伴的方式;指標根 +nu 與 -nu,正是貝索函數成對出現的原因。最後一個從整個階段帶過來的誠實提醒:即使弗羅貝尼烏斯成功了,它產生的級數也只收斂到最近的那個*另一個*奇點——它的收斂半徑由奇點的幾何決定,而非由你的努力決定。這個方法很強大,但它從不承諾封閉形式,只承諾一個收斂的局部描述。