目錄
- RC 低通濾波器電路
- 一階線性微分方程推導
- 拉普拉斯域表達(傳遞函數)
- 傳遞函數 H ( s ) H(s) H(s)
- 頻率響應(令 s = j ω s = j\omega s=jω)
- 幅頻特性:
- 相位特性:
- Bode 圖(線性系統頻率響應)
RC 低通濾波器電路
最常見的一階低通濾波器是 RC 電路,如下所示:
Vin ────[ R ]────┬──── Vout│=== C│GND
一階線性微分方程推導
-
根據電路的基爾霍夫定律(KCL)和電容的電壓-電流關系:
- 電容電流: i C = C ? d V out d t i_C = C \cdot \frac{dV_{\text{out}}}{dt} iC?=C?dtdVout??
- 電阻電流: i R = V in ? V out R i_R = \frac{V_{\text{in}} - V_{\text{out}}}{R} iR?=RVin??Vout??
-
因為串聯, i C = i R i_C = i_R iC?=iR?,所以:
C ? d V out d t = V in ? V out R C \cdot \frac{dV_{\text{out}}}{dt} = \frac{V_{\text{in}} - V_{\text{out}}}{R} C?dtdVout??=RVin??Vout??
- 整理為標準一階微分方程y’+p(x)y=q(x)形式:
d V out d t + 1 R C V out = 1 R C V in \frac{dV_{\text{out}}}{dt} + \frac{1}{RC} V_{\text{out}} = \frac{1}{RC} V_{\text{in}} dtdVout??+RC1?Vout?=RC1?Vin?
- 記 τ = R C \tau = RC τ=RC,即:
τ ? d y ( t ) d t + y ( t ) = x ( t ) \tau \cdot \frac{dy(t)}{dt} + y(t) = x(t) τ?dtdy(t)?+y(t)=x(t)
-
其中:
- y ( t ) = V out y(t) = V_{\text{out}} y(t)=Vout?:濾波器輸出;
- x ( t ) = V in x(t) = V_{\text{in}} x(t)=Vin?:濾波器輸入;
- τ = R C \tau = RC τ=RC:時間常數,單位為秒,𝜏表示系統響應變化的快慢程度。
- 注:對于一階系統(例如 RC 電路、熱傳導、液位控制、機械阻尼、速度響應等),系統對一個階躍輸入信號(如電壓從 0 變為 1V)的響應,在經過時間 𝜏 后,輸出達到最終穩態值的約 63.2%。表示系統輸出 y ( t ) y(t) y(t) 跟隨輸入 x ( t ) x(t) x(t),但有慣性,不能瞬間響應。 τ \tau τ 越小 → 響應快, τ \tau τ 越大 → 響應慢,更平滑。
拉普拉斯域表達(傳遞函數)
τ ? d y ( t ) d t + y ( t ) = x ( t ) \tau \cdot \frac{dy(t)}{dt} + y(t) = x(t) τ?dtdy(t)?+y(t)=x(t)
- 對其兩邊做拉普拉斯變換,假設 初始值為零(零初始條件),即 y ( 0 ) = 0 y(0) = 0 y(0)=0,得到:
τ ? L { d y ( t ) d t } + L { y ( t ) } = L { x ( t ) } \tau \cdot \mathcal{L} \left\{ \frac{dy(t)}{dt} \right\} + \mathcal{L}\{y(t)\} = \mathcal{L}\{x(t)\} τ?L{dtdy(t)?}+L{y(t)}=L{x(t)}
- 利用 一階導數的拉普拉斯變換 公式:
τ ( s Y ( s ) ? y ( 0 ) ) + Y ( s ) = X ( s ) \tau (sY(s) - y(0)) + Y(s) = X(s) τ(sY(s)?y(0))+Y(s)=X(s)
- 假設 y ( 0 ) = 0 y(0) = 0 y(0)=0 :
τ s Y ( s ) + Y ( s ) = X ( s ) \tau s Y(s) + Y(s) = X(s) τsY(s)+Y(s)=X(s)
- 提取 Y ( s ) Y(s) Y(s):
Y ( s ) ( τ s + 1 ) = X ( s ) Y(s)(\tau s + 1) = X(s) Y(s)(τs+1)=X(s)
傳遞函數 H ( s ) H(s) H(s)
H ( s ) = Y ( s ) X ( s ) = 1 τ s + 1 \boxed{H(s) = \frac{Y(s)}{X(s)} = \frac{1}{\tau s + 1}} H(s)=X(s)Y(s)?=τs+11??
- 或者換成 RC 表達:
H ( s ) = 1 R C s + 1 \boxed{H(s) = \frac{1}{RC s + 1}} H(s)=RCs+11??
頻率響應(令 s = j ω s = j\omega s=jω)
- 代入 s = j ω s = j\omega s=jω 就描述了系統對頻率為 ω \omega ω 的正弦輸入的穩態響應
H ( j ω ) = 1 1 + j ω τ H(j\omega) = \frac{1}{1 + j\omega \tau} H(jω)=1+jωτ1?
幅頻特性:
∣ H ( j ω ) ∣ = 1 1 + ( ω τ ) 2 |H(j\omega)| = \frac{1}{\sqrt{1 + (\omega \tau)^2}} ∣H(jω)∣=1+(ωτ)2?1?
相位特性:
arg ? H ( j ω ) = ? tan ? ? 1 ( ω τ ) \arg H(j\omega) = -\tan^{-1}(\omega \tau) argH(jω)=?tan?1(ωτ)
Bode 圖(線性系統頻率響應)
-
注:在電路筆記(信號):數字濾波電路的拉普拉斯變換與零極點分析(增益=零點距/極點距,相位=零點角-極點角)的用零極點圖分析電路中提供了一個RC高通的Bode圖估計示例。
-
注:繪制代碼如下
import numpy as np
import matplotlib.pyplot as plt# 參數設置
tau = 0.05 # 時間常數 (s)
w = np.logspace(0, 4, 500) # 頻率向量 (rad/s)
mag = 20 * np.log10(1 / np.sqrt(1 + (w * tau)**2)) # 幅度響應
phase = -np.degrees(np.arctan(w * tau)) # 相位響應# 截止頻率
wc = 1 / tau
fc = wc / (2 * np.pi)# 繪圖
plt.figure(figsize=(10, 6))# 幅度圖
plt.subplot(2, 1, 1)
plt.semilogx(w, mag)
plt.axvline(wc, color='red', linestyle='--', label=f'Cutoff @ {fc:.1f} Hz')
plt.title('Bode Plot of 1st-Order Low-Pass Filter')
plt.ylabel('Magnitude (dB)')
plt.grid(True, which='both')
plt.legend()# 相位圖
plt.subplot(2, 1, 2)
plt.semilogx(w, phase)
plt.axvline(wc, color='red', linestyle='--')
plt.xlabel('Frequency (rad/s)')
plt.ylabel('Phase (degrees)')
plt.grid(True, which='both')plt.tight_layout()
plt.show()