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

基于MATLAB设计巴特沃斯低通滤波器课程设计专业：XXXXXX姓名：XXX学号： XX指导老师：XXX2011年11 月26日通信系统仿真课程设计任务书院（系）：电气信息工程学院目录1 绪论 (1)1.1 引言 (1)1.2 数字滤波器的设计原理 (1)1.3 数字滤波器的应用 (2)1.4 MATLAB的介绍 (3)1.5 本文的工作及安排 (3)2 滤波器分类及比较 (4)2.1 滤波器的设计原理 (4)2.2 滤波器分类 (4)2.3 两种类型模拟滤波器的比较 (6)3 巴特沃斯低通滤波器 (7)3.1 巴特沃斯低通滤波器简介 (7)3.2 巴特沃斯低通滤波器的设计原理 (7)4 MATLAB仿真及分析 (11)4.1 MATLAB工具箱函数 (11)4.2 巴特沃斯低通滤波器的MATLAB仿真 (11)另附程序调试运行截图： (13)5.1 总结 (13)5.2 展望 (13)1 绪论1.1 引言凡是有能力进行信号处理的装置都可以称为滤波器。

滤波器在如今的电信设备和各类控制系统里面应用范围最广、技术最为复杂，滤波器的好坏直接决定着产品的优劣。

自60年代起由于计算机技术、集成工艺和材料工业的发展，滤波器发展上了一个新台阶，并且朝着低功耗、高精度、小体积、多功能、稳定可靠和价廉方向努力，其中小体积、多功能、高精度、稳定可靠成为70年代以后的主攻方向。

使以数字滤波器为主的各种滤波器得到了飞速的发展，到70年代后期，数字滤波器的单片集成已被研制出来并得到应用。

80年代，致力于各类新型滤波器的研究，努力提高性能并逐渐扩大应用范围。

90年代至现在主要致力于把各类滤波器应用于各类产品的开发和研制。

当然，对数字滤波器本身的研究仍在不断进行。

[1]滤波器主要分成经典滤波器和数字滤波器两类。

从滤波特性上来看，经典滤波器大致分为低通、高通、带通和带阻等。

本文主要对低通数字滤波器做主要研究。

1.2 数字滤波器的设计原理所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。

因此，数字滤波器的概念和模拟滤波相同，只是信号的形式和现实滤波方法不同。

正因为数字滤波器通过数值运算实现实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。

如果要处理的是模拟信号，可以通过A/DC 和D/AC ，在信号形式上进行匹配转换，同样可以使用数字滤波器对模拟信号进行滤波。

[2]大多数的数字滤波器都归类于选频滤波器,其频率响应函数）（ωj e H 如下：）（）（）（ωθωωj j j e e H e H = （1.1）式中，）（ωj e H 称为幅频特性函数；称为相频特性函数。

幅频特性反应的是信号从此滤波器通过后各个频率成分的振幅衰减情况，相频特性表示的是经过滤波器之后各个频率成分在时间上的延时情况。

因此，即使两个滤波器幅频特性相同，而相频特性不同，对相同的输入，滤波器输出的信号波形也是不一样的。

通常情况下幅频特性决定了选频滤波器的技术要求，因为巴特沃斯低通滤波器具有固定的相频特性，所以设计时对相频特性基本没有要求。

图1.1低通滤波器的技术要求图1.1是低通滤波器的幅频特性，p ω和s ω表示通带边界频率和阻带截止频率。

