最近在做希尔伯特黄变换,文献上说边际谱即是希尔伯特谱对时间的积分,那么如下简单代码是否能直接获得边际谱呢?
[imf,residual,info] = emd(X,'Interpolation','pchip'); %X是源信号,进行EMD分解
[hs,f,t]= hht(imf,fs); %对IMF进行希尔伯特谱计算
lms=sum(hs,2); %对希尔伯特谱时间积分
最近在做希尔伯特黄变换,文献上说边际谱即是希尔伯特谱对时间的积分,那么如下简单代码是否能直接获得边际谱呢?
[imf,residual,info] = emd(X,'Interpolation','pchip'); %X是源信号,进行EMD分解
[hs,f,t]= hht(imf,fs); %对IMF进行希尔伯特谱计算
lms=sum(hs,2); %对希尔伯特谱时间积分
由自带函数HHT得到hs的数据顺序是错的。时频谱矩阵相当于把时频谱行方向用频率切割,列方向用时刻切割,得出多个小方块,每一个方块对应的频率用中心频率表示,对应的时刻则记录数据的时刻,小方块里的数据则表示该时刻,该频率的能量值(振幅的平方)。原文在下面链接
————————————————
版权声明:本文为CSDN博主「jy-cheng」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/weixin_41406486/article/details/104775769