FIR数字滤波器的优化设计

1、FIR数字滤波器的优化设计学生：赵然然，信息系指导教师：李永全，信息系摘要FIR数字滤波器在保证幅度特性满足技术要求的同时，又很容易做到有严格的线性相位特性，故FIR数字滤波器用得很多，因而研究FIR数字滤波器的优化设计具有重要的理论意义。本文介绍FIR数字滤波器的原理，并在此基础上利用加权平方误差最小准则、等波纹切比雪夫逼近和最小二乘法对FIR滤波器进行优化设计。在每种设计方法中，都对原理、设计要求、设计步骤、选择特性、最优处理以及相关工具函数等进行了详细介绍。同时，采用MATLAB软件对FIR数字滤波器进行设计，简化了设计中繁琐的计算，而通过滤波前后信号的频谱图的对比，分析不同滤波器的滤波

2、效果，从而更快的得到优化设计方法。关键词FIR数字滤波器 优化设计方法加权平方误差最小准则 等波纹切比雪夫逼近 最小二乘法 Optimum Design for FIR Digital FilterStudent:Zhao Ranran, Information DepartmentInstructor：Li Yongquan, Information DepartmentAbstractFIR digital filter ensure the range characteristics and meet the technical requirements at the same time

3、,and it is easy to do a strict linear phase characteristic,so it is used much more.It is meaningful to research the optimum design for FIR digital filter. This paper introduced the principle of the FIR filter,and used weighed least2squares error criterion, equivalent ripple Chebyshev approximation a

4、nd the least squares method to carry on the optimized design to the FIR filter. Every kind of design method are explained in detail about its requirements, steps, select properties and optimal treatment. Designing the FIR filter by Matlab can simplify the complicated computation in simulation,and co

5、mparing the signals spectrum viewers and the sound files which have been generated to analyzing the filtering effect of different digital filters,and gaining the optimization techniques faster.KeywordsFIR digital filter optimization techniques weighed least2squares error criterionequivalent ripple C

6、hebyshev approximationleast squares method前言数字滤波器(digital filter)是由数字乘法器、加法器和延时单元组成的一种装置。MATLAB是第四代计算机语言,是目前公认的国际上最流行的科学与工程计算的软件工具。随着科技和新兴技术的发展，数字滤波在图像处理、语音识别和模式识别等 数字信号处理中占有越来越重要的地位。与模拟滤波器相比，数字滤波器可以满足滤波器幅度和相位特性的严格要求，可以克服模拟滤波器所无法克服的电压漂移、温度漂移和噪声等问题。有限冲激响应（FIR）滤波器可以保证严格的线性相位。同时由于其实现结构主要是非递归的，因此 FIR 滤波

7、器可以稳定工作。 FIR 滤波器被广泛用于各类数字信号处理系统中实现卷积、相关、自适应滤波、正交插值等处理,对于非实时系统和低速采样系统,FIR 滤波器的运算可在 CPU 或 DSP 处理器上采用软件实现。而由于它越来越重要的地位，本文将对其进行一些简要的介绍，并通过 MATLAB来实现 FIR 数字滤波器的优化，使FIR 数字滤波器设计精度提高，简化其设计过程，而通过这一系列的介绍和实验，让我们进一步的理解和掌握 FIR 数字滤波器。研究的目的和意义FIR数字滤波器的优化设计1 研究的目的和意义1.1 研究的目的 与模拟滤波器相比,数字滤波器除了具有数字信号处理固有优点外,还有滤波精度高、稳

8、定性好、灵活性强等优点。FIR数字滤波器是数字信号处理的基础，用来对信号进行过滤，检测与参数估计等处理，利用数字滤波器可改变信号中所含频率分量的相对比例或滤除某些频率分量，使其达到所需要的效果。它在通信、语言，图像、自动控制、雷达，军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。尤其在图像处理、数据压缩等方面取得了令人瞩目的进展和成就。FIR滤波器可以得到严格的线性相位,但它的传递函数的极点固定在原点,只能通过改变零点位置来改变性能,为了达到高的选择性,必须用较高的阶数,对于同样的滤波器设计指标,FIR滤波器要求的阶数可能比IIR滤波器高510倍。通过本次设计，完成FIR数字滤波器设

9、计的优化，达到提高设计精度 ,简化编程过程,降低滤波器成本等目的。1.2 研究的意义 大部分工程技术领域中都会涉及到信号的处理问题，其信号表现形式有电、磁、机械以及热、光、声等。信号处理的目的一般是对信号进行分析、变换、综合、估值与识别等。如何在较强的噪声背景下提取出真正的信号或信号的特征，并将其应用于工程实际是信号处理的首要任务。数字信号处理中一个非常重要且应用普遍的技术就是数字滤波。数字滤波器有FIR数字滤波器和IIR数字滤波器，IIR数字滤波器的设计方法是利用模拟滤波器成熟的理论及设计图表进行设计的，因而保留了一些典型模拟滤波器优良的幅度特性，但设计中只考虑了幅度特性，没考虑相位特性，所

