第七章2连续小波变换ppt课件_第1页
第七章2连续小波变换ppt课件_第2页
第七章2连续小波变换ppt课件_第3页
第七章2连续小波变换ppt课件_第4页
第七章2连续小波变换ppt课件_第5页
已阅读5页,还剩26页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、2 连续小波变换连续小波变换基本小波基本小波连续小波变换的定义连续小波变换的定义连续小波变换的性质连续小波变换的性质常用的基本小波常用的基本小波连续小波变换的逆变换连续小波变换的逆变换连续小波变换的再生核连续小波变换的再生核小波时频分析小波时频分析CWT的变换过程示例的变换过程示例连续小波变换的数值积分结果演示连续小波变换的数值积分结果演示连续小波变换的应用连续小波变换的应用基本小波基本小波定义:定义:2|( )|Cd| 则称)(t为基本小波,或母小波。2( )( )tL R,其Fourier变换为 ()满足设允许条件性质:性质:(1) 带通性质(2) 零均值和波动性(3) 局部化特性t)ab

2、-t(t)|a |b) )(a,(21 -dffW尺度因子平移因子归一化因子a,ba,b (W f )(a, b)f (t)(t)dtf (t), (t)连续小波变换的定义连续小波变换的定义定义:定义: 对任意,它的连续小波变换定义为:2( )( )f tL R-1/2a,bt-b(t)|a|(): aR, a0; bR a称为小波基函数一般可简记为连续小波变换的性质连续小波变换的性质 (4)自相似性:对应不同尺度参数a和不同的平移参数b的连续小波变换之间是自相似的 线性性: 一个多分量信号的连续小波变换等于各个分量 的小波变换之和(2)平移不变性:假设的小波变换为( )f t( , )fW

3、ab()f t,那么( ,)fW ab的小波变换为(3)伸缩共变性:假设的小波变换为( )f t的( , )fW ab,那么( )f ct1( , ),0fW cacbcc小波变换为(5)冗余性: 连续小波变换中存在信息表述的冗余连续小波变换的逆变换连续小波变换的逆变换 逆变换存在的可能性: 窗口宽度任意调节,在时域上或频域上能完全恢复出信号的信息。 连续小波变换结果有很大的冗余度。 以a,b(t)为变换核的连续小波变换的逆变换: 假设 f(t)、(t) L2(R):(t) (t),t(t)(t)b) )(a,( ba,ba,fdffWbaa(t)b) )(a,(C(t)2ba,1 -ddfW

4、f 2|( )| Cd|其中:常用的基本小波常用的基本小波 Haar小波小波101/2( )11/210ttt 其它/224( )sin/4iie 2. Daubechies小波小波D4尺度函数与小波尺度函数与小波 012345-0.4-0.200.20.40.60.811.21.4-2-10123-1.5-1-0.500.511.52D6尺度函数与小波尺度函数与小波 常用的基本小波常用的基本小波 3、双正交小波、双正交小波双正交双正交B样条小波样条小波5-3)、)、 (9-7小波滤波器小波滤波器bior2.2, bior4.4(7-5小波滤波器小波滤波器: 11 1 3 1 1,8 2 4

5、2 82h2,2nnnnphqh2022221222222232021438245182411644281 212qpqqqpqqpqqqpqqqq 13355533,164 16 2 164 162h常用于图形学中。其中尺度函数是一常用于图形学中。其中尺度函数是一个三次个三次B样条。样条。常用的基本小波常用的基本小波 4. Morlet小波小波20/2( )itttee20() /2( )2 e Morlet小波不存在尺度函数小波不存在尺度函数; 快速衰减但非紧支撑快速衰减但非紧支撑. Morlet小波是Gabor 小波的特例。 2221/421ti tg tetg t e Gabor 小波

6、Morlet小波1,5常用的基本小波常用的基本小波 5. 高斯小波高斯小波 2/212ttte 2/2i e ( ) t( ) 这是高斯函数的一阶导数,在信号与图象的边缘提取中具有重要的应用。这是高斯函数的一阶导数,在信号与图象的边缘提取中具有重要的应用。主要应用于阶梯型边界的提取。主要应用于阶梯型边界的提取。 特性:特性: 指数级衰减,非紧支撑;具有非常好的时间频率局部化;指数级衰减,非紧支撑;具有非常好的时间频率局部化; 关于关于0轴反对称。轴反对称。常用的基本小波常用的基本小波 6. Marr小波也叫墨西哥草帽小波)小波也叫墨西哥草帽小波) 这是高斯函数的二阶导数,在信号与图象的边缘提取

