训练bp神经网络连接权值的源代码.txt

var matlab:variant //通过"变体"调用接口是比较低效的,但很方便 begin //变体这种结构

Ir=0.05; %Ir为学习速率
err_goal=0.001;%期望最小误差值
max_epoch=10000;
X=[0.75 1 0 1;0.25 0 0.55 0;0 0 0.45 0;0 0 0 0;0 0 0 0;1 1 1 1;0.5 0 0.5 0.5;0.5 0 0.5 0.5;0 1 0 0;1 1 1 1;0 0 0 0;0 0 0 0;0 0 0 0;0 0 0 0;1 1 1 1];
T=[0 0 0 0;0 1 1 1;0 0 0 0;1 0 0 0];%提供4组15输入4输出训练集和目标集
[M,N]=size(X);q=10;[L,N]=size(T);
Wij=rand(q,M);
Wki=rand(L,q);
b1=zeros(q,1);b2=zeros(L,1);—随机给定隐含层、输出层偏值
for epoch=1:max_epoch
Oi=tansig(Wij*X,b1);
Ok=purelin(Wki*Oi,b2);
E=T-Ok;
deltak=deltalin(Ok,E);%计算输出层的delta
deltai=deltatan(Oi,deltak,Wki);%计算隐含层的deita
[dWki,db2]=learnbp(Oi,deltak,Ir);%调整输出层加权系数
Wki=Wki+dWki;b2=b2+db2;
[dWij,db1]=learnbp(X,deltai,Ir);
Wij=Wij+dWij;b1=b1+db1;
SSE=sumsqr(T-purelin(Wki*tansig(Wij*X,b1),b2));
if(SSE<err_goal) break;end
end
epoch %显示计算次数
X1=X;
Oi=tansig(Wij*X1,b1);%各隐含层输出
Ok=purelin(Wki*Oi,b2);%显示网络输出层的输出