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是对角阵,那么左乘跟右乘是一样的。