学科分类号（二级）510.8020
题 目 用MATLAB设计带通滤波器
姓 名 杨 洪 斌
学 号 *********
院、 系 物理与电子信息学院
专 业 物理学
指导教师 王六玲
职称（学历） 副教授
用 MATLAB设计带通滤波器
摘要：关于滤波器设计的传统思路，本文借助Matlab 软件中滤波器设计工具，探讨一种更为简洁、高效的设计方式。
传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。本文介绍了一种利用MATLAB信号处理工具箱（Signal Processing Toolbox）快速有效的设计由软件组成的常规数字滤波器的设计方法。本文给出了使用MATLAB语言进行程序设计和利用信号处理工具箱的FDA Tool工具进行界面设计的方法。利用MATLAB设计滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。 关键词：带通滤波器；MATLAB；设计工具
随着信息科学与计算技术的迅速发展，数字信号处理的理论与应用得到飞跃式发展，形成了一门极为重要的学科[1]。滤波技术是信号分析、处理技术的重要分支。无论是信号的获取、传输, 还是信号的处理和交换都离不开滤波技术,它对信号安全可靠和有效灵活地传递是至关重要的。故滤波技术是极为敏感与热门的课题, 对滤波器的研制也历来为各国所重视[3]。
1 MATLAB软件简介和滤波器概述
1.1 MATLAB软件简介
MATLAB软件具有很强的开放性和适应性。在保持内核不变的情况下，MATLAB可以针对不同的应用学科推出相应的工具箱(toolbox)，目前己经推出了图象处理工具箱、信号处理工具箱、小波工具箱、神经网络工具箱以及通信工具箱等多个学科的专用工具箱，极大地方便了不同学科的研究工作。国内已有越来越多的科研和技术人员认识到MATLAB的强大作用，并在不同的领域内使用MATLAB来快速实现科研构想和提高工作效率。
MATLAB提供了20类图像处理函数,涵盖了图像处理的包括近期研究成果在内的几乎所有的 技术方法,是学习和研究图像处理的人员难得的宝贵资料和加工工具箱。这些函数按其功能可分为:图像显示;图像文件I/O;图像算术运算;几何变换;图像登记;像素值与统计;图像分析;图像增强;线性滤波;线性二元滤波设计;图像去模糊;图像变换;邻域与块处理;灰度与二值图像的形态学运算;结构元素创建与处理;基于边缘的处理;彩映射表操作;彩空间变换;图像类型与类型转换。
1.2 滤波器
滤波器，顾名思义，其作用是对输入信号起到滤波的作用，数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器。
IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配。所以iir滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。
segg
FIR数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。
我们平时所要设计的数字滤波器，阶数和类型并不一定是完全给定的，很多时候都是要根据设计要求和滤波效果不断的调整，以达到设计的最优化。在这种情况下，滤波器的设计就要进行大量复杂的运算，单纯的靠公式计算和编制简单的程序很难在短时间内完成设计。利用MATLAB强大的计算功能进行计算机辅助设计，可以快速有效的设计数字滤波器，大大的简化了计算量，直观简便。
1.3 带通滤波器
带通滤波器是指能通过某一频率范围内的频率分量、但将其他范围的频率分量衰减到极低水平的滤波器，与带阻滤波器的概念相对。一个模拟带通滤波器的例子是电阻-电感-电容电路(RLC circuit)。这些滤波器也可以用低通滤波器同高通滤波器组合来产生。
除了电子学和信号处理领域之外，带通滤波器应用的一个例子是在大气科学领域，很常见
的例子是使用带通滤波器过滤最近3到10天时间范围内的天气数据，这样在数据域中就只保留了作为扰动的气旋。还有，许多音响装置的频谱分析器均使用此电路作为带通滤波器，以选出各个不同频段的信号，在显示上利用发光二极管点亮的多少来指示出信号幅度的大小。
带通滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。在高保真的信号处理，如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。
2 滤波器设计工具
四氧化锰
2.1 基本功能
滤波器设计工具以MATLAB的计算环境和信号处理工具箱(SPTool) 为基础。主要功能包括：仿真分析滤波器、FFT 和IFFT 计算、滤波器设计计算(包括自适应滤波器)、支持滤波器的频带变换。操作该工具可以设计出满足各种性能指标的滤波器，并可方便的查看它们的各性能分析图形。设计完成后，可以把其系数结果直接导出为MATLAB 变量保存记录。
2.2 基本操作快递系统
在安装了MATLAB的PC 机上运行该软件,在其命令窗口输入FDA Tool 便可打开该设计工具,如图1 所示:
裤带蛇图1 滤波器设计分析工具
Fig.1连通区域 filter design analysis tools
FDA Tool界面总共分两大部分，一部分是design filter，在界面的下半部，用来设置滤波器的设计参数，另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。design filter部分主要分为：
filter type（滤波器类型）选项，包括lowpass（低通）、highpass（高通）、bandpass（带通）、bandstop（带阻）和特殊的fir滤波器。
design method（设计方法）选项，包括iir滤波器的butterworth（巴特沃思）法、chebyshev type i（切比雪夫i型）法、 chebyshev type ii（切比雪夫ii型） 法、elliptic（椭圆滤波器）法和fir滤波器的equiripple法、least-squares（最小乘方）法、window（窗函数）法。
filter order（滤波器阶数）选项，定义滤波器的阶数，包括specify order（指定阶数）和minimum order（最小阶数）。在specify order中填入所要设计的滤波器的阶数（n阶滤波
器，specify order＝n-1），如果选择minimum order则matlab根据所选择的滤波器类型自动使用最小阶数。
frenquency specifications选项，可以详细定义频带的各参数，包括采样频率fs和频带的截止频率。它的具体选项由filter type选项和design method选项决定，例如bandpass（带通）滤波器需要定义fstop1（下阻带截止频率）、fpass1（通带下限截止频率）、fpass2 （通带上限截止频率）、fstop2（上阻带截止频率），而lowpass（低通）滤波器只需要定义fstop1、fpass1。采用窗函数设计滤波器 时，由于过渡带是由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。
magnitude specifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义wstop1（频率fstop1处的幅值衰减）、wpass （通带范围内的幅值衰减）、wstop2（频率fstop2处的幅值衰减）。当采用窗函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。
window specifications选项，当选取采用窗函数设计时，该选项可定义，它包含了各种窗函数。
3 带通滤波器设计
3.1设计要求
在小电流接地系统中注入83.3hz的正弦信号，对其进行跟踪分析，要求设计一带通滤波器，滤除工频及整次谐波，以便在非常复杂的信号中分离出该注入信号。参数要求：96阶fir数字滤波器，采样频率1000hz，采用hamming窗函数设计。
3.2设计工具
本例中，首先在filter type中选择bandpass（带通滤波器）；
（1）在design method选项中选择fir window（fir滤波器窗函数法），接着在window specifications选项中选取hamming；
（2）紫砂饮水机指定filter order项中的specify order＝95；
（3）由于采用窗函数法设计，只要给出通带下限截止频率fc1和通带上限截止频率fc2，选取fc1＝70hz，fc2＝84hz。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议