大家来帮我看看这个程序错在哪啊....谢谢了
我在学习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
=hhspectrum(imf);
=toimage(A,f);
disp_hhs(im,[],100);
colormap(flipud(gray))%绘出Hilbert谱,如图3
页:
[1]