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

tamujin · 发表于 2012-2-27 15:12

回复 1 # bcyangbc 的帖子

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

lixiaolong08 · 发表于 2012-2-29 16:59

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

hustlikaifu · 发表于 2012-3-4 18:41

我也下载了学习

danielrichie · 发表于 2012-3-12 15:00

谢谢楼主分享！

one23for · 发表于 2012-3-14 16:57

受益匪浅，谢谢

若凡 · 发表于 2012-3-16 10:53

本帖最后由若凡于 2012-3-16 10:56 编辑

[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�'])
  end
end

zhaowei_huang · 发表于 2012-3-18 23:24

谢谢高手们的无私奉献

yoooog · 发表于 2012-3-20 11:48

感谢楼主分享，学习下~

白云飞 · 发表于 2012-3-28 17:34

感谢楼主无私奉献

sjdwjt · 发表于 2012-3-29 13:37

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

lovme · 发表于 2012-3-29 16:27

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

vampire101 · 发表于 2012-4-2 20:56

学习了，不错啊

lixuanya · 发表于 2012-4-8 20:48

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

yghit08 · 发表于 2012-4-13 22:08

回复 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 · 发表于 2012-4-13 22:12

回复 90 # tamujin 的帖子

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

账号		自动登录	找回密码
密码			我要加入

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