小波分析ppt课件

1、1.1.小波和小波变换小波和小波变换(Wavelet and Wavelet (Wavelet and Wavelet Transform)Transform) x x x 几点约定几点约定:我们的讨论范围只是函数空间我们的讨论范围只是函数空间 L2R）；）；小写小写 是时间信号，大写是时间信号，大写 是其是其Fourier变换；变换；尺度函数总是写成尺度函数总是写成 （时间域和（时间域和 （频率（频率域）；域）；小波函数总是写成小波函数总是写成 （时间域和（时间域和 （频（频率域）。率域）。1.1 小波小波(Wavelet) 小波就是空间小波就是空间L2(R)L2(R)中满足下述条中满足下述

2、条件的函数或者信号件的函数或者信号 : : x Rdxx2 *2RdC 这时，这时， 也称为小波母函数，也称为小波母函数，(2) 称为容称为容许性条件。许性条件。 x （1）（2）连续小波连续小波 abxaxba 1,函数：函数：为由小波母函数为由小波母函数 生成的依赖生成的依赖于参数于参数a,b的连续小波，简的连续小波，简称为小波。称为小波。 x （3）注释注释注释：如果小波母函数注释：如果小波母函数 的的Fourier 变换变换 在原点在原点 是连续是连续 的，那么公式的，那么公式(2)阐明阐明 ， x 0 00 于是于是 x dxR 0这说明函数这说明函数 有波动的特点，公式有波动的特点

3、，公式(1)又说明函数又说明函数 有衰减的特点，因此，有衰减的特点，因此，称函数称函数 为为“小波小波”。 x x x 1.2 小波变换小波变换(Wavelet Transform)对于任意的函数或者信号对于任意的函数或者信号 ，其，其小波变换为小波变换为 RLxf2 RRbafdxabxxfadxxxfbaW 1,（4）性质性质这样定义的小波变换具有下列性质：这样定义的小波变换具有下列性质：Plancherel恒等式：恒等式： 22,RgfRadadbbaWbaWdxxgxfC 小波变换的逆变换公式：小波变换的逆变换公式： 22,1RbafadadbxbaWCxf （5）（6）性质性质吸收公

4、式：当吸收条件吸收公式：当吸收条件 0202 dd成立时，有吸收的成立时，有吸收的Plancherel恒等式恒等式 1220Cf x g x dxW abW abdbdaafg,（7）（8）性质性质吸收的逆变换公式吸收的逆变换公式 02,2adadbxbaWCxfbaf （9）1.3.二进小波和二进小波变换二进小波和二进小波变换(Dyadic Wavelet Transform) 如果小波函数如果小波函数 满足稳定性条件满足稳定性条件 x BAj 2 (10)则称则称 为二进小波，对于任意的整数为二进小波，对于任意的整数k,k,记记 x kkxxk2212 (11)逆变换逆变换对于任意的对于任

5、意的 ，其二进小波变换为：，其二进小波变换为： RLxf2 bfdxxxfbWkRkkkf2221 这时，逆变换公式是这时，逆变换公式是 kRkkkfdbbxbWxf22 （12）（13）重构小波重构小波其中其中 的的Fourier变换满足变换满足 x 122 kkk 称为二进小波称为二进小波 的重构小波，比如可取：的重构小波，比如可取： x kk22 （14）（15）设小波为设小波为 ，对于任意的整数，对于任意的整数k 和和j，记，记1.4. 正交小波和小波级数正交小波和小波级数(Orthonormal Wavelet) 构成空间构成空间 的标准正交基，则称的标准正交基，则称 是正是正交小波

6、。交小波。 x jxxkkjk 222, ZZjkjxxkkjk , ; 222, RL2 x 如果函数族如果函数族(16)(17)小波级数小波级数这时，逆变换公式就是小波级数这时，逆变换公式就是小波级数 kjjkjkxxf, (18) 其中小波系数其中小波系数 的算法是的算法是 jk, Rjkjkjkdxxxff, , (19)连续和离散统一连续和离散统一上的取值，因此，小波系数上的取值，因此，小波系数 实际上是实际上是信号信号f(x)的离散小波变换。其实，这也是小的离散小波变换。其实，这也是小波变换迷人的风采之一：波变换迷人的风采之一： baWf, jkk 2 , 2jk, 小波系数是信号

