希尔伯特-黄变换 工具箱非常好用

算法用于非平稳随机信号的希尔伯特-黄变换处理。并附有实际算例。已验证，非常好用。

% This is an EMD/EEMD program
%
%   function allmode=eemd（YNstdNE）
%
% INPUT:
%       Y: Inputted data;
%       Nstd: ratio of the standard deviation of the added noise and that of Y;
%       NE: Ensemble number for the EEMD
% OUTPUT:
%       A matrix of N*（m+1） matrix where N is the length of the input
%       data Y and m=fix（log2（N））-1. Column 1 is the original data columns 2 3 ...
%       m are the IMFs from high to low frequency and comlumn （m+1） is the
%       residual （over all trend）.
%
% NOTE:
%       It should be noted that when Nstd is set to zero and NE is set to 1 the
%       program degenerates to a EMD program.
%
% References can be found in the “Reference“ section.
%
% The code is prepared by Zhaohua Wu. For questions please read the “Q&A“ section or
% contact
%   zhwu@cola.iges.org
%

function allmode=eemd（YNstdNE）
xsize=length（Y）;
dd=1:1:xsize;
Ystd=std（Y）;
Y=Y/Ystd;

TNM=fix（log2（xsize））-1;
TNM2=TNM+2;
for kk=1:1:TNM2 
    for ii=1:1:xsize
        allmode（iikk）=0.0;
    end
end

for iii=1:1:NE
    for i=1:xsize
        temp=randn（11）*Nstd;
        X1（i）=Y（i）+temp;
    end

    %for jj=1:1:xsize
    %    mode（jj1） = Y（jj）;
    %end
    
    xorigin = X1;
    xend = xorigin;
    
    nmode = 0;
    while nmode <= TNM
        xstart = xend;
        iter = 1;
   
        while iter<=10
            [spmax spmin flag]=extrema（xstart）;
            upper= spline（spmax（:1）spmax（:2）dd）;
            lower= spline（spmin（:1）spmin（:2）dd）;
            mean_ul = （upper + lower）/2;
            xstart = xstart - mean_ul;
            iter = iter +1;
        end
        xend = xend - xstart;
   
   	    nmode=nmode+1;
        
        for jj=1:1:xsize
            mode（jjnmode） = xstart（jj）;
        end
    end
   
    for jj=1:1:xsize
        mode（jjnmode+1）=xend（jj）;
    end
   
    allmode=allmode+mode;
    
end

allmode=allmode/NE;
allmode=allmode*Ystd;

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

     文件       7763  2015-06-09 19:50  HHT\4.37第1组特征波形1.asv

     文件       7615  2015-05-25 15:06  HHT\4.37第1组特征波形1.txt

     文件      20426  2015-05-25 15:06  HHT\4.8-4.16HHT模拟对比.txt

     文件       2014  2015-01-23 15:58  HHT\eemd.m

     文件       2883  2015-01-23 19:25  HHT\eemd.txt

     文件       7554  2015-06-09 20:58  HHT\eemd_plot.m

     文件      22275  2007-06-27 16:12  HHT\emd.m

     文件      23135  2015-01-25 12:47  HHT\emd.txt

     文件       7596  2015-06-09 20:58  HHT\emd_plot.m

     文件       1545  2015-06-09 20:03  HHT\fmlin.m

     文件       2931  2015-06-09 19:44  HHT\instfreq.m

     文件       1123  2015-04-18 13:44  HHT\plot_hht_3d.m

     文件       1643  2015-04-16 10:08  HHT\plot_hht_3d.txt

     文件      20649  2015-06-09 20:02  HHT\Untitled1.m

     文件       1913  1998-04-13 13:11  HHT\工具箱1-时频分析\altes.m

     文件       2253  1998-04-13 13:11  HHT\工具箱1-时频分析\ambifunb.m

     文件       4821  1998-05-19 12:08  HHT\工具箱1-时频分析\ambifuwb.m

     文件       1124  1998-04-13 13:11  HHT\工具箱1-时频分析\amexpo1s.m

     文件       1132  1998-04-13 13:11  HHT\工具箱1-时频分析\amexpo2s.m

     文件       1158  1998-04-13 13:11  HHT\工具箱1-时频分析\amgauss.m

     文件       1105  1998-04-13 13:11  HHT\工具箱1-时频分析\amrect.m

     文件       1138  1998-04-13 13:11  HHT\工具箱1-时频分析\amtriang.m

     文件       1475  1998-06-12 12:03  HHT\工具箱1-时频分析\anaask.m

     文件       1564  1998-06-12 12:03  HHT\工具箱1-时频分析\anabpsk.m

     文件       1604  1998-06-12 12:03  HHT\工具箱1-时频分析\anafsk.m

     文件        901  1998-04-13 13:11  HHT\工具箱1-时频分析\anapulse.m

     文件       1721  1999-10-26 17:59  HHT\工具箱1-时频分析\anaqpsk.m

     文件       1382  1998-04-13 13:11  HHT\工具箱1-时频分析\anasing.m

     文件        905  1998-04-13 13:11  HHT\工具箱1-时频分析\anastep.m

     文件       5491  1998-05-19 13:33  HHT\工具箱1-时频分析\ATOMS.M

............此处省略375个文件信息