第2章小波变换的基础理论研究

1、精品文档第2章小波变换的基础理论研究2.1引言Fourier变换是信号处理的重要工具，它在语音、雷达、声纳、地震、图像、通信系统、自动控制、生物医学工程、机械振动、遥感遥测、电力系统等许多领域都得到了应用。但是Fourier变换反映的是信号的整体特性，不能得到信号的局部特性。小波变换是时间和频率的局部变换，能更加有效的提取和分析信号的局部特性。小波分析在许多领域，如信号分析、图像识别、计算机视觉、视频图像分析、数据压缩和传输、故障诊断等领域有着重要的应用。图像的小波变换及其压缩编码是当今图像压缩编码领域中的研究热点，Saprio1提出的图像内嵌小波零树EZW(EmbeddedZerotreeW

2、avelet)编码技术使该技术由理论走向实用。Said和Pearlman2在EZW基础上又给出了更为精细的基于分层树结构的集合分裂算法SPIHT(SetPartitioninginHierarchicalTree)。由于这种算法既降低了码率，又加快了算法的执行速度，因而得到了广泛的应用，在新标准JPEG20003,4中已有应用。然而在图像的小波压缩编码技术中，小波基的正则性、信号的边界处理对压缩率及失真的影响仍然是值得研究的问题5。本章在小波变换的理论基础上对这些问题做了一些分析和探讨。2.2小波变换的原理7-112.2.1 短时Fourier变换(STFT:ShortTimeFourierT

3、ransform)对于一些非平稳信号，如音乐信号、语音信号，图像信号等，它们的频域特性都是随时间变化的。对这一类信号用Fourier变换进行分析，仅能知道信号所含有的频率信息，但不能知道这些频率信息究竟出现在什么时段上，为了研究这些信号的局部形态，需要对信号进行二维时-频分析。二维时-频分析实际上就是依赖于时间的频谱特性。STFT首先是由Garbor提出的。考虑一个信号x(t),集中在一个局部点J假定通过一个窗函数g(t)(g(t)是在有限的时间区间内定义的)，通过窗口的x(t)g”(t-9的Fourier变换就是STFT7:STFTx(f)=£x(t)g(t-力ej2ftdt(2.

4、1)这就把一维信号x经STFT变换映射到二维的时-频平面(小f)上。STFT非常强的依赖于窗函数g(t)oSTFT的反变换7为：Mt)=,Jdf.JSTFXTf)g(t-加j2ftdp(2.2)(2.(1) 从两方面进行解释6,7。在二维时-频平面上，如图2.1所示，垂直方水平方向的横条表示在给定频率f处，用 f脉冲响应为g的带通滤波器对所有时间 的信号进行滤波。g(t)不但要求在时域是近 似有限宽的，而且在频域也要求是近似有 限宽的，如图2.2所示。STFT中的窗函数g(t)一旦确定，它的 时间窗大小和频率窗大小就确定了。时间 窗大小和频率窗大小决定了g的时间分辨率和频率分辨率。根据文献7,

5、设g(t)fSTFT(f f，2STFT(力2)图2.1 STFT的时-频平面向上的竖条表示在出寸刻，在窗函数g(t)确定的窗口范围内所含有的所有频率分量；另一方面，从子带分析的角度来看，精品文档的Fourier变换为G，定义滤波器g(t)的带宽M为:(2.3)八f2|G(f)|2df-2|G(f)|2df其中分母表示g(t)的能量。两个正弦信号，只有它们的频率差大于3时，通过滤波器g(t)才能将它们区分出来，所以以g(t)为窗口的STFT的频率分辨率由由决定。类似，时间的宽度3定义为：22(2.4)废_t|g(t)|dtIg(t)12dt其中分母同样表示g(t)的能量，两个脉冲信号只有它们在

6、时间上相距大于加时，通过滤波器G才能将它们区分出来。在(2.3),(2.4)式中，假定g(t)和G(f)的中心点在t=0和f=0处，如图2.2所示在STFT中，由于A和Af的固定不变，在整个时-频平面上只能采用相同的频率、时间分辨率，这是STFT的不足。因为对于非平稳信号，也许某一小时间段上，是以高频信息为主，我们希望用小的时间窗进行分析，而对长时间段上的低频信号，希望采用大时间窗进行分析。因此，对于一个时变的非平稳信号，很难找到一个好的时间窗来适合于不同的时间段。小波变换的引入弥补了STFT的不足。小波变换采用了可变带宽的窗函数(滤波器)，在低频端用窄带滤波器进行分析，而在高频端用宽带滤波器

7、进行分析，这就是所谓的相对带宽固定的滤波器组(即恒Q特性)6。如图2.3(b)所示。图2.3为STFT下的绝对带宽恒定的滤波器组。G(f)8fo(b)相对带宽恒定的滤波器组频谱图2.3STFT和小波变换滤波器组的频谱2.2.2连续小波变换(CWT:ContinuousWaveletTransform)2.2.2.1 一维连续小波变换设dt)为平方可积函数，Mt)wL2(R),若其Fouriera变换A(co)满足:(2.5)则称a(t)为一个基本小波或小波母函数，(2.5)式为小波函数的容许条件。由(2.5)式可知：A(=0=fo(t)dt=0(2.6)R即小波母函数的均值为零，那么它一定是正

8、负交替的。如Marr小波：.22一t/2223/2阳)=e(1t)宣3)=可2九3e对一个小波母函数0<t）,通过平移和伸缩构成一组小波基，记为%,0）1 t- T%,p(t) =4).a aa>0, R(2.7)参数a,吩别称为尺度因子和位移因子（常数1/西是用于能量的归一化）。Marr小波的基及其Fourier变换如图2.4所示。由于a,t为连续变换的，所以看,小）为连续小波基函数图2.4 Marr小波基的时域波形、频域波形信号x（t）的连续小波变换定义为7CWTx(a,Xf(t),«a,T(t)>(2.8)1,tT=丁仅。)蚁)dt.aRaCWT将一维信号映射

