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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
) |. }% `9 w5 P! j, q+ ?; Y求高手解答！！！
! t# y9 @- i! d* V3 g. z我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
3 m4 o; Z/ b5 w+ V: [- Aclc
clear
k=[1,2,3,4,5];%输入
for i=1:length(k)
    T(i)=sum(k(1:i)); %输出
end
error_goal=0.01;%目标误差
. |1 U4 n- a- p/ N7 Mmax_epoch=200;%最大训练次数
lr=0.5;%学习速率
  r: m6 `( }: C; v) m%% 权值阀值初始化
" G: }" h! G* U# U+ A" t5 HE=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
W11=rands(1);%权值初始化
# ?. X" V/ N9 ^+ |# J* OW2=rands(2,1);
4 D* b. I7 \) W% }( C9 q2 ]% L[W3,theta]=rands(1,2)
/ }0 g, A! q/ Q% DW21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
# S; g4 S  F% Z7 H    DW110=0;
$ G1 c3 k* z. @; ?. _    DW210=0;
: q7 H( r3 `1 J  p. |; p  {6 c    DW220=0;
    DW310=0;
    DW320=0;
    Dtheta0=0;
4 N3 `- z! C0 M, |" W    mc=0.5;
6 `3 B" Q5 D+ x/ O, ^. ]%% 循环迭代
- c! W* _/ N$ p* m* d for j=1:max_epoch
* a3 P4 e2 x5 {$ r! ^# y, I0 [: [    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
$ P: ^* s# ~1 U  j9 r4 u8 Tfor i=1:length(k)
   
7 c+ m( t1 `8 z* I   %% 网络输出计算
) b. n7 D3 k  F; ?    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
4 w% C9 q+ B9 m+ H0 @' K    LC1=LB*W21          %LC层输出 线性函数
    LC2=LB*W22    %LC层输出  线性函数
1 X! u! ]$ Y; ~, z; y5 l4 d+ E+ g    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
- {2 d/ A* E( q5 X2 w        
   
    %% 权值修正
    error=T(i)-ym  %计算误差
    E(i,j)=error;    %误差统计  
    delta3=error
' {1 k) u: \6 L( b% M. v4 N+ G    delta21=error*W31
. d! S/ P2 w9 b% P4 Z    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
    %修改权值
    DW11=lr*delta1*k(i)
/ k7 U: P' c) o    DW21=lr*delta21*LB
    DW22=lr*delta22*LB
# H% |3 `5 L4 @. _8 W    DW31=lr*delta3*LC1
    DW32=lr*delta3*LC2
    Dtheta=lr*delta3;
2 |* x% z0 K( g  j3 }8 N6 W$ X( K    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
    theta=theta+Dtheta+mc*Dtheta0
    DW110=DW11
; {6 G! V# _7 a" C# P    DW210=DW21
2 _! b8 }- O' t+ S) {; M) g- p    DW220=DW22
& o0 v; j' k* M8 F; F* c) U    DW310=DW31
" l. R! o. U6 u    DW320=DW32
    Dtheta0=Dtheta
8 j7 K; h6 p9 S  a  _end
6 m; l7 y- u  h. b. J) e0 ~if sumsqr(E(:,j))<error_goal
9 @' G. x) Y) P    disp('good')
4 B  y# }1 F' P: d# i* p: h    x=input('x=')
    LB=1/(1+exp(-W11*x));   %LB层输出
  Q* k- L1 I1 b% `# h4 }    LC1=LB*W21;           %LC层输出
; Z2 ]" K' \) s- n: B/ o    LC2=LB*W22;    %LC层输出
5 i8 O  E7 `* V   
    LD=W31*LC1+W32*LC2  %LD层输出
& ~+ ^: z+ b! e* K# p    ym=LD-theta %网络输出值
& S8 I: ?' U! ?! G7 Fbreak;
end
end  
% z2 v! {# Z& B3 X# wYY=0;
    for i=1:length(k)
! H+ [: A' G0 R' Z) w" y0 k    LB=1/(1+exp(-W11*k(i)));   %LB层输出
$ o: [+ ^0 o4 e; [! F' m    LC1=LB*W21;           %LC层输出
  T+ m" J& ]- T    LC2=LB*W22;    %LC层输出
    LD=W31*LC1+W32*LC2;  %LD层输出
+ S& E+ ^! d+ }# I2 O; P# E# ?8 V    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
   
/ _8 R- q: Z) ^5 S+ `    end
+ f3 T* [3 ^& a/ s, z( ?     YY=YY(2:end);
# E7 `' X' i( L: V    YY

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
数学中国社区分享快乐！
4 E& N: `8 f) z) e/ Z' u; b