10、设计的滤波器一般是某种确定的非线性相位特性。为了得到线性相位特性，对IIR滤波器必须另外加相位校正网络，使滤波器设计变得复杂，成本也高，又难以得到严格的线性相位特性。而FIR滤波器在保证幅度特性，满足技术要求的同时，很容易做到有严格的线性相位特性，故FIR数字滤波器用得很多。数字滤波技术是数字信号分析、处理技术的重要分支。无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传输是至第1页（共30页）FIR数字滤波器的优化设计关重要的。在所有的电子系统中，使用最多技术最复杂的要算数字滤波器了。数字滤波器的优劣直接决定产品的优劣，因而研究FIR数字滤波器的优化

11、设计具有重要的理论意义。2 FIR数字滤波器的基本概念2.1 滤波器简介在通信这个大背景之下，信号在传输的过程中会受到很多因素的影响，例如多径衰落，噪音干扰，远近效应，甚至只是简单的电平衰落。这些复杂的环境，会严重干涉信号的正常传输。对某些系统而言，如果不对自然传输的信号加以处理，会使得有效信息产生较大的误差。所以，如何保证信息传输的有效性，成为信息处理的一个关键问题。滤波器就是这个环节中非常重要的一个组成部分。它的作用就是阻隔干扰信号，无用信号，使之相对有用信号更大幅度的衰减，最大可能实现滤波，达到系统要求的理论理想效果。2.1.1 滤波器的实现滤波器可以是软件也可以是硬件。除去滤波的特别作

12、用不谈，其实它也只被看做一个信号处理的系统。站在系统的角度，在设计中需要考虑很多因素，例如稳定性等。设计的方法，也基于信号与系统的理论前提。伴随数字系统的优势逐日显露，实际应用中的滤波器大多也是数字滤波器，它可以用以下函数来概括：（1）简单来看，数字滤波器的功能，就是把输入序列x (n)通过一定的函数运算变换成输出序列y (n)，满足相关映射关系。实际设计中，需要先明确滤波器的使用要求，然后确定滤波器的技术指标，从而开始具体设计。归纳起来，滤波器实现的过程包括四个一般步骤：确定逼近函数：设计能够满足理想技术指标的转换函数。实现方程：将转换函数对应为滤波网络中对数字序列进行运算的方程或相应的系数

13、向量组。第2页（共30页）FIR数字滤波器的基本概念研究缺陷：研究实际中的非理想因素的影响，如采样值是否超出存储的有限字长等问题。并做出相应的改进或转变，以求最大程度的逼近理想效果。产品实现：使用一定的硬件设备，例如DSP处理器甚至是普通计算机，能够满足专用运算来构建滤波器，它的优点是可以进行实时的处理。当然也可以直接通过软件和专用的数字信号处理芯片来实现。这与实际背景有关。本文的研究讨论只涉及到滤波器逼近函数的相关内容，从理论上认识和探讨一系列逼近函数和对应理想滤波器的差距。利用MATLAB仿真平台实现设计分析，达到认识和总结的目的。2.1.2 滤波器的分类滤波器种类繁多，不同的角度也会造成

14、分类方法的不同。可以从功能、实现方法或者是设计方法上来分类等。大方向上，滤波器可以分为模拟滤波器和数字滤波器。在模拟滤波器的设计过程中可能更多的涉及到电阻，电容，电感等元器件的应用；他们分别针对模拟系统和数字系统。实际中数字滤波器的应用相对广泛。从滤波效果来分：能通过的信号频段可以分为低通（Lowpass）、高通（Highpass）、带通(Bandpass)和带阻滤波器（Bandstop）四种。从硬件组成和应用来看，滤波器种类繁多，比如有源滤波器，自适应滤波器、复数滤波器以及多维滤波器等。总的来讲还可以分为两大类，就是经典滤波器和现代滤波器。如果信号和噪声处于不同的频带，则滤波器只要具有较好的

15、滤波选频特性就可以达到理想的效果。但是如果信号噪声并没有从频率上区分开来，传统的选频滤波器就无能为力了。另一类通过对随机信号的统计特性进行滤波的现代滤波器就可以满足这种要求。利用自相关函数和功率谱估计出来的信号可能比原来具有更高的性噪比。维纳滤波器、卡尔曼滤波器、自适应滤波器都属于现代滤波器。本文介绍的滤波器无论是原理还是应用上都只涉及到传统经典滤波器。面对选频背景，主要从滤波的算法上寻优。从实现方法上来看，数字滤波器还可以分为IIR和FIR，即无限单位冲激响应滤波器和有限单位冲激响应滤波器。FIR滤波器具有严格的线性相位特征，且始终稳定，应用相对比较广泛。而IIR滤波器则用在相位要求不是很严

