数字信号处理基础知识(pdf 25页)

第二章 数字信号处理基础

对采集到的信号进行处理，除了传统的时域分析之外，各种各样的变换发挥了重要作用，

从最熟悉的傅立叶（Fourier）变换，到现在的小波（wavelet）变换，以及主成分分析（principal

component analysis），独立成分分析（independent component analysis）和稀疏成分分析（sparse

component analysis）。每一种变换都有其独特的视野，为信号的分析处理提供了不同的思路。

这里我们将介绍最基本的一种变换，在线性时不变系统中广泛使用的傅立叶变换以及频谱分

析，本章还是以离散数据为主，介绍离散傅立叶变换（DFT）的有关知识。

第一节 傅立叶变换及其意义（Fourier Transform）

傅立叶分析方法的建立有过一段漫长的历史，涉及到很多人的工作和不同物理现象的研

究。在近代欧拉、伯努利、傅立叶、狄里赫利等学者的努力完善下，建立了傅立叶分析方法，

他们主要是集中在连续时间信号的分析问题上。与此同时，对于离散时间信号的傅立叶分析

方法却有着不同的发展过程，用于处理离散数据以产生数值近似的有关内插、积分和微分等

方面的公式早在17世纪的牛顿时代就被研究过，从事时间序列的研究曾吸引了18、19世纪

包括高斯在内的许多著名科学家，从而为离散傅立叶变换提供了数学基础。

在20世纪60年代中期，库利（Cooley）和图基（Tukey）独立发表了一篇论文，也就

是快速傅立叶变换算法（FFT）。FFT是非常高效的算法，使得计算变换所需要的时间减少

了几个数量级，由于计算机速度的迅速提高，越来越多的连续时间信号被离散化，然后用计

算机进行处理。

2.1.1 傅立叶变换的意义及各种变换对

利用“三角函数和”的概念来描述周期性过程至少可以追溯到古代巴比伦人时代，三角

函数和也即是成谐波关系的正弦和余弦或周期复指数函数的和。这些成谐波关系的复指数函

数在LTI系统分析中变得十分有用：如果一个LTI系统的输入可以表示为周期复指数的线性

组合，则输出也一定能表示成这种形式，并且输出线性组合中的加权系数与输入中对应的系

