| 查看: 4239 | 回复: 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 ] |
回复此楼» 猜你喜欢
国自科送审了吗
已经有10人回复
-
博士招生
已经有5人回复
-
夜，静悄悄的
已经有8人回复
-
研究生做的很差，你们会让毕业吗？
已经有6人回复
-
26年博士申请自荐-电催化
已经有5人回复
-
2026博士或科研助理转27年博士
已经有5人回复
-
26年申博自荐-计算机视觉
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
湖南大学化工院招收2026年入学博士
+1/177
-
东南大学废弃碳资源利用团队（杰青/海优团队）博士后、科研助理及27级博士生招募
+1/128
-
东南大学废弃碳资源利用团队（杰青/海优团队）27级博士研究生和科研助理招生
+1/127
-
2026年江西师范大学人文地理与低碳经济方向拟招收1名“申请入学制”学术博士研究生
+2/108
-
生命学部的面地青本子开始分发了么？
+1/61
-
松山湖材料实验室-大连理工大学联合招收博士研究生
+1/35
-
北京航空航天大学拟招收2026博士研究生
+1/35
-
Google ai pro 会员分享
+1/28
-
智合健物课题组2026年博士生招生(5月6日—2026年5月13日)
+1/22
-
深圳大学材料学院“智能分子材料团队”招聘博士后（连续流微反应器方向优先）
+1/17
-
2026年西安交通大学博士招生！
+1/16
-
南京林业大学-国家级青年人才团队 招2026级博士 （合成化学）
+1/10
-
吉林大学任雷课题组招机器人和触觉传感方向2026年入学博士生
+1/9
-
注册小木虫周年大礼包领取不了
+1/5
-
招收2026年秋季入学博士生1名（河北工业大学/北京科技大学联合 增材制造/生物材料）
+1/4
-
港科大广州 - 量子计算全额奖学金博士招聘（申请制，2026和2027均有名额）
+1/4
-
类器官/器官芯片 华西医院有编制研究员招聘
+1/4
-
天津医科大学基础医学院博士后招聘
+1/4
-
香港科技大学 招生 全奖博士 -- 机器人/电子/材料
+1/3
-
南京大学医学院医学数字孪生方向博士后招聘
+1/2
5楼2013-12-29 14:19:21
6楼2014-02-28 13:49:04
2楼2013-05-06 14:24:54
祝一切顺利，科研成果累累！顺便抢金币！3楼2013-05-06 18:12:42
4楼2013-05-06 21:20:24
7楼2014-02-28 13:50:24
8楼2015-01-18 14:40:24
9楼2016-02-14 18:56:51
|
本帖内容被屏蔽 |
10楼2017-06-02 21:48:43
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
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 发自小木虫Android客户端 |
11楼2017-06-04 11:09:57