《小波分析基础2》由会员分享，可在线阅读，更多相关《小波分析基础2（82页珍藏版）》请在金锄头文库上搜索。
1、1一、认识小波1、预备知识 从数学的角度讲，小波是构造函数空间正交基的基本单元，是在能量有限空间L2(R) 上满足允许条件的函数，这样认识小波需要L2(R) 空间的基础知识，特别是内积空间中空间分解、函数变换等的基础知识。 从信号处理的角度讲，小波(变换)是强有力的时频分析(处理)工具，是在克服傅立叶变换缺点的基础上发展而来的，所以从信号处理的角度认识小波，需要傅立叶变换、傅立叶级数、滤波器等的基础知识。2(1.1) 一个信号从数学的角度来看，它是一个自变量为时间t的函数f(t)。因为信号是能量有限的，即满足条件(1.1)的所有函数的集合就形成L2(R) 图像是二维信号，同样是能量有限的。实际
2、上任何一幅数字图像都是从真实的场景中经过采样和量化处理后得到的。从数学上看，图像是定义在L2(R2)上的函数。3如图1所示的LENA图像f(x,y)，假设图像的大小是512x512，量化级是256，即xy42、L2(R)空间的正交分解和变换1 对f(t)L2(R)，存在L2(R) 的一组标准正交基gi(t)，t R，i=1,2,使得其中(1.2)(1.3)5对于给定信号f(t)，关键是选择合适的基gi(t)，使得f(t)在这组基下的表现呈现出我们需要的特性，但是如果某一个基不满足要求，可通过变换将函数转换到另一个基下表示，才能得到我们需要的函数表示。常用的变换2有：(1)K-L变换(2)Wal
3、sh变换(3)傅立叶变换(4)小波变换 如如图图所所示示是信号f(t)的傅立叶变换示意图。信号f(t)经傅立叶变换由时域变换到频域，基底不同得到大变换也不同。在信号处理中，有两类非常重要的变换即傅立叶变换和小小小小波波波波变变变变换换换换。目前，可简单地将小波理解为满足以下两个条件的特殊信号：(1)小波必须时振荡的；(2)小波的振幅只能在一个很短的一段区间上非零，即是局部化的。61、Daubechies小波一些著名的小波3：72、Coiflets小波3、Symlets小波84、Morlet小波5、MexicanHat小波6、Meyer小波SKIPSKIP9不是小波的例10RETURNRETUR
4、N113、傅立叶变换与时频分析4我们知道，任何复杂的周期信号f(t)可以用简单的调和振荡函数表示成如下形式：这就是著名的傅立叶级数，都是简单的调和振荡函数，直观讲都是正弦波。是函数f(t)的傅立叶系数，可由以下公式计算：(1.4)12于是，周期函数f(t)就与下面的傅立叶序列产生了一一对应，即从数学上已经证明了，傅立叶级数的前N项和是原函数f(t)在给定能量下的最佳逼近：(1.5)(1.6)(1.7)13对于L2(R)上的非周期函数f(t)，有(1.8)(1.9)称为f(t)的傅立叶变换，反变换公式为(1.10)14有了傅立叶变换，我们可以很容易地将时域信号f(t)转换到频域上，于是信号的频率
5、特性一目了然，并且与傅立叶级数一样，傅立叶变换将一段信号的主要低频能量都集中在频率信号的前面几项，这种能量集中性有利于进一步的处理。在过去200年里，傅立叶分析在科学与工程领域发挥了巨大的作用，但傅立叶分析也有不足，主要表现在以下两点：q傅立叶分析不能刻画时域信号的局部特性；q傅立叶分析对非平稳信号的处理效果不好。下面通过两个例子来说明这两点。15例1、歌声信号歌声是一种声音震荡的波函数，其傅立叶变换就是将这个波函数转化成某种乐谱。但遗憾地是，傅立叶变换无法反映信号在哪一时刻有高音，在哪一时刻有低音，因此结果是所有的音符都挤在了一起，如图所示。16小波变换有效地克服了傅立叶变换的这一缺点，信号
6、变换到小波域后，小波不仅能检测到高音与低音，而且还能将高音与低音发生的位置与原始信号相对应，如图所示。17例2、信号逼近：如图(a)和(b)是原始信号，其余的是逼近信号。1819因此我们需要这样一个数学工具：既能在时域很好地刻画信号的局部性，同时也能在频域反映信号的局部性，这种数学工具就是“小波”。从函数分解的角度，希望能找到另外一个基函数(t)来代替sint。(t)应满足以下三个特性：q任何复杂的信号f(t)，都能由一个母函数(t)经过伸缩和平移产生的基底的线性组合表示；q信号用新的基展开的系数要能反映出信号在时域上的局部化特性；q新的基函数(t)及其伸缩平移要比三角基sint更好地匹配非平
7、稳信号。历史上，Haar第一个找到了这样一个基函数，这就是非常著名但又及其简单的Haar小波。(1.11)20数学上已经证明：小波级数、信号的小波逼近构成L2(R)的一个正交基，通过规范化处理，(1.12)构成L2(R)的一个规范正交基。故任何一个能量有限信号f(t)L2(R) 可以分解为(1.13)(1.14)(1.15)21二、小波变换的定义及特点定义定义11函数(t)L2(R)称为基本小波，如果它满足以下的“允许”条件：(2.1)如果是连续的，易得：(2.2)22(t)又称为母小波，因为其伸缩、平移可构成L2(R)的一个标准正交基：同傅立叶变换一样，连续小波变换可定义为函数与小波基的内积
8、：将a,b离散化，令可得离散小波变换：(2.3)(2.4)(2.5)23 总结：小小小小波波波波即小区域的波，是一种特殊的长度有限、平均值为零的波形。它有两个特点：一是“小”，即在时域具有紧支集或近似紧支集；二是正负交替的“波动性”，也即支流分量为零。傅傅傅傅立立立立叶叶叶叶分分分分析析析析是是是是将将将将信信信信号号号号分分分分解解解解成成成成一一一一系系系系列列列列不不不不同同同同频频频频率率率率的的的的正正正正弦弦弦弦波波波波的的的的叠叠叠叠加加加加，同同同同样样样样小小小小波波波波分分分分析析析析是是是是将将将将信信信信号号号号分分分分解解解解为为为为一一一一系系系系列列列列小小小小波
9、波波波函函函函数数数数的的的的叠叠叠叠加加加加，而而而而这这这这些些些些小小小小波波波波函函函函数数数数都都都都是是是是由由由由一一一一个个个个母母母母小小小小波波波波函函函函数数数数经经经经过过过过平平平平移移移移和尺度伸缩得来的。和尺度伸缩得来的。和尺度伸缩得来的。和尺度伸缩得来的。(2.6)(2.7)24小波分析优于傅立叶分析的地方是，它在时域和频域同时具有良好的局部化性质。而且由于对高频成分采用逐渐精细的时域或频域取样步长，从而可以聚焦到对象的任何细节，所以被称为“数学显微镜”。小波分析广泛应用与信号处理、图像处理、语音识别等领域。25可以这样理解小波变换的含义：打个比喻，我们用镜头观
10、察目标信号f(t)，(t)代表镜头所起的所用。b相当于使镜头相对于目标平行移动，a的所用相当于镜头向目标推进或远离。由此可见，小波变换有以下特点：多尺度/多分辨的特点，可以由粗及细地处理信号；可以看成用基本频率特性为()的带通滤波器在不同尺度a下对信号做滤波。适当地选择小波，使(t)在时域上为有限支撑,()在频域上也比较集中，就可以使WT在时、频域都具有表征信号局部特征的能力。26小波变换的思想来源于伸缩和平移方法。v尺度伸缩对波形的尺度伸缩就是在时间轴上对信号进行压缩和伸展，如图所示。2728v时间平移时间平移就是指小波函数在时间轴上的波形平行移动，如图所示。29小波运算的基本步骤：(1)选
11、择一个小波函数，并将这个小波与要分析的信号起始点对齐；(2)计算在这一时刻要分析的信号与小波函数的逼近程度，即计算小波变换系数C，C越大，就意味着此刻信号与所选择的小波函数波形越相近，如图所示。30(3)将小波函数沿时间轴向右移动一个单位时间，然后重复步骤(1)、(2)求出此时的小波变换系数C，直到覆盖完整个信号长度，如图所示；31(4)将所选择的小波函数尺度伸缩一个单位，然后重复步骤(1)、(2)、(3)，如图所示；(5)对所有的尺度伸缩重复步骤(1)、(2)、(3)、(4)。32v尺度与频率的关系尺度与频率的关系如下：小尺度a压缩的小波快速变换的细节高频部分大尺度a拉伸的小波缓慢变换的粗部
12、低频部分33三、多分辨分析由母小波按如下方式的伸缩平移可构成L2(R)空间的标准正交基如何构造母小波呢？1989年，Mallat和Meyer提出了按多分辨分析的思想来构造母小波，其基本思想是：q现构造一个具有特定性质的层层嵌套的闭子空间序列VjjZ，这个闭子空间序列充满了整个L2(R)空间。q在V0子空间找一个函数g(t)，其平移g(t-k)k Z构成V0子空间的Riesz基。q对函数g(t)进行正交化，得到函数称为正交尺度函数(t)。q由(t)计算出小波函数(t)。1 1、多分辨分析、多分辨分析、多分辨分析、多分辨分析(MRA)(MRA)的概念的概念的概念的概念55(3.1)34Riesz基
13、定义令H是Hilbert空间，H中的一个序列gjjZ是Riesz基，如果它满足以下的条件：A和B分别称为Riesz基的上下界，Riesz基又称为稳定基。(3.2)(3.3)35定义定义1 空间L2(R )中的多分辨分析是指L2(R )中的满足如下条件的一个子空间序列36多分辨空间的关系可用下图来形象地说明。37如果g(t-k)kZ是V0的Riesz基，可通过正交化得到V0空间的函数(t)V0，使得(t-k)kZ构成V0空间的规范正交基。由伸缩性和平移不变性可知，j,k(t)j,kZ构成Vj空间的一个规范正交基。于是(3.4)(3.5)38注意：(t)并不是L2(R )空间的小波函数，而是与其紧
14、密相关的尺度函数，j,k(t)j,kZ称为尺度基，多分辨空间序列VjjZ称为尺度空间，在MRA意义下，可由尺度基导出小波基。由MRA的单调性可以看出：Vj是Vj+1的严格子空间，设Wj是Vj关于Vj+1的正交补(子空间)，即(3.6)39对于一幅图像，量化级数决定了图像的分辨率，量化级数越高，图像就越清晰，即图像的分辨率高。对于任意一幅图像，都可以用不同的量化空间来表示，细节比较丰富的部分用高分辨率来表示，细节比较单一的部分可用低分辨率来表示。我们可以将不同的量化级数构成的空间看成不同的多分辨空间Vj，显然这些量化空间是相互嵌套的，(3.7)40从图像处理的角度，多分辨空间的分解可以理解为图像
15、的分解，假设有一幅256级量化的图像，不妨将它看成量化空间Vj中的图像，则可理解为Vj空间中的图像有一部分保留在Vj-1空间中，还有一部分放在Wj-1空间，如图所示如图所示如图所示如图所示。与尺度函数的产生一样，若存在(t)W0，使得(t-k)kZ构成空间W0的一个规范正交基，则构成L2(R)空间的一个规范正交基。 称为小波基，(t)称为母小波。(3.8)SKIPSKIP41VjWj-1Vj-1RETURNRETURN42MRA非常抽象，但是它给出了构造小波的一般框架。在实践中很难通过小波空间直接构造小波，但通过MRA可推导出一个非常重要的关系：双尺度方程，通过求解该方程，使我们有可能求出尺度
16、函数和小波函数。2 2、双尺度方程、双尺度方程、双尺度方程、双尺度方程由前面的分析，我们知道：(3.9)(3.10)43方程(3.9)和(3.10)称为双尺度方程。由(t)的正交性可得：对双尺度方程两边取傅立叶变换，可得频域上的的双尺度方程：(3.12)(3.11)(3.14)(3.13)44(3.16)(3.15)从信号处理的角度，h是与(t)对应的低通滤波器，g是与(t)对应的高同滤波器，h，g既可以表示为时域上的离散序列形式hk，gkkZ，也可以表示为频域上的2周期函数h ()，g()。两者本质上是一样的。45若kN时，hk=0，这样的滤波器称为有限脉冲响应滤波器(FIR)，FIR滤波器
17、具有好的局部化特性。此时，(t)只在有限区间0，N上取值，所以(t)是紧支的，其支集supp=0，N，(3.9)式变为：(3.17)此时(t)也是紧支的。所以只要滤波器的长度是有限的，我们称对应的小波(t)是紧支小波。46由(3.13)式得：(3.18)(3.19)47结论：结论：结论：结论：只要找到满足双尺度方程(3.9)的序列hkkZ，通过公式(3.15)就可以计算出2周期函数h ()，再由公式(3.19)就可以计算出，经过傅立叶反变换，最终可得尺度函数(t)，有了尺度函数就可以计算出小波函数(t)。通过解双尺度方程(3.9)，我们希望得到满足MRA的尺度函数(t)，并最终构造出小波函数(
18、t)，但有两个问题必须解决：问题问题问题问题1 1：双尺度方程(3.9)是否有解？解的唯一性如何？问题问题问题问题2 2：双尺度方程(3.9)的解是否满足MRA？关于问题1，I.Daubechies和Lagarias7在1991年给出了证明。48解决问题2却是一件非常困难的事情。这里牵涉到尺度函数(t)与滤波器系数hkkZ之间的关系问题：q如果有一个L2(R)空间的尺度函数(t)，一定能构造出双尺度方程(3.9)，从而找到一组满足(3.9)的滤波器hkkZ；q反过来，如果有一组滤波器hkkZ满足某个双尺度方程，由此求解得到的函数却不一定是满足MRA的尺度函数，这样无法保证双尺度方程解的平移构成
19、L2(R)Riesz基若(t)是正交的，则相应的滤波器h有什么性质呢？定理定理定理定理1 133若(t)是正交的，则相应的滤波器hk必须满足条件：(3.20)(3.21)但是，如果hk仅仅满足(3.20)和(3.21)，并不能保证由双尺度方程构造出的函数(t)是正交尺度函数。(3.20)和(3.21)称为构造正交小波的必要条件必要条件必要条件必要条件。49仅有必要条件是不够的，即hkkZ除了满足条件(3.20)和(3.21)外，还应满足其他条件。S.Mallat4，W.Lawton6等都在这方面作出了重大的贡献，并给出了一些有意义的结论。下面给出W.Lawton的充分条件。定理x2设h()是F
20、IR滤波器，若满足若矩阵A的特征值1是非退化的，则(t-k)kZ是标准正交的。50算法：算法：算法：算法：构造紧支小波基步骤步骤步骤步骤1 1寻找满足双尺度方程(3.9)和(3.10)的滤波器hk,gkk0,1,N步骤步骤步骤步骤2 2利用公式(3.15)计算2周期函数h()；步骤步骤步骤步骤3 3验证h()是否满足条件通过傅立叶反变换求出(t)步骤步骤步骤步骤5 5验证矩阵A的特征值1是否非退化；步骤步骤步骤步骤6 6(t-k)kZ是正交的尺度函数，对应的紧支小波由公式(3.10)计算。步骤步骤步骤步骤4 4 4 4 计算513 3、小波与共轭镜像滤波器、小波与共轭镜像滤波器、小波与共轭镜像
21、滤波器、小波与共轭镜像滤波器44我们知道尺度函数和小波函数(t),(t)tR是在时域刻画信号的性质，对应的滤波器h(),g()R从频域上刻画信号的性质。实际上，(t),(t)tR大量的性质都可以由对应的h(),g()R从频域上反映出来，甚至离散小波变换都可以借助滤波器来实现，因此小波与滤波器具有紧密的关系。3.1 3.1 正交尺度函数产生共轭镜像滤波器正交尺度函数产生共轭镜像滤波器正交尺度函数产生共轭镜像滤波器正交尺度函数产生共轭镜像滤波器定定定定义义义义若尺度函数(t)是正交的，则它所对应的滤波器h()称为共轭镜像滤波器。共轭镜像滤波器。共轭镜像滤波器。共轭镜像滤波器。h()满足以下条件：5
22、2滤波器hkkZ称为低通滤波器。所谓低通是指：当信号f(t)被hkkZ作用后，其低频成分能被保留下来，而高频成分(=)却被滤掉了。对应的小波滤波器g()也是共轭镜像滤波器。共轭镜像滤波器。共轭镜像滤波器。共轭镜像滤波器。也满足条件(3.22)另外，由于(t-k)kZ与(t-k)kZ分别是V0空间和W0空间的规范正交基，而V0W0，则(3.23)公式(3.23)反映了低通滤波器h()和高通滤波器g()之间的关系。RETURNRETURN53S.Mallat4同时给出了这样的结论：若高通滤波器g()满足公式(3.22)(3.22)和和和和(3.23)(3.23)，则由公式产生的小波基(t-k)kZ
23、构成W0空间的规范正交基。因此当尺度函数(t)已经确定时，只要能找到一个满足公式(3.22)和(3.23)的g()，就一定能找到对应的小波(t)，但是这样的解并不是唯一的。例如可取(3.24)可以验证g()满足(3.22)和(3.23)，对应的共轭镜像滤波器为：(3.25)54因此当找到低通共轭镜像滤波器hkkZ后，利用公式(3.25)马上可得高通共轭镜像滤波器gkkZ。总总总总 结结结结 ： 在 一 个 MRA下 的 正 交 尺 度 函 数 和 小 波 函 数(t),(t)tR，产生一组共轭镜像滤波器h,g，满足：(3.26)公式(3.26)还有几个等价形式，下面以定理的形式给出。55定理定
24、理定理定理设h,g是由正交尺度函数和小波函数产生的共轭镜像滤波器，则以下几个条件等价：q在频域上(3.26)式成立；q在时域上以下公式成立：(3.27)q定义调制矩阵：(3.28)56则(3.29)3.2 3.2 利用共轭镜像滤波器实现快速正交小波变换利用共轭镜像滤波器实现快速正交小波变换利用共轭镜像滤波器实现快速正交小波变换利用共轭镜像滤波器实现快速正交小波变换44L2(R)空间的一个MRA产生了两个子空间：尺度空间VjjZ和小波空间WjjZ。j,kj,kZ和j,kj,kZ分别是两个空间的规范正交基，信号f(t)L2(R)在两个空间上都可以做正交投影：(3.30)57信号在小波空间的展开为(
25、3.31)但实践中不可能进行无穷次逼近，不妨设f(t)VJ，则因为所以表示从尺度2-J到2-j进行了(J-j)次小波分解(jJ)58实际计算时，可以一次一次地进行小波分解，然后递推实现(J-j)次小波分解，不妨记一次小波分解的尺度系数和小波系数为(3.32)59而因为代入(3.32)式得故60从而我们得到如下的递推公式：(3.33)61现在来求dj,k的递推公式，(3.34)而62因为代入(3.34)式得故63从而我们得到如下的递推公式：(3.35)64通过公式(3.33)和(3.35)，可以很快计算出尺度系数和小波系数cj,k,dj,k，这就是著名的Mallat算法：因此，只要确定VJ空间的
26、初始序列cJ,kkZ，就可以算出任意空间Vj(jJ)的所有尺度系数和小波系数。公式(3.33)和(3.35)称为离散小波变换的分解公式。65又由于Vj+1=VjWj，VjWj，因此Vj上的标准正交基与Wj上的标准正交基是相互正交的。它们共同构成Vj+1上的标准正交基，则Vj+1上的函数j+1,nj,nZ可以由这两个基共同表示：有前面的计算可知：故66从而这就是Mallat重构算法：67小波的应用1,4,8,9小波的应用主要是信号的处理，其中最典型的应用是小波图象压缩。另外，小波在诸如信号去噪、特征提取等多方面均有成功的应用。下面以图象去噪为例说明小波应用策略。小波的各种应用均可分为以下三步：1
27、）对原始信号作小波变换，将信号由空域变换到频域；2）对小波系数做相应处理；3）对处理后的小波系数做小波逆变换，还原原信号。68小小波波信信号号去去噪噪一一69小小波波信信号号去去噪噪二二70图图像像融融合合71小波图像去噪小波图像去噪8因为噪声信号多包含在具有较高频率的细节中，所以小波去噪首先对图像信号进行小波分解，可利用门限阈值对所分解的小波系数进行处理，然后对图像信号进行小波重构，抑制图像信号中的无用部分，恢复图像信号中的有用部分。具体步骤为：图像信号的小波分解：选择合适的小波及恰当的分解层次N，对目标图像进行N层的小波分解；对分解后的高频系数进行阈值量化：对于分解的每一层，选择恰当的阈值
28、，对该层高频系数进行阈值量化处理；重构图像：根据小波分解后的第N层近似的低频系数和经过阈值量化处理后的细节高频系数，重构图像。 7273参考文献1唐远炎,王玲.小波分析与文本文字识别,科学出版社,20042李弼程,彭天强,彭波.智能图像处理技术,电子工业出版社,20043 I. Daubechies, Ten Lectures on Wavelets. Philadelphia: SIAM,1992.4S.Mallat.Awavelettourofsignalprocessing.AcademicPress,USA,19985S.Mallat,“Atheoryformultiresolutio
29、nsignaldecomposition:Thewavelet representation,” IEEE Trans. Pattern Anal. Machine Intell., vol.11,pp.674693,1989.6 W. Lawton. Tight frames of compactly supported wavelets, J.Math.Phys.,31:18981901,1990.747I.Daubechies,J.C.Lagarias.Two-scaledifferenceoperationsI: existence and glogal regularity of s
30、olutions, SIAM J. Math.Anal.,22:13881410,1991.8 D. Donoho, “De-noisingby soft-thresholding,”IEEE Trans. Inform.Theory, vol.41,pp.613627,1995.9B.JAWERTH,etc.“anoverviewofwaveletbasedmultiresolution analyses,” SIAM REVIEW,vol. 36,No.33,pp.377-412,September,1994.7576简单实例例 Fourier变换close all, clear all,
31、t=0:0.001:1.3;x=sin(2*pi*50*t)+sin(2*pi*300*t);f1=x+3.3*randn(1,length(t);subplot(121);plot(f1);xlabel(时间);ylabel(幅值);title(原始信号);y=fft(f1,1024);p=y.*conj(y)/1024; f2=1000*(0:511)/1024;subplot(122); plot(f2,p(1:512);xlabel(频率); ylabel(功率谱密度); title(信号功率)77例 连续小波变换close all, clear all,t=0:0.01:1;x=3*
32、sin(100*pi*t)+2*sin(64*pi*t)+5*cos(72*pi*t);f=x+randn(1,length(t);coefs=cwt(f,1:0.2:3,db3,plot);xlabel(时间); ylabel(尺度);title(不同的尺度小波变换系数值)78例 计算尺度方程系数close all, clear all,w=db5;Lo_D,Hi_D,Lo_R,Hi_R=wfilters(w);subplot(221);stem(Lo_D);title(分解低通滤波器);grid;subplot(222);stem(Hi_D);title(分解高通滤波器);grid;sub
33、plot(223);stem(Lo_R);title(重构低通滤波器);grid;subplot(224);stem(Hi_R);title(重构高通滤波器);grid;79例 一维离散小波变换close all, clear all,x=sin(0.3*1:300); lx=length(x); w=db9;subplot(5,1,1); plot(x); title(原始信号);Lo_D,Hi_D,Lo_R,Hi_R=wfilters(w);ca1,cd1=dwt(x,w);x1=idwt(ca1,cd1,w,lx);subplot(5,2,5);plot(x1);title(idwt重构
34、造信号);errx1=x-x1;errx1max=max(abs(x-x1)subplot(5,2,6);plot(errx1);title(idwt重构造误差);axis(0,300,-2e-11,2e-11);80例 含噪图像的处理close all, clear all,load woman2;subplot(221); image(X); colormap(map);title(原始图像);axis squareinit=2055615777; randn(seed,init);x=X+28*randn(size(X);subplot(222);image(x);colormap(map);81title(含噪图像);axis squarec,l=wavedec2(x,2,sym4);n=1,2;p=19,19;nc=wthcoef2(h,c,l,n,p,s);nc=wthcoef2(v,c,l,n,p,s);nc=wthcoef2(d,c,l,n,p,s);xx=waverec2(nc,l,sym4);subplot(223);image(xx);colormap(map);title(消噪的图像);axis square82
