第5章 FIR滤波器设计1.ppt

中心点为a的偶对称无限长非因果序列。 为了得到一个因果的线性相位的滤波器，进行“砍头去尾” f图，幅值用H(0)归一化 思考：窗的长度对结果有什么影响？ 主瓣宽度：4*pi/N，第一个零点位置2*pi/N，第二个零点位置4*pi/N 最大旁瓣对应位置：3*pi/N （1）主瓣越宽过渡带就越宽。 （2）由于窗口函数旁瓣的作用，使幅频特性出现波动。旁瓣所包围的面积越大，通带波动越大，阻带衰减越小。 （3）增加截取长度，只能缩小窗口频谱的主瓣宽度，而不能改变旁瓣的相对值。旁瓣与主瓣的相对关系只决定于窗口函数的形状。因此，增加截取长度只能相应地减小过渡带宽度，而不能改变过冲值的大小。 根据对过渡带及阻带衰减的要求，选择窗函数原则是保证阻带衰减满足要求的情况，尽量选择主瓣窄的窗函数； 需使窗函数谱满足： 主瓣尽可能窄，以使设计出来的滤波器有较陡的过渡带； 第一副瓣面积相对主瓣面积尽可能小，即能量尽可能集中在主瓣，外泄少，使设计出来的滤波器的肩峰和余振小。 上面两个条件是相互矛盾的，实际应用中，折衷处理，兼顾各项指标 矩形窗频谱的主瓣宽度为4?/N,第一旁瓣比主瓣低13dB。 若要求线性相位，则要求hd(n)对称与w(n)均为(N-1)/2对称 这种改进的升余弦窗，能量更加集中与主瓣，第一旁瓣的峰值比主瓣小40dB。 其他的窗函数是各以一定主瓣加宽为代价，来换取某种程度的旁瓣抑制，而Kaiser窗则是全面的反映主瓣与旁瓣衰减之间的交换关系，可以在它们两者之间自由地选择它们的比重。 实际应用中，FIR低通滤波器的阻带衰减至少要达到50dB以上，因此，很少采用矩形窗和汗宁窗。 X=fft(h,128); subplot(1,2,2); plot(abs(X)); %画幅频特性 title(Magnitude Response,N=41); axis([0,60,0,1.1]); xlabel(ω in units of \pi/N ); ylabel(H(ω)); grid Matlab Program: 当采用矩形窗时 根据式（7.19）可求得滤波器的幅频特性为 利用上式可以求得 取不同值时， 相应的值，画出幅频特性曲线，如图7.6（a）中实线所示，虚线为理想低通滤波器的幅频特性。 当采用汉宁窗时 由表7.2可知 求的汉宁窗口序列为 加汉宁窗后，有限长脉冲序列为 将 值代入下式，求得滤波器的幅频特性为 画出加汉宁窗后滤波器的幅频特性曲线，如图7.6 (b)中实线所示。从图7.6的(a)和(b)两图可以看出，加汉宁窗后幅频特性的过冲和波动基本消失，但过渡带有所增宽。 MATLAB程序代码如下 N=25;wc=1; n=[0:1:N-1]; hd=ideal_lp(wc,N); %设计理想低通滤波器 w_bl=(hann(N)); %选择汉宁窗 w_bl=(rectwin(N)); %选择矩形窗 h=hd.*w_bl; %单位脉冲响应 [db,mag,pha,grd,w,H]=freqz_m(h,[1]); %幅度响应 plot(w,mag) b=(w=1) hold on plot(w,b,--) function [db,mag,pha,grd,w,H]=freqz_m(b,a) [H,w]=freqz(b,a,1000,whole); H=(H(1:1:501));w=(w(1:501)); mag=abs(H); db=(mag+eps)/max(mag); pha=angle(H); grd=grpdelay(b,a,w) （2）高通滤波器的设计 已知理想高通滤波器的频率响应为 其中， ，求得无限单位脉冲响应为 根据 ， ， ，并选用偶对称形式，即可求得13个 的值，且 当采用矩形窗时 当采用汉宁窗时 将计算出的 值代入式（7.19），可求出加矩形窗和汉宁窗 后的幅频特性，见图7.7。 图7.7 FIR高通滤波器幅频特性 (a) 加矩形窗 (b) 加汉宁窗 MATLAB程序代码如下 N=25;wc=2; n=[0:1:N-1]; hd=ideal_lp(pi