FIR滤波器毕业设计论文非常标准的格式

1、摘 要21世纪是数字化的时代，纵观当代通信的发展趋势，已成为引领通信变革的主潮流。通信是在数字化浪潮的背景下，在计算机技术的应用和信息技术的发展的结果。数字信号滤波器在各种数字信号处理中发挥着重要的作用，数字信号设计一直是数字信号处理领域的重要研究课题。近年来，数字信号技术在我国也得到迅速发展，不论是在科学技术研究，还是在开发等发面，其应用越来越广泛，并取得了丰硕的成果。本文主要介绍如何用窗函数法和雷米兹交换法设计FIR滤波器的的具体步骤与方法，以及相关数字信号处理的一些具体算法，并在MATLAB环境下进行仿真。根据仿真运行的结果来说明各项运行指标均达到设计要求。分析和比较两种设计方法性，以及

2、它们分别所适用的范围，通过设计表达这两种方法的运算简单、精度高、设计过程简单易行，适合于工程设计。关键词：FIR数字滤波器 线性相位 MATLAB仿真 窗函数 雷米兹法AbstractThe 21st century is the era of digital, looking back at the development trend of contemporary communications, has become the main changes leading communications trend. This is the communication in the digital

3、 wave of context, the application of computer technology and the development of information technology results. Digital signal filter in a variety of digital signal processing plays an important role, digital signal design digital signal processing has been an important research topic in the field.

4、In recent years, Digital Signal Technology in China has been developing rapidly, both in scientific and technological research and development is failing in its application more and more widely, and achieved fruitful results.Therefore, this artcle introduced the FIR filter may realize the strict lin

5、ear phase under the window function and Remez function, designs gronp of filters coefficients ,vses least squares method to optimize these coefficients .in the foundation which smallest two rides to unify the Lagrange law first to restrain the belt is smallest two rides transfers asks the condition

6、extreme value ,introduces Lagrange to leave the Lagrange function while the structure ,then carries on the solution Finally the full use data analysis carries on the simulation realization under the MATLAB environment .Thas may know,restraint least squares method designs the filter has the algorithm

7、 simply.Key words: FIR Digital Filter Minimum MATLAB simulation Linear phase Window function Remez function目 录摘要IAbstract.II目录.III绪论.11 数字滤波器的简介.21.1 数字滤波器的介绍.21.2 数字滤波器的原理.21.3 数字滤波器的设计.41.3.1 数字滤波器的设计过程.41.3.2 数字滤波器的设计方法.52 FIR数字滤波器的基本结构.62.1 FIR滤波器的基本结构.62.2 最大误差最小化准则.93 线性相位的FIR数字滤波器.123.1 线性相位的

8、概念.123.2 线性相位滤波器.133.3 线性相位FIR数字滤波器的设计方法.153.3.1 雷米兹交换法设计FIR数字滤波器.153.3.2 FIR数字滤波器的线性规划设计 184 线性相位FIR滤波器的仿真设计.204.1 信号处理工具箱中的最优设计函数.204.2 MATLAB设计FIR数字滤波器的方法.214.3 线性相位FIR滤波器的仿真设计224.4 线性相位FIR滤波器的窗函数法的仿真设计234.5 线性相位FIR滤波器的雷米兹交换法的仿真设计.254.6 结果分析28结论.30致谢.31参考文献.32绪 论随着信息时代和数字世界的到来，数字信号处理己成为当今一门极其重要的学

9、科和技术领域，数字信号处理在通信、雷达、军事、航空航天、语音、图像、自动控制、医疗和家用电器等众多领域得到了广泛的应用。数字滤波器是数字信号处理的重要基础，在对信号的滤波、检测及参数的估计等信号应用中，数字滤波器是使用最为广泛的一种线性系统，在研究信号的时候，首先必须考虑噪声的干扰对信号的传输影响，噪声是一切干扰信号的泛指，有的仅希望最大限度地去除噪声而已，有的希望在去除噪声时能让滤波器具有线形相位，有的则是强调滤波的实时性，在设计时针对一些情况，制定有针对性的滤波器，来改善信号的质量。本文共分为四章，前三章为设计数字滤波器的基础知识，数字滤波器从功能上分为低通滤波器(LPF)、高通滤波器(H

10、PF)、带通滤波器(BPF)、带阻滤波器(BSF)；然而数字滤波器的实现方法有很多种，例如FPGA、DSP、MATLAB等，本文是利用MATLAB进行数字滤波器的仿真设计。最后一章是用三种方法设计FIR数字滤波器。我们可以得出第一类线性相位滤波器可以用于实现低通、高通、带通和带阻等各种滤波特性；通过数据和图形分析得出在相同的滤波器抽样响应长度下，如果在一个频带内赋予了大的加权，那么这个频带内将获得大的衰减。因此，通过调整加权值，可得到不同的衰减，在通带和阻带都具有较好的性能。1 数字滤波器的简介1.1 数字滤波器的介绍数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR

