用matlab制作带通滤波器

1、学科分类号（二级）510.8020YU N N AN NORMAL U N IVER S ITY本科学生毕业论文（设计）题 目用MATLAB设计带通滤波器姓名杨洪斌学号074090092院、系物理与电子信息学院专业物理学指导教师王六玲职称（学历）副教授用MATLA破计带通滤波器摘要： 关于滤波器设计的传统思路，本文借助Matlab 软件中滤波器设计工具，探讨一种更为简洁、高效的设计方式。传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。本文介绍了一种利用 MATLA茜号处理工具箱(Signal Processing Toolbox )快速有效的设计由软件组成的常规

2、数字滤波器的设计方法。本文给出 了使用MATLAB©言进行程序设计和利用信号处理工具箱的FDATool工具进行界面设计的方法。利用 MATLA股计滤波器，可以随时对比设计要求和滤波器特性 调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。关键词：带通滤波器；MATLAB设计工具随着信息科学与计算技术的迅速发展，数字信号处理的理论与应用得到飞跃式发展，形成了一门极为重要的学科1。滤波技术是信号分析、处理技术的重要分支。无论是信号的获取、传输, 还是信号的处理和交换都离不开滤波技术, 它对信号安全可靠和有效灵活地传递是至关重要的。故滤波技术是极为敏感与热门的课题 , 对滤

3、波器的研制也历来为各国所重视3 。1 MATLA歆件简介和滤波器概述1.1 MATLA瞅件简介MATLAB件具有很强的开放性和适应性。 在保持内核不变的情况下，MATLAB可以针对不同的应用学科推出相应的工具箱(toolbox) ，目前己经推出了图象处理工具箱、信号处理工具箱、小波工具箱、神经网络工具箱以及通信工具箱等多个学科的专用工具箱，极大地方便了不同学科的研究工作。国内已有越来越多的科研和技术人员认识到MATLAB强大作用，并在不同的领域内使用MATLAB快速实现科研构想和提高工作效率。MATLAB供了 20类图像处理函数，涵盖了图像处理的包括近期研究成果在内的几乎所有的技术方法, 是学

4、习和研究图像处理的人员难得的宝贵资料和加工工具箱。这些函数按其功能可分为: 图像显示; 图像文件I/O; 图像算术运算; 几何变换 ; 图像登记; 像素值与统计; 图像分析; 图像增强 ; 线性滤波 ; 线性二元滤波设计;图像去模糊; 图像变换; 邻域与块处理; 灰度与二值图像的形态学运算; 结构元素创建与处理; 基于边缘的处理; 色彩映射表操作; 色彩空间变换; 图像类型与类型转换。1.2 滤波器滤波器， 顾名思义，其作用是对输入信号起到滤波的作用，数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。

5、数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应(IIR )滤波器和有限长冲激响应(FIR)滤波器。IIR 数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配。所以 iir 滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。FIR数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。我们平时所要设计的数字滤波器，阶数和类型并不一定是完全给定的，很多时候都是要根据设计要求和滤波效果不断的调整，以达到设计的最优化。在这种情况下， 滤

6、波器的设计就要进行大量复杂的运算，单纯的靠公式计算和编制简单的程序很难在短时间内完成设计。利用MATLAB大的计算功能进行计算机辅助设计，可以快速有效的设计数字滤波器，大大的简化了计算量，直观简便。1.3 带通滤波器带通滤波器是指能通过某一频率范围内的频率分量、但将其他范围的频率分量衰减到极低水平的滤波器，与带阻滤波器的概念相对。一个模拟带通滤波器的例子是电阻- 电感 - 电容电路 （RLC circuit） 。这些滤波器也可以用低通滤波器同高通滤波器组合来产生。除了电子学和信号处理领域之外，带通滤波器应用的一个例子是在大气科学领域， 很常见的例子是使用带通滤波器过滤最近3到 10天时间范围内

7、的天气数据，这样在数据域中就只保留了作为扰动的气旋。还有，许多音响装置的频谱分析器均使用此电路作为带通滤波器，以选出各个不同频段的信号，在显示上利用发光二极管点亮的多少来指示出信号幅度的大小。带通滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。在高保真的信号处理，如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。2 滤波器设计工具2.1 基本功能滤波器设计工具以MATLA的计算环境和信号处理工具箱（SPTool）为基础。主要功能包括：仿真分析滤波器、FFT 和 IFFT 计算、滤波器设计计算（包括自适应滤波器） 、支持滤波器的频带变换。操作

8、该工具可以设计出满足各种性能指标的滤波器，并可方便的查看它们的各性能分析图形。设计完成后，可以把其系数结果直接导出为 MATLABE量保存记录。2.2 基本操作在安装了 MATLAB勺PC机上运行该软件，在其命令窗口输入FDA Tool便可打开该设计工具, 如图 1 所示 :13Filter Design & Analysie Tool - urtitled.fdaFile Edit Analysis Targets View W'indow Help AR集岛 理分x D KJ曰口旧熟中田 .,Li 膻尸。7if3g riHaxStruciure： Direct-Form F

