复件hht变换的三种方法matlab ument

1、与大师的近距离接触听HHT创始人黄鄂授课有感从不在BBS上发文章，因为觉得自己没有什么独特而有价值的东西与大家分享。我是个HHT的初学者，似乎只有两三个月，有点信号处理的常识，太多不懂的地方，所以从这个论坛上收获不少。今天之所以写，是因为前几天天听了HHT 之父黄鄂的讲课，真的是感触颇深，一直拜读他的文章，没想到可以见到本人，还可以请教他一些问题，实在是幸运而不幸，不幸来源于自己实在研究的不够深入，以至于有了这样的机会，也没法问出什么尖锐的问题，有点不爽。黄鄂在我们学校一共讲了三天，由于某些内部的原因，我很晚才得到消息，只听了最后一天的课，不过就只这一天都让我受益匪浅，大师就是大师，现将大师的

2、一些观念整理如下（很不好意思的加进了我自己的理解，望大家见谅）：1 黄老师对于HHT的热爱，或者说对于自己所从事的研究的热爱，是最让我影响深刻，也最值得我们铭记在心的，光凭这点，即便忽略掉他名字后面长长的title，都值得我们对他的尊重。我一直不怎么相信名人简历里写的东西，国人忽悠的本事太利害了，可能是我经历得太多。俗话说，眼见为实，一点不假。黄老师说他总是在不停的思考，我信这不是刻意的强调，他在我们学校上课，从早9点讲到晚5点，中午休息3小时，这是我参加了的。晚上他还要思考我们学校项目的问题，我肯定他确实做到了，因为第二天，他就给出了建议，而且听者无不觉得有道理。所以这都是我所亲眼见到的实事

3、，绝不是因为看了他的简历，就对他盲目的崇拜。这位满头白发，但双眼散发着光芒的师者，是把自己提出的理论，把HHT当成自己的孩子一样，希望他能被大家接受，希望他能成才。黄老师说，每当他拿到一段数据，就想用HHT来分析，一直到分析出结果为止。他的这种不停滞的态度，真正的热情，太令我们感到惭愧了。我们以为从网上down几个程序，装模作样的改几行，画几个图，凑篇文章就以为自己在从事研究，再靠着这些东西去混硕士，混博士文凭。我一开始并没有写专业相关的东西，而是提了些老生常谈，或许这才是决定人与人差距的地方。2 黄老师对于Fourier 有着自己独到的见解。我们都知道，做信息处理，Fourier 是个万金油

4、，至少很多年的实践经验和很多过去的学者这样认为，以致于我们拿到一段信号，第一反应就是做FT变换。这点没有错，错在Fourier也不是万能的，明明用Fourier没有得到有用的结果，为什么不去思考这种方法的弊端。黄老师认为很多以前的说法都是基于Fourier就是真理的前提，例如我们常看到关于instantaneous frequency 的种种悖论 ，（基本都是参考cohen的）都是基于这个假设的，如果Fourier在某些情况下不成立，这些悖论也就不是悖论了。Fourier是个global, linear的分析工具，这点黄老师一直强调，具体我就不写了，大家可以参考他的论文。3 现在的市面上每年都

5、有大量信息处理，分析的文章、理论、方法出现，深究其原因，是因为信号，图像的处理并不像物理，化学等学科，有自己的完整而严谨的系统，也就是说他们并没有一条固定的路可以走，我感觉随机性很大。例如图像分割，往往是用这种方法也行，用那种也行，或者是同一方法，用在这幅图效果很好，换一幅就不能用了。信号处理也有同样的问题。4 基于3所提出的问题，我们做信号分析时要用到的一个强有力的工具就是我们的眼睛，这点黄老师在他的98年的文章中也有提到。好好的观察信号的原始图（一般是时间幅度图），思考曲线为什么会这样走，再去选择分析的工具。5 可能由于国内太急于求成的环境，很多人都是认死理。做信号处理的就只关心信号处理，

6、全然不顾这是关于什么变化的信号，这个什么的将是分析处理信号的突破口。打个比方，你想做旋转机械的故障检测，你拿到了采集到的转子的震动图，你就必须知道转子在正常的时候是什么样子，在发生故障时可能是什么样子，可能有哪几种故障，各是什么样的震动图。我相信这些知识的获得是可以请教搞机械的人，他们可能不懂信号处理，但没有这些常识，你就什么都做不了。6 见多才能识广。这点和5点其实是一个意思。还是那句话，要多看，多问，多想。特别是要向自己的非同行请教。你可能没有见过高性能的战斗机分行，但你知道鸟在分行的时候做的是什么动作，其实都可以从中收到启发。黄老师来我们学校参观了桥梁健康检测系统，他说看了这个以后，再去

