| 查看: 4920 | 回复: 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 ] |
回复此楼» 猜你喜欢
» 抢金币啦！回帖就可以得到:
华南师范大学招收2027级生物学博士研究生1名，研究方向：神经生物学
+5/1005
-
诚征女友
+1/460
-
电子科技大学电子科学与工程学院向全军教授2027年博士研究生招生
+2/452
-
云南大学丁军桥团队2027年博士招生
+1/88
-
山东征女友，坐标济南
+1/75
-
昆明 征女友
+1/74
-
香港城市大学（东莞）雷进程课题组博士研究生、研究助理招聘
+1/42
-
湘潭大学材料科学与工程学院新型微电子材料与器件课题组招收博士，2027年秋入学
+1/39
-
2027年上海大学王亮课题组招收电催化研究背景的博士生
+2/38
-
华南师范大学（211）- 光电科学与工程学院招收物理学（一流学科）2027级博士
+1/34
-
【招生】国家级人才王新才教授课题组招收2027级“申请-考核”制博士研究生
+1/30
-
中国科学院北京纳米能源与系统研究所杨亚课题组招收2027年凝聚态物理硕士和博士
+1/22
-
湖南师范大学化学化工学院-熊二虎教授课题组招聘青年教师和博士后（长期有效）
+1/19
-
中山大学柔性电子方向，博士招生/博士后招聘（材料、材料加工、化学、化工等背景）
+1/19
-
香港科技大学（广州）段圭粤课题组招收风能与流体力学方向博士生、研究助理、博士后
+1/14
-
南方科技大学田雷蕾课题组招收2027年博士生
+1/13
-
配制溶液时如果其他都正确，定容时溶液温度为30度，溶液浓度偏大还是偏小
+1/13
-
【2027级博士申请】双非硕・一作 SCI 3 篇（2篇CEJ 一区 Top）储能方向
+1/3
-
gpt官方代充
+1/3
-
东北大学招收2027级硕士、博士研究生（推免、报考）凝聚态物理/微电子/电子等相关方向
+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