IIR滤波器和FIR滤波器的特点.doc

前言 随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。 数字滤波是数字信号处理的重要内容，数字滤波器可分为IIR和FIR两大类。对于IIR数字滤波器的设计，需要借助模拟原型滤波器，再将模拟滤波器转化为数字滤波器对于FIR数字滤波器的设计，可以根据所给定的频率特性直接设计，采用的设计方法是窗函数法。根据IIR滤波器和FIR滤波器的特点，在MATLAB坏境下分别用双线性变换法设计IIR和用窗函数设计FIR数字滤波器，并对采集的语音信号进行分析，最后给出IIR和FIR对语音滤波的效果。数字滤波器(Digital Filter)是指输入、输出都是离散时间信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数字滤波器在数字信号处理中起着非常重要的作用，在信号的过滤、检测与参数的估计等方面，是使用最为广泛的一种线性系统。实现数字滤波器的方法有两种，一是采用计算机软件进行，就是把所要完成的工作通过程序让计算机来实现；二是设计专用的数字处理硬件。这个地方主要用到的就是第一种方法。即是用Mafiab提供的信号处理工具箱来实现数字滤波器。Matlab信号处理工具箱提供了丰富的设计方法，可以使得繁琐的程序设计简化成函数的调用，只要以正确的指标参数调用函数，就可以正确快捷地得到设计结果。滤波是信号处理的基础，滤波运算是信号处理中的基本运算，滤波器的设计也就相应成为数字信号处理的最基本问题之一。信号带有噪声或无用信号，滤波器的作用是将这些干扰成分滤除，也就是让特定频段的信号通过达到对信号筛选的效果。 式中，称为幅频特性函数；称为相频特性函数。幅频特性表示信号通过该滤波器后各频率成分振幅衰减情况，而相频特性反映各频率成分通过滤波器后在时间上的延时情况。一般选频滤波器的技术要求由幅频特性给出，其相频特性是确定的，所以在设计过程中，对相频特性一般不作要求。在低通滤波器的幅频特性中，和分别称为通带边界频率和阻带截止频率。从到称为过渡带，过渡带上的频响一半是单调下降的。通常，通带内和阻带内允许的衰减一般用分贝表式，通带内允许的最大衰减用表示，阻带内允许的最小衰减用表示。根据对阻带衰减及过渡带的指标要求，选择适当的窗函数设计滤波器。 3.3.5 汉宁窗函数法 汉宁(Hanning)窗，又称升余弦窗 其频率响应和幅度响应分别为： 3.3.6设计程序如下 （1） % 参数设置 f=8*10^3; %信号频率 f1=30*10^3; %噪声频率 l=100; %信号的长度 fs=100*10^3; %采样频率 （2） % 信号的生成 t=(0:l-1)/fs; s=0.5*sin(2*pi*f*t)+cos(2*pi*f1*t); figure(1); subplot(1,2,1); plot(fs*t(1:60),s(1:60)); title(时域波形); xlabel(时间); ylabel(幅度); grid on; nf=2^nextpow2(l); Y=fftshift(fft(s,nf)); %k=-fs/2+(0:nf-1)*fs/nf; k=fs*linspace(0,1,nf)-fs/2; subplot(1,2,2); plot(k,abs(Y)); %信号的频谱图 title(信号频谱); xlabel(频率); ylabel(幅度); grid on; 以上程序会产生时间范围为0-60的有限掺杂噪音的取样信号，信号幅度范围为-1.5——1.5。根据频谱分析，所产生信号会产生频率范围为-5*10^4——5*10^4，幅度为0-43的频谱。 时域波形，信号的频谱图如图3-1： 图3-1时域波形和频谱图 （3） % 滤波器设定 ws=10*10^3; %通带截止频率 ws1=20*10^3; %阻带起始频率 ds=50; %阻带衰减 wt=2*pi*ws/fs; %经采样后的通带截止角频率 wz=2*pi*ws1/fs; %经采样后的阻带截止角频率 wc=(wt+wz)/2; %归一化后的滤波器截止频率 n=ceil(6.6*pi/(wz-wt)); %利用7-3求滤波器长度 t=(n-1)/2; q=fir1(n-1,wc/pi,hanning(n));%滤波器时域函数 w=linspace(0,pi,512); h1=freqz(q,1,512); %进行512个点的傅