%Samp3_9
N=100;dt=1;   %數據個數與采樣間隔
n = 0:N-1;t=n*dt;   %時間序列值
x = cos(pi*t/2);    %正弦信號
k = 0:N-1; f =k/(N*dt);        %頻率序列值 
X = dfs(x,N);    %  求解正弦信號的Fourier變換
realX=real(X);    %取X的實部
imagX=imag(X);    %取X的虛部
Xrealshft=cirshftt(realX,10,N);  %將X的實部信號移位10
Ximagshft=cirshftt(imagX,10,N); %將X的虛部信號移位10
Y=idfs(Xrealshft+j*Ximagshft,N); %將移位后信號結合為復數進行逆變換
subplot(2,2,1),plot(t,real(Y));ylim([-1 1])  %Y的實部
xlabel('時間/s')
title('移位后信號頻率域的實部')
subplot(2,2,2),plot(t,imag(Y));   %繪出Y的虛部
xlabel('時間/s')
title('移位后信號頻率域的虛部')
subplot(2,2,3),plot(t,x.*real(exp(j*2*pi*n*10/N))); %繪x(n) 的實部
xlabel('時間/s')
title('x*exp(j*2*pi*n*10/N)的Fourier變換的實部')
subplot(2,2,4),plot(t,x.*imag(exp(j*2*pi*n*10/N)));%繪x(n) 的虛部
xlabel('時間/s')
title(' x*exp(j*2*pi*n*10/N) 的Fourier變換的虛部')