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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
( Y, R3 h5 L$ X( d1 L求高手解答！！！
我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
* P0 I9 ^! t& o% c% c- u但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
clc
# q  j9 C) E  }) ^+ |2 zclear
k=[1,2,3,4,5];%输入
4 e: L7 e; G" I: I/ }) ffor i=1:length(k)
0 M0 l$ R( t$ b4 s    T(i)=sum(k(1:i)); %输出
7 @' x( T$ ?8 t; J! A  _9 a5 i  m! `end
error_goal=0.01;%目标误差
" G8 y- _; C9 cmax_epoch=200;%最大训练次数
lr=0.5;%学习速率
%% 权值阀值初始化
; n" Q2 h1 J' N( R& q. F9 dE=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
W11=rands(1);%权值初始化
! E5 P: \' s+ uW2=rands(2,1);
[W3,theta]=rands(1,2)
# c4 s) ?( t& z3 jW21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
    DW110=0;
  b1 B; h. {1 D' C+ ?    DW210=0;
    DW220=0;
    DW310=0;
4 p( ~6 S0 }, F/ I; ?    DW320=0;
$ H# m% F0 z! y/ r    Dtheta0=0;
    mc=0.5;
%% 循环迭代
for j=1:max_epoch
: J  N) l% P  d4 p9 `# o    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
- A" F- f9 O/ a2 D# B! L3 wfor i=1:length(k)
   
   %% 网络输出计算
( q' J) M- Z: E' y    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
5 Z+ @  ^2 a  p, _/ U    LC2=LB*W22    %LC层输出  线性函数
    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
4 d6 U4 ], @+ u5 s) G  c        
   
    %% 权值修正
    error=T(i)-ym  %计算误差
/ q) c0 v( f9 ?- U  D    E(i,j)=error;    %误差统计  
    delta3=error
    delta21=error*W31
& f" ]0 f4 `9 U3 E1 e7 {    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
' E* O" q9 K8 z8 |2 {1 |# H    %修改权值
    DW11=lr*delta1*k(i)
) C( Q0 J# X! c% y    DW21=lr*delta21*LB
$ c7 T( U8 b1 U    DW22=lr*delta22*LB
7 j, Y# E/ B: o- X    DW31=lr*delta3*LC1
    DW32=lr*delta3*LC2
    Dtheta=lr*delta3;
7 G' N9 Z& n+ n0 |    W11=W11+DW11+mc*DW110
    W21=W21+DW21+mc*DW210
    W22=W22+DW22+mc*DW220
    W31=W31+DW31+mc*DW310
: @0 U* H$ @, B: _1 B* z; @    W32=W32+DW32+mc*DW320
    theta=theta+Dtheta+mc*Dtheta0
9 e  ~0 y3 z4 Z" K$ o    DW110=DW11
    DW210=DW21
    DW220=DW22
    DW310=DW31
, z2 v$ O8 |/ ?+ j, W1 ]    DW320=DW32
    Dtheta0=Dtheta
end
0 O. X+ K3 P3 `+ {0 v$ jif sumsqr(E(:,j))<error_goal
    disp('good')
" k% e  X  }9 a2 l7 ~    x=input('x=')
    LB=1/(1+exp(-W11*x));   %LB层输出
    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
   
; A; j$ I' b$ [. e/ M* e) U    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
break;
6 v& z8 U$ r# oend
- ^; E# t: {: o% }/ Pend  
YY=0;
$ O5 k8 t2 Z$ k3 `+ f6 P9 N    for i=1:length(k)
  B* G: i- U2 n. N4 T4 h    LB=1/(1+exp(-W11*k(i)));   %LB层输出
. I) T; p) e# {7 A/ M    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
    LD=W31*LC1+W32*LC2;  %LD层输出
    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
1 ~" R8 a9 b/ B   
) l  R6 x/ U# r. a- ~9 [    end
, }/ B* L' k$ D% u     YY=YY(2:end);
    YY

: L  M7 z  j, h

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
# |4 {* v3 X/ j  c, c2 t数学中国社区分享快乐！
! g4 T. ]; Z2 q: d4 o