2022年2022年连续小波与二进小波变换 .pdf

《2022年2022年连续小波与二进小波变换 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年连续小波与二进小波变换 .pdf（35页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第 2 章连续小波与二进小波变换信号处理的应用随处可见，当你用数码相机拍摄照片，当你听着MP3 音乐，你有没有想过，正是信号处理技术使你轻松的获得娱乐。信号处理的主要任务是将现有的信号处理技术进行总结和抽象，信号处理的任务是认识客观世界中存在的信号的本质特征，并找出规律。从不同的角度去认识、分析信号有助于了解信号的本质特征。信号的表示方式很多，时间形式和频率形式是最重要的两种形式。时间形式是基于传感器采样得到的信号强度数据。这种数据很直观。除了时间以外，频率是一种表示信号特征最重要的方式。频率的表示方法是建立在傅里叶分析（Fourier Analysis）基础之上的，由于傅里叶分析是一种全局的

2、变换，要么完全在时间域，要么完全在频率域，因此无法表述信号的时频局部性质，而时频局部性质恰好是非平稳信号最基本和最关键的性质。为了分析和处理非平稳信号，在傅里叶分析理论基础上，提出并发展了一系列新的信号分析理论：短时傅里叶变换（Short Time Fourier Transform ） 、小波变换等。短时傅里叶变换是一种单一分辨率的信号分析方法，它选择一个时频局部化的窗函数。短时傅里叶变换窗函数受到W.Heisenberg 不确定准则的限制，时频窗的面积有下界。这也就从另一个侧面说明了短时傅里叶变换窗函数的时间与频率分辨率不能同时达到最优。Gabor 变换是海森伯不确定准则下的最优的短时傅里

3、叶变换。高斯窗函数是短时傅里叶变换同时追求时间分辨率与频率分辨率时的最优窗函数。具有高斯窗函数的短时傅里叶变换就是 Gabor 变换。与短时傅里叶变换一样，Gabor 变换也是单一分辨率的。小波变换使用小波基函数，时频窗面积不变，但形状可改变。小波函数根据需要调整时间与频率分辨率，具有多分辨分析（Multiresolution Analysis）的特点，克服了短时傅里叶变换分析非平稳信号单一分辨率的困难。小波变换是一种时间-尺度分析方法，而且在时间、尺度（频率）两域都具有表征信号局部特征的能力，在低频部分具有较高的频率分辨率和较低的时间分辨率，在高频部分具有较高的时间分辨率和较低的频率分辨率，

