小波变换发展和应用幻灯片.ppt

《小波变换发展和应用幻灯片.ppt》由会员分享，可在线阅读，更多相关《小波变换发展和应用幻灯片.ppt（39页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、小波变换发展和应用第1页，共39页，编辑于2022年，星期六 1988年 Mallat 提出的多分辨度分析理论，统一了几个不相关的领域：包括语音识别中的镜向滤波，图象处理中的金字塔方法，地震分析中短时波形处理等。当在某一个分辨度检测不到的现象，在另一个分辨度却很容易观察处理。例如：第2页，共39页，编辑于2022年，星期六 小波分析是纯数学、应用数学和工程技术的完美结合。从数学来说是大半个世纪“调和分析”的结晶（包括傅里叶分析、函数空间等）。小波变换是20世纪最辉煌科学成就之一。在计算机应用、信号处理、图象分析、非线性科学、地球科学和应用技术等已有重大突破，预示着小波分析进一步热潮的到来。第3

2、页，共39页，编辑于2022年，星期六 “小波分析”是分析原始信号各种变化的特性，进一步用于数据压缩、噪声去除、特征选择等。例如歌唱信号：是高音还是低音，发声时间长短、起伏、旋律等。从平稳的波形发现突变的尖峰。小波分析是利用多种“小波基函数”对“原始信号”进行分解。第4页，共39页，编辑于2022年，星期六小波的时间和频率特性小波的时间和频率特性小波的时间和频率特性小波的时间和频率特性运用小波基，可以提取信号中的运用小波基，可以提取信号中的“指定时间指定时间”和和“指定频率指定频率”的变化。的变化。时间：提取信号中时间：提取信号中“指定时间指定时间”（时间（时间A A或时间或时间B B）的变化

3、。顾名思义，小）的变化。顾名思义，小波在某时间发生的小的波动。波在某时间发生的小的波动。频率：提取信号中时间频率：提取信号中时间A A的比较慢速变化，称较低频率成分；而提取信号中的比较慢速变化，称较低频率成分；而提取信号中时间时间B B的比较快速变化，称较高频率成分。的比较快速变化，称较高频率成分。时间A时间B第5页，共39页，编辑于2022年，星期六 第6页，共39页，编辑于2022年，星期六 参考：M.Vetterli,”Wavelets and Subband Coding“,Prentice Hall PTR,1995 p.11第7页，共39页，编辑于2022年，星期六小波的小波的小波

4、的小波的3 3 个特点个特点个特点个特点 小波变换，既具有频率分析的性质，又能表示发生的时间。小波变换，既具有频率分析的性质，又能表示发生的时间。有利有利于分析确定时间发生的现象。（傅里叶于分析确定时间发生的现象。（傅里叶变换只变换只具有频率分析的性具有频率分析的性质）质）小波变换的多分辨度的变换，有利于各分辨度不同特征的提小波变换的多分辨度的变换，有利于各分辨度不同特征的提取（图象压缩，边缘抽取，噪声过滤等）取（图象压缩，边缘抽取，噪声过滤等）小波变换比快速小波变换比快速FourierFourier变换还要快一个数量级。变换还要快一个数量级。信号长度为信号长度为MM时，时，FourierFo

5、urier变换变换（左）（左）和和小波变换（右）计算复杂性分小波变换（右）计算复杂性分别如下公式：别如下公式：第8页，共39页，编辑于2022年，星期六小波基表示发生的时间和频率小波基表示发生的时间和频率小波基表示发生的时间和频率小波基表示发生的时间和频率“时频局域性时频局域性”图解：图解：FourierFourier变换的基（上）小波变换基（中）变换的基（上）小波变换基（中）和时间采样基（下）的比较和时间采样基（下）的比较 傅里叶变换(Fourier)基小波基时间采样基第9页，共39页，编辑于2022年，星期六信号的时频分析：信号时频分析的重要性：-时间和频率是描述信号的两个最重要的物理量。

