| 查看: 4080 | 回复: 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 ] |
回复此楼» 猜你喜欢
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有11人回复
-
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了
已经有9人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有10人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有10人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有11人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
医学超声影像负责人招聘-中国科学院赣江创新研究院
+1/980
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘
+2/256
-
985教授征女友
+1/252
-
苏州国家实验室和中国科学技术大学联培博士招生
+1/188
-
华南师范大学（211）- 光电科学与工程学院 - 申请审核制（2026年4-5月份面试考核）
+2/106
-
澳洲西澳大学Dr Yiran Liu招全额奖学金和CSC奖学金博士生（3.8万澳币/年）
+1/105
-
石河子大学化学化工学院CJ学者领衔分子炼油团队招收博士、硕士。【接收调剂】
+1/89
-
江西理工大学联合中国科学院赣江创新研究院招收2026级博士研究生
+1/82
-
欢迎报考南京农业大学植物环境适应课题组课题组2026级博士生。
+1/75
-
好玩的不敢搞，能搞的不挣钱，能挣钱的我不会做
+1/67
-
感谢小木虫的缘分
+1/38
-
罗格斯大学纽瓦克校区(Rutgers-Newark) 招收 PHD，计算材料物理方向
+1/37
-
香港中文大学医学院 诚聘 研究助理教授 (医工结合/生物信息学方向)
+1/36
-
澳大利亚麦考瑞大学(Macquarie University)国际博士硕士全额奖学金-计算机-26年中开学
+1/34
-
华南师范大学（211）- 光电科学与工程学院 - 申请审核制（2026年4-5月份面试考核）
+2/18
-
【青岛大学】2026年生物与医药申请考核制博士生招生（含少数民族骨干人才）
+1/18
-
都放假了嘛？
+1/7
-
广东省环境科学研究院招聘高分辨质谱方向博士一名
+1/7
-
【东南大学博士后、科研助理招聘】
+1/5
-
美国密苏里大学“柔性电子”课题组诚招博士研究生
+1/1
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