9、IROrder:50Stable:YesSourceDesigned orCHlEer .l ResponsnwrwasKvHiohoass Bandstop'Differenti. -Design,Butter. 9 FEquirip. ReadyFilter Speci. 11 d I q Minim.-OptiDensity 2 0 FrequencyUhitApas 1Aston ®0图i滤波器设计分析工具Fig.1 filter design analysis toolsFDATool界面总共分两大部分，一部分是design filter ,在界面的下半部， 用来设

10、置滤波器的设计参数，另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。design filter 部分主要分为：filter type（滤波器类型）选项，包括 lowpass （低通）、highpass （高通）、bandpass （带通）、bandstop （带阻）和特殊的fir滤波器。design method （设计方法）选项，包括 iir 滤波器的butterworth （巴特 沃思）法、chebyshev type i （切比雪夫 i 型）法、chebyshev type ii （切比 雪夫ii型）法、elliptic（椭圆滤波器）法和fir滤波器的equiripple

11、 法、least-squares （最小乘方）法、window （窗函数）法。filter order （滤波器阶数）选项，定义滤波器的阶数，包括specify order （指定阶数）和 minimum order （最小阶数）。在specify order中填入所要设计的滤波器的阶数（n阶滤波器，specify order =n-1 ）,如果选择 minimumorder 则matlab根据所选择的滤波器类型自动使用最小阶数。frenquency specifications 选项，可以详细定义频带的各参数，包括采样 频率fs和频带的截止频率。它的具体选项由 filter type选项和d

12、esign method 选项决定，例如bandpass （带通）滤波器需要定义fstop1 （下阻带截止频率）、 fpass1 （通带下限截止频率）、fpass2 （通带上限截止频率）、fstop2 （上阻带截止频率），而lowpass （低通）滤波器只需要定义fstopl、fpassl。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。magnitude specifications 选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义 wstopl （频率fstopl处的幅值衰减）、wpass （通带范 围内的幅值衰

13、减）、wstop2 （频率fstop2处的幅值衰减）。当采用窗函数设计 时，通带截止频率处的幅值衰减固定为6db,所以不必定义。window specifications 选项，当选取采用窗函数设计时，该选项可定义，它包含了各种窗函数。3 带通滤波器设计3.1 设计要求在小电流接地系统中注入83.3hz 的正弦信号，对其进行跟踪分析，要求设计一带通滤波器，滤除工频及整次谐波，以便在非常复杂的信号中分离出该注入信号。参数要求：96阶fir数字滤波器，采样频率1000hz,采用hamming窗函数设计。3.2 设计工具本例中，首先在filter type 中选择bandpass （带通滤波器）；（

14、1）在design method选项中选择fir window （fir 滤波器窗函数法），接着在 window specifications 选项中选取hamming；（2） 指定 filter order 项中的 specify order =95;（ 3）由于采用窗函数法设计，只要给出通带下限截止频率fc1 和通带上限截止频率 fc2 ,选取 fc1 =70hz, fc2 =84hz。（ 4）设置完以后点击design filter 即可得到所设计的fir 滤波器。（ 5）通过菜单选项analysis 可以在特性区看到所设计滤波器的幅频响应、相频响应、零极点配置和滤波器系数等各种特性。设

15、计完成后保存为1.fda 文件。在设计过程中，可以对比滤波器幅频相频特性和设计要求，随时调整参数和滤波器类型，以便得到最佳效果。其它类型的fir 滤波器和iir 滤波器也都可以使用 fdatool 来设计。3.3 设计编程MATLABJ部有丰富的函数可供调用，使用者只需要根据自己的需要查到所需要的函数名, 那就可以直接调用使用了. 不需要自己把函数详细的编写出来. 使编写内部程序变的简单。MATLABJ号工具箱主要提供了以下几种窗函数，如表 1所示：程序法只以hanning 窗为示例 ， 其他窗函数的程序大体相同， 只是在 window=hanning（N） 和 hn = fir1（N-1,w

16、c, hanning （N） 处要调自己的窗函数即可。表1 MATLAB窗函数Tab.1 MATLAB window function矩形窗Boxcar巴特利特窗Barlett三角窗Triang布莱克曼窗Blackman海明窗Hamming汉宁窗Hanning凯塞窗Kaiser切比雪夫窗Chebwinhanning窗设计程序：wls = 0.2*pi;wlp = 0.35*pi;whp = 0.65*pi;wc = wlp/pi,whp/pi;B = wlp-wls;N = ceil(8/0.15);n=0:N-1;window= hanning(N);h1,w=freqz(window,1)