9、到二维时间-尺度平面上，在二维时间-尺度平面上,有利于信号特征的提取。从时-频分析的角度来看，如令(2.9)则CWT可看作是STFTo信号x（t）在某一尺度a,平移点三上的小波变换系数,实质上表征的是7位置上，时间段aN上经过中心频率为的，带宽为4相带通滤波器的频率分量大小。随着尺度a的变化，带通滤波器的中心频率及带宽都发生变化。当分析低频（对应大尺度）信号时，其时间窗增大，滤波器中心频率和带宽减小，而当分析高频（对应小尺度）信号时，其时间窗减小，滤波器中心频率和带宽增大，这正好符合实际问题中高频信号持续时间短，低频信号持续时间长的自然规律。而在STFT中，窗口是固定不变的，这正是两者的本质区

10、别。小波变换的逆变换为8,9:x(t)=d|JCWT(a,如aMt)d忑(2.10)Chr-ar2.2.2.2 二维连续小波变换设a(x,y)为一二维连续函数，满足容许条件：- RX w d(2.12)(2.13)图2.5二维Marr小波母函数的时域波形，尺度 a分别为0.5、1、2则o(x,y)可以作为小波母函数。庆(附,孙)为(x,y)的二维Fourier变换。维连续函数f(x,y)的小波变换定义为:CWT(a,3©)="f(x,y)%3©(x,y)dxdyRR其中，,彳,为(x,y)为二维小波基，专，及为两个方向上的位移，a为尺度1x1y百洵q,及(x,y)

11、=；-'队,)|a|aa二维Marr小波母函数o<x,y)=(1-x2-y2)e-(x24y2)/2，时域波形见图2.5维小波变换的逆变换为:f (x,y) = Jo 粤二Ch a容CWTf (a, a,T2)dcd 及当二维小波母函数是可分离型时，即a(x,y)=g(x)a2(y)，则它可简化为一维小波变换。在实际的图像小波变换中大都采用可分离的小波变换02.3多分辨率分析与Mallat算法2.3.1 小波变换参数的离散化由于连续小波变换CWTx(a,)变换域参数是连续的，从降低信息冗余9的角度和实际应用的角度来说，需要将尺度参数和位移参数离散化。一种最常用的方法就是将尺度按幕

