%Samp3_8
clf
N=128;dt=1;   %數據長度和采樣間隔
n=0:N-1; t=n*dt;  %時間序列
xn1=cos(0.14*pi*t);   %原始信號
subplot(3,2,1),plot(t,xn1);title('原始信號');   %繪制原始信號
xn2=cirshftt(xn1,20,N);     %循環位移20個時間單位
subplot(3,2,2),plot(t,xn2);  % 繪出循環移位后的信號
title('時移20個單位的信號')
Xk1=dfs(xn1,N);   %將原信號進行Fourier變換
magXk1=abs(Xk1);phaXk1=angle(Xk1);  %得到原信號的振幅和相位
k=0:length(magXk1)-1;
subplot(3,2,3),
plot(k/(N*dt),magXk1*2/N);   %繪制原信號的振幅譜
ylabel('振幅');        
title('原始信號的振幅譜');
subplot(3,2,4),plot(k/(N*dt),unwrap(phaXk1)),ylabel('相位角/rad')  
%繪制原信號的相位譜,unwrap為將信號解卷繞，即將相位角展開
title('原始信號的相位譜')
Xk2=dfs(xn2,N);   %繪制移位后信號的振幅譜
Xk2=Xk2./exp(-2*j*pi*k*20/N);  %將移位后的Fourier變換與 相乘
%如果沒有此語句，得出的相位譜會發生變化
magXk2=abs(Xk2);phaXk2=angle(Xk2);
k=0:length(magXk2)-1;
subplot(3,2,5),
plot(k/(N*dt),magXk2*2/N);
xlabel('頻率/Hz');ylabel('振幅');
title('移位后與exp(-j*2*pi*k*20/N)相乘的振幅譜');
subplot(3,2,6),plot(k/(N*dt),unwrap(phaXk2)),ylabel('相位角/rad')
xlabel('頻率/Hz')
title('移位后與exp(-j*2*pi*k*20/N)相乘的相位譜')