17、;figure(1);stem(window);axis(0 60 0 1.2);grid;xlabel('n');title('Hanning 窗函数')；figure(2);plot(w/pi,20*log(abs(h1)/abs(h1(1);axis(0 1 -350 0);grid;xlabel('w/pi');ylabel('幅度(dB)');title('Hanning 窗函数的频谱');hn = fir1(N-1,wc, hanning (N);h2,w=freqz(hn,1,512);巾gure(3

18、);stem(n,hn);axis(0 60 -0.25 0.25);grid；xlabel('n');ylabel('h(n)');title( ' Hanning窗函数的单位脉冲响应');figure(4);plot(w/pi,20*log(abs(h2)/abs(h2(1);grid；xlabel('w/pi');ylabel('幅度(dB)');title( ' Hanning窗函数的单位幅频响应');4仿真结果与分析4.1 仿真结果图2汉宁窗函数波形图Fig.2 hanning windo

19、w function waveform figure汉宁窗（Hanning）函数时域形式表示为：W（n）4+MRN（公式1）汉宁窗是典型的开余弦窗，在运行程序后得到Hanning窗函数波形类似余弦 函数波形，基本符合汉宁窗函数的时域表达式。Fig.3 hanning window function spectrum diagram利用傅利叶变换的调制特性，由公式 1可得汉宁窗的频谱函数为：W（ej8）=10.5WRg ）+0.25Wr «WrNr> +-2沱一）二产L:< N -3 I N -1川.N 1=W® e '2 ,（公式 2）NJ上式中，NR（

20、e）=WR侬e 产（公式3）当N远大于1时，上式可近似表示为：（公式4）W侬 卜 0.5WR 侬)+0.25.Nr,舒卜这三部分之和使旁瓣互相抵消，能量更集中在主瓣。在运行程序后得到的Hanning窗函数的频谱图中阻带衰减大约在-70db到-300db。主瓣宽度也比旁瓣宽度大，能量更集中在主瓣。图4汉宁窗设计带通滤波器的冲击响应图Fig.4 hanning window design bandpass filter graph of shock response在matlab环境下运行汉宁窗设计程序，可得到 Hanning窗函数的单位脉冲 响应图，知道带通滤波器单位脉冲响应为偶对称。25020

21、0160100S-画 6。0图5汉宁窗设计带通滤波器的幅频响应图Fig.5 hanning window design band-pass filter the amplitude-frequency response figure运行程序后得到的Hanning窗函数的单位幅频响应图，清楚看到在频率0.35冗-0.65冗的信号会平稳输出，而在通带之外的频率的信号会大大衰减。4.2 结果分析由仿真波形图可以得出以下几个方面的分析结果：由图 2 可以看出，汉宁窗是典型的升余弦窗，基本符合汉宁窗函数的时域表达式，所加的窗口是余弦函数（类似余弦函数）。由图 3 可以看出，汉宁窗函数的阻带衰减大约在-7

22、0dB 到 -300dB, 也就是说旁瓣幅度远小于主瓣幅度，能量全部集中在主瓣，主瓣宽度也有所增加，效果比较明显。由图 4 可以看出，汉宁窗设计的带通滤波器的单位冲击响应为偶对称，也就是说当N为偶数时，单位冲击响应序列为偶函数，这说明FIR带通滤波器相位是 严格的线性关系。由图 5 可以看出，汉宁窗设计的带通滤波器带通下限截止频率大约为0.35冗，带通上线截止频率大约为0.65冗，在通带频率之外的频率的信号的幅度将受到很大的衰减，以致信号不能通过滤波器。根据以上的结果分析可以得到：利用matlab 设计的带通滤波器，经仿真成像后，与理论实际相符，本次设计结果符合设计要求。利用 matlab 的

23、强大运算功能，基于 matlab 信号处理工具箱的数字滤波器设计法可以快速有效的设计滤波器，设计方便、快捷，极大的减轻了工作量。在设计过程中可以对比滤波器特性，随时更改参数，以达到滤波器设计的最优化, 可以反映实际情况。参考文献1 丁玉美， 高西全 . 数字信号处理M. 第三版 . 西安： 西安电子科技大学出版社,2008， 6.2 刘令普 . 数字信号处理M. 哈尔滨：哈尔滨工业大学出版社，2002， 4.3 邓重一 . 滤波技术的发展现状J. 中国仪器仪表.2004， 2( 2) ： 1-4.4 王世一 . 数字信号处理M. 北京：北京理工大学出版社，2005， 6.5 郑南宁，程洪. 数

24、字信号处理M. 北京：清华大学出版社，2007， 9.6 刘正士， 王勇， 陈恩伟等. 一种数字滤波器的设计方法及其应用J. 中国机械工程 .2006， 1， 17(1) ： 88-91.7 凌燮亭 , 秦巍 , 胡波 . 对数域积分电路的分析与应用J. 电子学报，2000， 28(2) ： 46-48.8 蔡理， 马西奎 . 一种基于平衡式对数域积分器的高阶滤波器J. 固体电子学研究与进展，2002， 22(1) ： 68-71.9 赵录怀 . 高阶对数域滤波器设计的一种新方法J. 电工技术学报，1999,14(5) ： 27-30.10 White S A. Applications of

25、 Dist ributed Arithmetic to Digital Sgnal ProcessingJ.IEEE ASSP Magazine,1989,6(3):4210.11 R. Mahesh, A. P. Vinod, Edmund M-K. Lai,etal.Filter Bank Channelizers for Multi-Standard Software Defined Radio ReceiversJ.Journal of Signal Processing Systems,2010,10,61(3):251-387.Design of bandpass filter using MATLABAbstract: Aiming at the t raditional way of the filter design in the exisiting situation of the teaching material about the process of dig