%嚴(yán)格按照BP網(wǎng)絡(luò)計(jì)算公式來設(shè)計(jì)的一個matlab程序，對BP網(wǎng)絡(luò)進(jìn)行了設(shè)計(jì)
%yyy，即在o（k）計(jì)算公式時，當(dāng)網(wǎng)絡(luò)進(jìn)入平坦區(qū)時（<0.0001）學(xué)習(xí)率加大，出來后學(xué)習(xí)率又還原
%v（ij）=v（ij）+deltv（ij）+a*dv（ij）；動量項(xiàng)
clear?all
inputnums=1;%輸入層節(jié)點(diǎn)
outputnums=1;%輸出層節(jié)點(diǎn)
hidenums=7;%隱層節(jié)點(diǎn)
maxcount=50000;%最大迭代次數(shù)
samplenum=19;%一個計(jì)數(shù)器，無意義
precision=0.00001;%預(yù)設(shè)精度
yyy=1.3;%yyy是幫助網(wǎng)絡(luò)加速走出平坦區(qū)
alpha=0.02;%學(xué)習(xí)率設(shè)定值
a=0.9;%BP優(yōu)化算法的一個設(shè)定值，對上組訓(xùn)練的調(diào)整值按比例修改
error=zeros（1maxcount+1）;%error數(shù)組初始化；目的是預(yù)分配內(nèi)存空間
errorp=zeros（1samplenum）;%同上
v=rand（inputnumshidenums）;%1*10；v初始化為一個1*10的隨機(jī)歸一矩陣；v表輸入層到隱層的權(quán)值
deltv=zeros（inputnumshidenums）;%1*10；內(nèi)存空間預(yù)分配
dv=zeros（inputnumshidenums）;%1*10;

w=rand（hidenumsoutputnums）;%10*1
deltw=zeros（hidenumsoutputnums）;%10*1
dw=zeros（hidenumsoutputnums）;%10*1

samplelist=1:0.5:10;%樣本輸入值
expectlist=samplelist.^（-1）;%期望輸出值
%expectlist=exp（-samplelist）;
[samplelistminpmaxp]=premnmx（samplelist）;
%[expectlistmintmaxt]=premnmx（expectlist）;
count=1;
while（count<=maxcount）%結(jié)束條件1迭代20000次
????c=1;
????while（c<=samplenum）
????????for?k=1:outputnums
????????????d（k）=expectlist（c）;%獲得期望輸出的向量
????????end
????????for?i=1:inputnums
????????????x（i）=samplelist（c）;%獲得輸入的向量數(shù)據(jù)
????????end
????????
????????%forward（）;
????????for?j=1:hidenums
????????????net=0.0;
????????????for?i=1:inputnums
????????????????net=net+x（i）*v（ij）;%輸入層到隱層的加權(quán)和
????????????end
????????y（j）=1/（1+exp（-net））;%輸出層處理f（x）=1/（1+exp（-x）），單極性sigmiod函數(shù)
???????????????%y（j）=1/net;%輸出層處理f（x）=1/（1+exp（-x）），單極性sigmiod函數(shù)
????????end
????????for?k=1:outputnums
????????????net=0.0;
????????????for?j=1:hidenums
????????????????net=net+y（j）*w（jk）;
????????????end
????????????if?count>=2&&error（count）-error（count-1）<=0.0001
???????????????o（k）=1/（1+exp（-net）/yyy）;%平坦區(qū)加大學(xué)習(xí)率
?????????????????%o（k）=1/（net/yyy）;%平坦區(qū)加大學(xué)習(xí)率
???????????else?o（k）=1/（1+exp（-net））;%同上
????????????????%else?o（k）=1/net;%同上
????????????end
????????end
????????
????????%BpError（c）反饋/修改
????????errortmp=0.0;
????????for?k=1:outputnums
????????????errortmp=errortmp+（d（k）-o（k））^2;%第一組訓(xùn)練后的誤差計(jì)算
????????end
????????errorp（c）=0.5*errortmp;%誤差E=∑（d（k）-o（k））^2?*?1/2?
????%end
????%backward（）;
????for?k=1:outputnums
????????yitao（k）=（d（k）-o（k））*o（k）*（1-o（k））;%輸出層誤差偏導(dǎo)
????end
????for?j=1:hidenums
????????tem=0.0;
????????for?k=1:outputnums
????????????tem=tem+yitao（k）*w（jk）;%為了求隱層偏導(dǎo)，而計(jì)算的求和
????????end
????????yitay（j）=tem*y（j）*（1-y（j））;%隱層偏導(dǎo)
????end
????
????%調(diào)整各層權(quán)值
????for?j=1:hidenums
????????for?k=1:outputnums
????????????deltw（jk）=alpha*yitao（k）*y（j）;%權(quán)值w的調(diào)整量deltw（已乘學(xué)習(xí)率）
????????????w（jk）=w（jk）+deltw（jk）+a*dw（jk）;%權(quán)值調(diào)整，這里的dw=deltw（t-1），實(shí)際是對BP算法的一個改進(jìn)措施--增加動量項(xiàng)目的是提高訓(xùn)練速度
????????????dw（jk）=deltw（jk）;
????????end
????end
????for?i=1:inputnums
????????for?j=1:hidenums
????????????deltv（ij）=alpha*yitay（j）*x（i）;%同上deltw
????????????v（ij）=v（ij）+deltv（ij）+a*dv（ij）;
????????????dv（ij）=deltv（ij）;
????????end
????end
????c=c+1;
????end%第二個while結(jié)束；表示一次BP訓(xùn)練結(jié)束
????double?tmp;
????tmp=0.0;
????for?i=1:samplenum
????????tmp=tmp+error

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

?????文件?????751104??2008-01-05?16:01??神經(jīng)網(wǎng)絡(luò)程序1\M.doc

?????文件???????4419??2008-01-04?16:57??神經(jīng)網(wǎng)絡(luò)程序1\xfenzhiyibp.m

?????文件????????225??2007-12-20?09:31??神經(jīng)網(wǎng)絡(luò)程序1\人工神經(jīng)網(wǎng)絡(luò)要求.txt

?????目錄??????????0??2009-10-01?14:23??神經(jīng)網(wǎng)絡(luò)程序1

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

???????????????755748????????????????????4