hc456.m

书籍中matlab第二部分代码

%《数字信号处理教程——MATLAB释义与实现》第四章例4.5.6程序hc456
% T,N,L变化时频点集的对应关系
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
T0=[0.1,0.05,0.05,0.1];
N0=[16,16,32,32];
for r=1:4
    T=T0(r);N=N0(r);                        % 设定T和N的值
    L=T*N;DOmega=2*pi/L;
    k=floor(-(N-1)/2:(N-1)/2);                 % 设定正负对称频率下标范围
    subplot(4,1,r),plot(k*DOmega,0.1*ones(1,N),'.')     % 画出频点集
    xlabel('频率')
    axis([-70,70,-1,1])
    set(gca,'Xtick',[-pi/T,pi/T]),grid on,         % 标注乃奎斯特频率范围
    switch r,                         % 用switch语句写子图标题
    case 1,title('T=0.1,N=16,L=1.6')
    case 2,title('T=0.05,N=16,L=0.8')
    case 3,title('T=0.05,N=32,L=1.6')
    case 4,title('T=0.1,N=32,L=3.2')
    otherwise
    end
end
set(gcf,'color','w')              % 置图形背景色为白