11、）滤波器和有限长冲激响应（FIR）滤波器。IIR系统易取得较好的通带和阻带衰减特性，一般要求H（z）阶次要高，即M要大。FIR系统有自己突出的优点：系统总是稳定的，易实现线性相位，允许设计多通带（或多阻带）滤波器，后两项都是IIR系统不易实现的。FIR数字滤波器的设计方法有多种，如窗函数设计法、频率采样法和Chebyshev逼近法等。随着MATLAB软件尤其是MATLAB的信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计实现最优化。1.2 数字滤波器的原理数字滤波器可分为FIR（有限脉冲响应）和IIR（无限脉冲响应）两种。IIR滤波器的系统函数是两个Z的多项式

12、的有理分式，而FIR滤波器的分母为1，即只有一个分子多项式。数字滤波器的理想幅频特性如图1-1所示。在0到的全部频段上，其幅值为1的区域为通带，其余为阻带，即其幅值为0。根据wc1和wc2取值不同可分为4种类型：（1） 低通滤波器，当wc1=0时；（2） 高通滤波器，当wc2=时；（3） 带通滤波器，当wc1及wc2如图1-1所示时；（4） 带阻滤波器，当0,wc1及wc2,1区间幅度为1，wc1,wc2区间幅度为0时。图1-1 理想幅频特性有些情况下，还对滤波器的相位特性提出要求，理想的是线性相位特性，即相移与频率成线性关系。实际的滤波器不可能完全实现理想幅频特性，必有一定误差，因此要规定适

13、当的指标。低通滤波器在0，的通带区，幅频特性会在1附近波动；在1的阻带区，幅频特性不会真等于零是一个大于零的值；在, 之间，为过渡区；这三个与理想特性的不同点就构成了滤波器的指标体系。即通带频率和通带波动，阻带频率和阻带衰减。在许多情况下，人们习惯用分贝为单位，定义通带波动为（分贝）阻带衰减为（分贝）。 (1-1) （1-1）对于带通滤波器，范围为,；对于带阻滤波器，应表为。其他复杂形状的预期特性通常也可由若干理想的幅频特性叠合构成。FIR数字滤波器最大的优点是容易设计成线性相位特性，并且具有稳定性。1.3 数字滤波器的设计(1) 确定技术指标 在设计一个滤波器之前，必须首先根据工程实际的需要

14、确定滤波器的技术指标。在很多实际应用中，数字滤波器常被用来实现选频操作。因此指标的形式一般在频域中给出幅度和相位响应。幅度指标主要以两种方式给出。第一种是绝对指标。他提供对幅度响应函数的要求，一般应用于FIR滤波器的设计。第二种指标是相对指标。他以分贝值的形式给出要求。(2) 逼近 确定了技术指标后，就可以建立一个目标的数字滤波器模型（通常采用理想的数字滤波器模型）。之后，利用数字滤波器的设计方法（窗函数法、频率采样法等），设计出一个实际滤波器模型来逼近给定的目标。 (3) 性能分析和计算机仿真 上两步的结果是得到以误差或系统函数或冲激响应描述的滤波器。根据这个描述就可以分析其频率特性和相位特

15、性，以验证设计结果是否满足指标要求；或者利用计算机仿真实现设计的滤波器，再分析滤波结果。 数字滤波器的设计过程（1）按照实际需要，确定滤波器的性能要求。通常是在频域中给定数字滤波器的性能要求。通带截止频率在通带内幅度响应以的误差接近于1，即（1-2）为阻带起始频率，在阻带内幅度响应小于的误差接近于零，即 （1-3） 为了使逼近理想低通滤波器的方法成为可能，还必须提供一宽度为的不为零的过滤频带。在这个频带内，幅度响应从通带平滑地下落到阻带。这里（）指的是数字域频率，或者说是沿单位圆周的相角变化。相位特性受到稳定性和因果性要求的限制（即要求系统函数的极点必须位于单位圆内部）。（2）寻找满足预定性能

