數學 · 常微分方程 · 01

微分方程Differential equations

把「正在怎麼變」,寫成尋找一個函數的問題。

先認識方程與解,再用切線、候選函數和逐步代回,看清楚什麼才算真正滿足一條變化規則。

系列目錄 · 01/26

同一個高度,斜率也相同嗎?

規則固定為 y′ = y。選一個候選函數,再移動觀察位置。

候選函數實際切線方程要求的方向(虛線)
下方正文提供完整的函數比較與驗解步驟。
點畫布或拉動滑桿;畫布也可用左右方向鍵。座標軸尺度已標示,斜率看數值,不直接量螢幕上的角度。
高度 y(要求的斜率)1
實際導數 y′1
差值 y′ − y0

01 / 從數走到函數

我們要找的,是一整條變化。

一般的方程式問「哪個數字對」;微分方程問「哪一整條曲線,從頭到尾都照這條規則走」。要找的答案是一個函數。

看到 2a + 1 = 7,我們找的是一個數:a = 3。把 3 代回去,左右相等,問題就解決了。

現在換一個問題:有一個量隨著 x 改變,它在每個位置的變化率,都等於當下的值。究竟哪個函數能一直遵守這條規則?

y′(x) = y(x)左邊:函數的變化率。右邊:函數當下的值。

這就是一個微分方程。它把未知函數及其導數放在同一個等式裡,描述兩者必須滿足的關係。這裡未知的是整個 y(x),不只是某個位置上的一個數。

函數告訴你結果

y(x) = eˣ
給定 x,就能算出 y。

微分方程告訴你規則

y′(x) = y(x)
還不知道完整函數,但知道變化率和函數值如何相連。

規則和結果的方向不同:從函數求導,可以檢查它遵守什麼規則;從微分方程出發,則想找出哪些函數能遵守規則。

02 / 符號有各自的位置

高度,和高度正在怎麼變。

y 是現在有多高,y′ 是現在正在變多快。這一頁的方程要求這兩個數字在每一點都剛好相等。

x 是自變數:我們沿著它觀察變化。y(x) 是應變數:它的值依 x 而定。在圖上,x 是橫座標,y 是高度。

y′(x) 讀作「y 的導數」,也可寫成 dy/dx。它描述 x 改變一點時,y 對這個改變的反應;圖像上對應該點的切線斜率。

Δy ≈ y′(x) Δx在函數可微且 Δx 足夠小時,用切線近似局部變化;這裡是近似號。

例如在某一點,若 y′ = 2,x 向右增加很小的 0.01,y 大約增加 0.02。若導數為負,曲線向右走時下降;導數為 0,切線是水平的。

想複習:導數怎麼從平均變化率出現?

先取兩點,相差 h。平均變化率是 [y(x+h)−y(x)]/h,對應兩點連線的斜率。

讓 h 越來越接近 0;若這個比值有極限,就定義為該點的導數:

y′(x) = limh→0 [y(x+h)−y(x)]/h

極限並非把 h 直接代成 0。那會得到沒有意義的除零式。它問的是:h 趨近 0 的過程中,比值趨近哪裡?

因此 y′ = y 的意思很具體:若某點高度為 2,方程就要求斜率為 2;高度為 −1,要求斜率為 −1;高度為 0,要求水平。

斜率是 Δy/Δx,不是高度。這個方程特別要求兩個數值相等,才把高度和斜率連起來。一般函數沒有這個性質。

03 / 逐步代回

eˣ 為什麼是解?

檢查一個函數是不是解,可以先求導數,再代回原式,確認它在所選區間內每一點都成立。

自然指數函數 eˣ 有一個性質:對 x 微分後仍是自己。這是本例用到的微積分結果;e ≈ 2.71828 是固定常數,x 才是變數。

  1. 先提出候選。這一步是指定要檢查的函數,還沒有證明。

    y(x) = eˣ
  2. 計算候選的導數。使用自然指數的微分規則。

    y′(x) = eˣ
  3. 代進微分方程左、右兩邊。

    左邊 y′(x) = eˣ右邊 y(x) = eˣ
  4. 比較。不論 x 取哪個實數,兩邊都相等。

    y′(x) − y(x) = 0

