%?This?is?an?EMD/EEMD?program
%
%???function?allmode=eemd（YNstdNE）
%
%?INPUT:
%???????Y:?Inputted?data;
%???????Nstd:?ratio?of?the?standard?deviation?of?the?added?noise?and?that?of?Y;
%???????NE:?Ensemble?number?for?the?EEMD
%?OUTPUT:
%???????A?matrix?of?N*（m+1）?matrix?where?N?is?the?length?of?the?input
%???????data?Y?and?m=fix（log2（N））-1.?Column?1?is?the?original?data?columns?2?3?...
%???????m?are?the?IMFs?from?high?to?low?frequency?and?comlumn?（m+1）?is?the
%???????residual?（over?all?trend）.
%
%?NOTE:
%???????It?should?be?noted?that?when?Nstd?is?set?to?zero?and?NE?is?set?to?1?the
%???????program?degenerates?to?a?EMD?program.
%
%?References?can?be?found?in?the?“Reference“?section.
%
%?The?code?is?prepared?by?Zhaohua?Wu.?For?questions?please?read?the?“Q&A“?section?or
%?contact
%???zhwu@cola.iges.org
%

function?allmode=eemd（YNstdNE）
xsize=length（Y）;
dd=1:1:xsize;
Ystd=std（Y）;
Y=Y/Ystd;

TNM=fix（log2（xsize））-1;
TNM2=TNM+2;
for?kk=1:1:TNM2?
????for?ii=1:1:xsize
????????allmode（iikk）=0.0;
????end
end

for?iii=1:1:NE
????for?i=1:xsize
????????temp=randn（11）*Nstd;
????????X1（i）=Y（i）+temp;
????end

????for?jj=1:1:xsize
????????mode（jj1）?=?Y（jj）;
????end
????
????xorigin?=?X1;
????xend?=?xorigin;
????
????nmode?=?1;
????while?nmode?<=?TNM
????????xstart?=?xend;
????????iter?=?1;
???
????????while?iter<=10
????????????[spmax?spmin?flag]=extrema（xstart）;
????????????upper=?spline（spmax（:1）spmax（:2）dd）;
????????????lower=?spline（spmin（:1）spmin（:2）dd）;
????????????mean_ul?=?（upper?+?lower）/2;
????????????xstart?=?xstart?-?mean_ul;
????????????iter?=?iter?+1;
????????end
????????xend?=?xend?-?xstart;
???
???	????nmode=nmode+1;
????????
????????for?jj=1:1:xsize
????????????mode（jjnmode）?=?xstart（jj）;
????????end
????end
???
????for?jj=1:1:xsize
????????mode（jjnmode+1）=xend（jj）;
????end
???
????allmode=allmode+mode;
????
end

allmode=allmode/NE;
allmode=allmode*Ystd;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????2928192??2018-08-06?10:55??EMD\130.mat

?????文件???????2011??2018-10-16?09:26??EMD\eemd.m

?????文件?????838457??2018-09-28?08:57??EMD\EMD.docx

?????文件??????22275??2015-05-29?10:37??EMD\emd.m

?????文件????1122115??2018-09-27?11:18??EMD\emd.png

?????文件?????905216??2018-09-27?10:00??EMD\emd.ppt

?????文件????????435??2015-05-31?11:20??EMD\exzample.m

?????文件???????2326??2015-05-31?15:55??EMD\hilbert.m

?????文件????????427??2015-05-31?16:04??EMD\hua_baoluo.m

?????文件???????1435??2015-05-31?11:24??EMD\hua_fft.m

?????文件???????2010??2018-10-16?16:01??EMD\innerguzhang.m

????..A..H.???????162??2018-09-27?10:30??EMD\~$EMD.docx

?????目錄??????????0??2018-10-16?15:38??EMD

-----------?---------??----------?-----??----

??????????????5825061????????????????????13