无限长带通滤波器.doc

1、本 科 毕 业 设 计 (论 文)无限长带通滤波器Infinite bandpass filter 年 月无限长带通滤波器设计摘要本文分别研究了在MATLAB环境下IIR数字滤波器的设计方法及实现方法，并进行图形用户界面设计，以显示本文所介绍滤波器的特性。在无限脉冲响应（IIR）数字滤波器设计中，先进行模拟滤波器的设计，然后进行模拟-数字滤波器转换，即采用脉冲响应不变法及双线性Z变化法设计数字滤波器，最后进行滤波器的频带转换。为方便分析者直观，形象，方便的分析滤波器的特性，创新的设计出了图形用户界面滤波器分析系统。整个系统分为两个界面，其内容主要包含四部分：System(系统)、Analysi

2、s(分析)、Tool(工具)、Help(帮助)。关键词：数字滤波器、MATLAB、无限脉冲响应、有限脉冲响应、图形用户界面Infinite bandpass filterAbstract In this thesis, designs of the Infinite Impulse Response digital filter (IIR) under MATLAB are studied. And the Graphical User Interfaces (GUI) to analysis the characteristics of filter is designed. The des

3、ign of IIR filters can be achieved through three steps: firstly, the design of analog low-pass filter; secondly, it is analog-to-digital filter conversion; lastly, it is the conversion of filter frequency band. In design of FIR filters, We design the Graphical User Interfaces (GUI) of the digital fi

4、lter analysis system makes sure of people can analysis the characteristics of the design directly and easily. The whole system divide into two graphical interfaces, it contains four parts: System, Analysis, Tool and Help. Keywords: Digital Filter, MATLAB, IIR, Graphical User Interfaces 目 录1 绪论12 IIR

5、滤波器设计Y2.1 模拟滤波器设计Y2.2 巴特沃斯滤波器Y2.3切比雪夫滤波器Y3 Y（略）X Y结论 Y致谢 Y参考文献Y附录Y1 绪论1.1 数字滤波技术数字滤波是数字信号处理的重要基础，数字信号处理主要是研究用数字或符号的序列来表示信号波形，并用数字的方式去处理这些序列，把它们改变成在某种意义上更希望的形式，以便估计信号的特征参量，或削弱信号中的多余分量和增强信号中的有用分量。数字滤波器在对信号的过滤、检测与参数估计等处理过程中，是使用最为广泛的一种线性系统。数字滤波器（Digital Filter，简称为DF）是指完成信号滤波处理功能的、用有限精度算法实现的离散时间线性非时变系统。数

6、字滤波器的数学运算通常有两种实现方式。一种是频域法，即利用FFT快速运算办法对输入信号进行离散傅立叶变换，分析其频谱，然后根据希望的频率特性进行滤波，再利用傅立叶反变换恢复出时域信号。这种方法具有较好的频域选择特性和灵活性，并且由于信号频率与所希望的频谱特性是简单的相乘关系，所以它比计算等价的时域卷积要快得多。另一种方法是时域法，这种方法是通过对离散抽样数据作差分数学运算来达到滤波目的的。数字滤波器的输入是一组（由模拟信号取样和量化的）数字量，其输出是经过数字变换的另一组数字量。数字滤波器具有稳定性高、精度高、灵活性大等突出优点。随着数字技术的发展，用数字技术实现滤波器的功能愈来愈受到人们的重

7、视，并得到了广泛的应用。数字信号处理学科的一项重大进展是关于数字滤波器设计方法的研究。60年代中期，开始形成关于数字滤波器的一套完整的正规理论。这一时期，提出了各种各样的数字滤波器结构，有的以运算误差最小为特点，有的则以运算速度高见长；出现了数字滤波器的各种逼近方法和实现方法，对递归和非递归两类滤波器作了全面的比较；统一了数字滤波器的基本概念和理论，对有限冲激响应（IIR）和无限冲激响应（FIR）的认识有了完整理论。70年代后，科学技术蓬勃发展，数字信号处理开始与大规模和超大规模集成电路技术、微处理技术等新工艺新技术结合起来，并引进计算机辅助设计方法，大大丰富了数字滤波器的分析与设计，各种新的

