小波变换PPT

多媒体技术教程第7章小波与小波变换林福宗清华大学计算机科学与技术系[e-mail]2008年9月2020年10月30日第7章小波与小波变换2/46第7章小波与小波变换目录7.1小波介绍7.1.1小波简史7.1.2小波概念7.1.3小波分析7.1.4小波定义7.2哈尔函数7.2.1哈尔基函数7.2.2哈尔小波函数7.2.3函数的规范化7.2.4哈尔基的结构7.3哈尔小波变换7.4规范化算法7.5二维哈尔小波变换7.5.1二维小波变换举例7.5.2二维小波变换方法2020年10月30日第7章小波与小波变换3/467.1小波介绍小波(wavelet)是什么在有限时间范围内变化且其平均值为零的数学函数具有有限的持续时间和突变的频率和振幅在有限的时间范围内，它的平均值等于零2020年10月30日第7章小波与小波变换4/467.1小波介绍(续1)部分小波许多数缩放函数和小波函数以开发者的名字命名，例如，Moret小波函数是Grossmann和Morlet在1984年开发的db6缩放函数和db6小波函数是Daubechies开发的图7-1正弦波与小波——部分小波2020年10月30日第7章小波与小波变换5/467.1小波介绍(续2)1807:JosephFourier傅立叶理论指出，一个信号可表示成一系列正弦和余弦函数之和，叫做傅立叶展开式小波简史小波变换(wavelettransform)是什么老课题：函数的表示方法新方法：Fourier－Haar－wavelettransform2020年10月30日第7章小波与小波变换6/467.1小波介绍(续3)wherecossinjtetjt()()()()ejtjtFftedtftF12只有频率分辨率而没有时间分辨率可确定信号中包含哪些频率的信号，但不能确定具有这些频率的信号出现在什么时候2020年10月30日第7章小波与小波变换7/467.1小波介绍(续4)1909:AlfredHaarAlfredHaar对在函数空间中寻找一个与傅立叶类似的基非常感兴趣。1909年他发现并使用了小波，后来被命名为哈尔小波(Haarwavelets)2020年10月30日第7章小波与小波变换8/467.1小波介绍(续5)1945:Gabor开发了STFT(shorttimeFouriertransform)(,)()where:()signal()=windo(wing)functionjtgSTFTstedtstgttSTFT的时间-频率关系图2020年10月30日第7章小波与小波变换9/467.1小波介绍(续6)1980：Morlet20世纪70年代，在法国石油公司工作的年轻地球物理学家JeanMorlet提出小波变换(wavelettransform，WT)的概念。20世纪80年代,开发了连续小波变换(continuouswavelettransform，CWT)1986：Y.Meyer法国科学家Y.Meyer与其同事创造性地构造出具有一定衰减性的光滑函数，用于分析函数用缩放(dilations)与平移(translations)均为2j(j≥0的整数)的倍数构造了L2(R)空间的规范正交基，使小波分析得到发展2020年10月30日第7章小波与小波变换10/467.1小波介绍(续7)1988：Mallat算法法国科学家StephaneMallat提出多分辨率概念，从空间上形象说明小波的多分辨率的特性，并提出了正交小波的构造方法和快速算法，称为Mallat算法[1]该算法统一了在此之前构造正交小波基的所有方法，其地位相当于快速傅立叶变换在经典傅立叶分析中的地位2020年10月30日第7章小波与小波变换11/467.1小波介绍(续8)小波理论与工程应用InridDaubechies于1988年最先揭示了小波变换和滤波器组(filterbanks)之间的内在关系[2]，使离散小波分析变成为现实RonaldCoifman和VictorWickerhauser等著名科学家在把小波理论引入到工程应用方面做出了极其重要贡献在信号处理中，自从StephaneMallat和InridDaubechies发现滤波器组与小波基函数有密切关系之后，小波分析在信号(如声音和图像)处理中得到极其广泛的应用2020年10月30日第7章小波与小波变换12/467.1小波介绍——小波分析小波分析/小波变换变换目的是获得时间和频率域之间的相互关系小波变换对一个函数在空间和时间上进行局部化的一种数学变换通过平移母小波(motherwavelet)获得信号的时间信息通过缩放母小波的宽度(或称尺度)获得信号的频率特性对母小波的平移和缩放操作是为计算小波的系数，这些系数代表局部信号和小波之间的相互关系对比傅立叶变换提供了频率域的信息，但丢失了时间域的局部化信息小波分析中常用的三个基本概念连续小波变换离散小波变换小波重构2020年10月30日第7章小波与小波变换13/467.1小波介绍——小波分析(续1)连续小波变换(continuouswavelettransform，CWT)傅立叶分析用一系列不同频率的正弦波表示一个信号一系列不同频率的正弦波是傅立叶变换的基函数小波分析用母小波通过移位和缩放后得到的一系列小波表示一个信号一系列小波可用作表示一些函数的基函数凡能用傅立叶分析的函数都可用小波分析小波变换可理解为用经过缩放和平移的一系列函数代替傅立叶变换用的正弦波用不规则的小波分析变化激烈的信号比用平滑的正弦波更有效，或者说对信号的基本特性描述得更好2020年10月30日第7章小波与小波变换14/467.1小波介绍——小波分析(续2)CWT的变换过程示例，见图7-3，可分如下5步1.小波ψ(t)和原始信号f(t)的开始部分进行比较2.计算系数C——该部分信号与小波的近似程度；C值越高表示信号与小波相似程度越高3.小波右移k得到的小波函数为ψ(t-k)，然后重复步骤1和2，……直到信号结束4.扩展小波，如扩展一倍，得到的小波函数为ψ(t/2)5.重复步骤1~4图7-3连续小波变换的过程2020年10月30日第7章小波与小波变换15/467.1小波介绍——小波分析(续3)连续小波变换用下式表示(,)()(,,)Cscalepositionftscalepositiontdt该式含义：小波变换是信号f(t)与被缩放和平移的小波函数Ψ之积在信号存在的整个期间里求和CWT变换的结果是许多小波系数C，这些系数是缩放因子(scale)和位置(position)的函数离散小波变换(discretewavelettransform，DWT)用小波的基函数(basisfunctions)表示一个函数的方法小波的基函数序列或称子小波(babywavelets)函数是由单个小波或称为母小波函数通过缩放和平移得到的缩放因子和平移参数都选择2j(j0的整数)的倍数，这种变换称为双尺度小波变换(dyadicwavelettransform)2020年10月30日第7章小波与小波变换16/467.1小波介绍——小波分析(续4)图7-5离散小波变换分析图DWT得到的小波系数、缩放因子和时间关系，见图7-5图(a)是20世纪40年代使用Gabor开发的短时傅立叶变换(shorttimeFouriertransform，STFT)得到的图(b)是20世纪80年代使用Morlet开发的小波变换得到的2020年10月30日第7章小波与小波变换17/467.1小波介绍——小波分析(续5)执行DWT的有效方法用Mallat在1988年开发的滤波器，称为Mallat算法[1]DWT的概念见图7-6。S表示原始的输入信号；通过两个互补的滤波器产生A和D两个信号图7-6双通道滤波过程A表示信号的近似值(approximations)，大的缩放因子产生的系数，表示信号的低频分量D表示信号的细节值(detail)，小的缩放因子产生的系数，表示信号的高频分量2020年10月30日第7章小波与小波变换18/467.1小波介绍——小波分析(续6)小波分解树与小波包分解树由低通滤波器和高通滤波器组成的树原始信号通过一对滤波器进行的分解叫做一级分解。信号的分解过程可以迭代，即可进行多级分解。小波分解树(waveletdecompositiontree)用下述方法分解形成的树：对信号的高频分量不再继续分解，而对低频分量连续进行分解，得到许多分辨率较低的低频分量，见图7-7小波包分解树(waveletpacketdecompositiontree)用下述方法分解形成的树：不仅对信号的低频分量连续进行分解，而且对高频分量也进行连续分解，这样不仅可得到许多分辨率较低的低频分量，而且也可得到许多分辨率较低的高频分量，见图7-82020年10月30日第7章小波与小波变换19/467.1小波介绍——小波分析(续7)图7-7小波分解树2020年10月30日第7章小波与小波变换20/467.1小波介绍——小波分析(续8)图7-8三级小波包分解树1332SAAADDADDD2020年10月30日第7章小波与小波变换21/467.1小波介绍——小波分析(续9)图7-9降采样过程注意：在使用滤波器对真实的数字信号进行变换时，得到的数据将是原始数据的两倍例如，如果原始信号的数据样本为1000个，通过滤波之后每一个通道的数据均为1000个，总共为2000个。于是,根据尼奎斯特(Nyquist)采样定理就提出了采用降采样(downsampling)的方法，即在每个通道中每两个样本数据中取一个，得到的离散小波变换的系数(coefficient)分别用cD和cA表示，见图7-92020年10月30日第7章小波与小波变换22/467.1小波介绍——小波分析(续10)小波重构重构概念把分解的系数还原成原始信号的过程叫做小波重构(waveletreconstruction)或合成(synthesis)，数学上叫做逆离散小波变换(inversediscretewavelettransform，IDWT)两个过程在使用滤波器做小波变换时包含滤波和降采样(downsampling)两个过程，在小波重构时也包含升采样(upsampling)和滤波两个过程，见图7-10升采样是在两个样本数据之间插入“0”，目的是把信号的分量加长，其过程见图7-112020年10月30日第7章小波与小波变换23/467.1小波介绍——小波分析(续11)图7-10小波重构方法图7-11升采样的方法2020年10月30日第7章小波与小波变换24/467.1小波介绍——小波分析(续12)重构滤波器滤波器关系到能否重构出满意的原始信号。在信号的分解期间，降采样会引进畸变，这种畸变叫做混叠(aliasing)。这就需要在分解和重构阶段精心选择关系紧密但不一定一致的滤波器才有可能取消这种混叠低通分解滤波器(L)和高通分解滤波器(H)以及重构滤波器(L'和H')构成一个系统，这个系统叫做正交镜像滤波器(quadraturemirrorfilters，QMF)系统，如图7-12所示图7-12正交镜像滤波器系统2020年10月30日第7章小波与小波变换25/467.1小波介绍——小波分析(续13)小波变换演示网址：~maziarz/Wavelets/要安装下面的插件MacromediaShockwaveplug-in该演示有声音解说2020年10月30日第7章小波与小波变换26/467.2哈尔函数7.2.1哈尔基函数基函数是一组线性无关的函数，可以用来构造任意给定的信号，如用基函数的加权和表示哈尔基函数(Haarb