idft.m

数字信号分析处理方面的好教材,有DSP的MATLAB源程序.

function [xk] = idft(Xp,N)
% 计算逆离散付利叶变换
% -----------------------------------
% [xk] = idft(Xp,N)
% xk = 在0 <= k <= N-1间的N点有限长度序列 
% Xp = 在0 <= k <= N-1间的DFT 系数数组 
%  N = DFT的长度
%
k = [0:1:N-1];                       % k的行向量
p = [0:1:N-1];                       % p的行向量
WN = exp(-j*2*pi/N);                 % Wn 因子
kp = k'*p;                           % 产生一个含kp值的N 乘 N维矩阵
WNkp = WN .^ (-kp);                  % IDFT 矩阵
xk = (Xp * WNkp)/N;                  % IDFT 的行向量