通带频率范围为0≤ω≤p ω，在通带(0,p ω)中要求1)(11≤<-ωδj e H ，阻带频率范围为s ω≤ ω≤π，在阻带(s ω,π)中要求2)(δω≤j e H 。

从p ω到s ω为过渡带，过渡带上的频响一般是单调下降的。

一般情况下用分贝数表示通带及阻带内允许的衰减，通带范围内允许的最大衰减为p α，阻带范围内允许的最小衰减为s α。

p α和s α在低通滤波器里分别用下式定义：dB H H lg20pj 0j p ）（）（ωα= （1.2） dB H H lg 20sj 0j s ）（）（ωα = （1.3） 从上式可以看出p α愈小，通带波纹与通带逼近误差愈小；s α愈大，阻带波纹越小与阻带逼近误差愈小；p ω和s ω之间的距离愈小，过渡带也随之变得更加狭窄。

所以通带边界频率p ω、阻带边界频率s ω、通带最大衰减p α、阻带最小衰减s α决定了低通滤波器的设计指标。

1.3 数字滤波器的应用数字乘法器、加法器及延时单元三者共同构成了数字滤波器。

其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。

由于电子计算机技术和大规模集成电路的发展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。

数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。

[3]数字滤波器在很多领域都得到了广泛的应用，如图像信号处理、语言信号处理、医学生物信号处理等等，其中使用最普遍的是线性时不变数字滤波器。

近年来电子技术发展的非常迅速，数字滤波器的使用范围也随之发生了翻天覆地的改变：从家用的收音机、电视机到航天用的测控设备；从矿井用的通信机到巡航导弹；从超市用的报警器到日常生活的手机，由于电子产品门类及使用频段的不断扩展，各种电子设备之间的干扰也日趋严重，因而数1δ-δ字滤波器不但是确保电子产品本身正常可靠工作的重要部件，而且是减少相互影响、确保正常工作环境的重要器件，因而，可以毫不夸张地说，在具有特定功能的电子产品中均有滤波器的踪迹可寻。

1.4 MATLAB的介绍数字信号处理最重要的部分之一就是数字滤波器的设计， MATLAB的工具箱函数里有非常丰富的相关设计指令，掌握其应用后可以大大提高工作的效率。

MATLAB名字由MATrix和LABoratory两词的前三个字母组合而成。

新墨西哥大学计算机科学系主任Cleve Moler教授在20世纪七十年代年后期，为减轻大学生编程压力，用FORTRAN设计出一组使用简单方便的接口用来调用LINPACK和EISPACK库程序，这就是MATLAB的雏形。

经几年的校际流传，在Little的推动下，由Little、Moler、Steve Bangert 合作，于1984年成立了MathWorks公司，并把 MATLAB正式推向市场。

从此开始采用C语言来编写MATLAB的内核，在原来数值计算能力的基础上还推出了数据图视功能。

MATLAB以商品形式出现后，仅短短几年，就以其良好的开放性和运行的可靠性，使原先控制领域里的封闭式软件包（如英国的UMIST，瑞典的LUND和SIMNON，德国的KEDDC）纷纷淘汰，而改以MATLAB为平台加以重建。

在时间进入20世纪九十年代的时候，MATLAB已经成为国际控制界公认的标准计算软件。

到九十年代初期，在国际上30几个数学类科技应用软件中，MATLAB在数值计算方面独占鳌头。

在欧美大学里，应用代数、数理统计、自动控制、数字信号处理、模拟与数字通信、时间序列分析、动态系统仿真等课程的教科书都把MATLAB作为内容。

这几乎成了九十年代教科书与旧版书籍的区别性标志。

在那里，MATLAB是攻读学位的大学生、硕士生、博士生必须掌握的基本工具。

在国际学术界，MATLAB已经被确认为准确、可靠的科学计算标准软件。

在许多国际一流学术刊物上，（尤其是信息科学刊物），都可以看到MATLAB的应用。

在设计研究单位和工业部门，MATLAB被认作进行高效研究、开发的首选软件工具。

如美国National Instruments公司信号测量、分析软件LabVIEW，Cadence公司信号和通信分析设计软件SPW等，或者直接建筑在MATLAB之上，或者以MATLAB为主要支撑。

又如HP公司的VXI硬件，TM公司的DSP ，Gage公司的各种硬卡、仪器等都接受MATLAB的支持。

[4]1.5本文的工作及安排本文主要工作安排如下：（1）对巴特沃斯低通滤波器的一些基础理论进行详细阐述；（2）对巴特沃斯低通滤波器、切比雪夫Ⅰ型和Ⅱ型滤波器、椭圆滤波器四种典型滤波器加以简单的阐述，然后把四种滤波器加以比较，最后得出巴特沃斯低通滤波器在实际应用中的种种优点；（3）通过双线性变换法设计巴特沃思低通滤波器，然后用MATLAB软件对其进行仿真。

2 滤波器分类及比较2.1 滤波器的设计原理滤波器是使信号中特定的频率成分通过的选频装置，从而达到大幅度衰减其它频率成分的目的。

通过其选频作用，测试装置时就可以分析频谱或滤除干扰噪声。

广义地讲，任何一种信息传输的通道（媒质）都可视为是一种滤波器。

因为，任何装置的响应特性都是激励频率的函数，都可用频域函数描述其传输特性。

因此，构成测试系统的任何一个环节，诸如机械系统、电气网络、仪器仪表甚至连接导线等等，都将在一定频率范围内，按其频域特性，对所通过的信号进行变换与处理。

2.2 滤波器分类⒈从选频作用对滤波器分类 ⑴ 低通滤波器在0～f 2范围内，幅频特性几乎呈水平 发展，低于f 2的频率几乎不衰减通过， 高于f 2的频率几乎不能通过。

⑵ 高通滤波器从f 1～∞，幅频特性几乎呈水平发展。

高于f 1的频率几乎不衰减通过，低于 f 1的频率几乎不能通过。

⑶ 带通滤波器在f 1～f 2之间是带通通滤波器的通频带。

高于f 1低于f 2的频率成分几乎不衰减通过，其它成分几乎不能通过。

⑷ 带阻滤波器频率f 1～f 2之间是傣族滤波器的阻带。

高于f 1低于f 2的频率几乎不能通过，其余频率几乎不衰减通过。

f图2.1 低通滤波器0 f1 f图2.2高通滤波器f1f2图2.3带通滤波器f图2.4 带阻滤波器低通滤波器和高通滤波器是组成滤波器最基本的两种形式，剩下的滤波器都能分解演变成这两种滤波器，如：低通滤波器同高通滤波器进行串联可以成为带通滤波器，低通滤波器同高通滤波器进行并联可以成为带阻滤波器。

⒉ 从“最佳逼近特性”的方面来分类 ⑴ 巴特沃斯滤波器从幅频特性提出要求，而不考虑相频特性。

因为巴特沃斯滤波器具有最大平坦幅度的特性，它的幅频响应如下：n2/11H ）（）（n ωωω+=(2.1)图2.5 巴特沃斯滤波器的频响特性（2）椭圆滤波器的设计在通带和阻带内椭圆滤波器同时具有等波纹幅频响应特性。

因为极点位置与经典场论中的椭圆函数具有一定关联，所以叫做椭圆滤波器。

同时十九世纪三十年代初科学家考尔对椭圆滤波器第一次进行了理论上的证明，它也叫做考尔滤波器。

椭圆滤波器的典型幅频响应特性曲线如图2.7和2.8所示。

由图2.7可见，椭圆滤波器通带和阻带波纹幅度固定时，阶数越高，过渡带越窄；由图2.8可见，当椭圆滤波器阶数固定时，通带和阻带波纹幅度越小，过渡带就越宽。

所以椭圆滤波器的阶数N 由通带边界频率p Ω、阻带边界频率s Ω、通带最大衰减p α和阻带最小衰减s α共同决定。