歡迎來到微分方程的世界！
在你至今的數學旅程中，你已經花了很多時間解 \(x\)。在這章節，我們將要升級了。我們不再只是求出一個數字，而是要解出一整個函數。微分方程 (Differential Equations, DEs) 是包含導數（如 \(\frac{dy}{dx}\)）的方程式。它們非常重要，因為它們描述了事物隨時間變化的方式——從病毒在人口中的傳播，到笨豬跳（bungee jumping）時跳躍者的擺動過程。別擔心如果起初看起來很難；我們會一步一步把它拆解開來！
1. 通解與特解
在我們進入解法之前，我們需要理解什麼是「解」。在進階數學中，我們通常處理兩種類型：
1. 通解 (General Solution)： 這是答案最靈活的版本。由於涉及積分，你的答案會包含任意常數 (arbitrary constants)（如 \(A\)、\(B\) 或 \(c\)）。它代表了一整個可能的曲線「族群」。
2. 特解 (Particular Solution)： 這是特定的答案。如果你得到邊界條件 (boundary conditions)（例如「當 \(x=0\) 時，\(y=5\)」），你可以將這些值代入通解中，求出那些常數的確切數值。
快速回顧：常數陷阱
常見錯誤： 在積分的當下忘記加上積分常數 (\(+c\))。如果你等到題目最後才想著要「補上」，你的代數運算就會出錯！
2. 一階微分方程：積分因子法
當你遇到形如：
\(\frac{dy}{dx} + P(x)y = Q(x)\)
的方程式時，我們使用一種巧妙的技巧，稱為積分因子 (Integrating Factor)。你可以把它想像成一個「魔法乘數」，能將方程式的左邊變成一個乘積法則（product rule）導數的結果。
逐步流程：
1. 標準化： 確保方程式的形式為 \(\frac{dy}{dx} + P(x)y = Q(x)\)。如果 \(\frac{dy}{dx}\) 前面有數字或 \(x\)，請先將整條方程式除以它！
2. 找出積分因子 (\(I\))： 使用公式：\(I = e^{\int P(x) dx}\)。
3. 相乘： 將標準化方程式中的每一項都乘以 \(I\)。
4. 化簡： 左邊現在會自動變成 \((I \times y)\) 的導數。我們將其寫為 \(\frac{d}{dx}(Iy)\)。
5. 積分： 對兩邊進行關於 \(x\) 的積分：\(Iy = \int I \cdot Q(x) dx\)。
6. 求解： 整理得出 \(y = \dots\)
範例：解 \(\frac{dy}{dx} + \frac{1}{x}y = 3\)。這裡 \(P(x) = \frac{1}{x}\)。積分因子為 \(e^{\int \frac{1}{x} dx} = e^{\ln x} = x\)。全部乘以 \(x\) 後得到 \(x\frac{dy}{dx} + y = 3x\)，化簡為 \(\frac{d}{dx}(xy) = 3x\)。
重點提示： 積分因子能「解鎖」方程式。永遠記住 \(e^{\ln f(x)} = f(x)\)；這個恆等式在這裡是你最好的朋友！
3. 二階齊次方程
現在我們來到包含二階導數的方程式：\(y'' + ay' + by = 0\)。這些稱為齊次 (homogeneous) 方程，因為它們等於零。
為了求解這些方程，我們使用輔助方程 (Auxiliary Equation, AE)。我們將導數替換為 \(m\) 的次方：
\(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 i\beta\))
解：\(y = e^{\alpha x}(A \cos \beta x + B \sin \beta x)\)
你知道嗎？ 情況 3 描述的是振盪 (oscillations)，例如擺動的鐘擺或振動的吉他弦！
4. 二階非齊次方程
如果方程式不等於零呢？\(y'' + ay' + by = f(x)\)。
總解由兩個部分相加而成：通解 = 互補函數 (Complementary Function, CF) + 特解 (Particular Integral, PI)。
1. 互補函數 (CF)： 假設方程式等於零（使用上述的 AE 方法）來求解。
2. 特解 (PI)： 觀察「目標」函數 \(f(x)\) 並「猜測」一個試驗解。
試驗解小抄：
- 如果 \(f(x)\) 是多項式（例如 \(x^2 + 3\)）：嘗試 \(y = \lambda x^2 + \mu x + \nu\)。
- 如果 \(f(x)\) 是指數函數（例如 \(e^{5x}\)）：嘗試 \(y = \lambda e^{5x}\)。
- 如果 \(f(x)\) 是三角函數（例如 \(\sin 2x\) 或 \(\cos 2x\)）：嘗試 \(y = \lambda \sin 2x + \mu \cos 2x\)。
重要技巧： 如果你為 PI 做的「猜測」已經包含在你的 CF 中，那麼它將無法運作。你必須將你的猜測乘以 \(x\)（甚至 \(x^2\)）來使它變得獨特！
重點提示： CF 代表系統自然運作的方式，而 PI 代表系統對外力反應的方式。
5. 微分方程建模
在進階數學中，我們使用運動學 (kinematics) 記號。與其使用 \(x\) 和 \(y\)，我們常使用位移 (\(x\)) 和時間 (\(t\))：
- 速度 \(v = \dot{x} = \frac{dx}{dt}\)
- 加速度 \(a = \dot{v} = \ddot{x} = \frac{d^2x}{dt^2}\)
我們應用牛頓第二運動定律 (\(F = ma\)) 來建立方程式。如果力取決於速度（例如空氣阻力），它就會形成一個微分方程！
簡諧運動 (SHM)
一個特殊的情況是 \(\ddot{x} = -\omega^2 x\)。這描述了一個被拉回中心點的物體。其解永遠是一個波：
\(x = A \cos(\omega t) + B \sin(\omega t)\) 或 \(x = R \sin(\omega t + \phi)\)。
阻尼 (Damping)
當你在 SHM 中加入「摩擦力」時，就會得到阻尼振盪 (Damped Oscillations)。這通常建模為 \(a\ddot{x} + b\dot{x} + cx = 0\)。
- 過阻尼 (Overdamping)： 系統緩慢返回平衡狀態而沒有回彈（「摩擦力」太大）。
- 臨界阻尼 (Critical Damping)： 系統以最快速度返回平衡狀態且沒有振盪。（用於汽車懸吊系統！）
- 欠阻尼 (Underdamping)： 系統來回彈跳，但彈跳幅度逐漸減小。
重點提示： 阻尼的類型取決於輔助方程的判別式 (\(b^2 - 4ac\))。
6. 線性系統（耦合方程）
有時你有兩個變數同時在變化，例如狐狸和兔子的數量。這會給你兩個互相「對話」的方程式：
\(\frac{dx}{dt} = ax + by + f(t)\)
\(\frac{dy}{dt} = cx + dy + g(t)\)
如何求解：
1. 消元法： 整理其中一個方程式，使一個變數（比如 \(y\)）成為主項。
2. 微分： 將該新方程式對 \(t\) 微分。
3. 代入： 將這些結果代回另一個方程式。你會得到一個關於單一變數 (\(x\)) 的二階微分方程。
4. 求解： 使用 CF + PI 方法解該二階微分方程，然後利用結果求出第二個變數。
重點提示： 耦合系統總是能塌縮成一個「巨型」方程式。保持你的代數運算整潔！
總結檢查清單
• 我會找出一階微分方程的積分因子嗎？
• 我知道輔助方程根的三種情況嗎？
• 我記得非齊次方程的總解要加上 PI 嗎？
• 我能區分欠阻尼與過阻尼系統嗎？
• 如果 PI 猜測失敗，我記得「乘以 \(x\)」的技巧嗎？
你一定做得到的！微分方程只是一套食譜。一旦你認出了這道「菜」（方程式的類型），只要按照步驟去烹調出答案即可！