12、级数进行离散化，即取a=a°m(mwZ,Z为整数基，a/1,常取为2)。位移的离散化，为了不丢失信息，要满足Nyquist采样定理。当尺度a增加一倍时(m加1),对应的滤波器带宽减小一半，采样频率可以降低一半，采样间隔增大一倍。因此，在尺度a=1(m=0)时，位移e的采样间隔设为To,则在尺度a=aom时的采样间隔为aomToo因此，在尺度和位移都离散化后，小波基函数Ota7可表示为：a,m_m2.t2nT0-m2_m.2Xm)=2o(2t-nTo)->om,n(t),m,n=Z.2在归一化(设To=1)情况下，上式为：n(t)=2-m/242-mtn)(2.15)任意函数x的

13、离散小波变换为：DWTx(m,n)二亡x(t)淅,n(t)dt(2.16)2.3.2 多分辨率分析(MRA:Multi-ResolutionAnalysis)多分辨率分析方法是Mallat在研究计算机视觉时提出的11,它的基本思想是将图像在不同尺度下分解，得到不同尺度下图像分解的结果，然后进行比较，从而得到一些有用的信息。2.3.2.1 多分辨率分析的数学描述9,10,112设函数幻)WL(R)(平方可积空间)，若其整数平移序列%(t)=淋-n)相互正交，即：<%(t),%,(t)A9nn')n,nZ(2.17)则由%(t)所张成的子空间称为尺度空间网，而函数口(t)称为尺度函数

14、(或生成元)。由(2.15)式可知，在尺度函数序列%(t)中由于m=0,因此，由%(t)所张成的子空间为零尺度空间，记做V。，而厮(t)即为V0的一组基根据泛函分析的理论12,任意函数f(t)wVo,可以由Vo的一组基的线性表出，即：f(t)=£cn%(t)(2.18)n其中：Cn=£f(t)生(t)dt(2.19)同样可得到尺度m#0下的尺度函数序列mn(t)=25'22-mt-n),由,n(t)所张成的子空间为m尺度空间，记为Vmo那么任意f(t)WVm可由%,n(t)线性表出：f(t)=£cm,ngn(t)(2.20)n由此，尺度函数必)在不同尺度下

15、的平移序列张成了一系列的尺度空间Vm,m-Zo随着尺度m的增大，函数.n(t)的宽度增大，且实际的平移间隔(2mTo)也变大，所以它的线性表达式(2.20)就不能表示函数的细微变换(小于该尺度下的变化)，因此，其张成的尺度空间只能包含大尺度的慢变信号，相反，随着尺度m的减小，函数0n(t)的宽度变小，实际的平移间隔(2mTo)也变小，则它的线性表达式可以表达函数的更细微的变化，因此，其张成的尺度空间所包含的函数增多，尺度空间变大。由此，可以给出多分辨率分析严格的数学描述：(1)在L2(R)中，存在一系列嵌套子空间Vm,mZ,二V2二V1二V0二V二V工二(2.21)这一系列嵌套子空间具有：逼近

16、性：riVm=(0)UVm=L2(R)(2.22)mEZm之伸缩性：f(t)Vm=f(2t)Vm_1(2.23)(2)存在函数Mt)WV0,使得与=Mt-n)n已构成Vo的正交基，即V0=spanf加一n),o(tn)如m)dt=8(m-n)(2.24)R若厮(t)ng是Vo的正交基，则qm,n(t)=2f242-mt-n)m,ng是Vm的正交基。2.3.2.2 小波函数10虽然多分辨率分析的一系列子空间逼近L2(R),但是，由于它们之间是互相包含的，不具有正交性，因此它们的基0,n(t)m,nEZ在不同尺度下不具有正交性，因而也就不能作为L2(R)的正交基。为了寻找一组L2(R)的正交基，有

17、必要引入Vm的正交补。设Wm是Vm的在Vm-1中的正交补空间，即Vm-Vm二WmWmVVm(2.25)那么，对任意mF,Wm和Wn都是正交的，由(2.21)、(2.22)式可得：2L(R)=二Wm(2.26)m.-Z因此，Wm是构成L2(R)的一系列正交子空间，且Wm=Vm-Wm,W0=V-V0(2.27)若g(t)wW0,则g(2-mt)Wmlm,亦即：g(t)W0=g(2t)Wm(2.28)若露(t)=B(t-n)n且是性的一组正交基，由(2.28)式对任意尺度mZ,阶n(t)=2,2H2-mt-n)nez一定是Wm的一组正交基，再根据(2.26)式，全体他,n(t)m,na构成L2(R)

