count_1126 发表于 2007-3-19 23:51

求助 EMD程序 下面代码怎么理解?

% boundary conditions for interpolations :
      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);
      % in case symmetrized parts do not extend enough
    if tlmin(1) > t(1) | tlmax(1) > t(1)
      if lsym == indmax(1)
      lmax = fliplr(indmax(1:min(end,NBSYM)));
      else
      lmin = fliplr(indmin(1:min(end,NBSYM)));
      end
      if lsym == 1
      error('bug')
      end
      lsym = 1;
      tlmin = 2*t(lsym)-t(lmin);
      tlmax = 2*t(lsym)-t(lmax);
    end   
   
    if trmin(end) < t(lx) | trmax(end) < t(lx)
      if rsym == indmax(end)
      rmax = fliplr(indmax(max(end-NBSYM+1,1):end));
      else
      rmin = fliplr(indmin(max(end-NBSYM+1,1):end));
      end
      if rsym == lx
      error('bug')
      end
      rsym = lx;
      trmin = 2*t(rsym)-t(rmin);
      trmax = 2*t(rsym)-t(rmax);
    end
         
    mlmax =m(lmax);
    mlmin =m(lmin);
    mrmax =m(rmax);
    mrmin =m(rmin);

[ 本帖最后由 zhangnan3509 于 2007-7-4 14:44 编辑 ]

zhangnan3509 发表于 2007-3-20 12:38

端点延拓

count_1126 发表于 2007-3-20 13:03

我知道是端点延拓呀 就是看不懂
拜托你解释下行吗
谢谢了

zhangnan3509 发表于 2007-3-20 14:24

在语法上没有什么问题 关键是你要画图 这样才能明白！看看这篇文章《希尔伯特——黄变换的端点延拓》

count_1126 发表于 2007-3-20 17:46

这个程序就是实现里面的算法吗?
我想找到这个程序的算法看下
烦告知
谢谢了!

zhangnan3509 发表于 2007-3-20 18:14

这篇文章里面有提到！你好好看看，如果你找不到这篇文章，请告诉我邮箱，晚上发送给你！:@)

count_1126 发表于 2007-3-22 18:41

这篇文章我有
那个镜像闭延拓看不明白
你能解释下吗

dudubear2009 发表于 2007-5-29 17:26

回复 #1 count_1126 的帖子

同问：其中的NBSYM， lsym是什么意思啊？
页: [1]
查看完整版本: 求助 EMD程序 下面代码怎么理解?