xibu1102007 发表于 2012-3-26 10:22

EMD alan 注释版本交流

本帖最后由 xibu1102007 于 2012-3-26 11:04 编辑

谁下过这个程序啊，里面有些程序看不懂，主要是停止条件那块，求指导

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

x= transpose(x(:));               %转置为行矩阵
imf = [];

while ~ismonotonic(x)               %当x不是单调函数，分解终止条件
   x1 = x;
   sd = Inf;                        %均值 无穷大的意思？
                                    %直到x1满足IMF条件，得c1
   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);

   
      x1 = x2;
   end
   imf{end+1} = x1;
   x         = x-x1;
end

imf{end+1} = x;



这是分解程序。停止条件那块的。

xibu1102007 发表于 2012-3-26 10:30

沙发自己来做吗，论坛有米有EMD群啊？没有的话我C一个，希望大家都能讨论{:{39}:}

xibu1102007 发表于 2012-3-26 10:36

请加我QQ 746091962.在线等啊

xibu1102007 发表于 2012-3-26 11:19

可以看到这个程序的停止条件是判断单调性和SD的

xibu1102007 发表于 2012-3-26 15:27

不是吧，高手都不屑于理我么{:{40}:}

xibu1102007 发表于 2012-3-26 15:27

{:{10}:}

lovesea 发表于 2013-3-29 17:45

这个程序没有处理端点效应吧？
页: [1]
查看完整版本: EMD alan 注释版本交流