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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
求高手解答！！！
$ r( c- O7 w; y( M我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
4 \" |7 L) N. W9 Uclc
clear
' ?& k7 ?! j+ ]' p3 W' hk=[1,2,3,4,5];%输入
! O' p& z2 U8 c! r4 S% x" vfor i=1:length(k)
( Y% [) L: t7 x- k) a    T(i)=sum(k(1:i)); %输出
8 x) @. q* M5 f# {  d7 Kend
error_goal=0.01;%目标误差
max_epoch=200;%最大训练次数
$ [0 T" `9 x+ u1 c* D4 Olr=0.5;%学习速率
3 c) W3 A8 \3 ]5 y. R3 Y- Y# Z%% 权值阀值初始化
E=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
W11=rands(1);%权值初始化
2 k! v5 t% X, [! B( pW2=rands(2,1);
0 X9 d$ d/ m  E/ I4 N[W3,theta]=rands(1,2)
1 V9 M) g! f; y+ }1 A, }9 SW21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
! w2 y4 Y8 H- P7 T& w% g6 r& H    DW110=0;
    DW210=0;
+ Z' S% s) C! i/ \$ j; R' a! }    DW220=0;
    DW310=0;
1 W$ s7 [3 I1 i+ n% ^    DW320=0;
* ]- g( N' \# c# |' L$ t0 m3 w0 _  _    Dtheta0=0;
5 g! U- ?" p8 m6 [* }) R6 L    mc=0.5;
& R9 f8 G4 {; Z4 k# f* x%% 循环迭代
1 n' ~& z# z. \# Q5 Q for j=1:max_epoch
    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
7 E' q" R) X; x, Q+ {2 T/ t& ^  Wfor i=1:length(k)
+ o  W) _1 s) H9 t   
   %% 网络输出计算
! y' j9 u; L/ \    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
! L6 L# z) y, v, E; O2 Z7 u    LC1=LB*W21          %LC层输出 线性函数
    LC2=LB*W22    %LC层输出  线性函数
) [& k( ^9 Q' E% l9 y    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
- e) }/ g. M8 J        
# Z; L) |4 C6 V2 ]   
4 v, R- i  B7 M9 E# s    %% 权值修正
    error=T(i)-ym  %计算误差
    E(i,j)=error;    %误差统计  
1 `. h5 D% D8 u( g$ p: Y) b    delta3=error
0 Z' W5 M# n3 _/ Q7 E# E    delta21=error*W31
    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
9 d( u  e( X$ Q" x+ U7 w; C, F    %修改权值
    DW11=lr*delta1*k(i)
    DW21=lr*delta21*LB
" |) X* ?, I: }8 {3 @7 F: f; W  ]    DW22=lr*delta22*LB
    DW31=lr*delta3*LC1
6 r5 a: l7 _8 J( P, ~! p    DW32=lr*delta3*LC2
; @, ]1 Z" U) _+ n" A0 R6 h    Dtheta=lr*delta3;
    W11=W11+DW11+mc*DW110
    W21=W21+DW21+mc*DW210
    W22=W22+DW22+mc*DW220
  f0 S& C& F# j6 K0 N2 T+ M: p    W31=W31+DW31+mc*DW310
$ s1 c6 j% @- D: N7 h9 k* ?! n# b. U    W32=W32+DW32+mc*DW320
    theta=theta+Dtheta+mc*Dtheta0
    DW110=DW11
/ Z# U! r) t) G4 P* l! |9 U    DW210=DW21
" x; V0 |" `1 i' L' ]    DW220=DW22
- N" z5 r. f* t. e# Z/ E) ^+ o    DW310=DW31
    DW320=DW32
) K9 b8 ^2 }+ E$ G4 e    Dtheta0=Dtheta
end
2 ^, ^$ ?5 T( d& f+ Z: K6 _if sumsqr(E(:,j))<error_goal
$ E# [( w) ^$ h: C+ D    disp('good')
    x=input('x=')
7 e$ l$ Y7 n; B0 y3 _; L: d    LB=1/(1+exp(-W11*x));   %LB层输出
    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
   
7 P0 j: \* L$ ]6 W: ]8 m# a    LD=W31*LC1+W32*LC2  %LD层输出
, c1 q) R* j4 v1 o    ym=LD-theta %网络输出值
; G6 ^+ m" @" n. V" u7 R' n8 a# ~break;
4 x! C- ?5 Z  O0 ]: Q( Z! k" send
end  
YY=0;
    for i=1:length(k)
    LB=1/(1+exp(-W11*k(i)));   %LB层输出
# z+ {( E. d' N. n    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
    LD=W31*LC1+W32*LC2;  %LD层输出
# j& ]1 G( \6 h0 d4 j, q    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
   
: b2 Y& n! q/ T. s. U) d+ m    end
/ e4 H# ?' o1 S% n0 i     YY=YY(2:end);
& x4 O3 N5 p' A& c1 O    YY

1 Q' e1 W* W- b7 e) _

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
6 l4 T1 k# B$ S; P4 _# ?数学中国社区分享快乐！