7、看其他的桥就不会不知从何入手。7 关于现有文献中提到的HHT没有一个rigorous mathematics foundation，这点在论坛中也常有人提到，黄老师给出了一句名言：不能因为我不知道身体的消化过程原理，就不吃饭了。（出自一位天才物理学家，因为黄老师用英文说的，我不记得名字了，似乎和Maxwell Equation 有关系）所以我们也不能因为不清楚HHT的 rigorous mathematics foundation，就不去用它。在黄老师看来，这正是HHT的优势所在，HHT是一个超adaptive的方法，我们不清楚它到底是怎么回事，所以可以把它用在所有的信号上。也正因为它的ada

8、ptive的特性，我们在使用它的时候，也要灵活多变。其实这点也是基于3点的。打个比方：用EMD分解出了一系列分量，C1-C10，有时候要全部都用上，有时候其中一个分量的图比较接近我们的预想，这时就只用这个分量去画hilbert spectrum。这其中有很大一部分是来源于经验。所以还是那句，多做，做多了，就有feeling了。黄老师说他本身数学不行，所以申请到了一笔经费，请到Princeton 应用数学系的Daubechies做带头人去寻找HHT的 rigorous mathematics foundation。Daubechies就是那个定义Daubechies正交小波的人，也是个牛人，没想

9、到居然是位有点搞怪的美人。相信在这样一群人的带领下，找到HHT的 rigorous mathematics foundation指日可待。8 最后，还有一点，搞研究，听大师布道，英文真的非常非常重要。黄老师给我们讲课，开始还注意尽量使用中文，随着讲课的深入，他情不自禁的不停飚English，希望English不太好的同学真的要恶补下才行。（同去的同学硬要补上的，我觉得如果一直都是看英文论文的人是不会有这方面的问题）唠唠叨叨的写了一大堆，也不知道有没有传递出黄老师讲课的精髓，对大家有没有帮助。送给搞EMD或者HHT但没有下载到完整程序的朋友下载完整程序的方法：（考虑到程序更新问题，建议使用第二、

10、三种方法，flandrin就分别在2006年2月、2007年3月进行了更新）1. 本论坛（自己搜索一下）已经过时2. （适用于2007年3月之前）登录以下网站： HYPERLINK http:/perso.ens-lyon.fr/patrick.flandrin/emd.html t _blank http:/perso.ens-lyon.fr/patrick.flandrin/emd.html，这是flandrin的一个关于EMD和HHT的网页，进入后就可以看见主要的几个m文件，对只搞EMD的朋友这些程序就够用了，如果要进行HHT后续分析，则需要点击Time-Frequency ToolBo

11、x链接下载工具箱，具体位置在Time-Frequency ToolBox页面的download处，下载后得到tftb-0.1.tar.gz压缩包，压缩包中有个mfiles文件夹，所有程序就齐全了。注：instfreq 函数就在tftb-0.1.tar.gz里面；至于extr函数，如果阁下是在2006年2月之后从flandrin网站下载的程序，那么根本不需要这个m文件，因为该函数已集成到此版本的emd.m中，如果阁下是在2006年2月之前下载的，那么的确需要这个m文件，如果当时没有下载到，那么请下载此版本的emd.m，然后单独把extr函数提取出来；至于程序看不懂的问题，请按照emd.m文件中的

12、帮助信息下载相关的两篇文章，就是Huang和Flandrin写的关于EMD和HHT的经典著作Huang98年经典著作的下载地址： HYPERLINK /technologies/hht.htm t _blank /technologies/hht.htm3. （适用于2007年3月之后）登录上述网站下载压缩包，里面有详细说明。具体介绍可以参阅以下链接的帖子： HYPERLINK /forum/viewthread.php?tid=38891&extra=page%3D1 t _blank EMD估计有新进展注：程序看不懂的问题，请按照emd.m文件中的帮助信息下载相关的两篇文章（帮助信息列出的

13、5篇文章中的前两篇），就是Huang和Flandrin写的关于EMD和HHT的经典著作；另外，instfreq函数的实现请参阅 refguide.pdf 和 tutorial.pdf 两个文件中的该函数解释，并根据给出的参考文献查找相关资料阅读，这两个pdf文件在 hhspectrum.m 帮助信息关于时频工具箱(Time-Frequency Toolbox)的下载地址中可以找到EMD之包络线总结首先我要说这个论坛真的很不错，高手云集，让我这个菜鸟也对EMD有了一定的认识。我从众多的帖子里搜集到很多关于EMD包络线的资料，在这里想与大家分享一下，可能有的地方不是很准确，请各位版主多多批评指正。大概有三种方法1、在 emddec.