声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2011|回复: 8

[HHT] 边际谱

[复制链接]
发表于 2008-11-20 11:23 | 显示全部楼层 |阅读模式

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

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

x
我的数据通过在线采集的.采样频率是22050;
通过一段时间,我得到很多的数据150000;然后我就选取1:15000数据再EMD,边际谱.
N=15000;
如果按下面做
[A,fa,tt]=hhspectrum(y);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;  
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');

那么得到的频率很不准确.(选择的数据多,频率就会变).所以我想请教下大家.像我这种情况怎么把频率准确转换为信号的实际频率.
谢谢
回复
分享到:

使用道具 举报

 楼主| 发表于 2008-11-20 11:25 | 显示全部楼层
我到论坛上也找了好多帖子.但好象我没这样情况.
发表于 2008-11-20 12:02 | 显示全部楼层
你的意思是把坐标轴上的归一化频率转化成实际频率吗
 楼主| 发表于 2008-11-20 12:23 | 显示全部楼层
是啊
那样边际谱图出来就可以看见有哪些频率成分的信号
发表于 2008-11-20 22:46 | 显示全部楼层
在disp—hhs这个程序中把imagesc(t,fliplr([0,0.5]),im,[inf,0])改成imagesc(t,fliplr([0,0.5*fs]),im,[inf,0])就可以了
 楼主| 发表于 2008-11-21 11:21 | 显示全部楼层
回复棒棒糖,你是用什么程序求边际谱的?我上面那求边际谱没有那个.
你能否贴个程序看看?
发表于 2008-11-26 09:29 | 显示全部楼层


function disp_hhs(varargin)
error(nargchk(1,3,nargin));
fs = 0;
inf = -20;
im = varargin{1};
t = 1:size(im,2);
switch nargin
  case 1
    %raf
  case 2
    if isscalar(varargin{2})
      inf = varargin{2};
    else
      t = varargin{2};
    end
  case 3
    if isvector(varargin{2})
      t = varargin{2};
      inf = varargin{3};
    else
      inf = varargin{2};
      fs = varargin{3};
    end
  case 4
    t = varargin{2};
    inf = varargin{3};
    fs = varargin{4};
end
if isempty(inf)
  inf = -20;
end
if inf > 0
  inf = -inf;
elseif inf == 0
  error('inf must be nonzero')
end
M=max(max(im));
% warning off
im = 10*log10(im/M);
% warning on
figure
if fs == 0
  imagesc(t,[0,0.5],im,[inf,0]);
  ylabel('normalized frequency')
else
  imagesc(t/fs,[0,0.5*fs],im,[inf,0]);
  ylabel('frequency')
end
set(gca,'YDir','normal')
xlabel('time')
title('Hilbert-Huang spectrum')
 楼主| 发表于 2008-11-26 09:53 | 显示全部楼层
你这是求'Hilbert-Huang spectrum'的吧!
发表于 2013-5-24 16:42 | 显示全部楼层
li_1025tian 发表于 2008-11-26 09:53
你这是求'Hilbert-Huang spectrum'的吧!

请问有无求实际信号边际谱程序?我用lz的算算例准确,但算实际复杂信号就只能得到一条中间机会为零的线
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-5 20:31 , Processed in 0.065072 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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