脉冲响应不变法设计IIR数字滤波器

《脉冲响应不变法设计IIR数字滤波器》由会员分享，可在线阅读，更多相关《脉冲响应不变法设计IIR数字滤波器（20页珍藏版）》请在装配图网上搜索。

1、第一章摘要之老阳三干创作创作时间：二零二一年六月三十日本设计采纳脉冲响应不变法设计IIR数字低通滤波器,所以在根据脉冲响应不变法设计出来的IIR数字滤波器会很好地重现原模拟滤波器的频率特性；数字滤波器的单元脉冲响应完全模仿模拟滤波器的单元脉冲响应，时域特性迫近好，但容易发生频谱混叠现象，只适合低通和带通滤波器的设计，不适合高通和带阻滤波器的设计.关键词：数字滤波器；脉冲响应不变法；频率混叠第二章引言数字滤波器可以满足滤波器对幅度和相位特性的严格要求，精确度高，有高度的可编程性，灵活性好，而且它没有电压飘移、温度漂移及噪声等问题，基本不受环境影响，稳定性好等.正是由于数字滤波器的以上优点,使得数

2、字滤波器广泛应用于语音处置、图像处置、模式识别、频谱分析、医学仪器等领域.第三章设计原理3.1数字滤波器数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处置来到达频域滤波的目的.可以设计系统的频率响应,让它满足一定的要求,从而对通过该系统的信号的某些特定的频率成份进行过滤,这就是滤波器的基来源根基理.如果系统是一个连续系统,则滤波器称为模拟滤波器.如果系统是一个离散系统,则滤波器称为数字滤波器.数字滤波实质上是一种运算过程,实现对信号的运算处置.输入数字信号（数字序列）通过特定的运算转酿成输出的数字序列,因此,数字滤波器实质上是一个完成特定运算的数字计算过程,也可以理解为

3、是一台计算机.描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则,使其依照这个规则完成对输入数据的处置.时域离散系统的频域特性：其中EeO、|分别是数字滤波器的输出序列和输入序列的频域特性(或称为频谱特性)Hej|是数字滤波器的单元取样响应的频谱，又称为数字滤波器的频域响应.输入序列的频谱区经过滤波后XejHej因此,只要依照输入信号频谱的特点和处置信号的目的，适被选择|HejL使得滤波后的XejHej波原理.,这就是数字滤波器的滤3.2脉冲响应不变法脉冲响应不变法是实现模拟滤波器数字化的一种直观而经常使用的方法.它特别适合于对滤波器的时域特性有一定要求的场所.脉冲响应

4、不变法的核心原理是通过对连续函数ha(t)等间隔采样获得离散序列ha(nT).令h(n)二ha(nT),T为采样间隔.它是一种时域上的转换方法一个模拟滤波器的传递函数可以用有理分式表达式暗示为:M(3-1)H伺a0aiS.aMSNb。bis.b“s通过反拉普拉斯变换我们就可以获得它的冲激相应:iha（t）LiHa（S）（3-2）脉冲响应不变法就是要保证脉冲响应不变，即：h（n）ha（nT）|（3-3）对上式的冲激相应序列阪作Z变换,就可以获得数字滤波器的传递函数：H（z）ZTha（nT）（3-4）一般来说，的分母多项式阶次总是年夜于分子多项式的阶次.假定lHa（S）|的没有多重极点，则式（1）

5、就可分解为：Ha(S)（3-5）式中皿囚均为复数，总是的极点.其拉氏变换为一般来说，亘的分母多项式阶次总是年夜于分子多项式的阶次.数字滤波器的传递函数H（z）经过合并简化，成为一般形式的有理分式传递函数1LC。C|ZClZdodiZ1.d“zN(3-6)在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中，我们已经知道12H(z)STTHa(sjTK)IKI(3-7)依照z=的关系，每一个s平面上宽度为L/T的水平条带将重迭映射到z平面上.因此脉冲响应不变法将s平面映射到z平面,不是一个简单的一一对应的关系.对高采样频率（T小）的情况，数字滤波器在频域可能有极高的增益.为此我们采纳H(z

6、)TAii11eSTz1(3-8)Nh(n)Tha(nT)TAeTnu(n)i1(3-9)在脉冲响应不变法设计中，模拟频率与数字频率之间的转换关系是线性的（丨*T）同时，它可以坚持脉冲响应不变，回=|ha(nT)L因此，这一方法往往用于低通时域数字滤波器设计及相应的模拟系统数字仿真设计.3.3巴特沃斯低通滤波器巴特沃斯滤波器的特点是同频带内的频率响应曲线最为平坦，没有起伏，而在组频带则逐渐下降为零.在振幅的对数对角频率的波特图上，从某一鸿沟见频率开始，振幅随着角频率的增加而逐渐减少，趋向于负无穷年夜.一阶巴特沃斯滤波器的衰减率为每倍频20分贝，二阶巴特沃斯滤波器的衰减率为每倍频12分贝，三阶的