6、时间和频率是描述信号的两个最重要的物理量。-信号的时域和频域之间具有紧密的联系。信号的时域和频域之间具有紧密的联系。信号时频分析的主要方法：第10页，共39页，编辑于2022年，星期六反映傅立叶变换缺点的一个例子：第11页，共39页，编辑于2022年，星期六傅立叶变换的缺点：用傅立叶变换提取信号的频谱需要利用信号的全部时域信息。傅立叶变换没有反映出随着时间的变化信号频率成分的变化情况。傅立叶变换的积分作用平滑了非平稳信号的突变成分。第12页，共39页，编辑于2022年，星期六解决傅立叶变换缺点的方法：第13页，共39页，编辑于2022年，星期六窗口傅立叶变换（Gabor变换）：窗口傅立叶变换的

7、定义：-假设假设 f f(t)(t)L L2 2(R)(R)，则以，则以g g(t)(t)作为窗函数的窗口傅立叶作为窗函数的窗口傅立叶变换定义为：变换定义为：窗口傅立叶变换的物理意义：窗口傅立叶变换的物理意义：-若若g g(t)(t)的有效窗口宽度为的有效窗口宽度为D Dt t，则，则WFWFg g(,b,b)给出的是给出的是f f(t)(t)在在局部时间范围局部时间范围b-Db-Dt t/2,b+D/2,b+Dt t/2/2内的频谱信息。内的频谱信息。-有效窗口宽度有效窗口宽度D Dt t越小，对信号的时间定位能力越强。越小，对信号的时间定位能力越强。第14页，共39页，编辑于2022年，星

8、期六窗口傅立叶变换的频域性质：问题的提出：-窗口傅立叶变换窗口傅立叶变换WFWFg g(,b,b)=)=f 给出的是信给出的是信号在时域上的处理信息，一个很自然的问题是窗口号在时域上的处理信息，一个很自然的问题是窗口傅立叶变换在频域上是怎样处理信号的？傅立叶变换在频域上是怎样处理信号的？假设f(t)的傅立叶变换为F()，g,b(t)的傅立叶变换为G,b()，则根据Parseval定理有：WFg(,b)=/(2)窗口傅立叶变换频域上的物理意义：-若若G G()的有效窗口宽度为的有效窗口宽度为D D，则，则WFWFg g(,b,b)给出的给出的是是F F()在局部频率范围在局部频率范围 -D-D

9、/2,/2,+D+D /2/2内的频内的频谱信息。谱信息。-有效窗口宽度有效窗口宽度D D 越小，对信号的频率定位能力越强。越小，对信号的频率定位能力越强。第15页，共39页，编辑于2022年，星期六窗口傅立叶变换的性能分析：问题的提出：-窗口傅立叶变换是否既具有强的时间定位能力，又窗口傅立叶变换是否既具有强的时间定位能力，又具有强的频率定位能力？具有强的频率定位能力？-选择什么样的窗函数才能使得窗口傅立叶变换具有选择什么样的窗函数才能使得窗口傅立叶变换具有好的性能？好的性能？解决问题的思想：-从物理意义上来看从物理意义上来看D Dt t和和D D 是矛盾的，因此先定义是矛盾的，因此先定义D

10、Dt t和和D D 后，再计算后，再计算D Dt t和和D D 的乘积用以作为判断窗口傅的乘积用以作为判断窗口傅立叶变换性能的依据。立叶变换性能的依据。第16页，共39页，编辑于2022年，星期六窗口傅立叶变换的性能分析：具体分析过程：-假设：假设：-定义：定义：第17页，共39页，编辑于2022年，星期六窗口傅立叶变换的性能分析：-计算计算D Dt t2 2DD 2 2：海森堡测不准原理第18页，共39页，编辑于2022年，星期六窗口傅立叶变换的性能分析：-等号成立条件：等号成立条件：结论：结论：-窗口傅立叶变换的时间分辨率和频率分辨率不可能同时提高，窗口傅立叶变换的时间分辨率和频率分辨率不

