《数字信号处理实验报告——用窗函数法设计FIR数字滤波器》由会员分享，可在线阅读，更多相关《数字信号处理实验报告——用窗函数法设计FIR数字滤波器（7页珍藏版）》请在金锄头文库上搜索。
1、实验四用窗函数法设计FIR数字滤波器实验目的(1) 掌握用窗函数法设计FIR数字滤波器的原理和方法。(2) 熟悉线性相位FIR数字滤波器特性。(3) 了解各种窗函数对滤波特性的影响。二、实验原理、滤波器的理想频率响应函数为Hd(ej)，则其对应的单位脉冲响应为(ej)ejndw窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼hd(n)。由于hd(n)往往是无限长序列，且是非因果的，所以用窗函数。w(n)将hd(n)截断，并进行加权处理：h(n)=hd(n)w(n)h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数H(ej)为用窗函数法设计的滤波器性能取决于窗函
2、数w(n)的类型及窗口长度N的取值。设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。三、实验内容1. MATALAB程序任务一N1=15;N2=33;b1=fir1(14,1/4,hanning(15);b2=fir1(32,1/4,hanning(33);H1,W=freqz(b1,1);H1_db=20*log10(abs(H1);magH1=abs(H1);phaH1=angle(H1);H2,W=freqz(b2,1);H2_db=20*log10(abs(H2);magH2=abs(H2);phaH2=angle(H2);figure(1);sub
3、plot(2,1,1);stem(b1);title(N=15时,汉宁窗得到的FIR滤波器脉冲响应);subplot(2,1,2);stem(b2);title(N=33时,汉宁窗得到的FIR滤波器脉冲响应);figure(2);subplot(3,1,1);plot(W/pi,magH1);title(N=15时汉宁窗得到的FIR滤波器的幅频特性)subplot(3,1,2);plot(W/pi,H1_db);title(N=15时汉宁窗得到的FIR滤波器的对数幅频特性)subplot(3,1,3);plot(W/pi,phaH1);title(N=15时汉宁窗得到的FIR滤波器的相频特性)
4、figure(3);subplot(3,1,1);plot(W/pi,magH2);title(N=33时汉宁窗得到的FIR滤波器的幅频特性)subplot(3,1,2);plot(W/pi,H2_db);title(N=33时汉宁窗得到的FIR滤波器的对数幅频特性)subplot(3,1,3);plot(W/pi,phaH2);title(N=33时汉宁窗得到的FIR滤波器的相频特性)任务二N=33;b1=fir1(32,1/4,boxcar(33);b2=fir1(32,1/4,hanning(33);b3=fir1(32,1/4,bartlett(33);b4=fir1(32,1/4,b
5、lackman(33);H1,W=freqz(b1,1);H1_db=20*log10(abs(H1);magH1=abs(H1);phaH1=angle(H1);H2,W=freqz(b2,1);H2_db=20*log10(abs(H2);magH2=abs(H2);phaH2=angle(H2);H3,W=freqz(b3,1);H3_db=20*log10(abs(H3);magH3=abs(H3);phaH3=angle(H3);H4,W=freqz(b4,1);H4_db=20*log10(abs(H4);magH4=abs(H4);phaH4=angle(H4);figure(1
6、);subplot(2,2,1);stem(b1);title(矩形窗得到的FIR滤波器脉冲响应)subplot(2,2,2);stem(b2);title(汉宁窗得到的FIR滤波器脉冲响应)subplot(2,2,3);stem(b3);title(三角窗得到的FIR滤波器脉冲响应)subplot(2,2,4);stem(b4);title(布拉克曼窗得到的FIR滤波器脉冲响应)figure(2);subplot(2,2,1);plot(W/pi,magH1);title(矩形窗得到的FIR滤波器幅频特性)subplot(2,2,2);plot(W/pi,magH2);title(汉宁窗得到
7、的FIR滤波器幅频特性)subplot(2,2,3);plot(W/pi,magH3);title(三角窗得到的FIR滤波器幅频特性)subplot(2,2,4);plot(W/pi,magH4);title(布拉克曼窗得到的FIR滤波器幅频特性)figure(3);subplot(2,2,1);plot(W/pi,H1_db);title(矩形窗得到的FIR滤波器相频特性)subplot(2,2,2);plot(W/pi,H2_db);title(汉宁得到的FIR滤波器相频特性)subplot(2,2,3);plot(W/pi,H3_db);title(三角窗得到的FIR滤波器相频特性)su
8、bplot(2,2,4);plot(W/pi,H4_db);title(布拉克曼得到的FIR滤波器相频特性)figure(4);subplot(2,2,1);plot(W/pi,phaH1);title(矩形窗得到的FIR滤波器相频特性)subplot(2,2,2);plot(W/pi,phaH2);title(汉宁窗得到的FIR滤波器相频特性)subplot(2,2,3);plot(W/pi,phaH3);title(三角窗得到的FIR滤波器相频特性)subplot(2,2,4);plot(W/pi,phaH4);title(布拉克曼窗得到的FIR滤波器相频特性)2. 实验波形图任务一g时溟
9、宁唐f早到的FIR滤锻器的对故幅豹特性任务二三*宙得到的FI磁试盟对数幅频特性10-X40DQIQ20.3040S0.60?030.91四、实验结论1. N的大小决定了窗谱的主瓣宽度，N越大，窗谱的主瓣宽度越大2. 最小阻带衰减只有窗行决定，不受N的影响，过渡带宽度与N和窗形都有关,N越大，过渡带宽越小3由实验可知滤波特性：布拉克曼窗汉宁窗三角窗矩形窗五、思考题(1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器？写出设计步骤。答：根据A，s_p计算出所要设计的滤波器的过渡带查表找到满足的阻带最小衰减A所需的窗函数s 利用已选的窗函数的过渡带宽除以A得到低通滤波器单位脉冲响应的长度N 根据hC)hd仏0求出h(n)n0求出响应函数H(ej(2)如果要求用窗函数法设计带通滤波器，且给定上、下边带截止频率为卩和叫,试求理想带通的单位脉冲响应hd(n)。答：由s+p算出，求出滤波器过渡带宽A-和滤波器长度N,c2c12最后由h(n)d1(n-t)sin(n-T)nht或h(n)=巳,ntdT求出h(n)，其中dT=