18、的一组正交基，阳)就是小波母函数，Wm是尺度为m的小波空间。小波空间与尺度空间是互补的，尺度空间之间是包含关系，而小波空间是正交关系。2.3.2.3 一维信号的多分辨率分析根据多分辨率分析的定义，由于V0=Vi份Wi,如果一维信号f(t)WV0,则f(t)可分解(投影)为Vi和Wi上的两部分，在Vi上的投影称为f(t)的近似部分，记为f1a(t),在Wi上的投影称为f(t)的细节部分，记为f1d(t)。如果Mt)是尺度函数，F(t)是小波函数，那么：fia(t)=Nai,nM2-it-n)(2.29)nfid(t)=£din«2-it-n)(2.30)n其中ai,n和di,

19、n分别为尺度分解系数和小波分解系数：ai,n=：二f(t),2-it-n)(2.3i)di,nxf(t),3(2-it-n)>(2.32)重构信号：f(t)=f1a(t)f1d(t)_-1_1=£ai,nM2tn)+£di,n*2t-n)(2.33)nn同样，当f(t)wV-i,V=V0Wo,则f(t)可分解(投影)为Vo和Wo上的两部分：f(t)=foa(t)fod(t)=£ao,nMt-n)+£d°,n阳-n)(2.34)nnao,n=：二f(t),Kt-n).(2.35)do,nKf(t),atn)A(2.36)以上的分析是在m=-

20、1,0,1的尺度空间进行的，在其它尺度空间上同样适用。2.3.2.4 一维二尺度方程尺度函数和小波函数在相邻两个尺度上的关系就是二尺度方程，它反映了相邻尺度空间和小波空间之间的内在联系。由多分辨率分析的定义可知，若Mt)和阳)分别为尺度空间Vo和小波空间Wo的正交基，由于VoV-1,W0UV-1,所以染)和阳)也必然包含在V-1中，而V-1的一组正交基为乜n(t)=21/2戏t-n)nwz,所以必)和阳)可以表示为:如)=21/2%h(n)o(2t-n)(2.37)n1/2B(t)=2'g(n)o<2tn)(2.38)n其中系数h(n)和g(n)分别为：h(n)k如)，a,n(t

21、)>=21/2/加)W2t-n)dt(2.39)Rg(n)K阳)，a,n(t)a=21/2J阳)W2tn)dt(2.40)R对任意相邻空间，(2.37)和(2.38)都成立10,此二式就称为二尺度方程。系数h(n)和g(n)称作滤波器系数，它们是由尺度函数项)和小波函数阳)决定的，与具体的尺度无关。2.3.3 Mallat算法11设V0是一尺度空间，%(t)=必-n)n-Z是它的正交基，那么任意f(t)wV0可以由«n(t)n.线性表出，即：f(t)=£a0,no(tn)n其中：a0,nKf(t),a(t_n)A为f(t)在V0上的尺度系数。由于V0=V19W1,将f

22、(t)在Vi和Wi上分解,即：-，-一1、，-1、f(t)=£ai,nM2t-n)+Zdi,n做2t-n)nn其中："n(t)=2'/22'tn)nx是Vi的正交基，§,n(t)=2/2氐2,tn)nsz是Wi的正交基。atn=<f(t),«(2-it-n)>,d1nxf(t),H2-itn)>分别为f(t)在Vi上的尺度系数和Wi上的小波系数。根据二尺度方程不难得到i0,ii：ai,n=，h(m-2n)ao,m="h(m)ao,2nm(2.4i)mmdi,n="g(m2n)a0,m="g(

23、m)a0,2nm(2.42)mm其中：h(m)、g(m)分别是分解时相应的低通滤波器和高通滤波器，a0,m是0尺度下尺度系数，ai,m是i尺度的尺度系数，dim是i尺度的小波系数，初始系数a0,m一般取为函数的采样值i0o对aim继续分解就可以得到函数的金字塔型分解算法，如图2.6所示图2.6维Mallat分解算法重构：a0,m=：h(m-2k)ai,k八(m-2k)di,kkk=h(m)ai,(m4)/2g(m)di,(m*)/2(2.43)mm一一一、其中：h(m)、(m)分别是重构时相应的低通滤波器和图通滤波器。上式中ai,(m-k)/2、di,(m-k)/2当(m-k)取奇数时取零,即

