基于matlab的iir数字滤波器的设计,基于MATLAB的IIR数字滤波器设计

2015·9(上)《科技传播》 138 科技创新 设计无限长单位脉冲滤波器(IIR)的最通用的方法是借助于现有模拟滤波器的设计方法。模拟滤波器设计的技术和方法应用非常广泛。因此设计数字滤波器可以充分利用这些丰富的资源来进行。 1 利用模拟滤波器来设计IIR带通滤波器的步骤 1)按照一定的变换将给出的数字滤波器的指标转换为模拟带通滤波器的技术指标。 2)根据转换后的技术指标设计模拟低通滤波器 G1(s)(G1(s) 是低通滤波器的传递函数 )。 　 3)将模拟低通滤波器 G1(s) 转换成模拟带通滤波器 G2(s) 4)再按照一定的变换规则将 G2(s) 转换成数字带通滤波器的传递函数。 2 基于 MATLAB 设计 IIR 带通滤波器的步骤设计带通数字滤波器的方法可以归纳为如图所示的方法如图 1 所示。 图1 数字滤波器设计的方法 首先设计一个模拟原型低通滤波器 ，然后通过频率变换成所需要的模拟带通滤波器 ，最后再使用双线性变换成相应的数字带通滤波器。 设计先构造巴特沃斯、切比雪夫 I 和椭圆模拟低通滤波器，然后将模拟低通滤波器转换成模拟带通滤波器，最后利用双线性变换将模拟带通滤波器转换成数字带通滤波器。 在 MATLAB 下设计 IIR 滤波器可使用巴特沃斯函数设计出巴特沃斯滤波器 ，使用 Cheby1 函数设计出契比雪夫 I 型滤波器 ，使用 ellipord 设计椭圆滤波器。IIR 滤波器设计时的阶数不是由设计者所决定 ，而是根据设计者输入的各个滤波器参数由软件设计出满足这些参数的最低滤波器阶数。IIR 滤波器要知道通带上下限截止频率、阻带上下限截止频率、通带纹波和阻带纹波 ，从而确定滤波器的阶数 N。 1)确定数字滤波器的性能指标 ：通带段上限截止频率为 8Hz，通带段下限截止频率为 4Hz，通带纹波为 3dB，下限截止频率为3Hz，上限截止频率为 9Hz，阻带 衰减为 15dB。 2)转换成模拟滤波器的指标 ：由模拟角频率和数字角频率的对应关系确定模拟滤波器的带上下限截止频 率 p  、阻带上下限截止频率 s  、通带纹波 p  和阻带纹波 s  。 3)巴特沃斯滤波器 巴特沃斯低通滤波器可用如下振幅的平方对频率的公式表示 ： (1) 其中 ， = 滤波器的阶数 ； = 截止频率 = 振幅下降为 -3 分贝时的频率 ； = 通频带边缘频率。 根据指标可以得到 lg lg s a N   其中 10 10 10 1 10 1 s p a      (2) 2 1010 1 p C   (3) 根据模拟滤波器的指标可以得出波特沃斯滤波器的阶数和 3dB 频率。 同样可以求出契比雪夫 I 型滤波器和椭圆滤波器的阶数和 3dB 频率并以此得到模拟滤波器的系统函数     12 1 () N Hp p p p p p p     (4) 其中 21 2 , 1,2, , kN Nj k p e k N    4)利用双线性变换法将模拟滤波器映射成数字滤波器。 冲激响应不变法会产生频谱混淆 ，是由于 S 平面到Z 平面的多值对应关系造成的。所以采用双线性换法来将模拟滤波器映射成数字滤波器 ，双线性变换法是 S 平面与 Z 平面一单值对应关系 ，S 平面的虚轴对应于 Z 平面单位圆的一周。 采用非线性频率压缩方法 ，将整个频率轴上的频率范围压缩到横带之间 ，通过标准变换关系再变换到整个Z平面。这样就消除了多值变换性，也就消除了频谱混淆，映射关系如图 2。 基于MATLAB的IIR数字滤波器