第7章有限脉冲响应数字滤波器的设计(少)课件

第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 l7.1 线性相位线性相位FIR数字滤波器的条件和特点数字滤波器的条件和特点 l7.2 利用窗函数法设计利用窗函数法设计FIR滤波器滤波器l7.3 利用频率采样法设计利用频率采样法设计FIR滤波器滤波器l7.5 IIR和和FIR数字滤波器的比较数字滤波器的比较10/3/20227.1 线性相位线性相位FIR DF 的条件和特点的条件和特点一.引言为何要设计为何要设计FIRFIR滤波器？滤波器？(1)语音处理，图象处理等应用要求线性相位，任意幅度。FIR DF具有严格的线性相位，同时可具有任意幅度特性。(2)FIR DF单位抽样响应是有限长的，故滤波器是稳定的。只要经过延时，任何非因果有限长序列都能变成因果。(3)FIR可以用FFT算法来实现过滤信号。IIR数字滤波器数字滤波器缺点缺点：优点优点：相位非线性能借助模拟滤波器已有成果设计10/3/202210/3/202210/3/2022二.线性相位的定义长度为N的h(n),传输函数为：Hg():幅度特性,():相位特性。线性相位:传输函数表达式中的()是的线性函数第二类线性相位：第一类线性相位：群时延:是滤波器平均延迟的一个度量第一、二类线性相位的群时延相同，为第一、二类线性相位的群时延相同，为10/3/2022三.线性相位的条件第一类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称。第二类线性相位的条件是：h(n)是实序列且对(N-1)/2奇对称。四.线性相位条件证明(1)第一类线性相位的条件证明若 ，则 证明：10/3/2022所以H(z)表示为：将z=e j代入上式，得到：10/3/2022由上式可知：群时延：10/3/2022(2)第二类线性相位的条件证明证明：若 ，则 所以H(z)表示为：10/3/2022 将z=e j代入上式，得到：由上式可知：群时延：10/3/2022五.线性相位FIR滤波器的幅度特性 由于h(n)的长度N取奇数还是偶数，对Hg()的特性有影响，因此，对于两类线性相位，下面我们分四种情况讨论其幅度特性的特点：（1）h(n)=h(N-1-n)，即h(n)为偶对称偶对称，N=奇数 （2）h(n)=h(N-1-n)，即h(n)为偶对称偶对称，N=偶数 （3）h(n)=-h(N-1-n)，即h(n)为奇对称奇对称，N=奇数 （4）h(n)=-h(N-1-n)，即h(n)为奇对称奇对称，N=偶数10/3/2022（1）h(n)=h(N-1-n)，即h(n)为偶对称偶对称，N=奇数h(n)幅度函数:上式中：h(n)对(N-1)/2偶对称，余弦项也对(N-1)/2偶对称，可以以(N-1)/2为中心，把两两相等的项进行合并，可得：10/3/2022 看出：cos(n)对于w=0,2皆为偶对称，所以幅度函数Hg()也对 =0,2皆为偶对称。且H(0),H(),H(2)都可不为零(只要h(N-1)/2)不为零)。所以从从0 2 范围内，无范围内，无任何约束，可以设计成任何一种滤波器。低通、高通、带通、任何约束，可以设计成任何一种滤波器。低通、高通、带通、带阻）带阻）10/3/2022（2）h(n)=h(N-1-n)，即h(n)为偶对称偶对称，N=偶数h(n)幅度函数:上式中：h(n)对(N-1)/2偶对称，余弦项也对(N-1)/2偶对称，可以以(N-1)/2为中心，把两两相等的项进行合并，可得：10/3/2022 看出：cos(n-1/2)在=时为零，且对 =奇对称，所以幅度函数Hg()也对 =奇对称，且在 =处有一零点，使Hg()=0。用n代替m,有：(所以这种情况不能用于设计所以这种情况不能用于设计=时，时，Hg()0的的滤波器。滤波器。如：如：高通、带阻高通、带阻)10/3/202210/3/2022（3）h(n)=-h(N-1-n)，即h(n)为奇对称奇对称，N=奇数h(n)幅度函数:上式中，利用h(n)和sin(*)的奇对称性合并相同的项得：负负得正10/3/2022用n代替m,有：看出：sin(nw)对于w=0,2处皆为0，即Hg(w)在w=0,2处必为零。也即H(z)在z=1处都为零。sin(nw)对w=0,2呈奇对称形式。(不能用于 =0,时，时，Hg()0的的滤波器滤波器设计，故不能用作低通、高通和带阻滤波器的设计。)10/3/2022（4）h(n)=-h(N-1-n)，即h(n)为奇对称奇对称，N=偶数h(n)幅度函数:上式中，利用h(n)和sin(*)的奇对称性合并相同的项得：用n代替m,有：10/3/2022 看出：sin(n-1/2)对于w=0,2处皆为0，即Hg(w)在w=0,2处必为零。也即H(z)在z=1处有一个零点,且对w=0,2呈奇对称，对w=偶对称(不能用于 =0时，时，Hg()0的的滤波器滤波器设计，故不能用作低通和带阻滤波器的设计。)10/3/2022第一一种情况：h偶、N奇，Hg()关于0，2偶对称。四种滤波器都可设计(cos(n)。第二二种情况：h偶、N偶，Hg()关于奇对称。可设计低、带通滤波器，不能设计高通和带阻(cos(n-1/2)。第三三种情况：h奇、N奇，Hg()关于0，2奇对称。只能设计带通滤波器，其它滤波器都不能设计(sin(nw)。第四四种情况：h奇、N偶，Hg()关于0，2奇对称，关于偶对称。可设计高通、带通滤波器，不能设计低通和带阻(sin(n-1/2)。六.（总结）线性相位FIR滤波器的幅度特性10/3/202210/3/2022四种FIR数字滤波器的相位特性只取决于h(n)的对称性，而与h(n)的值无关。幅度特性取决于h(n)。设计FIR数字滤波器时，在保证h(n)对称的条件下，只要完成幅度特性的逼近即可。10/3/2022七.线性相位FIR滤波器零点分布特点1.若 z=zi 是H(z)的零点,则 z=zi-1 也是其零点2.h(n)为实数，则零点共轭成对线性相位滤波器的零点是互为倒数的共轭对线性相位滤波器的零点是互为倒数的共轭对得：由h(n)=22,20,33,6,6,33,-20,2210/3/2022八.线性相位FIR滤波器网络结构设为N为偶数，则滤波器的系统函数为：要想确定滤波器网络结构，首先要得到系统的差分方程或系统函数。10/3/2022N为偶数时系统函数：由前述可得：第一类线性相位情况下第二类线性相位情况下N为奇数时系统函数：N为偶数时系统函数：N为奇数时系统函数：根据系统函数可画出如下的滤波器网络结构：10/3/2022图7.1.2 第一类线性相位网络结构10/3/2022图7.1.3 第二类线性相位网络结构 与与FIR滤波器的直接结构相比，利用线性相位的对称特滤波器的直接结构相比，利用线性相位的对称特性可节省一半左右的计算量！性可节省一半左右的计算量！10/3/2022第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 l7.1 线性相位线性相位FIR数字滤波器的条件和特点数字滤波器的条件和特点 l7.2 利用窗函数法设计利用窗函数法设计FIR滤波器滤波器l7.3 利用频率采样法设计利用频率采样法设计FIR滤波器滤波器l7.5 IIR和和FIR数字滤波器的比较数字滤波器的比较10/3/20227.2 利用窗函数法设计利用窗函数法设计FIR滤波器滤波器1.由理想的频率响应,得到理想的 ；2.由 得到因果、有限长的单位抽样响应 ;3.对 加窗得到较好的频率响应。理想频率响应一、思路与方法：一、思路与方法：表达式及特点？10/3/2022则其脉冲响应：无限长 非因果 偶对称解决方法：截短，移位 保留即：隐含着使用了窗函数理想低通的传输函数为：特点：10/3/2022于是：注意：是因果的，且是线性相位的，即即事先给一线性相位为了省去每次的移位，可以令：在通带内这样：10/3/2022于是：使用了矩形窗上式的的表达式及设计 的思路可推广到高通、带阻及带通滤波器，也可推广到其它特殊类型的滤波器。实际上，给定一个 ，只要能积分得到 ，即可由截短、移位的方法得到因果的、且具有线性相位的FIR滤波器 。为了构造一个长度为N的线性相位滤波器，只有将hd(n)截取一段，并保证截取的一段对(N-1)/2对称：即：即：10/3/2022图7.2.1 理想低通的单位脉冲响应及矩形窗窗函数法设计窗函数法设计FIR滤波器图示滤波器图示10/3/2022补充几种常用的理想滤波器补充几种常用的理想滤波器高通：令：相当于用一个截止频率在 处的低通滤波器（实际上是全通滤波器）减去一个截止频率在 处的低通滤波器。10/3/2022令：相当于用一个截止频率在 处的低通滤波器减去一个截止频率在 处的低通滤波器。带通：10/3/2022二、二、FIR低通滤波器设计步骤低通滤波器设计步骤 已知滤波器的阶数已知滤波器的阶数N，和理想低通的，和理想低通的c c 利用公式利用公式 求移位值求移位值利用公式利用公式 求理想低通的单位脉冲响应求理想低通的单位脉冲响应hd(n)利用公式利用公式 截短截短hd(n)得可实现滤波器响应得可实现滤波器响应h(n)例7.2.1 用矩形窗、汉宁窗和布莱克曼窗设计FIR低通滤波器，设N=11,c=0.2rad。h(n)=0,0.0468,0.1,0.1514,0.1871,1,0.1871,0.1514,0.1,0.0468,0解：解：10/3/2022三、关于对三、关于对 截短的讨论（幅度特性）截短的讨论（幅度特性）称为用矩形窗对称为用矩形窗对hd(n)进行处理进行处理h(n)的DTFT为：为矩形窗的幅度特性。其中：其中：将理想低通滤波器的频率响应 表示为：10/3/2022则理想低通的幅度特性：则的设计低通的幅度特性：可见：设计的滤波器的幅度特性是矩形窗函数的幅度特设计的滤波器的幅度特性是矩形窗函数的幅度特性与理想低通滤波器的幅度特性的卷积性与理想低通滤波器的幅度特性的卷积（过程见下图）10/3/2022理想低通幅度特性窗函数幅度特性10/3/2022加矩形窗处理后，对理想频率响应产生了以下两点影响：加矩形窗处理后，对理想频率响应产生了以下两点影响：使理想频率特性不连续点=c 处，形成了一个过过渡带渡带，过渡带的宽度等于矩形窗的频率响应RN()的主瓣宽度=4/N 在截止频率c的两边=c 2/N 处（即过渡带 的两边），H()出现最大的肩峰值最大的肩峰值，肩峰的两侧形成起伏振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡的快慢，则取决于RN()波动的快慢。以上两点就是对hd(n)截断后，在频域的反映，称为吉布斯效应，其直接影响滤波器的性能。10/3/2022若增加截取长度N，则在主瓣附近的窗的频率响应为：随着x加大(N加大)，函数曲线波动的频率加快(因为过零点间隔变小)，主瓣幅度加高，旁瓣幅度也同样加高，主瓣与旁瓣幅度相对值不变。减少吉布斯效应影响的方法：1、过渡带宽=4/N。由此增加窗长N，能减小过渡带宽2、增加窗长N，在减小过渡带宽同时，能改善振荡的情况吗？分析：因此，减少带内波动和加大阻带的衰减只能从窗形入手。10/3/2022四、窗函数介绍四、窗函数介绍 窗函数的使用在数字信号处理中是不可避免的。数据、频谱、自相关函数等都需要截短。对窗函数提出那几方面的要求？对加窗的理解，关键是要搞清楚使用窗函数后所产生的影响：一个域相乘，在另一个域是卷积。窗函数的要求：(1)窗谱主瓣尽可能的窄，以获得较陡的过渡带。(2)尽量减小窗谱的最大旁瓣的相对幅度，也就是能量尽量集中在主瓣，这样使肩峰和波纹减小，可以增大阻带的衰减。10/3/2022（1）矩形窗（2）三角形窗（Bartlett)10/3/2022不同窗函数的波形以及幅度谱图不同窗函数的波形以及幅度谱图10/3/2022当N1时，N-1N，（3）汉宁(Hanning)窗升余弦窗10/3/202210/3/2022(7.2.11)当N1时，可近似表示为（4）哈明(Hamming)窗改进的升余弦窗10/3/202210/3/2022(7.2.13)