1、BP網(wǎng)絡的MATLAB仿真,樣本,% bp1;例4.1 P=-1:0.1:1; P2=-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,r+);,P=-1:0.1:1;,P2=-1:0.1:1;,T=-0.96 -0.577 -0.0729 0.377 0.641 0.66 0.461 0.1336 . -0.

2、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,r+);,R,Q=size(P);S2,Q=size(T);S1=5; W1,B1=rands(S1,R); W2,B2=rands(S2,S1); b1=;b2=; b1=B1*ones(1,21); b2=B2*ones(1,21); a2=W2*tansig(W1*P2+b1)+b2; A2=purelin(a2); hold on plot(P,A2) hold off disp(按任一

3、鍵繼續(xù)) pause,R,Q=size(P);,S1=5;,S2,Q=size(T);,W1,B1=rands(S1,R);,rands(S1,R);可得到一個S1*R的矩陣，其元素為（1，1）均勻分布隨機數(shù)， w,b=rands(S1,R)可得到一個S1*R矩陣和S1*1列矢量b, 其元素為（1，1）均勻分布隨機數(shù)， W1,B1=rands(S1,R); W2,B2=rands(S2,S1); 得到 W1：S1R矩陣,輸入層的權值矩陣 W2：S2S1矩陣，輸出層的權值矩陣 及兩個列向量B1：S1*1 B2：S2*1,W2,B2=rands(S2,S1);,b1=;b2=;空矩陣,b1=B1*

4、ones(1,21);,b1=B1*ones(1,21); 把列向量B1的每一行擴展成121向量，數(shù)值重復，,b2=B2*ones(1,21);,a2=W2*tansig(W1*P2+b1)+b2;,tansig is a transfer function. Transfer functions calculate a layers output from its net input. 計算網(wǎng)絡輸出,A2=purelin(a2);,purelin Linear transfer function,網(wǎng)絡初始輸入輸出關系,創(chuàng)建并訓練網(wǎng)絡,net=newcf(minmax(P),5,1,tansig,purelin,traingd); %創(chuàng)建兩層前向回饋網(wǎng)絡 net.trainParam.epochs=7000; %初始化訓練次數(shù) net.trainParam.goal=9.5238e-004; % sse=0.02 %初始化誤差值 net.trainParam.lr = 0.15; net,tr=train(net,P,T); %