最近,由于科研工作需求,需要处理某一类波动信号,遂利用经验模态分解(EMD)、集合经验模态分解(EEMD)和希尔伯特-黄变换对上述信号进行处理,期望获取有用的本征模态函数(IMF)、频率信息和边际谱等。
在信号实施和处理时发现,网上大多代码均由较早时期的MATLAB完成,彼时还没有自带的HHT函数。而在MATLAB2018a及其以后的版本中,集成和自带了HHT函数。
请各位网友解答我的疑惑,不胜感激!!!
问题1:函数hhspectrum()返回的是否为归一化后的频率?为何?
函数instfreq中有如下注释
The result FNORMHAT lies between 0.0 and 0.5.
问题2:函数toimage()的原理是什么?为什么以各个IMF分量的瞬时幅值和瞬时频率为输入得到的im是N*N的矩阵(N表示信号长度)?
[im,tt,ff] = TOIMAGE(A,f,t,splx,sply) transforms a spectrum made of 1D functions (e.g., output of "hhspectrum") in an 2D image