iir matlab 带通,基于Matlab的带通IIR数字滤波器设计与仿真

1引言数字滤波技术是数字信号分析、处理技术的重要分支[1]。无论是信号的获取、传输,还是信号的处理和交换都离不开滤波技术,它对信号安全可靠和有效灵活地传输是至关重要的。在所有的电子系统和各类控制系统中,数字滤波器的优劣直接决定产品的优劣。然而,传统的数字滤波器的设计过程复杂,计算工作量大,滤波特性调整困难,影响了它的应用。因此,本文旨在研究一种设计简单、滤波特性调整方便、可移植性强的数字滤波器设计方法。2基于matlab的数字滤波器设计方法鉴于传统的数字滤波器的设计过程复杂,计算工作量大,滤波特性调整困难的缺点,考虑到利用MAT-LAB信号处理工具箱(SignalProcessingToolbox)[2-3]可以快速有效的设计由软件组成的常规数字滤波器;并且,利用MATLAB设计滤波器,可以随时对比设计要求和滤波器特性调整参数,直观简便,极大的减轻了工作量,有利于滤波器设计的最优化,故,本文首先从基于Matlab的数字滤波器设计方法入手。2.1滤波器模型选择设给定的数字带通滤波器的参数是:通带为45Hz~55Hz,低截止频率为40Hz,高截止频率为60Hz,通带内衰减不大于3db,阻带衰减大于80db。借助Matlab仿真,可以得到设计成不同类型的滤波器所需的阶数:表1采用不同类型滤波器实现所需的阶数滤波器类型最低阶数稳定性FIR滤波器Kaiser窗函数201稳定Blackman窗函数180稳定Equirpple86稳定IIR滤波器Butterworth滤波器26稳定Chebishev滤波器16稳定椭圆滤波器12稳定可见,FIR滤波器阶数过高,导致的直接结果是消耗的资源较多,成本增加。因此从成本角度考虑,采用IIR滤波器实现更为合适。从稳定性角度考虑[7-9]:FIR滤波器零极点全部在原点,因而从理论上以及实际的有限精度的运算中,都是稳定的。巴特沃思和切比雪夫滤波器均为全极点网络,椭圆函数滤波器在有限频率上既有零点又有极点,但此三种IIR滤波器都可以通过采用递归结构,使极点在z平面单位圆内以保证其稳定性。虽然递归结构在运算中的四舍五入处理有时会引起寄生振荡,但可以通过增加字长来避免。故IIR滤波器的稳定性也是可以保证的。在IIR滤波器中,巴特沃思和切比雪夫滤波器阻带衰减是频率的单调减函数,而椭圆滤波器幅值响应在通带和阻带都是等波纹的,在同样的性能要求下,比巴特沃思和切比雪夫滤波器所需的阶数低,且过渡带比较窄,可获得更为陡峭的衰减特性曲线。综上考虑,采用椭圆函数滤波器更为适宜。2.2滤波器系数确定用下面的matlab程序可得到满足给定条件的12阶椭圆滤波器的直接型表示:fs=200;wp=[4555]*2/fs;ws=[4060]*2/fs;rp=3;rs=80;Nn=512;[n,wn]=ellipord(wp,ws,rp,rs);[b,a]=ellip(n,rp,rs,wn);freqz(b,a,Nn,fs);[z,p,k]=ellip(n,rp,rs,wn);zplane(z,p);其中,a、b分别为分母与分子的系数,即得到的直接型表示为H(z)=mi=0biz-ini=0aiz-i(a0=1)(1)由于直接型具有一些共同缺点[4-5]:(1)系数对滤波器的性能控制作用不明显。(2)极点对系数的变化过于灵敏,易出现不稳定或较大误差。(3)运算的累计误差较大。并且在这个设计中,直接型表示时系数bi和ai的最大值与最小值相差9104倍,考虑到在乘法器上实现时引入的截断误差很大,因此由直接型传输函数表达式来实现椭圆函数滤波器并不实用,可将其分解为多