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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
求高手解答！！！
2 e" G9 p- g7 a0 W; d" r我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
& y8 t% e, E( I9 Y但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
clc
clear
; t: x4 j7 \8 ]( @7 p& Rk=[1,2,3,4,5];%输入
, X1 z5 k8 f  z9 p. z+ s0 W# ofor i=1:length(k)
5 _$ }% G" }; q8 x' S    T(i)=sum(k(1:i)); %输出
, A1 e, G) H$ r0 A& zend
error_goal=0.01;%目标误差
max_epoch=200;%最大训练次数
8 R  I5 `! T9 qlr=0.5;%学习速率
1 u! Y: y( H: z%% 权值阀值初始化
E=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
W11=rands(1);%权值初始化
W2=rands(2,1);
- W7 A( U% T4 E* L$ U, A6 B6 _* M[W3,theta]=rands(1,2)
W21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
    DW110=0;
    DW210=0;
    DW220=0;
    DW310=0;
    DW320=0;
3 L3 r, y$ j, _% K    Dtheta0=0;
/ ], G% r- I/ `) T9 D! O& D2 y    mc=0.5;
%% 循环迭代
for j=1:max_epoch
/ c' L! L( P4 B; K/ Z7 o  C! F    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
for i=1:length(k)
   
   %% 网络输出计算
    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
: d) T1 G( f1 H8 o$ ~/ U5 n8 `4 i7 S    LC2=LB*W22    %LC层输出  线性函数
    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
        
   
: f: b/ J  J% Y    %% 权值修正
    error=T(i)-ym  %计算误差
' h% Q4 S( f7 d& }0 k7 n    E(i,j)=error;    %误差统计  
    delta3=error
    delta21=error*W31
    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
    %修改权值
4 A+ P% \' Q8 q6 H7 h* d: e    DW11=lr*delta1*k(i)
2 X! U9 K% l% Z2 G; V2 c    DW21=lr*delta21*LB
6 D2 i' B0 T) {8 R3 k- O1 T    DW22=lr*delta22*LB
    DW31=lr*delta3*LC1
    DW32=lr*delta3*LC2
. k( [0 x, l: K  g# \    Dtheta=lr*delta3;
" l+ R; I! G: b. E    W11=W11+DW11+mc*DW110
; R5 K1 H9 @8 R4 G; ?    W21=W21+DW21+mc*DW210
    W22=W22+DW22+mc*DW220
$ k0 |6 @, \) N) A% c    W31=W31+DW31+mc*DW310
+ K% B, M, l- @# o8 D    W32=W32+DW32+mc*DW320
    theta=theta+Dtheta+mc*Dtheta0
    DW110=DW11
    DW210=DW21
- `8 t) g0 Q# f( ]3 R/ f- {    DW220=DW22
: y1 b: X! u% Q) q    DW310=DW31
    DW320=DW32
    Dtheta0=Dtheta
$ D6 z! c8 q5 }) Y$ \3 E6 \' @" bend
" b0 _: y2 h+ O4 n# O6 m& T7 [if sumsqr(E(:,j))<error_goal
    disp('good')
    x=input('x=')
1 a( |8 b4 A4 i3 B, d: ]; s    LB=1/(1+exp(-W11*x));   %LB层输出
    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
   
    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
5 y; I& {8 H5 y, j+ ybreak;
end
end  
YY=0;
- T3 a9 M" H0 V3 U8 p    for i=1:length(k)
, t5 q) q) I. y' ~  Q2 ]3 r4 J; m    LB=1/(1+exp(-W11*k(i)));   %LB层输出
    LC1=LB*W21;           %LC层输出
( ~7 l- h+ a1 _& ?9 L7 [( a) A    LC2=LB*W22;    %LC层输出
3 A- S' U0 N+ m6 s    LD=W31*LC1+W32*LC2;  %LD层输出
    ym=LD-theta;   %网络输出值
+ e, }  b+ S2 y9 @# E- V6 i    YY=[YY,ym];
1 m$ B( S9 E( q. x! R  d   
    end
7 O, R2 t/ i; P$ P; T) p     YY=YY(2:end);
    YY

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
% D1 G( a  \; r5 s. b$ J. _数学中国社区分享快乐！
; `1 M* P( [: [0 P/ r. [  m- u( m