%N=16有限序列的DFT、IDFT結果
%DFT
N=16;
n=0:1:N-1;          %時域采樣
xn=sin(n*pi/8)+sin(n*pi/4);
k=0:1:N-1;          %頻域采樣
Wn=exp(-j*2*pi/N);
Nk=n'*k;
WnNk=Wn.^Nk;
Xk=xn*WnNk;
subplot(2,1,1),stem(n,xn);
subplot(2,1,2),stem(k,abs(Xk));
pause

%IDFT
N=16;
n=0:1:N-1;          %時域采樣
k=0:1:N-1;          %頻域采樣
Wn1=exp(j*2*pi/N);
Nk=n'*k;
WnNk1=Wn1.^Nk;
xn1=(Xk*WnNk1)./N;
subplot(2,1,1),stem(k,abs(Xk));
subplot(2,1,2),stem(n,real(xn1));