% ***用于函數(shù)逼近的BP網(wǎng)絡(luò)設(shè)計(jì)***
% 設(shè)計(jì)一兩層BP網(wǎng)絡(luò)，其網(wǎng)絡(luò)的隱層各神經(jīng)元的激活函數(shù)為雙曲正切型，
% 輸出層各神經(jīng)元的激活函數(shù)為線性函數(shù)，隱層含有5個(gè)神經(jīng)元，
% 并且有如下21組輸入矢量和相對(duì)應(yīng)的目標(biāo)矢量：
p=-1:0.1:1;
t=[-0.96 -0.577 -0.0729 0.377 0.641 0.66 0.461 0.1336 -0.201 -0.434 -0.5 -0.393 -0.1647 0.0988 0.3072 0.396 0.3449 0.1816 -0.0312 -0.2183 -0.3201];
plot(p,t)
pause
%首先對(duì)網(wǎng)絡(luò)進(jìn)行初始化
[R,Q]=size(p);
[S2,Q]=size(t);
S1=5;
[w1,b1]=rands(S1,R);
[w2,b2]=rands(S2,S1);
%用輸入矢量p計(jì)算網(wǎng)絡(luò)的輸出
A2=purelin(w2*tansig(w1*p,b1),b2);
plot(p,A2)
pause
%下面定義網(wǎng)絡(luò)訓(xùn)練參數(shù)
disp_fqre=10;max_epoch=18000;err_goal=0.01;lr=0.01;
TP=[disp_fqre max_epoch err_goal lr];
[w1,b1,w2,b2,epochs,errors]=trainbp(w1,b1,'tansig',w2,b2,'purelin',p,t,TP)
%運(yùn)行上述程序后，可以返回訓(xùn)練后的權(quán)值，訓(xùn)練次數(shù)和誤差平方和
ploterr(errors)