DFT的matlab代码

以下是一个简单的DFT Matlab代码示例：

function X = myDFT(x)
% 计算信号的DFT
% 输入：x – 输入信号向量
% 输出：X - DFT值向量

N = length(x); % 信号长度
n = 0:N-1;
k = n';
WN = exp(-2*pi*1i/N); % 旋转因子
nk = n*k;
W = WN.^nk; % DFT矩阵
X = x*W; % DFT值向量

end

这个函数接受一个信号向量x作为输入，并返回其DFT值向量X。DFT计算通过矩阵乘法实现。代码中，我们首先计算了信号长度N，然后创建了一个包含0到N-1的向量n和一个包含0到N-1的向量k。接着，我们计算了旋转因子WN，并使用矩阵乘法将其提升到n和k的元素级别。这给出了一个DFT矩阵W，我们用它来计算信号的DFT值向量X。