16、格的场合。此处简要的分析一下IIR相对于FIR的优势。IIR滤波器相对而言，明显的特点就是响应无限。同第3页（共30页）FIR数字滤波器的优化设计时它的极点可以处于任意位置，会引起系统的不稳定。并且它有反馈回路。它的优势在于所要求的阶数比FIR的低，同时可以借助模拟滤波器的研究成果，直接查表查图，设计过程非常方便。2.1.3 采样定理在讨论的滤波器认识之前，简要介绍采样定理。重要性在于它是数字信号理论研究方法的理论基础。自然产生的信号一般都是连续的，如果选用数字系统（如计算机）来处理信号，就会涉及到数字信号和模拟信号的相互转化。在模数转化（A/D）之前，必须选择采样的周期和量化的电平数。选择错

17、误会产生严重误差并丢失有用信息。采样定理给出了正确选择采样周期T的准则。对频率为的连续正弦信号采样时，采样定理要求采样频率应大于的两倍： （2）连续信号可以看做是由一个或多个正弦信号组成，假设其最高频率是，如果采样频率是最高频率的两倍或两倍以上，则正弦信号可以通过等间隔的样值来唯一表示。当采样序列通过一个对高于的正弦信号有抑制作用的系统时，原始连续信号就可以由采样序列重建。由上可知最小采样频率是，采样频率的一半被定义为奈奎斯特频率（Nyquist frequency），也被叫做折叠频率，它的最小值其实也就是。奈奎斯特频率。如果抽样信号的频率小于最小抽样频率但大于奈奎斯特频率，经采样成为离散信号

18、后，就有小于奈奎斯特的信号混入，混入频率的出现就仿佛有两个信号存在，一个频率是，另一个是。在实际应用中，通常选择采样频率为信号频率最大值的4倍。 利用采样定理正确提取数字信号成为数字系统分析和实践的基础理论。2.2 FIR数字滤波器特点归纳FIR数字滤波器有以下几个特点： FIR属于数字滤波器，所以它具有数字滤波器所具备的一切通有优点。例如精度高，灵活性可靠性强，能够达到高性能的技术指标，便于时分复用，和大规模集成生产。第4页（共30页）FIR数字滤波器的基本概念在结构上，它属于非递归结构，没有输入到输出的反馈。FIR滤波器响应有限，利于编成，计算延迟相对较小。可以做到任意幅度的频响特性。因为

19、有限，所以可以采用快速傅里叶变换。大大提高运算效率。故而实践起来方便有效，成为一种非常常见的数字滤波器。FIR滤波器具有非常好的线性相位。FIR的极点只能在圆点，所以只能通过调节零点来调整滤波器性能。因为FIR滤波器响应有限，且极点在单位圆内，所以一定稳定。在一定的时延的背景下，任意的信号都可以变成有限长序列，所以总能用因果系统进行实现。FIR的阶数要求可能是IIR滤波器的510倍，在这一点上增加了一定延迟。假设FIR数字滤波器的单位脉冲响应h(n)为一个N点序列，0nN-1，则滤波器的系统函数为：H(z)=n=0N-1h(n)z-n (3)说明该系统有N-1阶极点在z=0处，有N-1个零点位

20、于z平面。FIR数字滤波器有直接型、级联型、线性相关型、快速卷积型和频率采样型等几种结构。2.3 FIR数字滤波器的网络结构2.3.1 直接型结构直接型结构又被称为卷积型或横截性。其系统的差分表达式为： (4)在结构中，输出y(n)是每个沿着这条链的抽头信号由相应的系数（脉冲响应）加权，然后将所得乘积相加得到。设计方便简单，所用的乘法器相对较少，使用最为普遍。直接型对其他参数的控制，例如，零点位置等，没有其他结构有优势。根据差分方程可以直接画出系统的结构信号流图如图2-1所示，该结构称为直接型结构，也称为卷积型、横截型结构。第5页（共30页）FIR数字滤波器的优化设计图1 FIR滤波器的直接型

21、结构直接型结构的优点：简单直观，乘法运算量较少。缺点：调整零点较难。 2.3.2 级联型结构 当需要控制滤波器的传输零点时，可将H(z)分解为实系数二阶因子的乘积形式：Hz=n=0N-1h(n)z-n=k=1N-12(0k+1kz-1+2kz-2) (5)上式中，Hz为h(n)的z变换，0k，1k，2k为实数，N-12表示取的整数部分。当N为偶数时，N-1为奇数，除有N-12个阶子系统外，还有一个一阶子系统，图2给出了N为奇数时的级联型结构。图2 FIR滤波器的级联型结构该结构的优点：调整零点比直接型方便。缺点：Hz中的系数比直接型多，因而需要的乘法器多。当Hz的阶次高时，也不易分解。 2.3

