IIR滤波器设计分析.ppt

无限冲击响应滤波器 IIR 设计 IIR数字滤波器在很多领域中都有着广阔的应用。与FIR数字滤波器相比，IIR数字滤波器可以用较低的阶数获得较高的选择性，而且所用存储单元少。经济效率高。 但相位非线性，且选择性越好，相位非线性越严重。 IIR滤波器的设计方法有：脉冲响应不变法，双线性变换法、MATLAB求系数法。 IIR数字滤波器的结构 任意阶的IIR滤波器可通过数学分解,表示为： 其中为如下的二阶形式： 一般有N≥M，称N为滤波器的阶数 这样就可以将任意阶的IIR滤波器通过若干二阶网络(也称为滤波器的二阶基本节)级联起来构成,其结构如图1所示。其中,代表第i级的二阶网络。 对于每一个二阶基本节,可以转置直接II型结构加以实现,如图2所示。 采用这种级联结构实现IIR滤波器的优点是每一个基本节只是关系到滤波器的某一对极点和一对零点,调整系数a0i、a1i、a2i,只单独地调整了滤波器第i对零点而不影响其他任何零、极点。同样,调整b1i、b2i系数、也只单独调整了第i对极点。因此,这种结构便于准确地实现滤波器的零、极点,也便于调整滤波器的频率响应性能。这种结构的另一个优点是存储单元需要较少,在硬件实现时,可以用一个二阶节进行时分复用,从而降低对FPGA硬件资源的要求。 IIR数字滤波器的设计 IIR数字滤波器的设计实际上是求解滤波器的系数ak和bk。 利用MATLAB信号处理工具箱中的滤波器设计和分析工具(FDATool)可以很方便地设计出符合应用要求的未经量化的IIR滤波器。 示例：IIR低通滤波器设计要求： 低通巴特奥斯滤波器在其通带边缘1kHz处的增益为-3dB，12kHz处的阻带衰减为30dB，采样频率为25kHz。 用MATLAB的FDA Tools求取滤波器系数，步骤: 1)双击Launch pad窗口中的Filter Design Toolbox项中的Filter Design Analysis Tool(FDA Tool) 2)在Design Filter选项卡中选择Filter Type类型为Lowpass 3)在Design Method里选择IIR(Butterworth) 4)在Filter Order中选择Minimum order即最小的阶数 5)Frequency specifications里的Units选择为Hz，输入Fs(采样频率)：25000，Fpass(通带截止频率)：1000，Fstop(阻带截止频率):12000 注意：如果要输入的是归一化频率（2π f/fs），则选择Units为Normalized(0 to 1) 6)输入Apass(通带增益): -3 , Astop(阻带增益)：30 7)设计好各参数后，点击最下面的按钮Design Filter，生成滤波器及系数 8)查看系数：菜单Analysis-Filter Coefficitions 得到的滤波器系数为 Numerator: (分子多项式系数) 0.3346 0.3346 Denominator: (分母多项式系数) 1.0 （注意:该项为a0,实际运算时用不着） -0.3308 还可以直接用MATLAB的命令输入方法求系数： 如本文采用巴特沃斯滤波器，故需要在工具箱中调用的两个函数buttord和butter的调用格式为：[N，Wc]=buttord(wp，ws，Rp，As)[b，a] =butter(N，Wc) 其中N为滤波器阶数：wp和ws分别为通带截止频率矢量和阻带截止频率矢量，单位为π，一般需要模拟频率指标对采样频率的一半作归一化；Rp和As分别为通带最大衰减和阻带最小衰减，单位dB；wc为3 dB边缘频率矢量；b和a即为IIR公式中的系数矢量。 自行根据以上IIR数字滤波器设计方法，基于MATLAB设计一个IIR带通滤波器，其中带通的中心频率为ωp0=0.5π，;通带截止频率ωp1=0.4π,ωp2=0.6π;通带最大衰减αp=3dB;阻带最小衰减αs=15dB;阻带截止频率ωs2=0.3π,ωs2=0.7π (FDA Tool 或者直接命令输入) 脉冲相应不变法设计和双线性变换法 求解滤波器的系数和 ，它是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。 1.用脉冲相应不变法设计IIR数字滤波器 利用模拟滤波器来设计数字滤波器，也就是使数字滤波器能