Eminem 发表于 2016-5-9 09:48

Rilling的关于HHT端点延拓的代码

<font color="#000000">if indmax(1) < indmin(1)
      if m(1) > m(indmin(1))
      lmax = fliplr(indmax(2:min(end,NBSYM+1)));
      lmin = fliplr(indmin(1:min(end,NBSYM)));
      lsym = indmax(1);
      else
      lmax = fliplr(indmax(1:min(end,NBSYM)));
      lmin = ;
      lsym = 1;
      end
    else

      if m(1) < m(indmax(1))
      lmax = fliplr(indmax(1:min(end,NBSYM)));
      lmin = fliplr(indmin(2:min(end,NBSYM+1)));
      lsym = indmin(1);
      else
      lmax = ;
      lmin = fliplr(indmin(1:min(end,NBSYM)));
      lsym = 1;
      end
    end
   
    if indmax(end) < indmin(end)
      if m(end) < m(indmax(end))
      rmax = fliplr(indmax(max(end-NBSYM+1,1):end));
      rmin = fliplr(indmin(max(end-NBSYM,1):end-1));
      rsym = indmin(end);
      else
      rmax = ;
      rmin = fliplr(indmin(max(end-NBSYM+1,1):end));
      rsym = lx;
      end
    else
      if m(end) > m(indmin(end))
      rmax = fliplr(indmax(max(end-NBSYM,1):end-1));
      rmin = fliplr(indmin(max(end-NBSYM+1,1):end));
      rsym = indmax(end);
      else
      rmax = fliplr(indmax(max(end-NBSYM+1,1):end));
      rmin = ;
      rsym = lx;
      end
    end
   
    tlmin = 2*t(lsym)-t(lmin);
    tlmax = 2*t(lsym)-t(lmax);
    trmin = 2*t(rsym)-t(rmin);
    trmax = 2*t(rsym)-t(rmax);</font>转自:http://blog.sina.com.cn/s/blog_5def5a660100c278.html

lov_fore 发表于 2016-5-9 10:02

这个程序是在做EMD分解时,自动调用的吗?

hi凌寒 发表于 2016-5-18 17:38

问下楼主Rilling的 emd工具箱中,有对端点进行处理过吗

airazor 发表于 2016-8-29 20:12

hi凌寒 发表于 2016-5-18 17:38
问下楼主Rilling的 emd工具箱中,有对端点进行处理过吗

rilling用的是镜像延拓

zhangzy 发表于 2016-12-14 12:50

这个程序直接运行不了必须得有工具箱对吗
页: [1]
查看完整版本: Rilling的关于HHT端点延拓的代码