傅里叶变换矩阵

1. matlab 的fft()函数是没有归一化的DFT

matlab 得到归一化DFT矩阵

function F=DFT(N) %xn为序列
n = [0:N-1]; %n的行向量，为1*N矩阵
k = [0:N-1]; %k的行向量,为1*N矩阵
Wn = exp(-j*2*pi/N); %常数
nk = n'*k; %将n倒置之后与矩阵k进行矩阵的代数运算，为N*N矩阵，此处发生了N*N次乘法运算
Wnnk = Wn.^nk; %将常数Wn与nk进行点幂运算,为N*N矩阵，此处发生了N*N次点幂运算
F = Wnnk / sqrt(N);

K =  64;
F=DFT(K);
h_K = randn(K, 1);
H_F = F*h_K;
H_fft = fft(h_K) ./ sqrt(K);%fft做傅里叶变换是非归一化的

得到的值有轻微差别，但看起来是一样的 

2.DFT矩阵F是对称的,即F = F的转置

3.矩阵A左乘DFT矩阵是对矩阵A的每一列做傅里叶变换，矩阵A右乘DFT矩阵是对矩阵A的每一行做傅里叶变换，如果A是对角阵，那么左乘跟右乘是一样的。