声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1763|回复: 7

[HHT] HHT边际谱不太了解

[复制链接]
发表于 2014-12-10 19:10 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我的程序如下,在进行HHT边际谱分析的时候,为什么采样频率对我的边际谱图有影响,那么我怎么才能知道实际的频率是多大呢,很困惑,请求大家帮忙。我把程序和图黏贴如下:
  1. clear;
  2. fs=10;
  3. y=xlsread('data3.xlsx','B1:B488');
  4. z=y;
  5. imf=emd(z);
  6. [A,fa,tt]=hhspectrum(imf);
  7. [E,tt1,Cenf]=toimage(A,fa,tt,length(tt));
  8. cemd_visu(z,1:length(z),imf);
  9. figure(2)
  10. for k=1:size(E,1)
  11. bjp(k)=sum(E(k,:))*1/fs;
  12. end
  13. plot(Cenf(1,:)*fs,bjp);
  14. xlabel('频率/ Hz');
  15. ylabel('幅值');
复制代码
fs=10时和fs=50时结果就不一样,虽然图形一样,但是数值大小却变了。。。

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2014-12-11 10:55 | 显示全部楼层
对信号进行EMD分解时.为了得到一个IMF,是通过不断从信号中减去其局部均值而得到的,而局部均值是通过对信号极值点作三次样条插值求包络平均得到的,因此信号极值点的准确与否直接影响着EMD的分解结果,为了对信号正确进行EMD分解,就要使采样得到的信号的局部极值点不能太过偏离原信号的局部极值点,一方面,极值点的偏离会造成信号幅度的减小也即能量的减小,另一方面,极值点相对于原信号的不规则偏离(对应非整数倍采样率的情况)会导致EMD分解产生虚假模式。
为了提高极值点拟合的准确性,一般经验式地把采样率至少取为被分析信号最高频率的4倍以上。所以,采样率影响EMD的分解结果,也会影响后续的计算结果。
发表于 2014-12-11 11:10 | 显示全部楼层
Sampling effects on the Empirical Mode Decomposition,看看这篇文献。
发表于 2014-12-11 16:29 | 显示全部楼层

    bjp(k)=sum(E(k,:))/N;
这个是正确的。你的程序第11行有问题,我也是在反复调试中发现了这个问题。建议你先将用模拟信号调试一下程序,这样容易发现问题。
 楼主| 发表于 2014-12-11 22:13 | 显示全部楼层
airazor 发表于 2014-12-11 16:29
bjp(k)=sum(E(k,:))/N;
这个是正确的。你的程序第11行有问题,我也是在反复调试中发现了这个问题。 ...

好的,谢谢你的意见。。我在思考思考
 楼主| 发表于 2014-12-11 22:13 | 显示全部楼层
shuihai707 发表于 2014-12-11 11:10
Sampling effects on the Empirical Mode Decomposition,看看这篇文献。

非常感谢。。我来找看看
 楼主| 发表于 2014-12-11 22:15 | 显示全部楼层
shuihai707 发表于 2014-12-11 10:55
对信号进行EMD分解时.为了得到一个IMF,是通过不断从信号中减去其局部均值而得到的,而局部均值是通过对信 ...

主要是因为一开始我的的只是一些数据点,对于这些点内部隐藏的信号的频率具体多大也不清楚啊,我就直接按照时间采样步长的倒数作为采样频率来做的。。所以也很困惑
 楼主| 发表于 2014-12-17 08:52 | 显示全部楼层
仍小呆 发表于 2014-12-11 22:13
好的,谢谢你的意见。。我在思考思考

对了,你这个N是多少呢,我的程序里没有说到N啊。。谢谢
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 21:51 , Processed in 0.054686 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表