小波分析论文

1、多分辨率分析和Mallat算法摘要：多分辨率是小波优于傅里叶变换的重要特性，使小波成为“数学显微镜”，多分辨率分析将小波基的构造统一起来。Mallat算法是计算离散栅格上小波变换的快速算法，它的出现使小波从理论上升的实践。本文主要对多分辨率分析和Mallat算法进行分析，并给出其应用仿真实例。关键词：多分辨率分析Mallat算法应用仿真实例1引言小波变换的多分辨率分析是建立在函数概念上的理论，多分辨率分析(MultiresolutionAnalysis，MRA)概念是由S.Mallat和Y.Meyer在前人大量工作的基础上于1986年提出的，从空间的概念上形象的说明了小波的多分辨率特性，随着尺

2、度由大到小变化，在各尺度上可以由粗到细的观察图像的不同特征。在大尺度时，观察到图像的轮廓，在小尺度的空间里，则可以观察图像的细节。1989年，Mallat在小波变换多分辨率分析理论与图像处理的应用研究中受到塔式算法的启发，提出了信号的塔式多分辨率分析与重构的快速算法称为马拉特(Mallat)算法。多分辨率分析和Mallat算法构成了小波基础理论的重要板块，多分辨率分析从理论上说明了小波的自适应变频特性，也是构造小波滤波器的基础，Mallat算法使小波从理论走向工程，其功绩相当于FFT在傅里叶变换中的作用。2小波基本概念定义2.1小波的定义。定义时域具有紧支集(持续时间短)且平均值为零(不含直流

3、分量)的函数为小波，通常用屮(t)表示。定义式为j"屮C)dt=0g通常，屮(t)称为基本小波函数、母小波函数或小波母函数，简称小波函数、母小波等。由于上述定义是十分宽泛的，并不是所有满足定义的函数屮(t)都是具有价值的小波，若不加限制条件其逆变换也不一定存在，所以小波必须施加以限制条件，称为容许条件。定义2.2小波容许条件。假设屮(t)为小波母函数，则其傅里叶变换屮()应满足厂fb(®)2刁一C=fdx<g屮R3上述条件被称为小波函数容许条件。定义2.3连续小波变换定义。假设f(t)gL2(R)，屮(t)是连续小波母函数，屮(t)a,b为连续小波基函数，则小波变换可

