FIR滤 波器设计 课件.docVIP

第7章 ＦIＲ滤波器设计 第六章我们介绍了无限冲激响应(IIR)滤波器得设计方法、其中最常用得由模拟滤波器转换为数字滤波器得方法为双线性变换法,因为这种方法无混叠效应,效果较好。但通过前面得例子我们瞧到,IIR数字滤波器相位特性不好(非线性,如图 6—1１、图６－1３、图6—１５　),也不易控制。然而在现代信号处理中,例如图像处理、数据传输、雷达接收以及一些要求较高得系统中对相位特性要求较为严格,这种滤波器就无能为力了、改善相位特性得方法就是采用有限冲激响应滤波器。本章首先对FIＲ滤波器原理及其使用函数作基本介绍,然后重点介绍窗函数法设计ＦIR滤波器,并对最优滤波器设计函数进行介绍。 7、1 FIR滤波器原理概述及滤波函数 7、1、１ FＩＲ滤波器原理及设计方法分类 根据第 6 章对数字滤波器得介绍,我们知道ＦIＲ滤波器得传递函数为: 　　 　　　 　 　 　 　 　 　(７-1) 可得FIR滤波器得系统差分方程为: 　　 　　 　 　 因此,ＦＩR滤波器又称为卷积滤波器。根据第 4 章中所描述得系统频率响应,FIR滤波器得频率响应表达式为: 　 　 　 　 　 　 (７—2) 信号通过ＦＩR滤波器不失真条件与(6-6)式所描述得相同,即滤波器在通带内具有恒定得幅频特性与线性相位特性。理论上可以证明(这里从略):当ＦIR滤波器得系数满足下列中心对称条件: 　 　 　　 (７-3) 时,滤波器设计在逼近平直幅频特性得同时,还能获得严格得线性相位特性。线性相位FＩR滤波器得相位滞后与群延迟在整个频带上就是相等且不变得。对于一个 Ｎ 阶得线性相位ＦIR滤波器,群延迟为常数,即滤波后得信号简单地延迟常数个时间步长。这一特性使通带频率内信号通过滤波器后仍保持原有波形形状而无相位失真、 本章主要介绍得FIR数字滤波器设计方法及 MAＴLＡB 信号处理工具箱提供得了FIR数字滤波器设计函数,见表7—1。由于篇幅所限,本章我们主要介绍窗函数法与最优化设计方法。 表7—１ ＦIR滤波器设计得主要方法 函数设计方法 说明 工具函数 窗函数法 理想滤波器加窗处理 fｉr1(单频带) , ｆiｒ2(多频带) ,　kaｉserord 最优化设计 平方误差最小化逼近理想幅频响应或Park－ＭcCｌｅllaｎ 算法产生等波纹滤波器 ｆirls , ｒemez,rｅｍｅｚord 约束最小二乘逼近 在满足最大误差限制条件下使整个频带平方误差最小化 firｃlｓ,fircls１ 升余弦函数 具有光滑、正弦过渡带得低通滤波器设计 Ｆircos 7、1、２　 FIＲ数字滤波器滤波函数 相对于ＩIR 滤波器得滤波函数,FIR数字滤波器滤波函数除了dｉmpulsｅ与ｄｓteｐ仅适用于IＩR滤波器外,其她各种函数可直接应用于FIR滤波器,只就是输入得分母多项式向量a=1。另外,MＡTLAB还提供了一个函数fｆtfｉlｔ,该函数利用效率高得基于FＦＴ算法实现对数据得滤波,该函数只适用于FIR滤波器,调用形式为: y=fｆtｆiｌt(b,x[,n］) 式中,ｂ为FIR滤波器得系数向量;x为输入数据;n为FFＴ长度,缺省时,函数选用最佳得FFT长度,y为滤波器得输出。该函数执行下面得操作: n＝length(x); y＝iｆft(fft(ｘ).＊fft(b,n)、/fｆｔ(a,ｎ)); 应注意,ｙ=fｆtfilt(b,x)等价于y=fiｌtｅｒ(b,a,x)。 ７。2　 FIR滤波器得窗函数设计 7.２.1 　窗函数得基本原理 ＦIR滤波器设计得主要任务就是根据给定得性能指标确定滤波器得系数b,即系统单位脉冲序列ｈ(n),它就是一个有限长序列。 ＦIR滤波器得理想频率响应,可写成复数形式得Fｏurier级数形式: 　　 　　 　　 　 　 　 　 　 (7-4) 式中,hd(n)就是对应得单位脉冲响应序列。这说明滤波器得频率响应与单位脉冲响应互为Ｆouｒｉeｒ变换对。因此其单位脉冲响应可由下式求得, 　 　　 　 　　 　 　 　(7—5) 求得序列后,通过z变换,可得到 　 　 　 　　 　 　 　 　 　 (7—6) 注意,这里为无限长序列,因此就是物理上不可实现得。如何变成物理上可实现呢？一个自然得想法就是只取其中得某些项,即只截取中得一部分,比如n＝０,…,N—１,Ｎ为正整数。这种处理相当于将,n＝-∞~∞与函数w(ｎ)相乘,w(n)具有下列形式: w(n)相当于一个