| 查看: 4278 | 回复: 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 ] |
回复此楼» 猜你喜欢
宿州学院学报
已经有5人回复
-
青B发送上会通知了吗
已经有14人回复
-
博士申请
已经有5人回复
-
4,4二甲基联苯干啥用，有懂得吗
已经有4人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有7人回复
-
论文撤稿了
已经有9人回复
-
化学专业申博
已经有5人回复
-
招收2026级博士生
已经有5人回复
-
医学类期刊求推荐
已经有6人回复
-
26/27申博自荐
已经有10人回复
» 抢金币啦！回帖就可以得到:
-
【专业EPR】：自由基(超氧、羟基、单线态氧)，空位缺陷；定量计算，VX: 761711562。
+1/90
-
【通知】北京信息科技大学仪器科学与光电工程学院招收博士研究生（2026）
+1/85
-
武汉工程大学（省属一本）招聘师资博士后以及人才引进教师（事业编）
+1/78
-
福州大学化工学院电子化学品团队博士招生，还有一个名额！
+1/75
-
大连工业大学杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
+1/34
-
【通知】北京信息科技大学仪器科学与光电工程学院招收博士研究生（2026）
+2/30
-
招聘青年教师（有编制）——南京邮电大学柔性电子全国重点实验室徐申课题组
+1/21
-
盤後看下來，資金沒退，也算是吃上了。
+1/10
-
西安交通大学陕西省科技创新团队招聘2027年免试保送生和考研生，2-4个名额
+1/9
-
北理工国家杰青团队招博士后
+1/9
-
深圳大学物理与光电工程学院/深圳先进光源研究院招收2026级联合培养博士生
+1/7
-
电子科技大学材料学院SFT创新中心招收准备考硕和读博的科研助理 理工医交叉方向
+1/6
-
张兆威教授课题组长期招聘科研助理
+1/5
-
易度质量流量计在燃气阀流通性测试中的应用
+1/5
-
地大武汉（211）招生光刻胶半导体方向博士生
+1/4
-
上海交通大学海优建组招开山弟子
+1/2
-
四川大学招聘联合培养学生/科研助理/2027届博士
+1/2
-
北航杭州国际校区招聘3D 打印，陶瓷，力学等方向博士后
+1/1
-
【有偿访谈招募】高才通来港后，你过得还好吗？
+1/1
-
【有偿访谈招募】高才通来港后，你过得还好吗？
+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