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

peng2v

我在学习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