資源簡(jiǎn)介
對(duì)非平穩(wěn)信號(hào)進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解,提取IMF,并進(jìn)行HILBERT變換。很好的程序

代碼片段和文件信息
%?clear;?
%?clc;?
%?P=[-1?-1?2?2?4;0?5?0?5?7];?T=[-1?-1?1?1?-1];???%利用minmax函數(shù)求輸入樣本范圍?
%?net?=?newff(minmax(P)[81]{‘tansig‘‘purelin‘}‘trainrp‘);???%5?隱含層神經(jīng)元個(gè)數(shù)
%?net.trainParam.show=50;
%?net.trainParam.lr=0.05;?
%?net.trainParam.epochs=300;
%?net.trainParam.goal=1e-5;?
%?[nettr]=train(netPT);?
%?net.iw{11}%隱層權(quán)值?
%?net.b{1}%隱層閾值??
%?net.lw{21}%輸出層權(quán)值?
%?net.b{2}%輸出層閾值???
%?sim(netP)?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%BP網(wǎng)絡(luò)的第一階段學(xué)習(xí)期(訓(xùn)練加權(quán)系數(shù)wki,wij)
%初始化
lr=0.05;err_goal=0.001;?????%lr為學(xué)習(xí)速率;err_goal為期望誤差最小值
max_epoch=10000;a=0.9;??????%max_epoch為訓(xùn)練的最大次數(shù);a為慣性系數(shù)?
Oi=0;Ok=0;??????????????????%置隱含層和輸出層各神經(jīng)元輸出初值為0
%提供兩組訓(xùn)練集和目標(biāo)值(3輸入,2輸出)
X=[1?1;-1?-1;1?1];T=[1?1;1?1];
%初始化wki,wij(M為輸入節(jié)點(diǎn)j的數(shù)量;q為隱含層節(jié)點(diǎn)i的數(shù)量;L為輸出節(jié)點(diǎn)k的數(shù)量)
[MN]=size(X);q=8;[LN]=size(T);?%N為訓(xùn)練集對(duì)數(shù)量
wij=rand(qM);wki=rand(Lq);
wij0=zeros(size(wij));wki0=zeros(size(wki));
for?epoch=1:max_epoch
????%計(jì)算隱含層各神經(jīng)元的輸出
????NETi=wij*X;
????for?j=1:N
????????for?i=1:q
????????????Oi(ij)=2/(1+exp(-NETi(ij)))-1;
????????end
????end
????%計(jì)算輸出層各神經(jīng)元輸出
????NETk=wki*Oi;
????for?i=1:N
????????for?k=1:L
????????????Ok(ki)=2/(1+exp(-NETk(ki)))-1;
????????end
????end
????%計(jì)算誤差函數(shù)
????E=((T-Ok)‘*(T-Ok))/2;
????if(E ????%調(diào)整輸出層加權(quán)系數(shù)
????deltak=Ok.*(1-Ok).*(T-Ok);
????w=wki;
????wki=wki+lr*deltak*Oi‘+a*(wki-wki0);
????wki0=w;
????%調(diào)整隱含層加權(quán)系數(shù)
????deltai=Oi.*(1-Oi).*(deltak‘*wki)‘;
????w=wij;
????wij=wij+lr*deltai*X‘+a*(wij-wij0);
????wij0=w;
end
epoch???????????????%顯示計(jì)算次數(shù)
%BP網(wǎng)絡(luò)的第二階段工作期(根據(jù)訓(xùn)練好的wki,wij和給定的輸入計(jì)算輸出)
X1=X;???????????????%給定輸入
%計(jì)算隱含層各神經(jīng)元輸出
NETi=wij*X1;
for?j=1:N
????for?i=1:q
????????Oi(ij)=2/(1+exp(-NETi(ij)))-1;
????end
end
%計(jì)算輸出層各神經(jīng)元的輸出
NETk=wki*Oi;
for?i=1:N
????for?k=1:L
????????Ok(ki)=2/(1+exp(-NETk(ki)))-1;
????end
end
Ok??????????????????%顯示網(wǎng)絡(luò)輸出層的輸出
?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件???????2145??2016-03-14?15:33??EMD\BPnetwork.m
?????文件???????1081??2016-03-07?20:16??EMD\emd.m
?????文件????????206??2016-03-07?20:17??EMD\FFTAnalysis.m
?????文件????????909??2016-03-07?20:15??EMD\findpeaks.m
?????文件????????235??2016-03-07?20:17??EMD\HilbertAnalysis.m
?????文件???????1688??2016-03-07?20:16??EMD\plot_hht.m
?????文件???????1180??2016-03-14?14:40??EMD\test.m
?????目錄??????????0??2016-03-24?15:08??EMD
-----------?---------??----------?-----??----
?????????????????7444????????????????????8
評(píng)論
共有 條評(píng)論