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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
- M, K! K- N- ?; q* b求高手解答！！！
我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
clc
0 A$ K# V/ A6 S6 g5 _0 Mclear
k=[1,2,3,4,5];%输入
for i=1:length(k)
8 Z6 _8 Z. Y% L9 c3 ]    T(i)=sum(k(1:i)); %输出
. C, }' |% W+ x/ N" oend
error_goal=0.01;%目标误差
" w+ R! D" w/ M8 Vmax_epoch=200;%最大训练次数
+ X7 l1 {0 R7 o1 X- b( b5 plr=0.5;%学习速率
/ Y1 g& I3 _2 _0 B) G%% 权值阀值初始化
/ x, H# W2 S$ ?- U  ^5 ]% N3 S5 OE=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
9 n' O  J  ^+ r7 X" d9 [. D1 m' OW11=rands(1);%权值初始化
W2=rands(2,1);
9 h3 z3 a$ e* U[W3,theta]=rands(1,2)
1 @- d6 \! j: u( fW21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
    DW110=0;
    DW210=0;
    DW220=0;
, S8 K$ K+ A; v* |) Q    DW310=0;
    DW320=0;
    Dtheta0=0;
2 v, }! m* V/ s% v% s) r4 z7 s    mc=0.5;
/ e( Q; E' k2 q) L3 `9 M%% 循环迭代
! h$ Z3 X: z* R; K5 i9 T6 ` for j=1:max_epoch
* O3 Z$ g- I+ r  h! k5 ^6 t    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
9 S  k' N/ L7 U3 C/ ?% O* K5 }for i=1:length(k)
   
   %% 网络输出计算
    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
) L( Z8 p# Q# P8 x2 ~; P    LC2=LB*W22    %LC层输出  线性函数
    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
        
   
    %% 权值修正
  b4 q$ c& T! n$ L. O    error=T(i)-ym  %计算误差
) k4 p3 i! |' S  \/ O& {  {4 _    E(i,j)=error;    %误差统计  
# D; I  k1 i* X    delta3=error
5 E# i7 x! g" H( `4 a    delta21=error*W31
    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
    %修改权值
  v" j+ K' n# A8 E0 Z    DW11=lr*delta1*k(i)
    DW21=lr*delta21*LB
" |) e9 t# }4 D9 Z( ~: F7 O    DW22=lr*delta22*LB
    DW31=lr*delta3*LC1
! y, Q, e, K; M- g" M    DW32=lr*delta3*LC2
! n9 u3 H: N8 j0 z6 R5 F4 Q    Dtheta=lr*delta3;
2 N6 @' V2 o" Q" ]/ T& a    W11=W11+DW11+mc*DW110
5 b2 d' P4 t. h    W21=W21+DW21+mc*DW210
    W22=W22+DW22+mc*DW220
( N% J- w2 Y' }1 z3 W. L    W31=W31+DW31+mc*DW310
    W32=W32+DW32+mc*DW320
* m9 X+ b! Y  o  I% G    theta=theta+Dtheta+mc*Dtheta0
    DW110=DW11
    DW210=DW21
/ j1 C6 M2 J$ `; T    DW220=DW22
    DW310=DW31
+ `* c5 z1 {6 R    DW320=DW32
    Dtheta0=Dtheta
end
if sumsqr(E(:,j))<error_goal
" ~3 ^6 X) o0 S. l2 f  u    disp('good')
! m8 A# a+ _& q+ Q) K* z    x=input('x=')
' R$ s' _& ~! P4 W+ S% t$ i    LB=1/(1+exp(-W11*x));   %LB层输出
& D, s" {  N# Y0 g7 _& G    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
* E( g0 K( U# C% h6 `   
    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
/ K& o% ]1 A& o9 n9 F+ k  e! n! Cbreak;
end
) t7 y4 P% Q+ }4 c6 Nend  
YY=0;
* w. z! q: O1 E: O3 T    for i=1:length(k)
; v9 ?& D4 N4 D1 N: G    LB=1/(1+exp(-W11*k(i)));   %LB层输出
    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
# A- N$ ~) |" V7 |4 `    LD=W31*LC1+W32*LC2;  %LD层输出
" F4 |$ `( e( v- x    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
   
    end
4 i1 D, A$ {. H- `. [" n" t     YY=YY(2:end);
0 J. d. o- l- n- N    YY

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
数学中国社区分享快乐！
: c! e1 Y  V' E! A# A; X  i, J