BP神经网络原理

BP神经网络原理2.1基本BP算法公式推导基本BP算法包括两个方面：信号的前向传播和误差的反向传播。即计算实际输出时按从输入到输出的方向进行，而权值和阈值的修正从输出到输入的方向进行。图2-1BP网络结构Fig.2-1StructureofBPnetwork图中：jx表示输入层第j个节点的输入，j=1,…,M；ijw表示隐含层第i个节点到输入层第j个节点之间的权值；i表示隐含层第i个节点的阈值；()x表示隐含层的激励函数；kiw表示输出层第k个节点到隐含层第i个节点之间的权值，i=1,…,q；ka表示输出层第k个节点的阈值,k=1,…,L；x表示输出层的激励函数；ko表示输出层第k个节点的输出。（1）信号的前向传播过程隐含层第i个节点的输入neti：1a1kiwijwLakaqi………………1xjxMx1okoLo输出变量输入变量输入层隐含层输出层1Miijjijnetwx（3-1）隐含层第i个节点的输出yi：1()()Miiijjijynetwx（3-2）输出层第k个节点的输入netk：111()qqMkkiikkiijjikiijnetwyawwxa（3-3）输出层第k个节点的输出ok：111()()()qqMkkkiikkiijjikiijonetwyawwxa（3-4）（2）误差的反向传播过程误差的反向传播，即首先由输出层开始逐层计算各层神经元的输出误差，然后根据误差梯度下降法来调节各层的权值和阈值，使修改后的网络的最终输出能接近期望值。对于每一个样本p的二次型误差准则函数为Ep：211()2LpkkkETo（3-5）系统对P个训练样本的总误差准则函数为：2111()2PLppkkpkETo（3-6）根据误差梯度下降法依次修正输出层权值的修正量Δwki，输出层阈值的修正量Δak，隐含层权值的修正量Δwij，隐含层阈值的修正量i。kikiwEw；kkEaa；ijijEww；iiE（3-7）输出层权值调整公式：kikkkkkikkkikiwnetnetooEwnetnetEwEw（3-8）输出层阈值调整公式：kkkkkkkkkknetonetEEEaanetaoneta（3-9）隐含层权值调整公式：iiiijijiijiiijnetynetEEEwwnetwynetw（3-10）隐含层阈值调整公式：iiiiiiiiiinetynetEEEnetynet（3-11）又因为：11()PLppkkpkkEToo（3-12）ikikywnet，1kkneta，ijijnetxw，1iinet（3-13）11()'()PLppkkkkipkiETonetwy（3-14）)(iiinetnety（3-15）'()kkkonetnet（3-16）所以最后得到以下公式：11()'PLppkikkkipkwTonety（3-17）11()'PLppkkkkpkaTonet（3-18）11()'()PLppijkkkkiijpkwTonetwnetx（3-19）11()'()PLppikkkkiipkTonetwnet（3-20）图2-2BP算法程序流程图Fig.2-2TheflowchartoftheBPalgorithmprogram2.2基本BP算法的缺陷BP算法因其简单、易行、计算量小、并行性强等优点，目前是神经网络训练采用最多也是最成熟的训练算法之一。其算法的实质是求解误差函数的最小值问题，由于它采用非线性规划中的最速下降方法，按误差函数的负梯度方向修改权值，因而通常存在以下问题：（1）学习效率低，收敛速度慢（2）易陷入局部极小状态结束参数初始化：最大训练次数，学习精度，隐节点数，初始权值、阈值，初始学习速率等批量输入学习样本并且对输入和输出量进行归一化处理计算各层的输入和输出值计算输出层误差E(q)E(q)ε修正权值和阈值YesNYes开始2.3BP算法的改进2.3.1附加动量法附加动量法使网络在修正其权值时，不仅考虑误差在梯度上的作用，而且考虑在误差曲面上变化趋势的影响。在没有附加动量的作用下，网络可能陷入浅的局部极小值，利用附加动量的作用有可能滑过这些极小值。该方法是在反向传播法的基础上在每一个权值（或阈值）的变化上加上一项正比于前次权值（或阈值）变化量的值，并根据反向传播法来产生新的权值（或阈值）变化。带有附加动量因子的权值和阈值调节公式为：)()1()1(kwmcpmckwijjiij)()1()1(kbmcmckbiii其中k为训练次数，mc为动量因子，一般取0.95左右。附加动量法的实质是将最后一次权值（或阈值）变化的影响，通过一个动量因子来传递。当动量因子取值为零时，权值（或阈值）的变化仅是根据梯度下降法产生；当动量因子取值为1时，新的权值（或阈值）变化则是设置为最后一次权值（或阈值）的变化，而依梯度法产生的变化部分则被忽略掉了。以此方式，当增加了动量项后，促使权值的调节向着误差曲面底部的平均方向变化，当网络权值进入误差曲面底部的平坦区时，i将变得很小，于是)()1(kwkwijij,从而防止了0ijw的出现，有助于使网络从误差曲面的局部极小值中跳出。根据附加动量法的设计原则，当修正的权值在误差中导致太大的增长结果时，新的权值应被取消而不被采用，并使动量作用停止下来，以使网络不进入较大误差曲面；当新的误差变化率对其旧值超过一个事先设定的最大误差变化率时，也得取消所计算的权值变化。其最大误差变化率可以是任何大于或等于1的值。典型的取值取1.04。所以，在进行附加动量法的训练程序设计时，必须加进条件判断以正确使用其权值修正公式。训练程序设计中采用动量法的判断条件为：mcmc95.00其它)1()(04.1*)1()(kEkEkEkE，E（k）为第k步误差平方和。2.3.2自适应学习速率对于一个特定的问题，要选择适当的学习速率不是一件容易的事情。通常是凭经验或实验获取，但即使这样，对训练开始初期功效较好的学习速率，不见得对后来的训练合适。为了解决这个问题，人们自然想到在训练过程中，自动调节学习速率。通常调节学习速率的准则是：检查权值是否真正降低了误差函数，如果确实如此，则说明所选学习速率小了，可以适当增加一个量；若不是这样，而产生了过调，那幺就应该减少学习速率的值。下式给出了一个自适应学习速率的调整公式：)()(7.0)(05.1)1(kkkk其它)(04.1)1()()1(kEkEkEkE，E（k）为第k步误差平方和。初始学习速率(0)的选取范围可以有很大的随意性。2.3.3动量-自适应学习速率调整算法当采用前述的动量法时，BP算法可以找到全局最优解，而当采用自适应学习速率时，BP算法可以缩短训练时间，采用这两种方法也可以用来训练神经网络，该方法称为动量-自适应学习速率调整算法。2.4网络的设计2.4.1网络的层数理论上已证明：具有偏差和至少一个S型隐含层加上一个线性输出层的网络，能够逼近任何有理数。增加层数可以更进一步的降低误差，提高精度，但同时也使网络复杂化，从而增加了网络权值的训练时间。而误差精度的提高实际上也可以通过增加神经元数目来获得，其训练效果也比增加层数更容易观察和调整。所以一般情况下，应优先考虑增加隐含层中的神经元数。2.4.2隐含层的神经元数网络训练精度的提高，可以通过采用一个隐含层，而增加神经元数了的方法来获得。这在结构实现上，要比增加隐含层数要简单得多。那么究竟选取多少隐含层节点才合适？这在理论上并没有一个明确的规定。在具体设计时，比较实际的做法是通过对不同神经元数进行训练对比，然后适当地加上一点余量。2.4.3初始权值的选取由于系统是非线性的，初始值对于学习是否达到局部最小、是否能够收敛及训练时间的长短关系很大。如果初始值太大，使得加权后的输入和n落在了S型激活函数的饱和区，从而导致其导数f’(n)非常小，而在计算权值修正公式中，因为)('nf，当f’(n)0时，则有0。这使得0ijw，从而使得调节过程几乎停顿下来。所以一般总是希望经过初始加权后的每个神经元的输出值都接近于零，这样可以保证每个神经元的权值都能够在它们的S型激活函数变化最大之处进行调节。所以，一般取初始权值在（-1，1）之间的随机数。2.4.4学习速率学习速率决定每一次循环训练中所产生的权值变化量。大的学习速率可能导致系统的不稳定；但小的学习速率导致较长的训练时间，可能收敛很慢，不过能保证网络的误差值不跳出误差表面的低谷而最终趋于最小误差值。所以在一般情况下，倾向于选取较小的学习速率以保证系统的稳定性。学习速率的选取范围在0.01-0.8之间。3BP神经网络的应用现给出一药品商店一年当中12个月的药品销售量（单位：箱）如下：205623952600229816341600187314871900150020461556训练一个BP网络，用当前的所有数据预测下一个月的药品销售量。有两种方法实现，一种是编写matlab程序，一种是使用nntool工具箱。3.1matlab程序实现我们用前三个月的销售量预测下一个月的销售量，也就是用1-3月的销售量预测第4个月的销售量，用2-4个月的销售量预测第5个月的销售量，如此循环下去，直到用9-11月预测12月份的销售量。这样训练BP神经网络后，就可以用10-12月的数据预测来年一月的销售量。实现程序如下：p=[205623952600;239526002298;260022981634;229816341600;163416001873;160018731478;187314781900;147819001500;190015002046;]t=[229816341600187314871900150020461556];pmax=max(p);pmax1=max(pmax);pmin=min(p);pmin1=min(pmin);fori=1:9%归一化处理p1(i,:)=(p(i,:)-pmin1)/(pmax1-pmin1);endt1=(t-pmin1)/(pmax1-pmin1);t1=t1';net=newff([01;01;01],[71],{'tansig','logsig'},'traingd');fori=1:9net.trainParam.epochs=15000;net.trainParam.goal=0.01;LP.lr=0.1;net=train(net,p1(i,:)',t1(i));endy=sim(net,[150020461556]');y1=y*(pmax1-pmin1)+pmin1;如果神经网络的训练函数使用trainlm，则仿真步骤会很少，但需要较大的系统内存。经预测，来年一月的销售量（y1）为1.4848e+003箱（每次运行后的结果可能不同）。3.2nntool神经网络工具箱的使用1)在matlab（7.4）命令窗口键入nntool命令打开神经网络工具箱。如图：2)点击Import按钮两次，分别把输入向量和目标输出加入到对应的窗口（[Inputs]和[Targets]）中，有两种可供选择的加入对象（点击Import后可以看见），一种是把当前工作区中的某个矩阵加入，另一种是通过.mat文件读入。3)点击[NewNetwork]按钮，填入各参