24、插零2.4图像信号的小波变换将一维信号的分析结果推广到二维图像信号上来，就有以下的结果2.4.1 二维多分辨率分析设Vmm字是一维L2(R)的多分辨率分析的尺度空间嵌套，o(t)是它的尺度函数。Vm的一组正交基为om,n(t)=2；2戏-mt-n),Wn是Vm在Vm-1中的正交补，即Vm=VmWm,Wm1Vm,而外的一组正交基为3,n(t)=252乱2mt-n)。则可通过一维多分辨率分析的张量积来定义二维的多分辨率分析。设：以乂二%(2.50)则Vm满足：uV1cV0UV/U(2.51)一Vm=L2(R2),-Vm40)(2.51)mFZm三Z令(x,y)=0(x)0),则(x,y)为Vm的尺

25、度函数，那么，Vm空间的一组正交基为：m,ni,n2(x,y)=2a(2mx-n1)of2my-n2)m,ni,n2Z(2.52)Vm-Vmj二Vm=(Vm二Wm)二(Vm二Wm)123叫二Wm1二Wm2二Wm3(2.53)其中：Wm1=Vm®Wm，Wm2=Wm®Vm，Wm3=Wm®Wm(2.54)iijii且Vm!Wm(i=1,2,3),Wm!WmJiwj,i,jw1,2,3,及WmVW(miWm!2)。在m尺度上，有四个正交的子空间Vm,Wm1,Wm2,Wm3,它们的母函数分别是：(x,y)=o(x)o(y)(2.55a)&x,y)=o(x)B(y)(

26、2.55b)"x,y)=gx)a(y)(2.55c)n(x,y)=B(x)B(y)(2.55d)由母函数的伸缩和平移可构成它们的正交基，类似于(2.52)式。由(2.53)式，任一函数f(x,y)CV0可投影到V；,W1,W12和W13空间上，则有:f(x，y)=ZZa1,nn21,n1,n2(x，y)d1nl,n24,n1,n2(x，y)nin2n1n2一"孤一1,ni,n2(X,y)一二二d；ni,n21,ni,n2(X,y)(2.56)nin2nin2其中：ai,n1m281门口20=1,2,3)为展开系数，分别为：ai,n1,n2Kf(X,y),-：J1,n1,n2

27、(X,y)(2.57a)1di,nin2=f(x,y),i,ni,n2(x,y)(2.57b)d1,ni,n2=<f(X,y),i,ni,n2(X,y)(2.57c).3.di,nin2=：二f(X,y),i,n1,n2(X,y)(2.57d)ainz为f(X,y)在尺度1上的近似系数，dimm为该尺度三个方向上的细节系数。2.4.2图像小波变换的算法2.4.2.1 二维二尺度方程在二维情况下，每个尺度上有四个子空间，所以就有四个二尺度方程，在0尺度和1尺度之间二尺度方程为：O(x,y)=江ho(ni,血)Oi,m,n2(X,y)(2.58a)nin2&x,y)=£&#

28、163;hi(ni,ri2)3,小股(x,y)(2.58b)nin2(x,y)三年(5,6)im/xy)(2.58c)nin2(x,y)f卜3(小币2)i,n1,n2(x,y)(2.58d)nin2其中hi(X,y)i=0,1,2,3二维滤波器系数。对任意两相邻子空间Vm和Vm由，上面的二尺度方程总是成立。在(2.50)式的条件下，可以得到：ho(ni,n2)k(x,y),i,m,n(x,y)>(2.59a)=°(x)o(y)2o(2x-n1)o(2y-n2)dxdy=h(ni)h(ri2)，/、产/、产z、一hi(ni,n2)(x,y),im,n(x,y)=自己<x)R

