声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1726|回复: 7

[HHT] 一个有关边际谱的问题

[复制链接]
发表于 2008-7-26 16:15 | 显示全部楼层 |阅读模式

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

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

x
前几天再论坛上看到一个有关边际谱的程序,有的语句不太懂,请高手指点一下
clear;
fs=1000;    %fs
为采样频率;
tspan=2;
t=1/fs:1/fs:tspan;
N=length(t);     %采样点数
y1=5*sin(2*pi*241*t);
y2=3*sin(2*pi*73*t);
y=[y1;y2;zeros(size(y1))]; %IMF

%%%%%%%%%%%%%
求边际谱
[A,fa,tt]=hhspectrum(y);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E)
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('
幅值');
[E,tt1]=toimage(A,fa,tt,length(tt));这句是什么意思,其中的E是不是指toimage.m中的im, E(k,:)又是什么,为什么
边际谱就是bjp(k)=sum(E(k,:))*1/fs
回复
分享到:

使用道具 举报

发表于 2008-10-15 13:41 | 显示全部楼层
同样的问题。E指的是什么?toimage。m的输出没有E啊?
E是指im还是ff呢?不是很明白。
边际谱是对hilbert谱的时间积分,程序怎么体现出来的?
不明白。
望高手指点一下
发表于 2008-10-19 22:42 | 显示全部楼层

自己给自己解答!呵呵。。

 楼主| 发表于 2008-10-23 11:09 | 显示全部楼层
我觉得边际谱就是bjp(k)=sum(E(k,:)) 指某一频率对应的振幅总和,而不是bjp(k)=sum(E(k,:))*1/fs ; E就是im
另外频率是f0=(0:size(im,1)-1)*fs/size(im,1)+1/16;
法国人的程序中im是用了一个很好的技术得出来的,可以看子程序 toimage和DISP_HHS:@D

[ 本帖最后由 wang-wei 于 2008-10-23 11:12 编辑 ]
发表于 2008-10-23 16:39 | 显示全部楼层

时频工具箱

振动论坛里面有的。你再找找。
发表于 2008-10-23 16:43 | 显示全部楼层
原帖由 wang-wei 于 2008-10-23 11:09 发表
我觉得边际谱就是bjp(k)=sum(E(k,:)) 指某一频率对应的振幅总和,而不是bjp(k)=sum(E(k,:))*1/fs ; E就是im
另外频率是f0=(0:size(im,1)-1)*fs/size(im,1)+1/16;
法国人的程序中im是用了一个很好的技术得出来的, ...

可以把你说的法国的那个im传上来么?
另外而不是bjp(k)=sum(E(k,:))*1/fs 是要积分用的啊。
 楼主| 发表于 2008-10-25 14:20 | 显示全部楼层
im和这的E是一样的,不过表示不同
[E,tt1]=toimage(A,fa,tt,length(tt));即[im,tt1]=toimage(A,fa,tt,length(tt));
发表于 2014-6-27 06:43 | 显示全部楼层
谢谢指点!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-13 13:01 , Processed in 0.064300 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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