7、衰减率为每分贝18分贝，如此类推，巴特沃斯滤波器的振幅对角频率单调下降，而且滤波器的结束越高，在组频带振幅衰减速度越快，其他滤波器咼阶的振幅对角频率图和低阶数的振幅对角频率有分歧的形状21(S2)N(3-10)上述函数的特点是等距离分布在半径为的圆上.因此,极点用下式暗示为k01,2,N1(3-11)Ha(S)|的暗示式：nHa(S)(sSk)k0(3-12)为了使设计公式和图表统一，将频率归一化.巴特沃斯滤波器采纳3dB截止频率归一化，归一化后的系统函数为cc(3-13)令Rjsc,称为归一化频率,0称为归一化复变量，这样巴特沃斯滤波器的归一化低通原型系统函数为Ga(PPk)k0(3-14)

8、式中，bU,为归一化极点，用下式暗示Pke匕牯|k0,1,2,N1(3-15)第四章设计过程利用模拟滤波器来设计IIR低通数字滤波器是经常使用的方法,称之为模拟一数字转换法.设计步伐禾U用在MATLAB设计IIR数字滤波器可分以下几(1) 步来实现按一定例则将数字滤波器的技术指标转换为(2) 模拟低通滤波器的技术指标；根据转换后的技术指标使用滤波器阶数函数(3) 确定滤波器的最小阶数N和截止频率Wc；利用最小阶数N发生模拟低通滤波原型；利用截止频率Wc把模拟低通滤波器原型转换(4) 成模拟带通原型；利用冲激响应不变法或双线性不变法把模拟滤波器转换成数字滤波器.4.2频率混叠现象数字滤波器的频率

9、响应是模拟滤波器频率响应的周期延拓.只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以内时，才华使数字滤波器的频率响应在折叠频率以内重现模拟滤波器的频率响应，而不发生混叠失真.即j1H(e)THajT|w|n(4-1)可是,任何一个实际的模拟滤波器频率响应都不是严格限带的(非理想)，变换后就会发生周期延拓分量的频谱交叠，即发生频率响应的混叠失真.这时数字滤波器的频响就分歧于原模拟滤波器的频响，而是有一定的失真.当模拟滤波器的频率响应在折叠频率以上处衰减越年夜、越快时，变换后频率响应混叠失真就越小.这时,采纳脉冲响应不变法设计的数字滤波器才华获得良好的效果图4.1脉冲响应不变法的频率混叠现象

10、总结以上，脉冲响应不变法的优点是频率坐标变换是线性的，即W=,如果不考虑频率混叠现象，用这种方法设计的数字滤波器会很好地重现原模拟滤波器的频率特性.另一个优点是数字滤波器的单元脉冲响应完全模仿模拟滤波器的单元脉冲响应,时域特性迫近好其缺点是会发生频谱混叠现象,使数字滤波器的频响偏离模拟滤波器的频响.脉冲响应不变法适合低通、带通滤波器的设计,不适合高通、带阻滤波器的设计第五章法式设计5.1设计巴特沃斯模拟滤波器5.1.1模拟滤波器设计法式clearFs=1000;采样频率为周期倒数Wp=200*pi;Ws=600*pi;设置归一化通带和阻带截止频率Ap=3;As=20;设置通带最年夜和最小衰减N

