idft.m

计DFT

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的行向量