基于MATLAB做巴特沃斯低通滤波器

1、基于MATLAB设计巴特沃斯低通滤波器课程设计专 业： XXXXXX姓 名： XXX学 号： XX指导老师： XXX2011年 11 月 26日15通信系统仿真课程设计任务书院（系）：电气信息工程学院目 录1 绪论 .11.1 引言 . 11.2 数字滤波器的设计原理 . 11.3 数字滤波器的应用 . 21.4 MATLAB的介绍 . 31.5 本文的工作及安排 . 32 滤波器分类及比较 .42.1 滤波器的设计原理 . 42.2 滤波器分类 . 42.3 两种类型模拟滤波器的比较 . 63 巴特沃斯低通滤波器 .73.1 巴特沃斯低通滤波器简介 . 73.2 巴特沃斯低通滤波器的设计原理

2、 . 74 MATLAB仿真及分析. 114.1 MATLAB工具箱函数 . 114.2 巴特沃斯低通滤波器的MATLAB仿真 . 11另附程序调试运行截图： . 135.1 总结 . 135.2 展望 . 1315通信系统仿真课程设计1 绪论1.1 引言凡是有能力进行信号处理的装置都可以称为滤波器。滤波器在如今的电信设备和各类控制系统里面应用范围最广、技术最为复杂，滤波器的好坏直接决定着产品的优劣。自60年代起由于计算机技术、集成工艺和材料工业的发展，滤波器发展上了一个新台阶，并且朝着低功耗、高精度、小体积、多功能、稳定可靠和价廉方向努力，其中小体积、多功能、高精度、稳定可靠成为70年代以后

3、的主攻方向。使以数字滤波器为主的各种滤波器得到了飞速的发展，到70年代后期，数字滤波器的单片集成已被研制出来并得到应用。80年代，致力于各类新型滤波器的研究，努力提高性能并逐渐扩大应用范围。90年代至现在主要致力于把各类滤波器应用于各类产品的开发和研制。当然，对数字滤波器本身的研究仍在不断进行。1滤波器主要分成经典滤波器和数字滤波器两类。从滤波特性上来看，经典滤波器大致分为低通、高通、带通和带阻等。本文主要对低通数字滤波器做主要研究。1.2 数字滤波器的设计原理所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。

4、因此，数字滤波器的概念和模拟滤波相同，只是信号的形式和现实滤波方法不同。正因为数字滤波器通过数值运算实现实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。如果要处理的是模拟信号，可以通过A/DC和D/AC，在信号形式上进行匹配转换，同样可以使用数字滤波器对模拟信号进行滤波。2大多数的数字滤波器都归类于选频滤波器,其频率响应函数H（ej）如下：H（e）=H（eejjj（） （1.1） 式中，H（ej称为幅频特性函数；称为相频特性函数。幅频特性反应的是信号从此滤波器通过后各个频率成分的振幅衰减情况，相频特性表示的是经过滤波

5、器之后各个频率成分在时间上的延时情况。因此，即使两个滤波器幅频特性相同，而相频特性不同，对相同的输入，滤波器输出的信号波形也是不一样的。通常情况下幅频特性决定了选频滤波器的技术要求，因为巴特沃斯低通滤波器具有固定的相频特性，所以设计时对相频特性基本没有要求。15通信系统仿真课程设计1-图1.1低通滤波器的技术要求图1.1是低通滤波器的幅频特性，p和s表示通带边界频率和阻带截止频率。通带频率范围为0p，在通带(0,p)中要求1-1<H(ej)1，阻带频率范围为s，在阻带(s,)中要求H(ej)2。从p到s为过渡带，过渡带上的频响一般是单调下降的。一般情况下用分贝数表示通带及阻带内允许的衰减

