声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 937|回复: 3

[HHT] 求边际谱的问题,频率不对应

[复制链接]
发表于 2008-9-25 09:48 | 显示全部楼层 |阅读模式

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

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

x
fs=1000;    %fs为采样频率;
N=1000;     %
采样点数
t=1/fs:1/fs:1;
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*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('
幅值');
求出的边际谱如下图,为什么得到的频率不是30和45呢?
untitled.jpg
回复
分享到:

使用道具 举报

发表于 2008-9-25 11:24 | 显示全部楼层
HHT中频率是指瞬时频率,我想问的是f=(0:N-3)/N*(fs/2)是指瞬时频率吗?好像不对吧,如果这样画边际谱肯定不对的。
发表于 2008-9-26 15:22 | 显示全部楼层

回复 楼主 tjucruiser 的帖子

E=flipud(E);
把这个语句去掉,新版的程序不需要这个!
 楼主| 发表于 2008-9-26 15:27 | 显示全部楼层
多谢您了,我刚试了,果然是这样。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 20:11 , Processed in 0.066019 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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