8、数字信号处理系统，也都能用专用数字硬件实时加以实现。相信在未来，随着电子仪器与电子技术应用系统朝着数字化、小型化、自动化以及多功能化等方向发展，包括数字滤波器在内的数字信号处理技术会有以惊人的速度进行飞跃式发展。1.1.1 滤波器原理滤波器，顾名思义，其作用是对输入信号起到滤波的作用。对于图1-1所示的LSI系统，其时域输入输出关系为： (1-1)h (n) x(n) y(n) 图1-1若y(n),x(n)的傅立叶变换存在，则输入输出的频域关系是： (1-2)当输入信号x(n) 通过滤波器系统h(n) 后，其输出 y(n) 中不再含有|c的频率成分，仅使|c的信号成分通过。因此，滤波器的形状不

9、同，其滤波后的信号结果也不一样。若滤波器的输入、输出都是离散时间信号，那么该滤波器的单位冲激响应h(n)也必然是离散的，这种滤波器称为数字滤波器（DF，Digital Filter）。当用硬件实现一个DF时，所需元件是延迟器、乘法器和加法器；而利用MATLAB软件时，它仅需线性卷积程序便可实现。而模拟滤波器（AF，Analog Filter）只能用硬件实现。因此DF比AF容易实现，且更容易获得理想的滤波性能。数字滤波器还具有以下优点：精度和稳定性高；系统函数容易改变；灵活性高；不存在阻抗匹配问题；便于大规模集成；可实现多维滤波。1.1.2 数字滤波器设计方法概述数字滤波器从功能上分类：可以分为

10、低通滤波器(LP，Low Pass)、高通滤波器(HP，High Pass)、带通滤波器(BP，Band Pass)、带阻滤波器(BS，Band Stop)。从滤波器的网络结构或者从单位脉冲响应分类：如同模拟滤波器的性能可由g(t)和G(s)来表征一样，数字滤波器的性能完全取决于h(n)和H(z)。因此，数字滤波器可以按照单位取样响应（或称脉冲响应，冲激响应等）h(n)的性质分为两类：有限脉冲响应(Finite Impulse Response)数字滤波器，简称FIR数字滤波器，它的h(n)序列长度是有限的；无限脉冲响应(Infinite Impulse Response)数字滤波器，简称II

11、R数字滤波器，它的h(n)序列长度是无限的，即当时，h(n)仍有效。数字滤波器的设计与实现，通常按下述步骤进行：1、根据不同用途提出数字滤波器的技术指标、性能要求。2、设计一个稳定的、因果的数学模型H(z)来逼近所要求的技术指标，并用有限精度的运算实现所设计的系统。本文将重点介绍此步骤。3、设计专用的数字硬件来实现这个数字模型，或者通过电子计算机运行程序软件予以实现。本文将侧重于软件实现。IIR滤波器和FIR滤波器的设计方法很不相同：IIR滤波器设计方法有两类，经常用到的一类设计方法是借助于模拟滤波器的设计方法进行的。其设计思路是：先设计模拟滤波器得到传输函数G(s)，然后将G(s)按某种方法

12、转换为数字滤波器的系统函数H(z)。这一类方法是基于模拟滤波器的设计方法相对比较成熟，它不仅有完整的设计公式，也有完整的图标供查阅，更可以直接调用MATLAB中的对应的函数进行设计。另一种是直接在频域或者时域中进行设计，设计时必须使用计算机辅助，直接调用MATLAB中的程序或函数即可设计。FIR滤波器不能采用由模拟滤波器设计进行转换的方法，而经常使用的是窗函数法和频率采样法。也可以借助计算机辅助设计软件采用切比雪夫等波纹逼近法进行设计。1.2 MATLAB软件简介MATLAB是英文Matrix Laboratory（矩阵实验室）的缩写，它是由美国MathWorks公司推出的用于数值计算和图形处

13、理的数学计算环境。在MATLAB环境下，用户可以集成的进行程序设计、数值计算、图形绘制、输入输出、文件管理等各项操作。它优秀的数值计算能力和卓越的数据可视化能力使其在同类软件中脱颖而出。MATLAB系统最初由Cleve Moler用FORTRAN语言设计，现在的MATLAB程序是由MathWorks公司用C语言开发的。它的第一版（DOS版本1.0）发行于1984年，经过多年的改进，版本不断升级，其所包含的工具箱功能也越来越丰富，应用越广泛。MATLAB语言之所以能如此迅速的普及，显示出如此旺盛的生命力，是由于它有着不同于其他语言的特点，归纳如下：1、简单易学：MATLAB不仅是一个开发软件，也