16、要求的离散时间线性系统。IIR函数是的有理函数。FIR滤波器的系统函数是的多项式。这样，滤波器的设计问题，变成了一个数字逼近问题，即用一个因果稳定系统函数去逼近给定的性能要求，以确定滤波器系数。（3）用有限精度的运算来实现设计的系统。包括选择运算结构，滤波器的系数，输入变量，中间变量，和输出变量。（4）通过模拟，验证所设计的系统是否符合给定性能要求。根据这步的结果决定是否对第二步和第三步作修改，以满足技术的要求。 数字滤波器的设计方法设计FIR数字滤波器的基本方法有窗函数法、频率取样法和等波动最佳逼近法，这些方法主要是针对选频型滤波器（低通、高通、带通和带阻滤波器）的设计，这种滤波器的设计指标

17、是类似的，典型的指标为通带波动和阻带衰减。在FIR数字滤波器的设计中，还会涉及微分器和希尔伯特变换器之类的系统，这类非选频型滤波器的设计也遵循以上方法，更完善的设计则是基于任意频域指标的。数字滤波器的设计方法很多，大多数方法都在计算的复杂性和满足设计滤波器的指标两个问题间取得折衷。FIR滤波器的设计法方法可以分为以下几种：（1）频率采样法，（2）窗函数法，（3）雷米兹交替算法等。2 FIR数字滤波器2.1 FIR滤波器的基本结构在讨论任何一种滤波器时，都要着重分析其系统函数，FIR滤波器的系统函数为： 。从该系统函数可看出，FIR滤波器有以下特点：1) 系统的单位冲激响应h(n)在有限个n值处

18、不为零；2) 系统函数H(z)在|z|0处收敛，极点全部在z=0处(稳定系统)；3) 结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中(例如频率抽样结构)也包含有反馈的递归部分。FIR滤波器基本分为以下几种类型(1) 横截型（卷积型、直接型）a) 一般FIR滤波器的横截型(直接型、卷积型)结构：若给定差分方程为： 。 则可以直接由差分方程得出FIR滤波器结构如下图所示：图2-1 差分方程得出FIR滤波器结构图这就是FIR滤波器的横截型结构，又称直接型或卷积型结构。b.线性相位FIR滤波器的横截型结构若h(n)呈现对称特性，即此FIR滤波器具有线性相位，则可以简化成横截型结构，下面分情

19、况讨论：N为奇数时线性相位FIR滤波器实现结构如图2-2所示：图2-2 奇数时线性相位FIR滤波器实现结构如图N为偶数时线性相位FIR滤波器实现结构如图2-3所示图2-3 偶数时线性相位FIR滤波器实现结构图我们知道IIR滤波器的优点是可利用模拟滤波器设计的结果，缺点是相位是非线性的，若需要线性相位，则要用全通网络进行校正，比较麻烦，而FIR滤波器的优点是可以方便地实现线性相位。(2) 级联型将H(z)分解为若干个实系数一阶或二阶因子相乘： （2-1）实现结构如图2-4所示：图2-4 级联型实现结构如图该结构图中有2L=M个延迟器，2L+1=M+1个乘法器，2L=M个加法器。分析H(z)及结构

20、图可以得出级联型的特点：每个基本节控制一对零点，便于控制滤波器的传输零点。系数比直接型多，所需的乘法运算多。(3) 频率取样型 若FIR滤波器的冲激响应为有限长(N点)序列h(n)，则有如图2-5所示的关系：图2-5 FIR滤波器的冲激响应为有限长序列关系图因此，对h(n)可以利用DFT得到H(k)，然后利用内插公式： 来表示系统函数，这就为FIR滤波器提供了另外一种结构：频率抽样结构其中级联的第一部分为： 这是一个梳状滤波器，它滤掉了频率 及其各次谐波。级联的第二部分为N个一阶网络并联而成，第k个一阶网络为：，2L+1=M+1个乘法器，2L=M个加法器它在单位圆上有一个极点： 由上叙的理论分

