关于Hilbert-Huang的matlab实现,材料汇总,比较杂...感谢所有网络上的贡献者们:)
核心:以下代码计算HHT边际谱及其对应频率
工具包要求:G-Rilling EMD Toolbox,TFTB Toolbox
附:黄锷先生课题组开发的工具包(可以在 这里 找到),这里并未用到。
% Empirical mode decomposition, resulting in intrinc mode functions.
% Without parameter 'MAXMODES' the process may be seriously delayed by
% decompose original signals into too many IMFs (not necessary, 9 is
% enough generally)
imfs = emd(oriSig, 'MAXMODES', 9);
% HHT spectrum: hhtS
[A, f, t] = hhspectrum(imfs);
[hhtS, ~, fCent] = toimage(A, f, t);
% Marginal hilbert spectrum: hhtMS, xf: correspondig frequency
for k = 1 : size(hhtS, 1)
hhtMS(k) = sum(hhtS(k, : )) * 1 / fs;
end
xf = fCent(1, :) .* fs;
G-Rilling EMD Toolbox工具包相关配置:了凡春秋
简单来说,设置好路径之后输入 install_emd 即可。