14、是一门编程语言。其语法规则与结构化高级编程语言(如C语言等)大同小异，而且使用更为简便。用MATLAB编程写程序犹如在演算纸上排列出公式与求解问题，所以被称为演算纸式科学算法语言。2、计算功能强大：MATLAB拥有庞大的数学、统计及工程函数，可使用户立刻实现所需的强大数学计算功能。由各领域的专家学者们开发的数值计算程序，使用了安全、成熟、可靠的算法，从而保证了最快的运算速度和可靠的结果。此外，MATLAB还有数十个工具箱，可解决应用中的大多数数学、工程问题。3、先进的可视化工具：MATLAB提供功能强大的、交互式的二维和三维绘图功能，可使用户创建富有表现力的彩色图形。可视化工具包括曲面渲染、线

15、框图、伪彩图、光源、图像显示、动画等。4、开放性、可扩展性强：M文件是可见的MATLAB程序，所以用户可以查看源代码。开放的系统设计使用户能够检查算法的正确性，修改已存在的函数，或者加入自己的新部件。5、特殊应用工具箱：MATLAB的工具箱加强了对工程及科学中特殊应用的支持。工具箱也和MATLAB一样是完全用户化的，可扩展性强。将某个或几个工具箱与MATLAB联合使用，可以得到一个功能强大的计算组合包，满足用户的特殊要求。MATLAB数字信号处理工具箱和滤波器设计工具箱专门应用于信号处理领域。工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标参数调用相应的滤波

16、器设计程序或工具箱函数，便可以得到正确的设计结果，使用非常方便。接下来我们将分别对IIR数字滤波器的设计实现进行分析讨论。在具体的设计过程中，将充分发挥MATLAB软件在数值计算、图像处理中的优势，并充分利用其强大的滤波器设计工具箱(Filter Design Toolbox)及图像用户界面工具辅助设计分析。2 IIR滤波器设计及其MATLAB实现IR数字滤波器具有无限持续时间脉冲响应，而所有的模拟滤波器一般都具有无限长的脉冲响应，因此它与模拟滤波器相匹配。于是，设计IIR数字滤波器可以先设计模拟低通滤波器G(s)，然后按一定规则将G(s)进行滤波器变换而得到数字域的IIR低通滤波器H(z)，

17、最后通过频带变换转换成其它频带滤波器。本章按按此思路将IIR数字滤波器设计分为三部分：设计模拟滤波器；进行滤波器转换，得到数字滤波器；最后再简单分析数字滤波器的频带转换问题。下面，将以此进行分别讨论。2.1 模拟滤波器设计及其MATLAB实现IIR数字滤波器技术设计最通用的方法是借助于模拟滤波器的设计方法。模拟滤波器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图表可供查询，因此充分利用这些已有的资源将会给数字滤波器的设计带来很大方便。滤波器设计最重要的是寻找一个稳定、因果的系统函数去逼近滤波器的技术指标，因此模拟滤波器的设计十分重要。实际中，有三种广泛应用的滤波器

18、，即巴特沃斯低通滤波器、切比雪夫低通滤波器，椭圆低通滤波器，下面，将分别予以介绍。2.1.1 巴特沃斯低通滤波器设计巴特沃斯滤波器拥有最平滑的频率响应，在截止频率以外，频率响应单调下降。在通带中是理想的单位响应，在阻带中响应为0。在截止频率处有3dB的衰减。巴特沃斯滤波器除了具有平滑单调递减的频率响应的优点之外，其过渡带的陡峭程度正比于滤波器的阶数。高阶巴特沃斯滤波器的频率响应近似于理想低通滤波器。巴特沃斯(Butterworth)低通滤波器是将巴特沃斯函数作为滤波器的传递函数，它的平方幅度函数为： (2-1)式中,为滤波器频率，为3dB截止频率，N表示滤波器的阶次。MATLAB信号处理工具箱

