应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 数学综合 » 求高手解释HHT变换求频率等问题
111/1返回列表
查看: 4150  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tony1230

银虫 (小有名气)

[交流] 求高手解释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 ]
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2013-05-06 13:09:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

niuyun1217

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你好，请问HHT变换求频率你会了吗？可以把程序发给我吗？谢谢，急用
赞一下(2人) 回复此楼
5楼2013-12-29 14:19:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rock_yang

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
关于频率的部分的话，把计算得到的频率乘以采样频率就是实际频率了，至于为什么是0.5，你应该懂的

hilbert变换是对窄带信号才有效的，所以才有了EMD分解将任意信号分解成多个窄带信号

所以如果对一个窄带信号做hilbert变换的话，去掉负频率，double正频率就是hilbert谱了
赞一下(2人) 回复此楼
6楼2014-02-28 13:49:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tony1230

银虫 (小有名气)
我把上面的频谱滤掉了最高频的两个信号，依次得到下面两个频谱图，但是高频部分为何还留有比较强的振幅，请问这种情况合理吗？怎么解释呢？谢谢

HHT_spectrum02.png



HHT_spectrum03.png
赞一下(1人) 回复此楼
2楼2013-05-06 14:24:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiedt

金虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
祝一切顺利，科研成果累累！顺便抢金币！
赞一下(1人) 回复此楼
3楼2013-05-06 18:12:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tony1230

银虫 (小有名气)
引用回帖:
3楼: Originally posted by xiedt at 2013-05-06 18:12:42
祝一切顺利，科研成果累累！顺便抢金币！

给别人留点吧，谢谢！
赞一下(1人) 回复此楼
4楼2013-05-06 21:20:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rock_yang

新虫 (初入文坛)
至于高频的部分，你还要看不同的IMF的能量掉高频部分是不一要的，不能只简单地说去就完事了，这个还是要参考边际谱
赞一下(1人) 回复此楼
7楼2014-02-28 13:50:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zwg19870824

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by tony1230 at 2013-05-06 14:24:54
我把上面的频谱滤掉了最高频的两个信号，依次得到下面两个频谱图，但是高频部分为何还留有比较强的振幅，请问这种情况合理吗？怎么解释呢？谢谢

HHT_spectrum02.png

HHT_spectrum03.png
...

值得称赞
回复此楼
8楼2015-01-18 14:40:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao又又

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你好，我也是新手一枚，我通过希尔伯特变换后得到许多负频，怎么处理
赞一下 回复此楼
9楼2016-02-14 18:56:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

allen121

禁虫 (初入文坛)
本帖内容被屏蔽
10楼2017-06-02 21:48:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蓝鸦

新虫 (小有名气)
★
小木虫: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tony1230 的主题更新
111/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 复试调剂 +3 呼呼？~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 泣血叩求调剂恩，愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 一志愿山东大学，总分327，英语二79，有论文，有竞赛，已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
关闭
请填处理意见
关闭