数字信号处理(曹成茂)第五章 有限长单位脉冲响应(fir)滤波器的设计方法.ppt
《数字信号处理(曹成茂)第五章 有限长单位脉冲响应(fir)滤波器的设计方法.ppt》由会员分享,可在线阅读,更多相关《数字信号处理(曹成茂)第五章 有限长单位脉冲响应(fir)滤波器的设计方法.ppt(158页珍藏版)》请在三一文库上搜索。
1、第五章 有限长单位脉冲响应(FIR)滤波器的设计方法,序言 5.1 线性相位FIR数字滤波器的特性,5.2 窗口设计法(时间窗口法),5.3 频率采样法,5.4 FIR数字滤波器的最优化设计 5.5 IIR与FIR数字滤器的比较,序言,FIR数字滤波器的差分方程描述 ,对应的系统函数,因为它是一种线性时不变系统,可用卷积和形式表示 ,比较、得:,FIR数字滤波器的特点(与IIR数字滤波器比较): 优点 :(1)很容易获得严格的线性相位,避免被处理 的信号 产生相位失真,这一特点在 宽频带信 号处理、阵 列信号处理、数据传输等系统中 非常重要; (2 )可得到多带幅频特性; (3 )极点全部在原
2、点(永远稳定),无稳定 性问题; (4 )任何一个非因果的有限长序列,总可以通过一 定的延时,转变为因果序列, 所以因果性总是 满足; (5)无反馈运算,运算误差小。,缺点:(1)因为无极点,要获得好的过渡带特性,需以较 高的阶数为代价; (2)无法利用模拟滤波器的设计结果,一般无解 析设计 公式,要借助计算机辅助设计程序完成。,5.1 线性相位FIR数字滤波器的特性,5.1.1 线性相位的条件 线性相位意味着一个系统的相频特性是频率的线性函数,即,式中为常数,此时通过这一系统的各频率分量的时延为一相同的常数,系统的群时延为,线性相位FIR滤波器的DTFT为,式中 H()是正或负的实函数。等式
3、中间和等式右边的实部与虚部应当各自相等,同样实部与虚部的比值应当相等:,将上式两边交叉相乘,再将等式右边各项移到左边,应用三角函数的恒等关系,满足上式的条件是,另外一种情况是,除了上述的线性相位外,还有一附加的相位,即,利用类似的关系,可以得出新的解答为,偶对称,奇对称,图1 线性相位特性,分四种情况,5.1.2 线性相位FIR滤波器的幅度特性,分四种情况 1 h(n) 偶对称,N为奇数 h(n)=h(N-1-n),5.1.2 线性相位FIR滤波器的幅度特性,令 ,则,令,则,由于 偶对称,因此 对这些频率也呈偶对称。,2h(n)偶对称,N为偶数 h(n)=h(N-1-n),令 ,则,或写为:
4、,由于 奇对称,所以 对 也为奇对称,且由于 时, 处必有一零点,因此这种情况不能用于设计 时 的滤波器,如高通、带阻滤波器。,3. h(n)奇对称,N为奇数,h(n)=-h(N-1-n),令 n=m+(N-1)/2,得:,所以,由于 点呈奇对称,所以 对这些点也奇对称。 由于 时, 相当于H(z)在 处有两个零点,不能用于 的滤波器设计,故不能用作低通、高通和带阻滤波器的设计。,4.h(n)奇对称,N为偶数,令,由于 在=0,处为零,所以H()在=0, 2处为零,即H(z)在z=1上有零点,并对=0,2呈奇对称。,四种线性相位FIR滤波器,四种线性相位FIR滤波器,相位特性,辅助序列,适用,
5、幅度特性及特点,N为奇数,N为偶数,N为奇数,N为偶数,低通 高通 带通 带阻,低通 带通,带通,高通 带通,四种线性相位FIR DF特性,参考表5.1 第一种情况 ,偶、奇,四种滤波器都可设计。,第二种情况,偶、偶,可设计低、带通滤波器,不能设计 高通和带阻。 第三种情况,奇、奇,只能设计带通滤波器,其它滤波器 都不能设计。 第四种情况,奇、偶,可设计高通、带通滤波器,不能设 计低通和带阻。,例1 N=5, h (0) = h (1) = h (3) = h (4) = -1/2, h (2) = 2,求幅度函数H ()。 解 为奇数并且h(n)满足偶对称关系 a (0) = h (2) =
6、 2 a (1) = 2 h (3) = -1 a (2) = 2 h (4) = -1 H () = 2 - cos- cos2 = 2- (cos+cos2),小结:,四种FIR数字滤波器的相位特性只取决于h(n)的对称性,而与h(n)的值无关。 幅度特性取决于h(n)。 设计FIR数字滤波器时,在保证h(n)对称的条件下,只要完成幅度特性的逼近即可。 注意:当H()用H()表示时,当H()为奇对称时,其相频特性中还应加一个固定相移。,5.1.3 线性相位FIR滤波器的零点特性,由该式可看出,若z=zi是H(z)的零点,则z=z-1i也一定是H(z)的零点。由于h(n)是实数,H(z)的零
7、点还必须共轭成对,所以z=z*i 及 z=1/z*也必是零点。 所以线性相位滤波器的零点必须是互为倒数的共轭对,即成四出现,这种共轭对共有四种 可能的情况: 既不在单位园上,也不在实轴上,有四个互为倒数的两组共轭 对, zi z*i 1/zi 1/z*i 图4.2(a) 在单位圆上,但不在实轴上,因倒数就是自己的共轭,所以有一对共轭零点, zi,z*i 图4.2(b) 不在单位圆上,但在实轴上,是实数,共轭就是自己,所以有一对互为倒数的零点, zi, 1/zi 图4.2(c) 又在单位圆上,又在实轴上,共轭和倒数都合为一点,所以成单出现,只有两种可能, zi=1或zi=-1 图4.2(d),同
8、样道理,对于第三种FIR滤波器, h(n)奇对称,N为奇数,因 所以z=1,z=-1都是H(z)的单根; 对于第四种滤波器, h(n)奇对称,N为偶数,H(O)=0,所以z=1是H(z)的单根。 所以,h(n)奇对称H(0)=0 N为偶数H()=0 线性相位滤波器是FIR滤波器中最重要的一种,应用最广。实际使用时应根据需用选择其合适 类型,并在设计时遵循其约束条件。,我们从幅度响应的讨论中已经知道,对于第二种FIR滤波器(h(n)偶对称,N为偶数), , 即 是 的零点,既在单位圆,又在实轴,所以,必有单根;,5.2 窗口设计法(时域),如果希望得到的滤波器的理想频率响应为 ,那么 FIR滤波
9、器的设计就在于寻找一个系统函数,频率响应 去逼近 ,逼近方法有三种: 窗口设计法(时域逼近) 频率采样法(频域逼近) 最优化设计(等波纹逼近) 时间窗口设计法是从单位脉冲响应序列着手,使h(n)逼近理想的单位脉冲响应序列hd(n)。我们知道hd(n)可以从理想频响 通过傅里叶反变换获得,但一般来说,理想频响 是分段恒定,在边界频率处有突变点,所以,这样得到的理想单位脉冲响应hd(n)往往都是无限长序列,而且是非因果的。但FIR的h(n)是有限长的,问题是怎样用一个有限长的序列去近似无限长的hd(n)。最简单的办法是直接截取一段 hd(n) 代替 h(n) 。这种截取可以形象地想象为h(n)是通
10、过一个“窗口”所看到的一段hd(n),因此 ,h(n)也可表达为h(n)和一个“窗函数”的乘积,即 h(n)=w(n) hd(n) 在这里窗口函数就是矩形脉冲函数RN(n),当然以后我们还可看到,为了改善设计滤波器的特性,窗函数还可以有其它的形式,相当于在矩形窗内对hd(n)作一定的加权处理。,设计步骤:,1)由定义,3)卷积,插值,一.矩形窗口法,则,以一个截止频率为 c的线性相位理想低通滤波器为例,讨论FIR的设计问题。 a. 对于给定的理想低通滤波器 ,计算,:低通滤波器的延时,理想特性的hd(n)和Hd(),这是一个以 为中心的偶对称的无限长非因果序列,如果截取一段n=0N-1的hd(
11、n)作为h(n),则为保证所得到的是线性相位FIR滤波器,延时 应为h(n)长度N的一半,即,其中,b.计算,c.计算 。 设 为窗口函数的频谱: 用幅度函数和相位函数来表示,则有,其线性相位部分 则是表示延时一半长度,对频响起作用的是它的幅度函数,理想频响也可以写成幅度函数和相位函数的表示形式 Hd(ej)=Hd()e-j 其中幅度函数为 两个信号时域的乘积对应于频域卷积,所以有,如果也以幅度函数和相位函数来表示 H(ej), 则实际FIR滤波器的幅度函数H()为 正好是理想滤波器幅度函数与窗函数幅度函数的卷积。,4个特殊频率点看卷积结果:,(1)=0时, H(0)等于,在-c, c内的积分
12、面积,因一般,故H(0)近似为,在-, 内的积分面积,(2)=c时,一半重叠, H(c)=0.5 H(0);,(3) =c 2/N时,第一旁瓣(负数)在通带外,出现正肩峰;,( 4) =c +2/N 时,第一旁瓣(负数)在通带内,出现负肩峰。,矩形窗的卷积过程,窗口函数对理想特性的影响: 改变了理想频响的边沿特性,形成过渡带,宽为 , 等于WR()的主瓣宽度。(决定于窗长) 过渡带两旁产生肩峰和余振(带内、带外起伏),取决于 WR()的旁瓣,旁瓣多,余振多;旁瓣相对值大,肩峰强 ,与 N无关。(决定于窗口形状) N增加,过渡带宽减小,肩峰值不变。 因主瓣附近 其中x=N/2,所以N的改变不能改
13、变主瓣与旁瓣的比例关系,只能改变WR()的绝对值大小和起伏的密度,当N增加时,幅值变大,频率轴变密,而最大肩峰永远为8.95%,这种现象称为吉布斯(Gibbs)效应。,改变窗函数的形状,可改善滤波器的特性,窗函数有 许多种,但要满足以下两点要求:,窗谱主瓣宽度要窄,以获得较陡的过渡带; 相对于主瓣幅度,旁瓣要尽可能小,使能量尽量集中在主瓣中,这样就 可以减小肩峰和余振,以提高阻带衰减和通带平稳性。 但实际上这两点不能兼得,一般总是通过增加主瓣宽度来换取对旁瓣的抑制。,肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减,所以对滤波器的性能有很大的影响。,5.2.2 几种常用的窗函数 1. 矩
14、形窗,上面已讲过,不再细述 2. 汉宁窗(升余弦窗) 利用的移位特性,汉宁窗频谱的幅度函数W()可用矩形窗的幅度函数表示为:,三部分矩形窗频谱相加,使旁瓣互相抵消,能量集中在主瓣,旁瓣大大减小,主瓣宽度增加1倍,为 。,0,-80,-60,-44,-20,0,矩形窗,Hanning窗,dB,3. 汉明窗(改进的升余弦窗) 它是对汉宁窗的改进,在主瓣宽度(对应第一零点的宽度)相同的情况下,旁瓣进一步减小,可使99.96%的能量集中在窗谱的主瓣内。 4. 布莱克曼窗(三阶升余弦窗) 增加一个二次谐波余弦分量,可进一步降低旁瓣,但主瓣宽度进一步增加,为 。增加N可减少过渡带。 频谱的幅度函数为:,窗
15、口函数的频谱 N=51,A=20lg|W()/W(0)|,四种窗函数的比较,(a)矩形窗;(b) 汉宁窗;(c)汉明窗;(d)布莱克曼窗,几种窗函数的基本参数,5.凯塞窗 以上四种窗函数,都是以增加主瓣宽度为代价来降低旁瓣。凯塞窗则可自由选择主瓣宽度和旁瓣衰减。,I0(x)是零阶修正贝塞尔函数,参数可自由选择,决定主瓣宽度与 旁瓣衰减。越大,w(n)窗越窄,其频谱的主瓣变宽,旁瓣变小。一般取 49。,=5.44 接近汉明 =8.5 接近布莱克曼 =0 为矩形,图2 凯塞窗函数,图1 零阶修正贝塞尔函数,I0(x),x,0,1,不同 值凯塞窗特性,当MN时, hM(n)hd(n) 零阶贝塞尔函数
16、,窗口设计法的主要工作是计算hd(n)和w(n),当 较为复杂时,hd(n)不容易由反傅里叶变换求得。这时一般可用离散傅里叶变换代替连续傅里叶变换,求得近似值: 令,过渡带宽 At阻带最小衰减,在窗口函数方面,MATLAB可计算以下的窗口函数: W=hanning(N) W=hamming(N) W=blackman(N) W=kaiser(N, beta) 变量是窗函数的长度,beta就是凯塞窗的参数,它可以利用以上讨论的公式计算,是返回的窗口函数,长度为。 对于凯塞窗,MATLAB还提供了一个根据滤波器设计指标计算窗口函数的参数的函数文件,其调用形式如下: M, Wc, beta, fty
17、pe=kaiserord(f, a, dev, fs) f为通带和阻带边界频率,最高为fs/2;a为相应频带的幅度值;dev为波动值;fs缺省为2。M是滤波器阶数;Wc边界频率;beta即;ftype是滤波器类型,低通为low。,5.3 频率采样法,工程上,常给定频域上的技术指标,所以采用频域设计更直接。 一、基本思想 使所设计的FIR数字滤波器的频率特性在某些离散频率点上的值准确地等于所需滤波器在这些频率点处的值,在其它频率处的特性则有较好的逼近。,内插公式,逼近误差,由 得到了H(z) 或 。 要讨论 与 的逼近程度 ,以及 与H(k)的关系? 由,令,,则,单位圆上的频响为:,这是一个内
18、插公式。,式中 为内插函数 令 则,内插公式表明: 在每个采样点上, 逼近误差为零,频响 严格地与理想频响的采样值 H(k)相等;,在采样点之间,频响由各采样点的内插函数延伸迭加而形成,因而有一定的逼近误差,误差大小与理想频率响应的曲线形状有关,理想特性平滑,则误差小;反之,误差大。在理想频率响应的不连续点附近, 会产生肩峰和波纹。 N增大,则采样点变密,逼近误差减小。,二.设计方法,1)确定,并得到,3)计算,2)计算,三、 约束条件,为了设计线性相位的FIR滤波器,采样值 H(k)要满足一定的约束条件。 前已指出,具有线性相位的FIR滤波器,其单位脉冲响应h(n)是实序列,且满足 ,由此得
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号处理曹成茂第五章 有限长单位脉冲响应fir滤波器的设计方法 数字信号 处理 曹成茂 第五 有限 单位 脉冲响应 fir 滤波器 设计 方法
链接地址:https://www.31doc.com/p-3057422.html