首页电子与通信技术信息处理技术数字信号处理 09 第6章 FIR数字滤波器的设计与MATLAB实现-2.pptx

09 第6章 FIR数字滤波器的设计与MATLAB实现-2.pptx

页数:20页文件大小:473.6 KB来源: 数字信号处理ppt课件

本课件承接第6章后半部分,先由频率响应的内插公式出发说明频率采样法的设计思想,给出线性相位滤波器的采样值约束条件与完整设计步骤,并以N=33低通设计为例分析过渡带与阻带衰减;随后介绍等波纹最佳逼近法,阐述加权误差函数、交错定理与Remez迭代算法的求解过程,最后从性能、工具与设计灵活度等角度比较FIR与IIR两类数字滤波器。

6.4频率采样法设计FIR滤波器,由于有限长序列h(n)和它自身的DFT是一一对应的,因此可以由频域的N个频域采样值H(k)来唯一确定有限长序列h(n)
待设计滤波器的频率响应,和系统函数 根据H(z)的内插公式,可以由N个频域采样值H(k)恢复出FIR滤波器的系统函数H(z)和频率响应式中
内插公式与频率采样点,从内插公式可以看出,在各个频率采样点上,滤波器的实际频率响应是严格地和理想频率响应数值相等,频率采样点之间的频率响应是由各采样点的加权内插函数的延伸叠加而形成
查看本课件的文字版(共 20 页,供检索与复制)

第1页:6.4频率采样法设计FIR滤波器

频率采样法设计FIR滤波器 由于有限长序列h(n)和它自身的DFT是一一对应的,因此可以由频域的N个频域采样值H(k)来唯一确定有限长序列h(n)。然后以作为实际滤波器的频率特性的采样值,即 设计思想:频率采样法是从频域出发,把给定的理想滤波器的频率响应加以等间隔采样式中h(n)为待设计的滤波器的单位采样响应。

第2页:待设计滤波器的频率响应

待设计滤波器的频率响应和系统函数 根据H(z)的内插公式,可以由N个频域采样值H(k)恢复出FIR滤波器的系统函数H(z)和频率响应式中。

第3页:内插公式与频率采样点

从内插公式可以看出,在各个频率采样点上,滤波器的实际频率响应是严格地和理想频率响应数值相等,频率采样点之间的频率响应是由各采样点的加权内插函数的延伸叠加而形成。理想频率响应特性变化越平缓,则内插值越接近理想值,逼近误差越小。

第4页:6.4.2线性相位滤波器的条件

频率采样法设计线性相位滤波器的条件设的频响为:对于第I类线性相位滤波器,h(n)偶对称,N为奇数幅度函数H(w)g关于偶对称,即如果采样值也用幅值和相位特性表示,则有式中:必须满足偶对称关系。

第5页:第II类线性相位滤波器条件

频率采样法设计线性相位滤波器的条件 第II类线性相位滤波器,h(n)偶对称,N为偶数,其幅度函数H(w)奇对称: 对于第III类及第IV类线性相位FIR滤波器,两类滤波器的相位特性满足:这时满足奇对称要求,即:相位特性满足:第III类及第IV类滤波器采样值H(k)的幅值分别满足奇对称和偶对称关系。

第6页:频率采样法设计步骤

  • 频率采样法设计线性相位FIR滤波器的步骤一般,频率采样法设计线性相位FIR滤波器步骤可分为以下几步: 根据设计要求选择滤波器的类型。
  • 根据线性相位的约束条件,确定和,进而得到H(k)。
  • 将代入内插公式得到所设计的滤波器的频率响应。
  • 例题:利用频率采样法,设计一个低通FIR数字滤波器,其理想频率特性是矩形的,即已知,采样点数为N=33,要求滤波器具有线性相位。

第7页:例题:N=33低通设计

例题求解解:(1)根据设计指标要求,采样点数N=33是奇数,滤波器为低通滤波器,因此选择第I种类型滤波器,满足(2)由通带截止频率,确定通带内的采样点数,即取整数k=8,应在通带内设置9个采样点。(3)根据,可得因此。

第8页:滤波器频率响应求解

(4)将H(k)代入到的内插公式,即得到所设计滤波器的频率响应。(5)滤波器幅度特性与对数幅频特性的结果图。

第9页:过渡带与阻带衰减分析

  • 由上图可以看出过渡带带宽为,最小阻带衰减为-20dB在大多数情情况下这一衰减往往是不令人满意的。
  • 为了改善频率特性,在阻带和通带交界处安排一个或几个不等于1的采样值,就可以使得它的有用频带内的纹波得以减小。
  • 比如在本例中,在处增加一个,则得到下图所示结果,这相当于加宽过渡带带宽,其宽度为,算出阻带最小衰减为。

