MATLAB从零开始实现离散傅里叶变换DFT

一、基础目标

离散傅里叶变换（DFT）是信号处理领域的基石，它将信号从时域转换到频域，让我们能观察信号的频率成分。下面我们深入探讨其核心原理，并在MATLAB中从零开始实现它。

二、DFT的核心原理

DFT的数学定义如下：对于一个长度为N的离散时间信号x[n]，其DFT变换X[k]由以下公式给出：

X [ k ] = ∑ n = 0 N − 1 x [ n ] ⋅ e − j 2 π N k n , k = 0 , 1 , 2 , . . . , N − 1 X[k]=\sum_{n=0}^{N-1}x[n]\cdot e^{-j\frac{2\pi}{N}kn},\quad k=0,1,2,...,N-1X[k]=∑n=0N−1​x[n]⋅e−jN2π​kn,k=0,1,2,...,N−1

其中：

• x [ n ] x[n]x[n]是时域信号的第n个采样值
• X [ k ] X[k]X[k]是频域信号的第k个频率分量（复数值）
• N NN是信号的采样点数
• j jj是虚数单位（在MATLAB中用1 i 1i1i或1 j 1j1j表示）

逆离散傅里叶变换（IDFT）则可以将频域信号恢复回时域：

x [ n ] = 1 N ∑ k = 0 N − 1 X [ k ] ⋅ e j 2 π N k n , n = 0 , 1 , 2 , . . . , N − 1 x[n]=\frac{1}{N}\sum_{k=0}^{N-1}X[k]\cdot e^{j\frac{2\pi}{N}kn},\quad n=0,1,2,...,N-1x[n]=N1​∑