小波分析理论PPT课件.ppt

单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,.,*,单击此处编辑母版标题样式,*,第,1,章小波分析的基本理论,第,1,章小波分析的基本理论,1.1,傅里叶变换到小波分析,1.2,常用小波函数介绍,1.3,连续小波变换,1.4,离散小波变换,1.5,矢量小波变换,1.6,多分辨分析与,Mallat,算法,1.7,提升小波变换,1.8,小波包分析,1,.,小波分析属于时频分析的一种。传统的信号分析是建立在傅里叶,(Fourier),变换的基础上的，但是，傅里叶分析使用的是一种全局的变换，即要么完全在时域，要么完全在频域，它无法表述信号的时频局域性质，而时频局域性质恰恰是非平稳信号最根本和最关键的性质。为了分析和处理非平稳信号，人们对傅里叶分析进行了推广乃至根本性的革命，提出并发展了幌盗行碌男藕欧治隼砺郏憾淌备道镆侗浠弧逼捣治觥,Gabor,变换、小波变换、,RandonWigner,变换、分数阶傅里叶变换、线性调频小波变换、循环统计量理论和调幅调频信号分析等。,2,.,其中，短时傅里叶变换和小波变换也是因传统的傅里叶变换不能够满足信号处理的要求而产生的。短时傅里叶变换分析的基本思想是：假定非平稳信号在分析窗函数,g,(,t,),的一个短时间间隔内是平稳,(,伪平稳,),的，并移动分析窗函数，使,f,(,t,),g,(,t,t,),在不同的有限时间宽度内是平稳信号，从而计算出各个不同时刻的功率谱。但从本质上讲，短时傅里叶变换是一种单一分辨率的信号分析方法,(,因为它使用一个固定的短时窗函数,),，在信号分析上还存在着不可逾越的缺陷。,3,.,小波变换是一种信号的时间,尺度,(,时间,频率,),分析方法，它具有多分辨率分析,(Multi,-,resolutionAnalysis),的特点，而且在时频两域都具有表征信号局部特征的能力，是一种窗口大小固定不变，但其形状可改变，时间窗和频率窗都可以改变的时频局部化分析方法。即,在低频部分具有较高的频率分辨率和较低的时间分辨率，在高频部分具有较高的时间分辨率和较低的频率分辨率,，很适合于探测正常信号中夹带的瞬态反常现象并展示其成分，所以被誉为分析信号的显微镜。,4,.,1.1.1,傅里叶变换傅里叶变换是众多科学领域,(,特别是信号处理、图像处理、量子物理等,),里的重要的应用工具之一。从实用的观点看，当人们考虑傅里叶分析的时候，通常是指,(,积分,),傅里叶变换和傅里叶级数。,1.1,傅里叶变换到小波分析,5,.,定义,1.1,函数,f,(,t,),L,1,(,R,),的连续傅里叶变换定义为,(1.1)F(,w,),的傅里叶逆变换定义为,(1.2),6,.,为了计算傅里叶变换，需要用数值积分，即取,f,(,t,),在,R,上的离散点上的值来计算这个积分。在实际应用中，我们希望在计算机上实现信号的频谱分析及其他方面的处理工作，对信号的要求是：在时域和频域应是离散的，且都应是有限长的。下面给出离散傅里叶变换,(Discrete Fourier Transform,，,DFT),的定义。,7,.,定义,1.2,给定实的或复的离散时间序列,f,0,，,f,1,，,，,f,N,1,，设该序列绝对可积，即满足，称,(1.3),为序列,f,n,的离散傅里叶变换，称,8,.,(1.4),为序列,X,(,k,),的离散傅里叶逆变换,(IDFT),。在式,(1.4),中，,n,相当于对时间域的离散化，,k,相当于频率域的离散化，且它们都是以,N,点为周期的。离散傅里叶变换序列,X,(,k,),是以,2,p,为周期的，且具有共轭对称性。,9,.,若,f,(,t,),是实轴上以,2,p,为周期的函数，即,f,(,t,)L,2,(0,，,2,p,),，则,f,(,t,),可以表示成傅里叶级数的形式，即,(1.5),傅里叶变换是时域到频域互相转化的工具，从物理意义上讲，傅里叶变换的实质是把,f,(,t,),这个波形分解成许多不同频率的正弦波的叠加和。这样我们就可将对原函数,f,(,t,),的研究转化为对其权系数，即其傅里叶变换,F,(,w,),的研究。从傅里叶变换中可以看出，这些标准基是由正弦波及其高次谐波组成的，因此它在频域内是局部化的。,10,.,在进行傅里叶变换时，如果能合理运用它的有关性质，运算将很方便。下面列出了傅里叶变换的一些常用性质。,11,.,1.,线性性质设,F,1,(,w,),和,F,2,(,w,),分别为,f,1,(t),和,f,2,(t),的傅里叶变换，,a,和,b,为常数，则有,a,f,1,(t),b,f,2,(t),a,F,1,(,w,),b,F,2,(,w,),(1.6),这个性质表明，函数线性组合的傅里叶变换等于各函数傅里叶变换的线性组合。傅里叶逆变换亦具有类似的性质。,12,.,2.,位移性质设,F,(,w,),为函数,f,(t),的傅里叶变换，则有,(1.7),该性质表明，时间函数,f,(,t,),沿,t,轴向左或向右位移,t,0,的傅里叶变换等于,f,(,t,),的傅里叶变换乘以因子或。傅里叶逆变换亦具有类似的位移性质。,13,.,3.,微分性质设,F,(,w,),为函数,f,(,t,),的傅里叶变换，,f,(,t,),表示函数,f,(,t,),的微分，则有,(1.8),该性质表明，一个函数的导数的傅里叶变换等于这个函数的傅里叶变换乘以因子,j,w,。由该性质可以导出一般的微分公式：,14,.,4.,积分性质设,F,(,w,),为函数,f,(,t,),的傅里叶变换，如果当,t,时，则有,(1.9),15,.,5.,乘积定理设,F,1,(,w,),和,F,2,(,w,),分别为,f,1,(,t,),和,f,2,(,t,),的傅里叶变换，则有,(1.10),其中，,f,1,(,t,),和,f,2,(,t,),为,t,的实函数；和分别为,F,1,(,w,),和,F,2,(,w,),的共轭函数。,16,.,6.,能量积分设,F,(,w,),为函数,f,(,t,),的傅里叶变换，则有,(1.11),该式又称为巴塞瓦,(Parseval),等式。,17,.,例,1-1,在某工程实际应用中，有一信号的主要频率成分是由,50 Hz,和,300 Hz,的正弦信号组成，该信号被一白噪声污染，现对该信号进行采样，采样频率为,1000 Hz,。通过傅里叶变换对其频率成分进行分析。解该问题实质上是利用傅里叶变换对信号进行频域分析，其,MATLAB,程序如下：,t,0,：,0.001,：,1.3,；,%,时间间隔为,0.001,说明采样频率为,1000 Hz,x,sin(2*pi*50*t),sin(2*pi*300*t),；,%,产生主要频率为,50 Hz,和,300 Hz,的信号,18,.,f,x,3.5*randn(1,，,length(t),；,%,在信号中加入白噪声,subplot(321),；,plot(f),；,%,画出原始信号的波形图,Ylabel(,幅值,),；,Xlabel(,时间,),；,title(,原始信号,),；,y,fft(f,，,1024),；,%,对原始信号进行离散傅里叶变换，参加,DFT,的采样点个数为,1024,p,y.*conj(y)/1024,；,%,计算功率谱密度,19,.,ff,1000*(0,：,511)/1024,；,%,计算变换后不同点所对应的频率值,subplot(322),；,plot(ff,，,p(1,：,512),；,%,画出信号的频谱图,Ylabel(,功率谱密度,),；,Xlabel(,频率,),；,title(,信号功率谱图,),；程序输出结果如图,1.1,所示。,20,.,图,1.1,21,.,从图,1.1(a),中我们看不出任何频域的性质，但从信号的功率谱图,(,图,1.1(b),中，我们可以明显地看出该信号是由频率为,50 Hz,和,300 Hz,的正弦信号和频率分布广泛的白噪声信号组成的，也可以明显地看出信号的频率特性。,22,.,虽然傅里叶变换能够将信号的时域特征和频域特征联系起来，能分别从信号的时域和频域观察，但不能把二者有机地结合起来。这是因为信号的时域波形中不包含任何频域信息，而其傅里叶谱是信号的统计特性。从其表达式中也可以看出，它是整个时间域内的积分，没有局部化分析信号的功能，完全不具备时域信息，也就是说，对于傅里叶谱中的某一频率，不能够知道这个频率是在什么时候产生的。这样在信号分析中就面临一对最基本的矛盾：时域和频域的局部化矛盾。,23,.,在实际的信号处理过程中，尤其是对非平稳信号的处理中，信号在任一时刻附近的频域特征都很重要。如柴油机缸盖表面的振动信号就是由撞击或冲击产生的，是一瞬变信号，单从时域或频域上来分析是不够的。这就促使人们去寻找一种新方法，能将时域和频域结合起来描述观察信号的时频联合特征，构成信号的时频谱。这就是所谓的时频分析法，亦称为时频局部化方法。,24,.,1.1.2,短时傅里叶变换由于标准傅里叶变换只在频域里有局部分析的能力，而在时域里不存在局部分析的能力，因此,Dennis Gabor,于,1946,年引入了短时傅里叶变换,(Short,-,time Fourier Transform),。短时傅里叶变换的基本思想是：把信号划分成许多小的时间间隔，用傅里叶变换分析每一个时间间隔，以便确定该时间间隔存在的频率。其表达式为,(1.12),25,.,其中，“*”表示复共轭；,g,(,t,),为有紧支集的函数；,f,(,t,),为被分析的信号。在这个变换中，,e,j,w,t,起着频限的作用，,g,(,t,),起着时限的作用,。随着时间,t,的变化，,g,(,t,),所确定的“时间窗”在,t,轴上移动，使,f,(,t,)“,逐渐”进行分析。因此,g,(,t,),往往被称为窗口函数，,S,(,w,，,t,),大致反映了时刻为,t,、频率为,w,时,f,(,t,),的“信号成分”的相对含量。这样，信号在窗函数上的展开就可以表示为在,t,d,，,t,d,、,w,e,，,w,e,这一区域内的状态，并把这一区域称为窗口，,d,和,e,分别称为窗口的时宽和频宽，表示了时频分析中的分辨率，窗宽越小则分辨率就越高。,26,.,很显然希望,d,和,e,都非常小，以便有更好的时频分析效果，但,海森堡,(Heisenberg),测不准原理,(Uncertainty Principle),指出，,d,和,e,是互相制约的，两者不可能同时都任意小,(,事实上，,，且仅当为,高斯函数,时，等号成立,),，变换如图,1.2,所示。,27,.,图,1.2,28,.,由此可见，短时傅里叶,(STFT),虽然在一定程度上克服了标准傅里叶变换不具有局部分析能力的缺陷，但它也存在着自身不可克服的缺陷，即当窗函数,g,(,t,),确定后，矩形窗口的形状就确定了，,t,、,w,只能改变窗口在相平面上的位置，而不能改变窗口的形状。可以说,STFT,实质上是具有单一分辨率的分析，若要改变分辨率，则必须重新选择窗函数,g,(,t,),。,因此，,STFT,用来分析平稳信号犹可，但对非平稳信号，在信号波形变化剧烈的时刻，主频是高频，要求有较高的时间分辨率,(,即,d,要小,),，而波形变化比较平缓的时刻，主频是低频，则要求比较高的频率分辨率,(,即,e,要小,),，而短时傅里叶不能兼顾两者。,29,.,1.1.3,小波分析小波分析方法是一种窗口大小,(,即窗口面积,),固定但其形状可改变，时间窗和频率窗都可改变的时频局部化分析方法。即在,低频部分具有较高的频率分辨率和较低的时间分辨率，在高频部分具有较高的时间分辨率和较低的频率分辨率,，所以被誉为数学显微镜。正是这种特性，使小波变换具有对信号的自适应性。,信号长度越长，频率分辨率越好,30,.,小波分析被看成调和分析这一数学领域半个世纪以来的工作结晶，已经广泛地应用于信号处理、图像处理、量子场论、地震勘探、语音识别与合成、音乐、雷达、,CT,成像、彩色复印、流体湍流、天体识别、机器视觉、机械故障诊断与监控、分形以及数字电视等科技领域。原则上讲，传统上使用傅里叶分析的地方，都可以用小波分析取代。小波分析优于傅里叶变换的地方是，它在时域和频域同时具有良好的局部化性质。,31,.,设,y,(,t,),L,2,(,R,)(,L,2,(,R,),表示平方可积的实数空间，即能量有限的信号空间,),，其傅里叶变换为,Y,(,w,),。当,Y,(,w,),满足允许条件,(Admissible Condition),：,(1.13),时，我们称,y,(,t,),为一个基本小波或母小波,(Mother Wavelet),。将母函数,y,(,t,),经伸缩和平移后，就可以得到一个小波序列。,32,.,对于连续的情况，小波序列为,(1.14),其中，,a,为伸缩因子；,b,为平移因子。对于离散的情况，小波序列为,(1.15),33,.,对于任意的函数,f,(,t,),L,2,(,R,),的连续小波变换为,(1.16),其逆变换为,(1.17),34,.,小波变换的时频窗口特性与短时傅里叶的时频窗口不一样。其窗口形状为两个矩形,b,a,D,y,，,b,a,D,y,(,w,0,D,Y,)/,a,，,(,w,0,D,Y,)/,a,，窗口中心为,(b,，,w,0,/,a,),，时窗和频窗宽分别为,a,D,y,和,D,Y,/,a,。其中，,b,仅仅影响窗口在相平面时间轴上的位置，而,a,不仅影响窗口在频率轴上的位置，也影响窗口的形状。这样小波变换对不同的频率在时域上的取样步长是调节性的：在低频时，小波变换的时间分辨率较低，而频率分辨率较高；在高频时，小波变换的时间分辨率较高，而频率分辨率较低，这正符合低频信号变化缓慢而高频信号变化迅速的特点。,35,.,这便是它优于经典的傅里叶变换与短时傅里叶变换的地方。从总体上来说，小波变换比短时傅里叶变换具有更好的时频窗口特性。,36,.,1.1.4,小波分析与傅里叶变换的比较小波分析是傅里叶分析思想方法的发展与延拓，它自产生以来，就一直与傅里叶分析密切相关，它的存在性证明，小波基的构造以及结果分析都依赖于傅里叶分析，二者是相辅相成的。两者相比较主要有以下不同点。,(1),傅里叶变换的实质是把能量有限信号,f,(,t,),分解到以,e,j,w,t,为正交基的空间上去；小波变换的实质是把能量有限信号,f,(,t,),分解到,W,j,(,j,1,，,2,，,，,J,),和,V,j,所构成的空间上去。,37,.,(2),傅里叶变换用到的基本函数只有,sin(,w,t,),、,cos(,w,t,),、,exp(j,w,t,),，具有唯一性；小波分析用到的函数,(,即小波函数,),则具有不唯一性，同一个工程问题用不同的小波函数进行分析有时结果相差甚远。小波函数的选用是小波分析应用到实际中的一个难点问题,(,也是小波分析研究的一个热点问题,),，目前，往往是通过经验或不断的试验,(,对结果进行对照分析,),来选择小波函数。,38,.,(3),在频域中，傅里叶变换具有较好的局部化能力，特别是对于那些频率成分比较简单的确定性信号，傅里叶变换很容易把信号表示成各频率成分的叠加和的形式，如,sin(,w,1,t,),0.345sin(,w,2,t,),4.23cos(,w,3,t,),。但在时域中，傅里叶变换没有局部化能力，即无法从信号,f,(,t,),的傅里叶变换,F,(,w,),中看出,f,(,t,),在任一时间点附近的性态。事实上，,F,(,w,)d,w,是关于频率为,w,的谐波分量的振幅，在傅里叶展开式中，它是由,f,(,t,),的整体性态所决定的。,(4),在小波分析中，尺度,a,的值越大相当于傅里叶变换中,w,的值越小。,39,.,(5),在短时傅里叶变换中，变换系数,S,(,w,，,t,),主要依赖于信号在,t,d,，,t,d,片段中的情况，时间宽度是,2,d,(,因为,d,是由窗函数,g,(,t,),唯一确定的，所以,2,d,是一个定值,),。在小波变换中，变换系数,W,f,(a,，,b),主要依赖于信号在,b,a,D,y,，,b,a,D,y,片段中的情况，时间宽度是,2,a,D,y,，该时间宽度是随着尺度,a,的变化而变化的，所以小波变换具有时间局部分析能力。,40,.,(6),若用信号通过滤波器来解释，小波变换与短时傅里叶变换不同之处在于：对短时傅里叶变换来说，带通滤波器的带宽,D,f,与中心频率,f,无关；相反，小波变换带通滤波器的带宽,D,f,则正比于中心频率,f,，即,(C,为常数,),亦即滤波器有一个恒定的相对带宽，称之为等,Q,结构,(Q,为滤波器的品质因数，且有,),。,41,.,与标准傅里叶变换相比，小波分析中所用到的小波函数具有不唯一性，即小波函数,y,(,x,),具有多样性。但小波分析在工程应用中的一个十分重要的问题是最优小波基的选择问题，这是因为用不同的小波基分析同一个问题会产生不同的结果。目前，主要是通过用小波分析方法处理信号的结果与理论结果的误差来判定小波基的好坏，并由此选定小波基。,1.2,常用小波函数介绍,42,.,根据不同的标准，小波函数具有不同的类型，这些标准通常有：,(1),y,、,Y,、,f,和,F,的支撑长度。即当时间或频率趋向无穷大时，,y,、,Y,、,f,和,F,从一个有限值收敛到,0,的速度。,(2),对称性。它在图像处理中对于避免移相是非常有用的。,(3),y,和,f,(,如果存在的情况下,),的消失矩阶数。它对于压缩是非常有用的。,(4),正则性。它对信号或图像的重构获得较好的平滑效果是非常有用的。,43,.,但在众多小波基函数,(,也称核函数,),的家族中，有一些小波函数被实践证明是非常有用的。我们可以通过,waveinfo,函数获得工具箱中的小波函数的主要性质，小波函数,y,和尺度函数,f,可以通过,wavefun,函数计算，滤波器可以通过,wfilters,函数产生。在本节中，我们主要介绍一下,MATLAB,中常用到的小波函数。,44,.,1.2.1,Haar,小波,Haar,函数是在小波分析中最早用到的一个具有紧支撑的正交小波函数，同时也是最简单的一个函数，它是非连续的，类似一个阶梯函数。,Haar,函数与下面将要介绍的,db1,小波函数是一样的。,Haar,函数的定义为,(1.18),45,.,图,1.3,Harr,小波函数,46,.,尺度函数为,(1.19),在,MATLAB,中，可以输入命令,waveinfo(,haar,),获得,Haar,函数的一些主要性质，如图,1.3,所示。,47,.,1.2.2,Daubechies(dbN),小波系,Daubechies,函数是由世界著名的小波分析学者,Inrid Daubechies,构造的小波函数，除了,db1(,即,haar,小波,),外，其他小波没有明确的表达式，但转换函数,h,的平方模是很明确的。,db,N,函数是紧支撑标准正交小波，,它的出现使离散小波分析成为可能,。,假设,，其中，,为二项式的系数，则有,(1.20),48,.,其中，,小波函数,y,和尺度函数,f,的有效支撑长度为,2,N,1,，小波函数,y,的消失矩阶数为,N,。,大多数,db,N,不具有对称性，对于有些小波函数，不对称性是非常明显的。,正则性随着序号,N,的增加而增加。,函数具有正交性。,N,是小波的阶数,49,.,在这里，我们画出,db4,和,db8,小波的尺度函数、小波函数、分解滤波器和重构滤波器的图形，如图,1.4,所示。,Daubechies,小波函数提供了比,Haar,组更有效的分析和综合。,Daubechies,系中的小波基记为,db,N,，,N,为序号，且,N,1,，,2,，,，,10,。在,MATLAB,中，可以输入命令,waveinfo(,db,),获得,Daubechies,函数的一些主要性质。,50,.,图,1.4,51,.,1.2.3,Biorthogonal(bior,N,r.,N,d),小波系,Biorthogonal,函数系的主要特性体现在具有线性相位性，它主要应用在信号与图像的重构中，通常的用法是采用一个函数进行分解，用另外一个小波函数进行重构。众所周知，如果使用同一个滤波器进行分解和重构，对称性和重构的精确性将成为一对矛盾，而采用两个函数，将有效地解决这个问题。,设函数用于信号分解，而函数,y,用于信号重构，则分解和重构的关系式为,(1.21),52,.,(1.22),另外，,与,y,之间具有二元性,(1.23),这样，利用,函数的特性，在信号分解时可以获得一些很好的分解性质,(,如振动、零力矩,),，而利用,y,的特性，在信号重构时又可获得一些很好的重构性质,(,如正则性,),。,53,.,Biorthogonal,函数系通常表示成,bior,N,r.,N,d,的形式：,Nr,1,Nd,1,，,3,，,5,Nr,2Nd,2,，,4,，,6,，,8,Nr,3Nd,1,，,3,，,5,，,7,，,9,Nr,4Nd,4,Nr,5Nd,5,Nr,6Nd,8,其中，,r,表示重构,(Reconstruction),；,d,表示分解,(Decom,position),。,54,.,在这里，我们画出,bior2.4,和,bior4.4,小波,(,分别用于分解与重构,),的尺度函数、小波函数、分解滤波器和重构滤波器的图形，如图,1.5,所示。,在,MATLAB,中，可输入命令,waveinfo(,bior,),获得该函数的主要性质。,55,.,图,1.5,56,.,1.2.4,Coiflet(coif,N,),小波系,Coiflet,函数也是由,Daubechies,构造的一个小波函数，它具有,coif,N,(,N,1,，,2,，,3,，,4,，,5),这一系列。,Coiflet,具有比,db,N,更好的对称性。从支撑长度的角度看，,coif,N,具有和,db3,N,和,sym3,N,相同的支撑长度；从消失矩的数目来看，,coif,N,具有和,db2,N,和,sym2,N,相同的消失矩数目。在这里，我们画出,coif3,和,coif5,小波的尺度函数、小波函数、分解滤波器和重构滤波器的图形，如图,1.6,所示。,在,MATLAB,中，可输入命令,waveinfo(,coif,),获得该函数的主要性质。,57,.,图,1.6,58,.,1.2.5,SymletsA(sym,N,),小波系,Symlets,函数系是由,Daubechies,提出的近似对称的小波函数，它是对,db,函数的一种改进。,Symlets,函数系通常表示为,sym,N,(,N,2,，,3,，,，,8),的形式。在这里，我们画出,sym4,和,sym8,小波的尺度函数、小波函数、分解滤波器和重构滤波器的图形，如图,1.7,所示。,在,MATLAB,中，可输入,waveinfo(,sym,),获得该函数的主要性质。,59,.,图,1.7,60,.,图,1.8,61,.,1.2.6,Morlet(morl),小波,Morlet,函数定义为,(1.24),它的尺度函数不存在，且不具有正交性。在,MATLAB,中，可输入,waveinfo(,morl,),获得该函数的主要性质，如图,1.8,所示。,62,.,1.2.7,Mexican,Hat(mexh),小波,Mexican Hat,函数为,(1.25),63,.,图,1.9,64,.,它是,Gauss,函数的二阶导数，因为它像墨西哥帽的截面，所以有时称这个函数为墨西哥帽函数,(,如图,1.9,所示,),。墨西哥帽函数在时间域与频域都有很好的局部化，并且满足,(1.26),由于它的尺度函数不存在，因此分析不具有正交性。在,MATLAB,中，可输入,waveinfo(,mexh,),获得该函数的主要性质。,65,.,1.2.8,Meyer,函数,Meyer,小波的小波函数,y,和尺度函数,f,都是在频域中进行定义的，是具有紧支撑的正交小波。,(1.27),66,.,其中，,v,(,a,),为构造,Meyer,小波的辅助函数，且有,v,(,a,),a,4,(35,84a,70a,2,20a,3,),a,0,，,1,(1.28),(1.29),在,MATLAB,中，可输入,waveinfo(,meyr,),获得该函数的主要性质，如图,1.10,所示。,67,.,图,1.10,68,.,1.2.9,Battle-Lemarie,小波,Battle-Lemarie,小波在,MATLAB,工具箱中不存在，但它也是我们常用到的一个小波函数。它具有两种形式，一种具有确定的正交性，一种不具有确定的正交性。当,N,1,时，尺度函数是线性样条函数；当,N,2,时，尺度函数是具有有限支撑的,B-,样条函数。更一般的情况，对于一个,N,次,B,样条小波，尺度函数为,(1.30),69,.,当,N,为奇数时，,k,0,；当,N,为偶数时，,k,1,。式,(1.30),可以用来构造滤波器。,它的双尺度关系为,(1.31),当,N,为偶数时，,f,是对称的，,x,1/2,；当,N,为奇数时，,f,是反对称的，,x,0,。,70,.,1.2.10,其他一些实数小波简介下面介绍几个,MATLAB,工具箱中的实数小波函数。,1.RbioNr.Nd,小波,RbioNr.Nd,函数是,reverse,双正交小波。在,MATLAB,中，可输入,waveinfo(,rbio,),获得该函数的主要性质。,71,.,2.Gaus,小波,Gaus,小波是从高斯函数派生出来的,(1.32),其中，整数,p,是参数，由,p,的变化导出一系列的,f,(,p,),，它满足如下条件,f,(,p,),2,1,(1.33),在,MATLAB,中，可输入,waveinfo(,gaus,),获得该函数的主要性质。,72,.,3.Dmey,小波,Dmey,函数是,Meyer,函数的近似，它可以进行快速小波变换。在,MATLAB,中，可输入,waveinfo(,dmey,),获得该函数的主要性质。,73,.,1.2.11,其他一些复数小波简介下面所介绍的几个小波复数函数，均在,MATLAB,工具箱中。,1.Cgau,小波,Cgau,函数是复数形式的高斯小波，它是从复数的高斯函数中构造出来的，其表达式为,(1.34),其中，整数,p,是参数，由,p,的变化导出一系列的,f,(,p,),，它满足如下条件：,74,.,f,(,p,),2,1,(1.35),在,MATLAB,中，可输入,waveinfo(,cgau,),获得该函数的主要性质。,75,.,2.Cmor,小波,Cmor,是复数形式的,morlet,小波，其表达式为,(1.36),其中，,f,b,是带宽参数，,f,c,是小波中心频率。在,MATLAB,中，可输入,waveinfo(,cmor,),获得该函数的主要性质。,76,.,3.Fbsp,小波,Fbsp,是复频域,B,样条小波，表达式为,(1.37),其中，,m,是整数型参数；,f,b,是带宽参数；,f,c,是小波中心频率。在,MATLAB,中，可输入,waveinfo(,fbsp,),获得该函数的主要性质。,77,.,4.Shan,小波,Shan,函数是复数形式的,shannon,小波。在,B,样条频率小波中，令参数,m,1,，就得到了,Shan,小波，其表达式为,(1.38),其中，,f,b,是带宽参数；,f,c,是小波中心频率。在,MATLABA,中，可输入,waveinfo(,shan,),获得该函数的主要性质。下面，我们将,MATLAB,工具箱中,15,个小波,(,或小波系,),的一些主要的性质加以对比，见表,1-1,。,78,.,表,1-1,MATLAB,工具箱中,15,个小波,(,或小波系,),的主要性质,79,.,续表,80,.,1.3.1,一维连续小波变换定义,1.3,设,y,(,t,),L,2,(,R,),，其傅里叶变换为,Y,(,w,),，当,Y,(,w,),满足允许条件,(,完全重构条件或恒等分辨条件,),：,(1.39),1.3,连续小波变换,81,.,时，我们称,y,(t),为一个基本小波或母小波,(MotherWavelet),。将母函数,y,(t),经伸缩,(Dilation),和平移,(Translation),后得：,(1.40),称为一个小波序列。其中，,a,为伸缩因子；,b,为平移因子。,82,.,对于任意的函数,f,(,t,),L,2,(,R,),的连续小波变换为,(1.41),其重构公式,(,逆变换,),为,(1.42),83,.,由于基小波,y,(,t,),生成的小波,y,a,，,b,(,t,),在小波变换中对被分析的信号起着观测窗的作用，因此,y,(,t,),还应该满足一般函数的约束条件,(1.43),故,Y,(,w,),是一个连续函数。这意味着，为了满足完全重构条件,(1.39),，,Y,(,w,),在原点必须等于,0,，即,84,.,为了使信号重构的实现在数值上是稳定的，除了完全重构条件外，还要求小波,y,(,t,),的傅里叶变换满足下面的稳定性条件：,(1.44),式中，,0AB1),，公式,存在几种扩展的可能性，一种可能性是选择小波,f,(,t,),L,2,(,R,n,),使其为球对称，其傅里叶变换也同样为球对称，,(1.46),94,.,并且其相容性条件变为,(1.47),对所有的,f,，,g,L,2,(R,n,),有,(1.48),这里，,，,，其中,a,R,，,a,0,且,b,R,n,。,95,.,公式,(1.42),也可以写为,(1.49),如果选择的小波,y,不是球对称的，则可以用旋转进行同样的扩展与平移。例如，在二维时，可定义,(1.50),96,.,这里，,a,0,，,b,R,2,，,，相容条件变为,(1.51),该等式对应的重构公式为,(1.52),97,.,对于高于二维的情况，可以给出类似的结论。在,MATLAB,中，系统只给出了关于一维和二维的小波变换函数，有关这些函数用法的详细说明见,2.3,、,2.4,两节。,98,.,1.4.1,一维离散小波变换在实际运用中，尤其是在计算机上实现，连续小波必须加以离散化。因此，有必要讨论一下连续小波,y,a,，,b,(,t,),和连续小波变换,W,f,(,a,，,b,),的离散化。需要强调指出的是，这一离散化都是针对连续的尺度参数,a,和连续平移参数,b,的，而不是针对时间变量,t,的，这一点与我们以前习惯的时间离散化不同，希望引起注意。,1.4,离散小波变换,99,.,在连续小波中，考虑函数,(1.53),这里，,b,R,，,a,R,，且,a,0,，,y,是容许的，为方便起见，在离散化中，总限制,a,只取正值，这样相容性条件就变为,(1.54),100,.,通常，把连续小波变换中尺度参数,a,和平移参数,b,的离散化公式分别取作,，,，这里,j,Z,，扩展步长,a,0,1,是固定值，为方便起见，总是假定,a,0,1(,由于,m,可取正也可取负，所以这个假定无关紧要,),。所以对应的离散小波函数,y,j,，,k,(,t,),即可写作,(1.55),101,.,而离散化小波变换系数则可表示为,(1.56),其重构公式为,(1.57)c,是一个与信号无关的常数。,102,.,然而，怎样选择,a,0,和,b,0,才能够保证重构信号的精度呢？显然，网格点应尽可能密,(,即,a,0,和,b,0,尽可能小,),，因为如果网格点越稀疏，使用的小波函数,y,j,，,k,(t),和离散小波系数,c,j,，,k,就越少，信号重构的精确度也就会越低。在,MATLAB,中，离散小波变换可以用,dwt,、,dwt2,函数实现。,103,.,1.4.2,二维离散小波变换为了将一维离散小波变换推广到二维，只考虑尺度函数是可分离的情况，即,f,(x,，,y),f,(x),f,(y),(1.58),其中，,f,(x),是一维尺度函数，其相应的小波函数是,y,(x),，下列三个二维小波基是建立二维小波的基础：,y,1,(x,，,y),f,(x),y,(y),y,2,(x,，,y),y,(x),f,(y),y,3,(x,，,y),y,(x),y,(y),104,.,它们构成二维平方可积函数空间,L,2,(,R,2,),的正交归一基：,(1.59),二维离散小波分解的过程如下：,从一幅,NN,的图像,f,1,(,x,，,y,),开始，其中上标指示尺度,N,是,2,的幂。对于,j,0,，,2,j,2,0,1,尺度，也就是原图像的尺度。,j,值的每一次增大都使尺度加倍，而使分辨率减半。,105,.,在变换的每一层次，图像都被分解为四个,1/4,大小的图像，它们都是由原图与一个小波基图像的内积后，再经过在行和列方向进行,2,倍的间隔抽样而生成的。对于第一个层次,(,j,1),，可写成,后续的层次,(,j,1),，依次类推，形成如图,1.12,所示的形式。,106,.,图,1.12,107,.,若将内积改写成卷积形式，则有,108,.,因为尺度函数和小波函数都是可分离的，所以每个卷积都可分解成行和列的一维卷积。例如，在第一层，首先用,h,0,(,x,),和,h,1,(,x,),分别与图像,f,1,(,x,，,y,),的每行作卷积并丢弃奇数列,(,以最左列为第,0,列,),。接着这个,N,(,N,/2),阵列的每列再和,h,0,(,x),和,h,1,(,x),相卷积，丢弃奇数行,(,以最上行为第,0,行,),。结果就是该层变换所要求的四个,(,N,/2)(,N,/2),的数组，如图,1.13,所示。,109,.,图1.13,110,.,重构过程与上述过程相似。在每一层，通过在每一列的左边插入一列零来增频采样前一层的四个阵列；接着用,h,0,(,x,),和,h,1,(,x,),来卷积各行，再成对地把这几个,N,/2,N,的阵列加起来；然后通过在每行上面插入一行零，将刚才所得的两个阵列的增频采样为,N,N,；再用,h,0,(,x,),和,h,1,(,x,),与这两个阵列的每列卷积。这两个阵列的和就是这一层重建的结果。重构过程如图,1.14,所示。,111,.,图,1.14,112,.,1.4.3,二进小波变换离散小波变换要求对尺度参数,a,和平移参数,b,进行离散化。为了使小波变换具有可变化的时间和频率分辨率，适应待分析信号的非平稳性，我们很自然地需要改变,a,和,b,的大小，以使小波变换具有“变焦距”的功能。换言之，在实际中采用的是动态的采样网格。最常用的是二进制的动态采样网格：,a,0,2,，,b,0,1,，每个网格点对应的尺度为,2,j,，而平移为,2,j,k,。由此得到的小波,y,j,，,k,(,t,),2,j,/2,y,(2,j,t,k,),(,j,，,k,Z,),(1.60),称为二进小波,(DyadicWavelet),。,113,.,二进小波对信号的分析具有变焦距的作用。假定一开始选择一个放大倍数,2,j,，它对应为观测到信号的某部分内容。如果想进一步观看信号更小的细节，就需要增加放大倍数，即减小,j,值；反之，若想了解信号更粗的内容，则可以减小放大倍数，即增大,j,值。在这个意义上，小波变换被称为数学显微镜。,114,.,定义,1.5,设函数,y,j,，,k,(,t,),L,2,(,R,),，如果存在两个常数,A,、,B,，且,0,A,B,，使得稳定性条件几乎处处成立，即,(1.61),则,y,j,，,k,(,x,),为一个二进小波。若,A,B,，则称为最稳定条件。而函数序列,叫做,f,的二进小波变换，其中,(1.62),115,.,上式相应的逆变换为(1.63)二进小波不同于连续小波的离散小波，它只是对尺度参数进行了离散化，而对时间域上的平移参量保持连续变化，因此二进小波不破坏信号在时间域上的平移不变量，这也正是它同正交小波基相比具有的独特优点。,116,.,1.4.4,二进正交小波变换定义,1.6,设,y,j,，,k,(,t,),L,2,(,R,),