G. Rilling EMD分解程序端点延拓问题

huangqiang · 发表于 2009-9-17 10:43

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

G. Rilling的EMD分解边界处理函数  function [tmin,tmax,zmin,zmax] = boundary_conditions(indmin,indmax,t,x,z,nbsym)
有这么两句：
% 将序列根据对称中心，镜像到两边
tlmin = 2*t(lsym)-t(lmin);  %左边镜像
tlmax = 2*t(lsym)-t(lmax);
trmin = 2*t(rsym)-t(rmin); %右边镜像
trmax = 2*t(rsym)-t(rmax);

但分解结果端点效应的处理右端比左端要好，这是为何？
fs=200;
t = 1/fs:1/fs:1;
a=cos(2*pi*100*t);
b=cos(2*pi*30*t);
c=cos(2*pi*10*t);
s = a+b+c;
OPTIONS.FIX = 50;  %设置迭代计算筛分次数
imf=emd(s,OPTIONS);
err1=a-imf(1,:);
err2=b-imf(2,:);
err3=c-imf(3,:);

figure(1)
subplot(411)
plot(t,s);
subplot(412)
plot(t,err1)
subplot(413)
plot(t,err2)
subplot(414)
plot(t,err3)

[ 本帖最后由 huangqiang 于 2009-9-17 11:07 编辑 ]

huangqiang · 发表于 2009-9-17 10:55

把左边镜像做一下小改动，结果好像稍微好一些
% 将序列根据对称中心，镜像到两边
tlmin = -t(lmin);
tlmax = -t(lmax);

trmin = 2*t(rsym)-t(rmin);
trmax = 2*t(rsym)-t(rmax);
不解，盼高人指点。

pasuka · 发表于 2009-9-17 11:13

原帖由 huangqiang 于 2009-9-17 10:43 发表

 登录/注册后可看大图

G. Rilling的EMD分解边界处理函数 function [tmin,tmax,zmin,zmax] = boundary_conditions(indmin,indmax,t,x,z,nbsym)
有这么两句：
% 将序列根据对称中心，镜像到两边
tlmin = 2*t(lsym)-t(lmin); %左边镜像
...

lz的信号太完美了，换个别的，估计就是另一种结果了

huangqiang · 发表于 2009-9-17 11:36

这个信号分解已经相当糟糕了，
但仍然可以看出右端效果较好。
fs=200;
t = 1/fs:1/fs:1;
a=cos(2*pi*80*t);
b=cos(2*pi*30*t);
c=cos(2*pi*10*t);

huangqiang · 发表于 2009-9-17 14:06

呵呵……已经知道问题出哪儿了。

七七seven · 发表于 2009-9-17 19:42

LZ可以分享一下啊
问题出在哪

excalibur1982 · 发表于 2009-11-5 19:58

问题出在哪？楼主分享一下啊

nkdtxf · 发表于 2010-1-4 17:27

就是啊，原因出在哪里呢？

吃书的老虎 · 发表于 2010-1-4 22:38

G. Rilling的端点延拓不仅仅是对称延拓，他是根据端点和第一个极值点的关系再决定延拓的方法。楼主的问题我估计是左右端点不是同一种情况导致的。因为信号是周期的，可以将端点取成相同的，延拓应该左右情况会相同。

zhjinctgu · 发表于 2011-12-24 10:57

哦，原来如此~

xieqiyu1063 · 发表于 2012-11-2 09:41

楼主太监了

lovesea · 发表于 2013-4-1 16:11

本帖最后由 lovesea 于 2013-4-1 16:18 编辑

楼主，问题出在哪了？分享下啊！是信号取的不对称的结果吧？

华电机械 · 发表于 2013-8-18 17:34

lz现在对端点处理进展的怎么样，，有没有做出很好的改进，效果怎么样，贴图讨论下

账号		自动登录	找回密码
密码			我要加入

[HHT] G. Rilling EMD分解程序端点延拓问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

本帖被以下淘专辑推荐:

评分