aprilcat 发表于 2010-1-20 16:24

回复 沙发 rhenus 的帖子

你的例子容易理解，但是运行时候几乎是死机状态。我做了修改，但是还是不能分解出你给出的结果图。
clear
t = 0:0.01:10;
x1 = sin(2*pi*t);
x2 = sin(2*pi*10*t);
sig = x1 + x2;
= emd(sig);
s=size(imf);
cn=s(1);
subplot(cn+1,1,1)%首先画出原始信号
plot(sig);
title('original sigal');
for k = 1:cn-1%循环画出imf分量
subplot(cn+1,1,k+1);
plot(imf(k,:));
title(['imf ',int2str(k)]);
end
subplot(cn+1,1,cn+1)%最后画出残余量
plot(r);
title('residual');

aprilcat 发表于 2010-1-20 16:32

回复 9楼 wains 的帖子

imf(:,k)形式改为imf(k,:);形式。

cloudy_time 发表于 2012-9-7 18:05

2楼程序有点问题，用这个试试：


t = 0:0.01:10;
x1 = sin(2*pi*t);
x2 = sin(2*pi*10*t);
sig = x1 + x2;
= emd(sig);
cn=size(imf,1);
%首先画出原始信号
subplot(cn+1,1,1)
plot(sig)
title('original signal')
%循环画出imf分量
for k=2:cn
subplot(cn+1,1,k)
plot(imf(k,:))
title(['imf',num2str(k-1)])
end
%最后画出残余量
subplot(cn+1,1,cn+1)
plot(imf(cn,:))
title('residual')

还是算了吧 发表于 2013-9-11 20:23

cloudy_time 发表于 2012-9-7 18:05 static/image/common/back.gif
2楼程序有点问题，用这个试试：




= emd(sig)；
请问里面的“r”代表什么意思

小笨蛋0516 发表于 2014-8-28 16:12

还是算了吧 发表于 2013-9-11 20:23
= emd(sig)；
请问里面的“r”代表什么意思

残余信号分量

归零 发表于 2015-11-14 09:54

小笨蛋0516 发表于 2014-8-28 16:12
残余信号分量

cn=size(imf,1);这一句什么意思？？{:{03}:}

筱儿 发表于 2016-5-14 20:58

nkdtxf 发表于 2009-10-29 09:01
clear
t = 0:0.01:10;
x1 = sin(2*pi*t);


emd_visu(s,t,imf,1);
这个函数您有吗，能发我一下吗
页: 1 [2]
查看完整版本: 求助EMD分解的简单例子