第四章_多分辨率分析与正交小波变换.ppt

多分辨率分析与正交小波变换,概述,多分辨率是小波分析中的最重要的概念之一，它从函数空间的高度研究函数的多分辨率表示将一个函数表示为一个低频成分与不同分辨率下的高频成分。更重要的是，多分辨率能够提供一种构造小波的统一框架，并且能够提供函数分解与重构的快速算法。,本章主要内容,多分辨率分析尺度函数和小波函数二尺度方程及多分辨率滤波器组二进正交小波变换的Mallat算法常见小波函数,1.多分辨率分析,定义：多分辨率分析（MultiresolutionAnalysis,MRA）是用小波函数的二进伸缩和平移表示函数这一思想的更加抽象复杂的表现形式，它重点处理整个函数集，而非侧重处理作为个体的函数。基本思想：将L2（R）用它的子空间Vj，Wj表示，其中Vj，Wj分别称为尺度空间和小波空间。,补充：直和,设E是线性空间，L1，L2，Ln是E的子空间，如果任一元素xE可以惟一表示成x=x1+x2+xn,其中xkLk(k=1,2,n),则称E是L1，L2，Ln的直和，记为：,尺度函数，j0，-1，-2，-3；k=0，1,2，,（这里暂对j和k的范围做了限制）形成了伸缩平移系统，其中j不同，张成了不同的子空间，如图：,k=0，1，7，张成了子空间；,k=0，3，张成了子空间；,k=0，1，张成了子空间；,k=0,张成了子空间。由图可知：,比喻,类似于人的视觉系统。例如：人在观察某一目标时，不妨设他所处的分辨率为j（或2j），观察目标所获得的信息是Vj，当他走近目标，即分辨率增加到j-1（或2j-1），他观察目标所获得的信息为Vj-1，应该比分辨率j下获得的信息更加丰富，即，分辨率越高，距离越近；反之，则相反。,在分辨率分析中，Vj称为逼近空间，我们把平方可积的函数f(t)L2(R)看成是某一逐级逼近的极限情况。每次逼近都是用一低通平滑函数（t）对f(t)做平滑的结果，在逐级平滑时平滑函数（t）也做逐级逼近，这就是多分辨率，即用不同分辨率来逐级逼近待分析函数f(t)。,见word,性质,尺度空间Vj具有以下递归嵌套关系：将Vj，Vj1相关联的关键性质是：如：f(t)Vj，则f(t/2)Vj-1,f(2t)Vj+1。位移不变性：函数的时移不改变其所属空间，即如果f(t)Vj，则f(t-k)Vj。,空间的剖分是完整的，即当j-，VjL2（R），包含整个平方可积的实变函数空间。当j+，Vj0，即空间最终剖分到空集为止。,V0中的任意函数f(t)均可表示为的线性组合，我们设P0f(t)代表f(t)在V0上的投影，则有：是线性组合的权重，其求法如下：我们称P0f(t)为f(t)在V0处的平滑逼近，也就是f(t)在j=0下的概貌，称为f(t)在分辨率j=0下的离散逼近。,我们刚才推导出但是，毕竟不等于，也即比对x(t)近似的好，但二者之间肯定有误差。这一误差是由(tk)和(tk)的宽度不同而产生的，因此，这一差别应是一些“细节”信号，我们记之为。这样，有P0 x(t)=P1x(t)+D1x(t)该式的含义是：x(t)在高分辨率基函数所形成的空间中的近似等于它在低分辨率空间中的近似再加上某些细节。现在我们来寻找D1x(t)的表示方法。,我们设D1f(t)代表f(t)在W1上的投影，有是线性组合的权重，其求法：,进行类推，可得：Pjf(t)是f(t)在Vj中得投影，是f(t)在分辨率j下得平滑逼近，称为f(t)在分辨率j下得离散逼近。Djf(t)是f(t)在Wj中得投影，反映了Pjf(t)和Pj-1f(t)之间的细节差异。就是。,我们把空间做逐级二分解产生一组逐级包含的子空间：,各空间内的结构做进一步分析：,（1）设V0中有低通平滑函数（t），它的整数移位集合是V0中的正交归一基。我们称为尺度函数，所以有：,（2）根据二尺度伸缩性，如果（t）V0，则（t/2）V1，而且，如果是V0中的正交归一基，则,（3）如果在子空间W0中能找到一个带通函数，其整数位移的集合构成W0中的正交归一基，我们根据二尺度的伸缩性，可得W1中的任意函数f(t)均可以表示为的线性组合。,多分辨率概念,1.平移不变性。2.单调性。3.伸缩性。4.逼近性。5.Riesz基存在性。,性质1说明，函数的时移不改变所属的空间，等效为,2.尺度函数和小波函数,2.1尺度函数及其空间定义：函数为尺度函数，若其经过整数平移k和尺度j上的伸缩，得到一个尺度和位移均可变化的函数集合：称每一个尺度j上的平移系列jk（t）所组成的空间Vj为尺度为j的尺度空间。,对于任意函数所以，尺度函数在不同尺度下其平移系列组成了一系列的尺度空间。,2.2小波函数及其小波空间,L2（R）的正交基就是把直和的子空间的正交基合并起来。所以L2（R）的标准正交基为：比较二进小波的函数形式。,尺度函数和小波函数性质：,（1）尺度函数（2）小波函数,（3）同一尺度下，因为WjVj，所以小波函数和尺度函数之间是正交的，即：,3.二尺度方程及多分辨率滤波器组,即,类推到Wj和Vj-1之间，得：上面二式就是二尺度差分方程，其中，h0k和h1k是线性组合的权重。由于是正交归一基，它们值为：,二尺度关系存在于任意相邻尺度j和j-1之间设H0()为h0k的傅立叶变换，H1()为h1k的傅立叶变换，它们都是以2为周期的周期函数。,3.2滤波器系数h0k和h1k的性质,（1）h0k和h1k的总和分别为（2）频域初值,（3）递推关系,（4）滤波器H0（），H1（）特性：前两个式子是设计H0（），H1（）的主要依据，第三个式子给出了H0（）与H1（）之间的内在联系。它在时域中的表达式为：,4.二进正交小波变换的Mallat算法,根据多分辨率理论，Mallat提出了小波分解与重构的快速算法，称为Mallat算法，其在小波分析中的作用相当于FFT在傅立叶分析中的作用。它标志着小波分析走上了宽阔的应用领域。,4.1Mallat算法的信号分解过程,在多分辨率分析中，我们得出一个重要结论：,可推得：我们称上式为离散平滑逼近，下式是离散细节信号。,分解算法图例,4.2Mallat算法的信号重建过程,由前面所以,是由它们重建得到的第j-1级离散平滑信号。G0(k)、G1(k)为：如果从设计滤波器的角度考虑，设输入信号为x(k)，重建输出信号为y(n)，我们将x(k)进行二插值，得x(k/2)，k为偶数，所以：,重构算法图例（是分解过程的相反过程，要经过插值和滤波）,Mallat算法得分解与重构比较：,（1）在分解算法中信号是先滤波后抽取，而在重建算法中是先插值后滤波。（2）