%Samp3_8
clf
N=128;dt=1;   %數(shù)據(jù)長(zhǎng)度和采樣間隔
n=0:N-1; t=n*dt;  %時(shí)間序列
xn1=cos(0.14*pi*t);   %原始信號(hào)
subplot(3,2,1),plot(t,xn1);title('原始信號(hào)');   %繪制原始信號(hào)
xn2=cirshftt(xn1,20,N);     %循環(huán)位移20個(gè)時(shí)間單位
subplot(3,2,2),plot(t,xn2);  % 繪出循環(huán)移位后的信號(hào)
title('時(shí)移20個(gè)單位的信號(hào)')
Xk1=dfs(xn1,N);   %將原信號(hào)進(jìn)行Fourier變換
magXk1=abs(Xk1);phaXk1=angle(Xk1);  %得到原信號(hào)的振幅和相位
k=0:length(magXk1)-1;
subplot(3,2,3),
plot(k/(N*dt),magXk1*2/N);   %繪制原信號(hào)的振幅譜
ylabel('振幅');        
title('原始信號(hào)的振幅譜');
subplot(3,2,4),plot(k/(N*dt),unwrap(phaXk1)),ylabel('相位角/rad')  
%繪制原信號(hào)的相位譜,unwrap為將信號(hào)解卷繞，即將相位角展開
title('原始信號(hào)的相位譜')
Xk2=dfs(xn2,N);   %繪制移位后信號(hào)的振幅譜
Xk2=Xk2./exp(-2*j*pi*k*20/N);  %將移位后的Fourier變換與 相乘
%如果沒有此語(yǔ)句，得出的相位譜會(huì)發(fā)生變化
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)相乘的相位譜')