| 查看: 4182 | 回复: 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 ] |
回复此楼» 猜你喜欢
材料与化工371求调剂
已经有9人回复
-
085600调剂
已经有9人回复
-
346分的生物与医药08600求调剂
已经有6人回复
-
304求调剂（085602，过四级，一志愿985）
已经有13人回复
-
328分调剂
已经有6人回复
-
求调剂
已经有9人回复
-
358求调剂
已经有7人回复
-
生物与医药调剂
已经有3人回复
-
367求调剂
已经有3人回复
-
工科08专硕机械275求调剂
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
武汉轻工大学化学与环境工程学院廖丽霞副教授招收硕士研究生调剂信息
+1/482
-
2026年浙江农林大学园艺学院硕士招收，接受硕士调剂学生
+1/190
-
铟锡化合物半导体材料与器件团队[推介]
+3/154
-
德州学院材料与化工专硕招收大量调剂生
+1/95
-
济南大学前沿交叉科学研究院招收2026级调剂生（化学工程085602）
+1/86
-
东莞理工土木工程和土木水利专业硕士接受调剂，名额充足！
+1/85
-
邵阳学院食品与化学工程学院 生物与医药专业接受调剂生
+2/54
-
武汉纺织大学学硕调剂
+2/52
-
2026年宁波工程学院杜俊毅副研究员招收硕士研究生（含调剂）
+1/49
-
重庆大学李延军课题组招聘有机方向科研助理
+1/27
-
北京理工大学（珠海）招收2026级考核制博士
+1/18
-
招硕士研究生！招硕士研究生！招硕士研究生！招生专业：化学
+1/13
-
中科大科学岛分院2026年工程博士申请考核制-博士招生，4月10日截止
+1/11
-
首都师范大学化学系田洋课题组专硕招调剂生
+1/10
-
申博的看过来，湖南农大2026年农药学 紧急 补招一名博士研究生，时间紧迫
+1/9
-
长沙理工大学热环境方向招收专硕调剂生
+1/9
-
徐工-环境工程学院-招收调剂硕士
+1/8
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/7
-
浙江宁波 生物与医药 接受生物学、园艺学、作物学、农业调剂
+1/5
-
福建理工大学招收专业代码08开头学术型硕士研究生
+1/5
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
