爆破波形的边际谱问题！

blusun_2008

我用G. Rilling 2007年写的版本求出爆破波形的HHT谱图，然后用论坛上破凰写的求边际谱的程序求爆破波形的边际谱图，但是得出的结果是消减的太快，不知是为什么呢，哪位高手可以进来帮忙，帮看看呀！
这是我的程序：
sig=rsac(a);
    t=sig(:,1);               ％信号的时域采样点 ＝36001
    pm=sig(:,2);           ％幅值
    N=length(t);
    fs=50;                      ％采样频率
    T=N*1/fs;
   [imf,ort,nbits] = emd(pm);
   emd_visu(pm,t,imf);
[A,f,tt] = hhspectrum(imf(1:end-1,:));
[im,tt,ff] = toimage(A,f,tt);
  
    disp_hhs(im,[],fs);
    figure;
  for k=1:size(im,1);                ％size（im）=[400  35999]
      bjp1(k)=sum(im(k,:))*1/fs;   
   end
   NN=size(im,1);
    f=(0:NN-1)/NN*(fs/3);
    plot(f,bjp1,'r');
    xlabel('频率HZ');
    ylabel('幅值');

吃书的老虎

图中看出，能量主要集中在0～1.5HZ。不知你正确的结果应该是什么样子？

blusun_2008

爆破的频率一般是在几十赫兹左右，为什么我画的这个这么低呢？这个爆破信号的HHT谱图主要频率是集中在5－15hz左右的。