blusun_2008 发表于 2008-11-2 18:41

爆破波形的边际谱问题！

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

    disp_hhs(im,[],fs);
    figure;
for k=1:size(im,1);                ％size（im）=
      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('幅值');

吃书的老虎 发表于 2008-11-3 20:26

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

blusun_2008 发表于 2008-11-4 15:29

爆破的频率一般是在几十赫兹左右，为什么我画的这个这么低呢？这个爆破信号的HHT谱图主要频率是集中在5－15hz左右的。
页: [1]
查看完整版本: 爆破波形的边际谱问题！