emd-希尔伯特变换包络谱

对imf做希尔伯特 fft，包括幅值和频率的 包络

function imf = emd（x）
% Empiricial Mode Decomposition （Hilbert-Huang Transform）
% EMD分解或HHT变换
% 返回值为cell类型，依次为一次IMF、二次IMF、...、最后残差

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;

% 是否单调
function u = ismonotonic（x）
u1 = length（findpeaks（x））*length（findpeaks（-x））;
if u1 > 0
    u = 0;
else
    u = 1;
end

% 是否IMF分量
function u = isimf（x）
N  = length（x）;
u1 = sum（x（1:N-1）.*x（2:N） < 0）;                     % 过零点的个数
u2 = length（find

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-04-09 15:08  EMD\
     文件         227  2012-10-08 22:32  EMD\FFTAnalysis.m
     文件         250  2012-10-08 22:33  EMD\HilbertAnalysis.m
     文件        1108  2012-10-08 22:32  EMD\emd.m
     文件        2501  2017-02-02 20:34  EMD\exe_emd_20170201.asv
     文件        2129  2017-04-09 15:08  EMD\exe_emd_20170201.m
     文件        2595  2017-04-05 09:00  EMD\exe_emd_20170201_moni.asv
     文件         925  2012-10-08 22:30  EMD\findpeaks.asv
     文件         925  2012-10-08 22:30  EMD\findpeaks.m
     文件        1445  2007-06-27 22:12  EMD\hhspectrum.m
     文件        3427  2005-10-13 23:15  EMD\instfreq.m
     文件        1721  2012-10-12 15:48  EMD\plot_hht.asv
     文件        1755  2017-02-04 12:04  EMD\plot_hht.m
     文件         595  2012-10-11 22:50  EMD\test_bianjipu.m
     文件        1330  2017-01-26 17:01  EMD\test_emd.asv
     文件         442  2012-10-13 09:01  EMD\test_emd_shu.asv
     文件         274  2012-10-11 22:51  EMD\youyige_bianjipu.m