2.3 Bernoulli 與 Riccati 方程#

本章導言#

線性方程可以系統求解（2.1 節），非線性方程一般不能。但在兩者之間，有一批「偽裝成非線性的線性方程」：透過一次變數代換就能化歸為線性。這一類方程在 1690–1730 年間被 Bernoulli 家族、Leibniz 與 Riccati 逐一馴服，成為微積分黃金時代的競技場。本節講述 Bernoulli 方程與 Riccati 方程的故事——前者展示「代換化歸」的典範，後者展示「一個特解開萬鎖」的深招，並以 Liouville 1841 年的否定結果收尾：化歸之路，終有盡頭。

Bernoulli 方程：1695 年的公開挑戰#

1695 年，Jacob Bernoulli（雅各布・伯努利，1654–1705） 在《教師學報》上提出：

\[ \frac{dy}{dx} + p(x)\, y = q(x)\, y^n \qquad (n \neq 0, 1), \]

並坦言此方程超出他的能力，公開徵解。這是一階方程中最著名的「擾動線性」類型：多了 \(y^n\) 項。

解法（Leibniz 密授，Johann Bernoulli 1696 年公開）：當 \(y \neq 0\) 時，除以 \(y^n\) ：

\[ y^{-n} y' + p\, y^{1-n} = q. \]

令 \(v = y^{1-n}\) 。鏈鎖律給 \(v' = (1 - n) y^{-n} y'\) ，故 \(y^{-n} y' = \frac{v'}{1 - n}\) ，方程化為

