微分方程在工程上的應用：從振動分析到電路設計
title: “微分方程在工程上的應用：從振動分析到電路設計” date: 2026-10-07 tags: [數學, 工程, 微積分]
前言：為什麼微分方程如此重要
微分方程（Differential Equation, DE）是描述「變化率」的數學語言。在工程世界中，幾乎所有的物理現象都涉及變化：物體的運動、電流的流動、熱量的傳遞、結構的變形……這些現象的共同特徵是，系統的狀態隨時間或空間的變化率，與系統本身的狀態息息相關。
微分方程正是捕捉這種關係的工具。從牛頓第二運動定律 $F = ma$ 到歐姆定律的電容版本 $i = C\frac{dV}{dt}$，從熱傳導的傅立葉定律到流體力學的納維-斯托克斯方程，微分方程無處不在。
對工程師而言，理解微分方程不只是為了「解題」，更是為了預測系統行為、設計控制器、優化性能。本文將從基礎概念出發，逐步深入，並透過具體的工程實例，展示微分方程如何成為工程師的得力工具。
核心概念解釋：由淺入深
1. 什麼是微分方程？
簡單來說，微分方程是包含未知函數及其導數的方程式。例如：
$$\frac{dy}{dx} = 2x$$
這個方程告訴我們，函數 $y$ 的變化率等於 $2x$。解這個方程，就是找到所有滿足這個條件的函數。
2. 微分方程的分類
微分方程可以依以下方式分類：
- 常微分方程（ODE）：未知函數只有一個自變數，如 $\frac{dy}{dx} + y = 0$
- 偏微分方程（PDE）：未知函數有多個自變數，如 $\frac{\partial u}{\partial t} = \alpha \frac{\partial^2 u}{\partial x^2}$（熱傳導方程）
- 階數（Order）：方程中最高階導數的階數，如 $\frac{d^2y}{dx^2} + 3\frac{dy}{dx} + 2y = 0$ 是二階微分方程
- 線性與否：若未知函數及其導數都是一次的，則為線性微分方程
3. 從淺入深：從一階到高階
一階微分方程是最簡單的形式：
$$\frac{dy}{dx} = f(x, y)$$
例如，放射性衰變遵循 $\frac{dN}{dt} = -\lambda N$，其中 $N$ 是原子數量，$\lambda$ 是衰變常數。
二階微分方程則描述加速度或曲率：
$$\frac{d^2y}{dx^2} + a\frac{dy}{dx} + by = 0$$
這類方程在振動系統、電路分析中極為常見。
4. 初始條件與邊界條件
微分方程的解通常有無限多個，需要額外的條件來確定唯一解：
- 初始條件（Initial Condition）：給定函數在某一點的值，如 $y(0) = 1$
- 邊界條件（Boundary Condition）：給定函數在邊界上的值，如 $y(0) = 0, y(L) = 0$
工程應用實例
應用一：機械振動分析
在機械工程中，振動分析是確保結構安全的關鍵。考慮一個質量-彈簧-阻尼系統：
$$m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = F(t)$$
其中：
- $m$ 是質量
- $c$ 是阻尼係數
- $k$ 是彈簧常數
- $F(t)$ 是外力
- $x(t)$ 是位移
這個二階常微分方程描述了質量塊的運動。當 $F(t) = 0$ 時，我們得到齊次方程，其解取決於阻尼比 $\zeta = \frac{c}{2\sqrt{mk}}$：
- 欠阻尼（$\zeta < 1$）：系統振盪並逐漸衰減
- 臨界阻尼（$\zeta = 1$）：系統最快回到平衡位置而不振盪
- 過阻尼（$\zeta > 1$）：系統緩慢回到平衡位置
工程師利用這個模型來設計汽車懸吊系統、建築物的抗震設計，以及精密儀器的防震平台。
應用二：電路分析
在電子工程中，電容和電感引入了微分方程。考慮一個簡單的 RLC 串聯電路：
$$L\frac{d^2q}{dt^2} + R\frac{dq}{dt} + \frac{1}{C}q = V(t)$$
其中：
- $L$ 是電感
- $R$ 是電阻
- $C$ 是電容
- $q(t)$ 是電荷
- $V(t)$ 是電源電壓
這個方程與機械振動方程形式完全相同，體現了數學的統一性。當 $V(t) = V_0 \sin(\omega t)$ 時，電路會發生諧振，這在無線電通訊中用於選擇特定頻率的信號。
應用三：熱傳導分析
在熱工程中，熱傳導由偏微分方程描述：
$$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2}$$
其中 $T(x,t)$ 是溫度，$\alpha$ 是熱擴散係數。這個方程告訴我們，溫度隨時間的變化率與溫度的空間曲率成正比。
工程師利用這個方程來設計散熱器、優化建築物的隔熱材料，以及預測電子元件的溫度分布。
例題與詳解
例題 1：放射性衰變
問題：某放射性元素的半衰期為 10 年。若初始質量為 100 克，求 30 年後的剩餘質量。
解答：
放射性衰變遵循 $\frac{dN}{dt} = -\lambda N$，其解為：
$$N(t) = N_0 e^{-\lambda t}$$
其中 $N_0 = 100$ 克。
半衰期 $T_{1/2} = 10$ 年，所以：
$$\frac{1}{2} = e^{-\lambda \cdot 10} \Rightarrow \lambda = \frac{\ln 2}{10}$$
30 年後：
$$N(30) = 100 \cdot e^{-\frac{\ln 2}{10} \cdot 30} = 100 \cdot e^{-3\ln 2} = 100 \cdot \left(\frac{1}{2}\right)^3 = 12.5 \text{ 克}$$
例題 2：RLC 電路的諧振
問題：一個 RLC 串聯電路，$L = 0.1$ H，$C = 10$ μF，$R = 5$ Ω。求諧振頻率 $\omega_0$。
解答：
諧振頻率由下式給出：
$$\omega_0 = \frac{1}{\sqrt{LC}} = \frac{1}{\sqrt{0.1 \times 10 \times 10^{-6}}} = \frac{1}{\sqrt{10^{-6}}} = 1000 \text{ rad/s}$$
對應的頻率為：
$$f_0 = \frac{\omega_0}{2\pi} = \frac{1000}{2\pi} \approx 159.15 \text{ Hz}$$
例題 3：熱傳導的穩態解
問題：一根長度為 $L$ 的金屬棒，兩端溫度分別維持在 $T(0) = 100°C$ 和 $T(L) = 0°C$。求穩態溫度分布。
解答：
穩態時 $\frac{\partial T}{\partial t} = 0$，所以：
$$\frac{d^2 T}{dx^2} = 0$$
積分兩次得：
$$T(x) = Ax + B$$
利用邊界條件：
- $T(0) = B = 100$
- $T(L) = AL + 100 = 0 \Rightarrow A = -\frac{100}{L}$
所以：
$$T(x) = 100\left(1 - \frac{x}{L}\right)$$
這是一個線性溫度分布，從 100°C 線性下降到 0°C。
總結與延伸思考
微分方程是工程語言的核心。從最簡單的一階方程到複雜的偏微分方程，它們讓我們能夠：
- 預測系統行為：給定初始條件，我們可以預測系統未來的狀態
- 設計控制器：通過調整參數（如阻尼比），我們可以優化系統響應
- 理解物理本質：微分方程揭示了不同領域之間的深層聯繫
延伸思考：
- 非線性微分方程（如 $\frac{dy}{dt} = y(1-y)$）在生態學、人口動力學中有重要應用
- 數值方法（如歐拉法、龍格-庫塔法）讓我們能夠求解無法解析求解的微分方程
- 現代控制理論中的狀態空間方法，本質上是將高階微分方程轉化為一階方程組
微分方程不僅是數學工具，更是工程師理解世界的窗口。掌握它，就掌握了改變世界的鑰匙。
本文涵蓋了微分方程的基礎概念及其在機械、電子、熱工程等領域應用。希望讀者能從中體會到數學與工程的緊密聯繫，並在實際工作中靈活運用這些工具。