关键字：IIR滤波器
对于用双线性变换法来设计数字滤波器而言，由低通指标开始，其设计过程如下：
先低通模拟频率→数字频率指标→频率预畸变，得模拟低通指标→进行模拟逼近，求模拟低通的传输函数→双线性变换，求得数字低通的传输函数。
如果想了解双线性Z变换下的高通滤波器设计，可参看楼顺天等编著“基于MATLAB的系统分析与设计——信号处理”3.7节，该书中给出了设计的方法，并提供了一些函数。
实际上双线性Z变换法设计IIR数字滤波器并不复杂，只是有点麻烦，花些功夫还是可以搞清的。它是基于模拟的规正的低通滤波器(截止频率为1)变换而来的，在模拟滤波器设计中要把规正的低通滤波器变成非规正的低通、高通、带通和带阻，也是进行一个变换，它的变换公式是用以下S1的关系式代替传递函数中的S：低通滤波S1/wu高通滤波wl/S1带通滤波(S1^2+wu*w)/(S1*(wu-wl))
带阻滤波S1*(wu-wl)/(S1^2+wu*wl)
例如低通滤波器，用S1/wu=>S.式中wu是低通滤波器的截止频率，或带通/带阻滤波器的高端截止频率；wl是高通滤波器的截止频率，或带通/带阻滤波器的低端截止频率。这样便构成了模拟滤波器的传递函数。
设计双线性数字滤波器首先是要把数字滤波器的频率指标进行卷绕(畸变)，因为在S平面中，w是对应于无穷大，但在数字采样以后，频率只在-ws/2至ws/2的范围内，正负无穷大映射到正负ws/2内，频率轴被扭曲了，数字和模拟的频率轴不是线性关系，而是w1=(2/T)tan(wT/2)，上提到的wu和wl都要按这公式计算成wu1和wl1.
其次是把(2/T)((1-z^(-1))/(1+z^(-1)))替代各类滤波器(低通、高通、带通和带阻)传递函数中的S1，再整理合并，便是z^(-1)的传递函数，它分子和分母的系数便是数字滤波器的系数。 |