11、可能同时提高，只能以一种分辨率的降低来换取另一种分辨率的提高。只能以一种分辨率的降低来换取另一种分辨率的提高。-以高斯函数作为窗函数相对来说综合效果最好。以高斯函数作为窗函数相对来说综合效果最好。第19页，共39页，编辑于2022年，星期六解决窗口傅立叶变换缺点的方法：问题的提出：-窗口傅立叶变换窗口没有自适应性，只适合分析所窗口傅立叶变换窗口没有自适应性，只适合分析所有特征尺度大致相同的信号，不适于分析多尺度信有特征尺度大致相同的信号，不适于分析多尺度信号和突变过程。号和突变过程。解决方法：-引入窗口变化机制，同时求各种窗口大小下的变换，引入窗口变化机制，同时求各种窗口大小下的变换，这样变换

12、系数中就同时包含各种特征尺度下信号的这样变换系数中就同时包含各种特征尺度下信号的信息。信息。第20页，共39页，编辑于2022年，星期六小波变换的分类：连续小波变换-时间、控制窗口大小的参数和时移参数都连续的小时间、控制窗口大小的参数和时移参数都连续的小波变换。波变换。离散参数小波变换-时间连续，控制窗口大小的参数和时移参数离散的时间连续，控制窗口大小的参数和时移参数离散的小波变换。小波变换。离散小波变换-时间、控制窗口大小的参数和时移参数都离散的小时间、控制窗口大小的参数和时移参数都离散的小波变换。波变换。第21页，共39页，编辑于2022年，星期六连续小波变换：连续小波变换的定义：-假设信

13、号假设信号 f f(t)(t)L L2 2(R)(R)，则它的连续小波变换定义为：，则它的连续小波变换定义为：尺度伸缩参数时间平移参数归一化因子第22页，共39页，编辑于2022年，星期六连续小波变换：连续小波变换的物理意义：-时域上的意义：数学显微镜（一组有效宽度不同的时域上的意义：数学显微镜（一组有效宽度不同的窗口傅立叶变换的汇集）窗口傅立叶变换的汇集）第23页，共39页，编辑于2022年，星期六连续小波变换：-频域上的意义：频域上的意义：若若f f(t)(t)的傅立叶变换为的傅立叶变换为F F()，a,ba,b(t)(t)的傅立叶变换为的傅立叶变换为a,ba,b()，则根据则根据Pars

14、evalParseval定理，有：定理，有：第24页，共39页，编辑于2022年，星期六连续小波变换：-“恒恒Q Q性质性质”：假设假设(t)(t)的中心为的中心为t t0 0，有效宽度为，有效宽度为D Dt t；()的中心为的中心为 0 0，有，有效宽度为效宽度为D D；则；则 a,ba,b(t)(t)提取的是提取的是f f(t)(t)在窗口在窗口b+atb+at0 0-aD-aDt t/2,/2,b+atb+at0 0+aD+aDt t/2|/2|中的性质，相应地从频域上说中的性质，相应地从频域上说a,ba,b()提取地是提取地是F F()在窗口在窗口 0 0/a-D/a-D/(2a),/

15、(2a),0 0/a+D/a+D/(2a)/(2a)中的性质，因此对于中的性质，因此对于小波来说时域窗口宽度和频域窗口宽度的乘积始终为小波来说时域窗口宽度和频域窗口宽度的乘积始终为D Dt tD D。第25页，共39页，编辑于2022年，星期六母小波的例子：Harr小波：第26页，共39页，编辑于2022年，星期六母小波的例子：Mexico草帽小波：第27页，共39页，编辑于2022年，星期六母小波的例子：Morlet小波：第28页，共39页，编辑于2022年，星期六连续小波变换的逆变换：连续小波变换逆变换存在的可能性：-窗口宽度任意调节，在时域上或频域上能完全恢复窗口宽度任意调节，在时域上或

16、频域上能完全恢复出信号的信息。出信号的信息。-连续小波变换结果有很大的冗余度。连续小波变换结果有很大的冗余度。以a,b(t)为变换核的连续小波变换的逆变换：-假设假设 f f(t)(t)、(t)(t)L L2 2(R)(R)：第29页，共39页，编辑于2022年，星期六连续小波变换的逆变换：-证明思路：证明思路：第30页，共39页，编辑于2022年，星期六母小波的容许条件：从逆变换公式可以看出母小波的容许条件为：第31页，共39页，编辑于2022年，星期六连续小波变换的逆变换的其他形式：以不等于a,b(t)的小波函数a,b(t)为变换核的连续小波变换的逆变换：互为对偶关系第32页，共39页，编

