基于双线性变换法的IIR数字高通滤波器设计要点

1、基于双线性变换法的IIR数字高通滤波器设计基于双线性变换法的IIR数字高通滤波器设计摘要随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的 学科和技术领域。在数字信号处理中起着重要的作用并已获得广泛应用的是数字 滤波器(DF, Digital Filter )。数字滤波器是一种用来过滤时间离散信号的数 字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。实现IIR滤波器的阶次较低，所用的存储单元较少，效率高，精度高，而且能够保留一些模拟滤 波器的优良特性，因此应用很广。 Matlab软件以矩阵运算为基础，把计算、可 视化及程序设计有机融合到交互式工作环境中，并且为数字滤波的

2、研究和应用提 供了一个直观、高效、便捷的利器。尤其是 Matlab中的信号处理工具箱使各个 领域的研究人员可以直观方便地进行科学研究与工程应用。本次主要利用matlab的信号处理工具箱设计一个数字高通 IIR滤波器，并 用此滤波器处理一段音频信号。此报告重点介绍了用双线性不变法设计IIR数字 滤波器的基本流程，比较了各种设计方法的优缺点，总结了模拟滤波器的性能特 征。最后以双线性不变法设计了一个高通巴特沃斯IIR数字滤波器，介绍了设计步骤，然后在Matlab环境下进行了仿真与调试，实现了设计目标。关键字IIR数字滤波器 巴特沃斯 高通 双线性变换MATLAB2AbstractWith the

3、information era and the advent of the digital world, digitalsignal processing has become one of the important fields of science and technology. In the digital signal processing plays an important role and has been widely applied in the digital filter . Digital filter is a discrete time signal is use

4、d to filter the digital system,through the mathematical processing of sampling data to reach the goal of frequency domain filterin. Implementation of IIR filter order is low, less storage unit, with high efficiency, high precision, and can keep some simulation characteristics of filter , so a very w

5、ide application. Matlab software based on matrix computation, the calculation,visualization and program design of organic integration to interactive environment, and provides an intuitive, efficient, convenient tool for the research and application of digital filter .The main use of signal processin

6、g toolbox of MATLAB to design a digital high-pass IIR filter, and an audio signal withthe filter . This report focuseson the basic processof using the bilinearinvariant IIRdigital filter design, comparesthe advantages anddisadvantages of various design methods, summarizes thecharacteristics of the a

7、nalog filter. Finally, by bilinear invariant to design a high Butterworth IIR digitalfilter, introducesthe design steps, and then in the Matlab environment for the simulation and debugging, achieves the design goal.Keywords IIR Butterworth high-pass bilinear transform MATLABdigital filter基于双线性变换法的II

8、R数字高通滤波器设计目录摘要 IAbstract II1 .数字滤波器 11.1 数字滤波器介绍 11.2 IIR 数字滤波器设计原理 21.2.2双线性变换法 32 .数字滤波器设计实现 62.1 典型模拟滤波器比较 62.2 设计步骤 62.3 程序流程图 72.4 音频信号部分程序 82.5 仿真结果 82.5.1 滤波器性能仿真 92.5.2 滤波器除噪性能 103 .总结 10致谢 11参考文献 12附 录 13基于双线性变换法的IIR数字高通滤波器设计1 .数字滤波器1.1数字滤波器介绍数字滤波器是具有一定传输选才¥特性的数字信号处理装置,其输入、输出均为数字信号，实质上

