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
本文为原创文章,转载请注明出处!