hc451.m

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

%《数字信号处理教程——MATLAB释义与实现》第四章例4.5.1程序hc451
% 连续信号频谱计算中的T,N,L选择过程
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% 连续信号频谱计算中的T,N,L选择过程
T0=[2,1,0.5,0.1,0.1,0.1];   % 各次计算拟采用的T,编成向量T0
L0=[10,10,10,10,20,40];     % 各次计算拟采用的N,编成向量N0
for r=1:6                   % 循环计算六次
    T=T0(r);N=L0(r)/T0(r);      % 根据计算顺序选用T及N
    D=2*pi/(N*T);              % 频率分辨率
    n=0:N-1;x=exp(-0.1*n*T);    % 给出样本序列x
    Xa=T*fftshift(fft(x));       % 求x的FFT，移到对称位置
    k=floor(-(N-1)/2:(N-1)/2);  % 乃奎斯特频率下标
    [r,Xa(1)]                    % 显示乃奎斯特边界处的频谱X
    subplot(3,2,r),plot(k*D,abs(Xa))    % 在位置r处绘图
    set(gca,'Xtick',[-pi/T,-0.5*pi/T,0.5*pi/T,pi/T]),grid on
    set(gcf,'color','w')              % 置图形背景色为白
    xlabel('\Omega'),ylabel('|Xa(\Omega)|')
    switch r                         % 用switch语句写子图标题
    case 1, title('T=2,N=5')
    case 2, title('T=1,N=10')
    case 3, title('T=0.5,N=20')
    case 4, title('T=0.1,N=100')
    case 5, title('T=0.1,N=200')
    case 6, title('T=0.1,N=400')
    otherwise 
    end    
    axis([-3,3,0,10])           % 统一坐标，以便比较
end
set(gcf,'color','w')