fgp742.m

digital signal processing-MATLAB

%《数字信号处理教程——MATLAB释义与实现》第七章演示程序fgp742
% 给定理想低通滤波器的13个频点图
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
y=[ones(1,3),zeros(1,8),ones(1,2)];     % 13个频点的幅度
nx=0:12;                                % 分成的频谱自变量序号
x=nx/13*2*pi;                           % 13点的频率值
subplot(2,1,1)                          % 绘理想频点频谱图
plot(x,y,'o')
line([0,pi/3,pi/3,5*pi/3,5*pi/3,2*pi],[1,1,0,0,1,1])  % 绘理想的连续频谱图
line([0,2*pi],[0,0])                    % 画水平基线
line([pi,pi],[-0.1,0.1])                % 画对称中心线
axis([0,2*pi,-0.2,1.2])                 % 定坐标范围
subplot(2,1,2)
theta=-12*pi*[0:12]/13;                 % 13个频点的相位
plot(x,theta),hold on                   % 绘理想的连续相位图
plot(x,theta,'o')                       % 绘理想频点相位图
axis([0,2*pi,-12*3.5,6])                % 定坐标范围
grid
set(gcf,'color','w')                    % 置图形背景色为白