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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
求高手解答！！！
4 d+ }5 q$ ^% d我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
6 ~7 A& z- K/ C+ }但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
1 z8 S0 v6 p7 w7 a8 Zclc
+ e% d4 R4 P3 z" d$ H: J( s" Nclear
k=[1,2,3,4,5];%输入
for i=1:length(k)
    T(i)=sum(k(1:i)); %输出
' Y5 H4 D; d8 b7 p& k' U& _end
error_goal=0.01;%目标误差
: E) a. y# }% Y5 c( m6 M. Rmax_epoch=200;%最大训练次数
lr=0.5;%学习速率
+ i% M7 ?; ?1 K, R%% 权值阀值初始化
9 f* t( V; m) G  f! Y" ZE=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
W11=rands(1);%权值初始化
. f7 ?+ h2 h% C/ y5 pW2=rands(2,1);
[W3,theta]=rands(1,2)
8 O! j5 T$ L; D. q( |, q( j9 X/ z, \W21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
    DW110=0;
    DW210=0;
7 @3 ~7 d6 Y2 d: q0 @) S! `    DW220=0;
    DW310=0;
: o! W4 w8 B6 c1 P/ i, p) O# X    DW320=0;
" o/ b2 ?- b# t    Dtheta0=0;
    mc=0.5;
%% 循环迭代
for j=1:max_epoch
7 G9 ~7 {! n4 e: t5 e# w, y- e    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
for i=1:length(k)
5 w$ u! Q0 d2 ~1 h& f( b   
) r9 N2 K7 I; Z' S+ f   %% 网络输出计算
    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
    LC2=LB*W22    %LC层输出  线性函数
6 t2 R: p3 g9 N9 Q  y    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
/ b4 l/ M; F- H" B4 [    ym=LD-theta   %网络输出值
        
   
2 w, C) J) |8 s' y' F    %% 权值修正
    error=T(i)-ym  %计算误差
$ C* a5 k( S/ |    E(i,j)=error;    %误差统计  
    delta3=error
    delta21=error*W31
/ A7 ]$ ~8 }- x! {    delta22=error*W32
& T5 {( x/ I$ g4 a5 p6 F    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
3 o# v% e" m3 }! C5 g    %修改权值
: |. w- J* O+ [    DW11=lr*delta1*k(i)
( V1 l1 x1 ]. p: W    DW21=lr*delta21*LB
    DW22=lr*delta22*LB
8 X/ h3 r/ p9 W3 b: ?    DW31=lr*delta3*LC1
    DW32=lr*delta3*LC2
6 r, Q) q. A- a$ k0 E- t1 S# z: O    Dtheta=lr*delta3;
    W11=W11+DW11+mc*DW110
( s3 F. Q! ~, i% y* ^    W21=W21+DW21+mc*DW210
    W22=W22+DW22+mc*DW220
    W31=W31+DW31+mc*DW310
/ s% w; ?# t2 f" T2 ~    W32=W32+DW32+mc*DW320
    theta=theta+Dtheta+mc*Dtheta0
, W5 Q. ?& D" C    DW110=DW11
% e$ Y7 Y+ w; V3 n6 D    DW210=DW21
' F4 B: `1 N/ }' p! J8 N- }2 n4 Z! S    DW220=DW22
) F, Z0 B" O3 q0 f- L! R2 b/ x    DW310=DW31
    DW320=DW32
, i+ l3 i+ }: r, P# ]% z    Dtheta0=Dtheta
% p! v1 B$ B3 G3 T1 L$ ^1 Pend
% u# z+ i1 W& E7 w1 m; B* Fif sumsqr(E(:,j))<error_goal
  V4 @& ^3 ~7 t  [$ l    disp('good')
    x=input('x=')
+ ?# n. e2 f; D- @; J& Y& R; V    LB=1/(1+exp(-W11*x));   %LB层输出
) x, `) `; }, H+ K8 P+ g% r    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
- P8 F5 ]) h8 ?7 C   
    LD=W31*LC1+W32*LC2  %LD层输出
/ q/ O# ]$ B  i% R9 T7 Y$ g3 ~% o    ym=LD-theta %网络输出值
; X: c6 A: f" `break;
: }, _; I! z# h' k& z4 n8 }/ wend
1 K' q3 l# D+ P. {! @$ Qend  
YY=0;
+ B5 d( b$ c# Q" c    for i=1:length(k)
    LB=1/(1+exp(-W11*k(i)));   %LB层输出
, t* l! b6 ~2 x/ t7 L- q+ X6 e    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
  Z8 N/ ^3 o1 h5 Q9 Q$ u: g    LD=W31*LC1+W32*LC2;  %LD层输出
    ym=LD-theta;   %网络输出值
+ J3 c' X, ^: k/ c4 H, D0 D    YY=[YY,ym];
   
8 H  w0 ?( s5 ^' S% k1 v2 p    end
4 U% F! H4 w0 `6 N4 s+ }$ f     YY=YY(2:end);
0 g3 T7 r: t; q* s, f    YY

. l  x) d$ ]+ Z; N# T

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
数学中国社区分享快乐！