4、很适合于探测正常信号中夹带的瞬间反常现象并展示其成分。所以，小波变换被称为分析信号的显微镜。小波变换不会“一叶障目，不见泰山”，又可以做到“管中窥豹，略见一斑”。2.1 连续小波变换的定义定义 1 设)()(2RLxf，)(t是基本小波或母小波函数，则dxaxxfaaWRf)(1),(2.1) 称为)(xf的小波变换。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换2式中0a，称为尺度因子，反映位移，其值

5、可正可负。上式还可以用内积的形式表示，)(),(xgxf的内积表示为Rdxxgxfxgxf)()()(),(则连续小波变换用内积表示为)(),(),(xxfaWaf其中axaxa1)(是基本小波的位移和尺度伸缩。(2.1)式中，x，a和均为连续变量，因此称为连续小波变换（continuous wavelets transform ，CWT） 。关于式 (2.1)有以下几点补充说明（1）核函数连续傅立叶变换和小波变换在变换形式上，都可以视为信号和和核函数的内积。连续傅立叶变换的核函数为xie，基本小波( )x为小波变换的核函数。小波变换的核函数的形式没有唯一确定， 因此可以根据的应用需求选取不同

6、的小波基，从而达到最佳的处理效果。另外基本小波可以是实函数或者复函数，当基本小波是实函数时为实小波变换，当基本小波是复函数时就是复小波变换。例如21/42/22( )(1)e3xxx（Mexhat 小波）xiTxeex02)(=xiexeTxTx00sincos22（Morlet 小波）（2）频率与尺度因子小波作为时频分析工具，其将信号从时域变换到频域。那么小波变换的频率是怎样表现的呢？在连续小波变换中，尺度因子a的作用是将基本小波)(x作伸缩， 可以直观的得到其与频率是相关的。a愈大ax愈宽，相对于( ) x，其信号频率降低了。对于一个连续时间有限的小波)(x与)(xa关系如图名师资料总结

7、- - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换3)( x)(x2,aax2,aax小波函数的位移与伸缩在不同尺度下，小波的持续时间随尺度因子a的增大而增宽，幅度则与a成反比减小，但波的波形形状保持不变。（3）内积与卷积小波变换可以视为信号和小波核函数的内积，式（2.1）并不是卷积，由于内积dxxxfxxfR)()()(),(卷积dxxxfdxfdxfxxfRRR)()()()()()()()(比较两式， 如果()x满

8、足)()(xx，即如果)(x是偶函数，则上面的积分计算结果没有差别。从这一点出发，有些学者（ Mallat ）直接按卷积来定义小波变换，他们所采用的定义是dxaxxfaaWRf)(1),(*如果)()(xx，则有名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换4),(1),(*aWaaWff（4） 尺度因子)(xa前加因子a1的目的是使不同a值下)(xa的能量保持相等， 即基本小波的能量为RdxxE2)(

9、而)(xa得能量为EdxaxadxaxaERRa2211应该指出，目前小波函数的定义还不是唯一的，也就是说axaxa1)(不是小波函数族的唯一定义，有些学者主张对小波函数采用如下定义axaxa1)(其优点是在不同尺度下，可以保持各)(xa的频谱中幅频特征大小一致。事实上设)(x的傅立叶变换是)(?，则)(1xa的傅立叶变换是)(?)(?aaaa可见与)(?相比，只有频轴比例变化，没有幅度变化。2.2 连续小波变换的性质由于小波变换对)(xf而言是以)(x为核函数的线性变换，因此具有以下特性性质 1. 线性性如 果)(xf的 连 续 小 波 变 换 是),(aWf，)(xg的 连 续 小 波 变

10、 换 是),(aWg则名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换5)()()(21xgkxfkxh的连续小波变换是),(),(),(21aWkaWkaWgfh这是线性变换的基本特性，根据式（2.1）小波变换的定义很容易证明。性质 2. 时移性质如 果)(xf的 连 续 小 波 变 换 是),(aWf， 则)(0 xxf的 连 续 小 波 变 换 是),(0 xaWf。也就是)(xf的时移对应于小波变

11、换的移。证明令)()(0 xxfxg则),(aWg=dxaxxxfaR)(10令0 xxy，即0 xyx，dydx，则上式代入为),(aWg=dyaxyyfaR0)(1=),(0 xaWf性质 3. 尺度转换如果)(xf的连续小波变换是),(aWf，则xf的连续小波变换是,aWf，其中0。证明令xfxg)(，则),(aWg=dxaxxfaR1=dyayyfaR)(xy=dyayyfaR)(=dyayyfaR)(1=,aWf此性质表明，当信号)(xf做某一倍数伸缩时，其小波变换将在a，两轴上作同一比名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - -

12、- - - - 名师精心整理 - - - - - - - 第 5 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换6例的伸缩，但是不发生失真变形，这是使小波成为“数学显微镜”的重要理论依据。性质 4. 交叉项性质由于连续小波变换是线性变换，满足叠加性，因此不存在交叉项，但是由它引申出的能量分布函数2),(aWf却仍然有交叉项。设)()()(21xfxfxf，则)cos(),(),(2),(),(),(211212222fffffffaWaWaWaWaW式中1f，2f分别是),(1aWf和),(2aWf的幅角。证明),(aWf简记为fW，则2fW=)(2121212f

13、fffffWWWWWW=12212122ffffffWWWWWW21ffWW和12ffWW共轭，因此上式可以简化为)cos(2211212222fffffffWWWWW由上式可见， 小波变换的交叉项只出现在1fW和2fW同时不为零的),(a处，也就是两者相互交叠的区域中。性质 5 连续小波变换的等效频域形式小波变换dxaxxfaaWRf)(1),(的等效频域表示为RideafaaWf)(?)(?2),(证明 由傅立叶变换与函数卷积的性质?)(?)(?)()(fxxf所以?)(?)(? )()(fxxf从而? axxfa*)(1=)(?)(?afa可见RideafaaWf)(?)(?2),(小波

14、时频窗口性质名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换7式（ 1.19）定义的变换与短时傅立叶变换相比较，小波变换采用尺度的伸缩和平移因子平移。同时引入时间平移与频率伸缩，能保证能够建立起具有时间、频率同时局部化的窗口函数。为了研究小波的时频局部化性质，首先讨论如下一个基本性质。对于函数 f (x)的傅里叶变换)(?f而言，满足：?i?()e()f xf?()|()xfaf aa（1.30）现在讨论

15、由式（ 1.19）定义的窗口函数的时频局部化性质。从时域角度来看，当a,b(t)作为窗口函数时，其中心t0与窗口宽a,b分别为,12220,220,2,21() |( )| d |1|( ) | d|aaaaattttttttRR-式（1.19）可以表示成卷积()fa的形式。 而从频率的角度来看，利用 Parseval等式，又有,1?( , )()()d2aW f af，因此利用1i2,?| e()aaa，得到频率窗口的中心0与宽度,a分别为,20,2,212220,21?()d?1?()()d?aaaaaRR（1.31）记 a=1， =0， 此时,( )( )att ， 而相应的时、 频窗口

16、参数分别记为*0t，以及*0，?。于是，可以建立下面等式：,00?0?0*,*,aatataaa（1.32）下面讨论式（1.32）的证明。由于两个等式的证明相似，因此，为节省篇幅，只证明式（1.32）中第一行的等式。事实上，直接计算有2200221|() | d() |( ) | d*1|( ) | d|() | dtttaxxxaatattxxtaaRRRR类似得到,a ba。由式（ 1.32）建立时 -频窗口满足：,0000,a ba ba ba btt名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - -

17、 - - - - 第 7 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换8?000*,*,ataataa?0a此时，窗口的时间宽度为2a，频率宽度为2?/ a，因此其面积为4?，与a和的选取无关。窗口的特点：当需要检测高频分量时，减少a的值，此时时间窗口自动变窄，而频率窗口自动变宽，此时为一时宽窄而频宽大的高频窗；而在检测低频分量时，增加a值，时间窗口自动变宽，频率窗口自动变窄，此时为一时宽大而频宽窄的低频窗。t2a?2/ a?2/ a?2/ a?2/ a2a2a?2/ a?2/ a2a2a2a小波时频窗口面积?4基本小波 (x)应该具有快速衰减性质，其振幅为正负

18、相间的震荡形式。特别地，将式（1.19）所定义的变换称之为小波变换，而相应的函数(t)称之为小波函数。2.3 小波变换的反演及对基本小波的要求定理 3.1 小波变换的内积定理（Moyal 定理）以基本小波)(x分别对)(1xf和)(2xf作小波变换，设)(1xf的连续小波变换为：)(),(),(11xxfaWaf)(2xf的连续小波变换为：)(),(),(22xxfaWaf其中axaxa1)(，则有)(),(),(),(2121xfxfcaWaWff（*）式中dc02)(?。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精

19、心整理 - - - - - - - 第 8 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换9（* ）式可以写成更明确的形式，左边的内积石对, a作双重积分，而且由于式中a以倒数形式出现，所以微分为daa21，从而可将（ *）式写为更具体的形式。dxxfxfcdadxfxxxfaRRaa)()()(),()(),(1212102（* ）左边第二个内积中的两个因子次序对调反映了取共轭。证明首先根据 Parseval等式的广义形式dgfxgxfR)( ?)(?21)(),(得到)(),(1xxfa=dfRa)(?)(?211（i）)(),(2xfxa=dfRa)(?)(

20、?212（ii）又axaxa1)(的傅立叶变换为)(?a=ieaa)(?（iii）)(?a=ieaa)(?（iv）将式（ iii） 、 （iv）代入（ i） 、 （ii） ，再把式（ i） （ii）代入（ * ）式左边，并由)(2)(Ride整理后得到左边 =dadaaffaR)(?)( ?)(?)(?121210=dffdaaaR)(?)(?)(?212102设02)( ?daaa积分存在，即0202)( ?)(?daaadaaa=02)( ?d=c则上式最后成为左边 =Rdffc)(?)(?2121=)(),(21xfxfc=右边名师资料总结 - - -精品资料欢迎下载 - - - - -

21、 - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换10由证明可见，内积定理的成立以c存在为条件，存在的条件可以明确地表示成dc02)(?小波与常用变换的区别在于没有固定的核函数，但也不是任何函数都可以用作小波变换的基本小波)(x。任何变换只有存在反变换才有意义，对于小波变换而言，只有基本小波满足“容许条件” ，反变换才存在。2.4.1 容许条件定理 2.2 （小波反演公式）当dc02)(?时，才能由小波变换),(aWf反演原函数)(xf，此时RafdadxaWacx

22、f)(),(11)(02=RfdadaxaaWac)(1),(1102dc02)(?，便是对基本小波)(x提出的容许条件。证明利用内积定理来证明令)()(1xfxf，)()(2yxxf由)(x的采样性)()(),(yfyxxf则Raadadyxxxxfayfc)(),()(),(1)(02=RfdadayaaWa)(1),(102即RfdadaxaaWacxf)(1),(11)(02此式成立条件为dc02)(?，也就是内积定理存在的条件。注意连续小波变换中取共轭，而反演公式中，不取共轭。另外由容许条件dc02)(?，可以推出，基本小波的一个必要条件就是0)0( ?。事实上，名师资料总结 - -

23、 -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换11对于积分d02)(?，由于包含奇点0，只有当0)0( ?时积分才有可能存在。基于此，对于基本小波)(x，其在时域上的表现形式事实上由于Rixdxex)()(?Rdxx)()0(?=0 说明基本小波函数为上下交替的震荡波，其在整个实域上的积分为0。这也是小波的消失矩特性，小波的消失矩越高，其局部性能越好。2.4.2 能量的比例性由 Moyal 公式可以引出小波Parseva

24、l 等式，即小波变换幅度平方的积分和函数的能量成正比。RfdadaWa202),(1=dxxfcR2)(证明令)()()(21xfxfxf，得)(),(),(xxfaWaf)(),(),(xfxaWafdxxfxfxfR2)()(),(代入内积公式可得RfdadaWa202),(1=dxxfcR2)(2.4 消失矩条件满足容许条件的)(x可以用作基本小波，但是实际上要求更高，必须对)(x施加“正规性条件”。以便?( )在频域上表现更好的局部性，也就是要求( , )fWa随a的减少而迅速减少。 在小波理论中， 正规性条件也被称为消失矩条件。这要求)(x的前n阶原点矩等于 0，且n越高局部性能越好

25、，即( )0pRxx dx，1,2,.,pn（i）从上式可以看出，当0p时也就是容许条件。说明基本小波函数是一个在x轴上下震动，且面积相等的函数。下面来分析一下消失矩条件的等价频域性质。定理 2.1 小波函数的消失矩越高，则小波的在频域上的局部分析能力越强（( , )fWa随a名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换12的减少而迅速减少） 。基本小波函数)(x具有n阶消失矩， 等价于?( )在0

26、有n阶零点，其中一阶零点是容许条件。即10?()()n，0? (0)0（ii）式（i）和式（ ii）就是消失矩条件。证明先证（ i） ，由于dxaxxfaaWRf)(1),(将( )f x在x处泰勒展开有2(1)(2)1()( , )( )( )( ).1!2!fRttxWafffdxaa只取上式前n项，得()01()( )!pnpRpxxfdxpaa=1()01( )!ppnpRpayfy dypa,xydxadya也就是132( )20( , )( )()!ppnnpfRpayWafy dyo ap=1313(1)( )222201( )( )( ).()1!nnnnfffM aM aM

27、ao an式中( )ppRMxx dx由容许条件00M，如果进一步有0pM，1,2,.,pn，则( , )fWa随a的减少，而不低于32na的速度减少。下证（ i）等价于（ ii）由于?( )是)(x的傅立叶变换，所以有?()( )ixRx edx由此有?()()( )ixRdixx edxd，故，0?( )|()( )Rdixx dxd222?( )()( )ixRdixx edxd，故，22202?()|()( )Rdixx dxd. 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - -

28、第 12 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换13?()()( )nnixnRdixx edxd，故，0?()|()( )nnnnRdixx dxd可见要求( )0pRxx dx，1,2,.,pn等效于要求0?()|0ppdd，1,2,.,pn。也就是要求10?()()n，0? (0)0这项要求消除了( )f x的多项式展开式中()pxpn各项在小波变换中的贡献，以便突出信号的高阶起伏和高阶导数种可能存在的奇点。即小波变换将反映信号中的高阶变化。基于消失矩条件的小波构造由Daubechies 提出，是小波构造的最重要的方法。在后面的小波构造的章节，我们将

29、详细介绍。以二维数据（图像）为例，高的消失矩可使矩阵变得更加稀疏，这样使得数据非常利于压缩，因此在信号压缩处理中，我们往往选择消失矩较大的小波基比如CDF9/7。在信号检测的应用中为了能够有效地检测奇异点，小波基的消失矩也必须具有足够的阶数，它与Lipschitz 指数密切相关。然而，突变信号的Lipschitz 指数一般在0 到 1 内，因此为了分析突变信号，消失矩的阶数也不能太高，过高的阶数将使分析结果模糊，也增加计算量。因此选择合适的消失矩的小波变换来进行信号处理，是非常重要的因素。2.5 重建核与重建核方程重建核方程是小波变换的另一个重要性质，它说明小波变换的冗余性，即a平面上各点的小

30、波变换值是相关的。00(,)a处的小波变换值00(,)fWa可以表示成半平面（aR，R）上其他各处小波变换值的贡献。0000201(,)( , )(, , )ffRWaWaKaad daa（）式中0000,1(, , )aaRKaaxx dxc=000111Rxxdxcaaaa=001,aaxxc名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换14K是小波00,aaxx的内积，它反映两者的相关程度，称为

31、重建核；（）式称为重建核方程。证明由小波变换的定义式和反变换式，有00(,)( )( )faRWaf xx dxRafdadxaWacxf)(),(11)(02将下式代入上一式得到0 0002011(,)( , )( )( )ffaaRRWaWaxx dx d daca则得到0000201(,)( , )(, , )ffRWaWaKaad daa记0000,1(, , )aaRKaaxx dxc=0 01,aaxxc重建核定理还可以写成0000201(,)( , ),ffRaWaWaWd daaaa实际上就是要证明00( )( )aaRxx dx=00,aWaa证明0 00001( )( )a

32、aRRxxxx dxdxaaaa=000()Rayaydyaa令xya=000()Ryaaydyaaa=00,aWaa名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换1500(,)fWa00(,)fWa连续小波变换的冗余性说明1、 容许条件是小波的必要条件，只有满足容许条件的的函数( )f x才可以充当基本小波，重建核方程是二元函数( , )F a能够成为某个函数小波变换( , )fWa的必要条件。 也

33、就是说不是所有的( , )F a都可以作为( , )fWa，( , )F a必须满足重建核方程。2、 重建核00(, , )Kaa反映的是0 0aaxx和的相关性。当0aa，0时00(, , )Kaa达到最大。如果( , )a偏离（00,a）时K衰减较快，两者的相关区域便愈小。如果0000(, , )(,)Kaaaa，此时a半平面内各点的小波变换值互不相关，小波变换所包含的信息才没有冗余，这就要求不同尺度及不同位移的小波相互正交。由于,a是连续变量，这种正交性是不存在的。因此，当( )f x被变换成( , )fWa后信息是有冗余的。连续小波变换可以直观的看，用一个平面去表示一条曲线。从重建核方

34、程，很容易看出，一个一维的信号变换到二维以后，存在的信息冗余是非常大的。因此连续变化的尺度因子(0,)a是不必要的。连续小波变换怎样来选取一部分尺度因子(0,)a（其中为指标集），使得函数( )f x能重构。与自然数集同构时，a成为一个点列， 这样的方式是不是能实现信号的重构呢？名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换16尺度因子离散化方案中，a的元素如何选取也是我们需要解决的问题。在实际应用中

35、，信号进行离散采样，其小波变换中的位移因子可以对应于信号的历史采样点的采样位置。当 a 取所有大于0 的实数值时，频率轴0,+将会被完全覆盖，同时也不难发现，这种覆盖存在严重的交叉覆盖因而产生覆盖的冗余。为了在完全覆盖的前提下避免冗余现象的出现。如何通过选取频率值a 的离散值，仍能实现原信号的精确重构这是小波理论的重要研究问题。设?为频域上的一个固定的窗口宽度，做如下的频域轴分解：?0,)2jj1?,2)j02345671a0kb1242.6 连续小波的计算计算连续小波变换的基本方法是数值近似积分法，这种方法包括矩形法和梯形法。设0a，且( ) t为实函数小波1.矩形数值积分法在给定的尺度下，

