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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 数学综合 » 求高手解释HHT变换求频率等问题
111/1返回列表
查看: 4608  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗？20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