9、是一个由有限精度算法实现的线性时不变离散系统。它的基 本工作原理是利用离散系统特性对系统输入信号进行加工和变换，改变输入序列 的频谱或信号波形，让有用频率的信号分量通过，抑制无用的信号分量输出。数字 滤波器和模拟滤波器有着相同的滤波概念，根据其频率响应特性可分为低通、高 通、带通、带阻等类型，与模拟滤波器相比，数字滤波器除了具有数字信号处理的 固有优点外,还有滤波精度高(与系统字长有关)、稳定性好(仅运彳T在0与l两个 电平状态)、灵活性强等优点。时域离散系统的频域特性：Y(ejw )=X(ejw H(ejw ),其中Y(ejw ), X(ejw )分别 是数字滤波器的输出序列和输入序列的频域

10、特性(或称为频谱特性),H(ejw建 数字滤波器的单位取样响应的频谱，又称为数字滤波器的频域响应。输入序列的 频谱X(ejw好过滤波后X(ejw H(ejw )。因此，只要按照输入信号频谱的特点和处 理信号的目的， 适当选择H(ejw ),使得滤波后的X(ejw H(ejw )满足设计的要求， 这就是数字滤波器的滤波原理。数字滤波器根据其冲激响应函数的时域特性， 可分为两种，即无限长冲激响 应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。IIR数字滤波器的特征 是，具有无限持续时间冲激响应，需要用递归模型NN来实现，其差分方程为：y n = = aixn-i二6yn-ii =0i

11、1M brZ“系统函数为：H z =N.二 akz"kg设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数 H(z),使其频 率响应H(z)满足所希望得到的频域指标，即符合给定的通带截止频率、阻带截止 频率、通带衰减系数和阻带衰减系数。1.2 IIR 数字滤波器设计原理IIR数字滤波器是一种离散时间系统，其系统函数为M< bkz"k2Nk1 -v akZ假设Me N,当gN时，系统函数可以看作一个IIR的子系统和一个（M-N）的FIR 子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数 ak和R ,它是数学上的一种逼近问题，即在规定意义上（通常采用最小

12、均方误差准则） 去逼近 系统的特性。如果在S¥面上去逼近，就得到模拟滤波器；如果在z平面上去逼近， 就得到数字滤波器。设计通数字滤波器通常可以归纳为如图所示的两种常用方法。方法1方法2图1-1数字带通滤波器设计的两种方法方法1:首先设计一个模拟原型低通滤波器，然后通过频率变换成所需要的 模拟高通滤波器，最后再使用冲激不变法或双线性变换成相应的数字高通滤波 器。方法2:先设计一个模拟原型低通滤波器，然后采用冲激响应不变法或双线 性变换法将它转换成数字原型低通滤波器，最后通过频率变换把数字原型低通滤 波器变换成所需要的数字高通滤波器。本课程设计采用第一种设计方法，先构造一个巴特沃斯模拟低

13、通滤波器， 然 后将模拟低通滤波器转换成模拟高通滤波器，最后利用双线性变换将模拟高通滤 波器转换成数字高通滤波器。1.2.2双线性变换法为了克服冲激响应法可能产生的频率响应的混叠失真，这是因为从S平面到 Z平面是多值的映射关系所造成的。 为了克服这一缺点，可以采用非线性频率压 缩方法，将整个频率轴上的频率范围压缩到-兀/T兀/T之间，再用z=esT转换 到Z平面上。也就是说，第一步先将整个 S平面压缩映射到S平面的-冗/ T冗 /T 一条横带里；第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性， 也就消除了频谱混

14、叠现象，映射关系如图 1-3j'JS平面15图1-3双线性变换的映射关系为了将S平面的整个虚轴j Q压缩到S1平面j。1轴上的-冗/丁到冗/丁段上,2 tanT可以通过以下的正切变换实现式中，T仍是采样间隔。当Q1由-冗/T经过0变化到冗/T时，。由-oo经过0变化到+oo,也即映射 了整个j。轴。将式（1-9）写成I，1T/2 _er. 1T/2 e eaj.1T/2 . ajiT/2 e e将此关系解析延拓到整个 S平面和S1平面，令j Q=s, j Q1=s1,则得2 e51T/2 eTT/22gT 2 1-e1Ts = sT/2sT72tanh 二-ttT e、eT 2 T 1

15、 e 1T再将S1平面通过以下标准变换关系映射到 Z平面 z=es1T从而得到S平面和Z平面的单值映射关系为:21 -z，T 1 zi Ts_2_T1 s22s T_-2s T这两个关系式是S平面与Z平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换2 1 -eJ s -T 1 eJ=j-tan - = jGT 12J 即S平面的虚轴映射到Z平面的单位圆。其次，将s=o + jw代入，得-j1.1 z:> I'.1T因此由此看出，当仃0时，|z|v1；当仃0时，|z|1。也就是说，S平面的左 半平面映射到Z平面的单位圆内，S平面的右半平面映射到Z平面的单位圆

16、外， S平面的虚轴映射到Z平面的单位圆上。因此，稳定的模拟滤波器经双线性变换 后所得的数字滤波器也一定是稳定的。双线性变换法优缺点：双线性变换法与脉冲响应不变法相比， 其主要的优点 是避免了频率响应的混叠现象。这是因为S平面与Z平面是单值的一一对应关系。 S平面整个j。轴单值地对应于Z平面单位圆一周，即频率轴是单值变换关系。 这个关系重写如下：上式表明，S平面上。与Z平面的成非线性的正切关系，如图1-4所示由图1-4看出，在零频率附近，模拟角频率。与数字频率之间的变换关系 接近于线性关系；但当。进一步增加时，增长得越来越慢，最后当。一时, 终止在折叠频率3 =冗处，因而双线性变换就不会出现由于

17、高频部分超过折叠但是双线性变换的这个特点是靠频率的严重非线性关系而得到的，图1-4所示。由于这种频率之间的非线性变换关系，就产生了新的问题。首先，一个线 性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了；其次，这种非线性关系要求模拟滤波器的幅频响应必须是分段 常数型的，即某一频率段的幅频响应近似等于某一常数 （这正是一般典型的低通、 高通、带通、带阻型滤波器的响应特性），不然变换所产生的数字滤波器幅频响 应相对于原模拟滤波器的幅频响应会有畸变，如图1-5所示。1Ha(jQ)i1 argH (eo图1-5双线性变换法幅应和相位特性的非线性映射argHa(jQ)

18、0H (ej 耳对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波 器，但是各个分段边缘的临界频率点产生了畸变，这种频率的畸变，可以通过频率的预防来加以校正。也就是将临界模拟频率事先加以畸变， 然后经变换后正好 映射到所需要的数字频率上。2.数字滤波器设计实现2.1 典型模拟滤波器比较1、Butterworth巴特沃斯滤波器：它具有单调下降的幅频特性；即最平幅度。2、Chebyshev切比雪夫滤波器：在通带或阻带等波纹，可提高选择性。3、Bessel贝塞尔滤波器：在通带内有较好的线性相位特性。4、Ellipse椭圆滤波器：其选择性相对前三种是最好的。此处选择巴特沃斯主要是想获得

19、最平稳的幅频响应。而不计较相位特性，而用双线性变换法也会将线性相位变为非线性相位。2.2 设计步骤根据以上IIR数字滤波器设计方法，下面运用双线性变换法基于MATLA般计一个IIR高通滤波器。(1)确定性能指标在设计高通滤波器之前，首先根据工程实际的需要确定滤波器的技术指标：通带截止频率fp=3200Hz;阻带截止频率fs=3000Hz;阻带最小衰减 As=100dB和通带 最大衰减Ap=1dB(2)把频率转化为数字角频率wp=2*pi*fp*T;ws=2*pi*fs*T;(3)频率预畸变以及模拟高通性能指标转换成模拟低通性能指标用Q =2/T*tan(w/2)对高通数字滤波器 H(z)的数字

20、边界频率预畸变，得到高通模拟滤波器H(s)的边界频率主要是通带截止频率 Wp阻带截止频率Ws勺转换。抽样频率 Fs=44100Hz。通带截止频率 Wp2=(2/T)*tan(wp/2)阻带截止频率 Ws2=(2/T)*tan(ws/2)(4)归一化处理Wp2 _1 .1p 1 ,/他1，p Wp2p pWs2s = Wp2，1% =一；ssAs110 10 -1N =qlg/lgUs )；2 10 10 -1由式2-3, 2-4, 2-5计算出N,查表可得模拟低通滤波器的阶数，从而由下式 确定模拟高通滤波器的参数H(z)=H( s) | 1 z=H aL (p) | zzs =4 1p -：

21、'c d11，z -1 _ z -H(z)= H (s) 一s =1 +z(5)模拟低通滤波器的构造借助巴特沃斯滤波器模型得到模拟低通滤波器的传输函数Ha(s) oN,Wn=buttord(wp2,ws2,Ap,As,'s');z,p,k=buttap(N);Bap,Aap=zp2tf(z,p,k);G1=tf(Bap,Aap);(6)模拟低通滤波器转换成模拟高通滤波器调用lp2bp函数将模拟低通滤波器转化为模拟高通滤波器。Bbs,Abs=lp2hp(Bap,Aap,Wn);(7)模拟高通滤波器转换成数字高通滤波器利用双线性变换法将模拟高通滤波器Ha(s)转换成数字高通

22、滤波器H(z)Bbz,Abz=bilinear(Bbs,Abs,FS);2.3 程序流程图首先确定性能指标，把频率转化为数字角频率，进而在进行频率预畸变，用 Q =2/T*tan(w/2)对高通数字滤波器 H(z)的数字边界频率预畸变，得到高通模拟 滤波器H(s)的边界频率主要是通带截止频率 Wp;阻带截止频率 Ws的转换。抽样 频率 Fs=44100Hz上述准备工作做好之后，就先把模拟高通性能指标转换成模拟低通性能指 标，然后设计模拟低通滤波器，借助巴特沃斯滤波器模型得到模拟低通滤波器的 传输函数Ha(s)。然后将模拟低通滤波器转化为模拟高通滤波器。最后利用双线 性变换法将模拟高通滤波器 H

23、a(s)转换成数字高通滤波器H(z)。将指标转换成归一化模拟低通滤波器的指标图2-1程序流程图2.4 音频信号部分程序(1)写入声音信号fnose =10;T=length(x)/fs;t=T/length(x):T/length(x):T; j=1:1:length(x);(2)产生噪声并加到声音中nose =sin(2*pi*fnose*T*j);y=x+5*nose'b,a=butter(N,Wn,'high');Y=fft(y)；2.5 仿真结果2.5.1 滤波器性能仿真10Frequency rad/sjw'源程序设计了模拟低通滤波器、模拟高通滤波器与

24、数字高通滤波器，对数字 高通滤波器的性能仿真如下：w' mu n 1一件用胪一门-由揩tl B II I I I I BI I EQ i H三 mu ii 1；=s i 旧照；tl B B I I BI I < ill it;i ii ：i i i;：:ia ii li i i i iii iJ a aais -l n b t ii - _ _ t - r_r rff an! !：! ! ! Wi：II ,III Ifl 1! ; <Nii>!iiiiiiiiE棚i；i Hii ii if -, r-n r r- - r- -i1Ui#三! 持Fil n li I I

25、 hii ： HiHii i a !(. i ； ； ； iinli i » i ； ； ；l tl C I I ；用加 ： ； ： ii； tl B I IIIIIUH到m刖n m ; eh i i i i n11 iiiii工 m I ; tl B I IIIII1 旧旧 J a a i I I I I 1 I I I IIIII!：Z：；：；：於l三 Bi3 iiiii旧的 ；3 E EIII BBII Itl I I IIIi i ： ：1:tiiiiii: j niia i i m-F F ta-« - p - -T-1-B-l-ll h切阊:t iiijlr I

26、 !1 II1 J I I III a j >曲广 n a a ti i i in!；：用:"笑皿髓徵:割蕾:强肘:招; aa-a a.9a raa-a B.9B m aa"a ad-B Aaa a-a-a s Br4ilFT* !4a S4 J 4-a n 4 J I ;ii!： :|/：:II:，i!:"工- r , r iiiri ； ； ；：!：!方语Hriiiiitf1in1iaDid1i a1EFrequeincy (rad/s)in""叫 rr 丁朗 rrr "下洌i ： n ；!： ! ! !：；;： ：!：卜卡

27、卜卜3«卜“卜141 +卜*+卜小HHKI->>I-CH嘴-由缠¥时罚货门安陶厂一二个飞厂一门二拼机 3 )II !tlII tl I I IIlit! BB9 rrrjr-irmrrrsir- r rtti：； -?rr：finr-rj- 11 BUB B J II I I II II : B J !( BUB BII ri ,!II ! )I ! TTT?二二二三中3 一 丁 ： 一 三”若一二二Ti'前二T ” Q丁丁门W河/ «» K* T*K<I> - 7 ; : :iili i ! II J宣浦二rri瀚一zr

28、 t f Tiii ri- - t ' r trrv ?"： i i i ki i ai i: AH国归741期辅W-精箝-!4讣照0-曰由- II I I II II I I II -丁!？!说一厂:丁砾一三污样r十际六?_L.J-.jJJ1 ： : ： IIIJ I I » I I 'i'；：I I I 1 I 1 11I tl I -十QY il I I Li i i j a gIlfFrequency radl/8)2.5.2 滤波器除噪性能为了实现滤波器的应用，程序中加入了有噪声的音频信号，通过对其滤波处 理，来显示数字高通滤波器的功能，

29、下面显示未加入噪声，加入噪声和滤波后的 频谱波形图。滤波前加噪信号的波形和频谱滤波后信号的波形和频谱3.总结这次课设下来，对设计高通数字滤波器的整个过程有了很好的掌握，懂得了设计滤波器的基本方法，对双线性变换法，巴特沃斯滤波器有了一定了解，同时 呢也熟悉了 MATLAB勺环境，巩固了相关知识。初步掌握了 MATLA断言在数字信 号处理中一些基本库函数的调用和编写基本程序等应用；熟悉了滤波器设计的一 般原理，对滤波器有了一个感性的认识；学会了数字高通滤波器设计的一般步骤； 加深了对滤波器设计中产生误差的原因以及双线性变换法优缺点的理解和认识。总之，使理论联系了实际，巩固并深化了对课本基本知识的认

30、识和理解，使理论得以升华。在做本次课程设计的过程中，我深深地感受到了自己所学到知识的有限， 明 白了只学好课本上的知识是不够的，要通过图书馆和互联网等各种渠道来扩充自 己的知识。在实验过程中我们曾经遇到过问题。 但是从中我们学习到了如何对待 遇到的困难，进一步培养了我们一丝不苟的科学态度和不厌其烦的耐心。所有的这些心得会对我以后的学习和工作有帮助作用， 忠心感谢学校给我们提供这次实 验机会。致谢本设计过程中得到老师的悉心指导，设计过程中老师多次帮我分析思路，开拓视角，在我遇到困难想要放弃的时候给予我最大的支持和鼓励。老师那严谨求实的治学态度、踏实坚韧的工作精神，将使我终身受益；另外同学和朋友的

31、热 心帮助也使我学到了不少的知识。再多华丽的语言也显得苍白，在此，谨向老师 和同学们致以诚挚的谢意和崇高的敬意。因为他们我对这次课程设计有了深刻的 印象，虽然只有短暂的一周，但在这期间，却让我受益匪浅，我深刻地认识到了“理论联系实际”这句话的重要性与真实性。参考文献1 程佩青，数字信号处理.北京：清华大学出版社，2007.22 景振毅，张泽兵，董霖.北京：中国铁道出版社，2008.123 从玉良.数字信号处理原理及其 MATLAB实现M.北京：电子工业出版社.2009.74 胡广书.数字信号处理理论、算法与实现M.北京：清华大学出版社.2003,85 刘波.MATLAB1号处理.北京：电子工业出版社，20066 施阳等.MATLA断言工具箱.西安：西北工业大学出版社，19997 丁玉美.数字信号处理.西安电子科技大学出