求助EMD分解的简单例子

aprilcat

你的例子容易理解，但是运行时候几乎是死机状态。我做了修改，但是还是不能分解出你给出的结果图。
clear
t = 0:0.01:10;
x1 = sin(2*pi*t);
x2 = sin(2*pi*10*t);
sig = x1 + x2;
[imf,r] = 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

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

cloudy_time

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


t = 0:0.01:10;
x1 = sin(2*pi*t);
x2 = sin(2*pi*10*t);
sig = x1 + x2;
[imf,r] = 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')

还是算了吧

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

                               
登录/注册后可看大图

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

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

小笨蛋0516

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

残余信号分量

归零

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

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

筱儿

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

emd_visu(s,t,imf,1);
这个函数您有吗，能发我一下吗