声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1103|回复: 1

[HHT] 边际谱问题

[复制链接]
发表于 2008-10-30 19:34 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
在用论坛里求边际谱的程序时,出现了一些小问题。我的程序是这样的
%N为信号的采样点数。
[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,2);
       bjp1(k)=sum(im(:,k))*1/fs*1/T;  %T=N*1/fs
   end
   f=(0:N-3)/N*(fs/2);
   plot(f,bjp1,'r');
   title('边际谱');
   xlabel('频率HZ');
   ylabel('幅值');
可是运行出来的结果不对,在不该出现有幅值的地方也出现了,对应频率出现的机率和HHT谱出现的不对应。im是二维的谱信号,bjp1(k)=sum(im(:,k))*1/fs*1/T这样可不可以求出谱值的?我也不懂,把k=1:size(im,2);改为k=1:size(im,1);结果又是不一样,都不知道是哪里出了问题,请高手们帮帮忙啊!!!
回复
分享到:

使用道具 举报

 楼主| 发表于 2008-10-30 20:34 | 显示全部楼层
该信号的HHT谱图和边际谱如附件所示,信号的取样率为50hz,N=36001,为什么在20-25hz这一部分,HHT谱图和边际谱图不对应呢,出现了没有值,这是为什么呀?sig为信号。
    t=sig(:,1);
    pm=sig(:,2);
    N=length(t);
    fs=50;
    T=N*1/fs;
   % imf = emd(pm);
   [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);
   figure;
for k=1:size(im,2);
       bjp1(k)=sum(im(:,k))*1/fs*1/T;
   end
   NN=size(im,2);
   f=(0:N-3)/N*(fs/2);
   plot(f,bjp1);
   title('边际谱');
   xlabel('频率HZ');
   ylabel('幅值');

HHT谱图

HHT谱图

边际谱图

边际谱图
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-2 21:26 , Processed in 0.068415 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表