当前位置:文档之家› 小波变换的理解

小波变换的理解

小波变换的理解
小波变换的理解

由于小波变换的知识涵盖了调和分析,实变函数论,泛函分析及矩阵论,所以没有一定的数学基础很难学好小波变换.但是对于我们工科学生来说,重要的是能利用这门知识来分析所遇到的问题.所以个人认为并不需要去详细学习调和分析,实变函数论,泛函分析及矩阵论等数学知识.最重要是的理解小波变换的思想!从这个意义上说付立叶变换这一关必需得过!因为小波变换的基础知识在付立叶变换中均有提及,我觉得这也就是很多小波变换的书都将付立叶分析作为其重要内容的原因.所以我认为学习小波应从<数字信号处理>中的付立叶分析开始.当然也可从<信号与系统>这本书开始.然后再看杨福生老师的小波变换书.个人觉得他的书最能为工科学生所接受.

2信号的分解

付立叶级数将周期信号分解为了一个个倍频分量的叠加,基函数是正交的,也就是通常所说的标准正交基.通过分解我们就能将特定的频率成分提取出来而实现特定的各种需要,如滤波,消噪等.付立叶变换则将倍频谱转换为了连续谱,其意义差不多.小波变换也是一种信号分解思想:只不过它是将信号分解为一个个频带信号的叠加.其中的低频部分作为信号的近似,高频部分作为信号的细节.所谓的细节部分就是一组组小波分量的叠加,也就是常说的小波级数.

3小波变换的时频分析思想

付立叶变换将信号从时域变换到了频域,从整体上看待信号所包含的频率成分.对于某个局部时间点或时间段上信号的频谱分析就无能为力了,对于我们从事信号的奇异性检测的人来说,付立叶变换就失去了意义(包括加窗付立叶变换).因为我们要找的是信号的奇异点(时域方面)和奇异点处所包含的频带(频域方面)也就是说需要一种时频分析方法.当然能有纯时域的分析方法更好!(据说数学形态学能达到这种效果).小波变换之所以可以检测信号的奇异点,正在于它的"小".因为用小的波去近似奇异信号要比正弦波要好的多.

4小波变换的实质

小波变换的公式有内积形式和卷积形式,两种形式的实质都是一样的.它要求的就是一个个小波分量的系数也就是"权".其直观意义就是首先用一个时窗最窄,频窗最宽的小波作为尺子去一步步地"量"信号,也就是去比较信号与小波的相似程度.信号局部与小波越相似,则小波变换的值越大,否则越小!当一步比较完成后,再将尺子拉长一倍,又去一步步地比

较,从而得出一组组数据.如此这般循环,最后得出的就是信号的小波分解(小波级数).当然这只是一种粗略的解释.

5连续小波变换,二进小波变换与离散小波变换的关系

当尺度及位移均作连续变化时,可以理解必将产生一大堆数据,作实际应用时并不需要这么多的数据,因此就产生了离散的思想.将尺度作二进离散就得到二进小波变换,同时也将信号的频带作了二进离散.当觉得二进离散数据量仍显大时,同时将位移也作离散就得到了离散小波变换!

6 MALLAT算法的意义

想必大家都注意到,小波变换是以内积或卷积的形式实现的,这给数值计算带来了不利之处,因为用计算机作数值积分其计算量大.MALLAT算法则解决了这一问题,它不涉及小波的具体形式,只是对系数进行操作!其计算也就是用高通及低通滤波系数与小波系数作卷积.因为作信号处理时,我们往往并不关心小皮的具体形式,更为关心小波系数.需提出的是该算法仅适用于正交小波如果小波不是正交的(如B样条小波)则算法失效!

7小波变换的模极大值及其意义

对于我们搞信号奇异性检测的人来说,小波变换最重要的应用就是用模极大值定值奇异点.我觉得模极大值可以从两个方面去理解:第一,从直观角度,上文已说明小波变换的实质就是一种度量波形相似程度的方法.信号与小波越相似,则小波系数越大.这也就可理解为出现了小波变换的模极大值.因为当信号出现奇异点时,或是间断点,或是一阶导数不连续点,其在各个尺度下都将必然出现大的小波系数.从而可以定位奇异点!第二个方面从小波的取法来看,当小波取为光滑函数一阶导数或二阶导数时,从公式可以推导出小波变换将出现模极大值点或是过零点.也就是很多书上说的模极大值检测和零交叉检测.这些可以查书看!

我只谈谈连续小波变换,对于离散的也有同样的argument。小波函数的dilation和translation是这样一个形式:1/\sqrt{|s|}\psi((x-u)/s),s是scale,u是该小波atom 的center。

由于根据定义,小波的积分是0,也就是说小波函数的傅立叶变换在零点为零。再有于小波函数的傅立叶变换一般是连续的(比如如果小波是属于L_1的),这样在0的一个小临域

里面,小波的傅立叶变换很接近零,这也就是说小波函数的傅立叶变换可以看成某个高通滤波器的transfer function,这样小波变换W(f)实际是在measure该函数f在u点附近的variation。从这个角度看的话,如果小波的宽度很大(对应尺度s很大),该函数在该小波的窗口下的variation就很大;如果小波的宽度小(对应尺度s小),则函数在该小波的窗口下的variation就相对比较小(除非信号是fractal,呵呵)。

小波情结

到了小波版很久,总觉得应该写些什么。这篇文章也就献给那些所有正在研究或即将研究小波的同学、老师和科研人员们。这是篇与技术无关的文章,撰写的是我对小波的感受。从我开始接触小波,研究小波,到迷恋小波的真实记录。因此,我把它起名叫小波情结。刚开始,接触小波的时候在研一。关于老师布置的从频域构造一种小波的作业开始。后来我才知道,这种小波本质上就是meyer小波。当时,就一个字,嫩。实际上就是对小波毫无所知。脑子里就是一叠的公式。正交条件,容许条件等一大堆,与概念理解相差甚远的东东。但,还是乐此不疲的编程。总想看看,我亲手缔造的小波长的是什么样,也有些略带孩子气的,想把它作为桌面和自己的酷酷头像之类的欲望。于是,十一的头三天,我基本上闭门造车。当时,我用的是matlab,也是我最后得到哭笑不得结果的直接的助手与帮凶。因为构造的过程的起始,我就把函数离散化了。紧接着就是平移,对乘,积分,抽取,插值,dsp里的一套trick 把我搞得叫苦不迭。程序也累计到了1000行左右。当时,最可恨的就是对点,由于dsp下标的1,2,3离散化,所以我也就用手指开始傻傻的算。连续几天晚上鏖战,终于在3号的晚上。通过IFFT后,美妙的波形出来了。注意,美妙和丑陋只有一步之遥。这是我的对小波的第一课体会。当我一看屏幕,疯了,彻底疯了。一个DELTA函数类似的波形,就在我眼前。心想:忙乎了三天,整了个DELTA函数出来。这难道就是回报吗?别急,小波是紧支撑的啊。概念上对头,一定是取点的问题。我便拿起MATLAB自带的照妖镜(放大镜)一看,呵呵,一个差强人意的波形就在我眼前了。我当时大喊一声,爽,那时已经凌晨2:00。