4、定义为:1gtbWT(a,b)=ff(t”*()dtgfaa离散小波定义即对尺度因子a和平移因子b分别离散化。3多分辨率分析和Mallat算法3.1多分辨率分析多分辨率分析首先并不是应用于小波分析上，它曾成功应用于地球物理学，S.Mallat和Y.Meyer在前人大量工作的基础上于1986年把多分辨率分析应用到小波上，统一了小波基的构造，有人说只要掌握多尺度就可以自己构造自己的小波。分析空间L2(R)中的多分辨率分析是指L2(R)中满足下列条件的一个序列.。jjez单调性：对任意jeZ，有e3逼近性：nv二£,Jv二L2(r)。jjjeZj=s伸缩性体现了尺度的变换、逼近正交小波函伸

5、缩性：f(t)eVof(2t)eV，jj1数的变化和空间的变化具有一致性。(平移不变性：对任意keZ，有2NteVn2Nt-keV。t丿Riesz基存在性：存在(t)eV。,>构成V的Riesz基。jkeZ3.2Mallat算法分解算法：假定信号f(x)eV，即jf(x)=Ya申(x)j，kj，kk其中系数，kez为已知。如果我们希望将它分解为在j和wji空间的两个分量之和，即f(x)二工ap(x)+工d屮(x)j1，kj1，kj1，kj1，k也就是由已知序列a分别求出j1级的近似序列a和/-1级细节序列d。j，kj1，kj1，k重构算法：当序列a和d为已知时，重构信号f(x)在V空间中

6、的分量A(x)就j,kj,kj+1j+1是利用这两个序列来计算序列a。YCp+dp)pj,lj,ljljlj+1,k八ji，申j+i,k=工a®,申；+工dj厂j,lj+1,k:l四：仿真程序与仿真结果4.1一维信号的Mallat算法仿真程序与仿真结果程序：%一维信号的分解与重构loadnoisblocx=noisbloc;l=length(x);subplot(3,l,l);plot(x);title('原始信号')%信号的一级分解与重构swa,swd=swt(x,1,'db5');subplot(3,l,2);plot(swa);title(

7、9;低频');subplot(3,l,3);plot(swd);title('高频');A0=iswt(swa,swd,'db5');err=norm(x-A0);figure(2);subplot(3,l,l);plot(x);title('原始信号');subplot(3,l,2);plot(A0);title('重构信号');subplot(3,l,3);plot(x-A0);title('误差信号');%信号的单支重构nulcfs=zeros(size(swa);A1=iswt(swa,nulcfs

8、,'db5');D1=iswt(nulcfs,swd,'db5');figure(3);subplot(2,1,1);plot(A1);title('重构低频');subplot(2,1,2);plot(D1);title('重构高频');%信号的多级分解swa,swd=swt(x,3,'db5');kp=0;figure(4);fori=1:3subplot(3,2,kp+1),plot(swa(i,:);title('低频系数:level',num2str(i);subplot(3,2,kp+2

9、),plot(swd(i,:);title('高频系数:level',num2str(i);kp=kp+2;end仿真结果：原始信号低频高频原始信号30重构低频20100ft;、:冷尹心f、2004006008001000卄忖心沁A呼Z1200-5重构高频020040060080010001200500低频系数:level150010001500100100-100500100015000心飞-100''050010001500高频系数:level3200-200500100015004.2二维维信号的Mallat算法仿真程序与仿真结果程序：%二维信号的分解与重

10、构loadbelmont2;nbcol=size(map,1);%图像的多尺度小波分解C,S=wavedec2(X,2,'bior3.7');%提取系数的低频和高频部分cA2=appcoef2(C,S,'bior3.7',2);%单支重构A2=wrcoef2('a',C,S,'bior3.7',2);A1=wrcoef2('a',C,S,'bior3.7',1);H1=wrcoef2('h',C,S,'bior3.7',1);V1=wrcoef2('v'

11、;,C,S,'bior3.7',1);D1=wrcoef2('d',C,S,'bior3.7',1);H2=wrcoef2('h',C,S,'bior3.7',2);V2=wrcoef2('v',C,S,'bior3.7',2);D2=wrcoef2('d',C,S,'bior3.7',2);%显示多尺度分解结果figurecolormap(map);subplot(2,4,l);image(wcodemat(Al,nbcol);title('

12、低频Al');subplot(2,4,2);image(wcodemat(H1,nbcol);title(水平咼频Hl');subplot(2,4,3);image(wcodemat(V1,nbcol);title('垂直高频V1');subplot(2,4,4);image(wcodemat(D1,nbcol);title('对角高频D1');subplot(2,4,5);image(wcodemat(A2,nbcol);title('低频A2');subplot(2,4,6);image(wcodemat(H2,nbcol);

13、title(水平咼频H2');subplot(2,4,7);image(wcodemat(V2,nbcol);title('垂直高频V2');subplot(2,4,8);image(wcodemat(D2,nbcol);title('对角高频D2');%重构原始图像X0=waverec2(C,S,'bior3.7');figure(2);colormap(map);subplot(2,1,1);image(X);title('原始图像')；subplot(2,1,2);image(X0);title('重构图像&

14、#39;)；结果：原始图像10015020050100150200250300重构图像10015020050100150200250300单支分解与重构：低频A1100150200100200300水平高频H1垂直高频V1对角高频D150100150200100200300100200300100200300低频A2100150200100200300水平高频H2垂直高频V250100150200100200300100200300对角高频D2100200300五：结束语通过对小波课程的学习掌握了一门时频分析工具，学会了对一维信号与二维信号的多尺度分解与重构。在学习小波课程中遇到了许多困难，对小波的理解也不够深刻，以后会进一步加强学习。参考文献：1S.Mallat.ATheoryforMultiresolutionSignalDecomposition:TheWaveletRepresentation.IEEETransactionsonPatternAnalysisandMachingIn