21、析基础可以得到FIR滤波器的频率抽样结构。FIR滤波器的频率抽样结构如图2-6所示：图2-6 FIR滤波器的频率抽样结构如图频率抽样结构的特点是它的系数H(k)就是滤波器在 处的响应，因此控制滤波器的频率响应很方便。频率抽样结构存在问题的问题是：在有限长情况下，系数量化后极点不能和零点抵消，使FIR系统不稳定。(4) 快速卷积结构若FIR滤波器的单位冲激响应h(n)是一个N1点有限长序列，输入x(n)是一个N2点有限长序列，那么输出y(n)是x(n)与h(n)的线性卷积，它是一个LN1+N2-1点的有限长序列。 我们知道，将x(n)补上LN2个零值点，将h(n)补上LN1个零值点，然后进行L点

22、圆周卷积，就可以代替原x(n)与h(n)的线性卷积。而圆周卷积可以用DFT和IDFT的方法来计算，这样我们得到FIR滤波器的快速卷积结构：图2-7 FIR滤波器的快速卷积结构图这里DFT和IDFT都将采用快速傅里叶变换算法，当N1和N2足够长时，比直接计算线性卷积要快得多。2.2 最大误差最小化准则在滤波器的设计中，通常情况下通带和阻带的误差要求是不一样的。等波纹最大误差最小化准则就是通过对通带和阻带使用不同的加权函数，实现在不同频段（通常指的是通带和阻带）的加权误差最大值相同，从而实现其最大误差在满足性能指标的条件下达到最小值。（1） 加权切比雪夫逼近误差及交错定理线性相位FIR根据单位抽样

23、响应h(n)的奇偶对称性以及h(n)的长度N的奇偶性，总共可以分为四种类型。尽管如此，FIR的频率响应依然可以采用如下的统一形式来表示： (2-2) 其中：为幅度函数，它是一个可正可负的纯实数。利用三角恒等式知识和交错定理可得： (2-3)在FIR的四种类型中：加权切比雪夫误差公式可定义为： (2-4)其中： 为加权误差，为逼近误差加权函数，为理想幅度函数，为实际滤波器幅度函数。将(2)式带入(3)式并令：经推导可得： （2-5）（2-4）公式也是最终的加权切比雪夫逼近误差函数公式。那么线性相位FIR DF的加权切比雪夫等波纹逼近问题实际上就是求解表达式的问题，从而使得在实行逼近的频率范围内的

24、最大绝对值达到最小。在此定义该最小值表达式为： （2-6）A为实行逼近的频带。为了求解(25)式，ParksMcCLLan把逼近理论中的交错点定理应用到滤波器设计中，从而得出了如下的交错定理：设是个r余弦函数的线性组合，即： （2-7）A是内所研究的一个闭子集， 是A上的一个连续函数，则在A内能够最佳并且唯一地逼近的充要条件是：加权切比雪夫逼近误差函数在A中至少r+1个极值点，即在A中存在共r+1个频率点，各频率点均满足关系式： （2-8）3 线性相位FIR数字滤波器3.1 线性相位的概念数字滤波器的频率响应可以由幅值表示。如果数字滤波器的相位响应满足条件，那么称该数字滤波器是线性相位。式中是

25、一个常量。如果是正数，那么该系统延迟信号；否则就是一个超前系统，此时相位是频率的线性方程，可以归为的形式，其中斜率为，截距。假设正弦曲线的频率为，即周期为，其中 因为一个周期对应于，所以相位改变对应于延迟。同样 , 如果，那么对应的延迟时间为 ，可见延迟时间和频率无关。因此如果数字滤波器满足线性相位条件，那么所有频率分成的延迟量是一样的。这意味着滤波后的输出只是输入信号的一个简单的延迟信号。从另一方面来讲，如果滤波器不具有线性相位，那么输入信号的不同频率成分延迟量是不同的，这将会导致输出信号的失真，在实际设计中通常要回避这种情况。如果系统的频率响应为.那么其幅值响应是1,相位响应是,表明是线性

26、相位的。假设,那么它只是起到延迟的作用，其相位响应只画出了区间内的相位，这也导致了原本是线性的相位出现了弯折。下面考虑给出了频率响应的情况，给定 （3-1）式中是实数。既然是实数，所以它只会影响输入信号的幅值大小，而仅仅使输入信号产生相移。如果。相位，那么系统是一个线性想一系统。如果，相位，那么这种情况下，延迟时间和频率是有关的。从上面给出的线性相位的定义的角度来说，该系统不是严格意义上的线性相位系统。但是可以将上面的式子写成。这样一来中括号里的函数就变成了线性相位，此时波纹不再失真，负号只要将波纹沿纵轴反转即可。但是如果会改变符号，那么波纹就可能失真。只改变水平轴附近的符号，即阻带内的符号，

