| 查看: 4221 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
求高手解释HHT变换求频率等问题
|
|||
|
利用HHT变换分解信号以及求瞬时频率等，但是为何不同IMF分量求出来的瞬时频率是一样的呢，搞不清楚。下面是matlab代码求频率的，高手看看是不是存在什么问题。 hx = hilbert(C); xr=real(hx);xi=imag(hx); sz=sqrt(xr.^2 + xi.^2); sx=atan(xi./xr); dt = diff(t); dx = diff(sx); sp = dx./dt; figure plot(t(1:N-1),sp) title(['InstantFreq imf ' num2str(i)]) end %计算HHT时频谱和边际谱 [A,infeq,tt] = hhspectrum(C,t,1,1); [im,tt1,Cenf] = toimage(A,infeq,tt,length(tt)); disp_hhs(im,tt1) %二维图显示HHT时频谱，E是求得的HHT谱 这里由两种瞬时频率，一种是根据HHT定义的根据相位时间导数求的（sp），还有下面的根据matlab内部函数instfreq求的瞬时频率（infeq）。sp有正有负，infeq严格 在0-0.5之间。不知道哪种方法是合理的，本人不是学这专业的，最近用到这个工具，但是有些搞不懂，请高手指点。 另外，上面还有一个求解中心频率的（Cenf），但是出来的结果是一个从0-0.5线性序列。而且出来只有一个序列，看EMD分解可看到有5个IMF分量。看到有文章说，中心频率就是每个IMF瞬时频率的平均值，如果是这样算的话倒是可以理解了。 下面三幅图分别是EMD分解的结果，HHT频谱图，还有个边际谱图。虽然做了但是看不懂太深的东西。请大家帮忙分析指点。  marginal_sp.png  EMD.png  HHT_spectrum.png [ Last edited by tony1230 on 2013-5-6 at 13:42 ] |
回复此楼» 猜你喜欢
析晶
已经有8人回复
-
本人42，博士刚毕业，现在找不到工作，怎么办？:(
已经有19人回复
-
河北省自然基金
已经有6人回复
-
博士申请
已经有5人回复
-
26级硕士毕业生求博导收留
已经有6人回复
-
有人投过CCC中国控制会议吗？
已经有3人回复
-
3,4-二羟基苯乙酮如何纯化？
已经有5人回复
-
国基评审
已经有10人回复
-
2026-博士申请
已经有4人回复
-
考研调剂
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
失业的第25天
+2/582
-
华侨大学全日制硕士博士（可不需要统考过了四级就行）-厦门校区-明天截止-急招
+1/137
-
华侨大学全日制硕士博士（特殊情况可在职读）-厦门校区-名额急招-4月底截止
+1/82
-
中国水产科学研究院黄海水产研究所张旭志课题组招聘 可穿戴传感器研究 科研助理
+1/80
-
二区期刊Molecules征稿，欢迎投稿
+1/79
-
上海交通大学招聘催化反应工程/能源过程系统优化控制方向博士后
+1/30
-
大湾区大学张国强课题组招聘光催化方向博士后
+1/30
-
华中科技大学杨辉教授团队招收： 科研助理、 博士后、2027博士/硕士研究生
+1/29
-
中山大学院士团队微纳器件、脑机接口方向博士招生、博后招聘
+1/23
-
美国博士后招聘，代谢的转录调控方向，干+湿结合实验室，欢迎Email咨询。
+1/15
-
爱尔兰 - 可持续金融+航空方向 全奖博士 (2个名额)
+1/12
-
湖南大学张旺课题组硕士/博士/博士后招聘
+1/11
-
明尼苏达大学博士后招聘（结构病毒生物学, Prof. Fang Li's lab)
+1/8
-
深圳理工大学梁国进课题组招收液流电池界面电化学方向博士生
+1/6
-
广东工业大学机电工程学院高伟教授招收2026级博士生
+1/4
-
急聘！中国农科院兰州牧药所兽药创新与耐药性科研团队聘科研助理
+1/3
-
暨南大学招超表面方向博士后
+1/3
-
中国科学院青岛生物能源与过程研究所科研助理/工程师招聘启事
+1/3
-
大连理工大学张硕课题组 2026 年秋季博士生招生启事（有机合成/糖化学方向）
+1/2
-
药学人员在疾控中心到底能不能晋升职称
+1/1
7楼2014-02-28 13:50:24
2楼2013-05-06 14:24:54
祝一切顺利，科研成果累累！顺便抢金币！3楼2013-05-06 18:12:42
4楼2013-05-06 21:20:24