36、 对待分析信号( )f t和分析小波( ) t按照StnT，skT进行采样，其中SnT为采样间隔，则采用矩形数值积分法，小波变换可以近似计算如下：()( ,)()ssfssnTnk TWa kTf nTaa简记为( , )( )fnTnkWa kf naa式中，sTT，( )f n表示信号( )f t的采样序列，( )n表示分析信号( ) t的采样序列。对每个给定的a值，依次求出不同k 值下的一组小波系数。在实际计算时，a也需要名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页

37、，共 35 页 - - - - - - - - - 第 1 章连续小波变换17取离散值，通常取2ja。梯形数值积分法同样的，对待分析信号( )f t和分析小波( ) t按照StnT，skT进行采样，其中SnT为采样间隔。为使计算的数值更加精确，可采用梯形法做近似积分。其计算公式如下：1( ,)()(1)2sfsssssnTnknkWa kTf nTTfnTTaaa式中，sTT，简记为1( , )( )(1)2fnTnknkWa kf nTfnTTaaa这说明，对于任何一个a值，可以通过计算离散信号序列与小波函数采样序列的卷积来计算小波系数。应用小波变换分析的步骤为“（1）计算每一尺度下各个离散

38、位置的小波变换，然后将结果显示。（2）从频域上来看，用不同的a值做处理， 相当于用不同中心频率的带通滤波器做处理；从时域上来看，表现为信号在各局部时段的处理结果。（3）两者结合起来就能将一定范围内的信号特征加以凸显。在 Matlab 小波工具箱中，用cwt()函数计算连续小波变换。例如，对Matlab 中提供的某一信号进行连续小波变换，尺度a分别为 12.12，10.24，16.38，1.2，2，4，5，8，10，小波函数用db3，则连续小波变换的系数的灰度表示图，在图形中，小波系数的大小( , )fWa使用灰度的深浅表示，颜色越深，则表示变换以后的系数越大。离散小波变换的尺度和位移参数只能在

