现代信号处理教程 - 胡广书(清华)

1

第1章 信号分析基础

1．1 信号的时－频联合分析

我们生活在一个信息社会里，而信息的载体就是我们本书要讨论的主题——信号。在我们身边以及在我们身上，信号是无处不在的。如我们随时可听到的语音信号，随时可看到的视频图像信号，伴随着我们生命始终的心电信号，脑电信号以及心音、脉搏、血压、呼吸等众多的生理信号。

对一个给定的信号，如)(t x ，我们可以用众多的方法来描述它，如)(t x 的函数表达式，

通过傅立叶变换所得到的)(t x 的频谱，即)(Ωj X ，再如)(t x 的相关函数，其能量谱或功率谱等。在这些众多的描述方法中，有两个最基本的物理量，即时间和频率。显然，时间和频率与我们的日常生活关系最为密切，我们时时可以感受到它们的存在。时间自不必说，对频率，如夕阳西下时多变的彩霞，音乐会上那优美动听的旋律以及在一片寂静中突然冒出的一声刺耳的尖叫等，这些都包含了丰富的频率内容。正因为如此，时间和频率也成了描述信号行为的两个最重要的物理量。 信号是变化着的，变化着的信号构成了我们周围五彩斑斓的世界。此处所说的“变化”，一是指信号的幅度随时间变化，二是指信号的频率内容随时间变化。幅度不变的信号是“直流”信号，而频率内容不变的信号是由单频率信号，或多频率信号所组成的信号，如正弦波、方波、三角波等。不论是“直流”信号还是正弦类信号都只携带着最简单的信息。

给定了信号)(t x 的函数表达式，或x 随t 变化的曲线，我们可以由此得出在任一时刻处

该信号的幅值。如果想要了解该信号的频率成分，即“在××Hz 处频率分量的大小”，则可通过傅立叶变换来实现，即

?∞

∞

-Ω-=Ωdt e t x j X t j )()(

（1.1.1a ） ?

∞

∞

-ΩΩΩ=

d e j X t x t j )()(21π

（1.1.1b ）

式中f π2=Ω，单位为弧度/秒，将)(Ωj X 表示成)

