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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
7 R3 r  \6 f3 K: j) p求高手解答！！！
3 x" Z7 W2 R4 z; z* k( B7 g我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
! [6 U! B7 r- m3 E5 o9 x6 Nclc
+ p) S0 f& u$ d8 O) \- {6 j7 |clear
k=[1,2,3,4,5];%输入
for i=1:length(k)
# S% r8 S$ f! }" M% v    T(i)=sum(k(1:i)); %输出
end
/ w5 t7 F- Z# U) R# `error_goal=0.01;%目标误差
max_epoch=200;%最大训练次数
lr=0.5;%学习速率
% N$ u* |, e1 Q) ]* u%% 权值阀值初始化
E=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
7 Q; k9 Y; ~4 ^' _& y1 \+ FW11=rands(1);%权值初始化
. ~. _. t/ ^* s5 B0 d& o. E/ HW2=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;
; F$ k2 T6 ?5 D    DW110=0;
    DW210=0;
    DW220=0;
, F2 o7 B+ G6 ]) y+ s7 K0 H0 g& w    DW310=0;
    DW320=0;
    Dtheta0=0;
1 ?. e: L! b! M4 V* G    mc=0.5;
( j4 n/ K: k, _* G%% 循环迭代
for j=1:max_epoch
    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
6 j# p3 G+ r$ q% ^" B! @for i=1:length(k)
- S% ]; H4 s& _/ o2 [( J7 i. z# d   
   %% 网络输出计算
    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
    LC2=LB*W22    %LC层输出  线性函数
+ o7 r/ p2 N2 e& B# `0 J. @    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
; E9 M; I1 Y- G; n! u        
   
% T+ \/ j* l4 w    %% 权值修正
' I$ H4 J$ F% U" g$ A# H    error=T(i)-ym  %计算误差
' [# E, f" x: Q. T7 {. P    E(i,j)=error;    %误差统计  
# d# G6 ~% {7 |# w2 I: [6 w6 `    delta3=error
    delta21=error*W31
    delta22=error*W32
" h& l4 u  s$ ]2 X" D' |    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
    %修改权值
$ O2 B) T$ v: s* s) S3 m    DW11=lr*delta1*k(i)
* Z1 _- U# x. K+ X' I    DW21=lr*delta21*LB
    DW22=lr*delta22*LB
3 o. Q7 M/ i) V5 ]0 L) Y" h' g    DW31=lr*delta3*LC1
, i7 A' D; D' g+ `- @) P2 Q" E' h    DW32=lr*delta3*LC2
9 N) s- H  v6 u# i, _    Dtheta=lr*delta3;
    W11=W11+DW11+mc*DW110
( h% e! F! o( D& Z    W21=W21+DW21+mc*DW210
    W22=W22+DW22+mc*DW220
/ }9 }+ E" s; \    W31=W31+DW31+mc*DW310
    W32=W32+DW32+mc*DW320
# T3 o0 l5 v# d% r# T8 D! o    theta=theta+Dtheta+mc*Dtheta0
7 t$ k* I) x' ?% W: I( W    DW110=DW11
    DW210=DW21
! t# ]/ T/ z5 ]3 E2 s    DW220=DW22
    DW310=DW31
" ~  f0 H0 U4 f    DW320=DW32
    Dtheta0=Dtheta
- E3 Q1 }: k. \end
! N7 z: r( ]" V, Z; ^if sumsqr(E(:,j))<error_goal
    disp('good')
    x=input('x=')
" I$ J, J3 Z# H: V' i    LB=1/(1+exp(-W11*x));   %LB层输出
    LC1=LB*W21;           %LC层输出
, z$ U5 h  v# e7 }5 I; B( E* i    LC2=LB*W22;    %LC层输出
   
    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
break;
( ?; F1 B" Q1 g0 \+ E( l% {. _end
end  
YY=0;
% r+ V' d- G; {' |. R! Z    for i=1:length(k)
    LB=1/(1+exp(-W11*k(i)));   %LB层输出
4 ^. X6 U' C9 r, t0 O, v  H    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
    LD=W31*LC1+W32*LC2;  %LD层输出
! M8 q% t, L( o9 W" s2 ~" l: Z    ym=LD-theta;   %网络输出值
: @* e" A* t! E4 K8 L! m7 }. `2 P7 z    YY=[YY,ym];
   
    end
     YY=YY(2:end);
    YY
0 _* j4 I: P( C$ y0 ]1 M1 g

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
& x9 U3 r5 P  ~- k' Z7 L0 q数学中国社区分享快乐！