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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
6 R- g' N, d( @6 |求高手解答！！！
  v* T1 j6 _) ]/ C$ u, u) c( }& Q我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
+ ?) a  _6 `/ k' B但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
- E8 q0 B* v2 q7 ^% ~* N4 U# jclc
  l% Z, ?# I& t3 J4 s: @. V1 y6 ~- ]8 o* Tclear
9 w* ?' J0 J3 g$ p" N5 [k=[1,2,3,4,5];%输入
; [9 w; r2 |9 A8 \for i=1:length(k)
    T(i)=sum(k(1:i)); %输出
end
2 p# O) m6 `8 N6 _$ }: yerror_goal=0.01;%目标误差
6 ^4 i8 J" s5 M* d: `max_epoch=200;%最大训练次数
2 x2 |: k+ S0 _lr=0.5;%学习速率
8 G% H* y" o. Z( O%% 权值阀值初始化
E=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
W11=rands(1);%权值初始化
W2=rands(2,1);
[W3,theta]=rands(1,2)
W21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
' b# v7 a+ w6 Q8 L! K" |    DW110=0;
: g* ^, }0 m  n( t8 g    DW210=0;
! Q8 X8 W; g) J3 U    DW220=0;
    DW310=0;
    DW320=0;
" x& `0 D- R( X' h' Z    Dtheta0=0;
! F2 d. E' L" j) o  C, B' U1 d3 \    mc=0.5;
%% 循环迭代
' f$ Y8 r) z4 W% A+ H1 H- { for j=1:max_epoch
' Q2 n% [, x/ D& K  Q2 n    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
4 C# F) f) H- c2 v. j9 mfor i=1:length(k)
/ j1 ]' E/ ?4 x" S   
   %% 网络输出计算
- F( R. b( k' _* s) I( b6 h' P' H    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
- }: |" Y$ H0 u    LC2=LB*W22    %LC层输出  线性函数
    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
        
   
    %% 权值修正
( P* S+ @* I/ p  c7 T    error=T(i)-ym  %计算误差
$ ^. n$ \1 N$ t6 |. X    E(i,j)=error;    %误差统计  
0 _$ R9 ~+ b( V: J    delta3=error
. \( J" \0 q) A& f$ H, e    delta21=error*W31
    delta22=error*W32
+ t! S: q+ y. \4 M7 K; r' T' f% [    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
, W4 y$ ~. a. u% v: A" }$ n    %修改权值
/ J2 H+ o; N+ h7 z3 B  v$ f4 U$ b    DW11=lr*delta1*k(i)
    DW21=lr*delta21*LB
1 Z' Z$ r/ q5 `" U' k) y) [8 p' g    DW22=lr*delta22*LB
    DW31=lr*delta3*LC1
, }) {( b7 I$ z, Z    DW32=lr*delta3*LC2
    Dtheta=lr*delta3;
* V" F# `& F0 N: c4 h% T3 `0 a    W11=W11+DW11+mc*DW110
    W21=W21+DW21+mc*DW210
. G+ }/ ?# A, Q- h; E1 e+ P# e0 X    W22=W22+DW22+mc*DW220
    W31=W31+DW31+mc*DW310
6 P( J4 r: j. y2 w1 I, E) V) e    W32=W32+DW32+mc*DW320
    theta=theta+Dtheta+mc*Dtheta0
    DW110=DW11
) W% h* V" M2 ?4 h7 C2 s5 W    DW210=DW21
    DW220=DW22
' A+ b: G& j' o- b( Q    DW310=DW31
    DW320=DW32
1 X$ q6 S- ^- p) s! S& k    Dtheta0=Dtheta
- X+ Y( W2 }9 e2 gend
' [2 a: i2 \9 b+ Y8 Aif sumsqr(E(:,j))<error_goal
4 ]( _( {  s! ^( a    disp('good')
    x=input('x=')
9 ^: ?4 Q* P9 |3 v& F    LB=1/(1+exp(-W11*x));   %LB层输出
$ m, u. s( y# K+ y, Q    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
6 X- q. F( w3 E3 n6 z   
    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
break;
end
0 M# a' p+ F, b' Y( n* |! @end  
/ ^3 b6 Q2 R$ o0 NYY=0;
    for i=1:length(k)
5 b: Y3 E) ?5 q/ _    LB=1/(1+exp(-W11*k(i)));   %LB层输出
7 {. D3 A# B: b( z    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
( {. ~' _- r0 O( q0 T9 e    LD=W31*LC1+W32*LC2;  %LD层输出
; b: J( a0 k; F) n& l6 ~. m    ym=LD-theta;   %网络输出值
8 c, _8 x8 w9 q! S9 |0 s& @    YY=[YY,ym];
   
    end
     YY=YY(2:end);
    YY

. }. L5 a$ M% k

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
8 O. O# A/ f& Y数学中国社区分享快乐！
7 @& R9 E2 \6 x% e% T: Q