模拟滤波器到数字滤波器的映射方法

模拟滤波器到数字滤波器的映射方法

将已经设计好的模拟滤波器,通过一组映射关系变换为数字滤波器.对于这样一个映射来说,必须确保

1. 映射前后频率对应
2. 因果稳定的系统在映射后仍为因果稳定系统

以下是具体方法:

冲激响应不变法

从时域角度出发,使滤波器得单位冲激响应逼近模拟滤波器的单位冲激响应.为了实现这一目的,我们可以对模拟滤波器的冲激响应$h(t)$进行采样,得到的离散序列即为滤波器的单位冲激响应.
$$h(nT)=h_a(t){|}_{t=nT}$$

系统函数

根据时域推出s域的变换过程,设模拟滤波器的系统函数有N个一阶极点$s_i$,那么其系统函数就可写为
$$H_a(s)=\sum _{k=1}^N\frac{A_k}{s-s_k}$$
变换到时域后,抽样所得离散序列为
$$h(n)=h_a(nT)=\sum _{k=1}^N{A}_k{e}^{s_{k}nT}u(n)$$
等式取z变换求得原模拟滤波器所映射的数字滤波的系统函数
$$H(z)=\sum _{k=1}^N\frac{A_k}{1-e^{s_{k}T}{z}^{-1}}$$
这样看似完成了模拟滤波器的映射,但实际上我们需要注意到一点,或者说我们在设计的时候没有考虑到一点:抽样是否等价.这既是冲激响应不变法的核心,也是所来诸多问题的元凶.

借助于抽样信号频谱与原信号频谱之间关系
$$H(e^{j\omega })=\frac{1}{T}\sum _{k=-\infty }^{\infty }{H}_a(j\frac{\omega -2k\pi }{T})$$
可以看出,抽样信号实际是原信号频谱的周期延拓且加权,所以在实际设计中,令$h(n)=T{h}_a(nT$