实验六 用窗函数设计FIR数字滤波器.doc

实验六 用窗函数法设计FIR数字滤波器 一，实验目的 掌握窗函数法设计FIR数字滤波器的原理及具体方法； 深入理解吉布斯现象，理解不同窗函数的特点。 二，实验内容 例6.1 利用fir1函数和矩形窗设计一个N=51，截止频率ω=0.5π的低通滤波器，画出幅频特性。 MATLAB程序： clear N=51;wc=0.5; h=fir1(50,wc,boxcar(N)) [H,W]=freqz(h,1) plot(W/pi,abs(H)); title(‘矩形窗振幅特性/dB’); xlabel(‘相对频率’);ylabel(‘H(w)’) 例6.2 利用fir1函数和布莱克曼窗设计一个N=51，截止频率为=0.3π，=0.4π的带通滤波器。 MATLAB程序： clear N=51;wc=[0.3,0.4]; h=fir1(50,wc,bandpass,blackman(N)) [H,W]=freqz(h,1) plot(W/pi,abs(H)); title(布莱克曼窗带通振幅特性/dB); xlabel(相对频率);ylabel(H(w)) 6.3.1窗函数法设计低通数字滤波器， （1）N=26，分别利用矩形窗，汉宁窗和布莱克曼窗设计该滤波器，且滤波器具有线性相位。绘出脉冲响应及滤波器的频率响应； （2）增加N，观察过渡带和最大肩峰值的变化。 MATLAB程序： (1)N=26,分别利用矩形窗汉宁窗和布莱克曼窗设计该滤波器 N=26;wc=0.4; n=0:25; h1=fir1(25,wc,boxcar(N)); [H1,W]=freqz(h1,1); h2=fir1(25,wc,hanning(N)); [H2,W]=freqz(h2,1); h3=fir1(25,wc,blackman(N)); [H3,W]=freqz(h3,1); subplot(331);stem(n,h1,.); title(矩形窗对应的脉冲响应h1[n]); xlabel(相对频率);ylabel(h1[n]); subplot(332);plot(W/pi,abs(H1)); title(矩 形 窗 设 计 滤 波 器 的 幅 度 响 应/dB) xlabel(相 对 频 率);ylabel(H1(w)) subplot(333);plot(W/pi,angle(H1)); title(矩 形 窗 设 计 滤 波 器 的 相 位 响 应) xlabel(相 对 频 率);ylabel(DEGREES); subplot(334),stem(n,h2,.); title(汉宁窗对应的脉冲响应h2[n]) xlabel(相对频率);ylabel(h2[n]); subplot(335); plot(W/pi,abs(H2)); title(汉 宁 窗 设 计 滤 波 器 的 幅 度 响 应/dB) xlabel(相 对 频 率);ylabel(H2(w)); subplot(336); plot(W/pi,angle(H2)); title(汉 宁 窗 设 计 滤 波 器 的 相 位 响 应) xlabel(相 对 频 率);ylabel(DEGREES); subplot(337); stem(n,h3,.); title(布 莱 克 曼 窗 对 应 的 脉 冲 响 应h3[n]) xlabel(相 对 频 率);ylabel(h3[n]); subplot(338); plot(W/pi,abs(H3)); title(布莱克曼窗设计滤波器的幅度响应/dB) xlabel(相对频率);ylabel(H3(w)); subplot(339); plot(W/pi,angle(H3)); title(布 莱 克 曼 窗 设 计 滤 波 器 的 相 位 响 应) xlabel(相 对 频 率);ylabel(DEGREES (2)增加N=51设计程序如下 N=51;wc=0.4; n=0:50; h1=fir1(50,wc,boxcar(N)); [H1,W]=freqz(h1,1); h2=fir1(50,wc,hanning(N)); [H2,W]=freqz(h2,1); h3=fir1(50,wc,blackman(N)); [H3,W]=freqz(h3,1); subplot(331);stem(n,h1,.); title(矩 形 窗 对 应 的 脉 冲 响 应h1[n]); xlabel(相 对 频 率);ylabel(h1[n]); subplot(332); plot(W/pi,abs(H1)); ti