声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1041|回复: 0

[HHT] 大家来帮我看看这个程序错在哪啊....谢谢了

[复制链接]
发表于 2009-11-2 14:58 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我在学习hilbert-huang方法,在做信号的hilbert谱的时候,遇到了些问题,下面是我在论坛上找到位哥哥的程序,用Grilling提供的开源程序(http://perso.ens-lyon.fr/patrick.flandrin/emd.html ),附件中是我得到的图,怎么hilbert谱基本没值啊,不该是15 10 0.1 5 这些值的一些有些波动的点吗?  
T=0:0.01:1;
x=2*sin(30*pi*T)+4*sin(20*pi*T).*sin(2*pi*1/10*T)+sin(2*pi*5*T);
plot(T,x)  %绘出仿真波形图,如图1
OPTIONS.DISPLAY=2;
OPTIONS.T=0:0.01:1;
imf=emd(x,OPTIONS);
emd_visu(x,0:0.01:1,imf)  %绘出各个IMF波形,如图2
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im,[],100);
colormap(flipud(gray))  %绘出Hilbert谱,如图3

原信号

原信号

IMF图

IMF图

得到的Hilbert谱图

得到的Hilbert谱图
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 06:26 , Processed in 0.077923 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表