29、y)2o(2x-ni)R2y-n2)dxdy=h(ni)g(n2)(2.59b)其中h(ni),g(nj为一维滤波器,见(2.39)、(2.40)式。同理：h2(ni,n2)=g(ni)h(n2)(2.59c)h3(ni,n2)=g(ni)g(n2)(2.59d)2.4.2.2 图像小波变换的分解与重构算法10,11根据二尺度方程(2.58)式及(2.57)式可以得到二维图像信号的小波分解算法为：ai,k1,k2h(ni)h(n2)a0,n12klm22k2(2.60a)n1n2I1d1,k1,k2=h(n1)g(n2)a0,n12k1,n22k2(2.60b)n1n212d1,k1,k2=、

30、g(n1)h(n2)a0,n12kl,n22k2(2.60c)n1n2I3d1R,k2=g(n1)g(n2)a0,n12kle22k2(2.60d)n1n2上式将0尺度上的近似系数a。分解为1尺度上的近似系数和三个不同方向上的细节系数。对1维尺度上的近似系数a1再进行分解，即形成了二维的金字塔分解算法。1重构：a0,n1,n2=h(k1)h(k2)a1,(n1-k1)/2,(n2-k2)/2+-h(k1)g(k2)d1(n1-k1)/2(n2-k2)/2k1k2k1收十二g(k1)h(k2)d1(n14)/2-k2)/2+-9(4)9«2)&1(*)/2,(n2抬)/2(2.

31、61)k1k2k1检2.5 小波变换的正则性132.5.1 小波变换的滤波器表示在小波变换算法中，滤波器起着重要的作用。一维信号的一级小波分解与重构可以用图2.7所示的双通道的滤波器组来表示。图2.7双通道的滤波器组设h(k)和g(k)的z变换分别为H(z)和G(z),出)和%)的z变换分别为才和弓，a0,k和a0,k分别记为a0(k)和？0(k),它们的z变换分别为Ao(z)和及,则：Ao(z)=1H(z)(z)G(z)(5(z)A(z)21+2H(-z)H(z)+G(-z)G(z)A(-z)(2.62)当滤波器满足以下条件时：H(z)-G(-z)(2.63)G(z)=H(-z)(2.64)

32、H(z)G(-z)-H(-z)G(z)-2z-n(2.65)则：A0(z)=zHAo(z)(2.66)对应时域为：?o(k)=a0(k-n)(2.67)此式表明，重构信号仅是输入信号的n步延迟，即实现了完全重构。由下面的公式可知n一定为奇数。而(2.63),(2.64),(2.65)式就是完全重构条件，对应到时域为：h(n)-(-1)g(n)g(n)=(-1)h(n)iLk2j=n1-(-1)jS(-1)kh(k)g(k-j)=t，JkPj#n2.5.2 正则性探讨5,13在图像的小波变换编码中，一般选择正交或双正交小波基。小波基对应滤波器组，不同的滤波器组在进行信号分解与重构时产生不同的效果

33、。小波基与滤波器组的对应关系见表2.113:表2.1:小波基与滤波器组的对应关系小波基滤波器组对称性相位系数数量正则性正交小波不对称(除一组系数，其他三基不保证线性相位组系数是该组系数的逆序及其调制任意阶正则性Haar小波)双正交小两组系数，其他两有限阶正则性对称确保线性相位组系数是这两组系(与滤波器长度波基数的逆序及其调制有关)正则性是函数光滑程度的反映，是选择小波变换滤波器应考虑的一个重要因素。定义1:设函数o(x)的N阶导数存在，对任意x,殳R,若有a(N)(xM)a(N)(x)|<丁产，其中：0<K<1,C是与x9无关的常数，那么称o(x)具有N+上阶正则性。根据定义

34、，若Mx)具有N+K阶正则性，则o(x)具有N阶连续导数。正则性越高，函数越光滑。越光滑的函数，在频域中能量越集中。尺度函数的正则性与小波函数的消失矩成对应关系，高阶正则性的尺度函数对应高阶消失矩的小波函数。当给出分解低通滤波器h(n)后，可以构造出尺度函数o(x严,15。定理1:设h(n)是有限长序列，n=0,1,L-1,H(z)是h(n)的z变换。定义H4(z)=H(z2i),H(0(z)=1。H(k(z)所对应的时域序列为h(k*n),则h(k*n)i=0的长度为(2k-1)(L-1)+1。证明：用数学归纳法k=1时，HC(z)=H(z),对应的时间序列为h0(n),结论成立。八八cL/

35、ck=2时，H(z)=H(z)H(z),H(z)=2h(n)zn,对应的时间序列长度nHD为2(L-1)十1。根据卷积定理，H(2b)的时间序列为H(z)与H(z2)时间序列之卷积，其长度为2(L-1)+1+L-1=3(L-1)+1，结论成立。设k=K时，结论成立。当k=K+1时：H"*z)=H(K(z)H(z2K)由于H(z2K)的时间序列长度为2K(L-1)+1,所以H(K*h)的序列长度为(2K-1)(L-1)+1+2K(L-1)+1=(2K*)(L-1)+1。结论成立。定理2:构造一个函数：f(k)(t)=2k/2h(k)(n),n/2k4t4(n+1)/2k,h(k)(n)

36、与H(k)(z)为一Z变换对。则：f(k)(t)=21/2Sh(m)f(k-1)(2t-m),f(k)(t)的支撑集m为0MtML-1。证明：由于H(k)(z)=H(k,)(z)H(z2kj,所以：f=2k/2£h(m)h(k/)(n-2k-1m),n/2k<t<(n+1)/2k。m又由于:f2(2")=2(kJ)/2h(kJ)(n),n/2k,«2t-m«(n+1)/2k。令：n1=n+2"m,则：f(kJ)(2t-m)=2(kJ)/2h(kJ)(n1-2k-1m),n1/2k<t<(n1+1)/2k0所以：f(k)(

37、t)=21/2£h(m)f(k-1)(2t-m)。m根据定理1,0<t<L-10当kT®时，若f(k)(t)收敛到一个连续函数o(x),那么h(n)就具有正则性，Mx)就是由h(n)生成的尺度函数。只有h(n)满足一定的条件16,才能保证f(k)(x)收敛到一个连续函数。实际的图像除小部分边沿之外，大部分是光滑的。小波分解时，低频系数包含原图像的大量信息，而高频系数除极少数系数外，大部分接近于零。基于小波变换的图像压缩编码技术正是利用了小波变换的这个特性。设f(t)是一个充分光滑的函数，f(t)可由t=0点的Taylor级数近似：f(t)：Pm(t)=f(0)f

38、'(t)f-(0)t2-f-(0)tM(2.68)2!M!其中：f(t)表示f(t)在t=0点的i阶导数。根据Taylor级数的性质，当乂*时f'(0)t0。M!由(2-36)式，小波展开系数do.xf(t),附-n)a.巴Pm(t)虞-n)dtMfi(0)-i=£.gtB(t-n)dt(2.69)i=0i!定理3:当小波函数B(t)具有M阶消失矩时，即：巴,ti黑)dt=0,i=0,1，,M-1,则：邕ti&t-n)dt=0,i=0,1，,M-1,VnwZ(整数集)。i证明：Jti第-n)dt=g(t+n)i第)dt=gZGjtjni-j阳)dtj-j=0i

