数字信号处理窗函数法设计fir滤波器,窗函数选取原则.pptVIP

小结 37 通过上面两点讨论可以看出： 【1】根据阻带最小衰减选取窗函数的类型（尽量减少时域截取时的突变）。 【2】根据过渡带宽度确定所选窗函数的长度N。 【1】最小阻带衰减只由窗形决定，不受N的影响。 【2】过渡带的宽度既和窗形有关，又和N有关,且随窗长N的增加而减小。 因此实际中对于窗函数的选取应该遵循以下几个原则： 举例说明窗函数 选取原则 3 38 例题 试设计一个低通滤波器，其中 ， ， ,若要求 ，则应该选取什么窗函数，选取的窗函数的长度变化对滤波性能有什么影响。 39 例题 40 解： 计算可得： ， 求得 要求 ，因此在这里选择用汉明窗设计FIR低通滤波器。 例题 41 ws=0.2*pi; wp=0.4*pi; wd=0.5*(wp+ws); B=wp-ws; N=ceil(6.6*pi/B); n=0:N-1; window=hamming(N); [h1,w]=freqz(window,1); figure; stem(window); xlabel(n); title(汉明窗函数); figure; plot(w/pi,20*log10(abs(h1)/abs(h1(1)))); grid; xlabel(w/pi); ylabel(幅度(dB)); title(汉明窗的频谱); hn=fir1(N-1,wd/pi,hamming(N)); [h2,w]=freqz(hn,1); figure; stem(n,hn); xlabel(n); ylabel(h(n)); title(汉明窗的单位脉冲响应); grid; figure; plot(w/pi,20*log10(abs(h2)/abs(h2(1)))); grid; xlabel(w/pi); ylabel(幅度(dB)); title(汉明窗低通滤波器的频谱); 例题 42 例题 43 例题 44 过渡带宽度为 ，符合设计指标要求。 通过matlab可得： 总结 4 45 总结 窗谱主瓣尽可能的窄，以获得较陡的过渡带。（措施：加大窗长即增加N）。? 1 以上讨论可看出，为了满足工程上需要,可以通过改变窗函数的形状来改善滤渡器的幅频特性，窗函数的选择原则是： 46 2 具有较低的旁瓣幅度，尤其是第一旁瓣的幅度,旁瓣的幅度下降的速率要快．有利于增加阻带衰减。?（措施：选择合适的窗函数） ?。? 但在实际设计过程中 以上几点难以同时满足．实际中选用的窗函数往往是它们的折中，要根据设计指标的要求选取窗函数的形状和长度，只要设计出的FIR滤波器满足设计指标要求即可。 47 谢谢 Thanks 请加入标题 Logo 5 stem(window3); xlabel(n); title(汉宁窗函数); figure; plot(w/pi,20*log10(abs(h3)/abs(h3(1)))); grid; xlabel(w/pi); ylabel(幅度(dB)); title(汉宁窗的频谱); window4=hamming(N); [h4,w]=freqz(window4,1); figure; stem(window4); xlabel(n); title(汉明窗函数); figure; plot(w/pi,20*log10(abs(h4)/abs(h4(1)))); grid; xlabel(w/pi); ylabel(幅度(dB)); title(汉明窗的频谱); window5=blackman(N); [h5,w]=freqz(window5,1); figure; stem(window5); xlabel(n); title(布莱克曼窗函数); figure; plot(w/pi,20*log10(abs(h5)/abs(h5(1)))); grid; xlabel(w/pi); ylabel(幅度(dB)); title(布莱克曼窗的频谱); beta=8; window6=kaiser(N,beta); [h6,w]=freqz(wind