小波分析完整 第三章小波变换ppt课件

.,第三章小波变换(WaveletTransform，WT),.,Fourier分析对于平稳信号的分析是非常有效的,但是由于其不具有分析局部时域信号的局部频谱特性,没有时频局部化功能;短时Fourier变换虽然具有时频分析能力,但是时频窗是固定的,不能自动适应频率变化的需要.,.,3-1自适应窗函数的设计,Fourier变化的表达式为:只能求取整个时域上的频率特性，不能反映局部频率特性。短时Fourier变化的表达式为:,.,令,，则：,可以反映局部频率特性，但是窗函数一经设定，没有自适应能力，不能满足低频部分需要时窗宽、频窗窄，高频部分需要时窗窄、频窗宽的要求。为此，定义窗函数的一般形式为：它是经过平移和放缩的结果。,.,当a1时，被压窄且振幅被拉高，含有表现高频分量的特征。,.,3-2小波、小波变换的定义和条件,1、小波变换的定义把对给定平方可积信号f(t)的积分变换：称为小波变换。满足一定条件的函数称为允许小波函数。又称为基本小波，或母小波。是母小波经时移和伸缩所产生的一族函数，我们称之为小波基函数，或简称小波基。b是时移因子，a是尺度因子。式中a,b和t均是连续变量，因此该式又称为连续小波变换（CWT),.,在上式中，b的作用是确定对x(t)分析的时间位置。尺度因子a的作用是把基本小波(t)作伸缩。由(t)变成(at)，当a1时，若a越大，则(at)的时域支撑范围（即时域宽度）较之(t)变得越小；反之，当a1时，a越小，则(at)的宽度越宽。这样，a和b联合越来确定了对x(t)分析的中心位置及分析的时间宽度。,.,.,小波函数的范数不变性：此式表明：经过平移与伸缩以后，其模量没有改变。在不同的尺度a时，ab(t)终能和母函数(t)有着相同的能量。小波函数的频域特性：此式表明，经过平移和伸缩以后得到的新函数的频域特性随参数a的变化而变化。,.,2、小波变化的回复公式推导任何一种变换应该是可逆的。为推导小波变换的回复公式，先得推出与Fourier变换中类似的乘积公式。在Fourier变换中，有公式：对于小波变换而言:,.,同理：,用到公式,.,式中：,，,.,补充：（1）函数特点：质量为1的质点均匀分布在x轴的区间，平均密度分布函数为：若将该质点放置在座标原点，则密度分布函数为,.,由此，引出函数概念：,.,（2）阶跃函数,.,就有：，因此：,.,允许小波的条件：（1），即小波具有快速衰减性（平方可积，能量有限）；（2）由的连续性和可推知：。从而有：,这一结论指出，(t)的取值必然是有正有负，也即它是振荡的。,.,以上给我们勾画出了作为小波的函数所应具有的大致特征，即(t)是一带通函数，它的时域波形应是振荡的。此外，从时频定位的角度，我们总希望(t)是有限支撑的，因此它应是快速衰减的。这样，时域有限长且是振荡的这一类函数即是被称作小波（wavelet）的原因。,.,3.3小波变换的自适应时频窗分析,前面已经提到，WFT不具有自适应时频窗的能力，而小波变换具有这种能力，下面具体分析其时频局部化效果。与WFT相比，可将看成窗函数，令为时窗中心，为时窗半径，为频窗中心，为频窗半径，有：,.,如式：将窗函数的以上各量分别记为:，,。,.,可以证明：,，,.,.,.,由以上关系不难看出：平移因子b使时窗中心发生同量的平移，不影响时窗的宽度；而放缩因子a既可使时窗中心发生变化，也使时窗半径压缩为1/a；同样缩放因子将频窗中心和频窗半径拉宽a倍。,.,总结如下：,.,.,连续小波变换的计算性质,1时移性质,若f(t)的CWT是Wf(a,b)，那么f(t)的CWT是Wf(a,ba)。,证明：令，g(t)=f(t)，则,.,连续小波变换的计算性质,2尺度转换性质,如果f(t)的CWT是Wf(a,b)，令g(t)=f(t)，则,证明：,该性质指出，当信号的时间轴按作伸缩时，其小波变换在a轴上同时要作相反比例的伸缩，但小波变换的波形不变。这是小波变换优点的又一体现。,.,连续小波变换的计算性质,3微分性质,.,4小波变换的内积（乘积）定理(Parseval等式）,如果令,该式更清楚地说明，小波变换的幅平方在尺度位移平面上的加权积分等于信号在时域的总能量，因此，小波变换的幅平方可看作是信号能量时频分布的一种表示形式。,.,5、叠加性质设的小波变换为，的小波变换为则,.,常用小波,由前面的讨论可知，作为一个小波的函数(t)，它一定要满足容许条件，在时域一定要是有限支撑的，同时，也希望在频域也是有限支撑的，当然，若时域越窄，其频域必然是越宽，反之亦然。在时域和频域的有限支撑方面我们往往只能取一个折中。此外，我们希望由母小波(x)形成的ab(t),是两两正交的，或是双正交的；进一步，我们希望(x)有高阶的消失矩，希望与(x)相关的滤波器具有线性相位，等等。我们可以根据上述要求对现已提出的大量的小波函数作一粗略地分类。在下面的分类中，第一类是所谓的“经典小波”，在MATLAB中把它们称作“原始（Crude）小波”。这是一批在小波发展历史上比较有名的小波；第二类是Daubecheis构造的正交小波，第三类是由Cohen，Daubechies构造的双正交小波。,.,经典类小波,1.Haar小波Haar小波来自于数学家Haar于1910年提出的Haar正交函数集，其定义是：,(t)的傅里叶变换是：,.,Haar小波有很多好的优点，如：（1）Haar小波在时域是紧支撑的，即其非零区间为（0，1）；（2）若取a=2j,jZ+,bZ，那么Haar小波不但在其整数位移处是正交的，即=0，而且在j取不同值时也是两两正交的，即=0（3）Haar波是对称的。我们知道，离散的单位抽样响应若具有对称性，则该系统具有线性相位，这对于去除相位失真是非常有利的。Haar小波是目前唯一一个既具有对称性又是有限支撑的正交小波；（4）Haar小波仅取1和1，因此计算简单。,.,2.Morlet小波,Morlet小波定义为,其傅里叶变换,它是一个具有高斯包络的单频率复正弦函数。考虑到待分析的信号一般是实信号，所以在MATLAB中将其式改造为：,并取0=5。该小波不是紧支撑的，理论上讲t可取+。但是当0=5，或再取更大的值时，(t)和()在时域和频域都具有很好的集中。Morlet小波不是正交的，也不是双正交的，可用于连续小波变换。但该小波是对称的，是应用较为广泛的一种小波。,.,.,3.Mexicanhat小波,该小波的中文名字为“墨西哥草帽”小波，又称Marr小波。它定义为,该小波是由一高斯函数的二阶导数所得到的，它沿着中心轴旋转一周所得到的三维图形犹如一顶草帽，故由此而得名。其波形和其频谱如图,helpmexihat,.,该小波不是紧支撑的，不是正交的，也不是双正交的，但它是对称的，可用于连续小波变换。由于该小波在=0处有二阶零点，因此它满足容许条件，且该小波比较接近人眼视觉的空间响应特征，因此它在1983年即被用于计算机视觉中的图像边缘检测。,.,4Gaussian小波,高斯小波是由一基本高斯函数分别求导而得到的，定义为：,该小波不是正交的，也不是双正交的，也不是紧支撑的。当k取偶数时(t)正对称，当k取奇数时，(t)反对称。上图给出了k=4时的(t)的时域波形及对应的频谱。,.,正交小波,目前提出的正交小波大致可分为四种，即Daubechies小波，对称小波，Coiflets小波和Meyer小波。这些正交小波和前面所讨论的“经典小波”不同，它们一般不能由一个简洁的表达式给出(t)，而是通过一个叫做“尺度函数（Scallingfunction）”的(t)的加权组合来产生的。尺度函数是小波变换的又一个重要概念。由后面的讨论可知，小波函数(t)，尺度函数(t)同时和一个低通滤波器H0(z)及高通滤波器H1(z)相关连，H0(z)和H1(z)可构成一个两通道的分析滤波器组。这些内容构成了小波变换的多分辨率分析的理论基础。因此，在讨论正交小波时，同时涉及到尺度函数(t)，分析滤波器组H0(z)、H1(z)，及综合滤波器组G0(z)，G1(z)。,.,1Daubechies小波,Daubechies小波简称db小波。它是由法国女学者IngridDauechies于90年代初提出并构造的。Daubechies对小波变换的理论做出了突出的贡献，特别是在尺度a取2的整数次幂时的小波理论及正交小波的构造方面进行了深入的研究，其代表作TenLecturesonWavelet（小波十讲）深受同行们的欢迎。,dbN中的N表示db小波的阶次，N=210.当N=1时，db1即是Haar小波。因此，前述的Haar小波应归于“正交小波”类。Daubechies计算出了N=210时的(t),h0,h1,g0及g1。在MATLAB5.3中，N的阶次还可以扩展。db小波是正交小波，当然也是双正交小波，并是紧支撑的。(t)的支撑范围在t=0(2N1)，(t)的支撑范围在(1N)N。小波(t)具有N阶消失矩，()在=0处具有N阶零点。但db小波是非对称的，其相应的滤波器组属共轭正交镜像滤波器组（CQMFB）。,.,2.对称小波,对称小波简记为symN，N=2,3,L,8，它是db小波的改进，也是由Daubechies提出并构造的。它除了有db小波的特点外，主要是(t)是接近对称的，因此，所用的滤波器可接近于线性相位。下图是N=4时的对称小波。,.,3.Coiflets小波,该小波简记为coifN，N=1,2,5.在db小波中，Daubechies小波仅考虑了使小波函数(t)具有消失矩（N阶），而没考虑尺度函数(t)。R.Coifman于1989年向Daubechies提出建议，希望能构造出使(t)也具有高阶消失矩的正交紧支撑小波。Daubechies接受了这一建议，构造出了这一类小波，并以Coifman的名字命名。coifN是紧支撑正交、双正交小波，支撑范围为6N1，也是接近对称的。(t)的消失矩是2N，(t)的消失矩是2N1。下图是N=4时的coif4小波。,.,4Meyer小波,Meyer小波简记为meyr，它是由Meyer于1986年提出的。该小波无时域表达式，它是由一对共轭正交镜像滤波器组的频谱来定义的。Meyer小波是正交、双正交的，但不是有限支撑的，但其有效的支撑范围在8，8之间。该小波是对称的，且有着非常好的规则性。下图给出了Meyer小波的尺度函数(t)和小波函数(t)。,.,双正交小波,两通道正交镜像滤波器组具有仿酋性质。满足这一条件的分析滤波器H0(z)和H1(z)是功率对称的，且h0(n)和h1(n)之间有着正交性，再是h0(n)，h1(n)，g0(n)，g1(n)有着同样的长度，都不是线性相位的。为了取得线性相位的滤波器组，我们需放弃H0(z)的功率互补性质。这也就放弃了h0(n)和h1(n)之间的正交性，代之的是双正交关系。由于离散小波变换最后是由两通道滤波器组来实现。因此，正交小波条件下的(t)，(t)和h0,h1,g0与g1都不具有线性相位（Haar小波除外）。为此，Daubechies和Cohen提出并构造了双正交小波，其目的是在放宽小波正交性的条件下得到线性相位的小波及相应的滤波器组。双正交滤波器组简称biorNr,Nd，其中Nr是低通重建滤波器的阶次，Nd是低通分解滤波器的阶次。在MATLAB中，Nr和Nd的可能组合是：Nr=1,Nd=1,3,5Nr=2,Nd=2,4,6,8Nr=3,Nd=1,3,5,7,9Nr=4,Nd=4Nr=5,Nd=5Nr=6,Nd=8,.,这一类小波自然不是正交的，但它们是双正交的，是紧支撑的，更主要的是它们是对称的，因此具有线性相位。分解小波(t)的消失矩为Nr1。下图给出的bior3.7的分解小波、尺度函数及重建小波和尺度函数。,.,3.4离散小波变换及其频带特性,将a、b取为离散型整数形式，即：，可将表示为：相应的小波变换表示为离散小波变换：,.,因为离散小波是由经整数倍缩放和整数k平移所生成，它同样满足：若满足允许小波的条件，则也是允许小波；若具有时频局部化能力，也具有时频局部