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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
3 i: f$ l' X+ g) D求高手解答！！！
我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
+ s, L. p. u; A但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
clc
clear
k=[1,2,3,4,5];%输入
! o7 [* A: X# Yfor i=1:length(k)
6 ?* m4 N+ C$ r7 Y3 Z  I0 e    T(i)=sum(k(1:i)); %输出
" R) S; f. n, F2 N- Jend
( F* O5 I! o2 s- \3 c- Cerror_goal=0.01;%目标误差
1 O9 U0 ?8 z. r- F, ^, e2 X2 }max_epoch=200;%最大训练次数
lr=0.5;%学习速率
9 n  @& b% p9 _1 n%% 权值阀值初始化
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;
    DW110=0;
0 d* ]( P/ p! ]    DW210=0;
    DW220=0;
    DW310=0;
& d; P; ~0 W; V7 c! n) n8 v5 g    DW320=0;
; |7 W& b6 j. Z& K    Dtheta0=0;
    mc=0.5;
7 f5 D! E* {. z) g%% 循环迭代
for j=1:max_epoch
    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
for i=1:length(k)
- n' s# W1 G' [3 U, f+ ?   
7 [- R+ `, S: X( u) Q0 _   %% 网络输出计算
& X7 _; L! Z! d! h& `    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
( |/ g2 \4 g& _  b    LC2=LB*W22    %LC层输出  线性函数
    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
: m$ ]' y+ ]/ g( Q0 H' Q    ym=LD-theta   %网络输出值
        
; N2 N! }: F4 W" D2 c   
    %% 权值修正
6 `( b& P/ d( _0 y$ I' d    error=T(i)-ym  %计算误差
+ o3 N1 e2 |$ S: ~    E(i,j)=error;    %误差统计  
    delta3=error
- t3 e, s3 D  l    delta21=error*W31
    delta22=error*W32
; D( b. n* o0 \1 C% Z# F3 }% q$ }    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
4 y: }; b( S- a3 p5 f    %修改权值
    DW11=lr*delta1*k(i)
    DW21=lr*delta21*LB
% j6 b; h1 ?+ ]% R  s    DW22=lr*delta22*LB
% S( e' E2 j' H) B    DW31=lr*delta3*LC1
    DW32=lr*delta3*LC2
/ C9 b; [/ i2 h* G* j- D    Dtheta=lr*delta3;
5 H6 J/ x! J9 @, Z    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
    DW210=DW21
    DW220=DW22
2 l! }7 B% P/ v0 B# ?) |' Z    DW310=DW31
" W. m. q$ |- I    DW320=DW32
    Dtheta0=Dtheta
end
* k8 B" j; }" N, \. [! ^+ n6 x0 o2 eif sumsqr(E(:,j))<error_goal
+ P4 k/ \; X# X    disp('good')
# N) ~$ x' }+ [    x=input('x=')
    LB=1/(1+exp(-W11*x));   %LB层输出
6 V/ [4 v  v  s* A/ D9 {0 a    LC1=LB*W21;           %LC层输出
5 Z$ {  v. D3 M1 L    LC2=LB*W22;    %LC层输出
! P  s. w/ B7 l( f% `   
    LD=W31*LC1+W32*LC2  %LD层输出
% Q" ?% e4 `9 r$ o9 r- R    ym=LD-theta %网络输出值
break;
end
end  
YY=0;
    for i=1:length(k)
    LB=1/(1+exp(-W11*k(i)));   %LB层输出
    LC1=LB*W21;           %LC层输出
0 f  w, u0 M  }; K    LC2=LB*W22;    %LC层输出
! m7 Y+ H5 _( ~; K+ n' S$ K5 D7 g    LD=W31*LC1+W32*LC2;  %LD层输出
, F2 X, u9 e9 |    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
   
6 d- \, B4 [: ~6 q! J0 V6 Z! a+ C    end
8 {, Y/ G- ^9 ~2 u3 t9 S& Z# K8 ~     YY=YY(2:end);
    YY

杨帆

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