小波分析教案第1章Word格式.docx
- 文档编号:14520253
- 上传时间:2022-10-23
- 格式:DOCX
- 页数:36
- 大小:812.80KB
小波分析教案第1章Word格式.docx
《小波分析教案第1章Word格式.docx》由会员分享,可在线阅读,更多相关《小波分析教案第1章Word格式.docx(36页珍藏版)》请在冰豆网上搜索。
1.1傅里叶变换
首先,让我们回顾一下传统的傅立叶变换
(1.1)
(1.2)
注意(1.1)式中对时间t的积分是在整个时间轴上,即为了获得信号中某一特定频率分量的信息,我们必须知道信号在整个时间过程中的变化情况。
也就是说,它在时域内是非局域的。
即使信号的持续时间是有限的,由于傅里叶变换的核函数在时域内是无限的,(1.1)式中也必须在信号的整个持续时间内积分,傅里叶分析也仍然是非局域的。
从上述分析可以看到,描述了信号的时域特征,而描述了信号的频域特征。
也就是说,要么在时域,要么在频域描述信号的特征。
但在许多实际问题中,我们关心的却是信号在局部范围内的特征:
例如音乐和语音信号中人们关心的是什么时刻演奏什么音符,发出什么样的音节;
对地震波的记录人们关心的是什么位置出现什么样的反射波;
图像处理中的边缘检测关心的是信号突变部分的位置,即纹理结构。
尤其对非平稳信号的处理中,信号在任一时刻附近的频率特性都很重要。
如柴油机缸盖表面的振动信号就是由撞击或冲击产生的,是一瞬变信号,仅从时域或频域上来分析都是不够的。
由于傅里叶分析完成上述任务存在其故有的困难,这就需要所谓的时间—频率分析方法,亦称为时频分析方法。
由于傅里叶变换在时域内是非局域的,不太适宜于暂态的、非平稳信号的分析。
针对暂态的、非平稳信号的分析,曾经出现过许多改进的办法,其中比较有成效的有Wigner—Ville分布和加窗傅里叶变换两种,前者是一种非线性二次型变换,与小波变换的概念相去甚远。
而D.Gabor于1946年提出的加窗傅里叶变换在非平稳信号的分析中起过很好的的作用,而且与小波变换有很多相似之处。
所以在学习小波变换之前,有必要了解一下加窗傅里叶变换。
1.2加窗傅里叶变换
传统的傅里叶变换是对信号在整个时间过程中变化情况的分析。
它在时域内是全局化的、非局域的。
其原因在于傅里叶变换的核函数正弦函数在时域内是无限的。
显然,将正弦函数乘以一个时域内衰减很快的函数之后作为核函数是一种合理的改进方法,这种衰减很快的函数我们称为窗函数,以表示。
一般来说,是一个实偶函数,且其能量主要集中在原点附近。
例如均值为零的高斯函数,就是一个常用的窗函数。
为了能实现对信号在整个时间过程中任意局部范围内变化情况的分析,应将窗函数应沿时间轴移动,得,其能量主要集中在附近。
对任意一个能量有限信号,其加窗傅里叶变换定义为
(1.3)
式中是一个实偶函数,譬如说高斯函数,它的能量主要集中在原点附近。
所以加窗傅里叶变换给出了信号在附近的频率信息。
由上式可以看到,加窗傅里叶变换确实是传统傅里叶变换的一种改进,即核函数由改进为
(1.4)
另一方面,我们也可以把加窗傅里叶变换理解为信号乘以窗函数后得到的时间函数的傅里叶变换。
由于窗函数的局域性,所以(1.3)式中的积分实质上只在附近的一段时间内进行。
也就是说,加窗傅里叶变换给出了信号在附近的一段时间内的频率信息。
现在我们可以看到,加窗傅里叶变换在时域内是局域化的。
也正是由于这个性质,加窗傅里叶变换又称为短时(Short-time)傅里叶变换。
根据希尔伯特空间中内积的定义,加窗傅里叶变换也可以写成内积形式:
(1.5)
应该注意,加窗傅里叶变换是频率和时间平移量的二元函数,所以信号的加窗傅里叶变换是对信号在时间—频率域内的联合分析。
而信号的傅里叶变换仅是频率的函数,所以传统的傅里叶变换仅是对信号在频率域内的分析。
为了更清楚地了解加窗傅里叶变换在时间—频率域的性质,我们有必要先分析其核函数(1.4)式的性质。
由于窗函数应选择为一平滑函数,所以其傅里叶变换为一低通滤波器。
那么
(1.6)
的傅里叶变换为
(1.7)
它显然是中心频率为的带通滤波器。
由(1.3)式可以得到在时间—频率平面内点的加窗傅里叶变换
(1.8)
其中横线表示取共轭。
由上式可以明显看出加窗傅里叶变换在时域内的局域性。
应用巴塞瓦定理,由上式可得:
(1.9)
由于仍是中心频率为的带通滤波器,所以可以解释为提取了信号在附近一段时间内在附近一段频率内的信息,也就是说,加窗傅里叶变换不仅在时域内是局域化的,而且在频域内也是局域化的。
或者说,加窗傅里叶变换是对信号在时间—频率域内进行联合分析。
为了更形象地说明这种联合分析,我们引入时—频窗的概念,如图1.1所示。
为了定量地说明时—频窗的性质,我们定义:
(1.10)
和分别称为和的标准差,他们描述了图1.1所示时—频窗的宽度。
对加窗傅里叶变换,其时—频窗的中心可以分别由(1.8)和(1.9)式中的变量和来调节。
但时—频窗的的宽度与窗的中心和无关,是固定不变的,如(1.10)所示。
为了提高加窗傅里叶变换的精度,我们自然希望时—频窗的宽度在两个方向都同时变窄。
但同时在时—频域内提高精度是有限制的,著名的测不准原理告诉我们:
(1.11)
当且仅当窗函数为高斯函数时,上式中的等号成立。
和量子力学中的测不准原理一样,我们不能无限地同时提高时—频率域内的测试精度。
这并不是由于测试仪器的限制,而是具有深刻的物理含义的。
讨论某一时刻的频率是多少是没有意义的。
从数学的角度看,加窗傅里叶变换是将任一空间中的函数映射为空间中的一个函数。
应用巴塞瓦定理可以证明:
(1.12)
上式说明,经过加窗傅里叶变换之后,范数保持不变(带有一固定的比例系数),故加窗傅里叶变换具有唯一的反演公式:
(1.13)
上式也称为加窗傅里叶变换的反变换。
应该注意,(1.12)隐含着如下意义:
经过加窗傅里叶变换之后,并未丢失信息,只不过是将时间信号所携带的信息在时间—频率域内用表示出来,也正是因为加窗傅里叶变换未丢失信息,所以可以通过反变换将信号完整地恢复出来。
1.3连续小波变换
加窗傅里叶变换实现了信号在时间—频率域内局域化的联合分析,但是其时间—频率窗的宽度却是固定不变的。
在加窗傅里叶变换中,我们当然可以选择不同类型的窗函数,或者是调整窗函数的参数,例如可以调整高斯函数的方差。
但是一旦窗函数的类型及参数确定之后,唯一可以调整的就是平移量,而变量只能改变窗沿时间轴的位置,不能改变时—频窗的宽度。
加窗傅里叶变换中的变量ω也仅仅是改变窗沿频率轴的位置,对时—频窗的宽度没有影响。
时—频窗的宽度是由窗函数的类型及其参数确定的。
这样加窗傅里叶变换就不能适应许多复杂的信号分析和处理。
上述这些性质都是由加窗傅里叶变换的核函数决定的。
下面我们将要详细讨论的连续小波变换也是信号在时间—频率域内局域化联合分析的一种方法。
首先我们将引入小波的概念,说明小波是一种持续时间很短的波。
在小波分析的所有文献著作中,均用表示小波,将按(1.16)式伸缩和平移而得到的一族函数称为分析小波。
以为核函数的积分变换,如(1.17)式所示,正是本节要讨论的小波变换。
需要特别强调的是,小波变换的时间—频率窗不是固定不变的,而是可以自适应调整的。
这是小波变换和加窗傅里叶变换的根本区别。
也正是由于这个特点,决定了小波分析在实际应用中的独特地位。
前面已提到,小波是一种持续时间很短的波。
但并不是任意持续时间很短的波都是小波。
小波必须满足以下允许条件。
设的傅里叶变换为,则允许条件表示为
(1.14)
满足上式的时间函数称为母小波或基小波,通常简称为小波。
要满足允许条件,的傅里叶变换在(直流)处的值必须为零,即=0。
由傅里叶变换的定义不难得到:
(1.15)
由上式可以清楚地看到,必须时正时负地波动,否则,的积分不会为零。
在实际应用中,我们要求时,速降为零。
也就是说,小波是持续时间很短的衰减振荡,它在时间域内是局域的。
正弦波也时正时负地波动,但是它是一种等幅的波动,在时间上是无限的,它显然不满足(1.15)式。
“小波”不是指其波动的幅度很小,而是持续时间很短。
在频率域,我们同样也要求时,速降为零,这也是允许条件(1.14)式对的约束。
再注意到,可见为一带通滤波器的频率特性。
这就是说,小波在频率域内也是局域的。
小波在时间域和频率域都是局域的,这是小波最重要的特性。
小波在时—频域内的局域性实际是其能量在时—频域内的集中性,主要是时域内集中于某一时刻附近。
将母小波伸缩和平移之后得到的函数族
(1.16)
称为分析小波。
式中a为伸缩参数或尺度参数,取正实数。
当a〉1时,沿时间轴方向拉伸。
当a〈1时,沿时间轴方向压缩。
因子是为保持伸缩之后能量不变,或者说范数不变,即。
b为平移参数,可以取任意实数。
实际上,加窗傅里叶变换中不同频率可理解为正弦函数疏密程度的不同,也可以理解为是沿时间轴的伸缩。
只不过加窗傅里叶变换中的伸缩和平移是分别针对正弦函数和窗函数的,而小波变换中伸缩和平移都是针对小波函数的。
图1.2中,中间波形为母小波,其对称中心位于0;
左面为a=2,b=-8时的;
右面为a=1/2,b=5时的。
图1.2分析小波
对任一能量有限信号,其连续小波变换定义为:
(1.17)
实际使用的小波函数中大部分都是实函数,这时上式中取共轭的符号也可以没有。
我们可以看到,连续小波变换和传统的傅里叶变换及加窗傅里叶变换其数学描述方法是类似的,都是取信号和核函数的积分。
这些变换也都可以解释为信号和核函数相关程度的度量。
而各种变换的区别在于选取的核函数不同。
从数学形式上看,傅里叶变换将一维时间函数映射为一维频率函数,所以傅里叶变换是对时间信号的频率分析。
加窗傅里叶变换将一维时间函数映射为二维函数,是对时间信号的时—频分析。
类似的,小波变换将一维时间函数映射为二维函数,也是对时间信号的时—频分析。
由于连续小波变换是按积分形式定义的,所以又称为积分小波变换。
以后如不特别声明,均认为小波是一个实函数。
由Parseval定理很容易得到
(1.18)
由(1.17)和(1.18)可以看出小波变换和加窗傅里叶变换一样,在时—频域内是同时局域化的。
我们要再一次强调,小波是一种持续时间很短的波。
这样(1.17)式中的积分区间主要由分析小波的持续时间决定,于是我们可以这样理解积分小波变换:
小波函数相当于一个时间域内的观测窗,窗的宽度和位置由参数a和b调整,我们正是在这一观测窗内对信号进行分析、提取我们某种感兴趣的信息。
也就是说,积分小波变换是在时间域内的一种局域化分析。
由于在时间观测窗内,小波函数是衰减振荡的,可以预料到积分小波变换是在此时间窗内提取信号的频率信息。
事实确实如此。
前面已指出,为一带通滤波器,假定带通滤波器的中心频率和带宽分别为和。
那么仍然是一个带通滤波器,只不过其中心频率移至/a,而其带宽相应地调整为/a。
那么(1.18)告诉我们,积分小波变换可以理解为信号通过一个带通滤波器,这就是说,连续小波变换在频率域内也是对信号的局域化分析。
到此为止,我们已经可以清楚地看到
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 分析 教案