一阶RC低通滤波电路数字化

1. 一阶RC电路

1.1 电路图

其中，$R$为滤波电阻，$C$为滤波电容，$U_i$为输入电压，$U_0$为输出电压。

1.2 时域表达式

$$U_0=U_i-RC\frac{d{U}_0}{dt}$$

1.3 传递函数

同时求拉氏变换，在零初始条件下可得：
$$U_0(s)=U_i(s)-RCs{U}_0$$
则传递函数：
$$G(s)=\frac{U_0(s)}{U_i(s)}=\frac{1}{1+RCs}$$

1.4 $Z$变化 & 差分方程

已知传递函数，可利用双线性变换、一阶前向差分和一阶后向差分等方式求取$Z$函数。
(1) 一阶后向差分
$$s=\frac{1-z^{-1}}{T}$$
其中，$T$为离散周期，也就是采样周期。
写成$Z$函数为：
$$Y(Z)=\frac{T}{T+RC}X(Z)+\frac{RC}{T+RC}Y(Z)Z^{-1}$$
写成差分方程为：
$$y(n)=\frac{T}{T+RC}x(n)+\frac{RC}{T+RC}y(n-1)$$
(2) 双线性变换
$$s=\frac{2}{T}\frac{1-Z^{-1}}{1+Z^{-1}}$$
写成$Z$函数为：
$$Y(Z)=\frac{T}{T+2RC}X(Z)+\frac{T}{T+2RC}{Z}^{-1}X(Z)+\frac{2RC-T}{T+2RC}{Z}^{-1}Y(Z)$$

1.5 一阶RC数字滤波器的基本算法

按照一阶差分方程的形式：
$$y(n)=a\ast x(n)+(1-a)\ast y(n-1)$$
其中，$a$是一个和采样周期及RC有关的参数：
$$a=\frac{T}{T+RC}$$
a为小于1的数，在实际可以通过右移的方式实现。
$$\frac{1}{a}=\frac{T+RC}{T}=1+\frac{RC}{T}$$
此时，根据RC电路截止频率公式计算可得：
$$f_L=\frac{1}{2\pi RC}=\frac{a}{2\pi T(1-a)}$$
根据RC电路延迟时间公式计算可得：
$$t_{\tau }=RC=(\frac{1}{a}-1)T$$