小波变换的本质148.pdf

《小波变换的本质148.pdf》由会员分享，可在线阅读，更多相关《小波变换的本质148.pdf（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、为了应付老板的的一个任务而收集了几篇相关文章!我是搞电力系统故障波形分析的，正上研二，导师定的方向是用小波 变换进行信号的消噪及波形奇异点检测.出于研究方向的需要从去年 年底开始接触小波.毕竟是工科出身，学起小波来觉得难度很大.不 夸张地说常有学不下去的感觉.硬着头皮看了一段时间，终于觉得有 点眉目，现将我从信号奇异性方面的理解写出来，请各位同仁批评指 正，并希望能对刚接触小波的朋友有点帮助！1学习小波变换所需的基础知识 由于小波变换的知识涵盖了调和分析，实变函数论，泛函分析及矩 阵论，所以没有一定的数学基础很难学好小波变换.但是对于我们工 科学生来说，重要的是能利用这门知识来分析所遇到的问题

2、.所以个 人认为并不需要去详细学习调和分析，实变函数论，泛函分析及矩阵 论等数学知识.最重要是的理解小波变换的思想！从这个意义上说付 立叶变换这一关必需得过！因为小波变换的基础知识在付立叶变换中 均有提及，我觉得这也就是很多小波变换的书都将付立叶分析作为其 重要内容的原因.所以我认为学习小波应从V数字信号处理 中的付 立叶分析开始.当然也可从V信号与系统这本书开始.然后再看杨 福生老师的小波变换书.个人觉得他的书最能为工科学生所接受.2信号的分解 付立叶级数将周期信号分解为了一个个倍频分量的叠加，基函数是正 交的，也就是通常所说的标准正交基.通过分解我们就能将特定的频 率成分提取出来而实现特定

3、的各种需要，如滤波，消噪等.付立叶变 换则将倍频谱转换为了连续谱，其意义差不多.小波变换也是一种信 号分解思想：只不过它是将信号分解为一个个频带信号的叠加.其中 的低频部分作为信号的近似，高频部分作为信号的细节.所谓的细节 部分就是一组组小波分量的叠加，也就是常说的小波级数.3小波变换的时频分析思想 付立叶变换将信号从时域变换到了频域，从整体上看待信号所包含的 频率成分.对于某个局部时间点或时间段上信号的频谱分析就无能为 力了，对于我们从事信号的奇异性检测的人来说，付立叶变换就失去 了意义（包括加窗付立叶变换）.因为我们要找的是信号的奇异点（时 域方面）和奇异点处所包含的频带（频域方面）也就是

4、说需要一种时 频分析方法.当然能有纯时域的分析方法更好！（据说数学形态学能 达到这种效果）.小波变换之所以可以检测信号的奇异点，正在于它 的小.因为用小的波去近似奇异信号要比正弦波要好的多.4小波变换的实质 小波变换的公式有 内积形式和卷积 形式，两种形式的实质都是一样 的.它要求的就是一个个小波分量的系数也就是权.其直观意义 就是首先用一个时窗最窄，频窗最宽的小波作为尺子去一步步地 量信号，也就是去比较信号与小波的相似程度.信号局部与小波越 相似，则小波变换的值越大，否则越小！当一步比较完成后，再将尺 子拉长一倍，又去一步步地比较，从而得出一组组数据.如此这般循 环，最后得出的就是信号的小波

5、分解（小波级数）.当然这只是一种 粗略的解释.5连续小波变换，二进小波变换与离散小波变换的关系 当尺度及位移均作连续变化时，可以理解必将产生一大堆数据，作实 际应用时并不需要这么多的数据，因此就产生了离散的思想.将尺度 作二进离散就得到二进小波变换，同时也将信号的频带作了二进离 散.当觉得二进离散数据量仍显大时，同时将位移也作离散就得到了 离散小波变换！6 MALLAT 算法的意义 想必大家都注意到，小波变换是以内积或卷积的形式实现的，这给数 值计算带来了不利之处，因为用计算机作数值积分其计算量 大.MALLAT 算法则解决了这一问题，它不涉及小波的具体形式，只 是对系数进行操作！其计算也就是

6、用高通及低通滤波系数与小波系数 作卷积.因为作信号处理时，我们往往并不关心小皮的具体形式，更 为关心小波系数.需提出的是该算法仅适用于正交小波如果小波不是 正交的（如E样条小波）则算法失效！7小波变换的模极大值及其意义 对于我们搞信号奇异性检测的人来说，小波变换最重要的应用就是用 模极大值定值奇异点.我觉得模极大值可以从两个方面去理解：第一,从直观角度，上文已说明小波变换的实质就是一种度量波形相似程度 的方法.信号与小波越相似，则小波系数越大.这也就可理解为出现 了小波变换的模极大值.因为当信号出现奇异点时，或是间断点，或 是一阶导数不连续点，其在各个尺度下都将必然出现大的小波系 数.从而可以

