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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
求高手解答！！！
0 U; u& \5 h. n1 I8 R; ]  I. y我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
" e) [9 J- m+ q. j1 r但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
clc
clear
k=[1,2,3,4,5];%输入
4 m  |/ B! |0 vfor i=1:length(k)
    T(i)=sum(k(1:i)); %输出
end
error_goal=0.01;%目标误差
max_epoch=200;%最大训练次数
5 e5 L8 `( g9 c' A8 l- Elr=0.5;%学习速率
%% 权值阀值初始化
* K; J/ v8 ~* e6 ME=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
W11=rands(1);%权值初始化
W2=rands(2,1);
: A/ W4 a: d2 O# ~% o3 I[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;
5 b- y; A/ v/ u9 E7 B3 W    DW210=0;
    DW220=0;
    DW310=0;
# I! {/ I1 r0 r  n) @; p( |    DW320=0;
6 X8 \5 h+ ?9 @% b6 e9 s' @    Dtheta0=0;
    mc=0.5;
1 s; E- X9 ?9 d, H%% 循环迭代
3 v8 p' V% }: Y& T for j=1:max_epoch
! b; D' C! M* |' f% ]8 c! {    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
for i=1:length(k)
! A' Z! q0 `+ ?  {2 A" A  p7 E   
   %% 网络输出计算
    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
8 Y, |) K1 i. Q' `. \6 \    LC1=LB*W21          %LC层输出 线性函数
    LC2=LB*W22    %LC层输出  线性函数
- d- A% |' Z  E/ x    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
    ym=LD-theta   %网络输出值
        
   
    %% 权值修正
! S. h) y5 h! E" L# F; ]    error=T(i)-ym  %计算误差
+ a' l" w: \+ I/ ^' D! T5 T* [    E(i,j)=error;    %误差统计  
    delta3=error
8 Z- a% c- [- S; A    delta21=error*W31
2 u0 v& [' I+ I    delta22=error*W32
    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
: k$ O3 @0 P  Z0 }! V( y2 U    %修改权值
1 N; o6 C8 Y$ M; E: u1 d* z. x" z    DW11=lr*delta1*k(i)
    DW21=lr*delta21*LB
    DW22=lr*delta22*LB
    DW31=lr*delta3*LC1
* N" R/ q* s7 z    DW32=lr*delta3*LC2
    Dtheta=lr*delta3;
    W11=W11+DW11+mc*DW110
    W21=W21+DW21+mc*DW210
4 V' K4 w4 A& o3 b, s    W22=W22+DW22+mc*DW220
    W31=W31+DW31+mc*DW310
: ]% h. M* k# C5 h' C    W32=W32+DW32+mc*DW320
' B0 F1 l9 }: m* x: h    theta=theta+Dtheta+mc*Dtheta0
    DW110=DW11
    DW210=DW21
% b! k& {6 G! Z    DW220=DW22
    DW310=DW31
, `4 `& Y2 U. Z; ^* j- q9 Z: t    DW320=DW32
    Dtheta0=Dtheta
end
2 B6 r0 S8 c: h/ X- iif sumsqr(E(:,j))<error_goal
6 J  y& w# f; }9 {; r& v    disp('good')
( P, n, j5 K& e0 X4 v9 ^    x=input('x=')
5 P. I4 p% t( S. |, A    LB=1/(1+exp(-W11*x));   %LB层输出
    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
   
1 S7 p. E/ D' N# r    LD=W31*LC1+W32*LC2  %LD层输出
$ g$ }" l# G3 t1 \, x    ym=LD-theta %网络输出值
break;
$ o9 b& G. X# j4 Y0 qend
end  
YY=0;
    for i=1:length(k)
$ ~$ u7 l5 H) P: d. b, I+ \" ~1 r    LB=1/(1+exp(-W11*k(i)));   %LB层输出
    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
    LD=W31*LC1+W32*LC2;  %LD层输出
0 E6 q/ r( {8 @0 \# y    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
   
    end
     YY=YY(2:end);
    YY
/ r# G. O9 W& Y) Q# ^4 M2 D  j

杨帆

既然你诚信诚意的推荐了，那我就勉为其难的听听吧！
7 N1 K8 r7 K3 B数学中国社区分享快乐！
- u: _; R" c# r