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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 数学综合 » 求高手解释HHT变换求频率等问题
111/1返回列表
查看: 4797  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