Clear;close all;
N=21;wc=pi/4;                          %理想低通濾波器參數(shù)
n=0:N-1;r=(N-1)/2;
hdn=sin(wc*(n-r))/pi./(n-r);      %計(jì)算理想低通單位脈沖響應(yīng) 
if rem(N,2)~=0 hdn(r+1)=wc/pi;end  %N為奇數(shù)時(shí)，處理n=r點(diǎn)的0/0型
wn1=boxcar(N);                     %矩形窗
hn1=hdn.*wn1';	    %加窗
wn2=hamming(N);	%hamming窗
hn2=hdn.*wn2';	%加窗
w=2*[0:511]/512;
hw1=fft(hn1,512);
hw2=fft(hn2,512);
subplot(2,2,1);                    %作圖
stem(n,hn1,'.');grid;
xlabel('n'); ylabel('h(n)');
title('矩形窗設(shè)計(jì)的h(n) ');        	%圖形劃分及標(biāo)注語(yǔ)句
subplot(2,2,2);
plot(w,20*log10(abs(hw1)));grid;
xlabel('w/pi'); ylabel('Magnitude(dB)');
title('幅頻特性（dB）');
subplot(2,2,3);
 stem(n,hn2,'.');grid;
xlabel('n'); ylabel('h(n)');
title('hamming窗設(shè)計(jì)的h(n) ');
subplot(2,2,4);
plot(w,20*log10(abs(hw2)));grid;
xlabel('w/pi'); ylabel('Magnitude(dB)');
title('幅頻特性（dB）');