声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1074|回复: 3

[HHT] 边际谱问题求助,请帮忙一下看哪里出问题了(急)

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

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

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

x
我是用法国人的hhspectrum和toimage做的,自己编写的边际谱程序,
当我用如下程序计算式
function bjp(x,Fs)
N=length(x);
fs=Fs;
imf=emd_online(x);
[A,f,tt]=hhspectrum(imf);
[im,tt] = toimage(A,f);
im=flipud(im);
for k=1:size(im,1)
    bjp(k)=sum(im(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
title('Hilbert边际谱')

出现如f和bjp长度不等错误。


我每个变量计算后的确长度不等。
改为如下程序
function bjp(x,Fs)
N=length(x);
fs=Fs;
imf=emd_online(x);
[A,f,tt]=hhspectrum(imf);
[im,tt] = toimage(A,f);
%im=flipud(im);
for k=1:size(im,2)
    bjp(k)=sum(im(:,k))*1/fs;

end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
title('Hilbert边际谱')


但计算出如下的边际谱

边际谱

边际谱
回复
分享到:

使用道具 举报

 楼主| 发表于 2009-9-22 15:13 | 显示全部楼层
并且计算出的Hilbert谱为如下,请高手看看是否正确,谢谢。附计算数据
Hilbert谱.jpg

shuju.mat

120.68 KB, 下载次数: 0

 楼主| 发表于 2009-9-22 15:18 | 显示全部楼层
附清晰地边际谱图
边际谱.jpg
 楼主| 发表于 2009-9-22 15:22 | 显示全部楼层
附程序,谢谢

emd_online.m

25.7 KB, 下载次数: 11

disp_hhs.m

1.92 KB, 下载次数: 6

hhspectrum.m

1.41 KB, 下载次数: 6

toimage.m

2.93 KB, 下载次数: 8

bjp.m

271 Bytes, 下载次数: 7

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-13 19:49 , Processed in 0.069742 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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