对一实测信号的处理，希望大家参与讨论

tamujin

回复 1 # bcyangbc 的帖子

楼主你好，刚才那个问题解决了，Hilbert-Huang时频谱画出来了。
我把程序里的
disp_hhs(im,tt,[],sf)
改为
disp_hhs(im)
就不报错了！

lixiaolong08

楼主太厉害了，下下来好后学习下

hustlikaifu

我也下载了学习

danielrichie

谢谢楼主分享！

one23for

受益匪浅，谢谢

若凡

[A,fa,tt]=hhspectrum(IMF(1:nIMF,:));    Hilbert-Huang变换，计算瞬时频率和瞬时幅值
运行了程序，出现一个错误
??? Output argument "A" (and maybe others) not assigned during call to "package_emd\EMDs\hhspectrum.m (hhspectrum)".
Error in ==> hhspectrum at 20
if nargin < 2
请教下是什么原因？！！！
function [A,f,tt] = hhspectrum(imf,t,l,aff)%其中hhspectrum是从网上下载的，并未改动
if nargin < 2
  t=1:size(imf,2);
end
if nargin < 3
  l=1;
end
if nargin < 4
  aff = 0;
end
lt=length(t);
tt=t((l+1):(lt-l));
for i=1:(size(imf,1)-1)
  an(i,:)=hilbert(imf(i,:)')';
  f(i,:)=instfreq(an(i,:)',tt,l)';
  A=abs(an(:,l+1:end-l));
  if aff
    disp(['mode ',int2str(i),' trait&#65533;'])
  end
end

zhaowei_huang

谢谢高手们的无私奉献

yoooog

感谢楼主分享，学习下~

白云飞

感谢楼主无私奉献

sjdwjt

楼主好厉害，向你学习！我也在用EMD分析地震动信号，只是效果不太理想啊

lovme

收藏了，虽然不懂，但楼主大牛！贡献很大！

vampire101

学习了，不错啊

lixuanya

顶！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！11

yghit08

回复 91 # tamujin 的帖子

原来的那个disp_hhs里面有一个bug，修改后就可以了。
直接用disp_hhs(im)得出的结果是归一化的，disp_hhs(im,tt,[],sf)得到的实际的频率。
关于disp_hhs的改动就是将error(nargchk(1,3,nargin))改成error(nargchk(1,4,nargin));
这也是在论坛上得到的，忘了谁，感谢他

yghit08

回复 90 # tamujin 的帖子

在后面给你解决了，请参考