%《數字信號處理教程——MATLAB釋義與實現》第四章例4.6.2程序hc462
%  連續信號頻譜計算中含直流和多個正余弦分量的算例
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
% Program 4.6.2 
wm=62.83;T=pi/wm;       % 選擇采樣周期T
N0=[8,20];          % 設定兩種N值      
for r=1:2                       % 作兩次循環計算
    N=N0(r);DOmega=2*pi/(N*T);    % 取N,并求出模擬頻率分辨率
    k=[0:N-1]+eps;            % 求頻率下標，偏移微量
    X=sin(0.275*k*DOmega)./sin(0.025*k*DOmega); % 算出頻譜Xd
    x=ifft(X);                 % 求Xd的IDFT，得出x(n)
    subplot(1,2,r),stem([0:N-1]*T,x,'.')    % 繪圖
    set(gcf,'color','w');              % 置圖形背景色為白
end