声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1257|回复: 1

[HHT] HHT的一点问题,请大侠解答

[复制链接]
发表于 2009-6-5 23:15 | 显示全部楼层 |阅读模式

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

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

x
小弟最近做毕设。matlab和HHT都不是很熟悉。遇到了很多问题,看论坛上的贴子,得到了很多启发。但是还有一些问题。希望大侠帮我看看
我用的程序是法国人的。
下面是我自己编的一段主程序,用来求HHT谱。但是得到的谱图,横坐标是数据点,而不是时间。纵坐标是归一化的频率,而不是实际频率。如图
1.jpg
下面是我的程序:
imf=emd(x);
a=length(imf);
imf=cell2mat(imf);%我的matlab2009得到的imf集是cell格式的,这里先转换为矩阵。
xx=reshape(imf,7000,a);
xx=xx';
[A,f,tt] = hhspectrum(xx,1:7000,2,1);
[im,tt,ff] = toimage(A,f);
disp_hhs(im)

因为我对hhspectrum   toimage disp_hhs这三个程序的一些参数不是很明白,所以老是调试不对。请大家帮我看看我这段程序对不对,应该怎样改才能得到实际频率和时间。
另外,我用破凰的程序求边际谱得到的结果是这种的:
2.jpg
在x=0处总有一个很大的振幅。
fs=350;    %fs为采样频率;
N=7000;     %采样点数
[A,fa,tt]=hhspectrum(xx);
[E,tt1]=toimage(A,fa,tt,length(tt));
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('幅值');
这是程序。不知道是我这程序错了,还是前面的EMD分解就错了。请大侠解答,不胜感激。

x.mat

24.19 KB, 下载次数: 3

回复
分享到:

使用道具 举报

 楼主| 发表于 2009-6-9 21:49 | 显示全部楼层

回复 楼主 lfine100 的帖子

没人解答么
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-13 19:56 , Processed in 0.066894 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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