第一次的经历,对我来说收获很丰。然后,第二次挑战,则是彻底改变我对小波是个深不可测的家伙的看法。这次作业,就是用刚才构造的小波,做消噪。我这次,又一次的想起,爱情格言:我心灵的古堡经不起你轻轻的一击。女生问:结果呢?回答:碎了。一个如此,不

精确的波形,怎么能消噪呢?而且,当时老师要用连续小波的方法。也就是内积求和的方法。我和同学,首先合作,用mathmatic 做了个好一点的波形。因为,除最后一步,反傅里叶变换外,其他都是解析的。然后,一个困扰我许久的问题产生了。一个函数可以由无穷多个小波的膨胀和伸缩叠加起来。那么,我把函数从-inf 到+inf 积分,假设函数有直流分量,所以积分不为零。但是小波,积分却为零。这不是矛盾吗?后来,也就是研二我才知道,有些时候积分后不可以交换。还有,其实有限的小波逼近,必须加上尺度函数才可以。但当时,我们只是采用了把小波的支撑取宽的办法解决了此问题。但,我由于不太喜欢这种方法的冗长和费时,所以想令辟蹊径。于是,mallet一个令我崇拜的算法,终于在我阅读超星的时候,跳在了我生命里。首先,便是看冗长的证明,勉强理解了。当看到滤波器组的解释后,我开始豁然开朗。这是我熟悉的dsp概念。因此,我花了一晚上,把这个算法彻底搞懂了。但概念的理解和程序的成功编制,还是有一小步,就是这一步,使无数英雄竟折腰。我的幸运之神便是MATLAB里的DEMO。那个里面,有一个详细的算法解释。并且从哪里我知道了些怪怪的函数。WKEEP(),DYADDOWN(),DYADUP()等等。而且,又一个问题,理论和实际差别产生了。

这个问题甚至现在,还困扰着很多的小波工作者。一个长度为100的信号,分解后理论上高频50,低频50。但用卷积算法,假设滤波器长度为10。因此总长度109,做抽取后长度55。多了5。这怎么办呢。我去问了很多老师,回答都一样。就是MATLAB里用的函数WKEEP()。把两头丢掉。当时我勉强接受了这个结果。但始终有个概念,小波变换就是正交变换,它和傅里也变换一样,一定

可以写成正交阵的形式。

第二次作业的完成,我的小波课结束了。但我的小波情结还在继续。关于,正交阵的猜想还在困扰着我。一本电磁场和小波结合的外文书籍,帮助了我。圆周卷积的概念,历历在目。是呀,卷积对着傅里叶变换,而圆周卷积对应着离散傅里叶变换。这就是连续与离散的区别和联系啊。于是我用db小波,构造了一个完全正交的矩阵。当我把这个矩阵和它的转置相乘的时候,单位阵出来了。那天,我高兴得流泪。最终,我把圆周卷积用快速傅里叶变换实现出来。今后的日子,我便觉得,思维的水再也关不住了。

步步为营,我实现了db小波的时域构造,采用矩阵特征向量法和casade理论两种解法,我都成功了。慢慢的我开始醉心于消失矩,开始懂得框架,开始懂得双正交。然后就是,PR条件,二代小波,小波插值,因子化,等等。于是,我也在研学一边和大家交流,一边阅读大量书籍和文献,而且实现里面的每一个例子和思想。

当我们还在觉得自己懂点小波的时候,美国人已把它用于指纹压缩,产生了巨大的经济和社会效益;当我们,还在对二代不屑一顾的时候,一个叫JPEG2000标准的东东,彻底给我们上了一课。当我们,还在国家著名期刊上,打着错误的提升公式的时候,当我们,还在为些不值一提的程序保密的时候,一个叫各相异性小波的东东又开始蠢蠢欲动。看看那些大师们吧,看看他们的态度,再看看我们,我们努力的够吗。你说看不懂文献,我就要问你,你看了一遍,十遍,还是一百遍呢?如果说你认为是高手,你是否写了超过10万行以上小波的代码,看了10本以上的书,100篇的文献,实现里面所有的例子和思想了呢。我们差得很远。

但是我们服气吗,我们认输了吗,我们不再努力了吗。什么时候有中国的JPEG2008呢,什么时候我们能毫无保留的进行坦诚的交流和无私的分享呢,什么时候我们把学术的铜臭拨掉,把做小波看成一次和上苍对话的机会,和真理的交锋呢。我始终在问自己这些问题。为关于学术的单纯的问题。我找到了答案。在研学上。因此,我毫无保留的帮助大家,同时也在修正自己。我开始变得勤奋,开始编每个需要的程序,而基本上不用MATLAB提供的任何函数,除非是概念性验证。

我和所有那些从事这小波事业的人们一样,为实现这上面单纯而坚定的疑问而不停奋斗。我和你们一样,是一个奔跑着,一个向着小波的科学和真理殿堂不辞辛劳的奔跑者。即使路上满是荆棘,即使我们会暂时的迷失方向,但我相信小波这朵最美丽的奇葩,会以它最美丽的身姿,最沁人的芬芳,指引着我们。让我们结伴而行吧!我是引路者,也是跟随者。和你们一样,怀揣梦想,一起努力。虽然汗流浃背,虽然荆棘满身,但成功终会来的。我悯悯驾信,我正在奔跑,像春天里的孩子。

小波变换与傅里叶变换的对比异同

