求助hilbert边际谱

peng2v · 发表于 2009-11-2 15:15

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

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

x

我在使用破凰的求Hilbert谱的时候，得到的谱线很诡异，使用的hhspectrum也是http://perso.ens-lyon.fr/patrick.flandrin/emd.html 这个网站下载的程序，实在是搞不明白...谢谢大家的帮助

clear;
fs=1000; %fs为采样频率;
N=1000;    %采样点数
t=1/fs:1/fs:1;
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*t);
y=[y1;y2;zeros(size(y1))]; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);  %标记1
[E,tt1]=toimage(A,fa,tt,length(tt));  %标记2
E=flipud(E);
for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');

peng2v · 发表于 2009-11-2 15:16

帖图上不明显，在横坐标为0.5的地方有一根谱线，不明白，为什么程序一样，得到的结果不一样啊.......

bcyangbc · 发表于 2009-11-2 20:51

http://forum.vibunion.com/forum/viewthread.php?tid=37688
这个帖子讨论得很清楚
我也是用这个程序

peng2v · 发表于 2009-11-3 10:08

谢谢您的回复，我把那个帖子基本看一遍了，EMD程序包有2005版和2007版两个版本，我用的是2007版的，程序也是拷的论坛中的原程序，为什么我做的边际谱图和大家的不一样啊，感觉很迷茫啊...仿真没出来，就不能往下走了，很着急啊，所以把程序步骤中的计算结果发上来了，希望各位大哥给看下，谢谢！！

peng2v · 发表于 2009-12-21 16:45

这个文件要在matlab7.1 环境下运行才能得到正确的答案，我在7.0环境下运行时就出现了错误，遇到同样问题的同学看看了... 呵呵

daxue123 · 发表于 2012-7-22 23:11

遇到过同样问题

tjujhy · 发表于 2012-12-1 11:47

学到了很多，谢谢！！！

华电机械 · 发表于 2013-5-18 17:16

bcyangbc 发表于 2009-11-2 20:51

登录/注册后可看大图

http://www.chinavib.com/forum/viewthread.php?tid=37688
这个帖子讨论得很清楚
我也是用这个程序

我的结果是这样定义

bzp小鹏 · 发表于 2014-6-6 09:46

这个程序不错，挺好的

账号		自动登录	找回密码
密码			我要加入

[HHT] 求助hilbert边际谱

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

本帖被以下淘专辑推荐:

回复楼主 peng2v 的帖子

回复板凳 bcyangbc 的帖子

问题解决

[HHT] 求助hilbert边际谱

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

本帖被以下淘专辑推荐:

回复 楼主 peng2v 的帖子

回复 板凳 bcyangbc 的帖子

问题解决

回复楼主 peng2v 的帖子

回复板凳 bcyangbc 的帖子