MATLAB实现HHT变换的代码+

本人在网上找的一些用MATLAB实现HHT变换的程序，还附加了一些仿真的测试程序。

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

x   = transpose（x（:））;%表示对矩阵或者数组x进行转置
imf = [];
while ~ismonotonic（x）%当x不是单调函数，分解终止条件
    x1 = x;
    sd = Inf;%Inf-infinite表示无穷大的意思。实际意思是当算出的结果大于某个数
             %（这个数很大，比如10的很多次方），则MATLAB认为就是无穷大了，并返回 inf
    while （sd > 0.1） || ~isimf（x1）%当标准偏差系数sd大于0.1或x1不是固有模态函数时，分量终止条件
        s1 = getspline（x1）;         % 极大值点样条曲线
        s2 = -getspline（-x1）;       % 极小值点样条曲线
        x2 = x1-（s1+s2）/2;          % 此处的x2为文章中的h
       
        sd = sum（（x1-x2）.^2）/sum（x1.^2）;%Cauchy收敛性检测 用于判断抽取过程何时终止
        x1 = x2;
    end
   
    imf{end+1} = x1;
    x          = x-x1;
end
imf{end+1} = x;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-03-28 14:19  CYP—test\
     目录           0  2019-03-11 16:21  CYP—test\sample1\
     文件         286  2019-01-09 18:59  CYP—test\sample1\FFTAnalysis.m
     文件         419  2019-01-18 09:20  CYP—test\sample1\HilbertAnalysis.m
     文件         914  2019-01-09 14:42  CYP—test\sample1\emd.m
     文件         743  2019-01-05 10:50  CYP—test\sample1\fftfenxi.m
     文件         913  2019-01-06 10:49  CYP—test\sample1\findpeaks.m
     文件         156  2019-01-07 19:16  CYP—test\sample1\getspline.m
     文件        2262  2019-01-17 18:51  CYP—test\sample1\hhspectrum.m
     文件        1409  2019-03-07 10:55  CYP—test\sample1\high_pass_filtertest.m
     文件         292  2019-01-06 16:35  CYP—test\sample1\isimf.m
     文件         174  2019-01-06 16:34  CYP—test\sample1\ismonotonic.m
     文件        1684  2019-01-18 15:40  CYP—test\sample1\plot_hht.m
     文件        2858  2019-01-21 18:42  CYP—test\sample1\plot_hht1.m
     文件         635  2019-01-10 09:53  CYP—test\sample1\plot_imf.m
     文件        1358  2019-01-17 18:45  CYP—test\sample1\test.m
     目录           0  2019-03-07 20:38  CYP—test\sample3\
     文件       50624  2019-01-10 10:07  CYP—test\sample3\1.txt
     文件      327312  2019-03-07 20:38  CYP—test\sample3\8hz.mat
     文件         286  2019-01-09 18:59  CYP—test\sample3\FFTAnalysis.m
     文件         433  2019-01-09 19:19  CYP—test\sample3\HilbertAnalysis.m
     文件       33204  2019-03-06 15:19  CYP—test\sample3\T4.mat
     文件       18623  2019-03-06 15:20  CYP—test\sample3\T5.mat
     文件       14204  2019-03-06 15:20  CYP—test\sample3\T6.mat
     文件        6132  2019-03-06 15:22  CYP—test\sample3\T7.mat
     文件       11766  2019-03-06 15:22  CYP—test\sample3\T8.mat
     文件      104294  2019-03-06 15:23  CYP—test\sample3\T9.mat
     文件        2823  2019-01-10 10:00  CYP—test\sample3\cemd_visu.m
     文件        1583  2019-01-10 10:00  CYP—test\sample3\disp_hhs.m
     文件       22275  2019-01-10 10:11  CYP—test\sample3\emd.m
     文件         914  2019-02-21 10:46  CYP—test\sample3\emd1.m
............此处省略45个文件信息