实验三FIR数字滤波器设计.pptxVIP

实验三FIR数字滤波器设计第1页/共42页第2页/共42页实验三 FIR数字滤波器设计第3页/共42页实验目的了解FIR数字滤波器的原理。熟悉FIR数字滤波器从设计到实现的具体过程及其应用。 第4页/共42页实验原理一个N阶FIR滤波器只使用当前的和过去N-1个输入样点来获得当前的输出： 其中h(k)为滤波器的参数序列，也即时域冲击响应，H(m)称为滤波器的频率响应。FIR滤波器设计即寻求一组参数h(k)，使其H(m)具有所需要的频率响应特性。第5页/共42页实验原理FIR数字滤波器的特点(与IIR数字滤波器比较)：很容易获得严格的线性相位，避免被处理的信号产生相位失真；可得到多带幅频特性；极点全部在原点(永远稳定)，无稳定性问题；任何一个非因果的有限长序列，总可以通过一定的延时，转变为因果序列， 所以因果性总是满足；无反馈运算，运算误差小。第6页/共42页实验原理FIR数字滤波器的缺点：因为无极点，要获得好的过渡带特性，需以较高的阶数为代价；无法利用模拟滤波器的设计结果，一般无解析设计 公式，要借助计算机辅助设计程序完成。第7页/共42页一个简单的FIR低通滤波器(a). 一个5点平滑滤波器的参数序列h(k)(b). h(k)的归一化的离散频率幅度响应|H(m)|(c). H(m)的相位响应第8页/共42页(a). 归一化频率幅度响应|H(m)|(b). H(m)的相位响应(c). 滤波器的幅度响应0Hz~fs/2Hz第9页/共42页实验原理如果希望得到的理想滤波器的频响为 ，那么 FIR 滤波器的设计就在于寻找一个传递函数 去逼近 ，逼近的方法有以下三种：窗口设计法(时域逼近)频率采样法(频域逼近)最优化设计(等波纹逼近)第10页/共42页实验原理窗口设计法是从滤波器的单位脉冲响应序列，即参数序列着手，使h(k)逼近理想滤波器的单位脉冲响应序列hd(n)。而hd(n)可以通过对理想滤波器频响进行付氏反变换获得： 一般来说，理想频响是分段恒定，在边界频率处有突变，所以得到的理想单位脉冲响应hd(n)往往无限长，并且非因果。第11页/共42页实验原理为了用一个有限长的序列h(k)去近似无限长的hd(n)，最容易想到的办法就是直接截取hd(n)最主要的一段作为h(k) 。这种截取可看成是hd(n)和一个“窗函数”的乘积： h(k)=w(k) hd(n) 这里窗口函数就是矩形脉冲函数RN(k)，为了改善设计滤波器的特性，窗函数还有很多其它形式，除了截取之外，还在矩形窗内对hd(n)作一定的加权处理。第12页/共42页设计过程 a. 求解hd(n)一个截止频率为ωc的线性相位理想低通滤波器的频响为 (为滤波器的延时常数) 则：第13页/共42页设计过程 b. 加窗截取得到h(k)hd(n)是一个以 为中心的偶对称的无限长非因果序列。为了保证所得到的是线性相位FIR滤波器，截取时应以 为中心。 这里用矩形窗WR(k)截取N点作为h(k)。截取后滤波器的延时为第14页/共42页第15页/共42页设计过程 c. 计算H(ejω)由于时域相乘对应着频域卷积，所以h(k)的频响H(ejω)为理想频响Hd(ejω)和窗函数的频响W(ejω)相互卷积的结果。对频响起作用的只有幅度谱函数，且可以证明H(ω)也是Hd(ω)和WR(ω)的卷积。 矩形窗的幅度谱函数为 理想频率幅度函数为 则：第16页/共42页 矩形窗的卷积过程（P95的图4.5来说明）第17页/共42页实验原理通过分析可以了解到窗函数对理想特性的影响：改变了理想频响的边沿特性，形成过渡带，宽为4π/N，等于WR(ω)的主瓣宽度。(取决于窗长N)过渡带两旁产生肩峰和余振(带内、带外起伏)，取决于WR(ω)的旁瓣，旁瓣多，余振多；旁瓣与主瓣相对值大，肩峰强。(与N无关，取决于窗口形状)N的改变不能改变主瓣与旁瓣的比例关系，只能改变WR(ω)的绝对值大小和起伏的密度，当N增加时，幅值变大，频率轴变密，而最大肩峰永远不变(例如矩形窗时永远为8.95%的波动)，这种现象称为吉布斯(Gibbs)效应。第18页/共42页0N=15N=31-10-21-30-4000.250.50.751用矩形窗设计的?c=p/2 (fs/4)FIR滤波器的幅度响应第19页/共42页实验原理肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响。它取决于窗口形状。改变窗函数的形状，可改善滤波器的特性，窗函数有许多种，但要满足以下两点要求：窗谱主瓣宽度要窄，以获得较陡的过渡带；相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就可以减小肩峰和余振，以提高阻带衰减和通带平稳性。但实际上这两点不能兼得，一般总是通过增加主瓣