1、电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系2022-12-827.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法7.1 引言引言7.4 FIR数字滤波器的计算机辅助设数字滤波器的计算机辅助设7.5 FIR数字滤波器的特点数字滤波器的特点计方法计方法电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系7.1 引言引言2022-12-83l IIR数字滤波器数字滤波器它们的相频特性往往是它们的相频特性往往是非线性非线性的,的,这就会使信号产生失真。这就会使信号产生失真。l IIR数字滤波器数字滤
2、波器采用采用递归型结构递归型结构,系统存在极点,系统存在极点,因此设计系统函数时,必须把所有的极点置于因此设计系统函数时,必须把所有的极点置于单位圆内,否则系统不稳定。单位圆内,否则系统不稳定。电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性一、定义及特点一、定义及特点定义:定义:冲激响应有有限个非零值,系统函数可表示为冲激响应有有限个非零值,系统函数可表示为10()()NnnH Zh n Z相比相比IIRIIR滤波器的优点:滤波器的优点:1 1、可以得到线性相位频率特性、可以得到线性相位频率特性,保证通带内的信号不发
3、保证通带内的信号不发生畸变生畸变。2 2、没有递归运算,因此不论在理论上或实际应用中,均没有递归运算,因此不论在理论上或实际应用中,均不会因有限字长效应所带来的运算误差使系统不稳定。不会因有限字长效应所带来的运算误差使系统不稳定。全零点模型全零点模型2022-12-84电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系二二、FIR滤波器滤波器的线性相位频率特性的线性相位频率特性1、线性相位频率特性对信号的影响、线性相位频率特性对信号的影响()1|jk jcH ee 设一滤波器具有线性相位特性,表示为:设一滤波器具有线性相位特性,表示为:设输入信号为设输入信号为x(n),其频谱为:,
4、其频谱为:)(jeX)()()(jjjeHeXeY 则输出信号的频谱如下:则输出信号的频谱如下:在滤波器的通带之内:在滤波器的通带之内:kjjjeeXeY1)()(进行傅立叶反变换,根据性质,可得通带内输出信号为:进行傅立叶反变换,根据性质,可得通带内输出信号为:()()y nx nk7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性2022-12-85电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系结论:结论:l线性相位滤波器线性相位滤波器不会改变输入信号的形状不会改变输入信号的形状,而只是将,而只是将信号在时域上延迟了信号在时域上延迟了k个采样点。个采样点。l如果
5、滤波器的相频特性不是线性的,则输出信号就不如果滤波器的相频特性不是线性的,则输出信号就不会同输入信号一样,甚至会产生严重会同输入信号一样,甚至会产生严重畸变畸变。这就是这就是希望得到线性相位的原因。7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性二二、FIR滤波器滤波器的线性相位频率特性的线性相位频率特性1、线性相位频率特性对信号的影响、线性相位频率特性对信号的影响2022-12-86电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性1、线性相位频率特性对信号的影响、线性相位频率特性对信号的影响2022-1
6、2-87电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系线性相位的因果线性相位的因果FIR系统应具有系统应具有中心对称特性中心对称特性,即,即:2、FIR系统系统为为线性相位线性相位系统需系统需满足的条件满足的条件()(1)h nh Nn N可以是偶数也可以是奇数,对称可以是偶对称也可可以是偶数也可以是奇数,对称可以是偶对称也可以是奇对称,因此中心对称特性分为四种情况:以是奇对称,因此中心对称特性分为四种情况:N为偶数,偶对称为偶数,偶对称 N为偶数,奇对称为偶数,奇对称 N为奇数,偶对称为奇数,偶对称 N为奇数,奇对称为奇数,奇对称 可以证明在每种对称情况下,可以证明在每种对称
7、情况下,FIR滤波器都满足线性滤波器都满足线性相位特性。相位特性。7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性2022-12-88电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系n为偶数,偶对称为偶数,偶对称()(1)h nh Nn10)()(NnnZnhZH()(1)h nh Nn由于:由于:120)1()()(NnnNnZZnhZH12(1)011112()()()2220()()()NjjnjNnnNNNNjjnjnnH eh neeh n eee7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性2022-12-89电子信息工程学院电子信息工程
8、学院 生物医学工程系生物医学工程系112()201()()2cos()2NNjjnNH eeh nn分析:分析:1()2N 1201()()2cos()2NnNHh nn7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性n为偶数,偶对称为偶数,偶对称相位响应,为线性相位响应,为线性幅度响应,呈现低通特性幅度响应,呈现低通特性2022-12-810电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系相应的网络结构流图如下:相应的网络结构流图如下:7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性2022-12-811n为偶数,偶对称为偶数,偶对称电子信息工程学院
9、电子信息工程学院 生物医学工程系生物医学工程系7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性2022-12-812电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性2022-12-813电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系3、FIR线性相位线性相位系统系统的零点分布的零点分布特征特征()(1)h nh Nn()()(1)nnnnH zh n zh Nnz变量代换变量代换1m N n(1)(1)()()()mNNmmmH zh m z zzh m z(1)1()()NH zzH z
10、 表明:若表明:若z=z0是系统的零点,则是系统的零点,则z=1/z0也是系统的零点也是系统的零点7.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性2022-12-814电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系*()()h nh n根据根据ZT的性质:的性质:*()()()()ZT h nZT h nHzH z表明:若表明:若z=z0是系统的零点,则是系统的零点,则 z=z0*也是系统的零点也是系统的零点结论:结论:若若z=z0是系统的零点,则是系统的零点,则01/zz0zz*0zz*01/zz四个极点同时存在。四个极点同时存在。7.2 FIR数字滤波器的线性
11、相位特性数字滤波器的线性相位特性2022-12-8153、FIR线性相位线性相位系统系统的零点分布的零点分布特征特征电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系2022-12-8167.2 FIR数字滤波器的线性相位特性数字滤波器的线性相位特性7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法7.1 引言引言7.4 FIR数字滤波器的计算机辅助设数字滤波器的计算机辅助设7.5 FIR数字滤波器的特点数字滤波器的特点计方法计方法电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系 设计设计FIR滤波器最直接的方法,就是寻求系统单位取样滤波器最直接的方法,就
12、是寻求系统单位取样响应响应h(n),逼近理想滤波器的单位取样响应,逼近理想滤波器的单位取样响应hd(n)。1、加窗法加窗法设计设计FIR滤波器滤波器的思路的思路理想滤波器是非因果、非稳定的。理想滤波器是非因果、非稳定的。)(nhn7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-817电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系 对其加以改造:先进行移位,然后再截取有限的长对其加以改造:先进行移位,然后再截取有限的长度(加窗)从而得到一个因果稳定序列:度(加窗)从而得到一个因果稳定序列:nNnmnhnhd其它 010 )()(可以视为一个无限长序列
13、可以视为一个无限长序列hd(n-m)和一个有限长窗序列和一个有限长窗序列w(n)的乘积,即:的乘积,即:10 )()()(Nnnwmnhnhd思考思考:为什么为什么m的取值要保证的取值要保证h(n)中心对称?中心对称?7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-8181、加窗法加窗法设计设计FIR滤波器滤波器的思路的思路这是一种时域的逼近这是一种时域的逼近电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系2、加窗后频率响应的变化、加窗后频率响应的变化 加窗的实质是两者在时域相乘,反映在频域就是频率加窗的实质是两者在时域相乘,反映在频域就是频率响应卷
14、积的关系。响应卷积的关系。7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-819电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系 从图中可以看出,理想低通加矩形窗处理产生过渡带和从图中可以看出,理想低通加矩形窗处理产生过渡带和通带起伏,称为通带起伏,称为吉布斯效应。吉布斯效应。7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-820电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系吉布斯效应产生原因:吉布斯效应产生原因:1.1.理想幅频特性的陡直的边沿被加宽,形成一个理想幅频特性的陡直的边沿被加宽,形成一个过渡
15、过渡带带,过渡带的过渡带的宽度取决于窗函数频宽度取决于窗函数频响的主瓣宽度响的主瓣宽度。2.2.在过渡带两侧附近产生起伏的在过渡带两侧附近产生起伏的肩峰和波纹肩峰和波纹,它是由窗它是由窗函数频响的函数频响的旁瓣引起的旁瓣引起的,旁瓣相,旁瓣相对值愈大起伏就愈对值愈大起伏就愈强。强。3.3.增加截取长度增加截取长度N N,将缩小窗函数的主瓣宽度,将缩小窗函数的主瓣宽度,但却不,但却不能减小旁瓣相对值。能减小旁瓣相对值。旁瓣与主瓣的相对值主要取决旁瓣与主瓣的相对值主要取决于窗函数的形状于窗函数的形状。因此,增加截取长度。因此,增加截取长度N N,只能减小,只能减小过渡带宽度,而不能改善滤波器通带内
16、的平稳性和过渡带宽度,而不能改善滤波器通带内的平稳性和阻带中的衰减。阻带中的衰减。7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-821电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系为得到良好的滤波器性能,尽可能要求窗函数:为得到良好的滤波器性能,尽可能要求窗函数:1.主瓣宽度主瓣宽度窄,以获得较陡的过渡带。窄,以获得较陡的过渡带。2.旁瓣相对值旁瓣相对值尽可能小,以改善通带平稳度和增大阻带尽可能小,以改善通带平稳度和增大阻带中衰减。中衰减。7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-8223、各种窗函数的特性各种
17、窗函数的特性以上是窗的两个典型特征!加窗设计实质就是根据过渡以上是窗的两个典型特征!加窗设计实质就是根据过渡带及阻带衰减选择合适的窗!带及阻带衰减选择合适的窗!电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系3、各种窗函数的特性各种窗函数的特性矩形窗矩形窗 nNnnw其它 010 1)(1、定义公式、定义公式2、频率响应、频率响应mjmjjeNeWeW)2/sin()2/sin()()(|)2/sin(|)2/sin(|)2/sin()2/sin(|)(|NNeWj7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-823电子信息工程学院电子信息工程学院
18、 生物医学工程系生物医学工程系3 3、主瓣的宽度、主瓣的宽度4N4 4、第一副瓣衰减、第一副瓣衰减00sin(/2)|sin(/2)|()|limsin(/2)|sin(/2)|jNNW eN3/N 11|()|sin(/2)|/2jW edBN1332log20/2log207.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-824矩形窗矩形窗 电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系10)12cos(1 21)(NnnjjeNneW101212)21(21NnnjNnjNnjeee10)12(10)12(10414121NnnNjNnnNjN
19、nnjeee)2sin()2sin(41)2sin()2sin(41)2sin()2sin(21222111212121NeNeNeNjNjNj12,1221NN其中:)10()12cos(1 21)(NnNnnw汉宁汉宁(Hanning)窗窗7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-825电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系|)12(41)12(41)(21|)(|NWNWWeWRRRj7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-826汉宁汉宁(Hanning)窗窗电子信息工程学院电子信息工程学院
20、 生物医学工程系生物医学工程系海明海明(Hamming)窗窗 在在+=1的条件下的条件下,如何选取如何选取、使得幅频响应的第二副使得幅频响应的第二副瓣最大值衰减到瓣最大值衰减到0(实际上做不到,只能接近为(实际上做不到,只能接近为0),这就),这就是是海明窗海明窗,结果为,结果为 前述的汉宁窗可以写成前述的汉宁窗可以写成)12cos()12cos(1 21)(NnNnnw)12cos(46.054.0)(Nnnw幅频响应为幅频响应为|)12(23.0)12(23.0)(54.0|)(|NWNWWeWRRRj7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-827电子
21、信息工程学院电子信息工程学院 生物医学工程系生物医学工程系勃莱克曼勃莱克曼(Blackman)窗窗 前两个窗是消去第一副瓣的影响前两个窗是消去第一副瓣的影响,此窗是把第一、此窗是把第一、第二副瓣的影响都消去第二副瓣的影响都消去,其序列为:,其序列为:幅频响应为幅频响应为)14cos(08.0)12cos(5.042.0)(NnNnnw)(42.0|)(|RjWeW)12(25.0)12(25.0NWNWRR|)(04.0)(04.0)14()14(NjRNjReWeW7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-828电子信息工程学院电子信息工程学院 生物医学工
22、程系生物医学工程系凯塞凯塞(Kaiser)窗窗 凯塞窗在保证线性相位的条件下利用第一类修正零凯塞窗在保证线性相位的条件下利用第一类修正零阶贝塞尔函数对序列值进行调整,调整参数做了大量试阶贝塞尔函数对序列值进行调整,调整参数做了大量试验,其最后的效果达到最佳或较佳,验,其最后的效果达到最佳或较佳,凯塞窗序列为凯塞窗序列为)21()21()21()(00220NINnNInwa 其中其中I0(.)是第一类修正零阶贝塞尔函数,是第一类修正零阶贝塞尔函数,a为调整参数为调整参数,一般为一般为 8)21(4Na7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-829电子信息工
23、程学院电子信息工程学院 生物医学工程系生物医学工程系7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-830窗函数窗函数主瓣宽度主瓣宽度过渡带宽过渡带宽旁瓣峰值旁瓣峰值衰减衰减(dB)阻带衰减阻带衰减(dB)矩形矩形-131(4/N)-13-21汉宁汉宁-312(4/N)-31-44汉明汉明-412(4/N)-41-53布莱克曼布莱克曼-573(4/N)-57-744种种窗窗函数的比较函数的比较电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系4、加窗设计、加窗设计FIR数字滤波器的步骤数字滤波器的步骤1.根据技术要求根据技术要求(关键是阻带衰减)关键是
24、阻带衰减)先确定窗函数形式先确定窗函数形式(n)。并根据取样周期。并根据取样周期T,确定相应的数字频率,确定相应的数字频率TppTss2.根据过渡带宽根据过渡带宽=s-p确定加窗宽度确定加窗宽度N。/4PN3.确定确定单位取样:冲激响应位移系数单位取样:冲激响应位移系数m=(N-1)/2)()()(sin)(nwmnmnnhc常取常取 c=p。4.计算计算FIR滤波器的频率响应,校验技术指标滤波器的频率响应,校验技术指标7.3 FIR数字滤波器的加窗设计方法数字滤波器的加窗设计方法2022-12-831电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系1.很容易获得严格的线性相位,
25、避免被处理的信号产生相很容易获得严格的线性相位,避免被处理的信号产生相位失真,这一特点在宽频带信号处理、阵列信号理、数位失真,这一特点在宽频带信号处理、阵列信号理、数据传输等系统中非常重要;据传输等系统中非常重要;2.极点全部在原点(永远稳定),无稳定极点全部在原点(永远稳定),无稳定 性问题;性问题;3.任何一个非因果的有限长序列,总可以通过一定的延时任何一个非因果的有限长序列,总可以通过一定的延时,转变为因果序列,转变为因果序列,所以因果性总是满足;所以因果性总是满足;4.无反馈运算,运算误差小。无反馈运算,运算误差小。7.5 FIR数字滤波器的特点数字滤波器的特点优点:优点:2022-1
26、2-832电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系1.因为无极点,要获得好的过渡带特性,需以较高的因为无极点,要获得好的过渡带特性,需以较高的阶数为代价;阶数为代价;2.无法利用模拟滤波器的设计结果,一般无解析设计无法利用模拟滤波器的设计结果,一般无解析设计 公式,要借助计算机辅助设计程序完成。公式,要借助计算机辅助设计程序完成。7.5 FIR数字滤波器的特点数字滤波器的特点缺点:缺点:2022-12-833电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系2022-12-834FIR数字滤波器数字滤波器IIR数字滤波器数字滤波器性能性能严格线性相位,成本高,
27、延严格线性相位,成本高,延时大时大存储单元少,运算量小,存储单元少,运算量小,经济高效经济高效结果结果幅频和相频特性幅频和相频特性幅频特性幅频特性稳定性稳定性极点全部在原点极点全部在原点永远稳定永远稳定有时不稳定有时不稳定阶数阶数高高低低结构结构非递归非递归递归递归工具工具只有计算程序可循只有计算程序可循可借助可借助AF的设计成果的设计成果算法算法可以可以FFT实现,实现,减少运算量减少运算量无快速运算方法无快速运算方法FIR和和IIR数字滤波器的比较数字滤波器的比较电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系人有了知识,就会具备各种分析能力,明辨是非的能力。所以我们要勤恳读书,广泛阅读,古人说“书中自有黄金屋。”通过阅读科技书籍,我们能丰富知识,培养逻辑思维能力;通过阅读文学作品,我们能提高文学鉴赏水平,培养文学情趣;通过阅读报刊,我们能增长见识,扩大自己的知识面。有许多书籍还能培养我们的道德情操,给我们巨大的精神力量,鼓舞我们前进。电子信息工程学院电子信息工程学院 生物医学工程系生物医学工程系