function [xn] = idft(Xk,N)% Computes Inverse Discrete Transform% -----------------------------------% [xn] = idft(Xk,N)% xn = N-point sequence over 0 <= n <= N-1% Xk = DFT coeff. array over 0 <= k <= N-1%  N = length of DFT%n = [0:1:N-1];                       % row vector for nk = [0:1:N-1];                       % row vecor for kWN = exp(-j*2*pi/N);                 % Wn factornk = n'*k;                           % creates a N by N matrix of nk valuesWNnk = WN .^ (-nk);                  % IDFT matrixxn = (Xk * WNnk)/N;                  % row vector for IDFT values