7、定位奇异点！第二个方面从小波的取法来看，当小波取 为光滑函数一阶导数或二阶导数时，从公式可以推导出小波变换将出 现模极大值点或是过零点 也就是很多书上说的模极大值检测和零交 叉检测这些可以查书看！我只谈谈连续小波变换，对于离散的也有同样的 argument。小波函 数的 dilation 和 translation 是这样一个形式：1/sqrt|s|psi（x-u）/s）,s 是 scale,u 是该小波 atom 的 center。由于根据定义，小波的积分是 0，也就是说小波函数的傅立叶变换在 零点为零。再有于小波函数的傅立叶变换一般是连续的（比如如果小 波是属于 L_1 的），这样在 0

8、的一个小临域里面，小波的傅立叶变换 很接近零，这也就是说小波函数的傅立叶变换可以看成某个高通滤波 器的 transfer function，这样小波变换 W（f）实际是在 measure 该 函数 f 在 u 点附近的 variation。从这个角度看的话，如果小波的宽 度很大（对应尺度 s 很大），该函数在该小波的窗口下的 variation 就很大；如果小波的宽度小（对应尺度 s 小），则函数在该小波的窗 口下的 variation 就相对比较小（除非信号是 fractal，呵呵）。小波情结 到了小波版很久，总觉得应该写些什么。这篇文章也就献给那些所有 正在研究或即将研究小波的同学、老师和

9、科研人员们。这是篇与技术 无关的文章，撰写的是我对小波的感受。从我开始接触小波，研究小 波，到迷恋小波的真实记录。因此，我把它起名叫小波情结。刚开始，接触小波的时候在研一。关于老师布置的从频域构造一种小 波的作业开始。后来我才知道，这种小波本质上就是 meyer 小波。当 时，就一个字，嫩。实际上就是对小波毫无所知。脑子里就是一叠的 公式。正交条件，容许条件等一大堆，与概念理解相差甚远的东东。但，还是乐此不疲的编程。总想看看，我亲手缔造的小波长的是什么 样，也有些略带孩子气的，想把它作为桌面和自己的酷酷头像之类的 欲望。于是，十一的头三天，我基本上闭门造车。当时，我用的是 matlab，也是我

10、最后得到哭笑不得结果的直接的助手与帮凶。因 为构造的过程的起始，我就把函数离散化了。紧接着就是平移，对乘，积分，抽取，插值，dsp 里的一套 trick 把我搞得叫苦不迭。程序也 累计到了 1000 行左右。当时，最可恨的就是对点，由于 dsp 下标的 1，2，3 离散化，所以我也就用手指开始傻傻的算。连续几天晚上鏖 战，终于在 3 号的晚上。通过 IFFT 后，美妙的波形出来了。注意，美妙和丑陋只有一步之遥。这是我的对小波的第一课体会。当我一看 屏幕，疯了，彻底疯了。一个 DELTA函数类似的波形，就在我眼前。心想：忙乎了三天，整了个 DELTA 函数出来。这难道就是回报吗？别 急，小波是紧

11、支撑的啊。概念上对头，一定是取点的问题。我便拿起 MATLAB 自带的照妖镜（放大镜）一看，呵呵，一个差强人意的波形 就在我眼前了。我当时大喊一声，爽，那时已经凌晨 2:00。第一次 的经历，对我来说收获很丰。然后，第二次挑战，则是彻底改变我对 小波是个深不可测的家伙的看法。这次作业，就是用刚才构造的小波，做消噪。我这次，又一次的想起，爱情格言：我心灵的古堡经不起你 轻轻的一击。女生问：结果呢？回答：碎了。一个如此，不精确的波 形，怎么能消噪呢？而且，当时老师要用连续小波的方法。也就是内 积求和的方法。我和同学，首先合作，用 mathmatic 做了个好一点 的波形。因为，除最后一步，反傅里叶

12、变换外，其他都是解析的。然 后，一个困扰我许久的问题产生了。一个函数可以由无穷多个小波的 膨胀和伸缩叠加起来。那么，我把函数从-inf 到+inf 积分，假设 函数有直流分量，所以积分不为零。但是小波，积分却为零。这不是 矛盾吗？后来，也就是研二我才知道，有些时候积分后不可以交换。还有，其实有限的小波逼近，必须加上尺度函数才可以。但当时，我 们只是采用了把小波的支撑取宽的办法解决了此问题。但，我由于不 太喜欢这种方法的冗长和费时，所以想令辟蹊径。于是，mallet 一 个令我崇拜的算法，终于在我阅读超星的时候，跳在了我生命里。首 先，便是看冗长的证明，勉强理解了。当看到滤波器组的解释后，我 开

