hc734.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第七章例7.3.4程序hc734
% 用三种窗函数设计低通滤波器
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear
wc=0.2*pi; N=11;
hd= ideallp(wc,N);		% 用wc=0.2*pi的理想低通作为逼近滤波器
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),w,abs(H2),':',w,abs(H3),'-.');% 绘幅特性
legend('矩形窗','汉宁窗','布莱克曼窗')
xlabel('\omega'),ylabel('|H(\omega)|')
subplot(1,2,2),                 % 绘分贝幅特性
plot(w,20*log10(abs(H1)),w,20*log10(abs(H2)),':',w,20*log10(abs(H3)),'-.');
legend('矩形窗','汉宁窗','布莱克曼窗')
xlabel('\omega'),ylabel('dB')
set(gcf,'color','w');                % 置图形背景色为白