資源簡介
這是一段eemd分解的程序,可以直接使用,可用于故障信號的分解,提取故障特征,進行故障判斷的一段程序。
代碼片段和文件信息
%?Y:?Inputted?data;
%?Nstd:?ratio?of?the?standard?deviation?of?the?added?noise?and?that?of?Y;
%?NE:?Ensemble?member?being?used
%?TNM:?total?number?of?modes?(not?including?the?trend)
%
function?allmode=ceemd(YNstdNETNM)
%?find?data?length
xsize=length(Y);
dd=1:1:xsize;
%?Nornaliz?data
Ystd=std(Y);
Y=Y/Ystd;
%?Initialize?saved?data
TNM2=TNM+2;
for?kk=1:1:TNM2
????for?ii=1:1:xsize
????????allmode(iikk)=0.0;
????end
end
for?iii=1:1:NE
%?adding?noise
????for?i=1:xsize
????????temp=randn(11)*Nstd;
????????X1(i)=Y(i)+temp;
????????X2(i)=Y(i)-temp;
????end
????%?sifting?X1
????xorigin?=?X1;
????xend?=?xorigin;
%?save?the?initial?data?into?the?first?column
????for?jj=1:1:xsize
????????mode(jj1)?=?xorigin(jj);
????end
????nmode?=?1;
????while?nmode?<=?T
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2603??2014-04-01?21:11??ceemd.m
評論
共有 條評論