第10页:6.5等波纹最佳逼近法

利用等波纹最佳逼近法设计FIR数字滤波器设希望设计的滤波器幅度函数为,实际设计滤波器的幅度函数为加权误差为,则:预先指定的误差加权函数,用来控制不同频段的逼近精度等纹波滤波器设计问题就是求滤波器幅度函数,使得在一组频率点(第I类线性相位)中的系数上的E(w)绝对值最小,即。

第11页:加权函数与频率组

设计一低通滤波器,频率组可以是通带和阻带内的频率,如图6.14所示.过渡带是不关心的区域,求加权误差最小时不予考虑,可以采用交错定理求这个最优化问题。

第12页:交错定理

交错定理即在上必须存在L+2个极值频率点这样且,有:设是区间内封闭子集的并集,对于正的加权函数H(w)能成为唯一使加权误差|E(w)|最大值最小的函数。其充要条件是:在上在上E(w)至少有L+2个交错值。

第13页:交错定理的方程组

由交错定理可以得到:式中交错定理,确定了最优滤波器必须有的极值频率的最少数目这些关于未知数a(0,a(1),….a(L)以及的方程可以写成下面的形式。

第14页:系数求解与误差

方程:求解上式,求出系数a(0),a(1),a(2),…,a(L)以及误差,从而由a(n)求得最佳滤波器的单位采样响应h(n)。

第15页:Remez算法步骤

Remez算法的计算步骤Remez算法:靠一次次迭代求得一组交错点即极值频率点,从而求出系数,而且在每一次求解极值频率的迭代过程中能够避免直接求解。具体步骤为:由(6-67)式计算:式中一般不可能是最优值①首先给出L+2个极值频率的初始估计值,通常在频率区间内等间隔地取L+2个初始频率点,即。

第16页:拉格朗日插值求幅度函数

利用拉格朗日(Lagrange)插值公式得到幅度函数Remez算法的计算步骤②把L+2个极值频率点代入(6-67)式,可求得,这就是第一次指定极值频率的频率的偏差值。即:(6-69)式中。

第17页:Remez算法的迭代修正

  • Remez算法的计算步骤若在某些频率点处,则说明初始猜测频率点偏离了真正的交错频率,需要进行修改,进行后面的步骤。
  • ③利用和步骤②相同的方法,把在各频率处使的局部极值点作为新的极值频率的猜测值,选定一组新的极值频率值。
  • 把代入所有频率上都有,说明恰好就是交错点组,工作可结束。
  • 得到误差函数E(w).如果④重复以上步骤,每次迭代后的都是递增的,最后收敛到也就最佳一致逼近,求得h(n)。
  • 最佳极值,此时。

第18页:例6-6 等波纹低通设计

例6-6利用等波纹逼近法设计一个线性相位FIR低通滤波器 指标要求为:通带截止频率fp=800Hz,通带衰减0.5dB。阻带截止频率fst=1000Hz,阻带最小衰减40dB,采样频率为fs=4000Hz. 滤波器幅频特性。

第19页:6.5 FIR与IIR滤波器比较

  • FIR数字滤波器与IIR数字滤波器的比较 首先,从性能上说,IIR数字滤波器可以用较低的阶数获得好的选择特性,所用的存储单元少,运算次数少,所以经济高效。
  • FIR数字滤波器可以得到严格的线性相位,而IIR数字滤波器的选择性越好,相位的非线性越严重。
  • 从结构上看,FIR数字滤波器采用非递归结构,不存在内部的反馈,因此系统稳定,不论在理论上还是实际的有限精度运算中都不存在稳定性问题。
  • 而IIR数字滤波器采用递归结构,系统极点位置必须在单位圆以内,否则,系统将不稳定。
  • FIR数字滤波器可以采用快速傅里叶变换(FFT)算法大大减少运算量,运算速度快很多,而IIR数字滤波器没有快速运算。

第20页:FIR与IIR的设计工具比较

  • 从设计工具看,IIR数字滤波器可以借助于模拟滤波器的成果,因此一般都有有效的封闭形式的设计公式可供准确计算,计算工作量比较小,对计算工具的要求不高。
  • 而FIR数字滤波器无法借助于模拟滤波器的设计方法,只有计算程序可循,对计算工具要求较高。
  • 另外,IIR滤波器虽然设计简单,但主要是用于设计具有分段常数特性的滤波器,如低通、高通、带通和带阻滤波器等,往往脱离不了模拟滤波器的局限性。
  • 而FIR滤波器则灵活得多,尤其能适应某些特殊的应用,应用于巴特沃斯、切比雪夫等逼近不可能达到预定指标的情况。
顶部目录