%《數(shù)字信號(hào)處理教程——MATLAB釋義與實(shí)現(xiàn)》第四章例4.5.4程序hc454
% 連續(xù)信號(hào)頻譜計(jì)算中截?cái)嗟恼嘞倚盘?hào)頻譜的算例
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
N0=[50,100,500,628];
for r=1:4
    N=N0(r);
    T=0.1;n=1:N;        % 原始數(shù)據(jù)
    DOmega=2*pi/(N*T);                  % 頻率分辨率
    x=cos(5*n*T);        %  生成有限長的正弦序列
    X=T*fftshift(fft(x));X(1)       % 求x的FFT，移到對(duì)稱位置
    k=floor(-(N-1)/2:(N-1)/2);      % 頻率下標(biāo)向量
    subplot(2,2,r),plot(k*DOmega,abs(X))         % 繪圖
    v=axis;axis([-10,10,v(3),v(4)])
    xlabel('\Omega'),ylabel('|X(\Omega)|')
    switch r                         % 用switch語句寫子圖標(biāo)題
    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');              % 置圖形背景色為白