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