需要填補的那道缺口
到此為止，本節的頭條已成第二天性：a y'' + b y' + c y = 0 的解填滿一個二維空間，所以你需要兩個獨立解才能寫出通解。對常係數，特徵方程一次就把兩個都交給你，還整齊地分成三種情形。但這份禮物綁在常係數上。一旦有係數隨 x 變化——像 x^2 y'' - 3x y' + 4y = 0，或物理中著名的 Bessel 與 Legendre 方程——猜 y = e^(rx) 便瓦解，你只能瞪著一個代數把戲碰不到的問題。
然而你往往能用別的法子弄到*一個*解——一次幸運的猜測、一段冪級數計算、一個已知的特殊函數，或者在常係數世界裡，重根所吐出的那唯一一個指數函數。麻煩在於，一個解只張成一條過原點的直線，而非整張平面：只有一半答案。降階法正是這樣一種方法：它拿那一個解，記作 y1，製造出第二個保證與它獨立的解 y2。從單一顆種子，長出那缺失的維度。
核心想法：讓常數變成函數
整個把戲一句話講完。你早已知道，對任意常數 c，c * y1 都是解——但常數倍數*不*獨立於 y1，對建立一組基毫無用處。那就放鬆這個常數：不用 c，改允許一個*函數* v(x)，並尋找形如 y2 = v(x) * y1(x) 的第二個解。期望在於，v 的自由度能讓你既滿足方程，又使 y2 真正有別於 y1。這一招——把常數提升為未知函數——正是本法的核心，也與你將在非齊次階梯遇到的參數變異法出於同一直覺。
為何如此粗率的猜測會奏效？因為把 y2 = v * y1 代入方程，會生出某種神奇的東西：一條關於 v 的新微分方程，其中純粹的 v 項（不含 v 任何導數的那一項）完全消去。它必然消去——其係數恰是原方程在 y1 處的值，而那正因 y1 是解才為零。剩下的只含 v'' 與 v'，絕不含 v 本身。於是若你令 w = v'，那條關於 v 的二階方程便化為關於 w 的*一階*方程。名字便由此而來：階數被降了，從二降到一。
配方，一步一步來
先把方程化成標準型——整條除一除，使 y'' 的係數為 1：y'' + p(x) y' + q(x) y = 0。其中的 p(x) 正是本法在意的部分。接著執行代換 y2 = v y1，看著階數下降。以下是這套程序的明晰清單。
- 把方程寫成標準型 y'' + p(x) y' + q(x) y = 0，並備好已知的解 y1。
- 代入 y2 = v y1。算出 y2' = v' y1 + v y1' 與 y2'' = v'' y1 + 2 v' y1' + v y1'，再代進方程。
- 整理同類項。v 的係數是 y1'' + p y1' + q y1 = 0（因為 y1 解此方程），於是純粹的 v 消失，只剩 v'' 與 v'。
- 令 w = v'。剩下的方程關於 w 是一階線性的——其實可分離變數——故解之：w = (1/y1^2) * exp(-∫ p dx)。
- 再積一次，從 w = v' 還原 v，然後組出 y2 = v y1。捨去任何相加常數（它只會重建出 y1 的倍數）與任何前導常數因子。
第四、五步可壓縮成一條值得記住的公式，即降階公式：y2 = y1 * ∫ [ exp(-∫ p dx) / y1^2 ] dx。你可以背它，但更健康的做法是用代換重新推導幾次，直到那次消去感覺是必然而非魔法。注意分母裡的 y1^2——這正是本法為何需要 y1 在你所工作的區間上不為零，一條值得揣在口袋裡的但書。
一個演算示意與一位隱藏的朋友
取一個 e^(rx) 把戲派不上用場的變係數例子：x^2 y'' - 3x y' + 4y = 0，x > 0。假設有人告訴你 y1 = x^2 可行——你心算就能驗證：y1' = 2x、y1'' = 2，而 x^2*(2) - 3x*(2x) + 4*(x^2) = 2x^2 - 6x^2 + 4x^2 = 0。好。現在化標準型就是除以 x^2：y'' - (3/x) y' + (4/x^2) y = 0，故 p(x) = -3/x。p 的積分為 -3 ln x，於是 exp(-∫ p dx) = exp(3 ln x) = x^3。則 w = x^3 / (x^2)^2 = x^3 / x^4 = 1/x，而 v = ∫ 1/x dx = ln x。因此 y2 = v y1 = x^2 ln x。
Equation: x^2 y'' - 3x y' + 4y = 0 (x > 0) Known: y1 = x^2 Standard form: y'' - (3/x) y' + (4/x^2) y = 0 -> p(x) = -3/x Factor: exp(-INT p dx) = exp(3 ln x) = x^3 Reduce: w = v' = x^3 / (x^2)^2 = 1/x Integrate: v = INT (1/x) dx = ln x Second solution: y2 = v*y1 = x^2 ln x General: y = c1*x^2 + c2*x^2 ln x
於是通解為 y = c1 * x^2 + c2 * x^2 ln x。那個對數絕不可能從任何指數猜測冒出來；降階法找到它，是因為它向方程提了個更聰明的問題。而 ln x 這個因子當下就敲定獨立性：x^2 ln x 不是 x^2 的常數倍數，故這對的 Wronskian 不為零，它們構成一組真正的基本解組。事實上這個保證早已內建於本法，這就引出藏在公式裡的那位朋友。
再看一眼 exp(-∫ p dx) 這一項。前幾篇裡你恰恰見過這個式子，那是 Wronskian 的 Abel 恆等式：W(x) = W(x0) * exp(-∫ p dx)。這絕非巧合。降階法之所以奏效，是因為 v' = W / y1^2，其中 W 是 y1 與你正在建造的 y2 的 Wronskian。既然 Abel 恆等式說 W 永不為零（只要起初不為零），v' 便永不為零，所以 v 確實非常數——而這*正是* y2 獨立於 y1 的陳述。本法不可能不小心交給你一個 y1 的倍數；獨立性在結構上就被保證了。
它何時咬人，又通往何處
對它的侷限要看得清楚。其一，你必須已經擁有 y1——降階法是一台*第二*解的機器，不是第一解的神諭，而對真正困難的變係數方程，連找到 y1 都可能要動用後面階梯的冪級數或特殊函數理論。其二，公式要除以 y1^2，故在 y1 = 0 之處會絆倒；你通常在 y1 保持定號的區間上工作，必要時再跨越零點分段拼接。其三，那個倖存的積分 ∫ [exp(-∫ p dx)/y1^2] dx 也許沒有初等的封閉形式——這沒關係。一個由誠實積分定義的解，依然是一個解。
也值得說清楚降階法*不是*什麼。它不是求解非線性方程的法門——驅動它的那次消去，完全依賴線性、齊次的結構，正因如此「y1 是解」才能讓一整塊項消失。塞進一個 (y')^2，把戲就分崩離析，正如第一篇裡疊加原理所遭遇的那樣。這項技術，在線性之內方能存活，離開便亡。
有了這一篇，齊次這一節就圓滿了：你能讀出解空間的維數，用 Wronskian 與 Abel 恆等式檢驗獨立性，透過特徵方程的三種風味解盡每一個常係數情形，而如今——即便面對變係數——也能從單一顆種子建出一組完整的基。同樣那個「讓常數變動」的想法，正是通往下一節的門扉：那裡右邊坐著一個強迫項，而參數變異法把這把戲化為求解非齊次方程的通用方法。與其說你正闔上一章，不如說你正磨利那把要帶進下一章的工具。