6、，通带范围内允许的最大衰减为p，阻带范围内允许的最小衰减为s。p和s在低通滤波器里分别用下式定义：H（ H（jpj0p=20lgdB （1.2） s=20lgH（ H（jsj0dB （1.3） 从上式可以看出p愈小，通带波纹与通带逼近误差愈小；s愈大，阻带波纹越小与阻带逼近误差愈小；p和s之间的距离愈小，过渡带也随之变得更加狭窄。所以通带边界频率p、阻带边界频率s、通带最大衰减p、阻带最小衰减s决定了低通滤波器的设计指标。1.3 数字滤波器的应用数字乘法器、加法器及延时单元三者共同构成了数字滤波器。其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。由于电子计算机技术和大规

7、模集成电路的发展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。3数字滤波器在很多领域都得到了广泛的应用，如图像信号处理、语言信号处理、医学生物信号处理等等，其中使用最普遍的是线性时不变数字滤波器。 近年来电子技术发展的非常迅速，数字滤波器的使用范围也随之发生了翻天覆地的改变：从家用的收音机、电视机到航天用的测控设备；从矿井用的通信机到巡航导弹；从超市用的报警器到日常生活的手机，由于电子产品门类及使用频段的不断扩展，各种电子设备之间的干扰也日趋严重，因而数2通信系统仿真课程设计字滤波器不但是确保电子产品本

8、身正常可靠工作的重要部件，而且是减少相互影响、确保正常工作环境的重要器件，因而，可以毫不夸张地说，在具有特定功能的电子产品中均有滤波器的踪迹可寻。1.4 MATLAB的介绍数字信号处理最重要的部分之一就是数字滤波器的设计， MATLAB的工具箱函数里有非常丰富的相关设计指令，掌握其应用后可以大大提高工作的效率。MATLAB名字由MATrix和LABoratory两词的前三个字母组合而成。新墨西哥大学计算机科学系主任Cleve Moler教授在20世纪七十年代年后期，为减轻大学生编程压力，用FORTRAN设计出一组使用简单方便的接口用来调用LINPACK和EISPACK库程序，这就是 MATLA

9、B的雏形。经几年的校际流传，在Little的推动下，由Little、Moler、Steve Bangert 合作，于1984年成立了MathWorks公司，并把 MATLAB正式推向市场。从此开始采用C语言来编写MATLAB的内核，在原来数值计算能力的基础上还推出了数据图视功能。MATLAB以商品形式出现后，仅短短几年，就以其良好的开放性和运行的可靠性，使原先控制领域里的封闭式软件包（如英国的UMIST，瑞典的LUND和SIMNON，德国的KEDDC）纷纷淘汰，而改以 MATLAB为平台加以重建。在时间进入20世纪九十年代的时候，MATLAB已经成为国际控制界公认的标准计算软件。到九十年代初期

10、，在国际上30几个数学类科技应用软件中， MATLAB在数值计算方面独占鳌头。在欧美大学里，应用代数、数理统计、自动控制、数字信号处理、模拟与数字通信、时间序列分析、动态系统仿真等课程的教科书都把MATLAB作为内容。这几乎成了九十年代教科书与旧版书籍的区别性标志。在那里，MATLAB是攻读学位的大学生、硕士生、博士生必须掌握的基本工具。在国际学术界，MATLAB已经被确认为准确、可靠的科学计算标准软件。在许多国际一流学术刊物上，（尤其是信息科学刊物），都可以看到MATLAB的应用。在设计研究单位和工业部门，MATLAB被认作进行高效研究、开发的首选软件工具。如美国National Instr

11、uments公司信号测量、分析软件LabVIEW，Cadence公司信号和通信分析设计软件SPW等，或者直接建筑在MATLAB之上，或者以MATLAB为主要支撑。又如HP公司的VXI硬件，TM公司的DSP ，Gage公司的各种硬卡、仪器等都接受MATLAB的支持。41.5 本文的工作及安排本文主要工作安排如下：（1） 对巴特沃斯低通滤波器的一些基础理论进行详细阐述；（2） 对巴特沃斯低通滤波器、切比雪夫型和型滤波器、椭圆滤波器四种典型滤波器加以简单的阐述，然后把四种滤波器加以比较，最后得出巴特沃斯低通滤波器在实际应用中的种种优点；（3） 通过双线性变换法设计巴特沃思低通滤波器，然后用MATLA

