本科毕业设计（论文）-数字滤波器设计.pdf

本科毕业设计（论文）数字滤波器设计目录摘要.1第1章绪论.21.1 数字滤波器的研究背景与意义.21.2 数字滤波器的应用现状与发展趋势.21.3 数字滤波器的实现方法分析.41.4 本章小结.4第2章 数字滤波器的概述.52.1 数字滤波器的基本结构.52.2 HR滤波器的基本结构.52.3 FIR滤波器的基本结构.72.4 数字滤波器的设计原理.82.5 滤波器的性能指标.92.6 HR数字滤波器的设计方法.92.7 FIR数字滤波器的设计方法.102.8 HR滤波器与FIR滤波器的分析比较.122.9 本章小节.13第3章 数字滤波器的算法设计.143.1 由模拟滤波器设计HR数字滤波器.143.2 巴特奥兹滤波器.143.3 切比雪夫滤波器.153.4 椭圆滤波器.173.5 用MATLAB设计数字滤波器.183.6 FDATool 界面.183.7 用FDATOOL进行带通滤波器设计.203.8 将系统函数由直接型化成级联型.223.9 二阶节系数的确定.233.10 系数转换成二进制码.233.11 木章小结.25第4章仿真过程.274.1 设定系统的仿真对象.274.2 系统对象滤波器设计方法.274.3 MATLAB程序仿真设计.284.3.1产生一个含有10Hz.30Hz和60Hz的混合正弦波信号X.284.3.2产生一个4阶IIR带通椭圆滤波器.294.3.3对混合正弦波信号X进行滤波.304.3.4绘出信号滤波前、后的幅频图.314.4 创建仿真模型图.344.4.1仿真模块参数设置.354.5 系统仿真运行.384.6 小结.39第5章总结.405.1 滤波器功能和性能总结.405.2 设计心得和体会.40第6章结束语.41参考文献.42附录.43译文.47外文原文51数字滤波器设计摘要本文分析了国内外数字滤波技术的应用现状与发展趋势，介绍了数字滤波器的基本结 构，在分别讨论了 HR与FIR数字滤波器的设计方法的基础上，指出了传统的数字滤波器 设计方法过程复杂、计算工作量大、滤波特性调整困难的不足，提出了一种基于Ma tla b 和Modelsim软件的数字滤波器设计方法，完成了高Q值50Hz带通HR滤波器的设计，达 到了通带45-55Hz,衰减小于3db,阻带40-60Hz,衰减大于80db的参数指标。文中深入 分析了该滤波器系统设计的功能特点、实现原理以及技术关键，阐述了使用MATLAB进行 带通滤波器设计及仿真的具体方法。最后把整个设计方案用VHDL语言进行了描述并在 Modelsim上仿真。Modelsim与Ma tla b的仿真结果对比说明该设计准确性好，可精确到小 数点后六位，稳定后误差小于万分之一；可移植性强，在实际应用中，可根据不同的阶 数、精度和速度等要求对IIR滤波器系数进行灵活的修改，以实现任意阶数的IIR滤波 器。因此，该设计方法可靠性好，效率高，极大的减轻了工作量，有利于滤波器设计的最 优化。关键词 带通滤波器；IIR；Ma t la b；VHDL；Modelsim仿真；AbstractThis pa per a n a lyzes the situa tion of a pplica tion a n d developmen t of digita l filter techn ology home a n d a broa d.Tt in troduces the ba sic structure of a digita l filter,discusses differen t design methods of FIR a n d IIR filter,a n d poin ts out tha t the tra dition a l design method of digita l filter is n ot on ly complex but a lso of hea vy workloa d,even a djustmen t of filterin g pa ra metrer is very difficult.So it brin gs forwa d a n other design method of digita l filter which ba sed on the Ma tla b softwa re a n d Modelsim softwa re.This pa per in troduces the design method of a high Q va lue ba n d-pa ss IIR filter which meets the given sta n da rd the pa ss ba n d is 45-55Hz,with a tten ua tion of less tha n 3db;the stop ba n d is 40-60Hz with a tten ua tion grea ter tha n 80db.It deeply a n a lyzes the design fea tures a n d prin ciples of the filter system a n d the key techn ica l in the design.Then it describs the use of Ma tla b in design a n d simula tion of the ba n dpa ss filte design.In the en d,the procedure wa s simula ted on Modelsim.Simula tion results of Modelsim a n d Ma tla b compa red to ma ke proof of the a ccura cy of the design.It is proved tha t the design ca n be a ccura te a t the sixth decima l pla ce,sta ble error is within on e ten thousa n dth;The design ha s good porta bility a n d widely pra ctica l a pplica tion.In differen t a pplica tion s,a ccordin g to differen t performa n ce requiremen ts of order,speed a n d a ccura cy,etc,the IIR filter coefficien ts ca n cha n ge flexibly to ma ke up of IIR filter of a rbitra ry order;In a ddition,VHDL procedure ca n a lso ma ke certa in optimiza tion s,a ccordin g to the differen t a ctua l situa tion,to a chieve grea ter speed or con serva tion of FPGA resources.Key words digita l filter IIR Ma tla b VHDL Modelsim simula tion1数字滤波器设计第1章绪论1.1 数字滤波器的研究背景与意义当今，数字信号处理(DSP：Digta l Sign a l Processin g)技术正飞速发展，它不但自成 一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设 紧密相连；它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化 的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号、雷达信号、通信信 号、导航信号、射电天文信号、生物医学信号、控制信号、气象信号、地震勘探信号、机 械振动信号、遥感遥测信号，等等。上述这些信号大部分是模拟信号，也有小部分是数字 信号。模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。大 多数情况下一维模拟信号的自变量是时间，经过时间上的离散化(采样)和幅度上的离散化(量化)，这类模拟信号便成为一维数字信号。因此，数字信号实际上是用数字序列表示的 信号，语音信号经采样和量化后，得到的数字信号是一个一维离散时间序列；而图像信号 经采样和量化后，得到的数字信号是一个二维离散空间序列。数字信号处理，就是用数值 计算的方法对数字序列进行各种处理，把信号变换成符合需要的某种形式。例如，对数字 信号经行滤波以限制他的频带或滤除噪音和干扰，或将他们与其他信号进行分离；对信号 进行频谱分析或功率谱分析以了解信号的频谱组成，进而对信号进行识别；对信号进行某 种变换，使之更适合于传输，存储和应用；对信号进行编码以达到数据压缩的目的，等 等。数字滤波技术是数字信号分析、处理技术的重要分支已引。无论是信号的获取、传 输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传输是至 关重要的。在所有的电子系统中，使用最多技术最复杂的要算数字滤波器了。数字滤波器 的优劣直接决定产品的优劣。1.2 数字滤波器的应用现状与发展趋势在信号处理过程中，所处理的信号往往混有噪音，从接收到的信号中消除或减弱噪音 是信号传输和处理中十分重要的问题。根据有用信号和噪音的不同特性，提取有用信号的 过程称为滤波，实现滤波功能的系统称为滤波器。在近代电信设备和各类控制系统中，数 字滤波器应用极为广泛，这里只列举部分应用最成功的领域。(1)语音处理语音处理是最早应用数字滤波器的领域之一，也是最早推动数字信号处理理论发展的 领域之一。该领域主要包括5个方面的内容：第一，语音信号分析。即对语音信号的波形 特征、统计特性、模型参数等进行分析计算；第二，语音合成。即利用专用数字硬件或在 通用计算机上运行软件来产生语音；第三，语音识别。即用专用硬件或计算机识别人讲的 话，或者识别说话的人；第四，语音增强。即从噪音或干扰中提取被掩盖的语音信号。第 五，语音编码。主要用于语音数据压缩，目前已经建立了一系列语音编码的国际标准，大 量用于通信和音频处理。近年来，这5个方面都取得了不少研究成果，并且，在市场上已 2数字滤波器设计出现了一些相关的软件和硬件产品，例如，盲人阅读机、哑人语音合成器、口授打印机、语音应答机，各种会说话的仪器和玩具，以及通信和视听产品大量使用的音频压缩编码技 术。(2)图像处理数字滤波技术以成功地应用于静止图像和活动图像的恢复和增强、数据压缩、去噪音 和干扰、图像识别以及层析X射线摄影，还成功地应用于雷达、声纳、超声波和红外信号 的可见图像成像。通信在现代通信技术领域内，儿乎没有一个分支不受到数字滤波技术的影响。信源编码、信道编码、调制、多路复用、数据压缩以及自适应信道均衡等，都广泛地采用数字滤波 器，特别是在数字通信、网络通信、图像通信、多媒体通信等应用中，离开了数字滤波器,几乎是寸步难行。其中，被认为是通信技术未来发展方向的软件无线电技术，更是以数字 滤波技术为基础。(4)电视数字电视取代模拟电视已是必然趋势。高清晰度电视的普及指日可待，与之配套的视 频光盘技术已形成具有巨大市场的产业；可视电话和会议电视产品不断更新换代。视频压 缩和音频压缩技术所取得的成就和标准化工作，促成了电视领域产业的蓬勃发展，而数字 滤波器及其相关技术是视频压缩和音频压缩技术的重要基础。(5)雷达雷达信号占有的频带非常宽，数据传输速率也非常高，因而压缩数据量和降低数据传 输速率是雷达信号数字处理面临的首要问题。告诉数字器件的出现促进了雷达信号处理技 术的进步。在现代雷达系统中，数字信号处理部分是不可缺少的，因为从信号的产生、滤 波、加工到目标参数的估计和目标成像显示都离不开数字滤波技术。雷达信号的数字滤波 器是当今十分活跃的研究领域之一。(6)声纳声纳信号处理分为两大类，即有源声纳信号处理和无源声纳信号处理，有源声纳系统 涉及的许多理论和技术与雷达系统相同。例如，他们都要产生和发射脉冲式探测信号，他 们的信号处理任务都主要是对微弱的目标回波进行检测和分析，从而达到对目标进行探 测、定位、跟踪、导航、成像显示等目的，他们要应用到的主要信号处理技术包括滤波、门限比较、谱估计等。(7)生物医学信号处理数字滤波器在医学中的应用日益广泛，如对脑电图和心电图的分析、层析X射线摄影 的计算机辅助分析、胎儿心音的自适应检测等。音乐数字滤波器为音乐领域开辟了一个新局面，在对音乐信号进行编辑、合成、以及在音 乐中加入交混回响、合声等特殊效果特殊方面，数字滤波技术都显示出了强大的威力。数 字滤波器还可用于作曲、录音和播放，或对旧录音带的音质进行恢复等。(9)其他领域5数字滤波器的应用领域如此广泛，以至于想完全列举他们是根本不可能的，除了以上 3数字滤波器设计几个领域外，还有很多其他的应用领域。例如，在军事上被大量应用于导航、制导、电子 对抗、战场侦察；在电力系统中被应用于能源分布规划和自动检测；在环境保护中被应用 于对空气污染和噪声干扰的自动监测，在经济领域中被应用于股票市场预测和经济效益分 析，等等。1.3 数字滤波器的实现方法分析数字滤波器的实现，大体上有如下儿种方法：(1)在通用的微型机上用软件来实现。软件可以由使用者自己编写或使用现成的。自IEEE DSP Comm.于1979年推出第一 个信号处理软件包以来，国外的研究机构、公司也陆续推出不同语言不同用途的信号处理 软件包。这种实现方法速度较慢，多用于教学与科研。(2)用单片机来实现。目前单片机的发展速度很快，功能也很强依靠单片机的硬件环境和信号处理软件可用 于工程实际，如数字控制、医疗仪器等。(3)利用专门用于信号处理的DSP片来实现。DSP芯片较之单片机有着更为突出的优点，如内部带有乘法器、累加器，采用流水线 工作方式及并行结构，多总线，速度快，配有适于信号处理的指令等，DSP芯片的问世及 飞速发展，为信号处理技术应用于工程实际提供了可能。1.4 本章小结数字滤波器精确度高、使用灵活、可靠性高，具有模拟设备所没有的许多优点，已广 泛地应用于各个科学技术领域，例如数字电视、语音、通信、雷达、声纳、遥感、图像、生物医学以及许多工程应用领域。随着信息时代数字时代的到来，数字滤波技术已经成为 一门极其重要的学科和技术领域。以往的滤波器大多采用模拟电路技术，但是，模拟电路 技术存在很多难以解决的问题，例如，模拟电路元件对温度的敏感性，等等。而采用数字 技术则避免很多类似的难题，当然数字滤波器在其他方面也有很多突出的优点，在前面部 分已经提到，这些都是模拟技术所不能及的，所以采用数字滤波器对信号进行处理是目前 的发展方向。4数字滤波器设计第2章数字滤波器的概述数字滤波器可以用查分方程、单位取样响应以及系统函数等表示。对于研究系统的实 现方法，即它的运算结构来说，用框图表示最为直接。一个给定的输入输出关系，可以用多种不同的数字网络来实现。在不考虑量化影响 时，这些不同的实现方法是等效的；但在考虑量化影响时，这些不同的实现方法性能上就 有差异。因此，运算结构是很重要的，同一系统函数”，运算结构的不同，将会影响系 统的精度、误差、稳定性、经济性以及运算速度等许多重要性能。HR(无限冲激响应)滤波 器与FIR(有限冲激响应)滤波器在结构上有自己不同的特点，在设计时需综合考虑。2.1 数字滤波器的基本结构作为线形时不变系统的数字滤波器可以用系统函数来表示，而实现一个系统函数表达 式所表示的系统可以用两种方法：一种方法是采用计算机软件实现；另一种方法是用加法 器、乘法器、和延迟器等元件设计出专用的数字硬件系统，即硬件实现。不论软件实现还 是硬件实现，在滤波器设计过程中，由同一系统函数可以构成很多不同的运算结构。对于 无限精度的系数和变量，不同结构可能是等效的，与其输入和输出特性无关；但是在系数 和变量精度是有限的情况下，不同运算结构的性能就有很大的差异。因此，有必要对离散 时间系统的结构有一基本认识。2.2 IIR滤波器的基本结构一个数字滤波器可以用系统函数表示为:M义理-H(z)=1-2华必 k=Jz)X(z)(2-1)由这样的系统函数可以得到表示系统输入与输出关系的常系数线形差分程为：N M()=Z(一 份+3(幻(2-2)k=0 k=0可见数字滤波器的功能就是把输入序列()通过一定的运算变换成输出序列火)。不 同的运算处理方法决定了滤波器实现结构的不同。无限冲激响应滤波器的单位抽样响应 久)是无限长的，其差分方程如(2-2)式所示，是递归式的，即结构上存在着输出信号到输 入信号的反馈，其系统函数具有(2-1)式的形式，因此在z平面的有限区间(0|z 8)有 极点存在。前面已经说明，对于一个给定的线形时不变系统的系统函数，有着各种不同的等效差 分方程或网络结构。由于乘法是一种耗时运算，而每个延迟单元都要有一个存储寄存器，因此采用最少常熟乘法器和最少延迟支路的网络结构是通常的选择，以便提高运算速度和 减少存储器。然而，当需要考虑有限寄存器长度的影响时，往往也采用并非最少乘法器和 延迟单元的结构。5数字滤波器设计TIR滤波器实现的基本结构有：(l)IIR滤波器的直接型结构；优点：延迟线减少一半，变为N个，可节省寄存器或存储单元；缺点：其它缺点同直接/型。通常在实际中很少采用上述两种结构实现高阶系统，而是把高阶变成一系列不同组合 的低阶系统(一、二阶)来实现。(2)IIR滤波器的级联型结构；特点：.系统实现简单，只需一个二阶节系统通过改变输入系数即可完成；极点位置可单独调整；运算速度快(可并行进行)；各二阶网络的误差互不影响，总的误差小，对字长要求低。缺点：不能直接调整零点，因多个二阶节的零点并不是整个系统函数的零点，当需要准确的 传输零点时，级联型最合适。(3)IIR滤波器的并联型结构。优点：简化实现，用一个二阶节，通过变换系数就可实现整个系统；极、零点可单独控制、调整，调整而、如只单独调整了第i对零点，调整份限 则单独调整了第，对极点；各二阶节零、极点的搭配可互换位置，优化组合以减小运算误差；可流水线操作。缺点：二阶阶电平难控制，电平大易导致溢出，电平小则使信噪比减小。6数字滤波器设计C、串联型图2-1、IIR滤波器的基本结构2.3 FIR滤波器的基本结构FIR滤波器的单位抽样响应为有限长度，一般采用非递归形式实现。通常的FIR数 字滤波器有横截性和级联型两种。FIR滤波器实现的基本结构有：(l)FIR滤波器的横截型结构表示系统输入输出关系的差分方程可写作：N-1 y(n)-h(m)x(n-m)(2-3)m=0直接由差分方程得出的实现结构如图2-2所示：“w hjl*G1中41:加3 图2-2、横截型(直接型、卷积型)若()呈现对称特性，即此FIR滤波器具有线性相位，则可以简化加横截型结构，下 面分情况讨论：图2-3、N为奇数时线形相位FIR滤波器实现结构 图2-4、N为偶数时线性相位FIR滤波器实现结构(2)FIR滤波器的级联型结构将分解成实系数二阶因子的乘积形式：NT TH(z)=E 瓯+blkz-+b2kz-2(2-4)N=0 k=l7数字滤波器设计这时FIR滤波器可用二阶节的级联结构来实现，每个二阶节用横截型结构实现。如图 所示：x(n)Mi图2-5、FIR滤波器的级联结构这种结构的每一节控制一对零点，因而在需要控制传输零点时可以采用这种结构。2.4 数字滤波器的设计原理数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(HR)滤 波器和有限长冲激响应(FIR)滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。这种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲 激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的 方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法 和Chebyshev逼近法等等。随着MATLAB软件尤其是MATLAB的信号处理工作箱的不 断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。数字滤波器设计的基本步骤如下：(1)确定指标在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很 多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给 出幅度和相位响应。幅度指标主要以两种方式给出。第一种是绝对指标。它提供对幅度响 应函数的要求，一般应用于FIR滤波器的设计。第二种指标是相对指标。它以分贝值的形 式给出要求。在工程实际中，这种指标最受欢迎。对于相位响应指标形式，通常希望系统 在通频带中具有线性相位。运用线性相位响应指标进行滤波器设计具有如下优点：只包 含实数算法，不涉及复数运算；不存在延迟失真，只有固定数量的延迟；长度为N 的滤波器(阶数为AM),计算量为N/2数量级。因此，本文中滤波器的设计就以线性相位 FIR滤波器的设计为例。(2)逼近确定了技术指标后，就可以建立一个目标的数字滤波器模型。通常采用理想的数字滤 波器模型。之后，利用数字滤波器的设计方法，设计出一个实际滤波器模型来逼近给定的 目标。(3)性能分析和计算机仿真上两步的结果是得到以差分或系统函数或冲激响应描述的滤波器。根据这个描述就可 以分析其频率特性和相位特性，以验证设计结果是否满足指标要求；或者利用计算机仿真 实现设计的滤波器，再分析滤波结果来判断。8数字滤波器设计2.5 滤波器的性能指标我们在进行滤波器设计时，需要确定其性能指标。一般来说，滤波器的性能要求往往 以频率响应的幅度特性的允许误差来表征。以低通滤波器特性为例，频率响应有通带、过 渡带及阻带三个范围。在通带内：l-ApS|H(*|W|G)|coc 加匕)|在阻带中：|H3)|0Ast 3*|3 回c其中C为通带截止频率，g,为阻带截止频率，4为通 弓 带误差，4,为阻带误差。与模拟滤波器类似，数字滤波器按频率特性划分为低|；通、高通、带通、带阻、全通等类型，由于数字滤波器的频入、-士不，率响应是周期性的，周期为2兀。各种理想数字滤波器的幅度丁频率响应如图所示：图2-6、低通滤波器频率响应幅度特性的容限图图2-7、各种理想数字滤波器的幅度频率响应2.6 IIR数字滤波器的设计方法目前，IIR数字滤波器设计最通用的方法是借助于模拟滤波器的设计方法。模拟滤波 器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图 表供查询，因此，充分利用这些已有的资源将会给数字滤波器的设计带来很大方便，IIR 数字滤波器的设计步骤是：(1)按一定规则将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标；(2)根据转换后的技术指标设计模拟低通滤波器8；(3)在按一定规则将转换为。若所设计的数字滤波器是低通的，那么上述设计工作可以结束，若所设计的是高通、9数字滤波器设计带通或者带阻滤波器，那么还有步骤：(4)将高通、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标，然后 按上述步骤(2)设计出模拟低通滤波器H(s),再由冲击响应不变法或双线性变换将“转 换为所需的。S-Z映射的方法有：冲激响应不变法、阶跃响应不变法、双线性变换法等。下面讨 论双线性变换法。双线性变换法M是指首先把S平面压缩变换到某一中介平面S1的一条横带(宽度为2 冗T,即从-冗T到Ji T),然后再利用z=的关系把si平面上的这条横带变换到整个Z 平面。这样s平面与z平面是一一对应关系，消除了多值变换性，也就消除了频谱混叠 现象。s平面到z平面的变换可采用。=ta n(&g)(2-5)2_ 皿e丁-qz(2-6)J 万+Jh令 jCl=s,=S有：s=亘=1 _l(2-7)e 2+0 2从si平面到z平面的变换，即z=/(2-8)代入上式，得到：s(2-9)一般来说，为使模拟滤波器的某一频率与数字滤波器的任一频率有对应关系，可引Q T入代定常数 c,O=cta n(q-)(2-10)1-z_ 1则 s=cr(2-11)1+z这种S平面与Z平面间的单值映射关系就是双线性变换。有了双线性变换，模拟滤 波器的数字化只须用进行置换。2.7 FIR数字滤波器的设计方法IIR滤波器的优点是可利用模拟滤波器设计的结果，缺点是相位是非线性的，若需 要线性相位，则要用全通网络进行校正。FIR滤波器的优点是可方便地实现线性相位。FIR滤波器单位冲激响应力(力的特点：其单位冲激响应()是有限长(#卸-1),系统函数为：N7“(z)=(2-12)w=0在有限Z平面有(N-1)个零点，而它的(N-1)个极点均位于原点z=0处。10数字滤波器设计Fir滤波器线性相位的特点：如果FIR滤波器的单位抽样响应(“)为实数，而且满足以下任一条件：偶对称奇对称h(n)=由(N-1-玲其对称中心在=(N-l)/2处，则滤波器具有准确的线性相位。窗函数设计法：一般是先给定所要求的理想滤波器频率响应由Ke%)导出%()，我们知 道理想滤波器的冲击响应幻()是无限长的非因果序列，而我们要设计的是年5)是有限长 的FIR滤波器，所以要用有限长序列()来逼近无限长序列为()，设：(2-13)常用的方法是用有限长度的窗函数武)来截取电3)即:h(n)=conhdn(2-14)这里窗函数就是矩形序列凡(力，加窗以后对理想低通滤波器的频率响应将产生什么样 的影响呢?根据在时域是相乘关系,在频域则是卷积关系：H(ejco)=忆(暧)唯 Q3切2%上万(2-15)其中,航13s为矩形窗谱，(暧)是FIR滤波器频率响应.通过频域卷积过程看的幅度函数(的起伏现象,可知,加窗处理后，对理 想矩形的频率响应产生以下几点影响：使理想频率特性不连续点处边沿加宽，形成一个过渡带，其宽度等于窗的频率响应的 主瓣宽度。(2)在截止频率的两边的地方即过渡带的两边，出现最大的肩峰值，肩峰的两侧形成起伏 振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡的多少，则取决于旁瓣的多少。(3)改 变M只能改变窗谱的主瓣宽度，改变3的坐标比例以及改变的绝对值大小，但不能改 变主瓣与旁瓣的相对比例(此比例由窗函数的形状决定)o(4)对窗函数的要求a、窗谱主瓣尽可能窄，以获取较陡的过渡带；b、尽量减小窗谱的最大旁瓣的相对幅度；即能量集中于主瓣，使肩峰和波纹减小，增大阻带的衰减。频率采样法：窗函数设计法是从时域出发，把理想的支功用一定形状的窗函数截取成有限长的()，来近似理想的l,Q 增加，A(Q2)快速减小。心”3)=；，喘4,幅度衰减相当于3db衰减点。振幅平方函数的极点HA-S)Ha(S)=1+(3-2)可见，Butter worth滤波器的振幅平方函数有2N个极点，它们均匀对称地分布在|S|二0c的圆周上。考虑到系统的稳定性，知DF的系统函数是由S平面左半部分的极点(SP3,SP4,SP5)组成的，它们分别为:Sp3=ce,Sp4=-,Sp5=系统函数为：“(S%)(SS”)(S%)令得归一化的三阶BF：Ha(s)=；-八 S3+2S2+2S+1如果要还原的话，则有H(s)=-(5/Q(.)3+2(5/)2+2(5/n(.)+1(3-3)(3-4)(3-5)(3-6)3.3 切比雪夫滤波器巴特奥兹低通滤波器的幅频特性随。的增加而单调下降，当N较小时，阻带幅频特 性卜降较慢，要想使其幅频特性接近理想低通滤波器，就必须增加滤波器的阶数，这就将 导致模拟滤波器使用的原件增多，线路趋于复杂。切比雪夫滤波器的的阻带衰减特性则有 所改善。特点：误差值在规定的频段上等幅变化。巴特沃兹滤波器在通带内幅度特性是单调下降的，如果阶次一定，则在靠近截止频率 Oc处，幅度下降很多，或者说，为了使通常内的衰减足够小，需要的阶次(N)很高，为了 克服这一缺点，采用切比雪夫多项式逼近所希望f灯。)|2。切比雪夫滤波器的在通带范围内是等幅起伏的，所以同样的通带衰减，其阶 数较巴特沃兹滤波器要小。可根据需要对通带内允许的衰减量(波动范围)提出要求，如要 求波动范围小于Idbo振幅平方函数为9 14。2)=|/w(yn)|-二-(3-7)1+/)15数字滤波器设计式中Q一有效通带截止频率一与通带波纹有关的参量，大,波纹大,0 loVn(x)N阶切比雪夫多项式，定义为cos(#cos-1 x)Ixl 1W 时必(%)归1|%|1,|%,乙(%)个如图3-1,通带内-Qc,随0/0。/,|匕。)一。(迅速趋于零)当0=0时，H 1+/cos7V a rccos(O)十 2(3-8)(3-9)(3-10)N为偶数，风(/。)忆0=1,cos2(Nq)=l(n iin),(3-11)1 H-2图3-2、切比雪夫滤波器的振幅平方特性有关参数的确定：a.通带截止频率，预先给定b.由通带波纹表为。，d=20 lg 卜*=201g一J一此(Q)L iS=101g(l+/)给定通带波纹值分贝数/加)后,可求。(3-13)(3-14)16数字滤波器设计C.阶数N一由阻带的边界条件确定。(。,力事先给定)。二。,时,风(4)!(3-15)3.4 椭圆滤波器。=。,时,风(4)!(3-16)1g 317)1+V2 幺 AN SJ318)%段 C)忖1时,VN(x)=cosh(N cosh x)(3-19)得nWWi/320)ar cosh(Qr/Q)特点：幅值响应在通带和阻带内都是等波纹的，对于给定的阶数和给定的波纹要求,椭圆滤波器能获得较其它滤波器为窄的过渡带宽，就这点而言，椭圆滤波器山是最优的,其振幅平方函数为依),W.gf 1+?媪口。(3-21)式中，RM。，L)为雅可比椭圆函数，L是一个表示波纹性质的参量。L)一 一席一._!i Yo+1 a图3-3、N=5时 的特性曲线由图可见，在归一化通带内(twowi),在(o,1)间振荡，而超过以后，在IA 8间振荡。L越大，也变大。这一特点使滤波器同时在通带和阻带具有任意衰 减量。下图为典型的椭园滤波器振幅平方函数：17数字滤波器设计N为奇数图3-4、椭圆滤波器的振幅平方函数图中&和A的定义与切比雪夫滤波器相同。当。c、Qs、和人确定后，阶次N的确定方法为：确定参数人=QJQs确定参量占二/8 TaTN二型犬匹、KCkJKVTk7式中K(k)=j 炉户)i为第一类完全椭圆积分。3.5用ma tla b设计数字滤波器传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的 应用。利用MATLAB信号处理工具箱(Sign a l Processin g Toolbox)心可以快速有效的设 计由软件组成的常规数字滤波器的设计方法。给出了使用MATLAB语言进行程序设计和利 用信号处理工具箱的FDATool工具进行界面设计的详细步骤。利用MATLAB设计滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于 滤波器设计的最优化。3.6 FDATool 界面FDATool(Filter Design&An a lysis Tool)是MATLAB信号处理工具箱里专用的滤波 器设计分析工具，MATLAB6.0以上的版本还专门增加了滤波器设计工具箱(Filter Design Toolbox)o FDATool可以设计几乎所有的基本的常规滤波器，包括FIR和HR的各种设计 方法。它操作简单，方便灵活工FDATool界面总共分两大部分，一部分是Design Filter,在界面的下半部，用来设 置滤波器的设计参数，另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种 特性。Design Filter部分主要分为：18数字滤波器设计Filter Type（滤波器类型）选项,包括Lowpa ss（低通）、Highpa ss（高通）、Ba n dpa ss（带通）、Ba n dstop（带阻）和特殊的FIR滤波器。Design Method（设计方法）选项，包括HR滤波器的Butterworth（巴特沃思）法、Chebyshev Type 1（切比雪夫I型）法、Chebyshev Type H（切比雪夫H型）法、Elliptic（椭圆滤波器）法和FIR滤波器的Equiripple法、Lea st-Squa res（最小乘方）法、Win dow（窗函数）法。Filter Order（滤波器阶数）选项，定义滤波器的阶数，包括Specify Order（指定阶 数）和Min imum Order（最小阶数）。在Specify Order中填入所要设计的滤波器的阶数（N 阶滤波器，Specify 0rder=N-l）,如果选择Min imum Order则MATLAB根据所选择的滤波 器类型自动使用最小阶数。Fren quen cy Specif ica tion s选项，可以详细定义频带的各参数，包括采样频率Fs 和频带的截止频率。它的具体选项由Filter Type选项和Design Method选项决定，例如 Ba n dpa ss（带通）滤波器需要定义Fstopl（下阻带截止频率）、Fpa ssl（通带下限截止频 率）、Fpa ss2（通带上限截止频率）、Fstop2（上阻带截止频率），而Lowpa ss（低通）滤波器 只需要定义Fstopl、Fpa sslo采用窗函数设计滤波器时，由于过渡带是由窗函数的类型 和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。Ma gn itude Specifica tion s选项，可以定义幅值衰减的情况。例如设计带通滤波器 时，可以定义Wstopl（频率Fstopl处的幅值衰减）、Wpa ss（通带范围内的幅值衰减）、Wstop2（频率Fstop2处的幅值衰减）。当采用窗函数设计时，通带截止频率处的幅值衰减 固定为6db,所以不必定义。Win dow Specifica tion s选项，当选取采用窗函数设计时，该选项可定义，它包含了 各种窗函数，在通带内的衰减为6dB。Filter Design a Analysis Tool-untitled,fdaFile Edit Analysis Targets View Window Help0 B国目国 必比X 亡1 E 日四国麴哈出厂圈Q国回 eStructure:Direct-Form FIROrder:50Stable:YesSource:DesignedStore Filter.Filter Manager.品窗趣随演国Rea_ Response Type_LowpassHighpassBandpassBandstopDifferentiator _ Design Method_HR ButterworthFIR Equiripple _ Filter Order_Specify order:Minimum orderOptions-Density Factor:20Magnitude Specifications.Units:dBApass 1Astop 8019数字滤波器设计图3-5、FDATool的操作页面通过菜单选项An a lysis可以在特性区看到所设计滤波器的幅频响应、相频响应、零 极点配置和滤波器系数等各种特性。设计完成后将结果保存为*.fda文件。在设计过程中，可以对比滤波器幅频相频特性和设计要求，随时调整参数和滤波器类 型,以便得到最佳效果。其它类型的FIR滤波器和IIR滤波器也都可以使用FDATool来设 计。3.7 用Fda tool进行带通滤波器设计给定的数字带通滤波器的参数是：通带为45Hz55Hz,低截止频率为40Hz,高截止 频率为60Hz,通带内衰减不大于3db,阻带衰减大于80db。借助Ma tla b仿真，可以得到设计成不同类型的滤波器所需的阶数：表3-1、采用不同类型滤