13、始豁然开朗。这是我熟悉的 dsp 概念。因此，我花了一晚上，把这 个算法彻底搞懂了。但概念的理解和程序的成功编制，还是有一小步，就是这一步，使无数英雄竟折腰。我的幸运之神便是 MATLAB 里的 DEMO 那个里面，有一个详细的算法解释。并且从哪里我知道了些怪 怪的函数。WKEEP(),DYADDOWN(),DYADUP 等。而且，又一个问题，理论和实际差别产生了。这个问题甚至现在，还困扰着很多的小波工 作者。一个长度为 100的信号，分解后理论上高频 50,低频 50。但用 卷积算法，假设滤波器长度为 10。因此总长度 109，做抽取后长度 55。多了 5。这怎么办呢。我去问了很多老师，回答

14、都一样。就是 MATLAB 里用的函数 WKEEP(。)把两头丢掉。当时我勉强接受了这个结果。但 始终有个概念，小波变换就是正交变换，它和傅里也变换一样，一定 可以写成正交阵的形式。第二次作业的完成，我的小波课结束了。但 我的小波情结还在继续。关于，正交阵的猜想还在困扰着我。一本电 磁场和小波结合的外文书籍，帮助了我。圆周卷积的概念，历历在目。是呀，卷积对着傅里叶变换，而圆周卷积对应着离散傅里叶变换。这 就是连续与离散的区别和联系啊。于是我用 db 小波，构造了一个完 全正交的矩阵。当我把这个矩阵和它的转置相乘的时候，单位阵出来 了。那天，我高兴得流泪。最终，我把圆周卷积用快速傅里叶变换实 现

15、出来。今后的日子，我便觉得，思维的水再也关不住了。步步为营，我实现了 db 小波的时域构造，采用矩阵特征向量法和 casade 理论两 种解法，我都成功了。慢慢的我开始醉心于消失矩，开始懂得框架，开始懂得双正交。然后就是，PR 条件，二代小波，小波插值，因子 化，等等。于是，我也在研学一边和大家交流，一边阅读大量书籍和 文献，而且实现里面的每一个例子和思想。当我们还在觉得自己懂点 小波的时候，美国人已把它用于指纹压缩，产生了巨大的经济和社会 效益；当我们，还在对二代不屑一顾的时候，一个叫 JPEG200C 标准 的东东，彻底给我们上了一课。当我们，还在国家著名期刊上，打着 错误的提升公式的时候

16、，当我们，还在为些不值一提的程序保密的时 候，一个叫各相异性小波的东东又开始蠢蠢欲动。看看那些大师们吧，看看他们的态度，再看看我们，我们努力的够吗。你说看不懂文献，我就要问你，你看了一遍，十遍，还是一百遍呢？如果说你认为是高 手，你是否写了超过 10 万行以上小波的代码，看了 10 本以上的书，100 篇的文献，实现里面所有的例子和思想了呢。我们差得很远。但 是我们服气吗，我们认输了吗，我们不再努力了吗。什么时候有中国 的 JPEG2008 呢，什么时候我们能毫无保留的进行坦诚的交流和无私 的分享呢，什么时候我们把学术的铜臭拨掉，把做小波看成一次和上 苍对话的机会，和真理的交锋呢。我始终在问自

17、己这些问题。为关于 学术的单纯的问题。我找到了答案。在研学上。因此，我毫无保留的 帮助大家，同时也在修正自己。我开始变得勤奋，开始编每个需要的 程序，而基本上不用 MATLA 提供的任何函数，除非是概念性验证。我和所有那些从事这小波事业的人们一样，为实现这上面单纯而坚定 的疑问而不停奋斗。我和你们一样，是一个奔跑着，一个向着小波的 科学和真理殿堂不辞辛劳的奔跑者。即使路上满是荆棘，即使我们会 暂时的迷失方向，但我相信小波这朵最美丽的奇葩，会以它最美丽的 身姿，最沁人的芬芳，指引着我们。让我们结伴而行吧！我是引路者，也是跟随者。和你们一样，怀揣梦想，一起努力。虽然汗流浃背，虽 然荆棘满身，但成功终会来的。我悯悯驾信，我正在奔跑，像春天里 的孩子。