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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
求高手解答！！！
$ I1 }: H1 \3 G# y  t' B我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
clc
clear
& H4 ~2 |2 [5 o9 F7 @k=[1,2,3,4,5];%输入
5 P$ n# b7 K3 A% _5 w+ rfor i=1:length(k)
0 q) T$ f& x2 H( j3 t1 l% [' G    T(i)=sum(k(1:i)); %输出
0 w: z* p% l# @3 v! Wend
error_goal=0.01;%目标误差
& B" [& ~4 ]% |( D: M! cmax_epoch=200;%最大训练次数
3 K1 E$ M9 E' R: L' z$ r! R" {lr=0.5;%学习速率
! N' M9 D! R$ n0 L0 l%% 权值阀值初始化
3 F! a) G3 U; O. @7 r9 v8 j( L, bE=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
+ G, e" n! @' r. jW11=rands(1);%权值初始化
( G1 r( _) v- T3 G2 ?W2=rands(2,1);
[W3,theta]=rands(1,2)
' m4 I. y" g" `" x3 IW21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
    DW110=0;
    DW210=0;
( P2 T* E! j  ~0 N+ T    DW220=0;
    DW310=0;
* E' h6 t- [' Z+ }/ _9 S& t# @# H    DW320=0;
8 }' ]7 V9 B1 G& L; \/ u9 ~. p    Dtheta0=0;
    mc=0.5;
%% 循环迭代
. V9 ~4 }2 M: {/ ]/ k$ ~' \ for j=1:max_epoch
    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
for i=1:length(k)
5 j7 b0 L7 W# q. ]- J6 p   
+ V: h  X- r2 T. W: |2 D   %% 网络输出计算
, }3 g8 H" m4 P    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
    LC2=LB*W22    %LC层输出  线性函数
    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
1 g: X  ]. u- R8 Q        
% k/ y/ u# r- y7 D* U9 W, H9 I   
! A& o: a5 o. k' ]    %% 权值修正
    error=T(i)-ym  %计算误差
, E$ ^; |0 \3 R( J% X    E(i,j)=error;    %误差统计  
( w3 X& f. B5 h' }  R8 {    delta3=error
) A# @3 u& A8 c    delta21=error*W31
9 n& E6 l4 H: \  }    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
    %修改权值
1 x- B. V/ G0 v! O    DW11=lr*delta1*k(i)
    DW21=lr*delta21*LB
+ e  P9 c' `' c) y% r    DW22=lr*delta22*LB
  B5 G* Y  z5 X1 u! `4 s    DW31=lr*delta3*LC1
    DW32=lr*delta3*LC2
    Dtheta=lr*delta3;
2 i; D5 r1 s9 s" y* d9 q1 c, O# D    W11=W11+DW11+mc*DW110
    W21=W21+DW21+mc*DW210
    W22=W22+DW22+mc*DW220
    W31=W31+DW31+mc*DW310
    W32=W32+DW32+mc*DW320
! v3 R1 ]+ p) I4 P    theta=theta+Dtheta+mc*Dtheta0
4 m5 t, Q$ [8 \  l5 t- o+ _    DW110=DW11
; {6 a  h6 z% \# q' }! u    DW210=DW21
& _+ R, ?" }& p. w    DW220=DW22
  K& [% F. F  _* p    DW310=DW31
/ }5 V, M) ]! n( Z    DW320=DW32
& `& c6 h6 f  h$ ?    Dtheta0=Dtheta
end
if sumsqr(E(:,j))<error_goal
    disp('good')
    x=input('x=')
' n* u* O9 I: i; |' n    LB=1/(1+exp(-W11*x));   %LB层输出
    LC1=LB*W21;           %LC层输出
3 O; ~4 Q" K) Z# _. O3 o5 R+ X    LC2=LB*W22;    %LC层输出
6 B6 C& j1 {1 s9 ~1 C   
; |  a) K% S, Q: b( e9 R    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
: ?3 N  b! g  \4 e; gbreak;
) p; r  Z, o+ C9 D# wend
# n0 x! F& {4 s1 _* Q0 L0 u$ Rend  
( h  Q8 M; i0 A/ _6 \% @# ]* CYY=0;
    for i=1:length(k)
( x$ ]2 M( R( y( B; d( `    LB=1/(1+exp(-W11*k(i)));   %LB层输出
    LC1=LB*W21;           %LC层输出
! c* y. }& F/ j; P$ \8 t9 G2 M) r    LC2=LB*W22;    %LC层输出
    LD=W31*LC1+W32*LC2;  %LD层输出
    ym=LD-theta;   %网络输出值
, Y5 P  B  n- G" ~2 A    YY=[YY,ym];
/ }! J5 Z* [; o! G0 l, f   
2 C+ u3 E4 A$ t    end
     YY=YY(2:end);
    YY

1 l* A; H5 }% b( a+ ^

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
) w3 n- U9 e6 J8 g, a' [2 h数学中国社区分享快乐！