实验之FIR数字滤波器设计实验内容1

实验内容1：用升余弦窗设计一线性相位低通FIR数字滤波器，截止频率。窗口长度N=15,33。要求在两种窗口长度情况下，分别求出，打印出相应的幅频特性及相频特性曲线,观察3dB带宽和20dB带宽。总结窗口长度N对滤波特性的影响。设计低通FIR数字滤波器时，一般以理想低通滤波特性为逼近函数，即：

A.程序清单：

1、N=15;a=0.5*(N-1);n=0:N-1;wn=pi/4;w=sin(wn*(n-a))./(pi忮氽阝另*(n-a));if rem(N,2)~=0;w(a+1)=wn/pi;end;w1=hamming(N);h=w.*w1';subplot(3,2,1)stem(n,h,'.')title('图一 N=15 hamming窗设计的h(n)','FontSize',10,'FontWeight','bold');axis([0,14,-0.1,0.3]);ylabel('h(n)');text(0,0.15,['Made by：邱英尤'],'FontSize',10,'FontWeight','bold');gridh1=fft(h,512);w2=2*[0:511]/512;subplot(3,2,3)plot(w2,20*log10(abs(h1)))title('图二 N=15 幅频曲线','FontSize',10,'FontWeight','bold');xlabel('数字频率w/pi');ylabel('幅度/db');gridsubplot(3,2,5)plot(w2,angle(h1))title('图三 N=15 相频曲线','FontSize',10,'FontWeight','bold');xlabel('数字频率w/pi');ylabel('相角/rad');grid;N=33;a=0.5*(N-1);n=0:N-1;wn=pi/4;w=sin(wn*(n-a))./(pi*(n-a));if rem(N,2)~=0;w(a+1)=wn/pi;end;w1=hamming(N);h=w.*w1';subplot(3,2,2);stem(n,h,'.');title('图四 N=32 hamming窗设计的h(n)','FontSize',10,'FontWeight','bold');axis([0,32,-0.1,0.3]);grid;h1=fft(h,512);w2=2*[0:511]/512;subplot(3,2,4);plot(w2,20*log10(abs(h1)))title('图五 N=32 幅频曲线','FontSize',10,'FontWeight','bold');xlabel('数字频率w/pi');ylabel('幅度/db');gridsubplot(3,2,6)plot(w2,angle(h1))title('图六 N=32 相频曲线','FontSize',10,'FontWeight','bold');xlabel('数字频率w/pi');ylabel('相角/rad');grid;

B.运行结果：