小波分析论文..docVIP

多分辨率分析和Mallat算法 摘要：多分辨率是小波优于傅里叶变换的重要特性，使小波成为“数学显微镜”，多分辨率分析将小波基的构造统一起来。Mallat算法是计算离散栅格上小波变换的快速算法，它的出现使小波从理论上升的实践。本文主要对多分辨率分析和Mallat算法进行分析，并给出其应用仿真实例。 关键词：多分辨率分析 Mallat算法 应用 仿真实例 1 引言 小波变换的多分辨率分析是建立在函数概念上的理论，多分辨率分析（Multiresolution Analysis，MRA）概念是由 S.Mallat和 Y.Meyer 在前人大量工作的基础上于1986年提出的，从空间的概念上形象的说明了小波的多分辨率特性，随着尺度由大到小变化，在各尺度上可以由粗到细的观察图像的不同特征。在大尺度时，观察到图像的轮廓，在小尺度的空间里，则可以观察图像的细节。1989年，Mallat在小波变换多分辨率分析理论与图像处理的应用研究中受到塔式算法的启发，提出了信号的塔式多分辨率分析与重构的快速算法称为马拉特（Mallat）算法。多分辨率分析和Mallat算法构成了小波基础理论的重要板块，多分辨率分析从理论上说明了小波的自适应变频特性，也是构造小波滤波器的基础，Mallat算法使小波从理论走向工程，其功绩相当于FFT在傅里叶变换中的作用。 2 小波基本概念 定义2.1 小波的定义。定义时域具有紧支集（持续时间短）且平均值为零（不含直流分量）的函数为小波，通常用表示。定义式为 通常，称为基本小波函数、母小波函数或小波母函数，简称小波函数、母小波等。 由于上述定义是十分宽泛的，并不是所有满足定义的函数都是具有价值的小波，若不加限制条件其逆变换也不一定存在，所以小波必须施加以限制条件，称为容许条件。 定义 2.2 小波容许条件。假设为小波母函数，则其傅里叶变换应满足 上述条件被称为小波函数容许条件。 定义 2.3 连续小波变换定义。假设 ，是连续小波母函数，为连续小波基函数，则小波变换可定义为: 离散小波定义即对尺度因子a和平移因子b分别离散化。 3多分辨率分析和Mallat算法 3.1 多分辨率分析 多分辨率分析首先并不是应用于小波分析上，它曾成功应用于地球物理学，S.Mallat和 Y.Meyer 在前人大量工作的基础上于1986年把多分辨率分析应用到小波上，统一了小波基的构造，有人说只要掌握多尺度就可以自己构造自己的小波。分析空间中的多分辨率分析是指中满足下列条件的一个序列。 单调性：对任意,有。 逼近性：，。 伸缩性：，伸缩性体现了尺度的变换、逼近正交小波函数的变化和空间的变化具有一致性。 平移不变性：对任意，有。 Riesz基存在性：存在，使得构成的Riesz基。 3.2 Mallat算法 分解算法：假定信号，即 其中系数{，}为已知。如果我们希望将它分解为在和空间的两个分量之和，即 也就是由已知序列{}分别求出级的近似序列{}和级细节序列{}。 重构算法：当序列{}和{}为已知时，重构信号在空间中的分量就是利用这两个序列来计算序列{}。由 四：仿真程序与仿真结果 4.1 一维信号的Mallat算法仿真程序与仿真结果 程序： % 一维信号的分解与重构 load noisbloc x=noisbloc; l=length(x); subplot(3,1,1);plot(x);title(原始信号) %信号的一级分解与重构 [swa,swd]=swt(x,1,db5); subplot(3,1,2);plot(swa);title(低频); subplot(3,1,3);plot(swd);title(高频); A0=iswt(swa,swd,db5); err=norm(x-A0); figure(2); subplot(3,1,1);plot(x);title(原始信号); subplot(3,1,2);plot(A0);title(重构信号); subplot(3,1,3);plot(x-A0);title(误差信号); %信号的单支重构 nulcfs=zeros(size(swa)); A1=iswt(swa,nulcfs,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); for i=1:3 subplot(3,2,kp+1),plot(swa(i,:)); title([低频系数:level,num2str