有限脉冲响应数字滤波器的设计 (2).pptVIP

*第62页，共121页，星期日，2025年，2月5日常用窗函数的频谱窗型变化，旁瓣减小，阻带衰减增大，主瓣宽度增加*第63页，共121页，星期日，2025年，2月5日理想低通加窗后的幅度特性N=51,ωc=0.5π*第64页，共121页，星期日，2025年，2月5日凯泽（Kaiser）窗：第一类变形零阶 贝塞尔函数β*第65页，共121页，星期日，2025年，2月5日凯塞窗参数对滤波器的性能影响*第66页，共121页，星期日，2025年，2月5日六种窗函数的基本参数阻带最小衰减只由窗形状决定过渡带宽则与窗形状和窗宽N都有关*第67页，共121页，星期日，2025年，2月5日注意事项：对于窗函数法，加大N只能减少过渡带的宽度，不能加大阻带的衰减，要加大阻带的衰减只能通过选择合适的窗函数来实现。基本原则：在保证阻带衰减的情况下，尽量选择主瓣窄的窗函数。*第68页，共121页，星期日，2025年，2月5日六、用窗函数法设计FIR滤波器的步骤给出希望设计的滤波器的频率响应函数Hd(ejω)。计算hd(n)如果Hd(ejω)不能用简单函数表示，可以用求和代替积分。根据允许的过渡带宽及阻带衰减，选定窗函数和N值，见书表7-2,7-3*第69页，共121页，星期日，2025年，2月5日5.将hd(n)与窗函数相乘得FIR数字滤波器的冲激响应h(n)。6.计算FIR数字滤波器的频率响应，并验证是否达到所要求的指标，若不能满足要求，根据情况重复4、5、6步。*第70页，共121页，星期日，2025年，2月5日例窗函数法设计FIR滤波器例7.2.1用窗函数法设计线性相位FIR低通滤波器，设N=11，ω=0.2?rad解：理想数字低通滤波器单位取样响应*第71页，共121页，星期日，2025年，2月5日要求设计的FIR数字滤波器的单位取样响应*第72页，共121页，星期日，2025年，2月5日用不同窗函数设计的FIR滤波器用矩形窗时过渡带最窄，而阻带衰减最小，布莱克曼窗过渡带最宽，但阻带衰减加大。为保证有同样的过渡带，必须加大窗口长度N*第73页，共121页，星期日，2025年，2月5日MATLAB函数window=boxcar(N+1);矩形窗b=fir1(N,Wc,ftype,Windows)b=fir1(10,0.2*pi,window);window=boxcar(11);*第74页，共121页，星期日，2025年，2月5日MATLAB实现window=boxcar(11);b=fir1(10,0.2*pi,window);[h1,w1]=freqz(b,1);plot(w1/pi,20*log10(abs(h1)));holdon;window=hanning(11);b=fir1(10,0.2*pi,window);[h2,w2]=freqz(b,1);plot(w2/pi,20*log10(abs(h2)));holdon;window=blackman(11);b=fir1(10,0.2*pi,window);[h3,w3]=freqz(b,1);plot(w3/pi,20*log10(abs(h3)));holdon;*第75页，共121页，星期日，2025年，2月5日解：1）求数字频率设计一个线性相位FIR低通滤波器， 给定抽样频率为 通带截止频率为 阻带起始频率为 阻带衰减不小于-50dB，幅度特性如图所示例7.2.2*第76页，共121页，星期日，2025年，2月5日2）求hd(n)*第77页，共121页，星期日，2025年，2月5日4）确定N值3）选择窗函数：由确定哈明窗（-53dB）*第78页，共121页，星期日，2025年，2月5日5）确定FIR滤波器的h(n)6）求，验证若不满足，则改变N或窗形状重新设计*第79页，共121页，星期日，2025年，2月5日例7.2.3线性相位FIR高通滤波器的设计其单位抽样响应：理想高通的频响：*第80页，共121页，星期日，2025年，2月5日例7.2.4线性相位FIR带通滤波器的设计其单位抽样响应：理想带通的频响：*第81页，共121页，星期日，2025年，2月5日例7.2.5线性相位FIR带阻滤波器的设计其单位抽样响应：理想带阻的频响：*第82页，共121页，星期日，2