17、辑于2022年，星期六小波变换的分类：连续小波变换-时间、控制窗口大小的参数和时移参数都连续的小时间、控制窗口大小的参数和时移参数都连续的小波变换。波变换。离散参数小波变换-时间连续，控制窗口大小的参数和时移参数离散的时间连续，控制窗口大小的参数和时移参数离散的小波变换。小波变换。离散小波变换-时间、控制窗口大小的参数和时移参数都离散的小时间、控制窗口大小的参数和时移参数都离散的小波变换。波变换。第33页，共39页，编辑于2022年，星期六尺度和时移参数的离散化：问题的提出：-连续小波变换中含有很多冗余信息，冗余信息不利连续小波变换中含有很多冗余信息，冗余信息不利于对信号的分析和处理。于对信号

18、的分析和处理。-连续小波变换的计算量也大。连续小波变换的计算量也大。-由于连续小波变换中有冗余信息，可能对尺度和时由于连续小波变换中有冗余信息，可能对尺度和时移参数进行离散化后仍可重构信号。移参数进行离散化后仍可重构信号。尺度和时移参数离散化要解决的问题：-尺度和时移参数要怎样离散化？尺度和时移参数要怎样离散化？-尺度和时移参数离散化后要想重构信号对小波函数尺度和时移参数离散化后要想重构信号对小波函数应有什么样的要求？应有什么样的要求？第34页，共39页，编辑于2022年，星期六尺度和时移参数的离散化：尺度和时移参数离散化的方法：-尺度参数的离散化：尺度参数的离散化：a=aa=a0 0j j,

19、j ,j Z Z （通常取（通常取a a0 0的值为的值为2 2，称为二进小波），称为二进小波）-时移参数的离散化：取决于尺度参数时移参数的离散化：取决于尺度参数 b=kab=ka0 0j j,j,k ,j,k Z Z第35页，共39页，编辑于2022年，星期六尺度和时移参数的离散化：离散化后的小波变换：怎样选择小波函数才能够重构信号：怎样选择小波函数才能够重构信号：-小波函数仍应满足连续小波变换中的容许条件。小波函数仍应满足连续小波变换中的容许条件。-小波函数的选择与离散化的程度有关系，离散化参数取小波函数的选择与离散化的程度有关系，离散化参数取样间隔很小时对小波函数的限制也小，而离散化参数

20、的样间隔很小时对小波函数的限制也小，而离散化参数的取样间隔很大是对小波函数的限制也会很大。取样间隔很大是对小波函数的限制也会很大。第36页，共39页，编辑于2022年，星期六尺度和时移参数的离散化：-重构信号小波函数应满足的条件（框架理论）：重构信号小波函数应满足的条件（框架理论）：对任意的对任意的 f f(t)(t)L L2 2(R)(R)，称，称 j,kj,k 为一个框架，如果存为一个框架，如果存在正参数在正参数A A和和B B（0 0 A A B B ），使得：），使得：分析小波合成小波第37页，共39页，编辑于2022年，星期六尺度和时移参数的离散化：框架的一种特殊情况紧框架：-定义：定义：A=BA=B的框架称为紧框架。的框架称为紧框架。-性质：性质：分析小波的对偶是它本身，类似于正交变换。第38页，共39页，编辑于2022年，星期六标准正交小波基：标准正交小波基的优点：-变换系数没有冗余，能够很好地反映信号的性质。变换系数没有冗余，能够很好地反映信号的性质。-标准正交小波基与它的对偶相同。标准正交小波基与它的对偶相同。-计算简单：计算简单：假设：假设：-以下使用的都是二进小波，即小波函数的形式为：以下使用的都是二进小波，即小波函数的形式为：第39页，共39页，编辑于2022年，星期六