4.3 非齊次方程：待定係數與參數變異法#
本章導言#
齊次世界（4.2 節）已經完全馴服；非齊次的世界多了一個「強迫項」 \(g(x)\) ——彈簧被週期性搖動、電路接入交流電、橋梁承受風載。1.2 節的結構定理說：通解 = 一個特解 + 齊次通解。於是全部工作濃縮為一個問題：找到一個特解。兩大兵器登場——待定係數法（undetermined coefficients，聰明的猜測）與參數變異法（variation of parameters，把齊次解的常數鬆綁為函數；Johann Bernoulli 1690s 萌芽、Euler 1739 與 Lagrange 1774–1778 系統化）。本節講述兩法的原理、分工與著名的共振陷阱。
待定係數法：猜得有紀律#
當 \(g(x)\) 是多項式、指數、正弦餘弦及其乘積時，特解會「長得像 \(g\) 」——因為這些函數的導數族是封閉的（ \(e^{kx}\) 求導還是 \(e^{kx}\) 系； \(\cos\beta x\) 求導在 \(\{\cos\beta x, \sin\beta x\}\) 內循環；多項式求導降次）。因此設特解為同族函數的線性組合，係數待定，代入比較。
共振修正（關鍵規則）：若猜測的形式與齊次解撞形（ \(g\) 的族與齊次解空間相交），代入會恆得零、解不出係數——必須**乘 \(x\) **（撞一次乘 \(x\) ，撞 \(s\) 次乘 \(x^s\) ， \(s\) 為 \(g\) 的族與齊次解的重疊重數）。
試探形式速查：
| \(g(x)\) | 試探特解 \(y_p\) |
|---|---|
| \(P_n(x)\) （ \(n\) 次多項式） | \(Q_n(x) = A_n x^n + \cdots + A_0\) |
| \(e^{kx}\) | \(A e^{kx}\) （ \(k\) 為齊次根時 \(A x e^{kx}\) 或 \(A x^2 e^{kx}\) ） |
| \(\cos\beta x\) 或 \(\sin\beta x\) | \(A\cos\beta x + B\sin\beta x\) （永遠兩項都放） |
| \(P_n(x) e^{kx} \cos\beta x\) | 同族組合，含 \(\sin\) ；撞形時乘 \(x^s\) |
方法論評語：待定係數法只對係數為常數、且 \(g\) 族封閉時可用——它是「猜」的紀律化版本。其合法性（為何總能猜中）的深層原因在算子理論： \(L = aD^2 + bD + c\) 作用於這些族時是三對角/封閉矩陣，可逆（除非撞形）——猜測法其實是有限維線性代數。
參數變異法：常數的解放#
對一般 \(g(x)\) （甚至一般變係數齊次部分 \(y'' + p y' + q y = g\) ），參數變異法萬能。設齊次基底 \(y_1, y_2\) （Wronskian \(W \neq 0\) ），把通解中的常數 \(C_1, C_2\) 換成待定函數 \(u_1(x), u_2(x)\) ：
聰明的附加條件：令 \(u_1' y_1 + u_2' y_2 = 0\) （消去二階導數項的麻烦），代入方程整理得
係數行列式恰是 \(W(x) \neq 0\) ！Cramer 法則：
積分即得 \(u_1, u_2\) 。注意三件事：
- 這是 2.1 節場景三「常數變易」的二階版：齊次解 \(C_1 e^{r_1 x} + C_2 e^{r_2 x}\) 的常數被「鬆綁」為函數，用來吸收強迫項。
- 公式只用到 \(y_1, y_2\) 與 \(W\) ——變係數也照用（與待定係數法的根本分工）。
- 附加條件 \(u_1' y_1 + u_2' y_2 = 0\) 不是唯一選擇，而是最省工的規範（gauge choice）——方程組有兩個未知數、只需要一個條件鏈，剩下自由度可以任意分配，選最方便的。
歷史細節：Johann Bernoulli 1697 年在等周問題通信中已用「變異常數」的想法；Euler 1739 年解非齊次方程時系統化；Lagrange 1774–1778 年在月理論（月球攝動）中把它鑄成「variation of constants」方法——天文計算的需求是這個方法的直接推手。它今日在天體力學中的名字是「攝動理論」。
一個具體例題#
例題 4.3.1：解 \(y'' - 2y' + y = e^{x}\) 。
解：
齊次部分：特徵方程 \(r^2 - 2r + 1 = (r - 1)^2 = 0\) ，重根 \(r = 1\) ：
待定係數法（撞形演示）： \(g = e^x\) 。天真猜 \(y_p = A e^{x}\) ——但 \(e^x\) 與 \(x e^x\) 都是齊次解，代入左端恆得 0（它們根本是齊次解！）——這就是撞形。修正：乘 \(x^2\) ，猜 \(y_p = A x^2 e^{x}\) （ \(s = 2\) ，因與二重根撞兩層）。代入：
要等於 \(e^x\) ， \(A = \frac12\) 。故 \(y_p = \frac{x^2}{2} e^{x}\) ，通解
參數變異法覆核：取 \(y_1 = e^x\) 、 \(y_2 = x e^x\) ， \(W = e^{2x} \neq 0\) 。公式：
積分： \(u_1 = -\frac{x^2}{2}\) （取常數為零，併入齊次項）、 \(u_2 = x\) 。故
（ \(u_1 y_1 + u_2 y_2\) 中的額外 \(x^2 e^x\) 項看似與待定係數法不符，實為齊次解的倍數——積分常數的歸併自由度。）
物理解讀：這是共振：強迫頻率與系統固有頻率一致（ \(e^x\) 撞上固有響應 \(e^x\) ），響應按 \(x^2 e^x\) 多項式增長——不是 \(\sin\) 共振的線性增長 \(x\sin\omega_0 x\) ，而是平方增長，因為撞了二重根。共振階數 = 撞形重數（7.1 節的完整物理）。
兩法分工對照表#
| 待定係數法 | 參數變異法 | |
|---|---|---|
| 適用範圍 | 常係數 + \(g\) 為封閉族 | 幾乎一切（需齊次基底） |
| 計算量 | 解小線性組（代數） | 兩次積分（可能很難積） |
| 失效情形 | \(g\) 族不封閉（如 \(\tan x, \ln x\) ） | \(y_1, y_2\) 無法顯式求出 |
| 風格 | 紀律化的猜測 | 通用的機械程序 |
| 淵源 | 算子理論（十九世紀） | 天體力學（Lagrange, 1774–78） |
實務口訣：能用待定係數就用（便宜）； \(g\) 長得怪就參數變異（可靠但費工）。
歷史啟示#
- 特解 + 齊次解 = 通解：非齊次理論把「尋找全部解」化為「尋找一個解」——找一個比找全部容易，這是線性結構的巨大紅利。
- 天文學是方法之母：參數變異法為了算月球攝動而生；今日 Kalman 濾波（2.3 節）、衛星軌道修正都是它的後代。觀測壓力驅動方法發明。
- 共振的教訓：撞形（共振）不是方法失效，而是物理在喊話——系統的固有頻率被踩中了。數學的「例外情況」常是物理的「重點情況」。
與其他章節及領域的關聯#
- 4.1 節：參數變異的公式分母是 Wronskian；Abel 恆等式使其可顯式計算。
- 4.2 節：待定係數的「撞形修正」完全由齊次根結構（特徵方程）決定。
- 7.1 節：共振的物理（能量累積、橋梁事故）。
- 8.2 節：Laplace 變換解非齊次 IVP——待定係數的「撞形修正」對應 Laplace 域的部分分式重根情形。
- 9.1 節：PDE 的非齊次問題（熱源、波源）用同樣的「特解 + 齊次」結構；Green 函數是參數變異的無窮維版。
- 天體力學：攝動理論 = 參數變異 + 級數展開；月球近地點計算（2.2 節 Clairaut 插曲）是歷史現場。
- 控制工程：系統對輸入 \(g\) 的響應 = 暫態（齊次解，由初值決定）+ 穩態（特解，由輸入決定）——「暫態死、穩態留」是工程師的本節語言。
練習#
- 用待定係數法解 \(y'' - y = 2 e^{x} + 3\) （注意 \(e^x\) 撞形）。
- 用參數變異法解 \(y'' + y = \sec x\) （ \(0 < x < \frac{\pi}{2}\) ）——此 \(g\) 無法用待定係數。
- 解 \(y'' + 4y = \cos 2x\) （純共振），寫出通解並描述解隨時間增長的行為。
- 對 \(y'' + p y' + q y = g\) ，證明：若 \(y_{p,1}\) 對應 \(g_1\) 、 \(y_{p,2}\) 對應 \(g_2\) ，則 \(y_{p,1} + y_{p,2}\) 對應 \(g_1 + g_2\) （非齊次疊加原理）。
- 思考題：工程上「共振危險」通常以 \(\sin\) 強迫（ \(y'' + \omega_0^2 y = F\sin\omega_0 t\) ，響應 \(\propto t\cos\omega_0 t\) ）討論。請解釋：為什麼有阻尼時共振響應不再無界，而成為有限的峰值 \(y_{\max} \sim \frac{F}{2\delta\omega_0}\) ？（可用 7.1 節的正規化語言預習。）
進階研讀：特解思想的五個深化場景#
場景一：算子逆與 Green 函數。 \(L[y] = g\) 的特解 = \(L^{-1}[g]\) 。 \(L^{-1}\) 是積分算子（4.1 節場景四）： \(y_p(x) = \int G(x, \xi) g(\xi)\, d\xi\) 。參數變異法的公式實際上顯式構造了 \(G\) ： \(G = \frac{y_1(\xi) y_2(x) - y_1(x) y_2(\xi)}{W(\xi)}\) 。待定係數是在函數族裡查表，Green 函數是把查表變成積分——後者可以推廣到 PDE 與邊值問題，成為數學物理的主幹。
場景二：Duhamel 原理。非齊次初值問題的解可寫成
其中 \(G\) 是脈衝響應（8.3 節的 \(\delta\) 輸入的響應）。任何強迫 = 脈衝的連續堆疊，響應 = 脈衝響應的捲積——Duhamel（1833）的原理把「求對任意輸入的響應」化為「求對單一脈衝的響應」。Laplace 變換（第 8 章）的捲積定理是它的代數化。
場景三：共振與正規模。多自由度系統（耦合振子）的強迫響應在正規模頻率處共振——結構工程的地基理論。Taco Bell 大樓倒塌式笑話之外，真實案例是 1940 年 Tacoma 海峽大橋（氣動自激，嚴格說是負阻尼失穩而非純共振——教科書常混談，7.2 節將釐清自激振盪與共振的差異）。「共振」一詞的精確使用，是工程師的道德責任。
場景四：參數變異的現代身份——規範自由度。條件 \(u_1' y_1 + u_2' y_2 = 0\) 是一種規範選擇（gauge fixing）：表示 \((u_1, u_2) \to (u_1 + a, u_2 - a y_1/y_2)\) 不改變 \(y_p\) ，我們用一個規範條件凍結冗餘。這與電磁學的規範變換（ \(A \to A + \nabla\chi\) ）同構——同一個「表示有冗餘，用條件凍結」的思想，從 Lagrange 的月球計算走到 Yang–Mills 理論。
場景五：非齊次的離散與控制。差分形式 \(y_{n+1} = a y_n + u_n\) 中 \(u_n\) 是控制輸入；「特解 = 穩態響應」在數位濾波器設計中是頻率響應分析（ \(z\) 域）。本節的線性結構（特解 + 齊次、Duhamel 捲積）是所有線性系統理論（濾波、控制、通訊）的共同骨架——一句話：線性系統的全部本領，都是本節兩個方法的無窮維擴充。
補充：Lagrange 的月球計算與方法誕生#
1772–1778 年間，Lagrange 為巴黎科學院月球理論競賽發展「常數變異法」：把二體軌道參數（橢圓要素）視為時間函數，攝動力導出其演化方程。這一視角（軌道固定、要素漂移）今日稱經典攝動理論，仍是太空任務設計（軌道轉移、引力助推）的標準工具。方法誕生的工資單：Lagrange 因月球理論工作獲得科學院大獎與腓特烈大帝的柏林科學院職位——微積分技巧的經濟回報，從來都很實在。
補充：待定係數法「為何總能猜中」的代數證明#
以 \(g = e^{kx}\) 為例。 \(L[e^{kx} x^n e^{kx}]\) 全部落入 \(\mathrm{span}\{x^j e^{kx}\}_{j \le n}\) ——因為 \(D\) 對此族封閉。故 \(L\) 限制在 \(n+1\) 維族上是矩陣 \(M\) ；只要 \(k\) 非特徵根， \(M\) 可逆，方程 \(M \mathbf{A} = \mathbf{e}_n\) 有唯一解——待定係數總能解出。 \(k\) 是特徵根（撞形）時 \(M\) 奇異，需加大族（乘 \(x^s\) ）直到可逆——「猜得中」= 有限維封閉 + 可逆性，整個方法是線性代數的偽裝。
名詞速查表（本節出現的核心術語）#
- **特解 \(y_p\) **：滿足非齊次方程的單一解。
- 通解結構： \(y = y_p + y_h\) （特解 + 齊次通解）。
- 待定係數法：依 \(g\) 的函數族設試探解、代入定係數。
- 撞形修正：試探形式與齊次解重疊時乘 \(x^s\) 。
- 參數變異法： \(y_p = u_1 y_1 + u_2 y_2\) ， \(u_i' = \mp \frac{y_i g}{W}\) 。
- 共振：強迫與固有頻率一致；響應含多項式因子。
- Duhamel 原理：任意輸入的響應 = 脈衝響應的捲積。
- Green 函數： \(L^{-1}\) 的積分核。
- 暫態 / 穩態：齊次解（衰減）與特解（持續）的工程語言。
- 攝動理論：參數變異在天體力學中的名字。
年表速查（本節涉及的關鍵年份）#
- 1697：Johann Bernoulli 的常數變異萌芽。
- 1739：Euler 系統化非齊次解法。
- 1774–1778：Lagrange 參數變異法與月球理論。
- 1828：Green 函數。
- 1833：Duhamel 捲積原理。
- 1847：Boole 算子法（待定係數的代數基礎）。
- 1940：Tacoma 海峽大橋事故（振動理論的公眾教訓）。
- 1960s：Kalman 濾波（參數變異思想的自動控制版）。
深化問答三則#
問：參數變異法中，為什麼可以「隨意」附加條件 \(u_1' y_1 + u_2' y_2 = 0\) ？ 答：因為表示有冗餘： \((u_1, u_2)\) 與 \(y_p\) 不是一一對應（加上齊次解的倍數不變）。我們有兩個未知函數、只需滿足一個（原）方程——剩一個自由度。附加條件把它凍結在最方便的規範上，使 \(y_p''\) 的展開中 \(u_1'', u_2''\) 項全消。自由的正確用法，是把它換成簡單。
問：待定係數法的撞形與物理共振，是巧合嗎？ 答：不是巧合，是同一件事的兩面。撞形 = 強迫項 \(e^{kx}\) （或 \(\sin\omega x\) ）落在齊次解空間中 = 系統的固有模式與輸入同頻。此時代數上可逆性崩壞、特解必須含 \(x\) 因子；物理上能量逐週期淨輸入、振幅線性（或平方）增長。代數的奇異性就是動力學的共振——線性理論最漂亮的一次數學與物理的重合。
問：如果齊次解根本求不出顯式（如變係數 Bessel 方程），參數變異還有用嗎？ 答：有用，只是形式化：公式 \(u_i' = \mp \frac{y_i g}{W}\) 仍成立（用級數表示的 \(y_1, y_2\) ），特解可寫為級數或積分表達——這正是特殊函數論（第 5 章）與 Green 函數法的日常。方法的效力不依賴「解得漂亮」，只依賴「解存在且線性無關」——4.1 節的結構定理在此兌現全部價值。
（本節完）