chuang.m

用MATLAB实现数字信号处理的程序

%设计FIR低通滤波器
wc=0.2*pi;
N=11;
tao=(N-1)/2;
n=[0:(N-1)];
m=n-tao+eps;
hd=sin(wc*m)./(pi*m);
wd1=boxcar(N)';b1=hd.*wd1;
wd2=hanning(N)';b2=hd.*wd2;
wd3=blackman(N)';b3=hd.*wd3;
[H1,w]=freqz(b1,1);
[H2,w]=freqz(b2,1);
[H3,w]=freqz(b3,1);
subplot(1,2,1),plot(w,abs(H1),'r-.',w,abs(H2),'g*',w,abs(H3),'b:');
title('窗函数'),xlabel('w'),ylabel('H(w)');
legend('矩形窗','汉宁窗','布莱克曼窗');
subplot(1,2,2),plot(w,20*log10(abs(H1)),'r-.',w,20*log10(abs(H2)),'g*',w,20*log10(abs(H3)),'b:');
title('幅频特性'),xlabel('w'),ylabel('dB');
legend('矩形窗','汉宁窗','布莱克曼窗');