39、二进栅格上离散取值，因此只能得到二进尺度和二进栅格上的小波系数，使得它在许多信号分析场合显得粗糙。利用连续小波变换可以计算任意尺度上和相位上的小波系数，比离散小波变换具有更大的灵活性。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换18三维显示信号的连续小波变换名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - -

40、- - - 第 18 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换192.7 常用小波介绍与标准傅里叶变换相比，小波分析中所用到的小波函数具有不唯一性，即小波函数( )x具有多样性。 但小波分析在工程应用中的一个十分重要的问题是最优小波基的选择问题，这是因为用不同的小波基分析同一个问题会产生不同的结果。目前，主要是通过用小波分析方法处理信号的结果与理论结果的误差来判定小波基的好坏，并由此选定小波基。根据不同的标准，小波函数具有不同的类型，这些标准通常有：(1)、和的支撑长度。即当时间或频率趋向无穷大时，、和从一个有限值收敛到0 的速度。(2)对称性。它在图像处理

41、中对于避免移相是非常有用的。(3)和如果存在的情况下)的消失矩阶数。它对于压缩是非常有用的。(4)正则性。它对信号或图像的重构获得较好的平滑效果是非常有用的。但在众多小波基函数(也称核函数 )的家族中，有一些小波函数被实践证明是非常有用的。我们可以通过waveinfo 函数获得工具箱中的小波函数的主要性质，小波函数和尺度函数可以通过 wavefun 函数计算，滤波器可以通过wfilters 函数产生。我们主要介绍一下MATLAB中常用到的小波函数。1 Haar 小波Haar 函数是在小波分析中最早用到的一个具有紧支撑的正交小波函数，同时也是最简单的一个函数，它是非连续的，类似一个阶梯函数。Ha

