巴特沃斯低通滤波器分析报告Word文档下载推荐.docx

1、 impulse response method, the bilinear transformation method.The design process are by the analog filter system function h （s） to transform the digital filter of the system function h （z）.Keywords: Digital filter MATLAB Butterworth Chebyshev Bilinear transformation methed 1.前言所谓数字滤波器，是指输入、输出均为数字信号，通

2、过数值运算处理改变输入信号所含频率成分的相比照例，或者滤除某些频率成分的数字器件或程序。因此，数字滤波器的概念和模拟滤波一样，只是信号形式和实现滤波方式不同。正因为数字滤波通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波公能。数字滤波器是一个离散时间系统。应用数字滤波器处理模拟信号时，首先须对输入模拟信号进展限带、抽样和模数转换。数字滤波器输入信号的抽样率应大于被处理信号带宽的两倍。数字滤波器的频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即二分之一抽样频率点呈镜像对称。为得到模拟信号，数字滤波器处理

3、的输出数字信号须经数模转换、平滑。数字滤波器有低通、高通、带通、带阻和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。如果数字滤波器的部参数不随时间而变化，那么称为时不变的，否那么为时变的。如果数字滤波器在某一给定时刻的响应与在此时刻以后的鼓励无关，那么称为因果的，否那么为非因果的。如果数字滤波器对单一或多个鼓励信号的响应满足线性条件，那么称为线性的，否那么为非线性的。应用最广的是线性、时不变数字滤波器。与模拟滤波器一样，数字滤波器从功能上看可分为低通、高通、带通和带阻几类。因为它们的单位脉冲响应是非因果且无限长的，所以实际上理想滤波器是不可能实现的。与模拟滤波器不同

4、的是数字滤波器的传输函数都是以2为周期的，滤波器的低通频带处于2的整数倍处，而高频频带处于的奇数倍附近。根据对不同信号的处理可将滤波器分为模拟滤波器和数字滤波器两种。模拟滤波器和数字滤波器的概念一样，只是信号的形式和实现滤波方法不同。数字滤波器是指输入输出都是数字信号的滤波器。数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。这种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。2、数字滤波器的概述2.1数字滤波器的分类按照不同的分类方法，数字滤波器有许多种类，但总起来

5、可以分成两大类:经典滤波器和现代滤波器。经典滤波器的特点是其输入信号中有用的频率成分和希望滤除的频率成分各占有不同频带，通过一个适宜的选频滤波器滤除干扰，得到纯洁信号，到达滤波的目的。但是如果信号和干扰的频谱相互重叠，那么经典滤波器不能有效地滤除干扰，最大限度的恢复信号，这时就需要现代滤波器。滤波器的分类有很多的方式，如按元件、信号处理的方式、通频带分类、运放电路的构成等进展分类。 2.1滤波器的分类 2.2数字滤波器的技术指标技术指标是设计滤波器重要的一局部。在一般情况下，滤波器需要明确它的幅度特性，才能更好的进展设计。比方低通滤波器的频率响应特性，它有过渡带、通带和阻带三围，如图。在通带：

6、1-AP|H（ej）|1|c在阻带上：|H（ej）|Ast Wst|Wc参数如下：Ap通带误差；c通带截止频率；Ast阻带误差；st阻带截止频率低通滤波器频率响应幅度特性的容限图模拟滤波器和数字滤波器根据它的频率特性，有低通，高通，带通，带阻等类型，数字滤波器的响应频率的周期是2。如下图不同频率的滤波器：数字滤波器的幅度频率响应2.3数字滤波器的原理数字信号分析在设计数字滤波器中是必不可少的。数字滤波器被广泛运用于，如录音机、DVD、音响等。随着科学技术的开展，模拟滤波器有很多的缺乏之处，所以逐渐被数字滤波器取代。而与模拟滤波相比拟，数字滤波器具有准确、系统稳定和灵活性高等优点。但是信号在的滤

