离散傅里叶变换

(重定向自離散傅立葉轉換

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

在形式上,变换两端(时域和频域上)的序列都是的,而实际上这两组序列都应当被认为是离散周期信号的主值序列。即使对有限长的离散信号作DFT,也应当将其看作其周期延拓的变换。在实际应用中通常采用快速傅里叶变换计算DFT。

離散傅立葉轉換,是連續傅立葉轉換在離散樣本上的類比,目前廣泛應用於信號處理、數值分析、數位通信、音訊處理等領域,在快速演算法問世,以及硬體設備的提升後,能夠更廣泛的應用在人們的日常生活當中,也是個非常重要的學問之一。

定义

对于 点序列 ,它的离散傅里叶变换(DFT)为

 

其中 自然对数底数 虚数单位。通常以符号 表示这一变换,即

 

离散傅里叶变换的逆变换(IDFT)为:

 

可以记为:

 

实际上,DFT和IDFT变换式中和式前面的归一化系数并不重要。在上面的定义中,DFT和IDFT前的系数分别为  。有时会将这两个系数都改成 

从连续到离散

连续时间信号 以及对应的连续傅里叶变换 都是连续函数。由于数位系统只能处理有限长的离散信号,因此必须将  都离散化,并且建立对应的傅里叶变换。

假设 时限于 ,再通过时域采样将 离散化,就可以得到有限长离散信号,记为 。设采样周期为 ,则时域采样点数 

 

它的傅里叶变换为

 

这就是 在时域采样后的连续傅里叶变换,也就是离散时间傅里叶变换,它在频域依然是连续的。

下面将频域信号转化为有限长离散信号。与对时域信号的处理类似,假设频域信号是带限的,再经过离散化,即可得到有限长离散信号。依据采样定理,频域采样若要能完全重建原信号,频域信号 应当带限于 。由于时域信号时限于 ,由采样定理以及时频对偶的关系,频域的采样间隔应为 。故,频域采样点数为:

 

即频域采样的点数和时域采样同为 ,频域采样点为 。 而DTFT在频域上采样:

 

 ,将其归一化,就得到前面定义的离散傅里叶变换。因此,DFT就是先将信号在时域离散化,求其连续傅里叶变换后,再在频域离散化的结果。

DFT与CFT

下面考察离散傅里叶变换与连续傅里叶变换(CFT,Continuous Fourier Transform)的关系。连续傅里叶变换

 

的采样为:

 

将这个积分以黎曼和的形式近似,有:

 

这一结论指出离散傅里叶变换确实是连续傅里叶变换在某种意义上的近似。不过必须注意以下两点:

  • 时域采样必须满足采样定理
  • 离散傅里叶变换的处理对象是时限的

为此,通常对连续信号的时域采样再做一次加窗处理(Windowing),这样就得到带限的有限长离散信号。

DFT与DTFT

离散时间傅里叶变换(DTFT)是在时域上对连续傅里叶变换的采样。DFT则是在频域上对DTFT的均匀采样。离散信号 (n=0,...,N-1)的DTFT为:

 

 在离散的频点 上采样

 

即为 的DFT。由于DTFT在频域是周期的,所以在DTFT频域上的均匀采样也应是周期的。 实际上是这个周期序列的主值序列。

栅栏效应

 点序列的DFT只能在有限的 个频点上观察频谱,这相当于从栅栏的缝隙中观察景色,对于了解信号在整个频域上的特性是不够的。为了观察到其他频率的信息,需要对原信号 做一些处理,以便在不同的频点上采样。

将原来在DTFT频域上的采样点数增加到 点,这样采样点位置变为 。则对应的DFT成为

 

若在序列 之后补上M-N个零,设为 ,则上式变为

 

因此将 补零再做DFT就可以得到 的DTFT在其他频率上的值,相当于移动了栅栏,因而能够从其他位置进行观察。

频谱分辨率

 点DFT的频谱分辨率是 栅栏效应一节指出可以通过补零观察到更多的频点,但是这并不意味着补零能够提高真正的频谱分辨率。这是因为 实际上是 采样的主值序列,而将 补零得到的 周期延拓之后与原来的序列并不相同,也不是 的采样。因此  是不同离散信号的频谱。对于补零至 点的 的DFT,只能说它的分辨率 仅具有计算上的意义, 并不是真正的、物理意义上的频谱。频谱分辨率的提高只能在满足采样定理的条件下增加时域采样长度来实现。

从空间的角度分析

周期为N的离散信号构成一个N欧几里得空间 。在这一空间上两个信号xy内积定义为

 

下面给出 上的一组正交基

 

将信号x在这组正交基上分解,得

 

 

此即为离散傅里叶变换。又

 

则有

 

此即为离散傅里叶变换的逆变换。

由此可知,在正交分解的角度上说,离散周期信号 的离散傅里叶变换 实质上是 在正交基 上的分量。而从线性变换的角度上说, 圆周卷积特征向量 则是对应的特征值

DFT与圆周卷积

根据卷积定理,离散信号x与y的圆周卷积对偶于频域上x与y离散傅里叶变换的乘积。下面揭示DFT与圆周卷积的内在关系。

对长为N的离散信号  ,如果要计算它们的卷积,就必须定义  0≤n<N之外的值。如果将  作周期为N的延拓,则有

 

如此,周期为N的圆周卷积:

 

卷积的结果仍然是以N为周期的离散信号。

前面指出, 是DFT的特征向量,实际上它也是圆周卷积的特征向量。定义x与y的圆周卷积算子:

 

 与y的圆周卷积为

 

显然,y的DFT

 

就是圆周卷积的特征值。

性质

完全性

离散傅里叶变换是可逆的线性变换

 

其中C表示复数集。即,任意N-维复向量都存在DFT和IDFT,而且其DFT和IDFT也是N-维复向量。

正交性

向量组exp(2πi kn/N)是N-维复数空间上的一组正交基:

 

其中δknKronecker delta

移位定理

时域信号序列 的相位移动 (其中 为整数)在频域反映为“循环移位”,即:频域信号序列 变为 ,其中下标是指将k-mN 取余。类似的,由对偶性,时域信号序列的循环移位对应于频域信号序列的相移:

 
 
且有 

周期性

上文中DFT与DTFT一节已经证明,离散序列的傅里叶变换是周期的。有限长序列 的离散傅里叶变换 可以被看作是它的周期延拓序列 的离散时间傅里叶变换 。由时频对偶性可知 也可以被看作它的周期延拓的主值。

上一节的移位定理隐含着DFT的周期性,因为DFT的幅度 不受输入信号的循环移位的影响,因为时移在频域对偶于相移,循环移位仅仅使DFT的相位产生平移。周期性的边界条件在DFT的许多应用中有重要作用。解差分方程时,就假设边界条件是满足周期性的,这是一个很有用的性质(参见应用)。

普朗歇尔定理与帕塞瓦尔定理

如果XkYk分别是xnyn的离散傅立叶变换,那么就有普朗歇尔定理

 

其中星号表示复共扼。帕塞瓦尔定理是普朗歇尔定理的一个特例:

 

应用

DFT在诸多多领域中有着重要应用,下面仅是颉取的几个例子。需要指出的是,所有DFT的实际应用都依赖于计算离散傅里叶变换及其逆变换的快速算法,即快速傅里叶变换

快速傅里叶变换

簡介

快速傅立葉轉換(Fast Fourier Transform,簡稱FFT)是一種用來高效計算離散傅立葉轉換(Discrete Fourier Transform,DFT)及其反轉換的演算法。傅立葉轉換是一種數學變換,用來將信號從時間域或空間域轉換到頻率域。這在信號處理、圖像處理、數據分析等領域有著廣泛的應用。傳統的DFT計算複雜度較高,約為 ,而FFT將其降為 ,大大提高了計算效率。高效的計算能力使其在各種實際應用中扮演著關鍵角色。理解和應用FFT不僅能解決具體的技術問題,還能幫助我們更深入地了解信號處理的基本原理和方法。FFT的出現和應用,無疑是數學和工程領域的一大突破。

常見用於執行離散傅立葉轉換的演算法:

1. 庫利-圖基演算法(Cooley-Tukey algorithm)

庫利-圖基演算法是一種用於快速計算離散傅立葉轉換的方法,通過分治法策略(Divide-and-conquer)將一個長度為N的DFT分解為多個較小長度的DFT。這使得它特別適合處理長度為2的冪次方的序列,但實際上也適用於其他因數分解形式的序列。庫利和圖基在1965年首次提出這個方法,但類似的概念早在1805年由高斯提出。

2.基數-4、8、16、……演算法(Radix-4, 8, 16, …. Algorithms)

可以視為庫利-圖基演算法的一種推廣,使用不只是以2作為基底,以其他數字做為基底的快速演算法。

3.互質因子算法(prime factor algorithm)

互質因子算法,又稱Good–Thomas算法,是一種用於快速傅立葉變換(FFT)的算法。它將大小為 的DFT轉換為 的二維DFT,前提是  必須是互質的。這種算法在此特殊情況下比傳統的Cooley–Tukey算法更有效。

4.格策爾演算法(Goertzel Algorithm)

此演算法在1958年被傑拉德 · 格策爾所提出。

此算法,使用數字濾波方法計算離散傅立葉轉換(DFT)係數和信號頻譜。修改後的格策爾演算法可以用來計算信號頻譜,而無需像DFT算法那樣涉及複雜的代數計算。

5. 啁啾-Z轉換(Chirp Z transform (CZT))

啁啾-Z轉換(Chirp-Z transform)是一種離散傅立葉轉換(DFT)的擴展算法,適用於當取樣頻率間隔與取樣時間間隔的乘積不等於信號的時頻分佈面積時的情況。它利用卷積來實現任意大小的DFT,從而加快計算速度,是一種基於快速傅立葉轉換(FFT)的高效算法。

此外,啁啾-Z變換也是一種用於評估非單位圓上的Z變換的算法。

6.威諾格拉德演算法(Winograd algorithm)

由以色列裔美國計算機科學家Shmuel Winograd於1976年提出的,是一種改進型的快速傅立葉轉換(FFT)算法,旨在進一步提高計算效率和速度。此演算法使用的加法數量與庫利-圖基演算法大致相同,但其所需的乘法數量只有其演算法的約20%。

频谱分析

前面指出,DFT是连续傅里叶变换的近似。因此可以对连续信号x(t)均匀采样并截断以得到有限长的离散序列 ,对这一序列作离散傅里叶变换,可以分析连续信号x(t)频谱的性质。前面还提到DFT应用于频谱分析需要注意的两个问题:即采样可能导致信号混叠和截断信号引起的频谱泄漏。可以通过选择适当的采样频率(见奈奎斯特频率)消减混叠。选择适当的序列长度并加窗可以抑制频谱泄漏。

数据压缩

由于人类感官的分辨能力存在极限,因此很多有损压缩算法利用这一点将语音、音频、图像、视频等信号的高频部分除去。高频信号对应于信号的细节,滤除高频信号可以在人类感官可以接受的范围内获得很高的压缩比。这一去除高频分量的处理就是通过离散傅里叶变换完成的。将时域或空域的信号转换到频域,仅储存或传输较低频率上的系数,在解压缩端采用逆变换即可重建信号。

解偏微分方程

离散傅里叶变换及其多维形式在偏微分方程的求解中也有应用。此时DFT被看作傅里叶级数的近似。傅里叶级数将函数在复指数einx上展开,这正是微分算子的特征方程:d/dx einx = in einx。因此,通过傅里叶级数的形式,线性常微分方程被转换为代数方程,而后者是很容易求解的。此时得到的结果是偏微分方程解的级数表示,只要通过DFT逆变换即可得到其一般表示。这种方法被称作谱方法或级数解法。

长整数与多项式乘法

目前长整数多项式乘法最快速的算法是基于离散傅里叶变换的。由于整数(或多项式)乘法是逐位(或逐项)乘累加的形式,因此整数(或多项式)乘积的数字(或系数)可以用乘数数字(或乘式系数)的卷积表示。利用卷积定理,只要将数字(或系数)序列通过离散傅里叶变换变到频域,就可以将逐个乘累加的卷积变为对位的乘法,从而减少计算量,再以一次逆变换便可以得到乘法结果。需要注意整数乘法还有进位的问题。下面以多项式乘法为例说明这一应用:

假设需要计算多项式乘法c(x) = a(xbx),这一乘积结果的各项系数的表达式实际上是线性卷积的形式。由于离散傅里叶变换对应于圆周卷积,因此需要将这两个乘式的系数按照次数升序排列,序列后“补零”,使它们的长度d大于两式最高项次数的和:d > deg(a(x)) + deg(b(x))。然后作圆周卷积:

 

其中c就是c(x)系数的向量。由于圆周卷积的DFT是乘积,有:

 

 

利用快速傅里叶变换,这一算法的运算复杂度为 

OFDM

OFDM(正交频分复用)在宽带无线通信中有重要的应用。这种技术将带宽分为N个等间隔的子载波,可以证明这些子载波相互正交。尤其重要的是,OFDM调制可以由IDFT实现,而解调可以由DFT实现。OFDM还利用DFT的移位性质,在每个帧头部加上循环前缀(Cyclic Prefix),使得只要信道延时小于循环前缀的长度,就能消除信道延时对传输的影响。

特殊例子

数论转换

數論轉換是一種可以計算摺積的快速演算法。計算摺積的快速演算法中最常用的一種是使用快速傅里葉變換,然而快速傅立葉變換具有一些實現上的缺點,舉例來說,資料向量必須乘上複數係數的矩陣加以處理,而且每個複數係數的實部和虛部是一個正弦餘弦函數,因此大部分的係數都是浮點數,也就是說,必須做複數而且是浮點數的運算,因此計算量會比較大,而且浮點數運算產生的誤差會比較大。

而在數論轉換中,資料向量需要乘上的矩陣是一個整數的矩陣,因此不需要作浮點數的乘法運算,更進一步,在模數為的情況下,只有種可能的加法與種可能的乘法,這可以使用記憶體把這些有限數目的加法和乘法存起來,利用查表法來計算,使得數論轉換完全不須任何乘法與加法運算,不過需要較大的記憶體與消耗較大的存取記憶體量。

雖然使用數論轉換可以降低計算摺積的複雜度,不過由於只進行整數的運算,所以只能用於對整數的訊號計算摺積,而利用快速傅立葉變換可以對任何複數的訊號計算摺積,這是降低運算複雜度所要付出的代價,也不一定所有信號都可以找到適合做數論轉換的質數。

 ,此運算是根據模運算 而取得的, 需要是質數 ,如此可以建構出有限域,並且存在   個可以整除  實數根。如此可以得到   為正整數。令   為第  單位根,則第   個單位根   可以表示成   。另外,再令   次方   模運算之循環個數。

數論轉換的公式

 

 

其中有幾點限制以及應該注意的事項:

1.M應是一個質數

2.N應該是M-1的因數

3. 是滿足    modM=1 的整數(當 N = M−1 時,   = M−1)。它也被稱為 N 的模 M 的逆元

4. 是 N 階的單位根。


舉個例子來說,當  

 

 

 

 

因為    的模運算為   ,因此可以判定此情況為   個次方一循環,所以   可以整除  。則數論矩陣可以表示成

 

在一些特殊的情況下,令   ,而   不是值數也是有意義的。像是 Fermat Number Transform ,以及 Mersenne Number Transform .

數論轉換的數學性質

前面提到數論轉換是一種離散傅立葉轉換的特例,因此數論轉換具有許多離散傅立葉轉換的數學性質。

若F[k]為f[n]的數論轉換結果,即f[n]→F[k](NTT)

1.正交性質:

矩陣中的每個行或列,彼此正交,即任兩個行或列的內積為0

 

證明:

  

  

2.線性性質:

NTT 符合線性轉換,即:

 

並且對於常數 c 有:

 

 ,則 

3.平移性質:

 ,則 ,且 

4.可逆性質:

NTT 變換是可逆的,存在反轉換 INTT 使得:

INTT(NTT(n))=n

5.對稱性質:

 具有對稱性質,即 ,則 的數論轉換 也會有對稱性質,即 的特性。

 ,則 的數論轉換 也會有 的特性。

6.摺積性質:

此性質在線性非時變系統當中非常重要,NTT的主要應用多與此性質有關。

若考慮兩個信號,  

 

 

數論轉換的優缺點

優點:

1.使用整數進行計算,計算較容易,也可以節省記憶體空間。

2.使用整數不會像浮點數一樣產生誤差,可以非常精確地計算。

3.很多性質都跟離散傅立葉轉換相同。

缺點:

1.只能對整數做計算,若信號為非整數,可能需要事先處理。

2.不容易找到root of unity,即不容易找到合適的對數轉換。

3.相對而言缺乏物理意義。

参阅

参考文献

  1. Oppenheim, Alan V.; Schafer, R. W.; and Buck, J. R., (1999). Discrete-time signal processing, Upper Saddle River, N.J. : Prentice Hall. ISBN 0137549202
  2. Mallat, S., A Wavelet Tour of Signal Processing, Second Edition, Academic Press, ISBN 0-12-466606-x
  3. Gill, J., Fourier Transform and Its Applications[1]页面存档备份,存于互联网档案馆))
  4. Jian-Jiun Ding, Advanced Digital Signal Processing class note,the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2024
  5. Duhamel, Pierre, and Martin Vetterli. "Fast Fourier transforms: a tutorial review and a state of the art." Signal processing 19.4 (1990): 259-299.