12、B软件对其进行仿真。3通信系统仿真课程设计2 滤波器分类及比较2.1 滤波器的设计原理滤波器是使信号中特定的频率成分通过的选频装置，从而达到大幅度衰减其它频率成分的目的。通过其选频作用，测试装置时就可以分析频谱或滤除干扰噪声。广义地讲，任何一种信息传输的通道（媒质）都可视为是一种滤波器。因为，任何装置的响应特性都是激励频率的函数，都可用频域函数描述其传输特性。因此，构成测试系统的任何一个环节，诸如机械系统、电气网络、仪器仪表甚至连接导线等等，都将在一定频率范围内，按其频域特性，对所通过的信号进行变换与处理。2.2 滤波器分类从选频作用对滤波器分类 低通滤波器在0f2范围内，幅频特性几乎呈水平发

13、展，低于f2的频率几乎不衰减通过，高于f2的频率几乎不能通过。 高通滤波器从f1，幅频特性几乎呈水平发展。高于f1的频率几乎不衰减通过，低于f1的频率几乎不能通过。0 图2.1 低通滤波器 f 带通滤波器在f1f2之间是带通通滤波器的通频带。高于f1低于f2的频率成分几乎不衰减通过，其它成分几乎不能通过。 带阻滤波器频率f1f2之间是傣族滤波器的阻带。高于f1低于f2的频率几乎不能通过，其余频率几乎不衰减通过。0 f1 f2 f图2.2高通滤波器f1 f 图2.3带通滤波器图2.4 带阻滤波器4通信系统仿真课程设计低通滤波器和高通滤波器是组成滤波器最基本的两种形式，剩下的滤波器都能分解演变成这

