实验四FIR数字滤波器的设计.docVIP

《数字信号处理》 （2010-2011学年第1学期 ） 论文成绩： 实验四 FIR数字滤波器的设计 学生姓名：闫春遐 所在院系：电子信息工程学院自动化系 年级专业：2008级自动化系 学 号 指导教师：王亮 完成日期：2010年11月6日 实验四 FIR数字滤波器的设计 一、实验目的 （1）掌握用窗函数法、频率采样法及优化设计法设计FIR滤波器的原理及方法，熟悉相应的MATLAB编程。 （2）熟悉线性相位FIR滤波器的幅频特性和相频特性。 （3）了解各种不同窗函数对滤波器性能的影响。 二、实验内容 （1）=45，计算并画出矩形窗、汉明窗、布莱克曼窗的归一化的幅度谱，并比较各自的主要特点。 解答： N=45; window=boxcar(N); [h1 w]=freqz(window,1); subplot(1,2,1);stem(window,.);xlabel(n);ylabel(幅值); subplot(1,2,2);plot(w/pi,20*log(abs(h1)/abs(h1(1))));grid;xlabel(归一化频率);ylabel(幅值);title(矩形窗归一化的幅度谱); window=hamming(N); [h2 w]=freqz(window,1); subplot(1,2,1);stem(window,.);xlabel(n);ylabel(幅值); subplot(1,2,2);plot(w/pi,20*log(abs(h2)/abs(h2(1))));grid;xlabel(归一化频率);ylabel(幅值);title(汉明窗归一化的幅度谱); window=blackman(N); [h3 w]=freqz(window,1); subplot(1,2,1);stem(window,.);xlabel(n);ylabel(幅值); subplot(1,2,2);plot(w/pi,20*log(abs(h3)/abs(h3(1))));grid;xlabel(归一化频率);ylabel(幅值);title(布莱克曼窗归一化的幅度谱); （2）=15，带通滤波器的两个带通边界分别是=0.3，=0.5。用汉宁窗设计此线性相位带通滤波器，观察它的实际3 dB和20 dB带宽。=45，重复这一设计，观察幅频和相位特性的变化，注意长度变化的影响。 解答： wc=[0.3 0.5]; N1=15; b1=fir1(N1,wc,bandpass,hanning(N1+1)); [h1 w]=freqz(b1,1); subplot(1,2,1);plot(w/pi,20*log(abs(h1)));grid;xlabel(归一化频率);title(用汉宁窗设计N=15); N2=45; b2=fir1(N2,wc,bandpass,hanning(N2+1)); [h2 w]=freqz(b2,1); subplot(1,2,2);plot(w/pi,20*log(abs(h2)));grid;xlabel(归一化频率);title(N=45); （3）分别改用矩形窗和布莱克曼窗，设计（2）中的带通滤波器，观察并记录窗函数对滤波器幅频特性的影响，比较三种窗的特点。 解答： wc=[0.3 0.5]; N1=15; b1=fir1(N1,wc,bandpass,boxcar(N1+1)); [h1 w]=freqz(b1,1); subplot(1,2,1);plot(w/pi,20*log(abs(h1)));grid;xlabel(归一化频率);title(用矩形窗设计N=15); N2=45; b2=fir1(N2,wc,bandpass,boxcar(N2+1)); [h2 w]=freqz(b2,1); subplot(1,2,2);plot(w/pi,20*log(abs(h2)));grid;xlabel(归一化频率);title(N=45); wc=[0.3 0.5]; N1=15; b1=fir1(N1,wc,bandpass,blackman(N1+1)); [h1 w]=freqz(b1,1); subplot(1,2,1);plot(w/pi,20*log(abs(h1)));grid;xlabel(归一化频率);title(用布莱克曼窗设计N=15); N2=45; b2=fir1(N2,wc,