39、=£Cijni-j自tjB(t)dt=0j/一一回到(2.69)式，当不)的消失矩越高时，d0,n越接近于零。2.6 滤波器组的线性相位与边界处理13当信号经过小波分解后，分解后的信号要发生相移；当滤波器具有线性相位时，它的相移是线性的；而当滤波器不具有线性相位时，它的相移是非线性的。对线性相移，在重构时可用与分解滤波器成正交的重构滤波器进行补偿，使重构的信号在相位上保证与原信号一致。小波变换是针对无限长信号的，对于有限长的信号，必须处理好边界问题。处理边界的原则是(1)保持信息不丢失(即恢复信号不能畸变)，(2)数据量不能增大。为此，已有了零延拓、对称延拓17及周期延拓。我们认为周

40、期延拓是更好的方法，结合本文提出的滤波方法，可以满足上面两点要求。设待处理的有限长信号为x(0),x(1)，x(N-1),长度为N;对x(n)做周期延拓。分解时向后延拓M(滤波器长度)点，xf(n)=x(-M),x(-1),x(0)，x(N-1),其中x(-1)=x(N-1),x(-2)=x(N-2),，对x,(n)做滤波，设滤波器为h0(0),h0(1)，h°(M-1),则输出：y(n)="h0(m)xf(n-m)n=0,1，,N-1(2.70)m表示为矩阵形式：y(0)'n%"x(N -1)-一 0t(0)%(0)h0(M -1f：0h0(0)hi(1

41、)h0(M-1)0idd. d3h0(0)hi(1)h0(M -1)0x(0)x(-1)_x(M)(2.71)y(n)的相位比x(n)滞后一个角度，长度仍保持N。重构时则相反，将输入信号y(n)向前延拓M点，滤波器采用重构滤波器g0(m),则重构信号为？(n),用矩阵表示为：一？(0) g0(M -1)g0(1)?(1)0g0(M -1)- + * I+ :?(N1)_ g0(M -1)g0(0) g00g0(0)+g0(1)y(0)0y：N-1)： y(N)Bg0(0)： yN+Mf-(2.72)重构信号？(n)的相位比y(n)超前其中：y(N)=y(0),y(N+1)=y(1)，。一个角度

