离散傅里叶变换

title: 定义.
content: 对于formula_1点序列formula_2,它的离散傅里叶变换(DFT)为 其中formula_4是自然对数的底数,formula_5是虚数单位。通常以符号formula_6表示这一变换,即  离散傅里叶变换的逆变换(IDFT)为: 可以记为: 实际上,DFT和IDFT变换式中和式前面的归一化系数并不重要。在上面的定义中,DFT和IDFT前的系数分别为formula_10和formula_11。有时会将这两个系数都改成formula_12。
title: 从连续到离散.
content: 连续时间信号formula_13以及对应的连续傅里叶变换formula_14都是连续函数。由于数字系统只能处理有限长的离散信号,因此必须将formula_15和formula_16都离散化,并且建立对应的傅里叶变换。  假设formula_13时限于formula_18,再通过时域采样将formula_13离散化,就可以得到有限长离散信号,记为formula_20。设采样周期为formula_21,则时域采样点数formula_22。 它的傅里叶变换为 这就是formula_13在时域采样后的连续傅里叶变换,也就是离散时间傅里叶变换,它在频域依然是连续的。  下面将频域信号转化为有限长离散信号。与对时域信号的处理类似,假设频域信号是带限的,再经过离散化,即可得到有限长离散信号。依据采样定理,时域采样若要能完全重建原信号,频域信号formula_14应当带限于formula_27。由于时域信号时限于formula_18,由采样定理以及时频对偶的关系,频域的采样间隔应为formula_29。故,频域采样点数为: 即频域采样的点数和时域采样同为formula_1,频域采样点为formula_32。 而DTFT在频域上采样: 令formula_34,将其归一化,就得到前面定义的离散傅里叶变换。因此,DFT就是先将信号在时域离散化,求其连续傅里叶变换后,再在频域离散化的结果。
title: DFT与CFT.
content: 下面考察离散傅里叶变换与连续傅里叶变换(CFT,Continuous Fourier Transform)的关系。连续傅里叶变换 的采样为: 将这个积分以黎曼和的形式近似,有: 这一结论指出离散傅里叶变换确实是连续傅里叶变换在某种意义上的近似。不过必须注意以下两点: 为此,通常对连续信号的时域采样再做一次加窗处理(Windowing),这样就得到带限的有限长离散信号。
title: DFT与DTFT.
content: 离散时间傅里叶变换(DTFT)是在时域上对连续傅里叶变换的采样。DFT则是在频域上对DTFT的均匀采样。离散信号formula_38(n=0...,N-1)的DTFT为: 对formula_40在离散的频点formula_41上采样 即为formula_15的DFT。由于DTFT在频域是周期的,所以在DTFT频域上的均匀采样也应是周期的。formula_44实际上是这个周期序列的主值序列。
title: 栅栏效应.
content: formula_1点序列的DFT只能在有限的formula_1个频点上观察频谱,这相当于从栅栏的缝隙中观察景色,对于了解信号在整个频域上的特性是不够的。为了观察到其他频率的信息,需要对原信号formula_38做一些处理,以便在不同的频点上采样。  将原来在DTFT频域上的采样点数增加到formula_48点,这样采样点位置变为formula_49。则对应的DFT成为 若在序列formula_38之后补上M-N个零,设为formula_52,则上式变为 因此将formula_38补零再做DFT就可以得到formula_38的DTFT在其他频率上的值,相当于移动了栅栏,因而能够从其他位置进行观察。
title: 频谱分辨率.
content: formula_1点DFT的频谱分辨率是formula_57。栅栏效应一节指出可以通过补零观察到更多的频点,但是这并不意味着补零能够提高真正的频谱分辨率。这是因为formula_38实际上是formula_13采样的主值序列,而将formula_38补零得到的formula_52周期延拓之后与原来的序列并不相同,也不是formula_13的采样。因此formula_63与formula_16是不同离散信号的频谱。对于补零至formula_48点的formula_66的DFT,只能说它的分辨率formula_67仅具有计算上的意义,formula_63并不是真正的、物理意义上的频谱。频谱分辨率的提高只能在满足采样定理的条件下增加时域采样长度来实现。
title: 从空间的角度分析.
content: 周期为N的离散信号构成一个"N"维欧几里得空间formula_69。在这一空间上两个信号"x"和"y"的内积定义为 下面给出formula_69上的一组正交基: 将信号"x"在这组正交基上分解,得 令 此即为离散傅里叶变换。又 则有 此即为离散傅里叶变换的逆变换。  由此可知,在正交分解的角度上说,离散周期信号formula_15的离散傅里叶变换formula_16实质上是formula_15在正交基formula_80上的分量。而从线性变换的角度上说,formula_80是圆周卷积的特征向量,formula_16则是对应的特征值。
title: DFT与圆周卷积.
content: 根据卷积定理,离散信号x与y的圆周卷积对偶于频域上x与y离散傅里叶变换的乘积。下面揭示DFT与圆周卷积的内在关系。  对长为N的离散信号formula_83与formula_84,如果要计算它们的卷积,就必须定义formula_85与formula_86在"0≤n<N"之外的值。如果将formula_83与formula_84作周期为N的延拓,则有  如此,周期为N的圆周卷积:  卷积的结果仍然是以N为周期的离散信号。  前面指出,formula_91是DFT的特征向量,实际上它也是圆周卷积的特征向量。定义x与y的圆周卷积算子: 则formula_91与y的圆周卷积为  显然,y的DFT 就是圆周卷积的特征值。
title: 性质.
title: 完全性.
content: 离散傅里叶变换是可逆的线性变换 其中C表示复数集。即,任意"N"-维复向量都存在DFT和IDFT,而且其DFT和IDFT也是"N"-维复向量。
title: 正交性.
content: 向量组exp(2π"i kn/N")是"N"-维复数空间上的一组正交基: 其中δ是Kronecker delta。
title: 移位定理.
content: 时域信号序列formula_98的相位移动formula_99(其中formula_100为整数)在频域反映为“循环移位”,即:频域信号序列formula_101变为formula_102,其中下标是指将"k-m"对"N" 取余。类似的,由对偶性,时域信号序列的循环移位对应于频域信号序列的相移:
title: 周期性.
content: 上文中DFT与DTFT一节已经证明,离散序列的傅里叶变换是周期的。有限长序列formula_98的离散傅里叶变换formula_101可以被看作是它的周期延拓序列formula_108的离散时间傅里叶变换formula_109。由时频对偶性可知formula_101也可以被看作它的周期延拓的主值。  上一节的移位定理隐含着DFT的周期性,因为DFT的幅度formula_111不受输入信号的循环移位的影响,因为时移在频域对偶于相移,循环移位仅仅使DFT的相位产生平移。周期性的边界条件在DFT的许多应用中有重要作用。解差分方程时,就假设边界条件是满足周期性的,这是一个很有用的性质(参见应用)。
title: 普朗歇尔定理与帕塞瓦尔定理.
content: 如果"X"和"Y"分别是"x"和"y"的离散傅立叶变换,那么就有普朗歇尔定理:  其中星号表示复共扼。帕塞瓦尔定理是普朗歇尔定理的一个特例:
title: 应用.
content: DFT在诸多多领域中有着重要应用,下面仅是颉取的几个例子。需要指出的是,所有DFT的实际应用都依赖于计算离散傅里叶变换及其逆变换的快速算法,即快速傅里叶变换。
title: 快速傅里叶变换.
content: 快速傅里叶变换(即FFT)是计算离散傅里叶变换及其逆变换的快速算法。按照DFT的定义计算一个长为n的序列的DFT需要的计算复杂度达到了formula_114,而同样长度FFT的计算复杂度仅为formula_115。由于DFT的逆变换可以由DFT表示,所以DFT逆变换的计算同样可以由FFT完成。FFT算法的提出,使DFT得到了广泛的实际应用。
title: 频谱分析.
content: 前面指出,DFT是连续傅里叶变换的近似。因此可以对连续信号x(t)均匀采样并截断以得到有限长的离散序列formula_116,对这一序列作离散傅里叶变换,可以分析连续信号x(t)频谱的性质。前面还提到DFT应用于频谱分析需要注意的两个问题:即采样可能导致信号混叠和截断信号引起的频谱泄漏。可以通过选择适当的采样频率(见奈奎斯特频率)消减混叠。选择适当的序列长度并加窗可以抑制频谱泄漏。
title: 数据压缩.
content: 由于人类感官的分辨能力存在极限,因此很多有损压缩算法利用这一点将语音、音频、图像、视频等信号的高频部分除去。高频信号对应于信号的细节,滤除高频信号可以在人类感官可以接受的范围内获得很高的压缩比。这一去除高频分量的处理就是通过离散傅里叶变换完成的。将时域或空域的信号转换到频域,仅储存或传输较低频率上的系数,在解压缩端采用逆变换即可重建信号。
title: 解偏微分方程.
content: 离散傅里叶变换及其多维形式在偏微分方程的求解中也有应用。此时DFT被看作傅里叶级数的近似。傅里叶级数将函数在复指数"e"上展开,这正是微分算子的特征方程:"d"/"dx" "e" = "in" "e"。因此,通过傅里叶级数的形式,线性常微分方程被转换为代数方程,而后者是很容易求解的。此时得到的结果是偏微分方程解的级数表示,只要通过DFT逆变换即可得到其一般表示。这种方法被称作谱方法或级数解法。
title: 长整数与多项式乘法.
content: 目前长整数或多项式乘法最快速的算法是基于离散傅里叶变换的。由于整数(或多项式)乘法是逐位(或逐项)乘累加的形式,因此整数(或多项式)乘积的数字(或系数)可以用乘数数字(或乘式系数)的卷积表示。利用卷积定理,只要将数字(或系数)序列通过离散傅里叶变换变到频域,就可以将逐个乘累加的卷积变为对位的乘法,从而减少计算量,再以一次逆变换便可以得到乘法结果。需要注意整数乘法还有进位的问题。下面以多项式乘法为例说明这一应用:  假设需要计算多项式乘法"c"("x") = "a"("x")·"b"("x"),这一乘积结果的各项系数的表达式实际上是线性卷积的形式。由于离散傅里叶变换对应于圆周卷积,因此需要将这两个乘式的系数按照次数升序排列,序列后“补零”,使它们的长度"d"大于两式最高项次数的和:"d" > deg("a"("x")) + deg("b"("x"))。然后作圆周卷积: 其中c就是"c"("x")系数的向量。由于圆周卷积的DFT是乘积,有: 则 利用快速傅里叶变换,这一算法的运算复杂度为formula_120。
title: OFDM.
content: OFDM(正交频分复用)在宽带无线通信中有重要的应用。这种技术将带宽分为"N"个等间隔的子载波,可以证明这些子载波相互正交。尤其重要的是,OFDM调制可以由IDFT实现,而解调可以由DFT实现。OFDM还利用DFT的移位性质,在每个帧头部加上循环前缀(Cyclic Prefix),使得只要信道延时小于循环前缀的长度,就能消除信道延时对传输的影响。
title: 特殊例子.
title: 数论转换.
content: 数论转换是离散傅里叶转换(DFT)的一个特例 formula_121,此运算是根据模运算 而取得的,formula_122需要是质数,如此可以建构出有限域,并且存在 formula_123 个可以整除 formula_124 的实数根。如此可以得到formula_125,formula_126 为正整数。令 formula_127 为第 formula_124 个单位根,则第 formula_123 个单位根 formula_130 可以表示成 formula_131 。另外,再令formula_1 为 formula_130 次方 formula_122 模运算之循环个数。  则数论矩阵为  formula_135  formula_136 与 formula_137 各为 Column 与 Row 的指数(index), 令 formula_138 与 formula_139 各代表Column 与 Row 的总数,则formula_140, formula_141.  举个例子来说,当 formula_142 则  formula_143  formula_144  formula_145  formula_146  因为 formula_147 经 formula_148 的模运算为 formula_149,因此可以判定此情况为 formula_150 个次方一循环,所以 formula_151,formula_1可以整除 formula_124。则数论矩阵可以表示成  formula_154  在一些特殊的情况下,令 formula_155,而 formula_100 不是值数也是有意义的。像是 Fermat Number Transform 中formula_157,以及 Mersenne Number Transform 中 formula_158.

Summary: 离散傅里叶变换(英语:Discrete Fourier Transform,缩写为英语:DFT),是傅里叶变换在时域和频域上都呈离散的形式,将信号的时域采样变换为其DTFT的频域采样。

Transformation de Fourier discrète: La transformation de Fourier discrète (TFD), outil mathématique, sert à traiter un signal numérique. Elle constitue un équivalent discret de la transformation de Fourier (continue) utilisée pour traiter un signal analogique.