idft.m

matlab应用实例。如：FIR、IIR数字滤波器设计等。

function [xn]=idft(xk,N)
%this program computes Inverse Discrete Fourier Transform
%its format is
%xn= IDFT(Xk,N)
%N is length of IDFT
%xk is DFT coeff. Array over 0<=k<=N-1
%xn is N-point sequence over 0<=k<=N-1

n=[0:1:N-1];%n的行向量
k=[0:1:N-1];%k的行向量
WN=exp(-j*2*pi/N);%旋转因子
nk=n'*k;%产生一个含nk值的N乘N维矩阵
WNnk=WN.^(-nk);%I DPT矩阵
xn=(xk*WNnk)/N;%IDFT系数的行向量