基于EMD分解·和希尔伯特变换（HHT）的完整 MATLAB代码.rar

Hilbert-Huang 变换是一种适用于分析非线性、非平稳信号的数据处理方法，它是由美籍华人 Huang 以及他的同事在 1998 年提出的，从本质上讲这种方法是要对一个信号进行平稳化处理，得到信号的时间-频率-能量特征。HHT 是近年来在信号处理领域中的一项重要突破。HHT 是分 EMD 和 Hilbert 变换两步来实现的，首先对非线性、非平稳信号进行 EMD 分解，逐级分解出原始信号中不同尺度的波动或变化趋势，这些具有不同特征尺度的一系列时间序列分量叫做本征模态函数(IMF)，接着对每个 IMF 分量进行 Hilbert 变换。对于 EMD 分解得到的每个分量都有着不同的频率成分，通过对各分量的 Hilbert 变换能够得到具有物理意义的瞬时属性参数。 Hilbert谱表示的是信号幅值在整个频率段上随时间和频率的变化规律，Hilbert边际谱表示信号幅值在整个频率段上随频率的变化情况，它相当于傅里叶谱，但比傅里叶谱具有更高的频率分辨率。Hilbert边际谱是通过对Hilbert谱积分得到的。

function disp_hhs（imtinf）

% DISP_HHS（imtinf）
% displays in a new figure the spectrum contained in matrix “im“
% （amplitudes in log）.
%
% inputs : - im : image matrix （e.g. output of “toimage“）
%          - t （optional） : time instants （e.g. output of “toimage“） 
%          - inf （optional） : -dynamic range in dB （wrt max）
%            default : inf = -20
%
% utilisation : disp_hhs（im） ; disp_hhs（imt） ; disp_hhs（iminf） 
%              disp_hhs（imtinf）

figure
colormap（bone）
colormap（1-colormap）;

if nargin==1
  inf=-20;
  t = 1:size（im2）;

end

if nargin == 2
  if length（t） == 1
    inf = t;
    t = 1:size（im2）;
  else
    inf = -20;
  end
end

if inf >= 0
  error（‘inf doit etre < 0‘）
end

M=max（max（im））;

im = log10（im/M+1e-300）;

inf=inf/10;


imagesc（tfliplr（（1:size（im1））/（2*size（im1）））im[inf0]）;
set（gca‘YDir‘‘normal‘）
xlabel（[‘time‘]）
ylabel（[‘normalized frequency‘]）
title（‘Hilbert-Huang spectrum‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        936  2005-07-11 17:30  希尔伯特变换（HHT）的 完整 MATLAB程序\disp_hhs.m

     文件       7987  2019-08-07 10:46  希尔伯特变换（HHT）的 完整 MATLAB程序\emd.m

     文件       1854  2019-08-07 10:46  希尔伯特变换（HHT）的 完整 MATLAB程序\emd_fmsin.m

     文件       9869  2004-06-04 06:49  希尔伯特变换（HHT）的 完整 MATLAB程序\emd_local.m

     文件       3120  2004-06-04 06:25  希尔伯特变换（HHT）的 完整 MATLAB程序\emd_n.m

     文件      26138  2004-06-04 06:49  希尔伯特变换（HHT）的 完整 MATLAB程序\emd_online.m

     文件        891  2004-06-04 06:51  希尔伯特变换（HHT）的 完整 MATLAB程序\emd_sampling.m

     文件       1262  2004-06-04 06:51  希尔伯特变换（HHT）的 完整 MATLAB程序\emd_separation.m

     文件        676  2005-07-11 17:33  希尔伯特变换（HHT）的 完整 MATLAB程序\emd_triang.m

     文件       2491  2004-06-04 06:50  希尔伯特变换（HHT）的 完整 MATLAB程序\emd_visu.m

     文件       1777  2005-07-11 17:29  希尔伯特变换（HHT）的 完整 MATLAB程序\extr.m

     文件       2063  2019-08-07 10:46  希尔伯特变换（HHT）的 完整 MATLAB程序\ex_online.m

     文件        316  2017-12-21 21:28  希尔伯特变换（HHT）的 完整 MATLAB程序\fftfenxi.m

     文件        845  2004-06-04 06:49  希尔伯特变换（HHT）的 完整 MATLAB程序\hhspectrum.m

     文件       4543  2019-08-15 13:16  希尔伯特变换（HHT）的 完整 MATLAB程序\HHT.m

     文件       2071  2004-06-04 06:25  希尔伯特变换（HHT）的 完整 MATLAB程序\hspec.m

     文件        335  2005-07-11 17:29  希尔伯特变换（HHT）的 完整 MATLAB程序\io.m

     文件       1889  2005-07-11 17:30  希尔伯特变换（HHT）的 完整 MATLAB程序\toimage.m

     文件        378  2005-07-11 17:34  希尔伯特变换（HHT）的 完整 MATLAB程序\triang.m

     目录          0  2019-08-15 15:16  希尔伯特变换（HHT）的 完整 MATLAB程序

----------- ---------  ---------- -----  ----

                69659                    21