27、此时阻带内的信号极大地衰减。所以信号通过一个频率响应系统时，通带内信号没有产生任何失真。这样的系统也常常称为线性相位系统。这里顺便要指出的是模拟滤波器不可能有线性相位特性，只可能在很小的一个频带内近似地认为是线性相位。3.2 线性相位滤波器可以很容易设计出满足线性相位特性的FIR滤波器，这使得FIR滤波器得到广泛应用。如前所述，FIR滤波器肯定是稳定的。另一方面，线性相位IIR滤波器的设计就不那么简单了，而且通常只能使它在一定的频带频率范围内满足线性相位性质。虽然如此，但是IIR也有比FIR优越的方面，那就是当IIR滤波器和FIR滤波器具有相同幅值响应时，前者所有的系数少很多。假设一个因果FI

28、R滤波器同（如下式中滤波器的长度为，为滤波器系数。以下表明滤波器的线性相位特性可以通过滤波器系数成偶数对称或奇数对称的。系数偶对称意味着，系数奇对称意味着。令，此时就有六个系数，假设这些系数是偶数对称的，如下面关系式，其中所以其脉冲响应关于点n对称。 （3-2）对于M为奇数的一般情况，有以下结论：（M是奇数，系数偶对称） （3-3）注意到求和号后每项都是实数，所以式（3-3）和（3-2）具有相同形式，所以FIR-1滤波器是线性相位的。令，此时就有七个系数，而且是对称的。此时的频率响应变为 （3-4）（M为偶数，奇数偶对称） （3-5）那么当奇数奇对称时，对于，假设系数如图所示，此时符号发生改变

29、，这导致了分子中的余弦函数被替代。下面给出了通用公式：（M为奇数，系数奇对称） （3-6）因为M是奇数，所以系数的中心应该是本身的负数，所以必须等于0。对于，通用公式应该为 （3-7） FIR-III和FIR-IV滤波器和等式的形式不同，所以不是线性相位的。这种滤波器的相位响应为。相应响应的滤波器常常称为广泛线性相位。 （3-8）或是这一情况的特例，此时的滤波器就变成前面讨论的线性相位了。广泛线性相位滤波器在许多场合有广泛应用，包括窄带滤波器以及通讯信号的解调。这些滤波器有固定的群延迟或时间延迟，其定义如下： （3-9）3.3 线性相位FIR数字滤波器的设计方法最优设计就是充分利用技术指标来进

30、行设计。误差容限设计低通滤波器，要求在频带内以最大误差逼近1，在频带内以最大误差逼近零位。我们将一要求表示为加权逼近误差函数的形式。并且使用最大误差最小化准则将其描述为切比雪夫逼近问题。最优线性相位FIR数字滤波器的设计就是要设法求得切比雪夫逼近的最优解的滤波器的系数。人们在寻求最优化设计上做了大量的工作。1970年发表了非线性方程的方法求解切比雪夫逼近的最优解。1971年出现了更好的拉格朗日内插多项式求解法。到了1973年又找到了雷米兹算法求解加权误差的方法。非线性方程解法及多项式内插法，之适用于设计那些误差极值点数目为最大可能性的滤波器，也即最多波纹滤波器。同时由于N，是固定的，所以滤波器

31、的频带边缘不能预先规定，需在最后的解求得以后，才能计算出来。它可用来设计任何最优（最大误差最小化）线性相位FIR滤波器。此外，目前还有线性规划技术设计方法，下面对雷米兹算法及线性规划技术设计法分别加以介绍。 雷米兹交换法设计FIR数字滤波器雷米兹交换算法是为了在N固定时，能控制和的需要而产生的。前面已将最优线性相位FIR滤波器的设计问题描述为切比雪夫逼近问题，逼近函数是r个独立的余弦函数之和。交替定理给出了加权逼近误差函数的一组必要充分条件，使逼近成为所需频率响应的唯一最好逼近。基于交替定理的最优FIR滤波器的设计程序的主要步骤：（1） 输入部分：规定所需要的频率响应为，加权函数和滤波器的长度

