hc329.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第三章例3.2.9程序hc329
% 求理想数字低通滤波器对应的脉冲序列
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear,close all
wc=input('wc= （范围从0到pi)');             % 低通滤波器的截止频率
n=[-10:10]+1e-10;
hd=sin(n*wc)./(n*pi);                       % 理论的反变换脉冲序列
subplot(1,2,1),
plot([-pi,-wc,-wc,wc,wc,pi],[0,0,1,1,0,0])  % 画出理想频率响应
xlabel('频率（1/秒）'),ylabel('幅度')
axis([-pi,pi,-0.1,1.1]),grid on
subplot(1,2,2),stem(n,hd),grid on           % 画出对应的脉冲序列
xlabel('n'),ylabel('序列')
axis([-10,10,-0.1*wc,0.4*wc])               % 确定适当的坐标比例尺
set(gcf,'color','w')                        % 置图形背景色为白