查看“︁离散傅里叶变换”︁的源代码

{| class="wikitable" style="float:right; width:320px; margin-left:1em;"
|+ style="font-weight:bold; font-size:1.2em;" | 数学词条：离散傅里叶变换
|-
! 英文名称
| Discrete Fourier Transform (DFT)
|-
! 核心定义
| 将有限长离散时间序列映射为等长离散频域序列的数学变换，是计算机进行频谱分析的基础
|-
! 核心本质
| 对信号的离散时间傅里叶变换（DTFT）在频域进行等间隔采样，揭示信号的离散频率成分
|-
! 核心公式
| X[k] = Σ x[n]·e^(-j2πkn/N) （将时域卷积转化为频域乘积）
|-
! 根本目标
| 让计算机能够处理和运算信号的频谱，为快速傅里叶变换（FFT）提供理论原型
|}

== 1 概述 ==
'''离散傅里叶变换'''（Discrete Fourier Transform，简称 DFT）是数字信号处理中最基础、最核心的数学工具。在现实世界中，计算机无法直接处理时间和幅度都连续的模拟信号，也无法处理无限长的离散序列。DFT 的出现完美解决了这一问题，它将一个长度为 N 的有限长时域离散序列 x[n]，转换为另一个长度为 N 的频域离散序列 X[k]。

DFT 不仅在数学上严密，更具有深刻的物理意义：它告诉我们，任何有限长的离散信号，都可以被分解为 N 个不同频率的复指数（正弦/余弦）信号的加权和。这些复数权重 X[k] 就代表了信号在对应离散频率点上的幅度与相位信息。

== 2 物理本质与核心原理 ==

DFT 的建立依赖于从连续到离散的完整数学推导链条，其本质是对连续信号的“双重离散化”。

=== 2.1 从 FT 到 DFT 的演变 ===
为了让计算机处理频谱，必须经历以下三个数学步骤：
# '''时域离散化'''：对连续模拟信号进行等间隔采样，得到离散时间序列。此时，信号的频谱会从连续谱变为以 2π 为周期的周期谱（即离散时间傅里叶变换 DTFT）。
# '''时域截断'''：计算机无法存储无限长的序列，因此必须截取有限长度 N 的数据（相当于乘以矩形窗）。
# '''频域离散化'''：对周期的 DTFT 频谱在 0 到 2π 区间内进行 N 点等间隔采样。这一过程最终导出了 DFT，使其在时域和频域上都是离散且有限的。

=== 2.2 DFT 的数学表达式 ===
设 x[n] 为长度为 N 的有限长序列，其 N 点 DFT 定义为：
<center><math>X[k] = \sum_{n=0}^{N-1} x[n] \cdot e^{-j\frac{2\pi}{N}kn}, \quad k = 0, 1, \dots, N-1</math></center>
其对应的逆变换（IDFT）为：
<center><math>x[n] = \frac{1}{N} \sum_{k=0}^{N-1} X[k] \cdot e^{j\frac{2\pi}{N}kn}, \quad n = 0, 1, \dots, N-1</math></center>
其中，e^(-j2πkn/N) 被称为'''旋转因子'''（Twiddle Factor），它是 DFT 计算的核心。

=== 2.3 隐含的周期性 ===
DFT 在数学上隐含了时域和频域的周期性。虽然我们在工程中通常只处理主值区间（0 到 N-1），但在进行圆周卷积、频域采样等操作时，必须将序列视为以 N 为周期的周期序列来处理。

== 3 核心性质与定理 ==

DFT 拥有一系列完美的数学性质，这些性质是设计数字滤波器和信号分析系统的理论基石：

# '''线性性质'''：多个信号线性组合后的 DFT，等于各信号 DFT 的线性组合。
# '''圆周移位性质'''：时域序列的圆周移位，对应频域乘以一个线性相位因子（幅度谱不变，仅相位改变）。
# '''帕塞瓦尔定理 (Parseval's Theorem)'''：信号在时域的总能量等于其在频域的总能量，这保证了信号经过 DFT 变换后能量守恒。其数学表达式为：
<center><math>\sum_{n=0}^{N-1}|x[n]|^2 = \frac{1}{N}\sum_{k=0}^{N-1}|X[k]|^2</math></center>
# '''圆周卷积定理'''：这是 DFT 最强大的性质。两个序列在时域的'''圆周卷积'''，等于它们 DFT 的乘积。即：
<center><math>x_1[n] \otimes x_2[n] \leftrightarrow X_1[k] \cdot X_2[k]</math></center>
这一性质使得利用 FFT 加速长序列卷积运算成为可能。

== 4 关键技术指标与工程参数 ==

在利用 DFT 对真实信号进行谱分析时，以下参数直接决定了分析的精度与效果：

