| 查看: 4426 | 回复: 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人回复
-
b口会评
已经有3人回复
-
2026年面上项目中了，2A+B, 会评顺利通过
已经有13人回复
-
生命口会评
已经有4人回复
-
大龄残疾硕士的一点执念
已经有21人回复
-
今年E04面上
已经有20人回复
-
网络举报学术不端成为了男女之间情杀的手段了，真荒唐。。。
已经有3人回复
-
Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联，生成B2Pin2杂质
已经有3人回复
-
刑法学论文投稿求助
已经有5人回复
-
Journal of Environmental Chemical Engineering
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
filecode有几个（（JTJC））
+1/563
-
27年入学博士招生（图像智能处理相关）
+1/180
-
热分析测试
+1/85
-
固态电池专家
+1/73
-
在深圳的河南人，男找女
+1/68
-
工材学部会评时间知道么？
+1/60
-
浙江大学衢州研究院肖成梁课题组招聘博士后及科研助理
+1/34
-
上海交通大学化工学院邱惠斌教授课题组招募科研助理、交流生
+1/34
-
厦门大学张桥保教授团队招收2027年硕士/博士研究生以及博士后（材料、化学方向）
+1/31
-
回首往事
+1/30
-
沈阳农业大学招聘微生物酶工程方向博士！！！！！！！！！！！！！
+1/29
-
华南理工大学电力学院雪映教授招收博士生（2027年秋季入学）
+2/26
-
暨南大学开明峰教授团队博士后招聘启事
+1/13
-
Urgent！知名外资仪器厂家急招Application Scientist（2名，售前，SH）
+1/9
-
2氯5甲基吡啶和2氯6三氯甲基吡啶大生产工艺需求
+1/6
-
圣路易斯华盛顿大学管建均教授课题组招聘博士后（2名）
+1/5
-
力学研究所某专项团队支撑岗位和特别研究助理岗位招聘启事
+1/5
-
清华大学韩军功教授课题组博士后招聘启事
+1/3
-
力学研究所某专项团队支撑岗位和特别研究助理岗位招聘启事(长期有效)
+1/2
-
科研测试平台求助
+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