7、小波系数是信号f(x)的小波变换的小波变换 在二在二进离散点进离散点(20)连续变换和离散变换形式统一；连续变换和离散变换形式统一；连续变换和离散变换都适合全体信号；连续变换和离散变换都适合全体信号； 2. 小波分析和时小波分析和时-频分析频分析(Time-Frequency Analysis )2.1 窗口窗口Fourier变换和变换和Gabor变换变换(Windowed Fourier Transform and Gabor Transform) D.Gabor在在1946年开创时年开创时-频分析的先河提出频分析的先河提出Gabor Transform一般的时一般的时-频分析是频分析是Wi

8、ndowed Fourier TransformShort-Time Fourier TransformWindowed Fourier Transform称为信号称为信号 的窗口的窗口Fourier变换，其变换，其中的函数中的函数 称为窗口函数，一般要求称为窗口函数，一般要求是：是： RfdxxixxgxfxS0000exp, RLxf2 RLxg2 1 Rdxxg具体地具体地(21)Gabor Transform 4exp212xxg RfdxxSF0000, D.Gabor取取(22)是是Gaussian函数，对应的变换称为函数，对应的变换称为Gabor变变换换(1946)。对于。对于G

9、abor变换，存在如下的频变换，存在如下的频率再分割公式：率再分割公式：(23)物理解释物理解释Gabor变换变换 是是信号信号 在在x=x0点点“附近的频率为附近的频率为 的的频率成分；频率成分；只要把信号只要把信号 在各个时间点在各个时间点“附近的频率为附近的频率为 的频率成分全部累加起来，理所当的频率成分全部累加起来，理所当 然就然就应 该 是 这 个 信 号 的 频 率 为应 该 是 这 个 信 号 的 频 率 为 的频率成的频率成 分；分；Gabor变换变换 可以认为是信号可以认为是信号f(x)的另一种等价描述的另一种等价描述(因为因为Fourier变换是信号的等价描述变换是信号的等

10、价描述) RLxf2 0 00, xSf RLxf2 0 0 00, xSf局限局限 Gabor变换没有变换没有“好的即可好的即可以构成标架或者正交基离散以构成标架或者正交基离散形式；形式； Gabor变换没有快速算法：比变换没有快速算法：比如没有类似于离散如没有类似于离散Fourier变变换之换之FFT的快速数值算法；的快速数值算法； 遗憾的是，遗憾的是，Gabor变换存在如下局限：变换存在如下局限： Appendix A Fig.1. Gabor变换的固定时变换的固定时-频窗口频窗口t00t1t12.2. 时时-频分析频分析(Time-Frequency Analysis）时时- -频分析

11、本质上是信号描述、分析和处频分析本质上是信号描述、分析和处理的一种方法，它给信号的理的一种方法，它给信号的“最优描述问题最优描述问题提供一种解决方案。提供一种解决方案。R.Balian(1981)R.Balian(1981)早在早在八十年代就清清楚楚地描述了这个问题：八十年代就清清楚楚地描述了这个问题：在通讯理论中，人们对于在给定的时在通讯理论中，人们对于在给定的时间内，把一个信号表示成间内，把一个信号表示成“每一个都每一个都同时具有足够确定的位置及频率的谐同时具有足够确定的位置及频率的谐波的叠加这种信号的描述方法极感波的叠加这种信号的描述方法极感兴趣兴趣 最优描述问题最优描述问题有用的信息总

12、是同时被所发射信号的频率有用的信息总是同时被所发射信号的频率特性与信号的时间结构所传递，最好的例特性与信号的时间结构所传递，最好的例子是演奏音乐；子是演奏音乐；把信号表成时间的函数其频率特征无法突把信号表成时间的函数其频率特征无法突出，而出，而FourierFourier分析又无法标定各个分量分析又无法标定各个分量发射的瞬时位置和持续时间；发射的瞬时位置和持续时间；“最优描述应该综合这两种描述的优点，最优描述应该综合这两种描述的优点，并用一个离散的刻画来表示，以适应信息并用一个离散的刻画来表示，以适应信息理论和计算机处理的需要。理论和计算机处理的需要。 Wigner分布函数分布函数Wigner

13、分布函数是信号时分布函数是信号时-频频分析的另一种具体的解决途径。信分析的另一种具体的解决途径。信号号f(x)的的Wigner分布函数是著名分布函数是著名理论物理学家理论物理学家E.P.Wigner在在1932年提出来的，定义是：年提出来的，定义是： RfdixfxfxW exp22,(24) 显然，这是一个实的二元函数显然，这是一个实的二元函数 。性质性质 22, RRRgfdxxgxfddxxWxW 22,xfdxWRf 2, FdxxWRf Wigner分布函数有如下性质：分布函数有如下性质：(25)(26)(27)Wigner分布函数的物理意义分布函数的物理意义Wigner分布函数的P

14、lancherel恒等式成立；Wigner分布函数 标明信号的瞬时频率的位置；Wigner分布函数 标明信号的瞬时位置的频率。 , xWf在能量的意义下，在能量的意义下，Wigner分布函数的分布函数的物理意义是：物理意义是： , xWf x x x Wigner分布函数理论的局限分布函数理论的局限Wigner分布函数的三个局限：分布函数的三个局限： Wigner分布函数分布函数 只记忆信号的部分信息；只记忆信号的部分信息；Wigner分布函数分布函数 没有有效的重建算法；没有有效的重建算法；Wigner分布函数分布函数 的的“瞬时是渐近意义的。瞬时是渐近意义的。 , xWf , xWf ,

15、xWf2.3. 小波的时小波的时-频分析频分析(Wavelets Time- Frequency Analysis) 小波变换是一种时小波变换是一种时-频描述，它的信息记频描述，它的信息记忆是完全的，是一种等价的变换描述，具忆是完全的，是一种等价的变换描述，具有独特的时有独特的时频分析性质。引入记号：频分析性质。引入记号： RgdxgxgxE22(28)中心中心半径半径 21222 RggdxgxgEx(29)对于对于 ，如果满足条件，如果满足条件:窗口函数及说明窗口函数及说明则称之为窗口函数，则称之为窗口函数， 和和 分别称为它的分别称为它的时间中心和时间半径，而时间中心和时间半径，而 和和

16、 分别称为分别称为它的谱中心和谱半径。它的谱中心和谱半径。 RLxg2 Rdxxxg2gEg GEG 0221 gxgxg阐明：中心和半径是下述分布的期望和均方差小波的时小波的时-频中心与半径频中心与半径 aEEaEbEbaba, aababa, 2.3.2.小波的时-频半径 2.3.1.小波的时小波的时-频中心频中心（29）（30）2.3.3. 小波的时小波的时-频窗频窗 aaEaaEaaEbaaEbEEEEbabababababababa , + , , , , (32)Appendix B Fig.2.小波在时小波在时-频相平面上的窗频相平面上的窗t00t12t12.3.4. 小波的时小

17、波的时-频特性频特性 小波时小波时-频窗的面积恒等于频窗的面积恒等于 ； 小波的时小波的时-频窗是时频窗是时-频相平面中的可频相平面中的可变的矩形；变的矩形； 小波时小波时-频窗的变化规律：频窗的变化规律： 4（1尺度参数尺度参数a增大时，小波的时窗变宽，增大时，小波的时窗变宽，同时，它的主频变低，频窗变窄；同时，它的主频变低，频窗变窄；（2尺度参数尺度参数a减小时，小波的时窗变窄，减小时，小波的时窗变窄，同时，它的主频变高，频窗变宽；同时，它的主频变高，频窗变宽； 小波的频率分辨率小波的频率分辨率小波分析具有固定的相对频率分辨率小波分析具有固定的相对频率分辨率 22EaaE(33) 主频变低

18、时，频窗变窄，频率分辨率提高；主频变低时，频窗变窄，频率分辨率提高；主频变高时，频窗变宽，频率分辨率降低；主频变高时，频窗变宽，频率分辨率降低；高频时出现较低的频率分辨率难题！）。高频时出现较低的频率分辨率难题！）。 小波的频带特性小波的频带特性 (1)小波变换处理频域的方式完全不同于经小波变换处理频域的方式完全不同于经典的典的Fourier变换，任何小波本质上都变换，任何小波本质上都是以频带的形式出现在频域中，这样避是以频带的形式出现在频域中，这样避免了许多理论和计算上的麻烦；免了许多理论和计算上的麻烦；(2)二进小波频域划分的特色：将参数二进小波频域划分的特色：将参数a按按二进方式离散化为

19、二进方式离散化为kka 2 x E 3选择二进小波选择二进小波 满足满足二进小波二进小波 的主频是的主频是二进小波的分频特性二进小波的分频特性 xbk,2 kk23 212 , 2kk kkk212 , 2 , 0(34)所在的频带是所在的频带是当当k取遍全体整数时，这些频带正好分离覆取遍全体整数时，这些频带正好分离覆盖正频轴，即盖正频轴，即这就是著名的二进小波频带划分技术。这就是著名的二进小波频带划分技术。2.4. 正交小波的时正交小波的时-频分析频分析Orthonormal Wavelets Time-Frequency Analysis对于正交小波对于正交小波 ， x kjjkjkxxf

20、, (35) 2, , ; Zjkxjk 其中系数是其中系数是是一个标准正交基，所以，对于任何信号是一个标准正交基，所以，对于任何信号f(X)f(X)，可以展开成小波级数：，可以展开成小波级数： Rjkjkdxxxf, (36)正交小波的吸收谱正交小波的吸收谱jk, baWf, jkk 2 , 2由小波变换的定义可知，正交小波由小波变换的定义可知，正交小波级数的系数级数的系数 正好是信号正好是信号f(x)的小波的小波变换变换 在二进离散点：在二进离散点：(37)上的取值。这说明：对于正交小波来说，任上的取值。这说明：对于正交小波来说，任何信号在二进离散点上的小波变换包含了它何信号在二进离散点上

21、的小波变换包含了它的小波变换的全部信息，所以的小波变换的全部信息，所以正交小波具有优美的谱吸收特点。正交小波具有优美的谱吸收特点。小波变换与小波变换与FourierFourier变换变换 Fourier变换：变换：对于任何信号对于任何信号f(x)，只有当它是时，只有当它是时间有限时，它的谱间有限时，它的谱F()(Fourier变换变换)才是频率吸收的；才是频率吸收的；反过来，只有当它是频域有限时，反过来，只有当它是频域有限时，f(x)才是时间吸收的才是时间吸收的;小波变换小波变换: :对于正交小波分析来说，任何信号的对于正交小波分析来说，任何信号的正交小波谱都是谱吸收的，即二维小正交小波谱都是

22、谱吸收的，即二维小波谱所包含的信息完全被二进离散点波谱所包含的信息完全被二进离散点上的谱吸收。上的谱吸收。一点评论一点评论正交小波变换谱的完全吸收性为小波变正交小波变换谱的完全吸收性为小波变换的理论分析、数值计算和各种应用提换的理论分析、数值计算和各种应用提供了极大的方便。同时，这些离散的小供了极大的方便。同时，这些离散的小波谱点，本质上意味着时波谱点，本质上意味着时-频分析中频谱频分析中频谱分析的频带统计意义下的区间），因分析的频带统计意义下的区间），因此，小波分析成功地实现了人们梦寐以此，小波分析成功地实现了人们梦寐以求的求的“频带信息的点处理方式频带信息的点处理方式”；在在(a,b)-W

23、(a,b)给出的二维小波谱空间，给出的二维小波谱空间，二进离散小波谱点的分布规律可以用二进离散小波谱点的分布规律可以用Appendix C Fig.3. 加以说明。加以说明。 Appendix C Fig.3.正交小波的点谱吸收特性正交小波的点谱吸收特性01234567891011121314150123456701230103. 正交小波和多分辨分析正交小波和多分辨分析(Orthonormal Wavelet and Multiresolution Analysis)多分辨分析：多分辨分析： 上的一列闭的上的一列闭的线性子空间线性子空间 和一个函数和一个函数 共同称为一个多分辨分析，如果它共

24、同称为一个多分辨分析，如果它们满足如下的五个要求：们满足如下的五个要求： RL2 ZkVk ; x 3.1. 多分辨分析多分辨分析(Multiresolution Analysis)多分辨分析多分辨分析 0 ZkkV RLVZkk2 12 kkVxfVxf2.唯一性公理：唯一性公理：3.稠密性公理：稠密性公理：4.伸缩性公理：伸缩性公理：(39)(40)(41)5.构造性公理：构造性公理： Zkkx ; (42)生成生成V0的标准正交基。其中的函数的标准正交基。其中的函数 称为称为尺度函数尺度函数(Scale Function)。 x 1.单调性公理：单调性公理：1 kkVV(38)图像的多分

25、辨分析图像的多分辨分析多分辨分析多分辨分析(Multiresolution Analysis)方法，在计算机科学和信号处方法，在计算机科学和信号处理中，特别是在图像分析中，通常称为理中，特别是在图像分析中，通常称为多尺度分析方法多尺度分析方法(Multiscale Analysis) ，在小波分析建立之前就已经得到了一些在小波分析建立之前就已经得到了一些理论研究和应用，这推动了小波变换理理论研究和应用，这推动了小波变换理论的产生和完善。实际上，信号论的产生和完善。实际上，信号f(x)在在子空间子空间Vk上的正交投影上的正交投影fk(x)是是 ZkjkjkkxCxf, 图像的多分辨分析续）图像的

26、多分辨分析续）正交投影正交投影fk(x)正好是原象正好是原象f(x)在一定的分辨率在一定的分辨率之下的模糊象，公式之下的模糊象，公式(40)阐明，当分辨率足够阐明，当分辨率足够高时，模糊象和原象重合，即高时，模糊象和原象重合，即 xfxfkk 因此，对因此，对fk(x)的分析实际是对原象的多种分的分析实际是对原象的多种分辨率的分析。多分辨分析的困难在于如何从辨率的分析。多分辨分析的困难在于如何从低分辨率的模糊象有效地添加恰当的细节，低分辨率的模糊象有效地添加恰当的细节，得到正确的高分辨率下的模糊象。这些问题得到正确的高分辨率下的模糊象。这些问题的研究都属于多分辨分析的范围。的研究都属于多分辨分

27、析的范围。 3.2. 小波构造小波构造(Y.Meyer and S.Mallat, 1988)称之为尺度方程。系数列称之为尺度方程。系数列 叫低通滤波系数。叫低通滤波系数。 ZkVk ; x ZlZkhk2 ; Zkkkxhx22 ZlZkhk2 ; 假设假设 和函数和函数 是一个多分辨分是一个多分辨分析，那么，必然存在一列析，那么，必然存在一列 系数，使得系数，使得(43)构造定理构造定理(Y.Meyer and S.Mallat, 1988)令令 ， 并构造并构造 kkkhg 11Zk Zkkkxgx22 (44) 是L2(R)的标准正交基 2, ; Zjkxjk 则有如下结论则有如下结论

28、:(45) 是Vk在Vk+1中的正交补kkkWVV 1 ZjxSpanWjkk ; , 构造定理的延伸结果构造定理的延伸结果0 11 pVWWWVpkpkkkk， 01ppkkWV 02ppkkWVRL kkWRL2(46)(47)(49) (48)4. 多分辨分析和金字塔算法多分辨分析和金字塔算法(Multiresolution Analysis and Pyramid Algorithms)4.0. 记号记号 （Notation）：）： RmkmkRmkmkdxxxfDdxxxfC, 分别表示信号的趋势和波动或者模糊象和细节分别表示信号的趋势和波动或者模糊象和细节(50)4.1. 小波分解

29、算法小波分解算法(Decomposition Algorithms of Wavelet) ZmmkjmjkZmmkjmjkCgDChC,2,1,2,1(51)4.2. 小波重建算法小波重建算法(Reconstruction Algorithms of Wavelet) ZmmkmjmkmjjkDgChC, 12, 12,(52) 4.3. 金字塔算法金字塔算法(Pyramid Algorithms) ZmmkmkkZmmkmkkxDxgxCxf, (53) 引入记号引入记号:它们的几何意义分别是原信号它们的几何意义分别是原信号 在在子空间子空间Vk和和WK上的正交投影，且它们是相上的正交投影

30、，且它们是相互正交的。由多分辨分析的意义可得互正交的。由多分辨分析的意义可得 RLxf2 xgxfxfkkk 1(54)4.3.1. 分解金字塔算法分解金字塔算法(Decomposition Pyramid Algorithms)信号的分解信号的分解(Decomposition of Signal) 121 + + + + kkkk mk mf xfxfxfxfx 121 kkk mk mgxgxgxgx 空间的分解空间的分解空间的分解空间的分解(Decomposition of The Subspace) 1+21121mkmkkkmkmkkkkWWWWVVVVV 系数的分解系数的分解系数的

31、分解系数的分解(Decomposition of The Coefficients) mkkkkmkHHkHkHkHkDDDDGGGGGCCCCC 3213214.3.2. 重建金字塔算法重建金字塔算法(Reconstruction Pyramid Algorithms)信号的重建信号的重建(Reconstruction of Signal) xxxgxxfxfxfxfxfmkmkkkmkmkkkk g g g + + + + 1+21121空间的重建空间的重建空间的重建空间的重建(Reconstruction of Subspace) m121121 kmkkkmkmkkkkWWWWVVVV

32、V系数的重建系数的重建mkkkmkHkHkHkDDDGGGCCCC 2121 系数的重建(Reconstruction of The Coeffients)信号的小波分解和合成算法信号的小波分解和合成算法有限数字信号的高低通滤波器有限数字信号的高低通滤波器矩阵分解算法矩阵分解算法矩阵合成算法矩阵合成算法有限数字信号的小波变换编码有限数字信号的小波变换编码数字信号小波编码数据量关系数字信号小波编码数据量关系小波应用基本模式小波应用基本模式数字图像二维小波编码数字图像二维小波编码数字图像二维小波重建数字图像二维小波重建数字图像的矩阵小波变换数字图像的矩阵小波变换5.Malvar小波小波(H.S.M

33、alvar 1987）（R.Coifman and Y.Meyer 1991)5.1 Malvar小波小波(H.S.Malvar 1987) 选择窗口函数选择窗口函数 满足如下要求：满足如下要求： x 3or xx 0 x xxx 2 , 10 x 122 xx 时时时时Malvar小波基构造小波基构造Malvar小波基是函数族小波基是函数族 ZlkorZlkxulk2 ,0 ,.,3 ,2 , 1 ; , 1,2,= , 12 2sin22 0=, 2 2 , 2 , 1 , 2 2cos22,kZlxklxkZllxkZlxklxxulk (55)阐明阐明容易验证，上述函数族构成容易验证，

34、上述函数族构成L2(R)的的标准正交基。一般称这个函数族的小波为标准正交基。一般称这个函数族的小波为Malvar小波。小波。Malvar小波和离散余弦变换小波和离散余弦变换(DCT)、离散正弦变换、离散正弦变换(DST)有许多相似之有许多相似之处，根本的差别在于，处，根本的差别在于，Malvar小波是真正小波是真正局部化了的离散余弦变换和离散正弦变换局部化了的离散余弦变换和离散正弦变换分析，同时，它还具有变换结果的递推数分析，同时，它还具有变换结果的递推数值算法。值算法。 让人们惊奇的是，物理学家让人们惊奇的是，物理学家K.Wilson和和数学家数学家I.Daubechies也得到了极其相似的

35、结果。也得到了极其相似的结果。但是，他们两人和但是，他们两人和Malvar的工作之间并没有的工作之间并没有必然的逻辑的关系。必然的逻辑的关系。K.Wilson的想法是，对的想法是，对于实数轴的长度是于实数轴的长度是2的等长划分，按照各个的等长划分，按照各个区间的奇偶变化，分别轮番使用离散余弦变区间的奇偶变化，分别轮番使用离散余弦变换 和 离 散 正 弦 变 换 进 行 信 号 分 析 ；换 和 离 散 正 弦 变 换 进 行 信 号 分 析 ；I.Daubechies的想法是，不仅如此，而且必须的想法是，不仅如此，而且必须加以局部化，局部化因子是同一个函数加以局部化，局部化因子是同一个函数 的

36、的2倍整数平移，只不过要求函数和它的倍整数平移，只不过要求函数和它的Fourier变换都是指数衰减的并使得前述函数变换都是指数衰减的并使得前述函数族构成族构成 的标准正交基的标准正交基 。 RL25.2 Malvar小波小波(R.Coifman and Y.Meyer 1991) 选择选择 和和 并构造窗口函数列并构造窗口函数列 满足满足: ,.3, 2, 1, 0 ; kAk ,.3, 2, 1, 0 ; 0 kk ,.3, 2, 1, 0 ; kxk 11101kkkkkkkkkAAlAAAAA 窗函数的构造窗函数的构造 kkkkkkkkkkkkkkkkkkkkkAAxAAxxAxorAx

37、xAxAx122111110 110实际上，函数实际上，函数 本质上是区间本质上是区间 的的特征函数的光滑化特征函数的光滑化 xk 1 , kkAAAppendix D Fig.4.窗函数的形状示意图窗函数的形状示意图Ak-1AkAk+1Ak+kAk-kAk+1-k+1k(t)k-1(t)第一类第一类MalvarMalvar小波基小波基第一类第一类Malvar小波为：小波为： 0,1,2,= , ,21cos2,lZkAxllxlxukkkklk (56) 第二类第二类MalvarMalvar小波基小波基第二类第二类Malvar小波基为小波基为 1,2,3,= , 12 sin2 0= , 2

38、 11,2,3,= , 2 cos2,lZkAxllxlxulZkxlxulZkAxllxlxukkkklkkklkkkkklk (57) 6. 小波包小波包(Wavelet Packets)(R.Coifman and Y.Meyer and M.V.Wickerhauser 1992) 设设 和和 是一个多分辨分析且是一个多分辨分析且(43)和和(44)成立。记成立。记 ZkVk ; x xxxx 106.1 正交小波包正交小波包(Orthonormal Wavelet Packets)正交小波包的定义正交小波包的定义递推定义的函数族递推定义的函数族 ,.3 , 2 , 1 , 0 ; m

39、xm ZnmnmZnmnmnxgxnxhx2222122 ZnnxxSpankmknkmmk ; 222,; (58)(59)k是整数，是整数，m是自然数。是自然数。称之为小波包。引入记号称之为小波包。引入记号正交小波包定理正交小波包定理正交小波包定理正交小波包定理(Coifman and Meyer and (Coifman and Meyer and Wickerhauser 92)Wickerhauser 92)空间构造空间构造 是是 的标的标准正交基准正交基空间关系空间关系 (60) (60)特殊空间关系特殊空间关系 Znxnkm ; ,; mk 1221 mkmkmk122 mkmk

40、kkV 0kkW 1正交小波包的空间分割正交小波包的空间分割 小波包实现小波空间的再分割小波包实现小波空间的再分割120220120201222122726252423121111. kkkklllllklklklkkkkkkkkkW6.2小波包和时小波包和时-频分析频分析(Wavelet Packets andits Time-Frequency Analysis)利用正交小波的构造定理可知，子空间利用正交小波的构造定理可知，子空间Wk是是Vk在在Vk+1中的正交补：中的正交补： kkkWVV 1同时，根据小波的时同时，根据小波的时-频分析特性，可得下列频分析特性，可得下列关系关系: WGkkkk2212 ,VHkkk0 , 22正交小波实现有限频带的二进分割正交小波实现有限频带的二进分割 1111212 ,22 , 02 , 0 kpklllpkkpkkpkkWWVV 11212 ,22 , 0 klllklklkWV正交小波实现全频域的二进分割正交小波实现全频域的二进分割 111022,22,0,0 klllkppkkWVRL kkkkkWRL122 ,2 , 0正交小波包对二进频带的等分割正交小波包对二进频带的等分割 120122212211112,122 ,22,1llllllnkklklklklkkl