| 查看: 4051 | 回复: 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 ] |
回复此楼» 猜你喜欢
求助一下有机合成大神
已经有4人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
请教限项目规定
已经有5人回复
-
最失望的一年
已经有16人回复
-
存款400万可以在学校里躺平吗
已经有33人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
浙江师范大学招收大气遥感、光学或人工智能方向的2026年物理学博士生1名
+2/278
-
专任教师招聘
+1/176
-
招收2026年春-申请考核博士生-仿生组装/生物能源方向
+1/174
-
海南大学专任教师招聘：海洋工程与海洋智能装备相关方向【长期有效】
+2/130
-
保障电池制造超纯环境：威格科技电池生产全线解决方案
+1/87
-
大连理工大学智能系统实验室优秀硕博研究生招生
+1/81
-
大连理工大学智能系统实验室优秀硕博研究生招生
+1/77
-
西交利物浦大学黄彪院士招收26年全奖博士生1名（工业智能方向）
+1/41
-
[长期合作招募] 同济大学肖倩老师团队诚邀港澳学者学术交流
+1/34
-
悉尼大学 AMME 机械工程 双ARC Future Fellows团队招收CSC博士生
+1/27
-
学生申博
+1/23
-
专科起点双非硕士，申博求建议
+1/21
-
中国科学院杭州医学研究所方晓红课题组-上海交通大学肖泽宇课题组联合博士后招聘启示
+2/20
-
沈阳征女友
+1/20
-
云南大学智能光电子课题组招收申请-考核博士研究生1名（光电、光学）
+1/16
-
南开大学齐迹课题组诚聘分子生物学、免疫学、有机分子合成相关方向的博士后、博士硕士
+1/14
-
中国科学院大连化学物理研究所-环境催化工程研究组（DNL 902组）事业编外项目聘用人员
+2/12
-
北京工业大学高靓教授课题组2026级博士研究生招生
+1/9
-
山东大学（青岛校区）招博士后（COF\MOF\催化\电池）
+1/3
-
中国农业科学院棉花研究所阴祖军研究员诚招植物生物学科研助理
+1/3
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