22、.3 线性相关型结构FIR滤波器的线性相位结构有偶对称和奇对称，不论h(n)为偶对称还是奇对称都有： 第6页（共30页）FIR数字滤波器的基本概念当N为偶数时，系统函数为： å-=-±= Hz=n=0(N2)-1h(n)z-n+z-N-1-n (6)当N为奇数时，系统函数为： Hz=n=0(N2)-1h(n)z-n+z-N-1-n+hN-12z-N-12 (7)对这两种情况，都可以用FIR直接型实现，其信号流图如图3所示。(a)N为偶数(b)N为奇数图3 线性相位型结构 这种结构在本质上是直接型，但乘法次数比直接型省了一半。 2.3.4 频率采样型结构系统函数在单位圆上作N

23、等分取样就是单位取样相应h(n)的离散傅里叶变换Hk。Hk与系统函数之间的关系可用内插公式表示：H(z)=1-zNk=0N-1Hk1-WN-kz-1 (8)这个公式为FIR数字滤波器提供了另外一种结构频率采样型结构。频率采样型结构是一种用系数将滤波器参数化的一种实现结构，这种结构由两部分级联构成即第7页（共30页）FIR数字滤波器的优化设计H(z)=1NHczk=0N-1Hkz (9)其中，级联的第一部分Hcz=1-z-N为由N节延时单元构成的梳状滤波器，是由N节延时单元组成的全零点网络。 令Hcz=1-z-N=0，即zkN=1=ej2k，则有 zk=ej2kN ，k=0，1，2，N-1 (1

24、0)即在单位圆上有N个等间隔的零点，其频率响应特性为 Hcej=1-e-jN=2je-jN2sinN2 (11)因而幅度响应为 Hcej=2sinN2 (12)级联的第二部分为 k=0N-1Hkz=k=0N-1Hk1-WN-kz-1 (13)它是由N个一阶网络并联组成，每个一阶网络在单位圆上有 一个极点 zpk=WN-k=ej2kN (14)即该网络在频率为=2kN处响应为无穷大，故等效于谐振频率为2kN的谐振器。一阶网络的极点正好与梳状滤波器的一个零点相抵消，从而使得在该频率上的频率响应等于H(k)。这样N个一阶网络的极点就和梳状滤波器的N个零点相互抵消，从而在N个频域采样点上的频率响应就分

25、别等于N个H(k)值。FIR数字滤波器的频率采样型结构如图4所示。图4 FIR滤波器的频率采样结构 第8页（共30页）FIR数字滤波器的基本概念频率采样结构的优点：(1)在频率采样点k，H(ejk)，只要调整H(k)就能有效地调整频响特性;(2)只要h(n)长度N相同，对于任何频响，其梳状滤波器部分和N个一阶网络部分完全相同，只是各支路增益H(k)不同。相同部分便于标准化、模块化。 缺点:(1)寄存器长度都是有限的，零、级点可能不能正好抵消，造成系统不稳 ;(2)当N很大时，其结构很复杂，需要的乘法器和延时单元很多。2.4 FIR数字滤波器的线性相关如果FIR滤波器的单位脉冲响应h(n)为实数

26、序列，并满足一定的对称性（可以是偶对称也可以是奇对称）。并且对称中心在处，则该滤波器就能保证具有准确的线性相位。面对取值N的不同，h(n)可以被划分为四种情况，它们分别对应了四种线性相位的滤波器。详细介绍如下。1.时域特点FIR的时域特点可以归结为以下两类：群时延定义为，是一个常数，所以在此将第一类和第二类线性相位的特征统称为恒定群时延特性。偶对称和奇对称相比，除了具有线性相位特征之外，还有的固定相移，它们共同构成了正交变换网络。2.频域特点将线性相位具体分为四类，并在频域上做出分析将特点归纳如下：第9页（共30页）FIR数字滤波器的优化设计综合以上特点可以将四种线性相位的FIR滤波器的特点总

27、结如下表1：表1 四种线性相位FIR滤波器的特性第10页（共30页）FIR数字滤波器的优化设计3 FIR数字滤波器的优化设计3.1 FIR数字滤波器的优化设计准则3.1.1 均方误差最小准则该准则下是希望使得误差能量最小。用表示理想的频率响应，表示实际设计中所得到的滤波响应，并且以来表示频率响应误差： 所以均方误差为：（15）将用冲击响应来表示可以得到：（16）将冲击带入可以得到：（17）从此式中可以观察得出：等式右边第二项的大小和设计值无关，只取决与给定的特性。并且要想让取得最小值，就必须使第一项求和值最小，即：也就是说满足：（18）上面这项表达式中刚好是矩形窗的设计结果。所以，满足均方误差

