untitled2.m

经改进的希尔伯特——黄分解

%t=linspace(0,512,1024)
%x=sin(2/30*pi*t)+sin(2/34*pi*t)
t=linspace(0,400,400)
x=cos(2*pi*0.05*t)+cos(2*pi*0.02*t)
imf=emd(x,t)

k=5
for v=1:k
subplot(k,1,v),plot(t,imf(v,:))
ylabel('IMF ')
v=v+1
end