跳轉到

1.1 從 Newton 的運動定律到微分方程#

本章導言#

1687 年，Isaac Newton（牛頓，1643–1727） 出版《自然哲學的數學原理》（Philosophiæ Naturalis Principia Mathematica），把宇宙寫成一行字：

\[ F = m a. \]

這行字看似靜態，實則是全書最動態的宣言：力決定加速度，也就是位置的二階導數。因此，知道力就等於知道一個關於位置函數 \(x(t)\) 的方程——

\[ m\,\ddot{x}(t) = F\big(t,\; x(t),\; \dot{x}(t)\big). \]

這正是微分方程（differential equation）：一個聯繫未知函數與其導數的方程。科學從此改變了提問方式——不再是「行星在哪裡？」而是「描述行星位置的函數滿足什麼方程？」本節將沿著歷史軌跡，說明微分方程如何從力學中誕生，並介紹解、初值條件等核心概念。

從幾何到力學：問題的第一次轉向#

十七世紀微積分誕生時，數學家問的是「已知曲線，求切線」。但很快出現了反問題：「已知切線性質，求曲線。」例如：

  • Galileo（伽利略，1564–1642） 在《兩門新科學》（1638）中指出，自由落體的速度與時間成正比： \(v = g t\) 。用導數語言， \(\dot{v} = g\) ——這是最簡單的微分方程。
  • 1638 年，Descartes 與 Fermat 等人討論「法線長等於横坐標某函數」的曲線，實質上是一階微分方程。
  • Huygens（惠更斯，1629–1695） 1659 年研究單擺時，已在不自覺中解二階方程。

真正的轉折是 Newton 的動力學。Newton 的第二定律把「力」翻譯成「導數」，於是一切運動問題都是微分方程問題。Newton 在《原理》中其實避開了微積分記號，用純幾何證明；他的流數法解方程的技巧，直到 1690 年代才在書信與論文中逐漸披露。

「微分方程」一詞的誕生#

術語「equatio differentialis」（微分方程）由 Leibniz（萊布尼茲，1646–1716） 與其學派在 1670–1680 年代的手稿中提出，並由 Johann Bernoulli（約翰・伯努利，1667–1748） 在 1690 年代推廣。里程碑年份如下：

年份 人物 事件
1676 Leibniz 手稿中首次系統研究 \(\dfrac{dy}{dx} = f(x)\) 型方程
1687 Newton 《原理》：以幾何形式處理天體運動方程
1690 Jacob Bernoulli 等時曲線問題，提出「積分方程」式論證
1691 Leibniz 發明變數分離法（separation of variables）
1693 Leibniz 致 Huygens 信中 明確寫出「微分方程」概念
1695 Jacob Bernoulli 提出後世所稱 Bernoulli 方程，1696 由 Johann Bernoulli（經 Leibniz 提示）解出

1690 年 Jacob Bernoulli 解「等時曲線」（isochrone curve：沿此曲線下滑，到達最低點所需時間與起點無關）時，寫出一個積分關係並「倒著微分」得到曲線——這被視為微分方程史上第一個公開求解的著名問題，也開啟了 Bernoulli 家族與 Leibniz 之間既合作又爭鬥的黃金三十年。

什麼是解？#

定義：設 \(n\) 為正整數。 \(n\) 階微分方程是形如

