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

ylabel(h(n));

text(0,0.15,,FontSize,10,FontWeight,bold);

grid

h1=fft(h,512);

w2=2*/512;

subplot(3,2,3)

plot(w2,20*log10(abs(h1)))

title(图二 N=15 幅频曲线,FontSize,10,FontWeight,bold);

xlabel(数字频率w/pi);ylabel(幅度/db);

grid

subplot(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();

grid;

h1=fft(h,512);

w2=2*/512;

subplot(3,2,4);

plot(w2,20*log10(abs(h1)))

title(图五 N=32 幅频曲线,FontSize,10,FontWeight,bold);

xlabel(数字频率w/pi);ylabel(幅度/db);

grid

subplot(3,2,6)

plot(w2,angle(h1))

title(图六 N=32 相频曲线,FontSize,10,FontWeight,bold);

xlabel(数字频率w/pi);ylabel(相角/rad);

grid;

B.运行结果：1/32/33/3滤波器