32、N。（2） 用公式表示逼近问题，即形成。（3） 用雷米兹多次交换算法，求逼近问题的解。（4） 计算滤波器的单位取样响应。第一步设计滤波器算法，表达所要求设计的滤波器的类型和必须满足的性能要求。第二步在前面切比雪夫加权逼近已提及。第三步用雷米兹算法求逼近问题的解。需要指出的是，在整个程序中，雷米兹算法是作为一个子程序出现的，在调用该子程序以前，主程序已完成了以下几点。 读输入数据（滤波器的技术规格） 根据滤波器的类型和长度确定了逼近函数cos的个数r 用密集的格点代替了频率区间。确定了两格点间的距离为因而总格点数等于（N+1）格点密度/2，并给所有下标格点赋上了标称频率值。调用了子程序EFF和W

33、ATE计算各格点频率上所要求的函数值和加权函数值.根据四种情况统一的公式将、，变成了。根据交替定理，建立了一组等间隔的极值频率初始值。等波纹的误差曲线是在多次迭代中形成的。雷米兹迭代计算是从（r+1）个极值频率的初始假设值开始的。第一次迭代的（r+1）个极值频率是按等间隔假设的，这些频率位于区间内，并且由于和是固定的，所以中的某一频率，即。假定这些频率点上的误差函数的数值为，其符号为正负交替。这就是说根据问题的原始要求，对于给定的一组极值频率，需要求以下方程式中 （3-10）计算出以后，确定出r个极值频率上的值 （3-11）利用拉格朗日内插公式得出式中 （3-12）要注意也可以内插到 （3-1

34、3）求出的内插值以后，在根据公式 (3-14)在密集的频率组上计算值。若在改频率组的某些频率上，则选择（r+1）个新的频率作为可能的极值点，新的频率就选所得误差曲线上那些峰值点频率，然后重新计算这些点上的误差函数值。作为这次迭代寻找新的r+1个误差最大点的比较标准，看在这些频率上计算值。为求那些峰值点，应在通带和阻带上把频率分点取得更密一些，以便在这些细分点上搜索出峰值点。如果在任一次迭代中的极值点多于（r+1）个，就保留值最大的（r+1）个频率作为下一次迭代的假设极值点。随着迭代次数的增加，极值频率的位置逐次向最佳位置上调整，一直重复到与其前一个值相同为止，最后一次迭代的结果对应于问题的解。

35、此时误差曲线上每个格点频率处的误差值都满足，r+1个极值频率处的，并且具有正负交替的符号，这标志着加权切比雪夫等波纹逼近误差已经形成，最佳逼近找到了。在每次迭代中都是将，定位极点频率中的两个频率。若最后得到的,并且规定加权函数为 （3-15）则在通带内值逼近，在阻带内值逼近于。由雷米兹迭代算法最后所求得值是要求的最小值。如果要求的和值是已知的，则计算滤波器时可以固定，改变，重复进行以上迭代计算，直至得出要求的和值。过度区宽度出现的局部极小值，曲线上的这些点相当于超波纹滤波器。极小值之间的所有点对应于按交替定理为最优的滤波器。第四步求滤波器的单位取样响应。在区间的个等间隔的频率点上计算值，利用I

36、DFT求得，最后，根据四种线性相位滤波器的不同情况，求出于a(n)相应的单位取样响应h(n)。于是满足预定要求的最优线性相位FIR滤波器被唯一确定。根据上述原理，已编出一个通用的线性相位FIR滤波器的设计程序。此程序的内部接口，用于设计多频带通带阻滤波器（包括低通、高通、带通和带阻滤波器），微分器，希尔伯特变换器。 FIR数字滤波器的线性规划设计FIR线性相位滤波器是对所有，误差函数的最大值极小化的滤波器。可以用一组线性不等式描述这个最大值极小化问题 （3-16）这里表示最大误差，F为要进行逼近的频带中的一组密集，为r个余弦函数的线性组合。将上式写成线性规划的形式为 （3-17）可以用线性规划

37、技术解上述不等式组，但要比雷米兹法慢得多。用线性规划法设计对频域和时域都有约束的滤波器是唯一简单的方法。4 线性相位FIR滤波器的仿真设计4.1 信号处理工具箱中的最优设计函数 要构成完整的最优等波动滤波器设计的程序，除了最小最大波动公式外，还有许多实际问题要考虑。（1） 滤波器长度N（或阶数M=N-1）如何确定？（2） 极值数目的确定。最优等波动滤波器的误差函数在上有（L+2）或（L+3）个极值。大多数等波动滤波器有（L+2）个极值。但是对于某些的组合，可能得到有（L+3）和极值的滤波器。此处L=floor(N-1)/2)。（3） 如何建立进行频率修正的算法，在程序中可以自动进行反复的迭代修