小波变换与傅里叶变换的对比、异同 一、基的概念 两者都是基,信号都可以分成无穷多个他们的和(叠加)。而展开系数就是基与信号之间的内积,更通俗的说是投影。展开系数大的,说明信号和基是足够相似的。这也就是相似性检测的思想。但我们必须明确的是,傅里叶是0-2pi 标准正交基,而小波是-inf到inf之间的基。因此,小波在实轴上是紧的。而傅里叶的基(正弦或余弦),与此相反。而小波能不能成为Reisz基,或标准稳定的正交基,还有其它的限制条件。此外,两者相似的还有就是PARSEVAL定理。(时频能量守恒)。 二、离散化的处理 傅里叶变换,是一种数学的精妙描述。但计算机实现,却是一步步把时域和频域离散化而来的。第一步,时域离散化,我们得到离散时间傅里叶变换(DTFT),频谱被周期化;第二步,再将频域离散化,我们得到离散周期傅里叶级数(DFS),时域进一步被周期化。第三步,考虑到周期离散化的时域和频域,我们只取一个周期研究,也就是众所周知的离散傅里叶变换(DFT)。这里说一句,DFT是没有物理意义的,它只是我们研究的需要。借此,计算机的处理才成为可能。所有满足容许性条件(从-INF到+INF积分为零)的函数,都可以成为小波。小波作为尺度膨胀和空间移位的一组函数也就诞生了。但连续取值的尺度因子和平移因子,在时域计算量和频域的混叠来说,都是极为不便的。用更为专业的俗语,叫再生核。也就是,对于任何一个尺度a和平移因子b的小波,和原信号内积,所得到的小波系数,都可以表示成,在a,b附近生成的小波,投影后小波系数的线性组合。这就叫冗余性。这时的连续小波是与正交基毫无关系的东西,它顶多也只能作为一种积分变换或基。但它的显微镜特点和相似性检测能力,已经显现出来了。为了进一步更好的将连续小波变换离散化,以下步骤是一种有效方法。第一步,尺度离散化。一般只将a二进离散化,此时b是任意的。这样小波被称为二进小波。第二步,离散b。怎么离散化呢?b取多少才合适呢?于是,叫小波采样定理的东西,就这样诞生了。也就是小波平移的最小距离(采样间隔),应该大于二倍小波基的最高频率(好像类似,记不清了)。所以b取尺度的整数倍就行了。也就是越胖的小波,对应频谱越窄,平移量应该越大,采样间隔越大。当然,第一二两步的频域理解,即在满足频域窗口中心是3倍的频域窗口半径的前提下,频域就在统计上是完美二分的。(但很多小波满足不了这个条件,而且频域窗口能量不?,所以只是近似二分的).这时的小波变换,称为离散二进小波变换.第三步,引入稳定性条件.也就是经过变换后信号能量和原信号能量有什么不等式关系.满足稳定性条件?后,也就是一个小波框架产生了可能.他是数值稳定性的保证.一个稍弱的稳定条件???,就是?

小波变换的基本原理

10.2小波变换的基本原理 地质雷达的电磁波信号和地震波信号都是非平稳随机时变信号,长期以来,因非平稳信号处理的理论不健全,只好将其作为平稳信号来处理,其处理结果当然不满意。近年来,随着科学技术的发展和进步,国内外学术界已将注意力转向非平稳随机信号分析与处理的研究上,其中非平稳随机信号的时频表示法是研究热点之一。在这一研究中,戈勃展开、小波变换、维格纳分布与广义双线性时频分布等理论发展起来,这些方法既可以处理平稳信号过程,也可以处理非平稳随机时变信号。 小波变换是上世纪80年代中后期逐渐发展起来的一种数学分析方法。1984年法国科学家J.M OLET在分析地震波的局部特性时首先使用了小波这一术语,并用小波变换对地震信号进行处理。小波术语的含义是指一组衰减震动的波形,其振幅正负相间变化,平均值为零,是具有一定的带宽和中心频率波组。小波变换是用伸缩和平移小波形成的小波基来分解(变换)或重构(反变换)时变信号的过程。不同的小波具有不同带宽和中心频率,同一小波集中的带宽与中心频率的比是不变的,小波变换是一系列的带通滤波响应。它的数学过程与傅立叶分析是相似的,只是在傅立叶分析中的基函数是单频的调和函数,而小波分析中的基函数是小波,是一可变带宽内调和函数的组合。 小波变换在时域和频域都具有很好的局部化性质,较好地解决了时域和频域分辨率的矛盾,对于信号的低频成分采用宽时窗,对高频成分采用窄时窗。因而,小波分析特别适合处理非平稳时变信号,在语音分析和图象处理中有广泛的应用,在地震、雷达资料处理中将有良好的应用前景。 下边就小波分析的基本原理、主要作用及在雷达资料处理中的应用三方面作以介绍。 10.2.1小波分析的基本原理 小波函数的数学表达

小波变换的几个典型应用

第六章小波变换的几个典型应用 6.1 小波变换与信号处理 小波变换作为信号处理的一种手段,逐渐被越来越多领域的理论工作者和工程技术人员所重视和应用,并在许多应用中取得了显著的效果。同传统的处理方法相比,小波变换取得了质的飞跃,在信号处理方面具有更大的优势。比如小波变换可以用于电力负载信号的分析与处理,用于语音信号的分析、变换和综合,还可以检测噪声中的未知瞬态信号。本部分将举例说明。 6.1.1 小波变换在信号分析中的应用 [例6-1] 以含躁的三角波与正弦波的组合信号为例具体说如何利用小波分析来分析信号。已知信号的表达式为 应用db5小波对该信号进行7层分解。xiaobo0601.m 图6-1含躁的三角波与正弦波混合信号波形 分析: (1)在图6-2中,逼近信号a7是一个三角波。 (2)在图6-3中细节信号d1和d2是与噪声相关的,而d3(特别是d4)与正弦信号相关。 图6-2 小波分解后各层逼近信号 图6-3 小波分解后各层细节信号 6.1.2 小波变换在信号降躁和压缩中的应用 一、信号降躁 1.工程中,有用信号一般是一些比较平稳的信号,噪声通常表现为高频信号。2.消躁处理的方法:首先对信号进行小波分解,由于噪声信号多包含在具有较高频率的细节中,我们可以利用门限、阈值等形式对分解所得的小波系数进行处理,然后对信号进行小波重构即可达到对信号的消躁目的。 小波分析进行消躁处理的3种方法: (1)默认阈值消躁处理。该方法利用ddencmp生成信号的默认阈值,然后利用wdencmp函数进行消躁处理。 (2)给定阈值消躁处理。在实际的消躁处理过程中,阈值往往可通过经验公式获得,且这种阈值比默认阈值的可信度高。在进行阈值量化处理时可利用函数wthresh。 (3)强制消躁处理。该方法时将小波分解结构中的高频系数全部置为0,即滤掉所有高频部分,然后对信号进行小波重构。方法简单,消躁后信号比较平滑,但易丢失信号中的有用成分。 小波阈值去噪方法是目前应用最为广泛的小波去噪方法之一。 3.信号降噪的准则: 1.光滑性:在大部分情况下,降噪后的信号应该至少和原信号具有同等的光滑性。

小波变换的原理及matlab仿真程序讲解学习

小波变换的原理及m a t l a b仿真程序

