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

yqxxqy

网络结构神经元数为输入：1，隐1层：1，隐二层：2，输出：1
* d! G) D% l4 k3 s3 y) Y求高手解答！！！
/ s2 w0 d: D) `我单步执行的结果是第一隐层输出饱和了。因为是为了后面的程序做准备，所以不想归一化。用工具箱的话是可以实现的。
但是从底层编的话该怎么改呢？望高手解答！心里那个急啊~~~~~~
* L1 H( r. I5 R) `clc
clear
: T: @* y) k! M( o' d. b- fk=[1,2,3,4,5];%输入
for i=1:length(k)
    T(i)=sum(k(1:i)); %输出
end
error_goal=0.01;%目标误差
max_epoch=200;%最大训练次数
lr=0.5;%学习速率
! {; r, C- n% ?. b%% 权值阀值初始化
+ v: i" Z% A3 v* [5 }E=zeros(length(k),max_epoch);%误差矩阵，便于查看误差
W11=rands(1);%权值初始化
4 Z: F5 _6 E2 |/ W4 ~, qW2=rands(2,1);
[W3,theta]=rands(1,2)
- S$ U& L+ X0 d$ VW21=W2(1)/4;W22=W2(2)/4;W31=W3(1)/4;W32=W3(2)/4;theta=theta/4;
5 c  ?: n1 Z! ^% @% |    DW110=0;
    DW210=0;
    DW220=0;
0 i3 l2 y3 [' {5 B8 Z    DW310=0;
    DW320=0;
    Dtheta0=0;
    mc=0.5;
%% 循环迭代
for j=1:max_epoch
) }) O0 Q7 P8 i4 I    lr=lr/sqrt(1+j)   %防止震荡，学习效率递减
for i=1:length(k)
9 I0 K0 t4 w4 c# d# `   
  t4 S$ u' `) C( C, [+ c# L   %% 网络输出计算
/ n1 W: {% ?2 m    LB=1/(1+exp(-W11*k(i)))   %LB层输出 sigmod函数
    LC1=LB*W21          %LC层输出 线性函数
    LC2=LB*W22    %LC层输出  线性函数
+ e" T* C/ m' E% x    LD=W31*LC1+W32*LC2  %LD层输出  线性函数
9 u8 S2 V& z4 n  \1 r$ C    ym=LD-theta   %网络输出值
        
   
    %% 权值修正
    error=T(i)-ym  %计算误差
, C* R% x, p, K; m    E(i,j)=error;    %误差统计  
: c% E4 S1 ^$ O. ?8 B0 }. I    delta3=error
    delta21=error*W31
    delta22=error*W32
/ r! b/ ]1 e& ]: [3 v6 r6 s    delta1=(delta21*W21+delta22*W22)*(LB*(1-LB))
    %修改权值
    DW11=lr*delta1*k(i)
    DW21=lr*delta21*LB
4 V4 B$ k" f0 ]  S8 G/ y- B    DW22=lr*delta22*LB
    DW31=lr*delta3*LC1
! {  |* n1 W) B5 O5 A    DW32=lr*delta3*LC2
6 g& P& a3 }* p" f4 l: E- n7 {; w* N* z( W    Dtheta=lr*delta3;
    W11=W11+DW11+mc*DW110
* V5 e3 N7 ~' f; N    W21=W21+DW21+mc*DW210
# y7 X% P9 ^+ ~  d5 H    W22=W22+DW22+mc*DW220
& K5 u( X3 H3 D* B4 j. H# q    W31=W31+DW31+mc*DW310
5 z. S; K4 b# S7 w$ U0 n! F    W32=W32+DW32+mc*DW320
2 V; W3 {' X) }9 E3 \0 c* q8 A9 D/ O    theta=theta+Dtheta+mc*Dtheta0
    DW110=DW11
1 f3 h! ~" S8 n: v% Z7 e6 J; w7 Q- T    DW210=DW21
    DW220=DW22
    DW310=DW31
, o5 L' q3 u/ ^    DW320=DW32
    Dtheta0=Dtheta
6 Z. f( Q; E8 ]0 a+ Q. V, }* x+ Mend
if sumsqr(E(:,j))<error_goal
    disp('good')
% k; H4 G# A( P- H% i9 ?    x=input('x=')
    LB=1/(1+exp(-W11*x));   %LB层输出
. M" k$ y$ h) C- n2 `' W7 F* B    LC1=LB*W21;           %LC层输出
; ~4 k* J  y7 Z    LC2=LB*W22;    %LC层输出
   
    LD=W31*LC1+W32*LC2  %LD层输出
    ym=LD-theta %网络输出值
break;
end
end  
  ~+ Y/ n" [* {9 K- H7 wYY=0;
    for i=1:length(k)
% O9 T$ f: U5 ~- x/ |, {/ k    LB=1/(1+exp(-W11*k(i)));   %LB层输出
    LC1=LB*W21;           %LC层输出
    LC2=LB*W22;    %LC层输出
    LD=W31*LC1+W32*LC2;  %LD层输出
    ym=LD-theta;   %网络输出值
    YY=[YY,ym];
  }9 h' S% t; Y, Q   
    end
. C2 v2 r$ X* v5 k     YY=YY(2:end);
    YY

. F. s. b; |% k2 e8 `8 o

杨帆

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