14、两种滤波器，如：低通滤波器同高通滤波器进行串联可以成为带通滤波器，低通滤波器同高通滤波器进行并联可以成为带阻滤波器。 从“最佳逼近特性”的方面来分类 巴特沃斯滤波器从幅频特性提出要求，而不考虑相频特性。因为巴特沃斯滤波器具有最大平坦幅度的特性，它的幅频响应如下：H（=1+（/n）2n(2.1)H(-100-300/c图2.5巴特沃斯滤波器的频响特性（2）椭圆滤波器的设计在通带和阻带内椭圆滤波器同时具有等波纹幅频响应特性。因为极点位置与经典场论中的椭圆函数具有一定关联，所以叫做椭圆滤波器。同时十九世纪三十年代初科学家考尔对椭圆滤波器第一次进行了理论上的证明，它也叫做考尔滤波器。椭圆滤波器的典型幅

15、频响应特性曲线如图2.7和2.8所示。由图2.7可见，椭圆滤波器通带和阻带波纹幅度固定时，阶数越高，过渡带越窄；由图2.8可见，当椭圆滤波器阶数固定时，通带和阻带波纹幅度越小，过渡带就越宽。所以椭圆滤波器的阶数N由通带边界频率p、阻带边界频率s、通带最大衰减p和阻带最小衰减s共同决定。5它的典型幅频响应特性如图2.7和2.8所示。5通信系统仿真课程设计图2.7 椭圆滤波器的典型幅频响应特性曲线图2.8 椭圆滤波器的典型幅频响应特性曲线2.3 两种类型模拟滤波器的比较图2.9和图2.10是巴特沃思和椭圆滤波器的频响特性曲线。调用MATLAB滤波器涉设计函数，很容易验证：当阶数相同时，对相同的通带

16、最大衰减和阻带最小衰减，巴特沃斯滤波器具有单调下降的幅频特性，过渡带最宽。比椭圆滤波器的过渡带宽。巴特沃斯滤波器的最大平坦幅度特性致使它在实际中应用最为广泛，所以本文主要对巴特沃斯低通滤波器进行研究。6通信系统仿真课程设计3 巴特沃斯低通滤波器3.1 巴特沃斯低通滤波器简介巴特沃斯滤波器是电子滤波器的一种，特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。这种滤波器最先由英国工程师斯替芬·巴特沃斯（Stephen Butterworth）在1930年发表在英国无线电工程期刊的一篇论文中提出的，可以构成低通、高通、带通和带阻四种组态，是目前最为流行的一类数字

17、滤波器 ,经过离散化可以作为数字巴特沃思滤波器 ,较模拟滤波器具有精度高、稳定、灵活、不要求阻抗匹配等众多优点 ,因而在自动控制、语音、图像、通信、雷达等众多领域得到了广泛的应用，是一种具有最大平坦幅度响应的低通滤波器。 873.2 巴特沃斯低通滤波器的设计原理 巴特沃斯低通滤波器的幅度平方函数H（j）用下式表示： aH（j）a22=1+（1c）2N （3.1）j=1；=c时，H（j=1/2,c是3dBN为滤波器的阶数。当=0时，H（aa截止频率。=c时，逐渐增大，幅度下降非常迅速。、N同幅度特性关系如图3.1所示。N决定了幅度下降速度，N越大，通带就越平坦，过渡带也随之变窄，阻带幅度同过渡带

18、下降的速度越迅速，总体频响特性同理想低通滤波器的实际误差越小。10.707 c2图3.1 、N同幅度特性关系11+（sjc）2N j）变成s的函数： 用s代替j，把幅度平方函数H（aH（s）H（-s）=aa （3.2）s=+j，此公式说明了幅度平方函数有2N个极点，极点sk可以用下面的公式来表达：7通信系统仿真课程设计j（2k+1）sk=（jc2N）=c 2jj2k+12N=c12k+1j（+22N（3.3）k=0,1,2，···，2N-1。2N个极点等间隔分布在半径为c的圆上，间隔是/N rad。如图3.2所示：图3.2 三阶巴特沃斯滤波器极点分布s），而右为形

19、成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成H（as）的表示式为 -s）半平面的N个极点构成H（。H（aaH（s）=acN-1N（3.4）k(s-sk=0)例如N=3,通过下式可以计算出6个极点s0=c2j3，s1=cj，s2=c4j3，s3=c4j35j3，s4=cj2，s5=c7j3当N=3时，6个极点中位于左半平面的三个分别为:s0=c2j3，s1=c j=-c，s2=cc32-j3cs）： 取s平面左半平面的极点s0，s1，s2组成H（aH（s）=a1N-1（s+c）（s-2j3c（3.5）（s-将H（s）=a（k=0sc-skc对3dB截止频率c归一化：）s）归一化后

20、的H（表示为： aH（s）=a1N-1（k=0sc-skc（3.6）8通信系统仿真课程设计令p=sc=jc=j，p称为归一化拉氏复变量。 =/c，称为归一化频率。经过归一化后巴特沃斯滤波器的传输函数为：H（p）=a1N-1 （3.7）k（p-pK=0）式中，pk=sk/c为归一化极点，sk为位于左半平面的极点用下式表示： pk= 12k+1j（+22N ，k=0,1， N-1 （3.8）将极点表示式(3.8)代入(3.7)式，得到的H（的分母是p的N阶多项式，用下式表p）a示：1Ha(p)=（3.9） 2N-1N b0+b1p+b2p+bN-1p+p下面来确定N由技术指标p，p s和s确定。在

21、定义p=-10lgH（jp （3.10）a22（3.11）p=-10lgH（jp aH（j=a1+（21c）2N （3.12）js，再将 中，将=p和=s分别代入(6.2.6)式中，得到H（jp和H（aa22H（jpa2js代入(3.13)和(3.14)式中，得到： 和H（a2p=-10lg1+(1p, )2N （3.13）cs=-10lg1+(1sc)2N （3.14）整理得：p2N1+（）=10p/10 （3.15） c1+（sc）2N=10s/10 （3.16）由(3.15)和(3.16)式得到：9通信系统仿真课程设计（ps）=N1010p/10s/10-1-1 （3.17） 令sp=s

22、/p，ksp=则N由下式表示： N=-lgksp1010p/10s/10-1-1 （3.18）lgsp （3.19）取大于等于N的最小整数。 关于3dB截止频率c，如果技术指标中没有给出，可以按照(3.18)式或(3.19)式求出， 1+（1+（pcsc）2N=10=10p/10 （3.20）2Ns/10 （3.21）由(3.20)式得到： c=（10p0.1p-1）2N （3.22） -1由(3.21)式得到： c=（10s0.1s-1）2N (3.23)lgk-1经过总结，巴特沃斯低通滤波器的设计步骤大致为： 通过p,p,s和s的值，用公式N=-根据公式pk=H（p）=a1N-1splgs

23、p算出滤波器的阶数N0。 j（12+2k+12N，k=0,1， N-1，求出归一化极点pk，将pk代入 p）中，得出归一化传输函数H（。 ak（p-pk=0）p）p）将H（去归一化。将p=s/c代入H（之中，从而得到实际的滤波器传输函数aaH（s）。 a10通信系统仿真课程设计4 MATLAB仿真及分析4.1 MATLAB工具箱函数MATLAB的信号处理工具箱提供了滤波器的函数 buttap、buttord、butter。由z,p,k = buttap(n)函数可设计出n阶巴特沃斯低通滤波器原型,其传递函数为H(s)=1(s-p(1)(s-p(2)(s-p(n)ks+bn-isnn-1 （4.

24、1） 所以事实上z为空阵。上述零极点形式可以化为： H(s)=+b1s+b0 （4.2）其中b0=cn，令c=1rad/s，得到巴特沃斯滤波器归一化结果,如表1所示。表1 n=18阶的巴特沃斯滤波器系数 n b7 b6 b5 b4 b3 b2 b1 b0buttord函数可在给定滤波器性能的情况下 ,选巴特沃斯滤波器的阶数 n 和截止频率c,从而可用butter函数设计巴特沃斯滤波器的传递函数。9n ,c = buttord (p,c,Rp,Rs,s) 可得到足性能的模拟巴特沃斯滤波器的最小阶数 n及截止频率c ,其中p为通带的拐角频率, s为阻带的拐角频率, p和s的单位均为rad/s; R

25、s为通带区的最大波动系数,Rp为Rs阻带区的最小衰减系数, Rs和Rp的单位都为dB。 b,a = butter (n,c,s) 可设计截止频率为c的n 阶低通模拟巴特沃斯滤波器,其传递函数为:H(s)= 10 b(1)s+b(2)ss+a(2)snnn-1B(s)A(s)=+b(n+1)n-1+a(n+1) (4.3)4.2 巴特沃斯低通滤波器的MATLAB仿真由于冲击响应不变法和阶跃响应不变法的不足 ,下面以双线性变换法设计一个巴特沃思低通滤波器。通过仔细研究MATLAB软件自带的信号处理工具箱下%:MATLABD IR%15通信系统仿真课程设计toolbox signal signal buttord.m和%MATLABD IR% toolbox signalsignalbutter.m ,可以清楚看到:MATLAB默认采用双线性变换法设计,其实完全可以采用更简洁的方法。8滤波器的性能指标如下:通带截止频率fp=2400Hz，通带最大衰减p=3dB，阻带截止频率fs=5000Hz，阻带最小衰减s=25dBMATLAB程序如下:%巴特沃斯低通滤波器f_N=16000; %采样率f_p=2400; f_s=5000; R_p