EEMD例子学习后感想

王跃峰

         近几日，有幸学习了一下台湾国立中央大学黄院士开发的EEMD工具。按照指导网上下载的指导手册（查看附件），前面4部分程序都能运行。如果换成了自己的数据一定要根据提取的IMF个数来适当修改指导手册命令，因为指导手册的命令默认的是生成8个IMF分量。 此外，还发现在指导手册最后一部分（分量显著性检验）中，总会提示undefined function“significance”，仔细检查发现下载的包里确实缺少了“significance.m”这个文件，但是这个文件也可以在国立中央大学的网上下载到，把这个code粘贴到你的文件夹里就可以运行这一步了。
     另外，我看网上大家都建议把下载的包解压到MATLAB安装路径下的toolbox文件夹内，我也是这么做的，但是我觉得不放进去也可以计算，大家可以尝试一下。
     附件还上传了一篇有助于大家理解EEMD的中文文章，写的挺好，大家可以辅助理解下。
     另外我还有个地方不明白：计算出的IMF分量的瞬时频率是一个过程线（见附件中的图片），如何最终确定为一个数值，是看峰值吗？还是根据几个峰值求平均，然后再求倒数，得到周期，不是很懂。有明白的请赐教一下。



来自群组: 南气院校友总会

李白

好东西！感谢楼主

aridlhy

好东西！感谢楼主

rceclx

学习一下！谢谢楼主！

surion

楼主您好~~求问为什么我把EEMD解压到了TOOLBOX以后，
输入中央大学的这部分代码以后提示是输入参数的数目不足，出不了图。我就是下载了指导手册里面的数据。
无法EEMD分解~  
先谢谢楼主的回答啦~

>> year=gsta(:,1);

>> inData=gsta(:,2);

>> rslt=eemd(inData,0,1);

>> plot(year,rslt(:,2));

>> hold on;

>> plot(year,rslt(:,3)-0.3);

>> plot(year,rslt(:,4)-0.6);

>> plot(year,rslt(:,5)-0.9);

>> plot(year,sum(rslt(:,6:8),2)-1.3,'r-');

>> hold off

>> set(gca,'yTickLabel',[]);

>> axis([1850 2010 -1.8 0.3]);

>> xlabel('year');

surion

啊我弄好了。set path以后这个第二步我画出来图了

thukd35

关于如何确定一个IMF的特征频率的方法，Huang在2006年有一个美国专利，COMPUTING FREQUENCY BY USING GENERALIZED ZERO-CROSSING APPLIED TO INTRINSIC MODE FUNCTIONS，不过我也没细看，不知道对你有没有帮助~

wuukaii

学习一下

wode_q_x

学习了,感谢楼主无私的分享!

浪子

跟原网上一模一样