mybp2.m

经典的神经网络matlab程序。

NNTWARN OFF
P=0:0.25:2*pi;
T=sin(3*P);

P2=0:0.025:2*pi;
%P2=-1:0.025:1;

[R,Q]=size(P);
%[S2,Q]=size(T);
[S3,Q]=size(T);
S1=8;
S2=10;
[W1,B1]=rands(S1,R);
%[W1,B1]=nwtan(S1,R);

[W2,B2]=rands(S2,S1);

[W3,B3]=rands(S3,S2);

%A2=purelin(W2*tansig(W1*P2,B1),B2);

disp_freq=10;
max_epoch=32000;
err_goal=0.05;
lr=0.01;
lr_inc=1.05;
lr_dec=0.7;
TP=[disp_freq max_epoch err_goal lr];
%TP=[disp_freq max_epoch err_goal lr lr_inc lr_dec err_ratio];
[W1,B1,W2,B2,W3,B3,epochs,TE]=trainbp(W1,B1,'tansig',W2,B2,'tansig',W3,B3,'purelin',P,T,TP)
%trainbp(W1,B1,'tansig',W2,B2,'purelin',P,T,TP)

figure;
A2=purelin(W3*tansig(W2*tansig(W1*P2,B1),B2),B3);
plot(P2,sin(3*P2),':',P2,A2,'-g');