\[ \frac{v'}{1 - n} + p\, v = q \quad\Longrightarrow\quad v' + (1 - n) p\, v = (1 - n) q, \]

——一階線性方程！用 2.1 節積分因子完成。整個方法的靈魂：把 \(y\) 的非線性冪，轉移為新變數 \(v\) 的線性身份。注意 \(y \equiv 0\) 是被除法排除的恆定解（ \(n > 0\) 時），須另行撿回。

歷史花絮：Leibniz 在 1695 年 12 月的私人信件中把代換法告訴 Johann Bernoulli；後者於 1696 年 3 月公開發表完整解法。Jacob 對「我出的題由弟弟經萊布尼茲指點解出」頗為不快——Bernoulli 兄弟的恩怨（1697 年最速降線之爭徹底決裂）是十八世紀數學最著名家庭劇。但學術本身是贏家：公開挑戰制催生了整套一階方程理論。

Riccati 方程：二次非線性的硬度#

1720 年代，威尼斯貴族數學家 Jacopo Francesco Riccati（里卡蒂，1676–1754）（其子 Vincenzo 亦為數學家，編輯了 Newton《原理》的義大利譯本）研究

\[ \frac{dy}{dx} = a(x) + b(x)\, y + c(x)\, y^2. \]

這是「最簡單的不可分離非線性方程」： \(y^2\) 項使它無法分離、也無法直接線性化。D'Alembert 1763 年命名之（今日通稱 Riccati 方程）。

核心定理（已知一個特解即可解）：設 \(y_1(x)\) 是一個特解。令

\[ y = y_1 + \frac{1}{u}. \]

代入（展開時 \(y_1^2\) 項因 \(y_1\) 是解而相消），整理後得

\[ u' = -\big(b + 2 c\, y_1\big) u - c, \]

——線性！解出 \(u\) 即得通解。兩層教訓：

一個特解 = 一把鑰匙：非線性方程的困難集中在「找到立足點」；一旦有特解，整個解族盡收囊中。這一招在二階線性方程中再現（4.1 節降階法）。
為何特解有效： \(y^2\) 項是唯一的「麻煩製造者」，而 \(y_1\) 恰好使 \(a + b y_1 + c y_1^2 = 0\) ——代入 \(y = y_1 + \frac1u\) 時，二次項以「 \((y_1 + \frac1u)^2 - y_1^2 = \frac{2 y_1}{u} + \frac{1}{u^2}\) 」的形式出現，其中 \(\frac{1}{u^2}\) 項恰好與 \(c\,\frac{y'}{\cdots}\) 匹配消去。讀者應親自展開驗證一次，體會「特解吸收非線性」的機制。

Riccati 方程與二階線性方程的等價：令 \(y = -\dfrac{u'}{c\, u}\) ，Riccati 方程化為 \(u\) 的二階線性方程（Riccati 變換）。這個雙向橋樑意義重大：

\[ \underbrace{y' = a + b y + c y^2}_{\text{一階非線性}} \;\Longleftrightarrow\; \underbrace{c\, u'' - (b + \tfrac{c'}{c})\, u' + a c\, u = 0}_{\text{二階線性（常係數情形：} c u'' - b u' + a c u = 0\text{）}}. \]

Riccati 方程解的性質（零點、爆破、振盪）直接翻譯為二階方程解的性質——第 5 章 Bessel、Legendre 方程與其 Riccati 伴侶的互換，正是二十世紀數學物理的常用武器（如量子力學中的 supersymmetric quantum mechanics，1980 年代）。

Liouville 的否定定理：化歸的邊界#

Riccati 方程看似可以「再努力一點」解出？Joseph Liouville（劉維爾，1809–1882） 在 1841 年證明：

對一般係數的 Riccati 方程，其解不能用初等函數與初等積分表達。

這是數學史上第一個「不可解性」的嚴格證明（採用他新創的、今稱 Liouville 理論的手段——後發展為微分代數，Ritt 與 Kolchin 1940–1950 年代完成）。否定結果的歷史意義不下於任何正結果：

它宣告「找到萬能積分法」的夢想死亡；
迫使數學家轉向級數解（第 5 章）、定性方法（第 6、7 章）與數值方法（第 10 章）；
其思想在 1880 年代被 Picard、Vessiot 發展為微分 Galois 理論，今日仍是活躍領域。

「解不出來」也是定理——這個觀念在 1841 年是革命性的。

一個具體例題#

例題 2.3.1：解 Bernoulli 方程

\[ \frac{dy}{dx} - \frac{y}{x} = -\frac{y^2}{x^2}, \qquad x > 0. \]

解：型如 \(y' + p y = q y^n\) ， \(p = -\frac1x\) ， \(q = -\frac{1}{x^2}\) ， \(n = 2\) 。

除以 \(y^2\) （並撿回恆定解 \(y \equiv 0\) ）：

\[ \frac{y'}{y^2} - \frac{1}{x y} = -\frac{1}{x^2}. \]

令 \(v = y^{1 - n} = y^{-1}\) ， \(v' = -y^{-2} y'\) ，即 \(\frac{y'}{y^2} = -v'\) ：

\[ -v' - \frac{v}{x} = -\frac{1}{x^2} \quad\Longrightarrow\quad v' + \frac{v}{x} = \frac{1}{x^2}. \]

線性方程，積分因子 \(\mu = e^{\int dx/x} = x\) ：

\[ (x v)' = \frac{1}{x} \quad\Longrightarrow\quad x v = \ln x + C \quad\Longrightarrow\quad v = \frac{\ln x + C}{x}. \]

代回 \(v = 1/y\) ：

\[ y = \frac{x}{\ln x + C}. \]

連同恆定解 \(y \equiv 0\) ，通解完備。檢驗： \(y \to \infty\) 當 \(\ln x + C \to 0\) ——解有爆破點 \(x = e^{-C}\) ，再次見證非線性的「有限壽命」現象（1.2 節）。

再看 Riccati 路線：此方程實為 Riccati 型（ \(a = 0\) ， \(b = -\frac1x\) ， \(c = -\frac{1}{x^2}\) ），且有特解 \(y_1 = \frac{x}{?}\) ……讀者可驗證 \(y_1 = \frac{x}{1 + \ln x + C_1}\) 並不「現成」；但另一路線 \(y = -\frac{x^2 u'}{(-1)\, u} = x^2 \frac{u'}{u}\) 的 Riccati 變換將其化為二階線性方程 \(x^2 u'' - x u' - \frac{1}{x^2}u = 0\) ，屬歐拉型方程（5.1 節進階）——可見同一方程有多重身份，會分類的人，選武器的自由度大。

歷史啟示#

家族與競爭塑造學科：Bernoulli 家族三代（Jacob、Johann、Daniel 與 Johann II）在 1690–1760 年間貢獻了當時一階方程理論的大半——競爭性家族 + 期刊挑戰制 = 高產出的數學生態。
化歸的系統化：分離（1691）→ 齊次代換（1691）→ 積分因子（1694）→ Bernoulli 代換（1696）→ Riccati 特解法（1720s）——四十年間歐陸建立了一整套「化為已知」的技術目錄。這份目錄至今仍在一階方程的教科書裡，原封不動。
否定結果的尊嚴：Liouville 1841 證明不可解，使數學家放下執念、開拓新戰場。科學的進步常以「放棄某類希望」為標誌。

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

2.1 節：Bernoulli 方程化歸的終點是線性方程與積分因子。
2.2 節：Riccati 方程一般非正合、無積分因子——正合理論的邊界案例。
3.2 節： \(y^2\) 項使 Lipschitz 條件在全域失效（但局部成立）——Riccati 解的爆破與此呼應。
4.1 節：「一個特解降階」在二階齊次線性方程中的版本（降階法）。
5.1 節：Liouville 不可解定理是冪級數解法登場的歷史理由——積分法窮盡之處，級數法接棒。
7.1 節：Riccati 變換把振幅方程化為相位的 Riccati 型方程（現代光學與量子力學中的 Ermakov–Pinney 方程，1883/1950）。
控制理論：矩陣 Riccati 方程（Rudolf Kalman, 1960）是濾波理論（Kalman 濾波）的核心，工程界每年執行億萬次——一個 1720 年代的方程在衛星導航系統中每天上崗。

練習#

解 \(y' + y = x\, y^3\) （ \(n = 3\) ），寫出代換與線性方程各步。
解 \(\dfrac{dy}{dx} = \dfrac{y}{x} + \dfrac{y^2}{x^2}\) ，並指出爆破行為。
已知 \(y_1 = x\) 是 \(y' = 1 + \dfrac{y}{x} + \dfrac{y^2}{x^2}\cdot 0 + \dfrac{y^2 - x^2}{2 x^2}\) 之簡化型 \(y' = 1 + \frac{y}{x} + \frac{y^2 - x^2}{2x^2}\) 的特解（請先驗證），用 Riccati 特解法求通解。
用 Riccati 變換 \(y = \dfrac{u'}{u}\) 把 \(y' + y^2 = 1\) 化為二階線性方程，並驗證 \(y = \tanh x\) （滿足 \(y' = 1 - y^2\) ，符號不同）與 \(y = \coth x\) 的角色差異。
思考題：Liouville 證明「一般 Riccati 方程不可初等積分」。請討論：這一否定結果如何改變了後世教科書的編排（為何第 5 章要講級數解、第 10 章要講數值法）？

進階研讀：非線性代換的五個深化場景#

場景一：代換的分類學。一階方程的主要代換： \(v = y/x\) （齊次，1.3）、 \(v = y^{1-n}\) （Bernoulli）、 \(v = y + \frac{1}{u}\) （Riccati 特解法）、 \(v = y'\) (降階)。可見代換的目的是改變未知數的身份而非自變數。Sophus Lie 把「所有能把方程化簡的變換」組織成群——李對稱分析（1880s）今日是尋找特殊代換的系統方法，取代了十九世紀前的「靈感式代換」。

場景二：爆破的定量研究。 \(y' = a + b y + c y^2\) （常係數）的解 \(y = \frac{y_+ (x - x_0)\, e^{(y_+ - y_-)x} \cdots}{}\) 型，爆破時間可由初值顯式算出。變係數時爆破時間只能估計（比較定理，Kamke 1930s）。金融上資產泡沫的「有限時間奇點」模型（Sornette 1990s–2000s 的對數週期冪律）直接借用 Riccati 式爆破數學——微分方程的奇點成為預測危機的語言。

場景三：Riccati 與量子力學。Schrödinger 方程 \(\psi'' + \big(E - V(x)\big)\psi = 0\) 經 \(\psi = e^{-\int W}\) （WKB 變換，1926 年 Wentzel–Kramers–Brillouin）化為 \(W' = W^2 - (E - V)\) ——Riccati 方程！WKB 近似（半經典物理的骨幹）本質上是 Riccati 方程的漸近分析。1983 年起的超對稱量子力學更把「找特解（基態）→ 生成整個解族」的思想工程化——Riccati 的 1720 年代技巧，是量子光學的日常工具。

場景四：Kalman 濾波的矩陣 Riccati。矩陣版 \(P' = A P + P A^T - P C^T R^{-1} C P + Q\) 源於 1960 年 Kalman 的最佳濾波理論：狀態估計誤差協方差 \(P\) 的演化。阿波羅登月的導航（1969）、GPS、無人機控制都內建此方程的即時求解。同一個二次非線性，在力學裡叫爆破，在資訊裡叫最佳估計。

場景五：微分代數的今天。Liouville 理論經 Ritt（1948）、Kolchin（1973）成為微分代數；Risch 演算法（1969–）讓電腦代數系統（Mathematica、SymPy）自動判斷「這個積分有沒有初等原函數」。你在電腦上按積分鍵得到的「不可積分」訊息，背後是 Liouville 1841 年定理的演算法化——否定定理最終變成了軟體。

補充：Bernoulli 家族速寫#

Jacob（1654–1705）：巴塞爾大學數學教授；機率論《推測術》（1713 出版，大數定律）；懸鏈線、等角螺線、Bernoulli 方程。墓碑刻對數螺線。
Johann（1667–1748）：Jacob 之弟、格羅寧根後巴塞爾教授；最速降線（1697）發起者之一；教出 Euler 與 L'Hôpital（1696 年 L'Hôpital 的微積分教科書實含 Johann 的有償供稿——史上第一本微積分教科書的「幽靈寫作」公案）。
Daniel（1700–1782）：Johann 之子；流體力學 Bernoulli 定律（1738）、弦振動疊加原理爭論（1753）、Gas kinetic 理論先驅；十次法蘭西科學院獎。
家族三代六位重要數學家，恩怨與成就同樣著名：Johann 曾因兒子 Daniel 的《流體力學》早於自己出版而將日期倒填——天才家族的內耗與傳承並行。

補充：Riccati 方程的可分離特例#

當係數皆常數（ \(a, b, c\) 常數）時，Riccati 方程可以初等解出：令 \(y = z + \frac{y_\pm}\) （ \(y_\pm\) 是代數方程 \(c y^2 + b y + a = 0\) 的兩根，設判別式 \(\geq 0\) ），化為

\[ z' = c\, z\, (z + y_+ - y_-), \]

可分離！解為指數型 \(z = \frac{y_+ - y_-}{K e^{-(y_+ - y_-)c x} - 1}\) 。可見「常係數二次非線性」其實溫馴；困難全部來自係數隨 \(x\) 變化——與線性情形（變係數需要級數解，第 5 章）形成平行結構。

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

Bernoulli 方程： \(y' + p y = q y^n\) ；代換 \(v = y^{1-n}\) 化線性。
Riccati 方程： \(y' = a + b y + c y^2\) ；一階非線性。
特解法：已知一特解時， \(y = y_1 + \frac1u\) 化為線性。
Riccati 變換： \(y = -\frac{u'}{c u}\) ，連接一階非線性與二階線性。
Liouville 不可解定理：一般 Riccati 方程不可初等積分（1841）。
微分代數：研究函數域上微分運算的代數結構；Ritt–Kolchin 理論。
WKB 變換： \(\psi = e^{-\int W}\) ，把線性二階方程化為 Riccati 型。
Kalman 濾波：矩陣 Riccati 方程驅動的最佳狀態估計。
爆破時間：解趨於無窮的有限時刻。
比較定理：用已知方程夾逼未知解的行為（Kamke）。

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

1695：Jacob Bernoulli 提出方程並公開徵解。
1695 年 12 月：Leibniz 密信 Johann Bernoulli 提示代換。
1696：Johann Bernoulli 公開解法。
1713：Jacob《推測術》出版（大數定律）。
1720s：Jacopo Riccati 研究二次方程。
1763：d'Alembert 命名「Riccati 方程」。
1841：Liouville 不可解定理。
1883：Ermakov 方程。
1926：WKB 近似。
1960：Kalman 濾波與矩陣 Riccati 方程。
1969：Risch 演算法。
1983：超對稱量子力學方法興起。

深化問答三則#

問：Bernoulli 代換 \(v = y^{1-n}\) 是怎麼被想到的？ 答：觀察非線性項 \(q y^n\) 與線性項 \(p y\) 的「冪差」是 \(n - 1\) 次除法可抹平的：除以 \(y^n\) 後，方程中只出現 \(y^{1-n}\) 與 \(y^{-n} y'\) ——後者恰是 \((y^{1-n})'\) 除以 \((1-n)\) 。代換的設計準則是「讓麻煩項與導數項共享同一個新變數」——這條準則也是 Riccati 特解法的內核（ \(\frac1u\) 的導數 \(\frac{-u'}{u^2}\) 提供 \(y^2\) 項的配對）。

問：既然 Liouville 說 Riccati 不可初等解，為什麼還要學它的解法？ 答：三個理由。（1）特例可解（常係數、已知特解、特定係數組合）——實際問題常落在特例中；（2）Riccati 變換把它接駁到二階線性方程，而後者有疊加、級數、Laplace 等全套武器；（3）定性理解（爆破、零點、漸近）不需要顯式解。「解」的意義早已從公式擴展到結構與行為。

問：一個特解「開萬鎖」的機制，能否推廣到高階或偏微分方程？ 答：可以，且形態多變。二階線性方程：一個齊次解 → 降階法（4.1 節）。Painlevé 方程（1900 年代，Painlevé 分類「解無流動奇點」的非線性二階方程）：已知特解未必解全部，但 Bäcklund 變換（1870s 起）可由一個解生成另一個解。PDE 中，Darboux 變換（1882）由 KdV 方程的一個解生成新解，是孤立子理論（1967 GGKM）的引擎。「由已知解繁殖新解」是非線性數學的普遍繁殖策略。

（本節完）