%?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?=?0;
????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;

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

?????文件???????7763??2015-06-09?19:50??HHT\4.37第1組特征波形1.asv

?????文件???????7615??2015-05-25?15:06??HHT\4.37第1組特征波形1.txt

?????文件??????20426??2015-05-25?15:06??HHT\4.8-4.16HHT模擬對比.txt

?????文件???????2014??2015-01-23?15:58??HHT\eemd.m

?????文件???????2883??2015-01-23?19:25??HHT\eemd.txt

?????文件???????7554??2015-06-09?20:58??HHT\eemd_plot.m

?????文件??????22275??2007-06-27?16:12??HHT\emd.m

?????文件??????23135??2015-01-25?12:47??HHT\emd.txt

?????文件???????7596??2015-06-09?20:58??HHT\emd_plot.m

?????文件???????1545??2015-06-09?20:03??HHT\fmlin.m

?????文件???????2931??2015-06-09?19:44??HHT\instfreq.m

?????文件???????1123??2015-04-18?13:44??HHT\plot_hht_3d.m

?????文件???????1643??2015-04-16?10:08??HHT\plot_hht_3d.txt

?????文件??????20649??2015-06-09?20:02??HHT\Untitled1.m

?????文件???????1913??1998-04-13?13:11??HHT\工具箱1-時頻分析\altes.m

?????文件???????2253??1998-04-13?13:11??HHT\工具箱1-時頻分析\ambifunb.m

?????文件???????4821??1998-05-19?12:08??HHT\工具箱1-時頻分析\ambifuwb.m

?????文件???????1124??1998-04-13?13:11??HHT\工具箱1-時頻分析\amexpo1s.m

?????文件???????1132??1998-04-13?13:11??HHT\工具箱1-時頻分析\amexpo2s.m

?????文件???????1158??1998-04-13?13:11??HHT\工具箱1-時頻分析\amgauss.m

?????文件???????1105??1998-04-13?13:11??HHT\工具箱1-時頻分析\amrect.m

?????文件???????1138??1998-04-13?13:11??HHT\工具箱1-時頻分析\amtriang.m

?????文件???????1475??1998-06-12?12:03??HHT\工具箱1-時頻分析\anaask.m

?????文件???????1564??1998-06-12?12:03??HHT\工具箱1-時頻分析\anabpsk.m

?????文件???????1604??1998-06-12?12:03??HHT\工具箱1-時頻分析\anafsk.m

?????文件????????901??1998-04-13?13:11??HHT\工具箱1-時頻分析\anapulse.m

?????文件???????1721??1999-10-26?17:59??HHT\工具箱1-時頻分析\anaqpsk.m

?????文件???????1382??1998-04-13?13:11??HHT\工具箱1-時頻分析\anasing.m

?????文件????????905??1998-04-13?13:11??HHT\工具箱1-時頻分析\anastep.m

?????文件???????5491??1998-05-19?13:33??HHT\工具箱1-時頻分析\ATOMS.M

............此處省略375個文件信息