生成离散傅里叶变换的矩阵的程序

生成N点的离散傅里叶变换的矩阵如下：

%% matlab离散傅里叶变换矩阵
%% Vsesion:1.0 Written by zhenhuaLiu@ 2021.11.16 HIT ATCI
function [WNnk] = dftWnkGenerate(N)

 
n = [0:1:N-1]; % row vector for n
k = [0:1:N-1]; % row vector for k
WN = exp(-1i*2*pi/N);%matlab提示要增强稳定性，使用li代替i
nk = n'*k;
WNnk = WN .^ nk;   %DFT matrix