dushudushu 发表于 2010-12-15 22:31

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

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

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

dushudushu 发表于 2010-12-15 22:37

怎么没有上传附件的选项呢

dushudushu 发表于 2010-12-15 22:40

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

数据

hitdely 发表于 2010-12-16 09:15

利用希尔伯特黄变换就行啊 还有一个编号的程序instfreq.m 都可以求解频率

hahaer 发表于 2010-12-17 22:47

我感觉楼主问的是时-频图，也就是希尔伯特谱图吧
这个例子论坛里有，搜索置顶的帖子

hivdg 发表于 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时频图');
页: [1]
查看完整版本: 怎么求对应IMF分量的频率呢