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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
求高手解答！！！
我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
  i- a3 G$ U3 {9 s! v但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
. G3 \* d  A4 n% u4 mclc
clear
$ i) o( m. j: k5 [k=[1,2,3,4,5];%输入
% ]! l* w5 L  o) D2 r8 Ifor i=1:length(k)
* E( G2 I2 _& I+ |. ~    T(i)=sum(k(1:i)); %输出
) K/ T# k$ x4 t4 Zend
: Q! I4 a. Y: V" Aerror_goal=0.01;%目标误差
max_epoch=200;%最大训练次数
1 f, u  v( H# ^+ o0 H( \lr=0.5;%学习速率
%% 权值阀值初始化
E=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
6 F0 r9 l  t; M$ h) S2 BW11=rands(1);%权值初始化
W2=rands(2,1);
[W3,theta]=rands(1,2)
W21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
7 K, }! T* ]: s% @2 K    DW110=0;
    DW210=0;
# R$ d, m; e% ]; f$ V    DW220=0;
0 @# o4 W  J" d4 v    DW310=0;
( T! a+ o% R; V5 ^! T; D0 D    DW320=0;
- W: [/ U) K/ _& t6 D# |    Dtheta0=0;
    mc=0.5;
# r9 D$ C5 F+ H1 G5 W" Y%% 循环迭代
, O+ c" n8 y% L7 O6 L7 J, V for j=1:max_epoch
6 l' S" V' `, l) M$ J    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
6 v7 n* T3 I) Z5 x: J1 H2 q" k0 Wfor i=1:length(k)
   
- J" X; w) R. H) Z2 v9 h   %% 网络输出计算
7 }/ l4 ]3 K  m& u; [    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
) c: N% G( c: t# ?% C0 c; h4 [    LC2=LB*W22    %LC层输出  线性函数
    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
        
   
+ y7 E7 c& O, G% o" o* v, g    %% 权值修正
5 ?  l1 c& \$ d  Y: K1 B    error=T(i)-ym  %计算误差
    E(i,j)=error;    %误差统计  
2 W, _' N( x" i+ l    delta3=error
( a6 N% w- z7 J1 o4 A2 F( g    delta21=error*W31
2 V1 k9 Y/ v( W8 y5 C" X    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
    %修改权值
& ?! b" l3 ]  z! ?/ t    DW11=lr*delta1*k(i)
0 S* z4 c& H2 j4 T6 P# D8 L) _    DW21=lr*delta21*LB
    DW22=lr*delta22*LB
    DW31=lr*delta3*LC1
( t# a; S7 l7 P" P$ b& z    DW32=lr*delta3*LC2
" k& h$ y/ x6 h( N2 B    Dtheta=lr*delta3;
    W11=W11+DW11+mc*DW110
2 y. _* ~3 {: t% b: u    W21=W21+DW21+mc*DW210
" T' ~7 x) @" ?; }' [, v    W22=W22+DW22+mc*DW220
    W31=W31+DW31+mc*DW310
8 S$ g$ Q8 a  F$ @2 S    W32=W32+DW32+mc*DW320
7 S; K. H2 d+ U3 Y9 g& t' u    theta=theta+Dtheta+mc*Dtheta0
+ o- d* K. V) G( R! F8 i9 K  Y    DW110=DW11
" k  u" r/ C0 Y- o    DW210=DW21
0 o, ]' n# _; W, a) [7 _    DW220=DW22
& Y  B5 |; o+ u    DW310=DW31
    DW320=DW32
    Dtheta0=Dtheta
: Q$ q9 L9 i! ?& [7 mend
if sumsqr(E(:,j))<error_goal
    disp('good')
9 F/ f) z3 u$ ]# i/ Z8 S% T' J0 Z( R    x=input('x=')
7 L+ g! y) Z" t& ^1 i) O9 q    LB=1/(1+exp(-W11*x));   %LB层输出
    LC1=LB*W21;           %LC层输出
1 s# O8 Z" x* z    LC2=LB*W22;    %LC层输出
   
    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
break;
& _& J; i4 v4 n5 T0 x: Wend
8 }# C! p( ?% F8 [" G9 U" P7 hend  
, w7 \( N7 ^9 o) \6 f4 D) i% Q1 LYY=0;
    for i=1:length(k)
    LB=1/(1+exp(-W11*k(i)));   %LB层输出
    LC1=LB*W21;           %LC层输出
1 \( ]! ~4 Y  |  D& C: p0 z+ A  ?    LC2=LB*W22;    %LC层输出
+ _" m! u7 O: f7 M1 U1 A    LD=W31*LC1+W32*LC2;  %LD层输出
/ z) t; X( P6 l7 C    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
   
    end
     YY=YY(2:end);
, M& U! q. V" b9 n' h, @* W3 @+ V: w    YY
4 e1 r% E- w% |9 d

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
+ P5 G6 W9 T9 \7 h- U! w* D0 y数学中国社区分享快乐！
  A4 Q  {/ j7 }4 X7 b) R" T8 q