tony1230

銀蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 511.9
帖子: 72
在線: 73.3小時(shí)
蟲號(hào): 1793983

[交流] 求高手解釋HHT變換求頻率等問(wèn)題

利用HHT變換分解信號(hào)以及求瞬時(shí)頻率等，但是為何不同IMF分量求出來(lái)的瞬時(shí)頻率是一樣的呢，搞不清楚。下面是matlab代碼求頻率的，高手看看是不是存在什么問(wèn)題。
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

%計(jì)算HHT時(shí)頻譜和邊際譜
[A,infeq,tt] = hhspectrum(C,t,1,1);
[im,tt1,Cenf] = toimage(A,infeq,tt,length(tt));
disp_hhs(im,tt1) %二維圖顯示HHT時(shí)頻譜，E是求得的HHT譜

這里由兩種瞬時(shí)頻率，一種是根據(jù)HHT定義的根據(jù)相位時(shí)間導(dǎo)數(shù)求的（sp），還有下面的根據(jù)matlab內(nèi)部函數(shù)instfreq求的瞬時(shí)頻率（infeq）。sp有正有負(fù)，infeq嚴(yán)格在0-0.5之間。不知道哪種方法是合理的，本人不是學(xué)這專業(yè)的，最近用到這個(gè)工具，但是有些搞不懂，請(qǐng)高手指點(diǎn)。

另外，上面還有一個(gè)求解中心頻率的（Cenf），但是出來(lái)的結(jié)果是一個(gè)從0-0.5線性序列。而且出來(lái)只有一個(gè)序列，看EMD分解可看到有5個(gè)IMF分量。看到有文章說(shuō)，中心頻率就是每個(gè)IMF瞬時(shí)頻率的平均值，如果是這樣算的話倒是可以理解了。

下面三幅圖分別是EMD分解的結(jié)果，HHT頻譜圖，還有個(gè)邊際譜圖。雖然做了但是看不懂太深的東西。請(qǐng)大家?guī)兔Ψ治鲋更c(diǎn)。

marginal_sp.png

EMD.png

HHT_spectrum.png

[ Last edited by tony1230 on 2013-5-6 at 13:42 ]

回復(fù)此樓

» 猜你喜歡

招收調(diào)劑碩士已經(jīng)有12人回復(fù)
296求調(diào)劑已經(jīng)有8人回復(fù)
招收博士1-2人已經(jīng)有3人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有11人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有9人回復(fù)
274求調(diào)劑已經(jīng)有8人回復(fù)
290求調(diào)劑已經(jīng)有4人回復(fù)
一志愿西南交通專碩材料355 本科雙非求調(diào)劑已經(jīng)有4人回復(fù)
295復(fù)試調(diào)劑已經(jīng)有6人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有8人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2013-05-06 13:09:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

niuyun1217

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46.5
帖子: 1
在線: 1.2小時(shí)
蟲號(hào): 2767947

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

你好，請(qǐng)問(wèn)HHT變換求頻率你會(huì)了嗎？可以把程序發(fā)給我嗎？謝謝，急用

贊一下(2人)

回復(fù)此樓

5樓2013-12-29 14:19:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

rock_yang

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 619.3
帖子: 49
在線: 17.8小時(shí)
蟲號(hào): 1348710

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

關(guān)于頻率的部分的話，把計(jì)算得到的頻率乘以采樣頻率就是實(shí)際頻率了，至于為什么是0.5，你應(yīng)該懂的

hilbert變換是對(duì)窄帶信號(hào)才有效的，所以才有了EMD分解將任意信號(hào)分解成多個(gè)窄帶信號(hào)

所以如果對(duì)一個(gè)窄帶信號(hào)做hilbert變換的話，去掉負(fù)頻率，double正頻率就是hilbert譜了

贊一下(2人)

回復(fù)此樓

6樓2014-02-28 13:49:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

tony1230

銀蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 511.9
帖子: 72
在線: 73.3小時(shí)
蟲號(hào): 1793983

我把上面的頻譜濾掉了最高頻的兩個(gè)信號(hào)，依次得到下面兩個(gè)頻譜圖，但是高頻部分為何還留有比較強(qiáng)的振幅，請(qǐng)問(wèn)這種情況合理嗎？怎么解釋呢？謝謝

HHT_spectrum02.png

HHT_spectrum03.png

贊一下(1人)

回復(fù)此樓

2樓2013-05-06 14:24:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiedt

金蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1299.7
帖子: 1904
在線: 628.6小時(shí)
蟲號(hào): 845974

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

祝一切順利，科研成果累累！順便搶金幣！

贊一下(1人)

回復(fù)此樓

3樓2013-05-06 18:12:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

tony1230

銀蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 511.9
帖子: 72
在線: 73.3小時(shí)
蟲號(hào): 1793983

引用回帖:

3樓: Originally posted by xiedt at 2013-05-06 18:12:42

祝一切順利，科研成果累累！順便搶金幣！

給別人留點(diǎn)吧，謝謝！

贊一下(1人)

回復(fù)此樓

4樓2013-05-06 21:20:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

rock_yang

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 619.3
帖子: 49
在線: 17.8小時(shí)
蟲號(hào): 1348710

至于高頻的部分，你還要看不同的IMF的能量掉高頻部分是不一要的，不能只簡(jiǎn)單地說(shuō)去就完事了，這個(gè)還是要參考邊際譜

贊一下(1人)

回復(fù)此樓

7樓2014-02-28 13:50:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

zwg19870824

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 14.8
帖子: 6
在線: 5.1小時(shí)
蟲號(hào): 3550375

引用回帖:

2樓: Originally posted by tony1230 at 2013-05-06 14:24:54
我把上面的頻譜濾掉了最高頻的兩個(gè)信號(hào)，依次得到下面兩個(gè)頻譜圖，但是高頻部分為何還留有比較強(qiáng)的振幅，請(qǐng)問(wèn)這種情況合理嗎？怎么解釋呢？謝謝

HHT_spectrum02.png

HHT_spectrum03.png
...

值得稱贊

回復(fù)此樓

8樓2015-01-18 14:40:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiao又又

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 14
帖子: 4
在線: 2.2小時(shí)
蟲號(hào): 4394812

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

你好，我也是新手一枚，我通過(guò)希爾伯特變換后得到許多負(fù)頻，怎么處理

贊一下

回復(fù)此樓

9樓2016-02-14 18:56:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

allen121

禁蟲 (初入文壇)

本帖內(nèi)容被屏蔽

10樓2017-06-02 21:48:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

藍(lán)鴉

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 666.9
帖子: 167
在線: 138.4小時(shí)
蟲號(hào): 4256922

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

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

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

11樓2017-06-04 11:09:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 317求調(diào)劑 +4	申子申申 2026-03-19	9/450	2026-03-20 11:08 by 申子申申
[考研] 265求調(diào)劑 +6	梁梁校校 2026-03-17	6/300	2026-03-20 08:59 by ZHANG0tao
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +3	Charlieyq 2026-03-19	3/150	2026-03-20 08:53 by xingguangj
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-19	3/150	2026-03-19 13:36 by houyaoxu
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 0817調(diào)劑 +3	沒(méi)有答案_ 2026-03-14	3/150	2026-03-19 09:51 by Xu de nuo
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無(wú)際的草原
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike