4.2 常係數齊次方程#

本章導言#

\[ a\, y'' + b\, y' + c\, y = 0. \]

最平凡的方程，最驚人的解法：**猜 \(y = e^{rx}\) **，代入得代數方程

\[ a r^2 + b r + c = 0, \]

解 \(r\) 就解了微分方程。微積分問題降級為中學代數——這一招是 1739 年前後 Euler 的傑作，也是整個微分方程理論中被使用次數最多的一頁。本節講述特徵方程的完整理論（相異實根、重根、複根三案例）、Euler 公式的歷史、以及這套代數化如何預示了現代算子理論。

猜測的邏輯：為什麼是指數函數#

指數函數 \(e^{rx}\) 的特殊性：求導只是乘上常數—— \(\frac{d}{dx} e^{rx} = r e^{rx}\) ，求 \(n\) 次導乘 \(r^n\) 。代入微分方程：

\[ a\, r^2 e^{rx} + b\, r e^{rx} + c\, e^{rx} = (a r^2 + b r + c)\, e^{rx} = 0. \]

由 \(e^{rx} \neq 0\) ，方程等價於 \(a r^2 + b r + c = 0\) 。「求導 = 乘 \(r\) 」把微積分的指數函數變成代數的數 \(r\) ——這個觀察的正式化就是二十世紀的算子記號： \(D = \frac{d}{dx}\) ，方程即 \((a D^2 + b D + c) y = 0\) ，特徵多項式 \(a r^2 + b r + c\) 是多項式 \(aD^2 + bD + c\) 在「 \(D = r\) 」的取值。Boole（1847）與十九世紀英國符號代數學派正式發展了算子演算——先把 \(D\) 當數用（合法嗎？），再補證明（合法，因 \(D\) 滿足與數相同的交換加乘律）。

歷史細節：Euler 在 1739 年致 Johann Bernoulli 的信與《積分法》（1768–1770 出版）中系統發展此法；複數根的處理依賴他 1748 年的

\[ e^{i\theta} = \cos\theta + i \sin\theta. \]

這條「數學中最美的公式」的原始動機之一，正是讓複根的微分方程有意義——微分方程是 Euler 公式的催生婆。

三案例的完整理論#

特徵方程 \(ar^2 + br + c = 0\) 的判別式 \(\Delta = b^2 - 4ac\) 分三案：

**案例一： \(\Delta > 0\) ，相異實根 \(r_1 \neq r_2\) **。兩解 \(e^{r_1 x}, e^{r_2 x}\) ，Wronskian \(= (r_2 - r_1) e^{(r_1 + r_2)x} \neq 0\) ——線性無關。通解：

\[ y = C_1 e^{r_1 x} + C_2 e^{r_2 x}. \]

**案例二： \(\Delta = 0\) ，重根 \(r_1 = r_2 = -\frac{b}{2a}\) **。只有一個解 \(e^{rx}\) ；第二解由降階法（4.1 節）取得：

\[ y_2 = e^{rx} \int \frac{e^{-\int (b/a) dx}}{e^{2rx}}\, dx = e^{rx} \int dx = x\, e^{rx}. \]

通解：

\[ y = (C_1 + C_2 x)\, e^{rx}. \]

「重根長出 \(x\) 因子」不是補丁，而是必然：解空間必須二維（4.1 節），缺的位置由 \(x e^{rx}\) 補上。物理解讀：重根是「臨界阻尼」——系統恰在不振盪與衰減的分界（7.1 節）。

**案例三： \(\Delta < 0\) ，複根 \(r = \alpha \pm i\beta\) **。形式解 \(e^{(\alpha \pm i\beta)x}\) ；由 Euler 公式拆實虛部：

\[ e^{\alpha x} \cos\beta x, \qquad e^{\alpha x} \sin\beta x. \]

通解：

\[ y = e^{\alpha x}\big(C_1 \cos\beta x + C_2 \sin\beta x\big). \]

複數在此不是權宜之計而是必需：振盪（ \(\cos\beta x\) ）與增長/衰減（ \(e^{\alpha x}\) ）的耦合，只有複指數能一筆寫盡。量子力學的波函數、電機工程的相量（phasor）分析，都住活在這個案例裡。

三案例的力學與電學對照#

對 \(m y'' + \gamma y' + k y = 0\) （質量-阻尼-彈簧），定義 \(\omega_0^2 = k/m\) 、 \(2\delta = \gamma/m\) ：

判別	力學名稱	解的行為	應用現象
\(\delta^2 < \omega_0^2\)	欠阻尼	\(e^{-\delta x}\) 乘振盪	減震器太軟、鐘擺
\(\delta^2 = \omega_0^2\)	臨界阻尼	\((C_1 + C_2 x) e^{-\delta x}\)	門關閉器、車用避震的設計目標
\(\delta^2 > \omega_0^2\)	過阻尼	兩個衰減指數	糖漿中的擺、老式電表

工程師的黃金設計：臨界阻尼——回到平衡最快且不振盪。汽車避震器的調校、儀表指針的穩定，都是在特徵方程的判別式上做文章。

高階推廣與範式#

\(n\) 階常係數齊次方程

\[ a_n y^{(n)} + \cdots + a_1 y' + a_0 y = 0 \]

由特徵多項式 \(a_n r^n + \cdots + a_0\) 完全解決：每個單根 \(r\) 貢獻 \(e^{rx}\) ；每個 \(k\) 重根 \(r\) 貢獻 \(e^{rx}, x e^{rx}, \ldots, x^{k-1} e^{rx}\) ；複根按共軛對出現、以 \(\cos/\sin\) 組合。代數基本定理（Gauss 1799 首證）保證 \(n\) 個根（計重數），恰好填滿 \(n\) 維解空間——微積分定理與代數定理在此嚴絲合縫地咬合。

一個具體例題#

例題 4.2.1：解初值問題

\[ y'' + 4 y' + 13 y = 0, \qquad y(0) = 2, \quad y'(0) = 0. \]

解：特徵方程 \(r^2 + 4r + 13 = 0\) ，判別式 \(\Delta = 16 - 52 = -36 < 0\) ：

\[ r = \frac{-4 \pm \sqrt{-36}}{2} = -2 \pm 3i. \]

故 \(\alpha = -2\) 、 \(\beta = 3\) ，通解

\[ y = e^{-2x}\big(C_1 \cos 3x + C_2 \sin 3x\big). \]

初值： \(y(0) = C_1 = 2\) 。求導：

\[ y' = -2 e^{-2x}(C_1\cos 3x + C_2 \sin 3x) + e^{-2x}(-3 C_1 \sin 3x + 3 C_2 \cos 3x), \]

\(y'(0) = -2 C_1 + 3 C_2 = 0 \Rightarrow C_2 = \frac{4}{3}\) 。故

\[ y = e^{-2x}\Big(2 \cos 3x + \frac43 \sin 3x\Big). \]

振幅形式： \(2\cos 3x + \frac43 \sin 3x = R \cos(3x - \phi)\) ， \(R = \sqrt{4 + \frac{16}{9}} = \frac{2\sqrt{13}}{3}\) ， \(\tan\phi = \frac{2}{3}\) 。

解讀：這是阻尼振盪——振幅 \(e^{-2x}\) 指數衰減，包絡半衰期 \(x = \frac{\ln 2}{2} \approx 0.35\) ；振盪頻率 3（每 \(\frac{2\pi}{3}\) 一個週期）。衰減速度與振盪速度的比 \(\frac{\delta}{\beta} = \frac{2}{3}\) ，稱為對數衰減率——測量相鄰兩個峰值的高度比 \(\frac{y(x)}{y(x + 2\pi/3)} = e^{4\pi/3}\) ，實驗上就能反推系統的阻尼參數：特徵方程的根，是可以從實驗數據讀出來的。

驗算： \(r = -2 + 3i\) 應滿足 \(r^2 + 4r + 13 = (4 - 12i - 9) + (-8 + 12i) + 13 = (4 - 9 + 13) + (-12i + 12i) = 0\) ✓。

歷史啟示#

降維打擊的極致：把微分方程化為代數方程，是人類智力史上最成功的「化歸」之一。它成功的根源是指數函數是求導算子的特徵函數（線性代數語言）——好的猜測背後是深的結構。
複數的實用化轉折：複數在 16 世紀（Cardano 1545）被稱為「詭辯量」，在 18 世紀因 Euler 公式與微分方程而成為工程語言。一門數學對象的地位，常由它解決的應用問題決定。
範式遺產：特徵值思想在線性代數（矩陣特徵值）、量子力學（能量本徵態）、譜理論中全面統治——微積分是它的第一戰場。

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

4.1 節：重根的第二解由降階法（Abel 恆等式）產生；線性無關由 Wronskian 保證。
4.3 節：非齊次方程的待定係數法建立在齊次解（本節）之上。
7.1 節：欠阻尼/臨界/過阻尼與共振的物理全貌。
6.1 節：系統版本——矩陣特徵值取代純量特徵值， \(\mathbf{y}' = A\mathbf{y}\) 的解為 \(e^{\lambda t}\) 乘特徵向量；本節是其 \(2\times 2\) 伴隨矩陣特例。
8.2 節：Laplace 變換把本節的「特徵方程」變成「分母多項式」——同一代數，另一種抵達方式。
電機工程：RLC 電路（1.2 節）的暫態響應分析，與本節三案例逐字對應（ \(\gamma \leftrightarrow R\) 、 \(m \leftrightarrow L\) 、 \(k \leftrightarrow 1/C\) ）。
控制理論：特徵根的實部符號 = 穩定性（Re \(r < 0\) 全部成立則穩定）——Routh–Hurwitz 判準（1876/1895）即「不解根判符號」的代數工具，6.3 節的 Lyapunov 思想的線性先驅。

練習#

解 \(y'' - 5y' + 6y = 0\) ，並寫出滿足 \(y(0) = 1, y'(0) = 0\) 的特解。
解 \(y'' + 6y' + 9y = 0\) （重根），驗證 \(y_2 = x e^{-3x}\) 確實是解（直接代入）。
解 \(y'' + 2y' + 5y = 0\) ，把通解寫成振幅-相位形式 \(R e^{-\alpha x}\cos(\beta x - \phi)\) 。
求所有使 \(y'' + by' + 4y = 0\) 臨界阻尼的 \(b\) 值；並說明 \(b = 3\) 與 \(b = 5\) 時解的定性差異。
思考題：RLC 電路 \(L q'' + R q' + \frac{q}{C} = 0\) 。請把「欠阻尼／臨界／過阻尼」翻譯成 \(R, L, C\) 的不等式，並解釋為什麼增大電阻 \(R\) 會讓系統從振盪走向不振盪——這與「阻尼耗能」的直覺一致嗎？

進階研讀：特徵值思想的五個深化場景#

場景一：算子多項式與交換性。 \(D\) 是線性算子，多項式 \(p(D)\) 可代入、可因式分解： \(y'' - 5y' + 6y = (D - 2)(D - 3) y = 0\) 。因式分解給出階梯解法：令 \(z = (D - 3)y\) ，解 \((D - 2) z = 0\) 得 \(z = Ce^{2x}\) ，再解 \((D - 3) y = Ce^{2x}\) （一階線性！）得通解。可交換性 \((D-2)(D-3) = (D-3)(D-2)\) 是整套技巧的合法性根據——這正是 Boole 1847 年符號演算的核心洞察。

場景二：矩陣指數。系統 \(\mathbf{y}' = A\mathbf{y}\) 的解是 \(\mathbf{y}(t) = e^{At} \mathbf{y}_0\) ，其中矩陣指數以冪級數定義（Lagrange 1800s 隱含，Poincaré/Putnam 等標準化）。對角化 \(A = P\Lambda P^{-1}\) 時 \(e^{At} = P e^{\Lambda t} P^{-1}\) ——特徵值 \(\lambda_i\) 扮演本節 \(r_i\) 的角色。矩陣若不可對角化（Jordan 塊），就長出 \(t e^{\lambda t}\) 因子——與重根案例完全同構。線性代數與微分方程，原來是一門課。

場景三：穩定性判準的代數化。 \(n\) 階系統穩定 \(\Leftrightarrow\) 全部特徵根實部為負。Routh–Hurwitz 判準（1876，Routh 的板球運動員身分是科學史趣聞）不求解、只查係數行列式符號；根軌跡法（Evans 1948）把增益變化時根的遷移路徑畫成圖。設計自動駕駛或飛控的人，每天都在移動特徵根。

場景四：複根與振盪的量子版。Schrödinger 方程 \(-\frac{\hbar^2}{2m}\psi'' + V\psi = E\psi\) 中，束縛態對應「特徵值 \(E\) 離散」，散射態對應連續譜——本節的「根決定行為」升級為「譜決定物理」。譜理論（Hilbert 1900s–，von Neumann 1932）的直覺原型就是特徵方程。

場景五：離散化的映照。把 \(y' = \lambda y\) 離散為 \(y_{n+1} = (1 + h\lambda) y_n\) ，穩定條件 \(|1 + h\lambda| \le 1\) ；把 \(y''\) 用中心差分則特徵方程變 \(r^2 - (2 + h^2 c) r + 1 = 0\) 。連續特徵值 \(r\) 與離散放大因子 \(g = e^{rh}\) 的對應，是數值分析（10.1 節）穩定性理論的骨架；實部為負的 \(r\) 在離散世界對應 \(|g| < 1\) ——兩種穩定性的字典。

補充：Euler 公式 1748 的歷史語境#

Euler 在《無窮分析引論》（Introductio in analysin infinitorum, 1748）寫下 \(e^{ix} = \cos x + i\sin x\) ，並視為「指數與三角的深層統一」。他當時對複數的態度已相當現代（「它們服從運算規則並給出正確結果」），但嚴格的複數理論（Argand 平面，1806；Gauss 的嚴格化）還在後頭。微分方程使用者比基礎理論早跑了半個世紀——應用數學的常態節奏。

補充：特徵方程的名稱由來#

「特徵」（characteristic/charactéristique）一詞由 Lagrange 在 1760 年代研究一階偏微分方程時引入（特徵線），後被借用到線性代數（特徵值，characteristic root——德文 Eigenwert 由 Hilbert 1904 定名）。今日「特徵值」「本徵值」「固有值」三種中文譯名並存——術語的翻譯史本身就是一部微縮交流史。

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

特徵方程： \(ar^2 + br + c = 0\) ；把微分方程化為代數。
特徵根： \(r_1, r_2\) ；每個根貢獻一個解。
重根與 \(x\) 因子： \(k\) 重根貢獻 \(x^j e^{rx}\) ， \(j = 0..k-1\) 。
Euler 公式： \(e^{i\theta} = \cos\theta + i\sin\theta\) （1748）。
欠阻尼 / 臨界阻尼 / 過阻尼：判別式的三種符號對應的力學行為。
振幅-相位形式： \(R\cos(\beta x - \phi)\) ； \(R, \phi\) 由初值定。
對數衰減率：相鄰峰值比的自然對數，反推阻尼參數。
算子記號： \(D = d/dx\) ； \(p(D) y = 0\) 的符號演算。
矩陣指數： \(e^{At} = \sum \frac{(At)^n}{n!}\) ；系統版通解。
Routh–Hurwitz 判準：不解根判特徵根實部符號的代數條件。

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

1545：Cardano 首遇複數（「詭辯量」）。
1739：Euler 特徵方程法（通信中成形）。
1748：Euler《無窮分析引論》，Euler 公式。
1760s：Lagrange「特徵」術語。
1799：Gauss 代數基本定理首證。
1847：Boole 算子符號演算。
1876：Routh 穩定性判準。
1895：Hurwitz 判準。
1904：Hilbert 定名 Eigenwert。
1948：根軌跡法。

深化問答三則#

問：猜 \(e^{rx}\) 是不是「碰運氣」？ 答：不是。 \(e^{rx}\) 是求導算子 \(D\) 的特徵函數： \(D e^{rx} = r e^{rx}\) ——正如矩陣 \(A\) 的特徵向量 \(v\) 滿足 \(Av = \lambda v\) 。常係數方程是 \(D\) 的多項式，特徵函數自然是它的「本徵解」。事實上，任何與 \(D\) 交換的算子，其特徵函數都是指數——這是量子力學動量算子本徵態為平面波的同一原理。猜測背後是譜理論。

問：重根為什麼恰好長出 \(x\) 因子，而不是 \(x^2\) 或別的？ 答：因為解空間必須是二維（4.1 節），重根只提供一維 \(e^{rx}\) ，缺的一維由降階法唯一地補為 \(x e^{rx}\) 。用 Jordan 塊語言：重根對應非對角化的 \(\begin{pmatrix} r & 1 \\ 0 & r\end{pmatrix}\) ，冪零部分 \(\begin{pmatrix} 0 & 1 \\ 0 & 0\end{pmatrix}\) 作用兩次即零——所以只長出一層 \(x\) 。重數與 Jordan 結構精確決定多項式因子， \(x^2 e^{rx}\) 屬於三重根。

問：複根時「取實虛部」合法嗎？ 答：合法，且有定理保證：若複值函數 \(y = u + iv\) 滿足實係數線性方程，則 \(u, v\) 分別滿足（對方程取實虛部即可，因係數是實數，共軛保持方程）。故 \(e^{(\alpha + i\beta)x}\) 的實部與虛部都是解，且線性無關（Wronskian \(= \beta e^{2\alpha x} \neq 0\) ， \(\beta \neq 0\) ）。實係數方程的解空間在共軛作用下封閉，複解成對出現、實解各取一半。

（本節完）