第6章FIR滤波器设计1新剖析.ppt

将 值代入下式，求得滤波器的幅频特性为 画出加汉宁窗后滤波器的幅频特性曲线，如图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,N)-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=2) ;hold on plot(w,b,--) （3）带通滤波器的设计 已知理想带通滤波器的频率响应为 其中， ，求得无限单位脉冲响应为 由上式可以看出，当 时，为低通滤波器，当 时，为高通滤波器。根据 ， ， ，并 选用偶对称形式，即可求得 共13个 的值。 当采用矩形窗时 当采用汉宁窗时 将计算出的 值代入式（7.19），可求出加矩形窗和 汉宁窗后的幅频特性，见图7.8。 （4）带阻滤波器的设计 已知理想带阻滤波器的频率响应为 其中， ， 求得无限单位脉冲响应为 根据 ， ， ，并选用偶对称形式，即可求得 共13个 的值。 当采用矩形窗时 当采用汉宁窗时 将计算出的 值代入式（7.19），可求出加矩形窗和汉宁 窗后的幅频特性，见图7.9 MATLAB程序代码如下 N=25;wc1=1;wc2=2; n=[0:1:N-1]; hd=ideal_lp(wc1,N)+ideal_lp(pi,N)-ideal_lp(wc2,N); %设计理想带阻滤波器 w_bl=(hann(N)); %选择汉宁窗 %w_bl=(rectwin(N)); %选矩形窗 h=hd.*w_bl;