数有关，如图2.1所示，表示输入或者激励，表示系统输出或者响应，表)(nx)(ny)(jw

eH示系统单位脉冲响应的频率响应。 )(nh

)(nx)(ny

)(jw

eH

njw

e

0njwjw

eeH

0

0)()(nx)(ny

)(jw

eH

njw

kkkeA∑

njwjw

kkkkeeHA)(∑

图2.1 成谐波关系的复指数信号的响应

在研究LTI系统时，复指数信号的重要性就体现在图2.1中：一个LTI系统对复指数信号的

响应也是同样一个复指数信号，不同的只是乘了一个复振幅因子，频率并没有发生

变化，由于是复数因子，就有了幅度和相位或者实部和虚部的变化。 )(

kjw

eH

表2.1简要地综合了连续和离散时间信号的傅立叶级数和傅立叶变换表达式，有时候为

了统一，也把周期信号的傅立叶级数表示利用单位脉冲序列或单位冲激函数表示成傅立叶变

换。

表2.1 各种信号的傅立叶级数和傅立叶变换对

傅立叶级数 傅立叶变换

∑+∞

−∞==

kktjw

keatx

0)(~

时域是连续周期的 ∫+∞

∞−=dtejwXtxjwt

)(

21

)(

π

时域是连续非周期的

时域

频域 ∫−

=

00)(~1

0Tktjw

kdtetx

Ta

频域是离散非周期的 ∫+∞

∞−−

=dtetxjwXjwt

)()(

频域是连续非周期的

∑

==

)()/2(~

)(~

NknNjk

keanxπ

时域是离散周期的 ∫=

ππ

2)(

21

)(dweeXnxjwnjw

时域是离散非周期的

时域

频域 ∑

=−

=

)()/2(

)(~1

~

NnnNjk

kenx

Naπ

频域是离散周期的 ∑+∞

−∞=−

=

njwnjw

enxeX)()(

频域是连续周期的

从表2.1可以发现傅立叶变换在LTI系统分析中的思想，就是把一个无论多复杂的输入

信号分解成复指数信号的线性组合，那么系统的输出也能通过图2.1的关系表达成相同复指

数信号的线性组合，并且在输出中的每一个频率的复指数函数上乘以系统在那个频率的频率响应值。系数{

称为信号傅立叶级数系数或频谱系数或线谱等；称为的

频谱；也称为的频谱。作为线性组合所取的形式从求和过渡到积分，就是利用

傅立叶的思想，一个非周期信号可以看成是周期无限长的周期信号，当周期增加时，基频

越小，成谐波关系的各分量在频率上越来越近，当周期变得无穷大时，离散的线谱就形成了

一个连续谱，也就从求和变成了积分。 }

ka)(~

tx)(jwX)(tx

)(jw

eX)(nx

0w

从表2.1中时域和频域的关系还能得到如下规律：时域的离散必然导致频域的周期化，

频域的离散必然导致时域的周期化。简单的说，就是一个域离散必然另外一个域周期，相反

的，如果一个域连续必然另外一个域是非周期的。掌握了这个规律，我们很快就能判断出一

个信号在频域的表现形式。

2.1.2 离散傅立叶变换（DFT）

表2.1中给的四对傅立叶级数和傅立叶变换对在理论上有重要的意义，但在实际中往往

难以实现，尤其在数字计算机上实现是不太现实的，例如计算机无法处理连续的周期的信号。

因此我们需要的是一种在时域和频域都离散、非周期的一对傅立叶变换对，这就是离散傅立

叶变换，简称（DFT）。离散傅立叶变换的导出有多种方法，比较方便同时物理意义也比较

清晰的是从离散时间傅立叶变换（DTFT）和从离散傅立叶级数（DFS）入手。DTFT变换

为：，时域是离散非周期的，但频域是连续周期的，对连续变量

均匀采样，也就是对单位圆进行N等分，取一个周期的结果即得：∑+∞

−∞=−

=

njwnjw

enxeX)()(w

10,)()(1

02

2−≤≤=∑−

=−

=NkenxeXN

nnk

Nj

k

Nwjwπ

π，这样频谱变量由连续量变成了离散变量。

从DFS到DFT更加明显，DFS对应的时域和频域都是离散周期信号，可以在这两个域中分

别取它们的主值，也就是限定在一个周期内，这样就得到了DFT变换对。具体给出如下： w

10,)()]([)(1

02

−≤≤==∑−

=−

NkenxnxDFTkXN

nkn

Njπ

(2-1)

10,)(1

)]([)(1

02

−≤≤==∑−

=NnekX

NkXIDFTnxN

knk

Njπ

(2-2)

式（2－1）称为正变换，式（2－2）称为反变换。注意这一对变换对中信号的长)(nx度为N，它的频谱点长也为N，则和具有唯一的映射对应关系。也有可能

给一个7点的时间序列，求该信号的四点DFT或者10点DFT，明显前者求得的频谱不能真

实反映出信号，而后者不存在混叠现象，能反映出信号的频谱。因此在求DFT时要注意到

底是求多少点长的DFT，只有时域样点数小于或等于频域样点数，频谱才是真实的反映。 )(kX)(nx)(kX

【例2-1】试计算常用信号和)(nRN)()4

cos(nRn

NNπ

的N点DFT。

解：表示一个矩形窗，在滤波器设计中经常使用，我们在对一段未知信号进行

分析时的截取实际上就相当于乘了一个矩形窗，在实际分析信号时要注意矩形窗对信号频谱

的影响。 )(nRN

10),(

1,2,1,0

110,

)()(

22

1

02

1−≤≤=

⎪⎪

⎩⎪⎪

⎨⎧

−==

−−=

==

−−−

=−

∑

NkkN

Nk

eekN

enRkX

k

NjkN

NjN

nkn

Nj

Nδ

πππ

L

∑∑∑−

=+−−

=−−−

=−

+==1

0)2(2

1

0)2(2

1

02

2

21

21

)4

cos()(N

nkn

NjN

nkn

NjN

nkn

Nj

eeen

NkXπππ

π

与的求法类似，可得： )(

1kX

10))],2(()2([

221

)2(

2)(1

0)2(2

2−≤≤−−+−=+−=∑−

=−+−

NkNkkN

ekN

kXN

nNkn

Nj

δδδπ

为了表示方便，一般用符合来表示正交序列集中的基

NWNj

eπ2

−

，即＝

NWNj

eπ2

−

。因此

离散傅立叶变换对也可表示为：

10,)()]([)(1

0−≤≤==∑−

=NkWnxnxDFTkXN

nnk

N (2-3)

10,)(1

)]([)(1

0−≤≤==∑−

=−

NnWkX

NkXIDFTnxN

knk

N (2-4)

NW

具有下列性质：

周期性， rNn

Nn

NWW+

=

共轭对称性， ∗−

=)(n

Nn

NWW

可约性，或 n

rNrn

NWW

/=n

Nrn

rNWW=