% EMD_TRIANG.M%% P. Flandrin, Mar. 13, 2003 - modified Mar. 2, 2006%% computes and displays EMD for the sum of 2 triangular % waveforms + 1 tone %% produces Figure 2 in%% G. Rilling, P. Flandrin and P. Gon峚lv弒% "On Empirical Mode Decomposition and its algorithms"% IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing% NSIP-03, Grado (I), June 2003N = 1024;% # of data samplest = 1:N;% triangular waveform 1p1 = fix(N/6);% periodx1 = triang(N,p1);% tonef0 = 0.03;% frequencyx2 = real(fmconst(N,f0));% triangular waveform 2p3 = 5;% periodx3 = triang(N,p3);x = x1 + 0.4*x2' + .4*x3;[imf,ort,nbits] = emd(x); % <-------- adapted to current emd.memd_visu(x,t,imf,1);figure(1)