function [xn]=idft(Xk,N)
%計算逆離散傅立葉變換
%------------------
%[xn]=idft(Xk,N)
%xn=在0<=n<=N-1間的N點有限持續時間序列
%Xk=在0<=k<=N-1的DFT系數數組
%N=DFT的長度
%
n=[0:1:N-1];        %n的行向量
k=[0:1:N-1];        %k的行向量
WN=exp(-j*2*pi/N);  %Wn因子
nk=n'*k;            %產生一個含nk值得N乘N維矩陣
WNnk=WN.^(-nk);     %IDFT矩陣
xn=(Xk*WNnk)/N;     %IDFT的行向量