wc=0.25*pi;N=33;
alpha=(N-1)/2;
n=0:1:N-1;
m=n-alpha+eps;
hd=sin(wc*m)./(pi*m);
wd1=hanning(N)';b1=hd.*wd1;
wd2=hamming(N)';b2=hd.*wd2;
wd3=hanning(N)';b3=hd.*wd3;
wd4=hanning(N)';b4=hd.*wd4;
[H1,w]=freqz(b1,1);
[H2,w]=freqz(b2,1);
[H3,w]=freqz(b3,1);
[H4,w]=freqz(b4,1);
subplot(3,1,1)
plot(w,abs(H1),':',w,abs(H2),':',w,abs(H3),':',w,abs(H4),':')
legend('漢寧窗低通濾波器','海明窗低通濾波器','布來克曼窗低通濾波器','凱澤窗低通濾波器')
subplot(3,1,2)
plot(w/pi,20*log10(abs(H1)),':',w/pi,20*log10(abs(H2)),':',w/pi,20*log10(abs(H3)),':',w/pi,20*log10(abs(H4)),':')
legend('漢寧窗幅度響應（dB）','海明窗幅度響應（dB）','布來克曼窗幅度響應（dB）','凱澤窗幅度響應（dB）')
subplot(3,1,3)
plot(n,b1,':',n,b2,':',n,b3,':',n,b4,':')
legend('漢寧窗h(n)','海明窗h(n)','布來克曼窗h(n)','凱澤窗h(n)')