38、正，直到要求的精度为止。把最优等波动算法和上述的工程问题组合在一起，才能形成完整的设计算法，交替定理保证最大最小逼近问题的解存在并且惟一，但它并没有说明最优得到这个解，既不知道阶数N(或L)，也不知道极值频率和最大误差。滤波器技术指标中给出了，因此需要设定N的值。凯泽提出了以个简单的公式来逼近N。（4-1）这种算法首先猜设（L+2）个极值频率，估计这些频率上的最大误差，接着按给定的各点，你和一个L阶多项式，然后在一个较细的网格上确定局部极大误差及其极值频率，由这些新频率点你和出一个新的L阶多项式，重复以上过程。一直进式系数d(n),并最后算出滤波器脉冲响应h(n)。由于N是近似的，最大误差可能

39、不等于。如果出现这种情况，需要增加N或减小N，再次用remez算法确定一个新的。重复此过程。这样就得到了等波动滤波器。在MATLAB中，实现算法的函数为remez，它最常用的句法为：它由几种调用形式：设计一个M阶FIR数字滤波器，它的频率响应在数组f和A中给定。长度为N 的数组h是返回的滤波器系数。数组f单位为，即。其中包含各边缘频率，这些频率必须以递增次序排列，从0开始，到1结束。数组A为各定频率上预期的幅度响应。F和A的长度必须相等且为偶数（频率成对出现），形成依次的关心频带和不关心频带。每个关心频带中所用的权函数等于1，这说明在每个频带中的容限是相同的。：与上述情形相似，数组的长度是f的

40、一半，它是每个关心频带的权函数。：与第一种情形类似，不过加了一个滤波器类型作为输入变元。此函数响应的设计数字希尔伯特变换器或数字微分器。：与上述类似，只是由数组指定每个频带的加权函数。正如在讨论Parks_McClellan算法中指出的，为了使用程序remez,首先必须假设滤波器的阶数，得到数组h中的滤波器系数后，还必须检查最小阻带衰减，并与给定的比较，然后增加（或减少）滤波器的阶数。重读此过程直到得到期望的。4.2 雷米兹法设计FIR数字滤波器的方法信号处理工具箱采用remez算法实现线性相位FIR数字滤波器的等波纹最佳一致逼近设计。与其他设计法相比，其优点是，设计指标相同时，使滤波器阶数最

41、低；或阶数相同时，使通带最平坦，阻带最小衰减最大；通带和阻带均为等波纹形式，最适合设计片段常数特性的滤波器。其调用格式如下：b=remez(N,f,m,w,ftype)其中，w和ftype可默认。b为滤波器系数向量，调用参数N,f,m的含义与函数fir2中类同，但这里有一点不同，期望逼近的频幅响应值位于f(k)与f(k+1)（k为奇数）之间的频段上，而f(k+1)与f(k+2)之间为无关区。w为加权向量，其长度为f的一半。W(k)为对m中第k个常数片段的逼近精度加权值，w值越大逼近精度越高。rtype用于指定滤波器类型。Remezord函数用于估算FIR数字滤波器的等波纹最佳一致逼近设计的最低

42、阶数N，从而使滤波器在满足指标的前提下造价最低。基本调用格式如下：N,fo,mo,w=remezord(f,m,dev,Fs)其返回参数供remez函数使用。设计的滤波器可以满足由参数f,m,dev和Fs指定的指标。F和m与remez中所用的类似，这里f可以是模拟频率(Hz)或归一化数字频率，但必须以0开始，以Fs/2(用归一化频率时为1)结束，而且其中省略了0和Fs/2两个频点。Fs为采样频率，省略时默认为2Hz。dev为各逼近频段允许的幅频响应偏差（波纹振幅）。remez函数可直接调用remezord返回的参数，使用格式如下：b=remez(N,fo,mo,w)。4.3 线性相位FIR滤波