7、中具有重要的应用。这是高斯函数的二阶导数,在信号与图象的边缘提取中具有重要的应用。主要应用于屋脊型边界和主要应用于屋脊型边界和Dirac边缘的提取。边缘的提取。 22/22( )(1)3ttte242/22 2( )3e 特性:特性: 指数级衰减,非紧支撑;具有非常好的时间频率局部化;指数级衰减,非紧支撑;具有非常好的时间频率局部化; 关于关于0轴对称。轴对称。( ) t常用的基本小波常用的基本小波 7. Meyer小波小波它的小波函数与尺度函数都是在频域中进行定义的。具体定义如下: 122324sin1 22333481 2433280 ,332cosivve 42335847020 0,1

8、v tttttt 121222 33242cos1 223340 3v t( ) 常用的基本小波常用的基本小波 8. Shannon小波小波 sin1/2sin21/21/2tttt /21, 20, ie 其它在时域,在时域,Shannon小波是无限次可微的,具有无穷阶消失矩,不是小波是无限次可微的,具有无穷阶消失矩,不是紧支的,具有渐近衰减性但较缓慢;在频域,紧支的,具有渐近衰减性但较缓慢;在频域,Shannon小波是频率小波是频率带限函数,具有好的局部化特性。带限函数,具有好的局部化特性。 t常用的基本小波常用的基本小波 9. Battle-Lemarie样条小波样条小波 2242224

