%《數字信號處理教程——MATLAB釋義與實現》第四章例4.5.5程序hc455
%  連續信號頻譜計算中含直流和多個正余弦分量的算例
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
for i= 1:2
    N=input('N= ');
    T=0.2;n=1:N;                    % 原始數據
    DOmega=2*pi/(N*T);                   % 頻率分辨率
    x=1+2*sin(0.2*pi*n*T)-3*cos(pi*n*T);        %  生成有限長的正弦序列
    X=T*fftshift(fft(x));X(1);      % 求x的FFT，移到對稱位置
    k=floor(-(N-1)/2:(N-1)/2);      % 頻率下標向量
    subplot(1,2,i),plot(k*DOmega,abs(X))         % 繪圖
    title('N= ,L=')
    xlabel('\Omega'),ylabel('|X(\Omega)|')
    v=axis;axis([-10,10,v(3),v(4)])  
end
set(gcf,'color','w');               % 置圖形背景色為白