声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 6399|回复: 27

[HHT] 怀疑instfreq.m程序的正确性

[复制链接]
发表于 2007-3-20 21:02 | 显示全部楼层 |阅读模式

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

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

x
用instfreq.m求瞬时频率发现第二个imf的瞬时频率比第一个imf大
谁用过这个程序吗?
有问题吗?
回复
分享到:

使用道具 举报

发表于 2007-3-20 21:13 | 显示全部楼层
原帖由 count_1126 于 2007-3-20 21:02 发表
用instfreq.m求瞬时频率发现第二个imf的瞬时频率比第一个imf大
谁用过这个程序吗?
有问题吗?



正常现象,强烈建议你阅读本版关于EMD的帖子!!!!!特别是这个帖子:emd方法的 几点不明 请高手予以解答  7楼
发表于 2007-3-20 21:49 | 显示全部楼层
eight你的意见是什么呢?关于我的问题?
发表于 2007-3-20 21:56 | 显示全部楼层
原帖由 zhangnan3509 于 2007-3-20 21:49 发表
eight你的意见是什么呢?关于我的问题?



我是回答楼主的问题,指引他阅读那个帖子7楼的以下内容:

1. c1 c2 ... cn 的确是严格按照频率从高到低产生的,不过这里有一个误区,其意思并不是说c1的频率一定比c2的高,正确的理解是c1中的某个局部的频率比c2中相同局部的频率要高,这也正好反映了EMD算法局部性强的本质所在,也跟黄的说法“相邻的分量可能包含相同时间尺度的振荡,但是相同时间尺度的振荡绝对不会出现在两个不同的IMF分量的同一个位置”一致
 楼主| 发表于 2007-3-20 22:03 | 显示全部楼层

读过了 不过还是疑问

Himf(1,:)=hilbert(imf(1,:));
[fnormhat1,t1]=instfreq(Himf(1,:)');         

Himf(2,:)=hilbert(imf(2,:));
[fnormhat2,t2]=instfreq(Himf(2,:)');


plot(t1,fnormhat1,'r',t2,fnormhat2,'b');

对一个简单信号:x=cos(t)+0.5*cos(t/2);
红色的整个在蓝色图形下面.
发表于 2007-3-20 22:10 | 显示全部楼层
我是说以前的emd程序和后来改过的emd相比 有一些变化,比如后来的有fixe 和 fixe_h 你觉得这两个中止条件很有必要吗?起到了什么作用?
发表于 2007-3-20 23:27 | 显示全部楼层
原帖由 zhangnan3509 于 2007-3-20 22:10 发表
我是说以前的emd程序和后来改过的emd相比 有一些变化,比如后来的有fixe 和 fixe_h 你觉得这两个中止条件很有必要吗?起到了什么作用?



这个很少留意,呵呵。不过貌似 flandrin 在 emd.m 中实现很多功能,其中之一是用一个 mask 信号来区分频率相近的信号(可参考:The Use of a Masking Signal to Improve Empirical Mode Decomposition),还有就是你提到的 fixe 和 fixe_h,这两个均代表停止准则的不同实现,粗略看了一下, fixe 应该是不断迭代直到信号单调(极值点个数小于等于2),而 fixe_h 应该是只检验IMF的第一个条件(过零点条件),忽略零均值条件
发表于 2007-3-21 08:50 | 显示全部楼层
多加的这两个fixe和fixe_h在处理过程中作用很明显吗?mask好像是用来消除模态混叠的
发表于 2007-3-21 09:19 | 显示全部楼层
现在的问题是c1的瞬时频率的值全部小于c2,不知道是为什么?
发表于 2007-3-21 09:32 | 显示全部楼层
不是吧 是局部对应部分吗?你说的太泛泛了,大家都不知道具体情况是什么!
发表于 2007-3-21 09:42 | 显示全部楼层
5楼的程序对吗?就是用那个程序画出的图,从c1的瞬时频率全部小于c2的瞬时频率,就是在每一个点处,c1都小于c2。
 楼主| 发表于 2007-3-21 10:13 | 显示全部楼层

各位高手说说呀

确实有5楼说的问题
急死了!
还有那个hhspectrum.m
好像也不能用啊
发表于 2007-3-21 10:30 | 显示全部楼层
原帖由 count_1126 于 2007-3-21 10:13 发表
确实有5楼说的问题
急死了!
还有那个hhspectrum.m
好像也不能用啊



你以为只有你一个人急吗?我也很急,不过不是急你的问题,是急自己的问题,不要以为个个都很有空,只有你很忙。遇到问题不会自己先动手解决吗?我n多问题都是自己解决的,不要养成太依赖别人的习惯,尤其是男孩子
发表于 2007-3-21 10:32 | 显示全部楼层
多加的这两个fixe和fixe_h在处理过程中作用很明显吗?


没用过,可以试试


mask好像是用来消除模态混叠的


是的,不过效果也不太好
发表于 2007-3-21 10:33 | 显示全部楼层
原帖由 xxf 于 2007-3-21 09:19 发表
现在的问题是c1的瞬时频率的值全部小于c2,不知道是为什么?



没有试过,最近也忙,迟点有时间再看看,或者等其他高手解答
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-25 02:08 , Processed in 0.065842 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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