9、12sin164sin241sin3 8sin8sin34441( )() ()222 ige t常用的基本小波常用的基本小波 连续小波变换的再生核连续小波变换的再生核 尺度和位移连续变化的连续小波基函数构成了一组非正交的过度完全基,小波展开系数之间的相关关系可表示为,1( , ; , )( )( )aaRKaatt dtC 1. CWT系数具有很大的冗余,计算量比较大2. 利用冗余性可以实现去噪和数据恢复的目的。重建核方程重建核方程daaKaWadaaWff),;,(),(),(000200重建核小波时频分析小波时频分析 小波分析能提供随频率改变的时小波分析能提供随频率改变的时-频窗频窗(自

10、适应窗口自适应窗口)。 假设 是任一基本小波,并且 与都是窗函数,它们的中心 与半径分别为 ,t和。不妨设和尺度 a都是正数。1( , )( )()b atafb atatbWTa bf tdtaa ,1( , ),( )()22ibfa baWTa bffead给出了 ft在频域窗 ( , )fWTa b内的局部化信息。 ,aaaa给出了 ft在时间窗 ( , )fWTa b内的局部化信息。 ,batabata 小波时频分析小波时频分析 内的局部化信息, /a( , )fWTa b若用作为频率变量,那么,那么给出了信号t在时间频率平面(平面中一个矩形的时间频率窗 ,batabataaaaa

11、即小波变换具有时频局部化特征。 120aa2a窗宽窗宽:面积面积:,a b的宽度是宽度的a倍。对信号 f t的高频成分需用尺度较小的分析小波。信号进行分析。 这时时间窗会自动变窄,并在高频区域对4 “恒Q性质”: 小波函数a,b(t)的频率具有带通特性,带通滤波特性的品质因数Q为*/*2/2aQa小波时频分析小波时频分析 各种变换的比较各种变换的比较 小波变换的特性小波变换的特性分解种类:时间分解种类:时间-尺度或时间尺度或时间-频率频率分析函数:具有固定震荡次数的时间有限的波。分析函数:具有固定震荡次数的时间有限的波。 小波函数的伸缩改变其窗口大小。小波函数的伸缩改变其窗口大小。变量:变量:

12、 尺度,小波的位置尺度,小波的位置 信息:窄的小波提供好的时间局部化及差的频率信息:窄的小波提供好的时间局部化及差的频率 局部化,宽的小波提供好的频率局部化局部化,宽的小波提供好的频率局部化 及差的时间局部化。及差的时间局部化。适应场合:非平稳信号适应场合:非平稳信号短时短时Fourier变换的特性变换的特性 分解种类:时间分解种类:时间-频率频率 分析函数:由三角震荡函数复合而成的时间有限的波分析函数:由三角震荡函数复合而成的时间有限的波 变量:频率,窗口的位置变量:频率,窗口的位置 信息:信息: 窗口越小,时间局部化越好,其结果是滤掉低频成分;窗口越小,时间局部化越好,其结果是滤掉低频成分

13、; 窗口越大,频率局部化越好窗口越大,频率局部化越好, 此时时间局部化较差此时时间局部化较差. 适应场合:次稳定信号适应场合:次稳定信号 Fourier变换变换的特性变换变换的特性分解种类:频率分解种类:频率分析函数:正弦函数,余弦函数分析函数:正弦函数,余弦函数 变量:变量: 频率频率信息:组成信号的频率信息:组成信号的频率 适应场合:平稳信号适应场合:平稳信号CWT的变换过程示例的变换过程示例CWT的变换过程示例,见图,可分如下5步小波 (t)和原始信号f(t)的开始部分进行比较 计算系数C该部分信号与小波的近似程度;C值越高表示信号与小波相似程度越高小波右移k得到的小波函数为 (t-k)

14、 ,然后重复步骤1和2,直到信号结束 扩展小波,如扩展一倍,得到的小波函数为 (t/2) 重复步骤14 图2-1 连续小波变换的过程 对于含噪声的周期性信号noissin进行连续小波分解,首先使用db4小波对比小波对信号进行不同尺度的分解,重新选择分解尺度后脚明显地显示出了信号的周期性特征,最后使用复小波cgau4对信号进行分解,分别得到信号的实部、虚部、模和相位。结果如图2-2、2-3和图2-4所示。连续小波变换的数值积分结果演示连续小波变换的数值积分结果演示01002003004005006007008009001000-1.5-1-0.500.511.5连续小波变换的数值积分结果演示连续

15、小波变换的数值积分结果演示图2-2 原始信号连续小波变换的数值积分结果演示连续小波变换的数值积分结果演示Absolute Values of Ca,b Coefficients for a = 1 2 3 4 5 .time (or space) bscales a2004006008001000 1 4 710131619222528313437404346Absolute Values of Ca,b Coefficients for a = 2 4 6 8 10 .time (or space) bscales a2004006008001000 2 10 18 26 34 42 50

16、58 66 74 82 90 98106114122图2-3 不同尺度下的连续小波分解连续小波变换的数值积分结果演示连续小波变换的数值积分结果演示Real part of Ca,b for a = 2 4 6 8 10 .time (or space) bscales a2004006008001000 2 10 18 26 34 42 50 58 66 74 82 90 98106114122Imaginary part of Ca,b for a = 2 4 6 8 10 .time (or space) bscales a2004006008001000 2 10 18 26 34 42

17、 50 58 66 74 82 90 98106114122Modulus of Ca,b for a = 2 4 6 8 10 .time (or space) bscales a2004006008001000 2 10 18 26 34 42 50 58 66 74 82 90 98106114122Angle of Ca,b for a = 2 4 6 8 10 .time (or space) bscales a2004006008001000 2 10 18 26 34 42 50 58 66 74 82 90 98106114122图2-4 信号的复连续小波分解 正弦信号在位置正

18、弦信号在位置500处有一间断点,由于间断点持续的处有一间断点,由于间断点持续的时间很短,因此在波形图中无法辨识出间断点来。如图时间很短,因此在波形图中无法辨识出间断点来。如图2-5。连续小波变换用于断点分析连续小波变换用于断点分析020040060080010001200-101正 弦 信 号 x102004006008001000120000.050.1间 断 信 号 x2020040060080010001200-101含 间 断 点 的 正 弦 信 号 x图2-5 原始信号 对该信号进行快速对该信号进行快速Fourier变换,得到信号的频谱图。观察变换,得到信号的频谱图。观察信号的频谱,

19、除了正弦信号的基频之外,无任何额外的信息。信号的频谱,除了正弦信号的基频之外,无任何额外的信息。如图如图2-6。 连续小波变换用于断点分析连续小波变换用于断点分析0501001502002503003504004505000100200300信 号 x1的 频 谱 图0501001502002503003504004505000100200300信 号 x的 频 谱 图图2-6 频谱图 选取选取db4小波,尺度小波,尺度a=1,2,3,4,5,6,7,8,9,10,对信号实施连续小波变换,得到小波系数分布图。在空间位对信号实施连续小波变换,得到小波系数分布图。在空间位置置500处间断点的信息一

20、目了然地显现出来。在尺度为处间断点的信息一目了然地显现出来。在尺度为1处,处,在整个时间范围空间位置内分布的小波系数是正弦信号在整个时间范围空间位置内分布的小波系数是正弦信号的基频信息。而空间位置的基频信息。而空间位置500处的间断点在时间上是局部的,处的间断点在时间上是局部的,小波系数集中发生在间断位置的各个尺度上,将间断信息凸小波系数集中发生在间断位置的各个尺度上,将间断信息凸现了出来。如图现了出来。如图2-7。连续小波变换用于断点分析连续小波变换用于断点分析 传统的快速Fourier变换对此却无能为力,究其原因在于间断分布的能量在整个信号中所占的比例太小,而Fourier变换并无时间局部化能力。连续小波变换用于断点分析连续小波变换用于断点分析信 号 x的 连 续 小 波 变 换time (or space) bscales a1002003004005006007008009001000 1 2 3 4 5 6 7 8 910图2-7 连续小波变换连续小波变换用于自相似性检测连续小波变换用于自相似性检测 直观上讲,小波分解可通过计算信号和小波之间的自相似指数(小波系数)得到。自相似指数大,则信号的自相似程度就高;

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论