基于小波变换的信号降噪研究 2 小波分析基本理论 设Ψ(t)∈L 2( R) ( L 2( R) 表示平方可积的实数空间,即能量有限的信号空间) , 其傅立叶变换为Ψ(t)。当Ψ(t)满足条件[4,7]: 2 () R t dw w C ψψ =<∞? (1) 时,我们称Ψ(t)为一个基本小波或母小波,将母小波函数Ψ(t)经伸缩和平移后,就可以得到一个小波序列: ,()( )a b t b t a ψ -= ,,0a b R a ∈≠ (2) 其中a 为伸缩因子,b 为平移因子。 对于任意的函数f(t)∈L 2( R)的连续小波变换为: ,(,),()( )f a b R t b W a b f f t dt a ψψ-=<>= ? (3) 其逆变换为: 211()(,)()f R R t b f t W a b dadb C a a ψ ψ+-= ?? (4) 小波变换的时频窗是可以由伸缩因子a 和平移因子b 来调节的,平移因子b,可以改变窗口在相平面时间轴上的位置,而伸缩因子b 的大小不仅能影响窗口在频率轴上的位置,还能改变窗口的形状。小波变换对不同的频率在时域上的取样步长是可调节的,在低频时,小波变换的时间分辨率较低,频率分辨率较高:在高频时,小波变换的时间分辨率较高,而频率分辨率较低。使用小波变换处理信号时,首先选取适当的小波函数对信号进行分解,其次对分解出的参

数进行阈值处理,选取合适的阈值进行分析,最后利用处理后的参数进行逆小波变换,对信号进行重构。 3 小波降噪的原理和方法 3.1 小波降噪原理 从信号学的角度看 ,小波去噪是一个信号滤波的问题。尽管在很大程度上小波去噪可以看成是低通滤波 ,但由于在去噪后 ,还能成功地保留信号特征 ,所以在这一点上又优于传统的低通滤波器。由此可见 ,小波去噪实际上是特征提取和低通滤波的综合 ,其流程框图如图所示[6]: 小波分析的重要应用之一就是用于信号消噪 ,一个含噪的一维信号模型可表示为如下形式: (k)()()S f k e k ε=+* k=0.1…….n-1 其中 ,f( k)为有用信号,s(k)为含噪声信号,e(k)为噪声,ε为噪声系数的标准偏差。 假设e(k)为高斯白噪声,通常情况下有用信号表现为低频部分或是一些比较平稳的信号,而噪声信号则表现为高频的信号,下面对 s(k)信号进行如图结构的小波分解,则噪声部分通常包含在Cd1、Cd2、Cd3中,只要对 Cd1,Cd2,Cd3作相应的小波系数处理,然后对信号进行重构即可以达到消噪的目的。

第9章小波变换基础