(|)(|ΩΩ?j e j X 的形式，即可得到

|)(|Ωj X 和)(Ω?随Ω变化的曲线，我们分别称之为)(t x 的幅频特性和相频特性。

如果我们想知道在某一个特定时间，如0t ，所对应的频率是多少，或对某一个特点的频

2

率，如0Ω，所对应的时间是多少，那么傅立叶变化则无能为力。

分析（1.1.1）式，对给定的某一个频率，如0Ω，那么，为求得该频率处的傅氏变换)(0Ωj X ，

（1.1.1a ）式对t 的积分仍需要从∞-到∞+，即需要整个)(t x 的“知识”。反之，如果我们要求出某一时刻，如0t 处的值)(0t x ，由（1.1.1b ）式，我们需要将)(Ωj X 对Ω从∞-至∞+作积分，同样也需要整个)(Ωj X 的“知识”

。实际上，由（1.1.1a ）所得到的傅氏变换)(Ωj X 是信号)(t x 在整个积分区间的时间范围内所具有的频率特征的平均表示。反之，（1.1.1b ）式也是如此，因此，傅立叶变换不具有时间和频率的“定位”功能。

前已述及，信号的幅度不但随时间变化，而且对现实物理世界中的大部分信号，其频率也随时间变化。实际上，在时域中愈是在较短时间内发生幅度突变的信号，其包含的信息就愈多。但由傅立叶变换)(Ωj X 看不出在什么时刻发生了此种类型的突变。现举两个例子说明这一概念。

例1．1．1 设信号x(n)由三个不同频率的正弦所组成，即

???

??=),sin(),sin(),sin()(321n n n n x ωωω 1

1102

211-≤≤-≤≤-≤≤N n N N n N N n

（1.1.2）

式中12312,ωωω>>>>N N N 。ω为圆周频率，s f f /2πω=，f 是信号的实际频率，s f 为抽样频率，所以ω的单位为弧度，Ω和ω的关系是[19] :

s s f f T /2πω=Ω=

（1.1.3）

)(n x 的波形如图1.1.1(a)所示，)(n x 的傅立叶变换的幅频特性|)(|ωj e X ，如图1.1.1(b)所示。

显然，|)(|ω

j e X 只给出了在21,ωω及3ω处有三个频率分量，给出了这三个频率分量的大小，但由此图看不出)(n x 在何时有频率1ω，何时又有2ω及3ω，即傅立叶变换无时间定位功能。

图1.1.1(c)是用我们后面所讨论的方法求出的)(n x 的联合时－频分布。该图是三维图形

的二维投影，在该图中，一个轴是时间，一个轴是频率。由该图可清楚地看出)(n x 的时间－频率关系。若将1.1.1(c)画成三维图，则如图1.1.1(d)所示。

3

例1.1.2 令 )exp()exp()(2n jn n j n x ωω== （1.1.4） 该信号称作线性频率调制信号，其频率与时间序号n 成正比，在雷达领域中，该信号又称作chirp 信号，图1.1.2(a)是其时域波形，127~0=n ，图1.1.2(b)是其频谱。显然，无论从时域波形还是从频域波形，我们都很难看出该信号的调制类型及其他特点。和图 1.1.1(c)一样，图1.1.2(c)也是)(n x 的时－频分布表示，由该图可明显看出，该信号的频率与时间成

图1.1.1 信号的时－频表示 （a ）信号x(n), (b) x(n)的频谱， (c) x(n)时－频分布的二维表示，(d) x(n)时－频分布的三维表示，

4

正比，且信号)(n x 的能量主要集中在时间－频率平面的这一斜线上。图1.1.2(d)是图1.1.2(c)的立体表示。

图1.1.2 chirp 信号的时－频表示. （a ）信号x(n), (b) x(n)的频谱，(c) x(n)时

－频分布的二维表示，(d) x(n)时－频分布的三维表示，

频率随时间变化的信号（如例1.1.2中的)(n x ）称为时变信号。文献[13]称这一类信号为“非平稳”信号，而把频率不随时间变化的信号称为“平稳”信号。此处的“平稳”和“不平稳”和随机信号中的“平稳随机信号”及“非平稳随机信号”的意义不同。平稳随机信号是指该

5

类信号的一阶及二阶统计特征（均值与方差）不随时间变化，其自相关函数和观察的起点无关，而非平稳信号的均值、方差及自相关函数均与时间有关，即是时变的。尽管这两类说法的出发点不同，但非平稳信号的频率实质上也是时变的，因此，把频率随时间变化的信号统称为“非平稳信号”并无大碍。但要说一个信号是“平稳信号”，则要具体说明所指的是频率不随时间变化的信号还是平稳随机信号。 由上述两例可以看出，傅立叶变换反映不出信号频率随时间变化的行为，因此，它只适合于平稳信号，而对频率随时间变化的非平稳信号，它只能给出一个总的平均效果。 现在，我们再从“分辨率”的角度来讨论傅立叶变换的不足。“分辨率”包含了信号的时域和频域两个方面，它是指对信号所能作出辨别的时域或频域的最小间隔（又称最小分辨细胞）。分辨能力的好坏一是取决于信号的特点，二是取决于所用的算法。对在时域具有瞬变的信号，我们希望时域的分辨率要好（即时域的观察间隔尽量短），以保证能观察到该瞬变信号发生的时刻及瞬变的形态。对在频域具有两个（或多个）靠得很近的谱峰的信号，我们希望频域的分辨率要好（即频域的观察间隔尽量短，短到小于两个谱峰的距离），以保证能观察这两个或多个谱峰。有关分辨率的讨论见文献[19]的第三章。 （1.1.1a ）式的傅立叶变换可以写成如下的内积形式：

><=

ΩΩt j e t x j X ),()(21

π

（1.1.5）

式中>

?

>=