{| class="wikitable" style="width:100%"
! 参数名称 !! 符号/单位 !! 核心定义与工程意义
|-
! 频率分辨率
! <math>\Delta f</math> (Hz)
! DFT 能够分辨的两个相邻频率分量的最小间隔。由采样率 <math>f_s</math> 和采样点数 N 决定：<math>\Delta f = f_s / N</math>。要提高分辨率，必须增加信号的观测时间（即增加 N）。
|-
! 栅栏效应
! Picket Fence Effect
! DFT 只能观察到离散频率点 <math>k\Delta f</math> 处的频谱，就像透过栅栏看风景。如果信号的真实频率不在这些离散点上，其能量会泄漏到邻近的频点，导致测量误差。
|-
! 混叠现象
! Aliasing
! 如果采样率 <math>f_s</math> 不满足奈奎斯特采样定理（<math>f_s \ge 2f_{max}</math>），高频信号会伪装成低频信号出现在 DFT 频谱中。
|-
! 频谱泄漏
! Spectral Leakage
! 对信号的非整数周期截断会导致频谱扩散。在 DFT 运算前对信号加窗（如汉宁窗），可以有效抑制旁瓣，减小泄漏。
|}

== 5 DFT 与 FFT 的关系 ==

DFT 是数学定义，而'''快速傅里叶变换'''（FFT）是计算 DFT 的高效算法。
* '''直接计算 DFT'''：根据定义式计算 N 点 DFT，需要大约 N² 次复数乘法。当 N 很大时（如 N=4096），运算量极其庞大，无法满足实时处理需求。
* '''利用 FFT 计算'''：FFT 利用旋转因子的对称性和周期性，将运算量降低到 Nlog₂N 次复数乘法。
* '''工程现状'''：在现代工程实践（如 MATLAB 的 fft() 函数、Python 的 numpy.fft）中，我们口头上说的“做个 FFT”，实际上就是在计算数学上的 DFT。

== 6 典型应用与实战场景 ==

DFT（通过 FFT 实现）的应用几乎覆盖了所有涉及信号处理的领域：

{| class="wikitable" style="width:100%"
! 应用领域 !! 典型实例 !! 核心作用与原理
|-
! 频谱分析
! 音频均衡器 / 振动监测
! 通过 DFT 将麦克风或传感器采集的时域波形转换为频谱，直观展示信号中包含哪些频率成分及其能量大小。
|-
! 快速卷积
! 长序列滤波 / 图像处理
! 利用圆周卷积定理，将时域中复杂的卷积运算转化为频域中简单的乘法运算，极大提升了 FIR 滤波器和图像卷积核的处理速度。
|-
! 通信系统
! OFDM (4G/5G/Wi-Fi)
! 在正交频分复用技术中，利用 IDFT（IFFT）将并行的数据流调制到多个正交子载波上，接收端再用 DFT（FFT）进行解调。
|-
! 数据压缩
! MP3 / JPEG
! 利用 DFT（及其变种如离散余弦变换 DCT）将信号能量集中到少数几个频域系数上，舍弃人眼/人耳不敏感的高频分量，实现数据压缩。
|}

== 7 核心设计准则与常见误区 ==

# '''线性卷积与圆周卷积'''：DFT 对应的是圆周卷积。如果想用 DFT 计算两个有限长序列的'''线性卷积'''，必须对序列进行补零，使其长度 L ≥ N₁ + N₂ - 1，否则会产生时域混叠（重叠相加/保存法）。
# '''频率分辨率的真相'''：单纯在数据后面补零做更长的 DFT，只能让频谱曲线看起来更光滑（减小栅栏效应），但'''不能'''提高真实的物理频率分辨率。真实的分辨率只取决于信号的有效采样时长。
# '''实数序列的 DFT'''：实际采集的信号通常是实数。实数序列的 DFT 具有共轭对称性（即 X[k] = X*[N-k]），因此只需要计算前 N/2 个点即可完整描述信号频谱，后一半是冗余的。

== 8 学科发展与历史溯源 ==

* '''18世纪-19世纪'''：傅里叶提出连续傅里叶级数和变换理论，为频域分析奠定了基础。
* '''20世纪中期'''：随着数字计算机的诞生，科学家需要一种能让计算机处理频谱的数学工具，离散傅里叶变换（DFT）的理论体系逐渐成熟。
* '''1965年'''：库利（Cooley）和图基（Tukey）提出了快速傅里叶变换（FFT）算法，解决了 DFT 计算量过大的瓶颈，使得 DFT 从纯理论走向了广泛的工程应用，标志着现代数字信号处理时代的正式到来。

== 9 常见物理常数与参考 ==
* '''旋转因子'''：<math>W_N = e^{-j2\pi/N}</math>，DFT 运算的核心复数基底。
* '''频率分辨率公式'''：<math>\Delta f = f_s / N</math>，其中 <math>f_s</math> 为采样率，<math>N</math> 为采样点数。
* '''MATLAB/Python 实现'''：现代编程环境中，dft() 函数通常被高度优化的 fft() 函数取代，两者在数学结果上是等价的。

== 10 参见 ==
* [[快速傅里叶变换]] (FFT)
* [[数字信号处理]] (DSP)
* [[离散时间傅里叶变换]] (DTFT)
* [[采样定理]]
* [[圆周卷积]]

[[Category:应用数学]]
[[Category:信号处理]]
[[Category:算法]]
[[Category:电子工程]]