%% DFT transform
% Tom Wada
% 2011/6/20
clear all;
M=16;
n=(0:M-1);
%Šî’êŠÖ�”
q=zeros(M,M);
r(1,:)=exp(-1j*2*pi*(1-1)*n/M);
for l=2:M
    r(l,:)=exp(-1j*2*pi*(l-1)*n/M);
end
%’¼Œðƒ`ƒFƒbƒN
InnerProd=zeros(M,M);
for l=1:M
    for n=1:M
        InnerProd(l,n)=r(l,:)*r(n,:)'/M;
    end
end


