歡迎來到微分方程的世界!
在本章中,我們將跨越簡單代數的範疇,開始探討事物是如何變化的。微分方程 (Differential Equation, DE) 就是包含變化率(導數)的方程式。為什麼這很重要?因為現實世界中的幾乎一切事物——從咖啡如何變涼,到高空彈跳者的擺動——都可以用這些方程式來描述。讀完這些筆記後,你將能夠解開這些難題,並預測系統隨時間演變的規律。如果起初覺得有點棘手,別擔心;我們會將其拆解成簡單、易於掌握的步驟!
1. 一階微分方程:積分因子法
你已經知道如何解像 \(\frac{dy}{dx} = 2x\) 這樣的簡單方程式。但如果方程式長成這樣呢:\(\frac{dy}{dx} + P(x)y = Q(x)\)?這時我們使用一個巧妙的技巧,稱為積分因子 (Integrating Factor)。
檢查「標準形式」
在做任何事之前,你的方程式必須完全符合這個樣子:
\(\frac{dy}{dx} + P(x)y = Q(x)\)
如果 \(\frac{dy}{dx}\) 前面有數字或 \(x\),你必須先將整個方程式除以該項!這是學生最常犯的錯誤。
循序漸進:積分因子法
- 找出積分因子,記作 \(I(x)\)。公式為:\(I(x) = e^{\int P(x) dx}\)。
- 將標準形式方程式中的每一項都乘以這個 \(I(x)\)。
- 方程式的左邊現在會神奇地變成 \((I(x) \cdot y)\) 的導數。你可以將方程式重寫為:\(\frac{d}{dx}(I(x) \cdot y) = I(x) \cdot Q(x)\)。
- 將等式兩邊同時對 \(x\) 積分。
- 整理得到 \(y\)。別忘了你的積分常數 (\(+C\))!
小知識速覽:
如果你的解包含 \(+C\),它被稱為通解 (General Solution)(它代表一族曲線)。如果你得到座標(初始條件)來求出 \(C\) 的特定值,它則被稱為特解 (Particular Solution)。
重點總結:積分因子法實際上是「還原」了乘積法則的微分過程。記住,計算 \(I(x)\) 之前一定要先整理成標準形式!
2. 二階齊次微分方程
現在我們要探討包含二階導數的方程式:\(a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = 0\)。因為等式右邊為零,我們稱之為齊次 (Homogeneous) 方程。
輔助方程式 (Auxiliary Equation)
為了求解,我們「猜測」其解的形式為 \(y = e^{mx}\)。這引導我們得出輔助方程式:
\(am^2 + bm + c = 0\)
你只需像解普通二次方程式那樣求解即可!你得到的根(解)類型決定了解的形式:
- 情況 1:兩個相異實根 (\(m_1\) 和 \(m_2\))
解:\(y = Ae^{m_1x} + Be^{m_2x}\) - 情況 2:重實根 (\(m\))
解:\(y = (A + Bx)e^{mx}\) - 情況 3:複數根 (\(m = \alpha \pm \beta i\))
解:\(y = e^{\alpha x}(A\cos(\beta x) + B\sin(\beta x))\)
你知道嗎?情況 3 發生的情境就是物體產生振盪時,例如吉他弦的振動!
重點總結:解二階齊次微分方程就像解一個二次方程,然後套用正確的答案「模板」一樣簡單。
3. 二階非齊次微分方程
如果方程式右邊不等於零呢?\(a\frac{d^2y}{dx^2} + b\frac{dy}{dx} + cy = f(x)\)。
完整解由兩部分組成:\(y = \text{互補函數 (CF)} + \text{特解 (PI)}\)。
如何求解:
- 求 CF:將方程式設為零,按上一節的方法求解。
- 求 PI:根據 \(f(x)\) 的形式進行「猜測」。
常見的 PI 猜測:
- 如果 \(f(x) = \text{多項式 (例如 } x^2)\),試猜 \(y = px^2 + qx + r\)。
- 如果 \(f(x) = e^{kx}\),試猜 \(y = \lambda e^{kx}\)。
- 如果 \(f(x) = \sin(kx)\) 或 \(\cos(kx)\),試猜 \(y = p\cos(kx) + q\sin(kx)\)。
避免常見錯誤:如果你的 PI 「猜測」已經存在於 CF 中,那是行不通的!你必須將猜測乘以 \(x\)(甚至 \(x^2\))來使其與 CF 區分開來。
重點總結:你可以將 CF 視為系統的「自然」行為,而 PI 則視為對外力的「反應」。
4. 微分方程建模:簡諧運動與阻尼
微分方程是運動學 (Kinematics) 的語言。記住這些術語:
- 位移 (Displacement) = \(x\)
- 速度 (Velocity) = \(\dot{x}\) 或 \(\frac{dx}{dt}\)
- 加速度 (Acceleration) = \(\ddot{x}\) 或 \(\frac{d^2x}{dt^2}\)
簡諧運動 (SHM)
標準的 SHM 方程式為 \(\ddot{x} = -\omega^2x\)。
其解永遠是:\(x = A\cos(\omega t + \phi)\) 或 \(x = p\cos(\omega t) + q\sin(\omega t)\)。
阻尼振盪 (Damped Oscillations)
在現實世界中,摩擦力會減緩運動。我們將其建模為:\(a\ddot{x} + b\dot{x} + cx = 0\)。
輔助方程式的根決定了它減速的方式:
- 過阻尼 (Over-damping):兩個實根。系統緩慢回到平衡狀態而不產生振盪。(例如:門弓器)。
- 臨界阻尼 (Critical damping):重根。回到平衡狀態最快且不產生振盪的方式。(例如:汽車懸吊系統)。
- 欠阻尼 (Under-damping):複數根。系統會來回擺動(振盪),但擺幅會越來越小。
重點總結:你對二次方程式運算的結果,能直接預測一個物理系統會是彈跳、緩慢回歸還是完美停留在平衡點。
5. 聯立一階微分方程
有時兩個變量相互依賴。例如,森林中的兔子數量 (\(x\)) 和狐狸數量 (\(y\))。這些被稱為耦合方程式 (coupled equations)。
如何求解:
- 從兩個方程式開始,例如 \(\frac{dx}{dt} = ax + by\) 和 \(\frac{dy}{dt} = cx + dy\)。
- 將第一個方程式再次微分,得到 \(\frac{d^2x}{dt^2}\) 項。
- 將第二個方程式代入這個新方程式中以「消去」\(y\)。
- 你最終會得到一個關於 \(x\) 的單一二階微分方程,然後就可以使用我們之前學到的方法來求解!
鼓勵的話:這算是本章的「最終 Boss」,但它其實只是把你已經學過的兩個東西結合起來:代入法和二階微分方程!
重點總結:聯立微分方程透過將其轉化為一個二階微分方程來求解。先解出其中一個變量,再利用它求出第二個變量。
總結檢查清單
- 你能識別一階線性微分方程並找到積分因子嗎?
- 你記得輔助方程式的三種情況嗎?
- 你能根據方程式右側選擇合適的特解 (PI) 嗎?
- 你能將二階微分方程的數學與過阻尼、欠阻尼、臨界阻尼聯繫起來嗎?
- 你能透過消去變量來解聯立微分方程嗎?