HHT变换的三种方法 Matlab

三种常用的在MATLAB上实现希尔伯特黄变换的仿真代码，还有一些简单的帮助文档

function imf = emd（x）
% Empiricial Mode Decomposition （Hilbert-Huang Transform）
% imf = emd（x）
% Func : findpeaks

x   = transpose（x（:））;
imf = [];
while ~ismonotonic（x）
   x1 = x;
   sd = Inf;
   while （sd > 0.1） | ~isimf（x1）
      s1 = getspline（x1）;
      s2 = -getspline（-x1）;
      x2 = x1-（s1+s2）/2;
      
      sd = sum（（x1-x2）.^2）/sum（x1.^2）;
      x1 = x2;
   end
   
   imf{end+1} = x1;
   x          = x-x1;
end
imf{end+1} = x;

% FUNCTIONS

function u = ismonotonic（x）

u1 = length（findpeaks（x））*length（findpeaks（-x））;
if u1 > 0 u = 0;
else      u = 1; end

function u = isimf（x）

N  = length（x）;
u1 = sum（x（1:N-1）.*x（2:N） < 0）;
u2 = length（findpeaks（x））+length（findpeaks（-x））;
if abs（u1-u2） > 1 u = 0;
else              u = 1; end

function s = getspline（x）

N = length（x）;
p = findpeaks（x）;
s = spline（[0 p N+1][0 x（p） 0]1:N）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2010-04-16 00:30  复件 HHT变换的三种方法 Matlab\
     目录           0  2019-01-09 20:14  复件 HHT变换的三种方法 Matlab\G Rilling\
     文件       86016  2009-04-23 18:20  复件 HHT变换的三种方法 Matlab\G Rilling\EMD程序使用.doc
     文件       39936  2010-04-13 12:37  复件 HHT变换的三种方法 Matlab\G Rilling\document.doc
     文件      905216  2009-04-22 23:23  复件 HHT变换的三种方法 Matlab\G Rilling\emd.ppt
     文件       55088  2009-04-22 23:19  复件 HHT变换的三种方法 Matlab\G Rilling\pack_emd.tar.gz
     文件      100530  2009-04-22 23:20  复件 HHT变换的三种方法 Matlab\G Rilling\pack_emd.zip
     文件     1697793  2009-04-22 23:26  复件 HHT变换的三种方法 Matlab\G Rilling\tftb-0.2.tar.gz
     文件         162  2009-04-22 23:22  复件 HHT变换的三种方法 Matlab\G Rilling\~$cument.doc
     文件       30208  2009-04-22 23:27  复件 HHT变换的三种方法 Matlab\G Rilling\~WRL2980.tmp
     文件       31744  2009-04-22 23:33  复件 HHT变换的三种方法 Matlab\G Rilling\~WRL3223.tmp
     文件       31232  2009-04-22 23:25  复件 HHT变换的三种方法 Matlab\G Rilling\~WRL3695.tmp
     文件      962513  2009-04-23 18:14  复件 HHT变换的三种方法 Matlab\G Rilling\归一化频率.mht
     文件      344727  2009-04-23 18:12  复件 HHT变换的三种方法 Matlab\G Rilling\最原始的EMD程序.mht
     文件      173863  2009-04-23 22:56  复件 HHT变换的三种方法 Matlab\G Rilling\正弦插值法解决emd的端点.mht
     文件      275752  2009-04-23 18:00  复件 HHT变换的三种方法 Matlab\G Rilling\镜像延拓.mht
     目录           0  2010-04-16 00:30  复件 HHT变换的三种方法 Matlab\Matlab File Exchange Center\
     文件      476047  2009-04-22 23:17  复件 HHT变换的三种方法 Matlab\Matlab File Exchange Center\MATLAB Central - File detail - Hilbert-Huang Transform.mht
     目录           0  2010-04-16 00:31  复件 HHT变换的三种方法 Matlab\Matlab File Exchange Center\plot_hht\
     文件      958087  2007-01-25 14:59  复件 HHT变换的三种方法 Matlab\Matlab File Exchange Center\plot_hht\HHT.pdf
     文件      102892  2007-04-28 15:57  复件 HHT变换的三种方法 Matlab\Matlab File Exchange Center\plot_hht\Hum.wav
     文件         893  2008-03-13 11:49  复件 HHT变换的三种方法 Matlab\Matlab File Exchange Center\plot_hht\emd.m
     文件         146  2007-08-26 00:42  复件 HHT变换的三种方法 Matlab\Matlab File Exchange Center\plot_hht\findpeaks.m
     文件         965  2008-04-24 13:36  复件 HHT变换的三种方法 Matlab\Matlab File Exchange Center\plot_hht\plot_hht.m
     目录           0  2010-04-16 00:30  复件 HHT变换的三种方法 Matlab\National Taiwan Central University\
     文件      541035  2009-04-23 18:38  复件 HHT变换的三种方法 Matlab\National Taiwan Central University\HHT MATLAB PROGRAM.mht
     文件         854  2009-04-23 18:37  复件 HHT变换的三种方法 Matlab\National Taiwan Central University\dist_value.m
     文件        2011  2010-04-12 21:36  复件 HHT变换的三种方法 Matlab\National Taiwan Central University\eemd.m
     目录           0  2019-01-06 10:33  复件 HHT变换的三种方法 Matlab\National Taiwan Central University\example and help\
     文件        4884  2008-10-21 19:09  复件 HHT变换的三种方法 Matlab\National Taiwan Central University\example and help\gsta.dat
     文件      161280  2008-10-21 19:32  复件 HHT变换的三种方法 Matlab\National Taiwan Central University\example and help\help.doc
............此处省略7个文件信息