%Samp3_6
clf     %清除圖形框的內(nèi)容
N=120;dt=1;   %設置最大點數(shù)
n=0:N-1; t=n*dt;   %給出時間序列
xn=cos(2*pi*0.24*[0:99])+cos(2*pi*0.26*[0:99]);  
xn=[xn,zeros(1,N-100)];%給出原始信號的值序列
Xk=dfs(xn,N);      %對原始信號進行Fourier變換
magXk=abs(Xk);phaXk=angle(Xk);   %求出Fourier變換的振幅和相位
subplot(2,2,1),plot(t,xn); xlabel('時間/s')  %繪出原始信號
title('原始信號(N=120)');
xx=idfs(Xk,N);      %Fourier逆變換
x=real(xx);         %取變換后的實部，如做實驗可以驗證其虛部為零
subplot(2,2,2),plot(t,x),xlabel('時間/s'),title('運用Fourier逆變換得到的合成信號')
k=0:length(magXk)-1;
subplot(2,1,2),plot(k/(N*dt),magXk*2/N); %繪出Fourier變換的振幅譜，頻率采用（3-17）
%采用真實振幅[（3-23）式]繪圖
xlabel('頻率/Hz');ylabel('振幅');
title(' X(k)振幅(N=120)');