7、波，抽样和转换方面，模拟模拟滤波器仍然充当着重要的角色。 数字滤波器是利用信号数字化、离散化的方法来加工信号的。并且用数字方法根据信号的要求进展转换。简而言之，数字滤波器就是一个离散系统，通过运算来处理输入信号，滤除高频干扰信号，从而得到有用的输出信号。模拟信号的数字处理如图2-3所示。模拟信号通过抗混叠滤波器的进展频带处理，然后通过AD转换和数字处理模块，再通过DA转换为模拟信号，最后，再用滤波器滤除这个模拟信号中的干扰信号。 2.3模拟信号的数字处理3.巴特沃斯低通滤波器的设计 3.1巴特沃思低通滤波器的设计原理巴特沃斯滤波器以巴特沃斯函数来近似滤波器的系统函数。巴特沃斯滤波器是根据幅频特

8、性在通频带具有最平坦特性定义的滤波器。巴特沃思滤波器的低通模平方函数表示： N=1,2，.下面归纳了巴特沃斯滤波器的主要特征a.对所有的N，。 b.对所有的N， c.是的单调下降函数。d.随着阶次N的增大而更接近于理想低通滤波器。如下列图3.1所示，可以看出滤波器的幅频特性随着滤波器阶次N的增加而变得越来越好，在截止频率c处的函数值始终为1/2的情况下，通带有更多的频带区的值接近于1；在阻带更迅速的趋近于零。 3.1巴特沃思低通滤波平方幅频特性函数3.2系统函数设巴特沃斯的系统函数为Has，那么： 3.3设计过程巴特沃思低通滤波技术指标关系式为 ap-20log|Ha（j）|，P ass其中：

9、p为通带边界频率，s为阻带边界频率。经过化简整理可得：取满足上式的最小整数N作为滤波器的阶数。再将N代入可得： 或查表求得归一化传输函数H（s）,令s/c代替归一化原型滤波器系统函数中的s,即得到实际滤波器传输函数。（1）根据技术指标p、s和s，求出滤波器的阶数。（2）求出归一化极点k,将k代入,得到归一化低通原型系统函数Ga（p）。（3）将Ga（p）去归一化。将p=s/c代入Gap,得到实际的滤波器系统函数c为3dB截止频率4.Matlab设计与实现设计一个通带截止频率fp=5kHz,通带最大衰减p=2dB,阻带截止频率fs=12kHz,阻带最小衰减s=30dB的巴特沃斯低通滤波器。 4.1

10、Matlab程序设计1Z,P,K=butter（N）2N,wc=buttord（wp,ws,Rp,As）3N,wc=buttord（wp,ws,Rp,As,s）4B,A=butter（N,wc,ftype）5B,A=butter（N,wc,ftype,s）设计的程序如下： wp=2*pi*5000; ws=2*pi*12000; Rp=2; As=30; N,wc=buttord（wp,ws,Rp,As,s）N = 5wc = 3.7792e+004 B,A=butter（N,wc,B = 1.0e+022 * 0 0 0 0 0 7.7094A = 0.0000 0.0000 0.0000

11、0.0000 0.0007 7.7094 k=0:511; fk=0:14000/512:14000; wk=2*pi*fk; Hk=freqs（B,A,wk）; subplot（2,2,1）; plot（fk/1000,20*log10（abs（Hk）; grid on xlabel（频率（kHz）; ylabel（幅度dB axis（0,14,-40,5） 4.2 Matlab的实现实现结果如下列图4.1、4.2、4.3 4.1程序输入 4.2运算结果4.3巴特沃斯运行输出的损耗函数 4.3 Matlab结果分析滤波器的损耗函数曲线如图4.3所示。由图可以看出，阻带刚好满足指标要求，通带指

12、标有充裕。5.总结使用Matlab设计的过程中，找到了需要设计的低通滤波器的函数与依据。依据设计的滤波器的wp,ws,Rp,Rs四个参数来确定输入信号的频率以实现过滤的目的。实验过程中也遇到很多问题，通过查找资料都得到了解决。通过这次学习，掌握了matlab的仿真和设计，对低通滤波器设计原理的加深了解与技术指标的设置。这次设计不仅让我对数字信号处理这门学科的理论知识加深理解，而且让我对MATLAB中的许多函数的运用更加灵活，还增强了我的解决问题的能力。参考文献1 高西全、丁玉美等编 数字信号处理 电子科技大学2008年8月第3版2 志涌、祖樱等编 MATLAB教程 航空航天大学 2015年1月第1版