42、ar 函数与下面将要介绍的db1小波函数是一样的。Haar 函数的定义为101/ 211120Hxx其它尺度函数为101( )0 xx其它在 MATLAB中，可以输入命令waveinfo( haar )获得 Haar 函数的一些主要性质，如图 1.3所示。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换202 Daubechies(dbN) 小波系Daubechies 函数是由世界著名的小波分析学者In

43、rid Daubechies 构造的小波函数，除了db1(即 haar 小波 )外，其他小波没有明确的表达式，但转换函数h 的平方模是很明确的。dbN函数是紧支撑标准正交小波，它的出现使离散小波分析成为可能。假设110( )NNkkkkP yCy，其中1NkkC为二项式的系数，则有222()(cos)(sin ()22NHP其中21i01( )e2NkkkHh。22()()1HH计算出2()H，利用 Riesiz引理，求出()H，最后我们可以利用1( )()2kkxH当然无限乘积可能没有解析表示式解，这也是 Daubechies 小波没有解析表达式的原因。具体的构造过程在后面的小波构造理论中会

44、详细介绍。Daubechies 小波特点（1）小波函数和尺度函数的有效支撑长度为2N1， 小波函数 y 的消失矩阶数为N。（2）大多数 dbN 不具有对称性，对于有些小波函数，不对称性是非常明显的。（3）正则性随着序号N 的增加而增加。（4）函数具有正交性。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 20 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换21在这里，画出db4 和 db8 小波的尺度函数、小波函数、分解滤波器和重构滤波器的图形，如下图所

45、示。Daubechies 小波函数提供了比Haar 组更有效的分析和综合。Daubechies 系中的小波基记为dbN， N 为序号，且N 1，2， ，10。在 MATLAB中， 可以输入命令waveinfo(db )获得 Daubechies 函数的一些主要性质。3Biorthogonal(biorNr.Nd)小波系Biorthogonal 函数系的主要特性体现在具有线性相位性，它主要应用在信号与图像的重构中，通常的用法是采用一个函数进行分解，用另外一个小波函数进行重构。众所周知，如果使用同一个滤波器进行分解和重构，对称性和重构的精确性将成为一对矛盾，而采用名师资料总结 - - -精品资料欢

46、迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 21 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换22两个函数，将有效地解决这个问题。设函数用于信号分解，而函数y 用于信号重构，则分解和重构的关系式为_,( )( ) dj kj kCs xxx,j kj kj kSC另外具有双正交性,0,0,( )( )d0( )( )d0j kj kkkxxxxxx这样，利用函数,( )j kx和,( )j kx的特性，在信号分解时可以获得一些很好的分解性质(如振动、零力矩 )，而利用,( )j k

47、x和,( )j kx的特性，在信号重构时又可获得一些很好的重构性质 (如正则性 )。Biorthogonal 函数系通常表示成biorNr.Nd 的形式：Nr Nd 1 1，3，5 2 2，4，6，8 3 1，3，5，7，9 4 4 5 5 6 8 其中， r 表示重构 (Reconstruction)；d 表示分解 (Decomposition) 。在这里，我们画出 bior2.4 和 bior4.4 小波 (分别用于分解与重构)的尺度函数、 小波函数、分解滤波器和重构滤波器的图形，如图1.5 所示。在 MATLAB中，可输入命令waveinfo( bior )获得该函数的主要性质。名师资料

48、总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 22 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换234Coiflet(coifN)小波系Coiflet 函数也是由Daubechies 构造的一个小波函数，它具有coifN(N1，2，3，4， 5)这一系列。 Coiflet 具有比 dbN 更好的对称性。从支撑长度的角度看，coifN 具有和 db3N 和sym3N 相同的支撑长度；从消失矩的数目来看，coifN 具有和 db2N 和 sym2N 相同的消失矩

49、数目。在这里，我们画出coif3 和 coif5 小波的尺度函数、小波函数、分解滤波器和重构滤波器的图形，如下图所示。在 MATLAB中，可输入命令waveinfo( coif )获得该函数的主要性质。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 23 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换245SymletsA(symN) 小波系Symlets 函数系是由Daubechies 提出的近似对称的小波函数，它是对db 函数的一种改进。Symlet

50、s 函数系通常表示为symN(N2，3， ，8)的形式。在这里，我们画出sym4 和 sym8 小波的尺度函数、小波函数、分解滤波器和重构滤波器的图形，如图1.7 所示。在 MATLAB中，可输入waveinfo( sym )获得该函数的主要性质。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 24 页，共 35 页 - - - - - - - - - 第 1 章连续小波变换256Morlet(morl) 小波Morlet 函数定义为2/2( )ecos(5 )xxCx它的尺度函数不