19、为低通模拟巴特沃斯滤波器的产生提供了函数buttap,其调用的格式为：z,p,k=buttap(N)，其中，z表示零点,p表示极点,k表示增益，N表示阶次。下面举例设计巴特沃斯低通模拟滤波器：%巴特沃斯低通模拟滤波器clear all;n=0:0.01:2;for i=1:4switch i case 1 N=2;case 2 N=5; case 3 N=10; case 4 N=20; endz,p,k=buttap(N); %函数buttap-设计巴特沃斯低通滤波器b,a=zp2tf(z,p,k); %函数zp2tf-零极点增益模型转换为传递函数模型H,w=freqs(b,a,n); %函

20、数freqs-求解模拟滤波器频率响应magH2=(abs(H).2; %函数abs-取模值函数 hold on %函数hold-控制是否保持当前图形plot(w,magH2) %函数plot-画二维线性图axis(0 2 0 1); %函数axis-控制坐标轴比例和外观endxlabel(w/wc);ylabel(|H(jw)|2);title(巴特沃斯低通模拟滤波器);grid on由此得出的巴特沃斯低通模拟滤波器的平方幅度响应曲线如图所示： 图2-1 巴特沃斯滤波器的幅频特性(基于MATLAB实现)由图2-1可清楚地分析出，巴特沃斯滤波器拥有平滑的频率响应，在截止频率以外，频率响应单调下降

21、。其过渡带的陡峭程度正比于滤波器的阶数，高阶巴特沃斯滤波器的频率响应近似于理想低通滤波器。由图2-1所示，对N=2阶滤波器和N=20阶滤波器的平方幅度响应进行比较后，证明了高阶巴特沃斯滤波器有着更好的幅度特性，更接近理想低通滤波器。2.2 模拟-数字滤波器变换及其MATLAB实现在设计了模拟低通滤波器后，就可以把它们变成数字滤波器了。这些变换均是复值映射，许多文献对此都有研究，根据数字滤波器所保持的模拟滤波器的不同特性，研究出不同的变换技术。其中，最重要的有两种：脉冲响应不变法（保持脉冲响应不变，又叫冲激响应不变法）和双线性Z变换法（保持系统函数不变）。2.2.1 脉冲响应不变法脉冲响应不变法

22、的设计原理是使数字滤波器的单位抽样响应序列h(n)，模仿模拟滤波器的脉冲响应g(t)。设系统传递函数为G(s)的模拟滤波器的单位脉冲响应为g(t),并将脉冲响应g(t)进行等间隔采样，使得数字滤波器的单位抽样响应h(n)刚好等于g(t)的采样值，即： (2-6)其中的为采样周期。G(s)是模拟滤波器的系统传递函数，又令H(z)是数字滤波器的系统传递函数。采样信号的拉式变换与相应的采样序列Z变换的映射关系为： (2-7)所以系统函数G(s)和H(z)的关系为： (2-8)式(2-8)的物理意义为首先将模拟滤波器的系统函数G(s)作周期的延拓，在经过式(2-7)的映射变换，映射到Z平面上，从而得到

23、数字滤波器的系统函数H(z)。且模拟和数字频率满足下列关系：=T。经过式(2-7)的映射，s平面的左半平面映射为Z平面的单位圆内，因此，一个因果的和稳定的模拟滤波器映射成因果的和稳定的数字滤波器。经过以上分析，按照脉冲响应不变法，通过模拟滤波器的系统传递函数G(s)，可直接求得数字滤波器系统函数H(Z),其设计具体步骤归纳如下：(1)利用=T（可由关系式推出），将数字滤波器指标，转换为模拟滤波器指标，(2)根据指标，来设计模拟滤波器G(s)(3)利用部分分式展开法，把G(s)展成 (2-9)(4)最后把模拟极点转换为数字极点，得到数字滤波器： (2-10)根据上述理论，将举例在MATLAB环境

24、下用函数实现脉冲响应不变法设计一数字低通滤波器。其函数为b,a=impinvar(c,d,T)，其中，b表示数字滤波器自变量为的分子多项式，a表示数字滤波器自变量为的分母多项式，c表示模拟滤波器自变量为s的分子多项式,d表示模拟滤波器自变量为s的分母多项式，T表示采样变换参数。%设模拟低通巴特沃斯滤波器，通带纹波为Rp=1dB,通带上限角频率p=0.2,阻带下限角频率s=0.3，带阻最小衰减s=15dB，根据该低通模拟滤波器，利用冲激响应不变法设计相应的数字低通滤波器wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;T=1; %采样周期%性能指标Rip=10(-Rp/20);Atn

25、=10(-As/20);OmgP=wp*T;OmgS=ws*T;N,OmgC=buttord(OmgP,OmgS,Rp,As,s); %选取模拟滤波器的阶数cs,ds=butter(N,OmgC,s); %设计出所需的模拟低通滤波器b,a=impinvar(cs,ds,T); %应用脉冲响应不变法进行转换 db,mag,pha,grd,w=freqz_m(b,a); %求得相对，绝对频响及相位、群延迟响应%绘制各条曲线subplot(2,2,1);plot(w/pi,mag);title(幅频特性);xlabel(w(/pi);ylabel(|H(jw)|);axis(0,1,0,1.1);s

26、et(gca,XTickMode,manual,XTick,0 0.2 0.3 0.5 1);set(gca,YTickMode,manual,YTick,0 Atn Rip 1);gridsubplot(2,2,2);plot(w/pi,db);title(幅频特性(dB);xlabel(w(/pi);ylabel(dB);axis(0,1,-40,5);set(gca,XTickMode,manual,XTick,0 0.2 0.3 0.5 1);set(gca,YTickMode,manual,YTick,-40 -As -Rp 0);gridsubplot(2,2,3);plot(w/

27、pi,pha/pi);title(相频特性);xlabel(w(/pi);ylabel(pha(/pi);axis(0,1,-1,1);set(gca,XTickMode,manual,XTick,0 0.2 0.3 0.5 1);gridsubplot(2,2,4);plot(w/pi,grd);title(群延迟);xlabel(w(/pi);ylabel(Sample);axis(0,1,0,12);set(gca,XTickMode,manual,XTick,0 0.2 0.3 0.5 1);grid%本设计用到的M文件函数functiondb,mag,pha,grd,w=freqz_

28、m(b,a)%滤波器幅值响应（绝对、相对）、相位响应及群延迟%Usage:db,mag,pha,grd,w=freqz_m(b,a) %500点对应0,pi%db 相对幅值响应； mag 绝对幅值响应； pha 相位响应； grd 群延迟响应；%w 采样频率； b 系统函数H(z)的分子项%a 系统函数H(z)的分母项H,w=freqz(b,a,500); %500点的幅频响应mag=abs(H);db=20*log10(mag/max(mag);pha=angle(H);grd=grpdelay(b,a,w);图2-7 脉冲响应不变法设计数字低通滤波器的特性(基于MATLAB实现)总结以上，

29、脉冲响应不变法的优点是频率坐标变换是线性的，即=T，如不考虑频率混叠现象，用这种方法设计数字滤波器会很好的重现原模拟滤波器的频率响应。另外一个优点是数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，时域逼近好。但其也具有很大的缺点，若抽样频率不高或其它原因将产生混叠失真，不能重现原模拟滤波器频率响应。所以，脉冲响应不变法适合低通、带通滤波器设计，不适合高通、带阻滤波器的设计。2.2.2 双线性Z变换法利用脉冲响应不变法设计数字滤波器时，由于=T的频率关系是根据推导的，所以是轴每隔2/T便映射到单位圆上一周，引起了频域混叠。为克服这一现象，人们找到了另一种映射关系： (2-11)此关系称

30、为双线性Z变换法。双线性Z变换法的基本思路是：首先将整个s平面压缩到s1平面的一条带宽为2/T（丛-/T到/T）的横带里，然后通过标准的变化关系将横带变换成整个Z平面上去，这样就得到s平面与Z平面间的一一对应的单值关系，整个过程如图2-8所示： j j1 jIm(Z) /T 0 0 0 1 -/T s平面 s1平面 Z平面图2-8 双线性Z变换法的映射关系由式(2-11)得 (2-12)及 (2-13) (2-14)式(2-11)及式(2-12)给出了s和z之间的映射关系，而式(2-13) 和式(2-14)给出了和之间的映射关系，但这是一种非映射关系，双线性Z变换法正是利用了正切函数的非线性特

31、点，把整个j轴压缩到了单位圆的一周上。在MATLAB中，双线性Z变换可通过bilinear函数实现,其调用格式为：Bz,Az=bilinear(B,A,Fs)，其中B,A为模拟滤波器的传递函数G(s)的分子分母多项式系数分量，而Bz,Az为数字滤波器的传递函数H(Z)的分子分母多项式的系数分量。%使用双线性Z变换设计一低通数字滤波器，fp=100HZ,fs=300HZ，p=3dB,s=20dB,抽样频率为Fs=1000HZfp=100;fs=300;Fs=1000;rp=3;rs=20; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Fs=Fs/Fs %使Fs为一wap=tan(wp

32、/2);was=tan(ws/2);n,wn=buttord(wap,was,rp,rs,s)z,p,k=buttap(n);bp,ap=zp2tf(z,p,k)bs,as=lp2lp(bp,ap,wap)bz,az=bilinear(bs,as,Fs/2)H,w=freqz(bz,az,256,Fs*1000);图2-9 双线性Z变换设计数字低通滤波器的特性(基于MATLAB实现)2.3 频带变换及其MATLAB实现除了低通数字滤波器之外，实际中还常常需要高通、带通、及带阻数字滤波器。设计这三种滤波器的方法有很多，例如基于模拟滤波器转换法的IIR数字滤波器设计，基于直接数字域法的IIR数字滤

33、波器设计等等。下面将具体介绍基于滤波器转换法的IIR数字滤波器设计直接由模拟低通滤波器转换成数字高通、带通或带阻滤波器。2.3.2 模拟低通滤波器转换成数字带通滤波器若已知模拟低通滤波器的系统传递函数为G(s)，则模拟低通滤波器s平面到数字带通滤波器的Z平面的变换公式为： (2-18)频率变换公式为： (2-19)其中： (2-20)式中为下截止频率，为上截止频率。从而可得数字带通滤波器的传递函数表达式为： (2-21)接下来，基于MATLAB利用巴特沃斯模拟滤波器举例，设计一数字带通滤波器，要求带通上下截止频率为0.4,0.3,带通内衰减不大于3dB，阻带上下起始频率为0.5，0.2，阻带内

34、衰减不小于18dB。clear all;Wp=0.3*pi,0.4*pi;Ws=0.2*pi,0.5*pi;Ap=3;As=18;N,wn=buttord(Wp/pi,Ws/pi,Ap,As); %计算巴特沃斯滤波器阶次和截至频率b,a=butter(N,wn,bandpass); %频率变换法设计巴特沃斯带通滤波器b0,B,A=dir2cas(b,a) %数字带通滤波器级联型db,mag,pha,grd,w=freqz_m(b,a); %数字滤波器响应subplot(211);plot(w/pi,mag);title(带通数字滤波器幅频响|H(ejOmega)|)subplot(212);p

35、lot(w/pi,db);title(带通数字滤波器幅频响(db)其MATLAB运行结果为：b0 =0.0213B = 1.0000 -2.0000 1.0000 1.0000 2.0000 1.0000A = 1.0000 -0.6249 0.7840 1.0000 -1.0053 0.8061所以由巴特沃斯模拟滤波器设计的数字带通滤波器的传递函数表达式为：图2-11 基于巴特沃斯模拟滤波器设计的数字带通滤波器的幅频响应曲线2.4 小结这一章主要是用MATLAB语言进行IIR滤波器的设计和实现。IIR滤波器的设计步骤分为三步，即模拟低通滤波器设计，模拟-数字滤波器变换，滤波器的频带变换。在模拟低通滤波器的设计中，主要讨论了三种设计方法；在模拟-数字滤波器变换中，讨论了两种变换方法，即脉冲响应不变法和双线性Z变换法；在频带变换的实现中，主要以巴特沃斯滤波器为例并结合MATLAB信号处理工具箱中提供的几个相关函数来进行分析设计。整个设计过程都是在理论分析的基础上，用MATLAB语言来进行编程设计，并最终通过具体滤波器指标来加以实现的。