![Rank: 6](images/default/star_level3.gif) ![Rank: 6](images/default/star_level2.gif)
- UID
- 864567
|
对于用双线性变换法来设计数字滤波器而言,由低通指标开始,其设计过程如下:
先低通模拟频率--->数字频率指标--->频率预畸变,得模拟低通指标---->进行模拟逼近,求模拟低通的传输函数---->双线性变换,求得数字低通的传输函数.
如果想了解双线性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)的传递函数,它分子和分母的系数便是数字滤波器的系数。 |
|