数字信号处理材料.pptVIP

当 很复杂时，便不能按下式求 ： 是 的IDTFT。 一般采用IDFT计算： 其中： 当MN时： 由经验公式求得参数β和滤波器长度N。 其中： /At 为阻带最小衰减； △?为过渡带宽。 例2 用凯塞窗设计一FIR低通滤波器，低通边界频率ωc=0.3π，阻带边界频率 ωr=0.5π，阻带衰减At不小于50dB。 解：首先求解hd(n)，根据指标要求其边界频率应为: 式中a=(N-1)/2。 △ω=ωr-ωc=0.2π N=(50-8)/(2.285×0.2π)≈30 β=0.1102×(50-8.7)=4.55 wn=kaiser(30,4.55); nn=[0:1:29]; alfa=(30-1)/2; hd=sin(0.4*pi*(nn-alfa))./(pi*(nn-alfa)); h=hd.*wn; [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1)),r); axis([0,1,-80,10]); grid; xlabel(归一化频率/\pi); ylabel(幅度/dB); h= 对于矩形窗： ?阻带最小衰减为 20lg(8.95%) = -21dB， ?通常达不到指标的要求。 ?若要加大阻带衰减，只能改变窗函数形状。 （1）主瓣宽度尽可能小，以获取较窄的过渡带； 对窗函数要求： （2）与主瓣的幅度相比，旁瓣应尽可能小，以减少滤波器的通带衰减，增大阻带衰减。 但以上两个要求很难同时满足。 研究各种典型的窗函数特性是窗函数设计法的重要内容。 下面我们将会看到，加合适的窗函数后，能使H(?)的波纹变小，但其过渡带也将变宽。 因此说，用窗函数法设计FIR数字滤波器，其通带和阻带性能的改善是以加宽过渡带为代价的。 5.2.3 几种常用的窗函数及其特点 （1）矩形窗 前面已讨论过，现归纳如下： 其频率响应为： （2） 汉宁(Hanning)窗（升余弦窗） 我们可以求出： 当N1时，N-1≈N，其频响幅度特性函数近似为： 三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大大减小，但主瓣宽度增加1倍，为8π/N。 (a)三部分频谱相加 汉宁窗 (b)相加结果 汉宁窗 （3）汉明(Hamming)窗（改进的升余弦窗） 当N1时，其频响的幅度特性函数为： 可将99.963%的能量集中在窗谱的主瓣内。 与汉宁窗相比，主瓣宽度相同，为8π/N，但旁瓣幅度更小，旁瓣峰值小于主瓣峰值的1%。 所以人们常用hamming窗而不用hanning窗设计FIR滤波器。 （4）布莱克曼（Blackman）窗 （二阶升余弦窗） 其频响的幅度特性函数为： 这样可以得到更低的旁瓣，但是主瓣宽度却不得不进一步加宽到矩形窗的三倍（12π/N）。 四种窗口的比较 (1) ①为矩形窗， ②为汉宁窗， ③为汉明窗， ④为布莱克曼窗 四种窗函数的频谱（N=51 ） 四种窗口的比较 (2) 四种窗函数在同一指标下设计的滤波器的频率特性 （N=51,ωc=0.5π） 四种窗口的比较(3) 窗函数的主瓣宽度: 决定H(?)的过渡带； 窗函数的旁瓣： 决定通带波动和阻带波动。 四种窗口的比较(4) 设计参数的选择： （1）首先要满足阻带最小衰减； （2) 在满足（1）的条件下，选择最小 过渡带宽。 （5） 凯塞窗 凯塞窗可以自由选择主瓣宽度和旁瓣衰减。 I0(x)是零阶贝塞尔函数，β越大，则ω(n)窗越窄，频谱的旁瓣越小，但主瓣变宽也相应增加。 改变β值，就可以对主瓣变宽和旁瓣衰减进行选择。 一般选择4β9，此时旁瓣幅度与主瓣幅度的比值有3.1%变为0.047%。 β=0，相当于矩形窗。 当n=(N-1)/2，凯塞窗处于最大值，即： 从n=(N-1)/2这一点向两边变化时，ω(n)逐渐减小，最边上两点 参数β越大，ω(n)变化越快。 不同β值的凯塞窗特性 窗口设计法中的主要的计算工作： 先计算 和?(n)， 然后：h(n)= ? ?(n) 5.2 窗口设计法 5.2.1 窗口设计法的基本原理 我们要设计一个滤波器，其频率响 应为 ； ?找出与之相对应的理想滤波器，其频率 响应为 ； ?方法：用 去逼近 。 具体做法就是在时域中： 其中: 即: ? 是理想滤波器的频响， 其hd(n)是无限长的，而且是非因果的。 ? 而h(n)是有限长的。 ? 因此必须将hd(n