EMD

clear all;
close all;
clc;
% Time Domain 0 to T
T = 1000;
fs = 1/T;
t = (1:T)/T;
freqs = 2*pi*(t-0.5-1/T)/(fs);
x2=(cos(200*pi*t));
x3=sin(64*pi*t).*(t<0.5);
x=x2+x3+0.1*randn(size(x2));
figure;
plot(t,x);
imf=emd(x,0.1,100);
figure
m=size(imf,2);
for i =1:m
    subplot(m,1,i)
    plot(imf(:,i))
end

本文为原创文章，转载请注明出处！

打赏作者