%Samp3_7
clf
N=100; dt=1;   %數據長度為100，采樣間隔為1s
n=0:N-1;t=n*dt;  %給出時間序列
xn1=cos(2*pi*0.24*t);   %第一個振動
xn2=cos(2*pi*0.12*t);   %第二個振動
Xk1=dfs(xn1,N);    %第一個振動的Fourier變換
Xk2=dfs(xn2,N);    %第二個振動的Fourier變換
magXk1=abs(Xk1);phaXk1=angle(Xk1);  %第一個振動的振幅、相位
magXk2=abs(Xk2);phaXk2=angle(Xk2); %第二個振動的振幅、相位
k=0:length(magXk1)-1;
subplot(3,1,1),plot(k/(N*dt),magXk1*2/N);  %繪制第一個振動的振幅譜
ylabel('振幅');
title('第一個振動的Fourier變換');
k=0:length(magXk2)-1;
subplot(3,1,2),plot(k/(N*dt),magXk2*2/N);
ylabel('振幅');
title('第二個振動的Fourier變換');
Xk=dfs(xn1+xn2,N);  %兩個振動合成的Fourier變換
magXk=abs(Xk);phaXk=angle(Xk);   %合成振動的振幅和相位
k=0:length(magXk)-1;
subplot(3,1,3),plot(k/(N*dt),magXk*2/N);  %繪制合成振動的振幅和相位
xlabel('頻率/Hz');ylabel('振幅');
title('合成振動的Fourier變換');