hc454.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第四章例4.5.4程序hc454
% 连续信号频谱计算中截断的正余弦信号频谱的算例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N0=[50,100,500,628];
for r=1:4
    N=N0(r);
    T=0.1;n=1:N;        % 原始数据
    D=2*pi/(N*T);                  % 频率分辨率
    x=cos(5*n*T);        %  生成有限长的正弦序列
    X=T*fftshift(fft(x));X(1)       % 求x的FFT，移到对称位置
    k=floor(-(N-1)/2:(N-1)/2);      % 频率下标向量
    subplot(2,2,r),plot(k*D,abs(X))         % 绘图
    v=axis;axis([-10,10,v(3),v(4)])
    xlabel('\Omega'),ylabel('|X(\Omega)|')
    switch r                         % 用switch语句写子图标题
    case 1, title('(a)N=50,L=5')
    case 2, title('(b)N=100,L=10')
    case 3, title('(c)N=500,L=50')
    case 4, title('(d)N=628,L=62.8')
    otherwise 
    end    
end
set(gcf,'color','w');              % 置图形背景色为白