一阶低通滤波器完整推导笔记

### 一阶低通滤波器的数学推导 #### 连续时间域中的传递函数 在一阶低通滤波器的设计中，考虑一个RC电路作为物理模型。该系统的微分方程可以表示为： \[ \tau\frac{dV_{out}(t)}{dt} + V_{out}(t) = V_{in}(t), \] 其中 \( \tau = RC \)，\( R \) 是电阻值，\( C \) 是电容值[^3]。 为了得到频率响应特性，采用拉普拉斯变换方法求解上述微分方程，在零初始条件下可得: \[ H(s)=\frac{V_{out}(s)}{V_{in}(s)}=\frac{\omega_c}{s+\omega_c}, \] 这里引入角频率 \( \omega_c=1/\tau \) 表示截止频率处的时间常数关系。 #### 离散化处理 当转换到离散时间领域时，则需应用双线性变换法将连续形式映射至Z平面内。具体操作如下所示： 设采样周期 T_s ，则有 s 域与 z 域之间的对应关系为 \[ s=\frac{2}{T_s}\cdot\frac{(z-1)}{(z+1)}. \] 代入之前的模拟原型H(s): \[ H(z)=\frac{k_0}{1-1], \] 系数 k_0 和 k_1 可依据实际需求设定，通常取作： \[ k_0=\frac{T_s}{T_s+2R_C},\quad k_1=\frac{-2R_C+T_s}{T_s+2R_C}. \] ```matlab % MATLAB代码实现 function [b,a] = lpfilter(fc,Ts) % fc: 截止频率(Hz); Ts: 采样间隔秒(Sample Time) omega_c = 2*pi*fc; alpha = (Ts*omega_c)/(2+Ts*omega_c); b = [alpha]; % Filter numerator coefficients a = [1, -(1-alpha)];% Filter denominator coefficient end ```