楼主,我是刚接触HHT,可是按照你的第三个程序(3.EEMD的一些问题)运行,怎么结果跟你的不一样啊(EEMD的分解图)?
clear;clc
t=1:1000;
t1=t/100*2*pi;
a1=sin(t1);
t2=t/10*2*pi;
b1=linspace(0,0,1000);
for i=250:350
b1(i)=0.2*sin(t2(i));
end
for i=750:850
b1(i)=0.2*sin(t2(i));
end
x=a1+b1; %x是原信号
subplot(311);plot(t,a1);
subplot(312);plot(t,b1);
subplot(313);plot(t,x);
plot(x)
imf=emd(x);
emd_visu(x,t,imf)
%eemd
imf_eemd=eemd(x,0.1,100);
figure
subplot(511);plot(imf_eemd(:,1))
subplot(512);plot(imf_eemd(:,2))
subplot(513);plot(imf_eemd(:,3))
subplot(514);plot(imf_eemd(:,4))
subplot(515);plot(imf_eemd(:,5))
figure
subplot(511);plot(imf_eemd(:,6))
subplot(512);plot(imf_eemd(:,7))
subplot(513);plot(imf_eemd(:,8))
subplot(514);plot(imf_eemd(:,9))
subplot(515);plot(imf_eemd(:,10))
EEMD分解(2)
|