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

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 编辑

=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 = 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

{:{23}:}楼主好厉害，向你学习！我也在用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 的帖子

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

查看完整版本: 对一实测信号的处理，希望大家参与讨论