马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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 |