声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1819|回复: 5

[HHT] 怎么求对应IMF分量的频率呢

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

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

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

x
本帖最后由 dushudushu 于 2010-12-15 22:36 编辑

对一个信号用EMD分解后,得出了IMF分量,怎么作出每个分量的时频图呢,横坐标为时间,纵坐标为频率的,工具箱里有这个函数吗,新手,牛人帮忙看下啊,我把数据放到附件里,采样间隔为0.005s

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

 楼主| 发表于 2010-12-15 22:37 | 显示全部楼层
怎么没有上传附件的选项呢
 楼主| 发表于 2010-12-15 22:40 | 显示全部楼层
本帖最后由 dushudushu 于 2010-12-15 22:41 编辑

数据
发表于 2010-12-16 09:15 | 显示全部楼层
利用希尔伯特黄变换就行啊 还有一个编号的程序instfreq.m 都可以求解频率
发表于 2010-12-17 22:47 | 显示全部楼层
我感觉楼主问的是时-频图,也就是希尔伯特谱图吧
这个例子论坛里有,搜索置顶的帖子
发表于 2010-12-23 15:50 | 显示全部楼层
imf1=emd(y);
imf(i)=imf1(i,:);%你要求的第i个的IMF分量的时频谱
l=1;%默认时间截断头尾各一个点,共截断2个点
ttt=1:length(imf);
tt=ttt((l+1):(length(imf)-1));%截断头尾2个时间点后的截断时间,对应Hilber谱的时间序列
an=hilbert(imf')';
f=instfreq(an',tt,l)';
A=abs(an((l+1):(length(imf)-1)));
plot(tt/fs,f*fs);%fs是你的采样率
ylabel('瞬时频率');xlabel('时间');title('IMF时频图');
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 13:44 , Processed in 0.077396 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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