Hilbert-Huang Transform: matlab 希尔伯特-黄变换: matlab实现

转载

于 2016-07-14 17:51:00 发布 · 9.2k 阅读

108 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/minks/p/5671246.html

文章标签：

#matlab

本文汇总了关于Hilbert-Huang Transform在MATLAB中的实现，主要介绍了使用G-Rilling EMD Toolbox进行希尔伯特谱分析的步骤和代码示例。讨论了HHT相比于小波等时频分析方法的独特之处，强调了HHT在处理非平稳信号时的优势，如Hilbert边际谱和包络谱的概念，并提供了相关学习笔记和资料链接。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

关于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;