现代信号处理教程 - 胡广书(清华)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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]称这一类信号为“非平稳”信号,而把频率不随时间变化的信号称为“平稳”信号。此处的“平稳”和“不平稳”和随机信号中的“平稳随机信号”及“非平稳随机信号”的意义不同。平稳随机信号是指该