11、,Wc=buttord（Wp,Ws,Ap,As,s）;%调用butter函数确定巴特沃斯滤波器阶数B,A=butter(N,Wc,s);%调用butter函数设计巴特沃斯滤波器W=(0:pi:1000*pi);%指定一段频率值hs=freqs(B,A,W);%计算模拟滤波器的幅频响应plot(W/pi/2,abs(hs)/abs(hs(1);%绘出巴特沃斯模拟滤波器的扶贫特性曲线gridon;title（巴特沃斯模拟滤波器）;xlabel（频率/Hz）;ylabel（归一化幅值）;5.1.2用分贝显示幅值clearFs=1000;%采样频率为周期倒数Wp=200*pi;Ws=600*pi;%阻

12、带截止频率Ap=3;As=20;%夜和最小衰减N,Wc=buttord(Wp,Ws,Ap,As,butters);函数确定巴特沃斯滤波器阶数s);B,A=butter(N,Wc,数设计巴特沃斯滤波器W=(0:pi:1000*pi);hs=freqs(B,A,W);器的幅频响应hs0=abs(hs)/abs(hs(1)hs1=20*log10(hs0)plot(W/pi/2,hs1);%器的幅频特性曲线gridon;设置归一化通带和设置通带最年%调用调用butter函指定一段频率值计算模拟滤波绘出巴特沃斯模拟滤波采样频率为周设置归一化通带和设置通带最年）;%调用调用butter函title（巴特

13、沃斯模拟滤波器）;xlabel（频率/Hz）;ylabel（分贝/db）;设计巴特沃斯数字滤波器5.2.1数字滤波器设计法式clearFs=1000;%期倒数Wp=200*pi;Ws=600*pi;%阻带截止频率Ap=3;As=20;%夜和最小衰减N,Wc=buttord（Wp,Ws,Ap,As,sbutter函数确定巴特沃斯滤波器阶数B,A=butter（N,Wc,s）;%数设计巴特沃斯滤波器W=(0:pi:1000*pi);%指定一段频率值D,C=impinvar(B,A,Fs);%响应不变法调用脉冲W=(0:0.001*pi:pi)Hz=freqz(D,C,W);%返回频率响应plot(

14、W/pi,abs(Hz)/abs(Hz(1);%绘出巴特沃斯数字低通滤波器的扶贫特性曲线gridon;title（巴特沃斯模拟滤波器）;xlabel（频率/Hz）;ylabel（归一化幅值）;5.2.2用分贝显示幅值clearFs=1000;%采样频率为周期倒数Wp=200*pi;Ws=600*pi;%设置归一化通带和阻带截止频率Ap=3;As=20;%设置通带最年夜和最小衰减N,Wc=buttord(Wp,Ws,Ap,As,s);%调用butter函数确定巴特沃斯滤波器阶数调用butter函B,A=butter(N,Wc,s);%数设计巴特沃斯滤波器指定一段频率值计算出模拟滤W=(0:pi:

15、1000*pi);%hs=freqs(B,A,W);%波器的幅频响应D,C=impinvar(B,A,Fs);%调用脉冲响应不变法W=(0:0.001*pi:pi);Hz=freqz(D,C,W);%返回频率响应Hz0=abs(Hz)/abs(Hz(1);绘出巴特沃斯数字低通滤波Hz1=20*log10(Hz0);plot(W/pi,Hz1);%器的幅频特性曲线gridon;title（巴特沃斯数字滤波器）;xlabel（频率/Hz）;ylabel（分贝/db）;法式中涉及到的函数介绍Buttord功能：用于计算巴特沃斯数字滤波器的阶数N和3dB截止频率wc.N,wc二buttord（wp,w

16、s,ap,as）调用参数wp,ws分别为数字滤波器的通带、阻带截止频率的归一化值，要求：Owwp1,0ws1.1暗示数字频率pi.ap,as分别为通带最年夜衰减和组带最小衰减（dB）.当wswwp时,为高通滤波器;当wp和ws为二元矢量时,为带通或带阻滤波器,这时wc也是二元向量.N,wc作为butter函数的调用参数.N,Qc=buttord(Qp,Qs,ap,as,s)用于计算巴特沃斯模拟滤波器的阶数N和3dB截止频率Qc.Qp,Qs,Qc均为实际模拟角频率.Butter功能：利用butter函数直接设计各种形式的数字滤波器(也可以设计模拟滤波器)B,A=BUTTER(N,Wn),设计一个

17、阶数为n,频率为Wn的低通滤波器;B,A=BUTTER(N,Wn,ftype)可以设计高通,带阻滤波器,其中ftype参数的形式可以指定何种滤波器,ftype为high时，设计一个阶数为n,频率为Wn的高通滤波器;ftype为stop时,获得滤波器阶数为2*n,频率范围为Wn=W1,W2的带阻滤波器；Z,P,K=BUTTER(.),若返回值是三个,则分别是滤波器的零点,极点和增益,其他函数的参数同上面.Plot(X,Y)功能：绘制线性二维图形当X,Y均为实数向量时,而且为维数相同,X=X(i),Y=Y(i),则plot(X,Y)先描述点(X(i),Y(i),然后依次画线;当X,Y均为复数向量时

18、,则不考虑虚数部份;当X,Y均为实数矩阵时,而且为维数相同,plot依次依照对应的列画出线,矩阵有几列就有几条线;当X,Y一个为向量,一个为矩阵时,而且向量的维数即是矩阵的行数或者列数,则把矩阵依照向量的方向分解为几个向量,在与向量配对分别画图,矩阵分解几个向量就有几条线.Gridon2维/3维绘图设置网格线grid功能：来设置以后坐标系的网格线的开/闭.gridon功能：对以后坐标系添加主要网格线.gridoff功能：对以后坐标系移除主要/细节网格线.grid(axes_handle,.)输入：axes_handle坐标系句柄功能：对指定的坐标系设置网格线是否显示.gridminor功能：对

19、以后坐标系设置添加细节网格Impinvar功能：用“脉冲响应不变法”将模拟滤波器离散化BZ,AZ=impinvar(B,A,Fs)把具有B,A模拟滤波器传递函数模型转换成采样频率为Fs(Hz)的数字滤波器的传递函数模型BZ,AZ.采样频率Fs的默认值为Fs=1.Freqz、FreqsFreqz是计算数字滤波器的频率响应的函数Freqs是计算模拟滤波器的频率响应的函数H,F=freqz(B,A,N,Fs)其中B/A提供滤波器系数B为分子A为分母(b0+b1ZA-1+.)/(aO+a1ZA-1+.)N暗示选取单元圆的上半圆等间距的N个点作为频响输出;Fs为采样频率,该参数可以省略H为N个点处的频率

20、响应复值输出向量,其模即为频响幅值曲线幅值20log10(abs(H)DB,其幅角angle(H)即为频响相位曲线相位值.F为与第N点处对应的频率值f(Hz),如果Fs参数省略时,则频率值w为rad/sample,w=2*pi*f/FsH,F=freqz(B,A,N,whole)其参数意义如上,只是调用时角频率的范围由0-pi扩展到了0-2*pi,F返回的是离散系统频率响应H(exp(j*w)在0-2*pi范围内N个频率等分点的对应样值的列向量.freqz(B,A,N),调用时其实不返回系统频率响应的样值,而是直接绘制系统的幅频和相频特性曲线.但要注意,该调用绘制的是对数幅频特性曲线freqz

21、(B,A,N,whole),调用时绘制0-2*pi频率范围内系统对数幅频和相频特性曲线freqs的功能与freqz类似.Filter功能：FILTER是一维数字滤波器Y=FILTER(B,A,X),输入X为滤波前序列,Y为滤波结果序列,B/A提供滤波器系数,B为分子,A为分母整个滤波过程是通过下面差分方程实现的：a(1)*y(n)=b(1)*x(n)+b(2)*x(n-1)+.+b(nb+1)*x(n-nb)-a(2)*y(n-1)-.-a(na+1)*y(n-na)Y,Zf=FILTER(B,A,X,Zi),输入X为滤波前序列,Y为滤波结果序列,B/A提供滤波器系数,B为分子,A为分母，并输

22、入Zi指定X的初始状态,Zf为最终状态矢量第六章运行波形巴特沃斯模拟滤波器运行波形图如下图6.2巴特沃斯模拟低通滤波器幅频响应(分贝)如图6.1、6.2所示,模拟低通滤波器的幅频响应与所给参数基秘闻符.300Hz处的幅值衰减为20db,100Hz处的幅值衰减为1db.通带的指标与要求略有收支,这是因为在法式中计算的时候调用的公式保证了阻带衰减,而通带的指标有富余.数字滤波器的幅频响应图横轴为数字域频率,数字域频率为模拟角频率对采样频率的归一化频率.法式设计中，横轴为数字域频率除以n后的值，所以横轴为0-1.图6.3巴特沃斯数字低通滤波器归一化幅频响应图6.4巴特沃斯数字低通滤波器幅频响应（分贝

23、）如图6.3、6.4所示,数字滤波器的幅频响应与模拟滤波器的幅频响应基本一致.数字滤波器仍然是因果稳定的,数字滤波器的幅频响应模仿模拟滤波器的幅频响应,未发生频谱混叠,到达设计要求.第七章心得体会通过这一周数字信号处置的课程设计,熟悉了MATLAB勺运行环境，了解更多有关于Matlab软件的知识,初步掌握了MATLAB言在数字信号处置中一些基本库函数的调用和编写基本法式等应用；通过对数字低通滤波器的设计让我熟悉了滤波器设计的一般原理,对滤波器有了一个感性的认识；学会了数字低通滤波器设计的一般步伐；加深了对脉冲响应不变法的理解和认识.在设计的过程中,我也认识到了自己所学知识的缺乏.以前上课都是学

24、一些基本的工具,自以为知识掌握的很熟练,实际只是了解了最概况的工具,好多理论只是不能很好的应用于实践,现在运用学到得的工具做出有实际应用价值的工具,对所学知识点进一步的理解,并进行系统化.这也让我再次认识到知识是无尽的只有不竭的充分自己、完善自己的知识理论体系,才华够更好的胜任自己以后的工作.设计过程中知识的缺乏也让我更加坚定了终身学习的决心.第八章参考文献1薛年喜MATLAB在数字信号处置中的应用（第二版）清华年夜学出书社，20082谢平王娜林洪彬信号处置原理及应用机械工业出书社,20093吴湘淇肖熙郝晓莉信号系统与信号处置的软硬件实现电子工业出书社2002年4周浩敏.王睿.测试信号处置技术.北京航空航天年夜学出书社,2005年创作时间：二零二一年六月三十日