BP网络底层程序，求高手解答~~~~~

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
6 i. T. ~5 Z  M$ p求高手解答！！！
我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
4 s& n* f6 X* v! E. u7 f但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
clc
5 e6 Y0 `. c' I+ ]% v8 ]. Gclear
k=[1,2,3,4,5];%输入
for i=1:length(k)
" K2 b& q' U" s, D5 ?    T(i)=sum(k(1:i)); %输出
end
error_goal=0.01;%目标误差
% L7 f4 g7 f" y; k6 y) t* amax_epoch=200;%最大训练次数
lr=0.5;%学习速率
%% 权值阀值初始化
/ }$ J) k6 J6 m" V4 rE=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
7 x0 s4 i1 a; y) f( C& E$ kW11=rands(1);%权值初始化
W2=rands(2,1);
" e6 A1 t* O- l# A[W3,theta]=rands(1,2)
W21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
    DW110=0;
. ?1 n; }4 Q, L' I% L; f    DW210=0;
# K/ }4 ~6 D5 Z0 V9 n7 u    DW220=0;
    DW310=0;
  c3 ^1 o5 {0 f4 ?& L* B    DW320=0;
    Dtheta0=0;
    mc=0.5;
%% 循环迭代
for j=1:max_epoch
    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
for i=1:length(k)
6 \) i1 x' H! r9 R   
0 m, \& a. w# g% G* |; I   %% 网络输出计算
/ a( ]; i( Q! ?. N# A+ f    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
    LC2=LB*W22    %LC层输出  线性函数
    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
! |4 O0 H1 o+ `, W  [3 J; N! p    ym=LD-theta   %网络输出值
        
# i; ]( w& z; r, [& a% W, [, _% y   
" X; \3 `# c+ h    %% 权值修正
9 B; \3 \8 ^9 @% V8 _2 g    error=T(i)-ym  %计算误差
5 g% r. g, v( \# t: \    E(i,j)=error;    %误差统计  
    delta3=error
    delta21=error*W31
* L( ~5 s& W% b! j3 N$ M    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
$ [( o0 y  ~! W3 x    %修改权值
    DW11=lr*delta1*k(i)
& h0 c% L4 {+ @3 x+ X6 S    DW21=lr*delta21*LB
    DW22=lr*delta22*LB
    DW31=lr*delta3*LC1
' I7 e# r6 L  I+ E  B    DW32=lr*delta3*LC2
    Dtheta=lr*delta3;
0 \, u; b1 i) I1 {9 \    W11=W11+DW11+mc*DW110
    W21=W21+DW21+mc*DW210
9 Y9 O1 S, g$ a8 I* w, t. d    W22=W22+DW22+mc*DW220
    W31=W31+DW31+mc*DW310
5 J5 T3 T; e5 A, L5 N3 \! q    W32=W32+DW32+mc*DW320
    theta=theta+Dtheta+mc*Dtheta0
    DW110=DW11
9 }! O, `, E( b# S2 `    DW210=DW21
    DW220=DW22
    DW310=DW31
    DW320=DW32
    Dtheta0=Dtheta
end
( d# l, H* L8 R$ R+ ~9 zif sumsqr(E(:,j))<error_goal
    disp('good')
    x=input('x=')
    LB=1/(1+exp(-W11*x));   %LB层输出
. K; d% u. X7 A    LC1=LB*W21;           %LC层输出
; Y+ n3 }1 M2 e. P# p' n    LC2=LB*W22;    %LC层输出
   
    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
break;
( R$ F' O' M. W: Q4 i6 v  F# z' P- ^+ jend
. o, \4 m/ w" _3 `$ m$ X. _$ Vend  
( h" P0 y' \- m0 W1 dYY=0;
! A) \% Q! A& v2 K# a0 L, O    for i=1:length(k)
; h3 k0 Q7 H; s0 n    LB=1/(1+exp(-W11*k(i)));   %LB层输出
7 v# K: M9 c( \7 E3 `    LC1=LB*W21;           %LC层输出
+ \, s0 t% e8 f: a- [1 f    LC2=LB*W22;    %LC层输出
    LD=W31*LC1+W32*LC2;  %LD层输出
    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
/ p; a+ y: D' T% x. W   
    end
     YY=YY(2:end);
    YY

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
数学中国社区分享快乐！