《基于Matlab软件的数字滤波器设计方法.doc》由会员分享，可在线阅读，更多相关《基于Matlab软件的数字滤波器设计方法.doc（14页珍藏版）》请在三一办公上搜索。
1、 摘要本文分析了国内外数字滤波技术的应用现状与发展趋势，介绍了数字滤波器的基本结构，在分别讨论了IIR与FIR数字滤波器的设计方法的基础上，指出了传统的数字滤波器设计方法过程复杂、计算工作量大、滤波特性调整困难的不足，提出了一种基于Matlab软件的数字滤波器设计方法，完成了高Q值50Hz带通IIR滤波器的设计, 达到了通带45-55Hz，衰减小于3db，阻带40-60Hz，衰减大于80db的参数指标。文中深入分析了该滤波器系统设计的功能特点、实现原理以及技术关键，阐述了使用MATLAB进行带通滤波器设计及仿真的具体方法。最后把整个设计方案用VHDL语言进行了描述并在Modelsim上仿真。M
2、odelsim与Matlab的仿真结果对比说明该设计准确性好，可精确到小数点后六位，稳定后误差小于万分之一；可移植性强，在实际应用中，可根据不同的阶数、精度和速度等要求对IIR 滤波器系数进行灵活的修改，以实现任意阶数的IIR 滤波器。因此，该设计方法可靠性好，效率高，极大的减轻了工作量，有利于滤波器设计的最优化。关键词 带通滤波器；IIR；Matlab；VHDL；Modelsim仿真;AbstractThis paper analyzes the situation of application and development of digital filter technology hom
3、e and abroad. It introduces the basic structure of a digital filter, discusses different design methods of FIR and IIR filter, and points out that the traditional design method of digital filter is not only complex but also of heavy workload, even adjustment of filtering parametrer is very difficult
4、. So it brings forwad another design method of digital filter which based on the Matlab software and Modelsim software. This paper introduces the design method of a high Q value band-pass IIR filter which meets the given standard - the pass band is 45-55Hz, with attenuation of less than 3db; the sto
5、p band is 40-60Hz with attenuation greater than 80db. It deeply analyzes the design features and principles of the filter system and the key technical in the design. Then it describs the use of Matlab in design and simulation of the bandpass filte design. In the end, the procedure was simulated on M
6、odelsim. Simulation results of Modelsim and Matlab compared to make proof of the accuracy of the design. It is proved that the design can be accurate at the sixth decimal place, stable error is within one ten thousandth; The design has good portability and widely practical application. In different
7、applications, according to different performance requirements of order, speed and accuracy, etc, the IIR filter coefficients can change flexibly to make up of IIR filter of arbitrary order; In addition, VHDL procedure can also make certain optimizations, according to the different actual situation,
8、to achieve greater speed or conservation of FPGA resources.Key wordsdigital filter IIR Matlab VHDL Modelsim simulation 第一章 绪论1.1 数字滤波器的研究背景与意义当今，数字信号处理(DSP：Digtal Signal Processing)技术正飞速发展，它不但自成一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设紧密相连；它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能
9、化和网络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号、雷达信号、通信信号、导航信号、射电天文信号、生物医学信号、控制信号、气象信号、地震勘探信号、机械振动信号、遥感遥测信号等。上述这些信号大部分是模拟信号，也有小部分是数字信号。模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。大多数情况下一维模拟信号的自变量是时间，经过时间上的离散化(采样)和幅度上的离散化(量化)，这类模拟信号便成为一维数字信号。因此，数字信号实际上是用数字序列表示的信号，语音信号经采样和量化后，得到的数字信号是一个一维离散时间序列；而图像信号经采样和量化后，得到的数字信号是一个二维离散
10、空间序列。数字信号处理，就是用数值计算的方法对数字序列进行各种处理，把信号变换成符合需要的某种形式。例如，对数字信号经行滤波以限制他的频带或滤除噪音和干扰，或将他们与其他信号进行分离；对信号进行频谱分析或功率谱分析以了解信号的频谱组成，进而对信号进行识别；对信号进行某种变换，使之更适合于传输，存储和应用；对信号进行编码以达到数据压缩的目的，等等。数字滤波技术是数字信号分析、处理技术的重要分支。无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传输是至关重要的。在所有的电子系统中，使用最多技术最复杂的要算数字滤波器了。数字滤波器的优劣直接决定产品的优劣。1
11、.2 数字滤波器的应用现状与发展趋势在信号处理过程中，所处理的信号往往混有噪音，从接收到的信号中消除或减弱噪音是信号传输和处理中十分重要的问题。根据有用信号和噪音的不同特性，提取有用信号的过程称为滤波，实现滤波功能的系统称为滤波器。在近代电信设备和各类控制系统中，数字滤波器应用极为广泛，这里只列举部分应用最成功的领域。（1）语音处理语音处理是最早应用数字滤波器的领域之一，也是最早推动数字信号处理理论发展的领域之一。该领域主要包括5个方面的内容：第一，语音信号分析。即对语音信号的波形特征、统计特性、模型参数等进行分析计算；第二，语音合成。即利用专用数字硬件或在通用计算机上运行软件来产生语音；第三
12、，语音识别。即用专用硬件或计算机识别人讲的话，或者识别说话的人；第四，语音增强。即从噪音或干扰中提取被掩盖的语音信号。第五，语音编码。主要用于语音数据压缩，目前已经建立了一系列语音编码的国际标准，大量用于通信和音频处理。近年来，这5个方面都取得了不少研究成果，并且，在市场上已出现了一些相关的软件和硬件产品，例如，盲人阅读机、哑人语音合成器、口授打印机、语音应答机，各种会说话的仪器和玩具，以及通信和视听产品大量使用的音频压缩编码技术。（2）图像处理数字滤波技术以成功地应用于静止图像和活动图像的恢复和增强、数据压缩、去噪音和干扰、图像识别以及层析X射线摄影，还成功地应用于雷达、声纳、超声波和红外信
13、号的可见图像成像。（3）通信在现代通信技术领域内，几乎没有一个分支不受到数字滤波技术的影响。信源编码、信道编码、调制、多路复用、数据压缩以及自适应信道均衡等，都广泛地采用数字滤波器，特别是在数字通信、网络通信、图像通信、多媒体通信等应用中，离开了数字滤波器,几乎是寸步难行。其中，被认为是通信技术未来发展方向的软件无线电技术，更是以数字滤波技术为基础。（4）电视数字电视取代模拟电视已是必然趋势。高清晰度电视的普及指日可待，与之配套的视频光盘技术已形成具有巨大市场的产业；可视电话和会议电视产品不断更新换代。视频压缩和音频压缩技术所取得的成就和标准化工作，促成了电视领域产业的蓬勃发展，而数字滤波器及
14、其相关技术是视频压缩和音频压缩技术的重要基础。（5）雷达雷达信号占有的频带非常宽，数据传输速率也非常高，因而压缩数据量和降低数据传输速率是雷达信号数字处理面临的首要问题。告诉数字器件的出现促进了雷达信号处理技术的进步。在现代雷达系统中，数字信号处理部分是不可缺少的，因为从信号的产生、滤波、加工到目标参数的估计和目标成像显示都离不开数字滤波技术。雷达信号的数字滤波器是当今十分活跃的研究领域之一。（6）声纳声纳信号处理分为两大类，即有源声纳信号处理和无源声纳信号处理，有源声纳系统涉及的许多理论和技术与雷达系统相同。例如，他们都要产生和发射脉冲式探测信号，他们的信号处理任务都主要是对微弱的目标回波进
15、行检测和分析，从而达到对目标进行探测、定位、跟踪、导航、成像显示等目的，他们要应用到的主要信号处理技术包括滤波、门限比较、谱估计等。（7）生物医学信号处理数字滤波器在医学中的应用日益广泛，如对脑电图和心电图的分析、层析X射线摄影的计算机辅助分析、胎儿心音的自适应检测等。（8）音乐数字滤波器为音乐领域开辟了一个新局面，在对音乐信号进行编辑、合成、以及在音乐中加入交混回响、合声等特殊效果特殊方面，数字滤波技术都显示出了强大的威力。数字滤波器还可用于作曲、录音和播放，或对旧录音带的音质进行恢复等。（9）其他领域数字滤波器的应用领域如此广泛，以至于想完全列举他们是根本不可能的，除了以上几个领域外，还有
16、很多其他的应用领域。例如，在军事上被大量应用于导航、制导、电子对抗、战场侦察；在电力系统中被应用于能源分布规划和自动检测；在环境保护中被应用于对空气污染和噪声干扰的自动监测，在经济领域中被应用于股票市场预测和经济效益分析，等等。1.3 本章小结 数字滤波器精确度高、使用灵活、可靠性高，具有模拟设备所没有的许多优点，已广泛地应用于各个科学技术领域, 例如数字电视、语音、通信、雷达、声纳、遥感、图像、生物医学以及许多工程应用领域。随着信息时代数字时代的到来，数字滤波技术已经成为一门极其重要的学科和技术领域。以往的滤波器大多采用模拟电路技术，但是，模拟电路技术存在很多难以解决的问题，例如，模拟电路元
17、件对温度的敏感性，等等。而采用数字技术则避免很多类似的难题，当然数字滤波器在其他方面也有很多突出的优点，在前面部分已经提到，这些都是模拟技术所不能及的，所以采用数字滤波器对信号进行处理是目前的发展方向。 第二章 设计方案2.1 IIR滤波器的基本结构一个数字滤波器可以用系统函数表示为： (2.1-1) 由这样的系统函数可以得到表示系统输入与输出关系的常系数线形差分程为： (2.1-2) 可见数字滤波器的功能就是把输入序列x(n)通过一定的运算变换成输出序列y(n)。不同的运算处理方法决定了滤波器实现结构的不同。无限冲激响应滤波器的单位抽样响应h(n)是无限长的，其差分方程如(2-2)式所示，是
18、递归式的，即结构上存在着输出信号到输入信号的反馈，其系统函数具有(2-1)式的形式，因此在z平面的有限区间(0z)有极点存在。前面已经说明，对于一个给定的线形时不变系统的系统函数，有着各种不同的等效差分方程或网络结构。由于乘法是一种耗时运算，而每个延迟单元都要有一个存储寄存器，因此采用最少常熟乘法器和最少延迟支路的网络结构是通常的选择，以便提高运算速度和减少存储器。然而，当需要考虑有限寄存器长度的影响时，往往也采用并非最少乘法器和延迟单元的结构。IIR滤波器实现的基本结构有：(1)IIR滤波器的直接型结构；优点：延迟线减少一半，变为N 个，可节省寄存器或存储单元；缺点：其它缺点同直接I型。 通
19、常在实际中很少采用上述两种结构实现高阶系统，而是把高阶变成一系列不同组合的低阶系统(一、二阶)来实现。(2)IIR滤波器的级联型结构；特点： 系统实现简单，只需一个二阶节系统通过改变输入系数即可完成； 极点位置可单独调整； 运算速度快(可并行进行)； 各二阶网络的误差互不影响，总的误差小，对字长要求低。 缺点：不能直接调整零点，因多个二阶节的零点并不是整个系统函数的零点，当需要准确的传输零点时，级联型最合适。(3)IIR滤波器的并联型结构。优点： 简化实现，用一个二阶节，通过变换系数就可实现整个系统； 极、零点可单独控制、调整，调整1i、2i只单独调整了第i对零点，调整1i、2i则单独调整了第
20、i对极点； 各二阶节零、极点的搭配可互换位置，优化组合以减小运算误差； 可流水线操作。 缺点： 二阶阶电平难控制，电平大易导致溢出，电平小则使信噪比减小。 a、直接型 b、并联型c、串联型图2.1-3、IIR滤波器的基本结构2.2数字滤波器的设计原理数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。这种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计
21、方法有多种，如双线性变换法、窗函数设计法、插值逼近法和Chebyshev逼近法等等。随着MATLAB软件尤其是MATLAB的信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。数字滤波器设计的基本步骤如下：（1）确定指标在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给出幅度和相位响应。幅度指标主要以两种方式给出。第一种是绝对指标。它提供对幅度响应函数的要求，一般应用于FIR滤波器的设计。第二种指标是相对指标。它以分贝值的形式给出要求。在工程实际中，这
22、种指标最受欢迎。对于相位响应指标形式，通常希望系统在通频带中具有线性相位。运用线性相位响应指标进行滤波器设计具有如下优点：只包含实数算法，不涉及复数运算；不存在延迟失真，只有固定数量的延迟；长度为N的滤波器(阶数为N-1)，计算量为N/2数量级。因此，本文中滤波器的设计就以线性相位FIR滤波器的设计为例。（2）逼近确定了技术指标后，就可以建立一个目标的数字滤波器模型。通常采用理想的数字滤波器模型。之后，利用数字滤波器的设计方法，设计出一个实际滤波器模型来逼近给定的目标。（3）性能分析和计算机仿真上两步的结果是得到以差分或系统函数或冲激响应描述的滤波器。根据这个描述就可以分析其频率特性和相位特性
23、，以验证设计结果是否满足指标要求；或者利用计算机仿真实现设计的滤波器，再分析滤波结果来判断。2.2.1 滤波器的性能指标我们在进行滤波器设计时，需要确定其性能指标。一般来说，滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。以低通滤波器特性为例，频率响应有通带、过渡带及阻带三个范围。在通带内： 1- AP| H(ej)| 1 |c在阻带中： |H(ej)| Ast st |c其中c 为通带截止频率, st为阻带截止频率，Ap为通带误差, Ast为阻带误差。图2.2-1、低通滤波器频率响应幅度特性的容限图与模拟滤波器类似，数字滤波器按频率特性划分为低通、高通、带通、带阻、全通等类型，由于
24、数字滤波器的频率响应是周期性的，周期为2。各种理想数字滤波器的幅度频率响应如图所示：图2.2-2、 各种理想数字滤波器的幅度频率响应2.3. IIR数字滤波器的设计方法 目前，IIR数字滤波器设计最通用的方法是借助于模拟滤波器的设计方法。模拟滤波器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图表供查询，因此，充分利用这些已有的资源将会给数字滤波器的设计带来很大方便，IIR数字滤波器的设计步骤是：（1）按一定规则将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标；（2）根据转换后的技术指标设计模拟低通滤波器H(s)；（3）在按一定规则将H(s)转换为H(z)。若
25、所设计的数字滤波器是低通的，那么上述设计工作可以结束，若所设计的是高通、带通或者带阻滤波器，那么还有步骤：（4）将高通、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标，然后按上述步骤(2)设计出模拟低通滤波器H(s)，再由冲击响应不变法或双线性变换将H(s)转换为所需的H(z)。s - z 映射的方法有：冲激（脉冲）响应不变法、阶跃响应不变法、双线性变换法等。下面讨论冲激（脉冲）响应不变法和双线性变换法。一般来说，在要求时域冲激响应能模仿模拟滤波器的场合，一般使用冲激响应不变法。冲激响应不变法一个重要特点是频率坐标的变换是线性的，因此如果模拟滤波器的频响带限于折叠频率的话，则通过
26、变换后滤波器的频响应可不失真的反映原响应与频率的关系。所谓冲激（脉冲）响应不变法就是使数字滤波器的脉冲序列h（n）等于模拟滤波器的脉冲响应ha(t)的采样值，即h(t)=ha（t）|t=Nt=ha(nT)式中，T为采样周期。 因此数字滤波器的系统函数H(z)可由下式求得： H(z)=Zh(n)=Zha(nT)Z-表示对-的内容进行Z变换，Z变换的内容请参考相应的数字信号处理教材。如果已经获得了满足性能指标的模拟滤波器的传递函数Ha(s),求与之对应的数字滤波器的传递函数H（z）的方法是：（1）求模拟滤波器的单位脉冲响应ha(t)。 ha(t) =L-1Ha（s）式中，L-1Ha（s）表示对Ha
27、（s）的Laplace变换内容请参考高等数学的积分变换或信号处理教材。（2）求模拟滤波器单位冲激响应h(n)的采样值，即数字滤波器冲激响应序列h(n)。（3）对数字滤波器的冲激响应h(n)进行z变换，得到传递函数H(z)。 由上述方法推论出更直接地由模拟滤波器系统函数Ha(s)求出数字滤波器系统函数H(z)的步骤是： （1）利用部分分式展开将模拟滤波器的传递函数H(s)展开成： Ha(s)=双线性变换法是指首先把s 平面压缩变换到某一中介平面s1 的一条横带(宽度为2T,即从- T到T) ,然后再利用的关系把s1平面上的这条横带变换到整个z 平面。这样s 平面与z 平面是一一对应关系, 消除了
28、多值变换性, 也就消除了频谱混叠现象。s 平面到z 平面的变换可采用 (2.3-1) (2.3-2)令 ，有： (2.3-3)从s1 平面到z 平面的变换,即 (2.3-4)代入上式，得到： (2.3-5) 一般来说，为使模拟滤波器的某一频率与数字滤波器的任一频率有对应关系，可引入代定常数c， (2.3-6)则 (2.3-7) 这种s 平面与z 平面间的单值映射关系就是双线性变换。有了双线性变换，模拟滤波器的数字化只须用进行置换。 第三章 设计实现3.1 IIR反陷波器的实现及Simulink仿真分析 FDATool滤波器设计分析工具箱，他涵盖了 信号处理工具箱中所有的滤波器设计方法。利用这个
29、工具箱可以使复杂的滤波器最佳设计方法变得简单，用户只需在FDATool图形工具箱选择相应的滤波器设计方法，并填入滤波器设计指标，即可为用户生成所需的滤波器频率响应特性。根据技术指标要求，采用间接法设计一个带通的一阶IIR数字滤波器。模拟滤波器原型为切比雪夫1型，阻带上限截止频率为80Hz,下限截止频率为120Hz，阻带波纹为1Db,其幅频与相频特性分别如图：Simulink是MATLAB中一个专门用于对动态系统进行建模、仿真和分析的软件包，它支持连续离散及两者混合的线性和非线性系统，也支持具有多种采样速率的多速率系统。MATLAB/SIMULINK中嵌入了Embeddedtarget for
30、TI C2000 DSP,它是针对C2000系列的开发工具，从而可以在实际系统制作出来之前，预先对系统进行仿真与分析，并可以对系统做适当的实时修正或者按照仿真的最佳效果来调试及整定系统的参数，以提高系统的性能，减少系统设计过程中反复修改的时间，实现高效率地开发系统的目标。通过调试Simulink中的功能模块，可以构成仿真数字滤波器特性的结构框。对上述设计的IIR带通滤波器进行Simulink仿真分析，其仿真模型如图 所示。输入的两个信号为x1(n)=sin(30n)、x2(n)=5sin(100n)和x3(n)=sin(300n),经过加法器进行叠加，其中x2(n)为有用信号，x1(n)和x3(n)为干扰信号。经过所设计的IIR带阻数字滤波器滤波。滤波前后的时域波形分别经各自的示波器模块显示。仿真结果分别如图 所示滤波前的仿真图 滤波后的仿真结果图2.5 本章小结 第四章 总结4.1 脉冲响应与双线性变换法的比较4.2 设计心得和体会
