采用頻率取樣設計法設計FIR數字低通濾波器，滿足以下指標
 ， 
 ， 
a．取N = 20，過渡帶沒有樣本。
b．取N = 40，過渡帶有一個樣本，T = 0.39。
c．取N = 60，過渡帶有兩個樣本，T1 = 0.5925，T2 = 0.1099




 N=20;alpha=(N-1)/2;l=0:N-1;wl=(2*pi/N)*l;
Hrs=[1,1,1,zeros(1,15),1,1];   %取樣樣本
Hdr=[1,1,0,0];wdl=[0,0.25,0.25,1];    %理想幅度響應
k1=0:floor((N-1)/2);k2=floor((N-1)/2)+1:N-1;
angH=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)];
H=Hrs.*exp(j*angH);
h=real(ifft(H,N));
w=[0:500]*pi/500;
H=freqz(h,l,w);
[Hr,wr,type,tao]=amplrespond(h);
subplot(221);
plot(wdl,Hdr,wl(1:11)/pi,Hrs(1:11),'o');
axis([0,1,-0.1,1.1]);grid;
title('取樣樣本');xlabel('頻率 (單位:pi)');ylabel('Hr(k)');
subplot(222);
stem(l,h,'.');
axis([0,N-1,-0.1,0.3]);grid;
title('實際脈沖響應');xlabel('n');ylabel('h(n)');
subplot(223);
plot(wr/pi,Hr,wl(1:11)/pi,Hrs(1:11),'o');
axis([0,1,-0.2,1.2]);grid;
title('振幅響應');xlabel('頻率 (單位:pi)');ylabel('Hr(w)');
subplot(224);
plot(w/pi,20*log10((abs(H)/max(abs(H)))));
axis([0,1,-50,5]);grid;
title('振幅響應: dB');xlabel('頻率 (單位:pi)');ylabel('分貝');