%?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;（添加噪聲的標(biāo)準(zhǔn)方差與原始數(shù)據(jù)標(biāo)準(zhǔn)方差的比值）
%???????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
%???zwu@fsu.edu
%
%



function?allmode=eemd（YNstdNE）
xsize=length（Y）;
dd=1:1:xsize;
Ystd=std（Y）;
Y=Y/Ystd;Y=Y/Ystd;
%標(biāo)準(zhǔn)差歸一化，這樣設(shè)定加入噪聲額標(biāo)準(zhǔn)差才有意義，最后還要結(jié)果還要乘上?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;

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件????????2156??2012-04-11?08:35??Copy_of_eemd.m
?????文件?????????727??2012-04-17?20:10??Data_Partition.m
?????文件?????????854??2008-03-27?22:46??dist_value.m
?????文件?????8205843??2012-06-13?22:33??doppler_data.mat
?????文件????????2172??2012-04-11?08:25??eemd.asv
?????文件????????2011??2008-03-27?21:10??eemd.m
?????文件????????8980??2012-04-26?09:25??EEMD_ANN.m
?????文件?????????374??2011-11-25?15:35??eemd_emd_quzao.asv
?????文件?????????403??2011-11-25?15:36??eemd_emd_quzao.m
?????文件????????2187??2008-03-27?21:11??extrema.m
?????文件????????1020??2012-04-17?22:19??Hilbert_FFT_Plot.m
?????文件????????1825??2008-03-27?21:11??ifndq.m
?????文件????????1047??2012-04-15?16:24??Node_Index_generate.m
?????文件????????1094??2012-04-10?21:54??read?me.txt
?????文件????????2945??2008-03-27?22:46??significance.m
?????文件????????1398??2012-04-18?22:50??zhu.m
?????文件?????1216738??2011-10-29?07:22??大型旋轉(zhuǎn)機(jī)械非平穩(wěn)振動(dòng)信號(hào)的EEMD降噪方法.pdf
?????文件??????886664??2011-04-09?15:35??集合經(jīng)驗(yàn)?zāi)Ｊ椒纸庠诓裼蜋C(jī)機(jī)械故障診斷中的應(yīng)用.pdf
?????文件?????1544371??2011-04-09?15:34??集合經(jīng)驗(yàn)?zāi)Ｊ椒纸庠谛D(zhuǎn)機(jī)械故障診斷中的應(yīng)用.pdf