Hilbert分析和总结分析和总结.pdf

Hilbert 关于 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 即可。 Matlab 关于 EMD 分解后希尔伯特谱分析（相关函数均隶属 G-Rilling EMD Toolbox ） hhspectrum 函数说明（ 8 楼：老老的学生） % [A,f,tt] = HHSPECTRUM(imf,t,l,aff) % Input: %- imf : matrix with one IMF per row % 将 emd 分解得到 的 IMF 代入就可以， 就是你的程序中写的 c 变量， 不用加最 后一行的趋势项 %- t : time instants % 瞬时时间或持续时间 ？？（写 [1:信号长度 ] 就可以，真实的时间可以根据采样率转换 %- l : estimation parameter for instfreq % 瞬时频率的 估计参数 ？？（写 1 就可以，决定瞬时频率估计时的边界 从第几个点开始 %- aff : if 1, displays the computation evolution % 显示计 算进程选项，不想显示写 0 就可以 % % Output: % - A : amplitudes of IMF rows % - f : instantaneous frequencies % - tt : truncated time instants % 截止时间 ？？（截 断时间，返回的是瞬时频率对应的时间，要比原来信号的时 间按短，由前面的 l 值决定） % % calls: % - hilbert : computes the analytic signal % - instfreq : computes the instantaneous frequency % 瞬时频率 % % E