4.txt

采用频率取样设计法设计FIR数字低通滤波器

采用频率取样设计法设计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('分贝');