一维小波变换.doc

一维小波变换 　 现在可以正式定义若干密切相关的小波变换：一般小波序列展开、离散小波变换和连续小渡变换。它们在傅里叶域的对应部分分别是傅里叶序列展开、离散傅里叶变换和连续傅里叶变换。在7.4节，将定义一种计算效率很高的称做快速小波变换的离散小波变换。 　 一.小波序列展开 首先根据小波φ(x)和尺度函数φ(x)为函数f(x)∈L2(R)定义小波序列展开。据式(7.2.27)，可写出： (7.3.1) 其中j0是任意开始尺度，cj0和dj(k)分别是式(7.2.12)和式(7.2.21)中αk的改写。cj0(k)通常称为近似值或尺度系数；dj(k)称为细节或小波系数。这是因为式(7.3.1)的第一个和式用足度函数提供了f(x)在尺度j0的近似[除非f(x)∈Vj0，此时为其精确值]。对于第二个和式中每一个较高尺度的j≥j0，更细分辨率的函数(一个小波和)被添加到近似中以获得细节的增加。如果展开函数形成了一个正交基或紧框架(通常情况下是这样)，基于式(7.2.5)和式(7.2.9)的展开系数计算如下： (7.3.2) 和 (7.3.3) 如果展开函数是双正交基的一部分，上式中的φ和φ项要分别由它们的对偶函数和代替。 　 例7.7 y=x2的哈尔小波序列展开 考虑图7.13(a)中显示的简单函数： 使用哈尔小波——见式(7.2.14)和式(7.2.30)——和开始尺度J0=0，式(7.3.2)和式(7.3.3)可以被用来计算下述展开系数： 将这些值代入式(7.3.1)，可以得到小波序列展开： 上述展开中的第一项用c0(0)生成待展开函数的V0子空间近似值。该近似值如图7.13(b)所示，是原始函数的平均值。第二项使用d0(0)通过从W0子空间添加一级细节来修饰该近似值。添加的细节及V1的结果近似值分别如图7.13(c)和(d)所示。其他级 别的细节由子空间W1的系数d1(0)和d1(1)给出。 该附加细节如图7.13(e)所示，V2的结果近似值如图7.13(f)所示。注意，展开函数现在已经接近原始函数了。越高的尺度(或细节的级别越高)被叠加，近似值越变得接近函数的精确表示，它的极限是j→∞。 ?二.离散小波变换 与傅里叶序列展开相似，前节的小波序列展开将一个连续变量函数映射成一系数序列。如果待展开函数是一个数字序列，如连续函数f(x)的抽样值，得到的系数就称为f(x)的离散小波变换(DWT)。在这种情况下，式(7.3.1)到式(7.3.3)中定义的序列展开变成DWT变换对： (7.3.5) (7.3.6) 对于j≥j0还有： ? (7.3.7) 这里，f(x)，φjo,k(x)和φj,k(x)是离散变量x=0，1，2，·，M一1的函数。例如，对于某个x0，△x，x=0，1，2,…，M-1，f(x)=f(x0+x△x)。通常冷j0=0并选择M是2的幂(即M=2J)，对x=0，l，2，，M -1，j=0，1，2，…，J-l和k=0，1，2，…，2j-1求和。对于哈尔小波，变换中使用的离散尺度和小波函数(即基函数)与7.l.3节中的MM哈尔变换矩阵的行相对应。变换本身由M个系数组成，最小尺度是0，最大尺度是J-1。出于7.3.1节和例7.6中解释的原因，式(7.3.5)和式(7.3.6)中定义的系数分别称为近似值和细节系数。 式(7.3.5)到式(7.3.7)中的Wφ(j0，k)和Wφ(j，k)对应于前节的小波序列展开中的cj0(k)和dj(k)(这些变量的变化并不是必须的，但它为下一节的连续小 波变换中的标准符号提供了铺垫)。注意，序列展开中的积分变成了求和，而曾在4.2.1节的DFT中出现的归一化因子1/展开和反展开表达式中都有出现。该因子也可以在正展开和反展开表达式中以1／M的形式出现。最后，应该记住，式(7.3.5)到式(7.3.7)只对正交基和紧框架有效。对于双正交基，式(7.3.5)和式(7.3.6)中的φ和φ项必须由它们的对偶函数和代替。 　 例 7.8 计算一维离散小波变换 为说明式(7.3.5)和式(7.3.7)的用途,考虑四点的离散函数:f(0)=1,f(1)=4,f(2)=-3和f(3)=0.因为M=4.J=2且由于j0=0,对x=0,1,2,3(且对于j=0,k=0,或者对于j=1,k=0,1)求和.将使用哈尔尺度和小波函数,并假定f9x)的4个采样值分布在基函数的支撑区上,基函数的值为1.将4个采样点代入式(7.3.5),可得: 因为对于x=0,1,2,3,φ0,0(x)=1.注意,这里采用的是哈尔尺度函数对于j=0且k=0的均匀空间采样.该值对应于7.1.3节的哈尔变换矩阵H4的第一行.继续使用式(7.3.6)和相似间隔的采样点φj,k(x)(它对应于H4的第2,3,4行),可得 因此,这种简单的四点采样函数