43、器的仿真设计设计一个FIR线性相位系统h=3,-1,-5,4,6,4,-5,-1,3, 程序及运行结果如下：functionA,w,type,tao=amplres(h)h=3,-1,-5,4,6,4,-5,-1,3;M=length(h);tao=(M-1)/2;L=floor(M-1)/2);n=1:L+1;w=0:500*2*pi/500;if all(abs(h(n)-h(M-n+1)1e-10) A=2*h(n)*cos(M+1)/2-n)*w)-mod(M,2)*h(L+1) type=2-mod(M,2)elseif all(abs(h(n)+h(M-n+1)1e-10)&(h(

44、L+1)*mod(M,2)=0) A=2*h(n)*sin(M+1)/2-n)*w); type=4-mod(M,2);else errorendh=3,-1,-5,4,6,4,-5,-1,3;M=length(h);n=0:M-1;subplot(2,1,1),stem(n,h);ylabel(h(n) );xlabel(n);subplot(2,1,2),plot(w/pi,A);ylable(A);xlabel(pi);MATLAB命令窗显示：type=1图4-1 原始与输出信号4.4 线性相位FIR滤波器的窗函数法的仿真设计这种方法也称为傅立叶级数法。其设计是在时域进行的，先用傅氏反变

45、换求出理想数字滤波器的单位抽样响应hd(n)，然后时域移位并加时间窗w(n)对其截断，从而求得FIR 滤波器的单位抽样响应h(n)；（4-2）在设计过程中，将无限长序列变为有限长是通过时域加矩形窗乘积实现数据的截断的。时域乘积对应了频域卷积，从而对频响特征发生的改变。常见的窗函数有：矩形窗、三角形(Bartlertt)窗、汉宁(Hanning)窗。海明(Hamming)窗、布拉克曼(Blackman)窗、凯泽(kaiser)窗等，下面介绍几种常用的窗函数： （4-3）设计一个窗函数线性相位FIR滤波器，程序及运行结果如下：omegac=0.37;N=81;m=(N-1)/2;n=0:2*m+1

46、0;h=omegac/pi*sinc(omegac*(n-m)/pi);w=ones(1,N) zeros(1,length(n)-N);hd=h.*w;omega=-pi:2*pi/300:pi;Hd=freqz(hd,1,omega);plot(omega,abs(Hd);图4-2 输出信号4.5 线性相位FIR滤波器的雷米兹交换法的仿真设计用函数firs和remez设计一个50阶FIR多通滤波器，滤波器理想频率响应为：f=0 0.1 0.15 0.25 0.3 0.4 0.45 0.55 0.6 0.7 0.75 0.85 0.9 1 , a=1 1 0 0 1 1 0 0 1 1 0

47、0 1 1;将设计的滤波器的幅频响应和理想滤波器幅频响应进行比较。%Samp clf;n=20; f=0 0.3 0.4 1; a=1 1 0 0; b=firls(n,f,a); h,w1=freqz(b); bb=remez(n,f,a); hh,w2=freqz(bb); figure(1)plot(w1/pi,abs(h),w2/pi,abs(hh),r:); xlabel(归一化频率);ylabel(振幅);legend(firls,remez); grid on;figure(2)fs=50;t=0:1/fs:2; f1=5;f2=15; x1=sin(2*pi*f1*t)+8.*

48、cos(2*pi*f2*t);subplot(2,1,1),plot(t,x1)title(原始信号)y1=filter(b,1,x1);y2=filter(bb,1,x1);subplot(2,1,2),plot(t,y1,t,y2,r:)legend(firls,remez); title(输出信号)xlabel(时间/s) 图4-3 幅频特性图 4-4 原始与输出信号%Samp clf;n=50; f=0 0.1 0.15 0.25 0.3 0.4 0.45 0.55 0.6 0.7 0.75 0.85 0.9 1;a=1 1 0 0 1 1 0 0 1 1 0 0 1 1 b=firl

49、s(n, f,a); h,w1=freqz(b); bb=remez(n, f,a); hh,w2=freqz(bb); figure(1)plot (w1/pi,abs(h),b.,w2/pi, abs(hh),r:,f,a,m);xlabel(归一化频率);ylabel(振幅);legend(firls,remez,理想特性); figure(2)impz(bb,1),title (脉冲响应) xlabel(样本数);ylabel(幅度)图 4-5 幅频特性 图 4-6 脉冲响应4.6 结果分析首先由图4-1可见，这是一个第一类线性相位滤波器。滤波器系数N为9（奇数），且h(n)=h(N-1-n),幅度特性关于对称，在和处可以取任何值；可以用于实现低通、高通、带通和带阻等各种滤波特性。程序运行结果见图4-1和图4-3。比较两种设计方法所设计的滤波器的幅频响应可见，用firs设计的滤波器