所以 y(x)=eˣ 在整個實數範圍上滿足方程,是它的一個解。這個結論來自代回後的恆等式;拖動滑桿只是把它呈現出來。

「是解」必須連同區間來說

一個解要在考慮的區間內有方程所需的導數,並在每一個點滿足等式。若公式在某點沒定義,就不能假裝它跨過那個點仍然是同一段解。

本例的 eˣ 在所有實數都有定義且可微,所以不用切分區間。後面遇到分母、對數或發散的解時,區間會成為求解的一部分。

04 / 一個容易看錯的地方

碰巧相等,不等於一直遵守。

某一點剛好對,不算是解。方程要的是整段路上每一點都對;找到一個不對的點,就足以否決它。

再看候選 y=x+1。它是一條斜率固定為 1 的直線,因此 y′=1。

y=x+1 在不同位置的驗解
x高度 y導數 y′是否相等
011這一點相等
121不相等

在 x=0,它和 eˣ 恰好有相同的高度與切線。然而往右走,直線的斜率維持 1,方程要求的斜率卻跟著高度變大。它只在一點符合,不是 y′=y 在任何非零長度區間上的解。

y′ − y = 1 − (x+1) = −x只有 x=0 時差值為 0。

回到互動,選「y=x+1」,先停在 x=0,再移向右邊。金色實際切線和藍色要求方向會分開。

兩種看起來很簡單的常數函數

y=1 的導數是 0,所以方程會要求 0=1,不成立。y=0 的導數也是 0,但此時要求的是 0=0,處處成立。

零函數也是一個完整的解。解不一定要上升、振盪,甚至不一定要彎曲。

05 / 方程描述什麼

變化規則,需要寫清楚前提。

同一條式子可以描述很多真實情況,但你得先講清楚每個符號代表什麼、單位是什麼。

如果一個量的變化率和它目前的大小成正比,就可以先用 dy/dt=ky 建模。k 是比例常數;正值、負值與零,描述不同的變化規則。

當 t 有時間單位時,k 必須有「每單位時間」的單位,兩邊才一致。本頁的 y′=y 是無因次的數學例子,相當於比例係數取 1;不是說任何真實系統都照這個速度成長。

建立微分方程、求出解、檢查模型是否貼近現實,是三件相關但不同的工作。數學解完全正確,也可能建立在不適用的物理假設上。

常微分方程和偏微分方程有何不同?

當未知函數只有一個自變數,例如 y(x),使用普通導數的方程稱為常微分方程,英文縮寫 ODE。

若溫度同時依位置與時間變化,函數會有多個自變數;描述不同方向變化的方程可能含偏導數,稱為偏微分方程,縮寫 PDE。本段系列先處理 ODE。

06 / 把想法用一次

先判斷,再展開核對。

先自己判斷,再展開推理,看看哪一步和你的想法不同。

1. y=3eˣ,也是 y′=y 的解嗎?

看完整推理

3 是常數,微分時保留:y′=3eˣ。右邊的 y 也是 3eˣ,每個實數 x 都成立。因此它也是解;方程沒有只允許 eˣ 這一條曲線。

2. y=eˣ+2,還是同一個方程的解嗎?

看完整推理

常數 2 的導數為 0,所以 y′=eˣ。但 y=eˣ+2,兩者始終差 2,不是解。把函數「乘上常數」和「加上常數」是不同操作。

3. 圖上檢查了很多點,都符合,就算證明了嗎?

看完整推理

有限個點只能提供那些點的證據。區間內還有其他點;要證明這個明確候選是解,可以像上面那樣求導、代回,確認得到在整個區間成立的恆等式。只要找到一個不成立的點,就能否決它在包含該點的區間上是解。

現在我們已找到 eˣ、3eˣ 和 0。同一條變化規則,為什麼留下這麼多解?下一頁用「初始條件」把一整族曲線縮成一條。

核對與延伸

來源與範圍

這三頁以一個自變數的常微分方程為主,使用無因次的數學例子。畫面只顯示有限區域;曲線離開畫面,不代表解在那裡終止。

圖像幫助觀察,代回運算才檢查例子是否滿足方程。存在與唯一性的一般理論、更多解析解法,留在後續主題詳細展開。