第9章 小波变换基础 9.1 小波变换的定义 给定一个基本函数)(t ψ,令 )(1)(,a b t a t b a -= ψψ (9.1.1) 式中b a ,均为常数,且0>a 。显然,)(,t b a ψ是基本函数)(t ψ先作移位再作伸缩以后得到的。若b a ,不断地变化,我们可得到一族函数)(,t b a ψ。给定平方可积的信号)(t x ,即 )()(2R L t x ∈,则)(t x 的小波变换(Wavelet Transform ,WT )定义为 dt a b t t x a b a WT x )()(1),(-= ? *ψ ??==? * )(),()()(,,t t x dt t t x b a b a ψψ (9.1.2) 式中b a ,和t 均是连续变量,因此该式又称为连续小波变换(CWT )。如无特别说明,式中及以后各式中的积分都是从∞-到∞+。信号)(t x 的小波变换),(b a WT x 是a 和b 的函数, b 是时移,a 是尺度因子。)(t ψ又称为基本小波,或母小波。)(,t b a ψ是母小波经移位和 伸缩所产生的一族函数,我们称之为小波基函数,或简称小波基。这样,(9.1.2)式的WT 又可解释为信号)(t x 和一族小波基的内积。 母小波可以是实函数,也可以是复函数。若)(t x 是实信号,)(t ψ也是实的,则 ),(b a WT x 也是实的,反之,),(b a WT x 为复函数。 在(9.1.1)式中,b 的作用是确定对)(t x 分析的时间位置,也即时间中心。尺度因子 a 的作用是把基本小波)(t ψ作伸缩。我们在1.1节中已指出,由)(t ψ变成)(a t ψ,当1 >a 时,若a 越大,则)(a t ψ的时域支撑范围(即时域宽度)较之)(t ψ变得越大,反之,当1

小波变换详解

基于小波变换的人脸识别 近年来,小波变换在科技界备受重视,不仅形成了一个新的数学分支,而且被广泛地应用于模式识别、信号处理、语音识别与合成、图像处理、计算机视觉等工程技术领域。小波变换具有良好的时频域局部化特性,且其可通过对高频成分采取逐步精细的时域取样步长,从而达到聚焦对象任意细节的目的,这一特性被称为小波变换的“变聚焦”特性,小波变换也因此被人们冠以“数学显微镜”的美誉。 具体到人脸识别方面,小波变换能够将人脸图像分解成具有不同分辨率、频率特征以及不同方向特性的一系列子带信号,从而更好地实现不同分辨率的人脸图像特征提取。 4.1 小波变换的研究背景 法国数学家傅立叶于1807年提出了著名的傅立叶变换,第一次引入“频率”的概念。傅立叶变换用信号的频谱特性来研究和表示信号的时频特性,通过将复杂的时间信号转换到频率域中,使很多在时域中模糊不清的问题,在频域中一目了然。在早期的信号处理领域,傅立叶变换具有重要的影响和地位。定义信号(t)f 为在(-∞,+∞)内绝对可积的一个连续函数,则(t)f 的傅立叶变换定义如下: ()()dt e t f F t j ωω-? ∞ -∞ += (4-1) 傅立叶变换的逆变换为: ()()ωωπ ωd e F t f t j ? +∞ ∞ -= 21 (4-2) 从上面两个式子可以看出,式(4-1)通过无限的时间量来实现对单个频率

的频谱计算,该式表明()F ω这一频域过程的任一频率的值都是由整个时间域上的量所决定的。可见,式(4-1)和(4-2)只是同一能量信号的两种不同表现形式。 尽管傅立叶变换可以关联信号的时频特征,从而分别从时域和频域对信号进行分析,但却无法将两者有效地结合起来,因此傅立叶变换在信号的局部化分析方面存在严重不足。但在许多实际应用中,如地震信号分析、核医学图像信号分析等,研究者们往往需要了解某个局部时段上出现了哪个频率,或是某个频率出现在哪个时段上,即信号的时频局部化特征,傅立叶变换对于此类分析无能为力。 因此需要一种如下的数学工具:可以将信号的时域和频域结合起来构成信号的时频谱,描述和分析其时频联合特征,这就是所谓的时频局部化分析方法,即时频分析法。1964年,Gabor 等人在傅立叶变换的基础上引入了一个时间局部化“窗函数”g(t),改进了傅立叶变换的不足,形成窗口化傅立叶变换,又称“Gabor 变换”。 定义“窗函数”(t)g 在有限的区间外恒等于零或很快地趋于零,用函数(t )g -τ乘以(t)f ,其效果等同于在t =τ附近打开一个窗口,即: ()()()dt e t g t f G t j f ωττω-+∞ ∞--=?, (4-3) 式(4-3)即为函数f(t)关于g(t)的Gabor 变换。由定义可知,信号(t)f 的Gabor 变换可以反映该信号在t =τ附近的频谱特性。其逆变换公式为: ()()()ττωτωπ ωd G t g e d t f f t j ,21 ? ?+∞ ∞ --- = (4-4) 可见()τω,f G 的确包含了信号(t)f 的全部信息,且Gabor 窗口位置可以随着 τ的变化而平移,符合信号时频局部化分析的要求。 虽然Gabor 变换一定程度上克服了傅立叶变换缺乏时频局部分析能力的不

小波变换-完美通俗解读

小波变换和motion信号处理(一) 这是《小波变换和motion信号处理》系列的第一篇,基础普及。第二篇我准备写深入小波的东西,第三篇讲解应用。 记得我还在大四的时候,在申请出国和保研中犹豫了好一阵,骨子里的保守最后让我选择了先保研。当然后来也退学了,不过这是后话。当时保研就要找老板,实验室,自己运气还不错,进了一个在本校很牛逼的实验室干活路。我们实验室主要是搞图像的,实力在全国也是很强的,进去后和师兄师姐聊,大家都在搞什么小波变换,H264之类的。当时的我心思都不在这方面,尽搞什么操作系统移植,ARM+FPGA 这些东西了。对小波变换的认识也就停留在神秘的“图像视频压缩算法之王”上面。 后来我才发现,在别的很广泛的领域中,小波也逐渐开始流行。比如话说很早以前,我们接触的信号频域处理基本都是傅立叶和拉普拉斯的天下。但这些年,小波在信号分析中的逐渐兴盛和普及。这让人不得不感到好奇,是什么特性让它在图象压缩,信号处理这些关键应用中更得到信赖呢?说实话,我还在国的时候,就开始好奇这个问题了,于是放狗搜,放毒搜,找遍了中文讲小波变换的科普文章,发现没几个讲得清楚的,当时好奇心没那么重,也不是搞这个研究的,懒得找英文大部头论文了,于是作罢。后来来了这边,有些项目要用信号处理,不得已接触到一些小波变换的东西,才开始硬着头皮看。看了一

些材料,听了一些课,才发现,还是那个老生常谈的论调:国外的技术资料和国真TNND不是一个档次的。同样的事情,别人说得很清楚,连我这种并不聪明的人也看得懂; 国的材料则绕来绕去讲得一塌糊涂,除了少数天才没几个人能在短时间掌握的。 牢骚就不继续发挥了。在这个系列文章里,我希望能简单介绍一下小波变换,它和傅立叶变换的比较,以及它在移动平台做motion detection的应用。如果不做特殊说明,均以离散小波为例子。考虑到我以前看中文资料的痛苦程度,我会尽量用简单,但是直观的方式去介绍。有些必要的公式是不能少的,但我尽量少用公式,多用图。另外,我不是一个好的翻译者,所以对于某些实在翻译不清楚的术语,我就会直接用英语。我并不claim我会把整个小波变换讲清楚,这是不可能的事,我只能尽力去围绕要点展开,比如小波变换相对傅立叶变换的好处,这些好处的原因是什么,小波变换的几个根本性质是什么,背后的推导是什么。我希望达到的目的就是一个小波变换的初学者在看完这个系列之后,就能用matlab或者别的工具对信号做小波变换的基本分析并且知道这个分析大概是怎么回事。 最后说明,我不是研究信号处理的专业人士,所以文中必有疏漏或者错误,如发现还请不吝赐教。 要讲小波变换,我们必须了解傅立叶变换。要了解傅立叶变换,我们先要弄清楚什么是”变换“。很多处理,不管是压缩也好,滤波也好,图形处理也好,本质都是变换。变换的是什么东西呢?是基,也就是

详解傅里叶变换与小波变换

详解傅里叶变换与小波变化 希望能简单介绍一下小波变换,它和傅立叶变换的比较,以及它在移动平台做motion detection的应用。如果不做特殊说明,均以离散小波为例子。考虑到我以前看中文资料的痛苦程度,我会尽量用简单,但是直观的方式去介绍。有些必要的公式是不能少的,但我尽量少用公式,多用图。另外,我不是一个好的翻译者,所以对于某些实在翻译不清楚的术语,我就会直接用英语。我并不claim我会把整个小波变换讲清楚,这是不可能的事,我只能尽力去围绕要点展开,比如小波变换相对傅立叶变换的好处,这些好处的原因是什么,小波变换的几个根本性质是什么,背后的推导是什么。我希望达到的目的就是一个小波变换的初学者在看完这个系列之后,就能用matlab或者别的工具对信号做小波变换的基本分析并且知道这个分析大概是怎么回事。 要讲小波变换,我们必须了解傅立叶变换。要了解傅立叶变换,我们先要弄清楚什么是”变换“。很多处理,不管是压缩也好,滤波也好,图形处理也好,本质都是变换。变换的是什么东西呢?是基,也就是basis。如果你暂时有些遗忘了basis的定义,那么简单说,在线性代

数里,basis是指空间里一系列线性独立的向量,而这个空间里的任何其他向量,都可以由这些个向量的线性组合来表示。那basis在变换里面啥用呢?比如说吧,傅立叶展开的本质,就是把一个空间中的信号用该空间的某个basis的线性组合表示出来,要这样表示的原因,是因为傅立叶变换的本质,是。小波变换自然也不例外的和basis有关了。再比如你用Photoshop去处理图像,里面的图像拉伸,反转,等等一系列操作,都是和basis的改变有关。 既然这些变换都是在搞基,那我们自然就容易想到,这个basis的选取非常重要,因为basis的特点决定了具体的计算过程。一个空间中可能有很多种形式的basis,什么样的basis比较好,很大程度上取决于这个basis服务于什么应用。比如如果我们希望选取有利于压缩的话,那么就希望这个basis能用其中很少的向量来最大程度地表示信号,这样即使把别的向量给砍了,信号也不会损失很多。而如果是图形处理中常见的线性变换,最省计算量的完美basis就是eigenvector basis了,因为此时变换矩阵T对它们的作用等同于对角矩阵(Tv_n= av_n,a是eigenvalue)。总的来说,抛开具体的应用不谈,所有的basis,我们都希望它们有一个共同的特点,那就是,容易计算,用最简单的方式呈现最多的信号特性。 好,现在我们对变换有了基本的认识,知道他们其实就是在搞基。当然,搞基也是分形式的,不同的变换,搞基的妙处各有不同。接下来先看看,傅立叶变换是在干嘛。

小波变换基本原理

第五章 小波变换基本原理 问题 ①小波变换如何实现时频分析?其频率轴刻度如何标定? —尺度 ②小波发展史 ③小波变换与短时傅里叶变换比较 a .适用领域不同 b.STFT 任意窗函数 WT (要容许性条件) ④小波相关概念,数值实现算法 多分辨率分析(哈尔小波为例) Daubechies 正交小波构造 MRA 的滤波器实现 ⑤小波的历史地位仍不如FT ,并不是万能的 5.1 连续小波变换 一.CWT 与时频分析 1.概念:? +∞ ∞ --ψ= dt a b t t S a b a CWT )( *)(1),( 2.小波变换与STFT 用于时频分析的区别 小波 构造? 1910 Harr 小波 80年代初兴起 Meyer —小波解析形式 80年代末 Mallat 多分辨率分析—WT 无须尺度和小波函数—滤波器组实现 90年代初 Daubechies 正交小波变换 90年代中后期 Sweblews 第二代小波变换

3.WT 与STFT 对比举例(Fig 5–6, Fig 5–7) 二.WT 几个注意的问题 1.WT 与)(t ψ选择有关 — 应用信号分析还是信号复原 2.母小波)(t ψ必须满足容许性条件 ∞<ψ=? ∞ +∞ -ψdw w w C 2 )( ①隐含要求 )(,0)0(t ψ=ψ即具有带通特性 ②利用ψC 可推出反变换表达式 ??+∞∞-+∞ ∞-ψ -ψ= dadb a b t b a CWT a C t S )(),(11 )(2 3.CWT 高度冗余(与CSTFT 相似) 4.二进小波变换(对平移量b 和尺度进行离散化) )2(2)()(1 )(2 ,22,,n t t a b t a t n b a m m n m b a m m -ψ=ψ?-ψ= ??==--ψ dt t t S n CWT d n m m m n m )(*)()2,2(,,?+∞ ∞ ---ψ=?= 5.小波变换具有时移不变性 ) ,()() ,()(00b b a C W T b t S b a C W T t S -?-? 6.用小波重构信号 ∑∑ ∑∑+∞-∞=+∞ -∞ =+∞-∞=+∞ -∞ =ψψ= m n m n n m n m n m n m t d t d t S )(?)(?)(,,,,正交小波 中心问题:如何构建对偶框架{} n m ,?ψ

连续小波变换的概念

连续小波变换的概念swt,cwt,dwt 1。连续小波的概念。就是把一个可以称作小波的函数(从负无穷到正无穷积分为零)在某个尺度下与待处理信号卷积。改变小波函数的尺度,也就改变了滤波器的带通范围,相应每一尺度下的小波系数也就反映了对应通带的信息。本质上,连续小波也就是一组可控制通带范围的多尺度滤波器。 2。连续小波是尺度可连续取值的小波,里面的a一般取整数,而不像二进小波a取2的整数幂。从连续小波到二进小波再到正交离散小波,其实就是a、b都连续,a不连续、b连续,a、b都不连续的过程。操作他们的快速算法也就是卷积(快速傅里叶),多孔(a trous),MALLAT。在MATLAB里,也就是CWT,SWT,DWT。SWT称平稳小波变换、二进小波变换、或者非抽取小波变换。3。从冗余性上:CWT>SWT>DWT,前面两个都冗余,后面的离散小波变换不冗余。 4。从应用上:CWT适合相似性检测、奇异性分析;SWT适合消噪,模极大值分析;DWT适合压缩。 5。操作。就是在某个尺度上得到小波的离散值和原信号卷积,再改变尺度重新得到小波的离散值和原信号卷积。每一个尺度得到一个行向量存储这个尺度下的小波系数,多个尺度就是一个矩阵,这个矩阵就是我们要显示的时间-尺度图。 6。显示。“不要认为工程很简单”。我的一个老师说过的话。小波系数的显示还是有技巧的。很多人画出的图形“一片乌黑”就是个例子。第一步,一般将所有尺度下的小波系数取模;第二步,将每个尺度下的小波系数范围作映射,映射到你指定MAP的范围,比如如果是GRAY,你就映射到0-255;第三步,用IMAGE命令画图;第四步,设置时间和尺度坐标。MATLAB是个很专业的软件,它把这些做的很好,但也就使我们懒惰和糊涂,我是个好奇心重的人就研究了下。里面有个巧妙的函数把我说的(1,2)两个步骤封装在了一起,就是WCODEMAT,有兴趣的同学可以看看。 希望大家深入研究小波。 这里,还有要说的是,小波目前理论的热点: 1。不可分的小波或者具有可分性质的方向性小波; 2。XLET: CONTOURLET, WEDGELET, SHEARLET, BANDELET, RIDGELET, CURVELET; PLATELET. 3。多分辨率分析+多尺度几何分析的结合,才真正是我们所需要的。比如小波域的WEDGELET等等。 最后,几点建议: 1。理论研究和实际应用不同,工程上很多问题小波并不是最好的,在做项目的时候大家要实际情况,实际对待。

小波变换与傅里叶变换的对比异同

小波变换与傅里叶变换 的对比异同 IMB standardization office【IMB 5AB- IMBK 08- IMB 2C】

小波变换与傅里叶变换的对比、异同 一、基的概念 两者都是基,信号都可以分成无穷多个他们的和(叠加)。而展开系数就是基与信号之间的内积,更通俗的说是投影。展开系数大的,说明信号和基是足够相似的。这也就是相似性检测的思想。但我们必须明确的是,傅里叶是0-2pi标准正交基,而小波是-inf到inf之间的基。因此,小波在实轴上是紧的。而傅里叶的基(正弦或余弦),与此相反。而小波能不能成为Reisz基,或标准稳定的正交基,还有其它的限制条件。此外,两者相似的还有就是PARSEVAL 定理。(时频能量守恒)。 二、离散化的处理 傅里叶变换,是一种数学的精妙描述。但计算机实现,却是一步步把时域和频域离散化而来的。第一步,时域离散化,我们得到离散时间傅里叶变换(DTFT),频谱被周期化;第二步,再将频域离散化,我们得到离散周期傅里叶级数(DFS),时域进一步被周期化。第三步,考虑到周期离散化的时域和频域,我们只取一个周期研究,也就是众所周知的离散傅里叶变换(DFT)。这里说一句,DFT是没有物理意义的,它只是我们研究的需要。借此,计算机的处理才成为可能。所有满足容许性条件(从-INF到+INF积分为零)的函数,都可以成为小波。小波作为尺度膨胀和空间移位的一组函数也就诞生了。但连续取值的尺度因子和平移因子,在时域计算量和频域的混叠来说,都是极为不便的。用更为专业的俗语,叫再生核。也就是,对于任何一个尺度a和平移因子b的小波,和原信号内积,所得到的小波系数,都可以表示成,在a,b附近生成的小波,投影后小波系数的线性组合。这就叫冗余性。这时的连续小波是与正交基毫无关系的东西,它顶多也只能作为一种积分变换或基。但它的显微镜特点和相似性检测能力,已经显现出来了。为了进一步更好的将连续小波变换离散化,以下步骤是一种有效方法。第一步,尺度离散化。一般只将a二进离散化,此时b 是任意的。这样小波被称为二进小波。第二步,离散b。怎么离散化呢b取多少才合适呢于是,叫小波采样定理的东西,就这样诞生了。也就是小波平移的最小距离(采样间隔),应该大于二倍小波基的最高频率(好像类似,记不清了)。所以b取尺度的整数倍就行了。也就是越胖的小波,对应频谱越窄,平移量应该越大,采样间隔越大。当然,第一二两步的频域理解,即在满足频域窗口中心是3倍的频域窗口半径的前提下,频域就在统计上是完美二分的。(但很多小波满足不了这个条件,而且频域窗口能量不,所以只是近似二分的).这时的小波变换,称为离散二进小波变换.第三步,引入稳定性条件.也就是经过变换后信号能量和原信号能量有什么不等式关系.满足稳定性条件后,也就是一个小波框架产生了可能.他是数值稳定性的保证.一个稍弱的稳定条件,就是

小波变换学习心得

小波变换学习心得 第一章什么是小波变换 1从傅里叶变换到小波变换 1.1 短时傅里叶变换 为了克服傅里叶变换中时域和频域不能兼容的缺点,短时傅里叶变换把一个时间信号变为时间和频率的二维函数,它能够提供信号在某个时间段和某个频率围的一定信息。这些信息的精度依赖于时间窗的大小。短时傅里叶变换的缺点是对所有的频率成分,所取的时间窗大小相同,然而,对很多信号为了获得更精确的时间或频率信息,需要可变的时间窗。 1.2 小波变换 小波变换提出了变换的时间窗,当需要精确的低频信息时,采用长的时间窗,当需要精确的高频信息时,采用短的时间窗,图1.3 给出了时间域信号、傅里叶变换、短时傅里叶变换和小波变换的对比示意图。 由图1.3看出,小波变换用的不是时间-频率域。而是时间-尺度域,尺度越大,采用越大的时间窗,尺度越小,采用越短的时间窗,即尺度与频率成反比。 1.2 连续小波变换 小波是一个衰减的波形,它在有限的区域里存在(不为零),且其均值为零。图1.4是一个Daubechies小波(db10)与正弦波的比较。 正弦波:随时间无限振动的光滑波形,小波变换:尖锐变化而且是无规则的波形。因此小波能更好的刻画信号的局部特性。 在数学上,傅里叶变换的公式为

()()j t F f t e dt ωω+∞ --∞ =? 连续小波变换(Continue Wavelet Transform )的数学表达式 ()(),,a b a b CWT f t t dt ψ+∞ -∞ =? ()12 ,a b t b t a a ψψ--?? = ??? 式中,()t ψ为小波;a 为尺度因子;b 为平移参数。图1.6是小波变换的示意图。由图看出,小波变换给出了在各个时刻信号是由哪些尺度的小波构成。 小波中的尺度因子的作用是将小波在保持完全相似条件下“拉伸”或者“压缩”,图1.7给吃了尺度因子的“拉伸”和“压缩”作用。 小波中的平移参数,是简单地将波形沿时间轴平移。

小波分析考试题及答案

一、叙述小波分析理论发展的历史和研究现状 答:傅立叶变换能够将信号的时域和特征和频域特征联系起来,能分别从信号的时域和频域观察,但不能把二者有机的结合起来。这是因为信号的时域波形中不包含任何频域信息,而其傅立叶谱是信号的统计特性,从其表达式中也可以看出,它是整个时间域内的积分,没有局部化分析信号的功能,完全不具备时域信息,也就是说,对于傅立叶谱中的某一频率,不能够知道这个频率是在什么时候产生的。这样在信号分析中就面临一对最基本的矛盾——时域和频域的局部化矛盾。 在实际的信号处理过程中,尤其是对非常平稳信号的处理中,信号在任一时刻附近的频域特征很重要。如柴油机缸盖表明的振动信号就是由撞击或冲击产生的,是一瞬变信号,单从时域或频域上来分析是不够的。这就促使人们去寻找一种新方法,能将时域和频域结合起来描述观察信号的时频联合特征,构成信号的时频谱,这就是所谓的时频分析,亦称为时频局部化方法。 为了分析和处理非平稳信号,人们对傅立叶分析进行了推广乃至根本性的革命,提出并开发了一系列新的信号分析理论:短时傅立叶变换、时频分析、Gabor 变换、小波变换Randon-Wigner变换、分数阶傅立叶变换、线形调频小波变换、循环统计量理论和调幅—调频信号分析等。其中,短时傅立叶变换和小波变换也是因传统的傅立叶变换不能够满足信号处理的要求而产生的。 短时傅立叶变换分析的基本思想是:假定非平稳信号在不同的有限时间宽度内是平稳信号,从而计算出各个不同时刻的功率谱。但从本质上讲,短时傅立叶变换是一种单一分辨率的信号分析方法,因为它使用一个固定的短时窗函数,因而短时傅立叶变换在信号分析上还是存在着不可逾越的缺陷。 小波变换是一种信号的时间—尺度(时间—频率)分析方法,具有多分辨率分析(Multi-resolution)的特点,而且在时频两域都具有表征信号局部特征的能力,使一种窗口大小固定不变,但其形状可改变,时间窗和频率窗都可以改变的时频局部化分析方法。小波变换在低频部分具有较高的频率分辨率和较低的时间分辨率。在高频部分具有较高的时间分辨率和较低的频率分辨率,很适合于探测正常信号中夹带的瞬态反常现象并展示其成分,所以被誉为分析信号的显微镜。 小波分析最早应用在地震数据压缩中, 以后在图像处理、故障诊断等方面取得了传统方法根本无法达到的效果. 现在小波分析已经渗透到了自然科学、应用

小波变换 完美通俗解读2

这是《小波变换和motion信号处理》系列的第二篇,深入小波。第一篇我进行了基础知识的铺垫,第三篇主要讲解应用。 在上一篇中讲到,每个小波变换都会有一个mother wavelet,我们称之为母小波,同时还有一个father wavelet,就是scaling function。而该小波的basis函数其实就是对这个母小波和父小波缩放和平移形成的。缩放倍数都是2的级数,平移的大小和当前其缩放的程度有关。 还讲到,小波系统有很多种,不同的母小波,衍生的小波基就完全不同。小波展开的近似形式是这样: 其中的就是小波级数,这些级数的组合就形成了小波变换中的基basis。和傅 立叶级数有一点不同的是,小波级数通常是orthonormal basis,也就是说,它们不仅两两正交,还归一化了。 我们还讲了一般小波变换的三个特点,就是小波级数是二维的,能定位时域和频域,计算很快。但我们并没有深入讲解,比如,如何理解这个二维?它是如何同时定位频域和时域的? 在这一篇文章里,我们就来讨论一下这些特性背后的原理。 首先,我们一直都在讲小波展开的近似形式。那什么是完整形式呢?之前讲到,小波basis的形成,是基于基本的小波函数,也就是母小波来做缩放和平移的。但是,母小波并非唯一的原始基。在构建小波基函数集合的时候,通常还要用到一个函数叫尺度函数,scaling function,人们通常都称其为父小波。它和母小波一样,也是归一化了,而且它还需要满足一个性质,就是它和对自己本身周期平移的函数两两正交: 另外,为了方便处理,父小波和母小波也需要是正交的。可以说,完整的小波展开就是由母小波和父小波共同定义的。

几种时频分析综述1——傅里叶变换和小波变换

几种时频分析方法综述1——傅里叶变换和小波变换 夏巨伟 (浙江大学空间结构研究中心) 摘要:传统的信号理论,是建立在Fourier 分析基础上的,而Fourier 变换作为一种全局性的变化,其有一定的局限性。在实际应用中人们开始对Fourier 变换进行各种改进,小波分析由此产生了。小波变换与Fourier 变换相比,是一个时间和频域的局域变换因而能有效地从信号中提取信息,通过伸缩和平移等运算功能对函数或信号进行多尺度细化分析(Multiscale Analysis ),解决了Fourier 变换不能解决的许多困难问题。本文对傅里叶变换和小波变换进行了详细介绍,并用算例分析指出了两者的差别。 关键词:傅里叶变换;小波变换;时频分析技术; 1 傅里叶变换(Fourier Transform ) 1 2/201 22/0()()()()1()()()(::::)N j nk N ft N ft j nk N n H T h kT e H f h t e d DFT FT IFT IDFT t NT k h t H f e dt h nT H e N NT ππππ--∞ --∞∞--∞?=??=??????????→????=?=??? ∑??∑离散化(离散取样) 周期化(时频域截断) 2 小波变换(Wavelet Transform ) 2.1 由傅里叶变换到窗口傅里叶变换(Gabor Transform(Short Time Fourier Transform)/) 从傅里叶变换的定义可知,时域函数h(t)的傅里叶变换H(f )只能反映其在整个实轴的性态,不能反映h (t )在特定时间区段内的频率变化情况。如果要考察h(t)在特定时域区间(比如:t ∈[a,b])内的频率成分,很直观的做法是将h(t)在区间t ∈[a,b]与函数[][]11,t ,()0,t ,a b t a b χ?∈?=? ∈??,然后考察1()()h t t χ傅里叶变换。但是由 于1()t χ在t= a,b 处突然截断,导致中1()()h t t χ出现了原来h (t )中不存在的不连 续,这样会使得1()()h t t χ的傅里叶变化中附件新的高频成分。为克服这一缺点, D.Gabor 在1944年引入了“窗口”傅里叶变换的概念,他的做法是,取一个光滑的函数g(t),称为窗口函数,它在有限的区间外等于0或者很快地趋于0,然后将窗口函数与h(t)相乘得到的短时时域函数进行FT 变换以考察h(t)在特定时域内的频域情况。 22(,)()()()()(,)ft f ft f STFT ISTF G f h t g t e dt h t df g t G f e d T ππτττττ +∞ --∞ +∞+∞ -∞ -∞ =-=-??? ::

小波变换去噪基础地的知识整理

1.小波变换的概念 小波(Wavelet)这一术语,顾名思义,“小波”就是小的波形。所谓“小”是指它具有衰减性;而称之为“波”则是指它的波动性,其振幅正负相间的震荡形式。与Fourier变换相比,小波变换是时间(空间)频率的局部化分析,它通过伸缩平移运算对信号(函数)逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号分析的要求,从而可聚焦到信号的任意细节,解决了Fourier变换的困难问题,成为继Fourier变换以来在科学方法上的重大突破。有人把小波变换称为“数学显微镜”。 2.小波有哪几种形式?常用的有哪几种?具体用哪种,为什么? 有几种定义小波(或者小波族)的方法: 缩放滤波器:小波完全通过缩放滤波器g——一个低通有限脉冲响应(FIR)长度为2N和为1的滤波器——来定义。在双正交小波的情况,分解和重建的滤波器分别定义。 高通滤波器的分析作为低通的QMF来计算,而重建滤波器为分解的时间反转。例如Daubechies和Symlet 小波。 缩放函数:小波由时域中的小波函数 (即母小波)和缩放函数 (也称为父小波)来定义。 小波函数实际上是带通滤波器,每一级缩放将带宽减半。这产生了一个问题,如果要覆盖整个谱需要无穷多的级。缩放函数滤掉变换的最低级并保证整个谱被覆盖到。 对于有紧支撑的小波,可以视为有限长,并等价于缩放滤波器g。例如Meyer小波。 小波函数:小波只有时域表示,作为小波函数。例如墨西哥帽小波。 3.小波变换分类 小波变换分成两个大类:离散小波变换 (DWT) 和连续小波转换 (CWT)。两者的主要区别在于,连续变换在所有可能的缩放和平移上操作,而离散变换采用所有缩放和平移值的特定子集。 DWT用于信号编码而CWT用于信号分析。所以,DWT通常用于工程和计算机科学而CWT经常用于科学研究。 4.小波变换的优点 从图像处理的角度看,小波变换存在以下几个优点: (1)小波分解可以覆盖整个频域(提供了一个数学上完备的描述) (2)小波变换通过选取合适的滤波器,可以极大的减小或去除所提取得不同特征之间的相关性 (3)小波变换具有“变焦”特性,在低频段可用高频率分辨率和低时间分辨率(宽分析窗口),在高频段,可用低频率分辨率和高时间分辨率(窄分析窗口) (4)小波变换实现上有快速算法(Mallat小波分解算法) 另: 1) 低熵性变化后的熵很低; 2) 多分辨率特性边缘、尖峰、断点等;方法, 所以可以很好地刻画信号的非平稳特性 3) 去相关性域更利于去噪; 4) 选基灵活性: 由于小波变换可以灵活选择基底, 也可以根据信号特性和去噪要求选择多带小波、小波包、平移不变小波等。 小波变换的一个最大的优点是函数系很丰富, 可以有多种选择, 不同的小波系数生成的小波会有不同的效果。噪声常常表现为图像上孤立像素的灰度突变, 具有高频特性和空间不相关性。图像经小波分解后可得到低频部分和高频部分, 低频部分体现了图像的轮廓, 高频部分体现为图像的细节和混入的噪声, 因此, 对图像去噪, 只需要对其高频系数进行量化处理即可。 5.小波变换的科学意义和应用价值

相关主题
文本预览
相关文档 最新文档