有限脉冲响应数字滤波器.ppt

数字滤波器（FIR DF）的设计,有限脉冲响应,FIR DF的定义：如果一个DF的输出y(n)，仅取决于有限个过去的和现在的输入x(n)，则称这种DF为FIR DF。FIR DF的系统转移函数为：h(n)各个样点值与滤波器的各个系数对应相等。,8.1 线性相位FIR数字滤波器的条件和特点,FIR DF有以下特点：1、h(n)是有限长的 它永远是稳定的 2、如果对h(n)提出一些约束条件，很容易使H（Z）具有线性相位。,一、线性相位条件FIR DF的系统函数为：令 代入，得：将 表示成 其中 称为幅度特性，为可正可负的实函数 为相位特性,如果相位()满足：()=-,为常数 则称 具有线性相位 或：如果()满足下式：()=0-,0是起始相位，为常数。也称 具有线性相位特性 严格地说，此时()不具有线性相位，但由于满足群时延是一个常数，即所以也称()=0-为近似线性相位,下面讨论h(n)、0要满足什么条件，可使 具有线性相位：两式相除得：,所以，如果h(n)是以 为中心作偶对称（即 h(n)=h(N-1-n)）,那么 就必须是以 为中心作奇对称。这等效地要求：0=0,=.反之，如果h(n)是以 为中心作奇对称（即 h(n)=-h(N-1-n)）,则要求 是以 为中心的偶对称。这等效地要求：,所以，满足第一类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称，即h(n)=h(N-n-1)。满足第二类线性相位的条件是：h(n)是实序列且对(N-1)/2奇对称，即h(n)=-h(N-n-1),相对于N为奇数和偶数，线性相位FIR DF的h(n)具有四种形式，它们对应了不同类型的滤波器。1、偶对称h(n)=h(N-1-n)N为奇数 2、偶对称h(n)=h(N-1-n)N为偶数3、奇对称h(n)=-h(N-1-n)N为奇数 4、奇对称h(n)=-h(N-1-n)N为偶数,二、线性相位FIR DF幅度特性Hg()的特点1、h(n)=h(N-n-1),N=奇数设 N=2M+1,则h(n)的对称中心为M=,除h(M)外其余各项满足：h(M-n)=h(M+n),1nM,由于式中 项对=0,2皆为偶对称，因此幅度特性的特点是对=0,2是偶对称的。相位特性：显然，它是的线性函数。可以实现所有滤波特性,2)h(n)=h(N-n-1),N=偶数 设N=2M,则h(n)的对称中心为,h(n)的对称关系可写为：h(M-n)=h(M-1+n),1nM,可知：Hg()对=点呈奇对称,且在=处有一零点。对于高通和带阻不适合,3)h(n)=-h(N-n-1),N=奇数 设N=2M+1，则h(n)的对称中心为M=,h(M)=0，除h(M)外其余各项满足：h(M-n)=-h(M+n),1nM可证明：Hg()对=0和=均呈奇对称。只能实现带通滤波器,4)h(n)=-h(N-n-1),N=偶数N=2M，对称中心为M-1/2。h(M-1+n)=-h(M-n)关于=0、=2奇对称，=偶对称，不能实现低通、带阻,四种波形的幅度特性和相位特性如表所示：,三、系统函数H（Z）的零极点分布 令m=N-n-1,则有可看出，H（Z-1）的零点也是H（Z）的零点，反之亦然。,一般情况下,如果 是H（Z）的零点,则:也是H（Z）的零点.,设 H（Z）的一个零点为：、取不同的值,处于不同的位置1、,处于单位圆内2、,在实轴上3、,在单位圆上4、，在单位圆和实轴的交点上。,在第一种情况下，H（Z-1）的零点 也是H（Z）的零点，它与 是以单位圆为镜象对称的。因为h(n)一般都是实数，所以H（Z）的复数零点为共轭成对的。即 也是H（Z）的零点。所以如果H（Z）有一个零点，那么、都是H（Z）的零点，它们构成一个四阶系统，其系统函数H（Z）为：,在第二种情况下：,它无共轭零点存在，但有镜象零点 所以它们可构成一个二阶系统：在第三种情况下：,它无镜象零点，但有共轭零点，,它们可构成一个二阶系统：,在第四种情况下：既无镜象零点，又无共轭零点是一个简单的一阶系统这样，一个具有线性相位的FIR DF，其系统函数可表达为上述各式的级联。即：,8.2 利用窗函数法设计FIR DF 一、窗函数法设计FIR DF 设所希望设计的滤波器传输函数为Hd(ej),则其DTFT变换对为:：是与其对应的单位脉冲响应。由 可求出:,一般Hd(ej)是矩形频率特性,所以hd(n)是非因果的，且hd(n)从，物理上无法实现。但由此可得到一个逼近Hd(ej)的方法。即：将hd(n)截短为有限项,设为N项，则：,窗函数序列的形状及长度的选择很关键,为窗函数,其波形如图所示：,中心点在 的偶对称无限长非因果序列,信号特点：,为了构造一个长度为N的线性相位滤波器，将hd(n)截短为N长，即：,hd(n)必须是对称的，取对称中心,如图所示：,图8.2.1 理想低通的单位脉冲响应及矩形窗,由h(n)求得H（Z）：对应的频响特性为：上述设计方法由于所设计的线性相位FIR DF的h(n)是由Hd(ej)的傅立叶级数的系数hd(n)，截短后得到的，所以称为傅立叶级数法。同时又可将hd(n)截短的过程视为hd(n)乘以矩形窗口序列，又称为矩形窗口法。,加窗截断的影响：,取矩形窗函数：,则：,其中：,的波形如图所示,信号特点：,有主瓣和旁瓣，主瓣宽度为,正是这些主瓣和旁瓣的影响产生了吉伯斯现象。该现象引起通带内和阻带内的波动性，尤其使阻带的衰减小，从而满足不了技术上的要求。,也表示为：,则：,卷积过程如图所示,2、，一半重叠，,1、，H（0）值可近似看作 的全部积分面积,图8.2.2 矩形窗对理想低通幅度特性的影响,3、，最大旁瓣在外，卷积结果出现最大肩峰值,4、,卷积结果达到最负值，出现负的肩峰,小结：加窗处理后对原理想低通Hd(ej)的影响：(1)在理想特性不连续点=c附近形成过渡带。过渡带的宽度近似为4/N（WR()主瓣宽度）(2)通带内增加了波动，最大的峰值在 处。阻带内产生了余振，最大的负峰在 处。以上两点就是对hd(n)用矩形窗截断后在频域的反映。称为吉伯斯现象。.,增加截取长度N，则矩形窗幅度谱：,1、增加N，主瓣宽度变窄,可得出：,2、当x增大（N增大）时，主瓣幅度增大但同时旁瓣幅度也增加，保持主瓣和旁瓣幅度相对值不变,总结：N增大，的幅度波动并没有改善，如矩形窗时最大肩峰值比H（0）高8.95%，最大负峰比0值小8.95%,为了减小吉伯斯现象，可用一些旁瓣较小的窗口来代替矩形窗口,几种常用的窗函数设 h(n)=hd(n)w(n)w(n)表示窗函数。1.矩形窗(Rectangle Window)WR(n)=RN(n)其频率响应为其主瓣宽度为4/N，第一副瓣比主瓣低13dB,2.三角形窗(Bartlett Window)其频率响应为 其主瓣宽度为8/N，第一副瓣比主瓣低26dB,3.汉宁(Hanning)窗升余弦窗 当N1时，N-1N，其主瓣宽度为8/N,能量更集中在主瓣中。如图所示,图8.2.3 汉宁窗的幅度特性,4.哈明(Hamming)窗改进的升余弦窗其频域函数WHm(e j)为其幅度函数WHm()为当N1时，可近似表示为这种改进的升余弦窗能量更加集中在主瓣中，主瓣的能量约占99.96%,图7.2.4 常用的窗函数,表7.2.2 六种窗函数的基本参数,用窗函数法设计FIR DF 的步骤：,1、根据技术要求确定待求滤波器的单位脉冲响应,如果 复杂，可对 从 采样M个点，采样值为，则：,根据频率采样定理：,当M足够大时 是 的有效逼近,2、根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗口宽度N，设要求的过渡带宽为，则,3、计算滤波器的单位脉冲响应,4、验证技术指标是否满足要求,如矩形窗，A=4,例：设计一个N=13的线性相位FIR DF，使其幅频特性接近理想低通滤波器，理想低通滤波器的截频fc=100Hz，取样频率fs=1000Hz解：理想低通DF：,8.3 用频率采样法设计FIR滤波器,一、基本设计思想,n=0,1,2,N-1,求其Z变换,得系统函数：,同样，由频域内插公式利用这N个频域采样值Hd(k)也可求得FIR滤波器的系统函数H(z),基本思想：使所设计的FIR数字滤波器的频率特性在某些离散频率点上的值准确地等于所需滤波器在这些频率点处的值，在其它频率处的特性则有较好的逼近。,内插公式,二、线性相位约束条件,为了设计线性相位的FIR滤波器，采样值 要满足一定的约束条件。,n=0,1,2,N-1,已知：,对其取共轭，得：,因为 为实函数，则有,可求得：,将 表示为：,的取值应使 具有线性相位,可求得：,式中：,注意：当N为偶数时，由于，故：,所以，用频率采样法设计高通和带阻滤波器时，N不能取偶数。,三、逼近误差分析,由上述方法求得 现分析 与 的逼近程度,已知频域采样内插公式:,式中,()是内插函数,由()的幅度谱可看出,在各频率采样点=2k/N，k=0,1,2,N-1上,(-2k/N)=1,内插公式表明：在各采样点上，,逼近误差为零，频率响应 严格地与理想频响的采样值Hd(k)相等；在采样点之间，频率响应由各采样点的内插函数延伸迭加而形成，因而有一定的逼近误差，误差大小与理想频率响应的曲线形状有关，理想特性越平滑，则内插值越接近理想值，逼近误差越小，如图a所示；反之，如果采样点之间的理想频率特性变化越陡，则内插值与理想值的误差就越大，因而在理想频率特性的不连续点附近，就会产生肩峰和起伏。使阻带衰减减小。如图b所示.N增大，则采样点变密，逼近误差减小。,图a,图b,改进措施:,在频率响应间断点附近区间内插入一个或几个过渡采样点，使不连续点变成缓慢过渡，即人为地加一过渡带。如下图所示。这样就减小了频带边缘的突变，减小了通带和阻带的波动，因而增大了阻带最小衰减。,（a）一点过渡带;(b)二点过渡带;(c)三点过渡带,例:用频率采样法设计一个带通数字滤波器,其通带频率是500Hz700Hz,采样频率为fs=3300Hz,使用阶次N=33,解:对应的数字通带频率:,得:,H(ej)的幅频特性及衰减特性分别如图所示:显然其通带及阻带内都有较大的波纹.,增加两个过渡点,令,的幅值为0.5,重新求出H(ej)的幅频特性及衰减特性如图中红线所示:,显然,特性得到了较大的改善,8.4 FIR滤波器和IIR滤波器的比较,从性能上说，IIR滤波器可以用较少的阶数获得很高的选择特性，这样一来，所用存储单元少，运算次数少，较为经济而且效率高。但是这个高效率的代价是以相位的非线性得来的。FIR滤波器可以得到严格的线性相位。但是，如果需要获得一定的选择性，则阶数比较高,成本也高，信号延时较大。如果按相同的选择性和相同的相位线性要求的话，那么，IIR滤波器就必须加全通网络来进行相位校正，因此同样要大大增加滤波器的节数和复杂性。所以如果相位要求严格一点，那么采用FIR滤波器不仅在性能上而且在经济上都将优于IIR。,从结构上看，IIR必须采用递归型结构，极点位置必须在单位圆内;否则,系统将不稳定。此外，在这种结构中，由于运算过程中对序列的四舍五入处理，有时会引起微弱的寄生振荡。相反，FIR滤波器主要采用非递归结构，不论在理论上还是在实际的有限精度运算中都不存在稳定性问题，运算误差也较小。此外，FIR滤波器可以采用快速傅里叶变换算法，在相同阶数的条件下，运算速度可以快得多。,从设计来看，IIR滤波器可以借助模拟滤波器的成果，一般都有有效的封闭函数的设计公式可供准确的计算。又有许多数据和表格可查，设计计算的工作量比较小，对计算工具的要求不高。FIR滤波器设计则一般没有封闭函数的设计公式。窗口法虽然仅仅对窗口函数可以给出计算公式，但计算通阻带衰减等仍无显式表达式。一般，FIR滤波器设计只有计算程序可循，因此对计算工具要求较高。,此外，IIR滤波器主要是用于设计具有片段常数特性的滤波器，如低、高、带通及带阻等，往往脱离不了模拟滤波器的格局。而FIR滤波器则要灵活的多，尤其是频率采样设计法更容易适应各种幅度特性和相位特性的要求，可以设计出理想的正交变换、理想微分、线性调频等各种重要网络。因而有更大适应性。,从以上简单比较我们可以看到IIR滤波器与FIR滤波器各有所长，在实际应用时要从多方面考虑来加以选择。从使用要求来看，如对相位要求不敏感的语言通讯等，选用IIR较为合适。而对图像信号处理、数据传输等以波形携带信息的系统，一般对线性相位要求较高，这时采用FIR滤波器较好。当然,在实际设计中,还应综合考虑经济上的要求以及计算工具的条件等多方面的因素。,FIR DF的窗函数设计法,函数Fir1()采用经典窗函数法设计FIR DF,调用格式：,n：FIR DF的阶数,Wn：滤波器的截止频率，01,窗函数，缺省时，自动取哈明窗,b:为滤波器的系数向量,FIR具有下列形式:,Matlab应用,为窗函数,列向量,其长度为N+1.,Matlab提供的窗函数有:boxcar,hanning,hamming,bartlett,blackmann,kaiser等,例：用窗函数法设计一个线性相位FIR LP DF，并满足性能指标：,阻带衰减不小于40dB，所以取汉宁窗,wp=0.5*pi;ws=0.66*pi;ap=3;as=40;wdelta=ws-wp;N=ceil(8*pi/wdelta)Wn=(0.55+0.66)*pi/2;b=fir1(N,wn/pi,hanning(N+1)Freqz(b,1,512),