42、，从而使？(n)与x(n)保持同相位，长度为N2.7 仿真研究2.7.1 滤波器正则性与小波变换的关系正则性与小波变换的关系。选用不同正则性的滤波器组对标准图像Lena(256M256M8bit)及mosaic(256256x8bit)图像进行一次分解与重构，如图2.8所示，结果见表2.2,其中：a。为原图像的像素灰度值矩阵，a1为低频分量，dk图2.8图像的一次小波分解(k=1,2,3)为三个方向上的高频分量12822562ra="(a1,i,j)/(a0,i,j)i,ji,j128k22562rd=£(di,j)/Z(a0,i,j)。i,ji,j图2.9、图2.10为重

43、构图像。重构时只保留低频分量，并对其进行整数量化,舍去高频分量。表2.2正则性与小波变换的关系(Lena图像/mosaic图像)_4一正则性rard1rd2,3rdPSNR(db)波长度01二2二3r(h0)2.707,.7070.000.9941.0035.0020.000328.77aH/(g0)2.707,.7070.000.6250.1250.1250.12507.27交(h0)9.853,.377,-.111,-.024,.0381.068.9966.0012.0006.000233.25正/双(g0)7.788,.418,-.041,-.0651.701.8001.0632.063

44、2.063211.562.7.2 滤波器的线性相位与信号的边界问题选输入信号为正弦信号x(n)=sin(sn),频率切=兀/18,y1为分解后的低频信号,y为重构信号(舍去了高频信号)。(a)原图像(b)Haar小波(c)双正交小波图2.9Lena图像的小波分解与重构(a)原图像(b)Haar小波(c)双正交小波图2.10Mosaic图像的小波分解与重构Case1选用双正交小波滤波器组h0=0.0663,-0.1989,-0.1547,0.9944,0.9944,-0.1547,-0.1989,0.0663g0=0.0000,0.0000,0.1768,Case2:选用止交小波滤波器组0.53

45、03,0.5303,0.1768,0.0000,0,0000h0=-0.0106,0.0329,0.0308,-0.1870,-0.0280,0.6309,0.7148,0.2304g0=0.2304,0.7148,0.6309,-0.0280,-0.1870,0.0308,0.0329,-0.0106仿真结果见图2.11和图2.122 y0'0 一-1 _ _ _ _051015 (a) 20253035 n 400-2 051015 (b) 20253035 n 402 y一0 -'r-2 051015 20253035 n 40图2.11双正交小波分解与重构(a)输入信号

46、，(b)分解信号，(c)重构信号-2051015 (c) 20253035 n 40图2.12正交小波分解与重构(a)输入信号，(b)分解信号，(c)重构信号本章小结:小波变换是一种变尺度的时-频联合分析方法，它更适合于对非平稳信号的处理，能够分析到信号的局部形态，它的变尺度特性一一低频部分采用大的时间尺度，高频部分采用小的时间尺度，与人类的视觉系统相吻合。本章对小波变换的理论基础做了分析和归纳，特别对小波变换中的两个关键技术一一小波的正则性与信号的边界处理做了进一步的分析和探讨，从而得出如下结论：1 .正则性是我们选择小波基时要考虑的一个重要因素。对一般较为平滑的图像，如Lena图像，经过小波分解后，低频分量较大，而高频分量较小；而且正则性阶数越高，低频分量越大，高频分量越小。对于图像的压缩编码来说，压缩比大，恢复图像失真小。但对纹理较多的图像，如mosaic图像，小波分解后，其低频分量减小，高频分量增大，因而压缩比减小，恢复图像的失真加大。2 .对有限长信号，经过小波分解后，必引起相位的偏移，为了减小相位偏移所引起的相位失真，在分解与重构时应采用向前、向后的周期延拓；采用本文所提出分解与重构算法，可以使信号在小波变换前后的相位偏移得到完全的补偿。经过周期延拓，可以使数据量不扩