idft.m

数字信号处理及MATLAB实现 具体的例子与码源 很能帮助初学者学习

function [xn] = idft(Xk,N)
% 计算逆离散付利叶变换
% -----------------------------------
% [xn] = idft(Xk,N)
% xn = 在0 <= k <= 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 的行向量