\[ F\big(t,\, y,\, y',\, y'',\, \ldots,\, y^{(n)}\big) = 0 \]

的關係式，其中 \(y = y(t)\) 是未知函數。若函數 \(y(t)\) 在區間 \(I\) 上 \(n\) 階可微，代入後使方程在 \(I\) 上恆成立，則稱 \(y(t)\) 為該方程在 \(I\) 上的解（solution）。

例： \(y'' + y = 0\) 有解 \(y = \sin t\) 、 \(y = \cos t\) ；驗證： \((\sin t)'' = -\sin t\) ，故 \(-\sin t + \sin t = 0\) 。✓

通用解與任意常數#

一階方程的解通常含一個任意常數： \(y' = y\) 的解是 \(y = Ce^t\) 。 \(n\) 階方程的「通解」（general solution）則含 \(n\) 個獨立任意常數：

\[ y'' + y = 0 \quad\Longrightarrow\quad y = A\cos t + B\sin t. \]

這個「階數 = 常數個數」的經驗法則，最早由 Leibniz 學派在 1690 年代察覺，其嚴格證明要到十九世紀（存在唯一性理論，見本書第 3 章）才完成。要從通解中選出一條具體曲線，需要 \(n\) 個初始條件（initial conditions）：

\[ y(t_0) = y_0,\quad y'(t_0) = y_1,\quad \ldots,\quad y^{(n-1)}(t_0) = y_{n-1}. \]

在力學中這完全自然：知道初始位置與初始速度（加上二階方程），運動就唯一確定——這正是 Laplace（拉普拉斯，1749–1827） 1814 年著名的「決定論妖精」論斷的數學基礎：若知宇宙每一質點此刻的位置與速度，未來與過去皆可計算。

從 Newton 第二定律到具體方程：三個原型#

（一）自由落體： \(m\ddot{x} = -mg\) （取向上為正，重力向下），即 \(\ddot{x} = -g\) 。積分兩次得 \(x(t) = x_0 + v_0 t - \frac{1}{2}g t^2\) 。Galileo 由實驗猜到的規律，Newton 理論「算」了出來。

（二）有阻力落體：設阻力 \(\propto\) 速度： \(m\ddot{x} = -mg - k\dot{x}\) 。這是線性方程，可解（見第 2、4 章）。答案揭示「終端速度」（terminal velocity） \(v_t = -mg/k\) ：時間夠久後速度趨於定值——純數學推理預言了跳傘經驗。

（三）單擺：擺角 \(\theta\) 滿足精確方程

\[ \ddot{\theta} + \frac{g}{L}\sin\theta = 0. \]

它是非線性的（ \(\sin\theta\) ）。小幅擺動時以 \(\sin\theta \approx \theta\) 線性化，得 \(\ddot{\theta} + \frac{g}{L}\theta = 0\) ，週期 \(T = 2\pi\sqrt{L/g}\) ——Huygens 1659 年的結果。「線性化」從此成為數學處理世界的核心策略，其功過貫穿全書（見 6.2、7.3 節的教訓）。

Newton 的月地檢驗：一次偉大的方程求解#

1666 年，Newton 由 Kepler 第三定律推算「月球向地球墜落的加速度」應為 \(g/3600\) （因月球距地心約為地表半徑的 60 倍，而反平方律給 \(1/60^2\) ）。他比較月球每秒「墜落」的距離與地表落體每秒墜落 16 英尺，兩者相符。這是史上第一次：一個微分方程（反平方力下的運動）同時解釋蘋果與月亮。Newton 當時並未發表，直到 1679 年與 Edmond Halley（哈雷，1656–1742） 的通信才重拾此線，並在 Halley 力促下於 1687 年出版《原理》。Halley 甚至自費支付出版費——一本改寫人類宇宙觀的書，是靠天文學家的薪水印出來的。

一個具體例題#

例題 1.1.1：驗證 \(y(t) = C_1 e^{2t} + C_2 e^{-2t}\) （ \(C_1, C_2\) 為任意常數）是方程 \(y'' - 4y = 0\) 的解；並求滿足初值條件 \(y(0) = 3\) 、 \(y'(0) = 4\) 的特解。

解：求導兩次：

\[ y' = 2C_1 e^{2t} - 2C_2 e^{-2t}, \qquad y'' = 4C_1 e^{2t} + 4C_2 e^{-2t} = 4\big(C_1 e^{2t} + C_2 e^{-2t}\big) = 4y. \]

故 \(y'' - 4y = 4y - 4y = 0\) ，對任意 \(C_1, C_2\) 恆成立。✓

再解初值問題。 \(t = 0\) 時 \(e^{0} = 1\) ，故

\[ \begin{cases} C_1 + C_2 = 3,\\ 2C_1 - 2C_2 = 4 \;\Rightarrow\; C_1 - C_2 = 2. \end{cases} \]

兩式相加得 \(2C_1 = 5\) ，即 \(C_1 = \frac52\) ， \(C_2 = \frac12\) 。特解為

\[ y = \tfrac52 e^{2t} + \tfrac12 e^{-2t}. \]

歷史評註：這類「指數函數自救」的現象——指數函數求導不變、只乘係數——正是 1740 年代 Euler（歐拉，1707–1783） 發現常係數方程通解的鑰匙（見 4.2 節），而它的力學源頭就是本節的 \(F = ma\) 。

歷史啟示#

  • 提問方式的革命：從「靜態的形」到「動態的律」。Newton 之後，科學的第一語言是微分方程。
  • 理論與觀測的迴圈：Galileo 測量、Kepler 歸納、Newton 統一——微分方程是三者匯流的產物，不是孤立的數學發明。
  • 決定論的世界觀：「通解 + 初始條件」的結構，塑造了十八、十九世紀的機械宇宙觀；而二十世紀混沌理論（7.3 節）恰恰要動搖它「可預測」的另一半。

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

  • 1.2 節：把方程分類（階數、線性與否），是解方程前的第一個動作。
  • 1.3 節：變數分離法是 Leibniz 1691 年的原創技巧，處理最簡單的一類方程。
  • 第 3 章：「階數 = 常數個數」何時成立？答案藏在 Picard 疊代與存在唯一性定理中。
  • 第 4、7 章：單擺的線性化 \(\ddot\theta + \frac{g}{L}\theta = 0\) 是振盪理論的原型。
  • 第 6 章：非線性 \(\sin\theta\) 不能永遠當線性看待——相平面方法還它本來面目。
  • 物理學：Maxwell 方程（1865）、Schrödinger 方程（1926）、Einstein 場方程（1915），全是微分方程；本書的方法是它們的入門磚。
  • 微積分：本書假設讀者已熟悉導數與積分；微分方程是微積分的「反句法」——不是給函數求導，而是由導數還原函數。

練習#

  1. 驗證 \(y = t^2\) 與 \(y = t^2 + C\) 中哪些是 \(y' = 2\sqrt{|y|}\) （ \(t > 0\) ）的解，並指出「非唯一性」現象（此現象在 3.2 節將獲得解釋）。
  2. 由 \(F = ma\) 出發，寫出「質量 \(m\) 的物體掛在彈力常數 \(k\) 的彈簧下，忽略阻力」的運動方程，並指出其階數與線性性。
  3. 驗證 \(y = A\cos t + B\sin t\) 對任意常數 \(A, B\) 都是 \(y'' + y = 0\) 的解，並解初值問題 \(y(0) = 2\) 、 \(y'(0) = -1\) 。
  4. 自由落體方程 \(\ddot{x} = -g\) 有通解 \(x = x_0 + v_0 t - \frac12 g t^2\) 。說明此式如何「內含」Galileo 的三條實驗定律（位移 \(\propto t^2\) 、速度 \(\propto t\) 、與質量無關）。
  5. 思考題：Laplace 的決定論斷言「知初始條件即知未來」。請指出兩個可能使它在實務上失效的理由（提示：一個與 3.2 節的條件有關，一個與 7.3 節的現象有關）。

進階研讀：從運動定律出發的五個深化場景#

場景一：為什麼是二階？ Newton 定律用加速度（二階導）而非速度（一階導）表達，因為力依賴位置與速度但「伽利略相對性」要求慣性運動（ \(v\) 常數）無需力。若定律寫成 \(F = m\dot{x}\) ，慣性運動將需要力，宇宙觀全毀。這說明方程的階數不是數學偏好，而是物理對稱性的要求。

場景二：Newton 為何用幾何寫《原理》？ 當時微積分邏輯地位脆弱（Berkeley 之攻擊尚在後頭），且幾何是公認的嚴格語言。Newton 的《原理》通篇無 \(\dot{x}\) 記號——最偉大的微分方程書，不寫微分方程。這個弔詭提醒我們：記號的缺乏會阻礙發現的傳播；Leibniz 記號在十八世紀席捲歐陸，正因它「可操作」。

場景三：Halley 的問題清單。1684 年 Halley 問 Newton「反平方力下行星軌道為何是橢圓」，Newton 立答「我算過」。但《原理》實際寫作僅花了十八個月，其中還發明了變分思想的雛形。科學史家 Laplace 評：「《原理》是人類心靈最偉大的作品。」讀者可反思：一個好問題（Halley 的提問）對理論誕生的催化作用，不亞於天才本身。

場景四：決定論的邊界。Laplace 之妖需要「無限精確的初始條件」與「無限算力」。Poincaré 在 1890 年研究三體問題時發現：初始條件的微小誤差會指數放大（今稱混沌，7.3 節）。決定論（同樣的初值給同樣的未來）與可預測性（我們能算出未來）是兩回事——微分方程理論花了兩百年才分清這對概念。

場景五：從一個方程到一門學科。1690 年代的「微分方程」只是微積分習題集的一章；到 1750 年代，它已擁有專門技巧（積分因子、級數解）、專門期刊戰場（柏林、彼得堡科學院院刊）與專門人物（Euler 一人貢獻逾半數已知解法）。今日它更分裂為常微分、偏微分、動力系統、隨機微分等分支——本書是這座大廈的門廳。

補充：Newton 與 Leibniz 記號的長期競爭#

Newton 記號 \(\dot{y}\) 、 \(\ddot{y}\) 在英國沿用至十九世紀初，劣勢明顯：只適合以時間為自變數的力學。Leibniz 記號 \(\dfrac{dy}{dt}\) 暗示「可當分數操作」，催生了變數分離（把 \(\frac{dy}{dx}\) 的「分母」搬過去）——這個在嚴格定義下非法、在實踐中屢試不爽的操作，是 1691 年後十年間一大批解法的來源。英國數學因堅守 Newton 記號而與歐陸隔絕半世紀，史稱「數學的英倫冬眠」；1810 年代 Babbage、Peacock 等劍橋分析學會推動改用 Leibniz 記號，才重新接軌。記號不只是衣服，而是器官。

補充：等時曲線問題的骨架#

Jacob Bernoulli 1690 年的問題可化為（今用記號）：求曲線使 \(\displaystyle\int \frac{ds}{\sqrt{y}} = \text{常數}\) 與 \(y\) 無關。他對時間變元微分，得到一階方程後斷言解為擺線（cycloid）。Huygens 曾證明擺線是等時擺的軌道——兩條完全不同的路（幾何的與動力學的）指向同一曲線，令同代人大為震撼，也催生了 1697 年的最速降線問題（Johann Bernoulli 出題，Newton 匿名一夜解出）——變分法自此萌芽，而變分法反過來為微分方程提供 Euler–Lagrange 方程（1744）。一個方程，三代數學。

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

  • 微分方程：聯繫未知函數與其導數的方程。
  • 階數（order）：方程中出現的最高階導數。
  • 解：代入方程使其恆成立的可微函數。
  • 通解：含 \(n\) 個獨立任意常數的 \(n\) 階方程解族。
  • 初始條件：在某一點處指定的 \(y, y', \ldots, y^{(n-1)}\) 之值。
  • 初值問題（IVP）：方程 + 初始條件；第 3 章證明其解何時存在且唯一。
  • 線性化：以導數近似取代非線性項（如 \(\sin\theta \approx \theta\) ）。
  • 終端速度：有阻力落體中速度趨近的常數極限。
  • 流數（fluxion）：Newton 對導數的稱呼。
  • 等時曲線：沿其下滑時間與起點無關的曲線，即擺線。

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

  • 1638：Galileo《兩門新科學》——落體規律；Descartes–Fermat 切線之爭。
  • 1659：Huygens 解等時擺。
  • 1666：Newton 的月地檢驗（瘟疫年間）。
  • 1676：Leibniz 手稿研究一階方程。
  • 1679：Newton–Halley 通信，重拾天體力學。
  • 1687：《原理》出版。
  • 1690：Jacob Bernoulli 等時曲線問題。
  • 1691：Leibniz 變數分離法。
  • 1695–1696：Bernoulli 方程提出與求解。
  • 1744：Euler《尋求具有某種極大或極小性質的曲線的技巧》，變分法問世。
  • 1814：Laplace《概率的哲學隨筆》，決定論之妖。

深化問答三則#

問：微分方程與「函數方程」有何不同？ 答：函數方程（如 \(f(x+y) = f(x)f(y)\) ）聯繫函數值之間的關係；微分方程聯繫函數值與其導數。後者因導數是局部極限，天然攜帶「局部資訊決定整體行為」的性格——這正是存在唯一性定理（3.2 節）的直觀內涵，也是「初值決定未來」的數學根據。

問：為什麼要驗證「解含 \(n\) 個任意常數」而不是直接相信？ 答：因為有反例氛圍：一階方程 \(y'^2 = 4y\) 有含常數的解 \(y = (t - C)^2\) ，也有不含任意常數的奇解（singular solution） \(y = 0\) 。奇解不屬於通解族——十九世紀 Leibniz 學派後繼者與 Lagrange（1774） 已注意到此現象。可見「階數 = 常數個數」需要限制條件，這些限制正是第 3 章的主題。

問：如果 \(F = ma\) 這麼成功，為什麼還需要非線性研究？ 答：因為大多數 \(F\) 本身就是非線性函數：三體引力、流體阻力 \(\propto v^2\) 、擺的 \(\sin\theta\) 。線性方程是特例、是近視；非線性才是常態。二十世紀動力系統（Poincaré 1890 起）與混沌（Lorenz 1963）證明：非線性方程能產生線性世界裡根本不存在的現象——週期自我維持、突變、混沌。非線性不是困難，是寶藏。

（本節完）