28、最小准则的是矩形窗。矩形窗虽然过渡带最窄，但是根据前面分析，由于吉布斯效应（Gibbs）效应，窗谱的肩峰过大，造成所设计出的滤波器通带起伏不均匀，而阻带的衰减则过小，并不能满足设计要求。3.1.2 最大误差最小化准则该准则又被称作为加权切比雪夫等波纹逼近或雷米兹算法。在该准则中采用的方法为误差函数加权的方法。它的主要思想是希望达到不同频带（比如通带和阻带）加权误差的最大值相等的。第11页（共30页）具体分析，假设加权函数为，得到加权逼近误差函数可以定义为：。由于不同，所以值也可不同，根据不同频带的公差要求的松紧程度合理的选取不同加权值。只需使得各个频带上的即最大值要求一致即可。设计线性相位的加

29、权切比雪夫等波纹逼近问题还可以抽象为一组系数，使得在逼近的各个频带上，的最大值达到最小。对此帕克斯（Parks）和麦克莱伦（McClellan）引进了逼近理论的一个定理，得出交错定理。交错定理是切比雪夫最优逼近算法基础。定理明确指出，最优线性相位FIR滤波器加权逼近函数至少应该有（r+1）个极值（r是用于逼近函数的余弦函数的个数。具体而言的极值包括了以下2种情况：的极值点（大多数情况的极值同时也是的极值）单有极值点（不属于）两种极值点的和就是极值点的最大数目。定理中还提到，在不同频带上解决逼近问题，误差函数可以在每个频率的端点上得到一个极点，这些点一般都不是的极点，当然和除外，此时可能有极点。

30、比如在第1种线性相位中，最多有极值点r+2=(N+5)/2。同理可以推导出四种线性相位的FIR滤波器余弦数目（r）及极点数目以及的极值点数关系归纳在下表中：表2 四种线性相位参数特点类型极值1. N为奇数，h(n)偶对称2. N为偶数，h(n)偶对称3. N为奇数，h(n)奇对称4. N为偶数，h(n)奇对称对于低通响应来讲，在第1种情况中用r+2=(N+5)/2个极值实现的滤波器比交错原理要求的r+1=(N+3)/2个极值还多一个，故把这类滤波器叫做最多波纹滤波器。第12页（共30页）在雷米兹算法（remez）中可以设计任何最优线性相位FIR滤波器。但是在设计之前需要预先知道极值点数的最大值

31、数目。目前，它是一种最为实用的最优化算法。总结起来就是给定单位抽样响应长度N、带通和阻带的截止频率进行设计。算法由加权切比雪夫逼近来描述。逼近函数是由r个独立的预先函数之和，再利用交错定理给出最大误差最小化最优逼进等波纹逼近。3.2 基于加权平方误差最小准则的FIR数字滤波器的优化设计3.2.1 基于均方误差最小准则的FIR 数字滤波器的优化设计FIR数字滤波器的优化设计问题 ,是在满足给定要求下求出 FIR 数字滤波器的冲激响应 ,使h ( n)的宽度最小(或者说过渡带最窄) 。常用的均方误差最小准则,是以误差的能量最小为判据。若以 E( j)表示逼近误差谱,即E(j) = HD (j) -

