快速傅里叶变换(1).pdf

2.3快速傅立叶变换问题

1)问题背景

在数值电路的传输中，为了避免信号干扰，需要把一个连续信号x(t)先通

过取样离散化为一列数值脉冲信号x(0),x(1),……,然后再通过编码送到传

输电路中。如果取样间隔很小，而连续信号的时间段又很长，则所得到的数值脉

冲序列将非常庞大。因此，传输这个编码信号就需要长时间的占用传输电路，相

应地也需要付出昂贵的电路费用。

那么能否经过适当处理是使上述的数值脉冲序列变短，而同时又不会丧失有

用的信息？的经过研究，人们发现，如果对上述数值脉冲序列作如下的变换处理：

N1

X(n)x(k)e2nki/N,n0,1,...,N1,i1

k0(1)

则所得到的新序列X(0),X(1),……将非常有序，其值比较大的点往往集中在

某一很狭窄的序列段内，这将非常有利于编码和存储，从而达到压缩信息的目的。

公式（1）就是所谓的离散傅立叶变换，简称DFT。现在我们来分析一下计

算DFT所需要的工作量。如果我们不考虑公式（7.1）中指数项的运算，那么计

算其每一个点X(n)需要N次复数乘法和N-1次的复数加法。显然当N很大时，

这个工作量也非常巨大。正是由于这个原因，使得DFT的应用范围在过去很长的

时间里受到了严格的限制。注意到公式（1）是非常有规律性的，那么能否利用

这种规律性来降低DFT的计算时间？

1965年，凯莱和塔柯的提出了一种用于计算DFT的数学方法，大大减少了

DFT的计算时间，同时又特别适用于硬件处理，这就是所谓的快速傅里叶变换，

简称FFT。鉴于DFT的数据结构可以通过傅立叶变换的离散化获得，亦可通过三

角插值得到，而本质上又同连续傅里叶分析有着极为密切的关系。下面我们从傅

立叶级数级数和傅立叶积分入手，导出DFT结构的来源和FFT的工作原理。

2）傅立叶变换

如果x(t)是定义在整个实轴上的实值或复值函数，则其傅立叶变换可由下式

给出：



X(f)x(t)e2nift/Tdt,i1（2）

若对任意参数f，上述积分都存在，则（2）式确定了一个函数X(f)，称为x(t)

的傅立叶变换。如果已知X(f)则利用如下的傅立叶逆变换，还可复原x(t)：



x(t)X(f)e2nift/Tdf,i1（3）

若x(t)和X(f)同时满足（2）、（3）式，则称他们是一个傅立叶变换对，记

为x(t)X(f)。。通常X(f)是一个复函数，因此可以写成如下两部分：

X(f)R(f)I(f)i，i1（4）

式子中R(f)，I(f)分别是X(f)的实部和虚部。将上式表示为指数形式：

X(f)X(f)ei(f)i1（5）

其中

22

X(f)R(f)I(f)

（6）

1I(f)

(f)tg()

R(f)

工程技术中，常将x(t)看成一时间信号，相应的空间，称为时间域和空域；

将其傅立叶变换X(f)看成频率函数，相应的空间称为频域。X(f)称为x(t)的

傅立叶谱，而(f)称为其相角，这在物理上是有良好背景的。的譬如此频率的

的含义可以这样来理解：应用欧拉公式可将指数项表示成正弦-余弦的形式，如

果把（2）式解释成离散项和的极限，则显然X(f)是包含了无限项正弦-余弦的