bp1.m

来自「该程序代码实现BP网络函数逼近」· M 代码 · 共 34 行

%Bp网络的训练
P=rands(10);%大家可以自己设
T=rands(10);
[R,Q]=size(P);
[S2,Q]=size(T);
S1=10;%可以自己设定，但必须满足网络要求
[W1,B1]=rands(S1,R);
[W2,B2]=rands(S2,S1);
A1=tansig(W1*P,B1);
A2=purelin(W2*A1,B2);
E=T-A2;
err_goal=0.2;
max_epoch=40;
SEE=sumsqr(E);
for epoch=1:max_epoch
    if SEE<err_goal
        epoch=epoch-1;
        break
    end
lr=0.5;
D2=deltalin(A2,E);
D1=deltatan(A1,D2,W2);
[dW1,dB1]=learnbp(P,D1,lr);
[dW2,dB2]=learnbp(A1,D2,lr);
W1=W1+dW1;
B1=B1+dB1;
W2=W2+dW2;
B2=B2+dB2;
end
W1
B1
W2
B2