32、 H(j) (19)HD (j)为给定滤波器的频率响应, H( j)是根据所寻找的 h ( n)求得的频率响应 Hj=n=0N-1h(n）e-jn (20)均方误差为: 2=1202HD（ej）- H（ej）2d=1202E(j)2d （21）均方误差最小准则就是使2最小。按矩形窗口法进行设计所产生的均方误差比任何其他设计方法所产生的都小,但矩形窗口的旁瓣电平太平,致使阻带衰减太小,常常不能满实际要求。因此,均方误差最小准则不适于用作 FIR 数字滤波器的优化设计,采用 “最大误差最小准则” ,即:minmax| E( j) | A , N 个抽样值 (22)用改变 N 个抽样值(或 N 个

33、h ( n)值) ,使在A 域内最大的绝对误差最小,这个 A 域包括通带和阻带。为了统一使用最大误差最小准则,采用加权函数的形式,即min0WE2d (23)(23)式中 W () 为预先给定的误差加权函数。采用随机抽样最小二乘法进行设计,将 FIR数字滤波器看成一个线性系统,从而可通过辩识系统参数的方法来设计滤波器的系数。第13页（共30页）3.2.2 FIR数字滤波器优化设计的算法实现设计 FIR 数字滤波器的算法如下:(1)根据设计要求确定 HD() 、H () 、 滤波器的类型(系数奇或偶对称)及滤波器的阶数 N ;(2)选定采样点个数 L ,置参数初值(0) = 0及递推最小二乘算法

34、的协方差矩阵初值 P (0) =0 ,置 t = 1 ;(3)按输入输出数据产生方法产生 x ( t )及 y( t ) ;( t ) = ( t - 1) +Pt-1x(t)1+xtPt-1x(t)yt-xt(t-1) (24)P(t)=P(t-1)-P(t-1)x(t)xtPt-11+xtPt-1x(t) (25)(4)若 t < L ,则 t = t + 1 并返回(3) ,否则到(5) ;(5)根据(4)(20)式计算 h (0) , h (1) , , h( N - 1) ;(6) H( Z) = h (0) + h (1) z-1+ h (2) z-2+ h ( N - 1)

35、 z-(N-1)，结束。3.3 基于等波纹切比雪夫逼近准则的最优化设计3.3.1 等波纹切比雪夫逼近准则在滤波器的设计中 ,通常情况下通带和阻带误差要求是不一样的。等波纹切比雪夫逼近准则就是通过对通带和阻带使用不同的加权函数 ,实现在不同频段(通常指的是通带和阻带)的加权误差最大值相同 ,从而实现其最大误差在满足性能指标的条件下达到最小值。3.3.2 加权切比雪夫逼近误差及交错定理线性相位FIR 数字滤波器根据单位抽样响应 h ( n) 的奇偶对称性以及 h( n) 的长度 N 的奇偶性 ,总共可以分为四种类型。尽管如此 ,FIR数字滤波器的频率响应依然可以采用如下的统一形式来表示:H（ej）

36、=e-jN-1/2ej/2kHd (26)其中: k 0 ,1 , Hd() 为幅度函数,它是一个可正可负的纯实数。 利用三角恒等式知识和交错定理可得:第14页（共30页）H() = Q() P() (27)在FIR数字滤波器的四种类型中:Q() , P() 的表达式可参阅文献。 则加权切比雪夫误差公式可定义为:E() = W() Hd - H() (28)其中: E() 为加权误差, W() 为逼近误差加权函数, Hd为理想幅度函数, H() 为实际滤波器幅度函数。 将(27)式代入(28)式并令:W () = W() Q () , Hd =Hd/ Q()经推导可得:E() = W () H

37、() - P() (29)(29)式也是最终的加权切比雪夫逼近误差函数公式。 那么线性相位 FIR数字滤波器的加权切比雪夫等波纹逼近问题实际上就是求解 P() 表达式的问题,从而使得在实行逼近的频率范围内 E() 的最大绝对值达到最小。 在此定义该最小值表达式为:E() = minmax | E(ej) | (30)A 为实行逼近的频带。 为了求解(30)式, Parks -McClellan把逼近理论中的交错点定应用到滤波器设计中,从而得出了如下的交错定理:设 P() 是 r个余弦函数的线性组合,即：P() =n=0r-1a(n)cosn (31)A 是0内所研究的一个闭子集, Hd是 A

38、上的一个连续函数,则 P() 在 A 内能够最佳，并且唯一地逼近 Hd的充要条件是:加权切比雪夫逼近误差函数 E() 在A中至少存在r +1个极值点,即在 A中存在1<2 < L <r+1共 r + 1个频率点,各频率点均满足关系式:E(i) = - E(i+1) = ±E() (32)| E(i) | = maxA E() i= 1 ,2 r (33)3.3.3 Remez 算法Remez算法是由 Parks和McClellan等人在1972年推导出来的。 它是将FIR数字滤波器的五个参数( N ,1 , 2 ,p , s) 中的 N , p , s和1/2固定,

39、而视1 (或2) 为变量的一种迭代方法。 尽管 Herrmann 等人在1971年也推导出来了将 N和1 , 2固定,而将p , s设为变量的另一种迭代方法,但由于前一种算法最灵活且最有效。 因而成为最优化设计的主要方法。 该方法求解过程为:第15页（共30页）FIR数字滤波器的优化设计(1)求解:首先在滤波器的通带和阻带内等间隔地取 r + 1个频率点k ( k = 0 ,1 r)作为交错点的初始值, 然后利用交错定理中P() =n=0r-1a(n)cosn的表达及的解析式求解出满足下式的值: WHd-P() =-1k k = 0 ,1 , L , r (34)(2)求解 P() :利用求解

40、出来的值和预先假定的r +1个频率点求出 P(i) (其中 i = 0 ,1r - 1)的值,然后根据拉格朗日插值公式求出 P() 的最终表达式。(3)求解 E() :将求得的 P() 代入下式判定其是否满足不等式| E() | =| W () Hd() - P() | < 的要求。 若满足要求,则说明已经获得了最优解;若在某些频率点不满足要求,则需要将这些频率点作为新的极值点重新计算,经过反复的迭代直到在所有的频率点上都满足不等式的要求。 这时的值就是最终所求的值。 这样便获得了最佳逼近。3.3.4 基于等波纹切比雪夫逼近准则的FIR数字滤波器的最优化设计步骤(1)给出所需的频率响应

41、Hd , 加权函数W() 以及滤波器的单位抽样响应 h ( n) 的长度 N。(2)由(1)中给定的参数来形成所需的 Hd, W , P() 的表达式。(3)根据 Remez算法,求解逼近问题。(4)根据求得的 P() 表达式,利用傅立叶逆变换计算出单位抽样响应 h(n) 的表达式即可获解。3.4 基于最小二乘法的 FIR 数字滤波器的优化设计3.4.1 设计思想RLS 思想是选择使得各次估计误差平方和为最小即Q=t=1net2=t=1nYt-'Xt2。换言之,Q 函数最小二乘估计就是使误差函数 Q 为最小那个， 这个量可以理解为误差能量.因此,最小二乘估计使误差能量为最小.结合 RL

42、S 对 FIR 滤波器进行优化设计,在频域计算用 H（ej）表示实际得到滤波器频响,用Hd（ej）表示要求频响,以 E（ej）表示频响误差,即 E（ej）=Hd（ej）- H（ej） （35）则误差能量：E（ej）2=Hd（ej）- H（ej）2 （36）均方误差： e2=12-Hd（ej）- H（ej）2d （37）第16页（共30页）设计目的是选择一组 h(n）F-1H（ej）使得e2最小.先将公式(35)中Hd（ej）和H（ej）分别用它们冲激响应表示Hdej=-hd(n)e-jn，Hej=n=0N-1h(n）e-jn。由于用 FIR 滤波器来逼近,故 h(n)长度是有限的.将它们代入(

43、35)可得Eej=n=0N-1hdn-h(n）e-jn+其他nhd(n)e-jn (38)按照帕塞瓦公式有e2=n=0N-1hdn-h(n）2+其他nhdn2 (39)由此式可看出,等式右边第二个求和项只取决于给定特性hdn，它和设计值 h（n）无关，故是一个常数，要使e2最小,就必须使第一项求和式最小,即希望hdn-h(n）=0，0nN-1. 在这一条件下 ，就有e2=min(e2) ，也就是说要满足 h （n）hdn 0nN-10 其他n，这个式子恰好是矩形窗的结果。3.4.2 推导低通、高通、带通、带阻 h（n）表达式(1)低通 h（n）表达式假设Hdej=1 0c0 c , 窗函数Rn

44、=1 -N-12nN-120 其他n 则 h（n）sin(n-N-12)c(n-N-12) 其中 N 值由矩形窗过渡带宽决定，N4V=2Vf(2) 高通 h(n)表达式 h（n）l- c n=N-12-sinc(n-N-12)(n-N-12) nN-12 其中 N值为 N4V=2Vf , 为过渡带宽。（3）带通 h(n)表达式 h（n）sin2(n-N-12)(n-N-12)-sin1(n-N-12)(n-N-12) 其中 N 值仍由过渡带宽决定.第17页（共30页）(4)带阻 h(n)表达式h（n）1+1+2 n=N-12sin1n-N-12-sin2n-N-12n-N-12 n=N-123

45、.4.3 FIR 数字滤波器幅度特性推导由于 h（n）hdn n，所以根据 FT 性质，时域乘积，在频域卷积。因此实际 FIR 滤波器幅度函数为 H（）H（ej）为H（）12-Hd（）WR-d 下面结合几个关键频率点,根据上式说明卷积过程：(1)当 0 时,响应 H（0）可以近似看作 从 到 WR全部积分面积，我们用 H（0）进行归一化。(2)当 =c时，Hd（）正好与WR-一半重叠,因此卷积结果正好等于 一半,即H(c)H(0)=0.5。(3)当 =c-2N 时，响应值 H(c-2N ) ，整个WR（）主瓣在Hd（）通带以内，因此卷积得到最大值H(c-2N )=max，出现正肩峰。(4)当

46、=c+2N 时，响应值 H(c+2N ) ， WR（）主瓣全部在通带Hd（）之外,通带内旁瓣负面积大于正面积,因此卷积达到最小值H(c+2N )min，出现负肩峰 。(5)当 >c+2N 时，WR（）左尾旁瓣扫过通带 ，因此Hd（）围绕零点波动。(6)c+2N<<c-2N 时，WR（）主瓣和左右旁瓣扫过通带Hd（），所以围绕 1 波动。4 用MATLAB来实现FIR滤波器的优化设计4.1 MATLAB简介MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术

47、计算语言和交互式环境，主要包括Matlab和Simulink两大部分。第18页（共30页）用MATLAB来实现FIR数字滤波器的优化设计4.1.1 MATLAB基本功能MATLAB一种工程计算的高级语言，由美国的MathWorks公司发布，主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple并称为三大数学软

48、件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连 接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中加入了对C,FORTRAN,C+,JAVA支持，可以直接调用，此外，许多的MATLAB爱好者都编

49、写了一些经典的程序，用户可以直接进行下载就可以用。4.1.2 MATLAB的特点MATLAB最重要的特点就是它的易扩展性。每个MATLAB用户都可成为对其有贡献的作者之一，并可创造自己的应用程序。在MATLAB近几年的发展过程中，有许多科学家、数字家和工程师等开发了一些新的有价值的应用程序，所有程序完全不需要使用底层代码来编写，从而极大的促进了MATLAB的发展。4.1.3 MATLAB的优势 工作平台编程环境十分友好 编程语言简单易用 数据的计算处理能力十分强大 图像处理能力强大 模块集合工具箱应用广泛 程序的接口和发布平台很实用第19页（共30页）FIR数字滤波器的优化设计 可以开发用户界

50、面4.2 优化设计的MATLAB实现在优化设计的MATLAB实现中，程序中经常使用remez函数，这种函数的使用方法为： b=remez(n,f,a,w,ftype)n为待设计滤波器的阶数；f是一个向量，它是一个0到1的正数a是一个向量，指定频率段的幅度值；w对应于各个频段的加权值函数的返回值b是设计出的滤波器的系数组成的一个长度为n+1的向量4.2.1 利用Remez函数设计等波纹低通滤波器设计要求：通带截频0.5，阻带截频0.6 , 采样频率2000Hz阻带衰减大于等于40dB，通带波纹0.1710和阻带波纹0.01程序参见附录中的1-(1)利用Remez函数设计等波纹低通滤波器图5 等波

51、纹低通滤波器的增益响应从参考程序及图5以得到所设计出滤波器的参数如下：滤波器的采样频率为2000Hz，滤波器的阶数为22滤波器的通带截频0.5，阻带截频0.6，过渡带宽均为0.1第20页（共30页）用MATLAB来实现FIR数字滤波器的优化设计阻带衰减为40dB，通带波纹为0.1710，阻带波纹为0.01对比设计要求与所设计出滤波器的参数可知，其各项参数均满足设计指标，所设计出的滤波器即为设计所要求的滤波器。图6 信号滤波前的时域图和频域图图7 信号滤波后的时域图和频域图从图6和图7的图像中可以看到：输入信号是由两个不同频率的正弦信号叠加而成，信号频域图中位于滤波器通带内的频率分量保留了下来，

52、位于滤波器阻带内的频率分第21页（共30页）FIR数字滤波器的优化设量被滤除，滤波器的效果符合设计要求。4.2.2 利用Remez函数设计等波纹带通滤波器设计要求：通带截频0.3、0.6，阻带截频0.2、0.7 阻带衰减大于等于40dB通带波纹0.1710和阻带波纹0.01采样频率2000Hz程序参见附录中的1-（2）利用Remez函数设计等波纹带通滤波器图8 等波纹带通滤波器的增益响应从参考程序及图8可以得到所设计出滤波器的参数如下：滤波器的采样频率为2000Hz，滤波器的阶数为22通带截频0.3、0.6，阻带截频0.2、0.7，过渡带宽均为0.1阻带衰减为40dB，通带波纹为0.1710，

53、阻带波纹为0.01对比设计要求与所设计出滤波器的参数可知，其各项参数均满足设计指标，所设计出的滤波器即为设计所要求的滤波器。第22页（共30页）用MATLAB来实现FIR数字滤波器的优化设计图9 信号滤波前的时域图和频域图图10 信号滤波后的时域图和频域图从图9和图10的图像中可以看到：输入信号是由四个不同频率的正弦信号叠加而成，信号频域图中位于滤波器通带内的频率分量保留了下来，位于滤波器阻带内的频率分量被滤除，滤波器的效果符合设计要求。4.2.3 利用Remez函数设计等波纹带阻滤波器第23页（共30页）FIR数字滤波器的优化设计设计要求：阻带截频0.3、0.6，通带截频0.2、0.7 阻带衰减大于等于15dB通带波纹0.01和阻带波纹0.1710采样频率2000Hz程序参见附录二中的1-（3）利用Remez函数设计等波纹