深度神经网络之反向传播算法

1.DNN反向传播算法简介

回顾我们前面学到的监督问题，通常会遇到这种情况，假如有m$m$个训练样本，分别为{(x1,y1),(x2,y2),(x3,y3),...,(xm,ym)}$\{(x_1,y_1),(x_2,y_2),(x_3,y_3),...,(x_m,y_m)\}$，其中x$x$为输入变量，特征维度为n_in，y为输出向量，特征维度为n_out。现在我们利用这m个训练样本来训练模型，当有测试样本(xtest,?)$(x_{test},?)$时，需要我们能够预测出ytest$y_{test}$向量的输出。

现在对应到我们的DNN模型之中，即输入层有n_in个神经元，输出层有n_out个神经元，再加上一些含有若干个神经元的隐含层。此时我们需要找到所有隐含层和输出层所对应的线性系数矩阵W、偏倚向量b，希望通过DNN对所有的训练样本计算后，计算结果能够等于或很接近样本输出，当有新的测试样本数据时，能够有效预测样本输出。但怎样找到合适的线形系数矩阵W和偏倚变量b呢?

回顾我们前面学习的机器学习之Logistic回归、机器学习之SVM支持向量机等机器学习算法，很容易联想到，我们可以用一个合适的损失函数来度量训练样本的输出损失。然后对损失函数优化，求损失函数最小化的极值，此时对应的线性系数矩阵W，偏倚变量b便是我们希望得到的结果。深度神经网络中，损失函数优化极值求解的过程，通常是利用梯度下降法迭代完成的。当然也可以利用其他的迭代方法，比如牛顿法或拟牛顿法。梯度下降算法以前在机器学习之线形回归中有过详细介绍，有兴趣可以回顾一下。

对DNN损失函数用梯度下降法进行迭代优化求极小值的过程，便是我们的反向传播算法(Back Propagation,BP)。

2.DNN反向传播算法数学推导

进行DNN反向传播算法之前，我们需要选择一个损失函数，来度量计算样本的输出和真实样本之间的损失。但训练时的计算样本输出怎么得到呢？

初始时，我们会随机选择一系列W,b，然后利用神经网络之前向传播算法中介绍到的al=σ(zl)=σ(Wlal−1+bl)$a^l=\sigma (z^l)=\sigma (W^l{a}^{l-1}+b^l)$，计算输出层所对应的aL$a^L$，此时的aL$a^L$便是DNN计算样本的输出。为专注DNN反向传播算法的推导，我们选择较为简单的损失函数，为此我们使用最常见的均方差来度量损失。

即对于每个样本，我们期望能够最小化下式，其中aL$a^L$和y$y$为特征维度的n_out的向量，||S||2$||S|{|}_2$为S的L2范数。

J(W,b,x,y)=12||aL−y||22$$J(W,b,x,y)=\frac{1}{2}||a^L-y|{|}_2^2$$

通过损失函数，我们能够用梯度下降法来迭代求解每一层的W，b。首先计算的是输出层，其中输出层的W，b满足下式

aL=σ(zL)=σ(WLaL−1+bL)$$a^L=\sigma (z^L)=\sigma (W^L{a}^{L-1}+b^L)$$

J(W,b,x,y)=12||aL−y||22=12||σ(WLaL−1+bL)−y||22$$J(W,b,x,y)=\frac{1}{2}||a^L-y|{|}_2^2=\frac{1}{2}||\sigma (W^L{a}^{L-1}+b^L)-y|{|}_2^2$$

然后对WL,bL$W^L,b^L$分别求偏导，其中符号⊙$\odot$表示Hadamard积，对于两个维度的向量A(a1,a2,a3,...,an)T$A(a_1,a_2,a_3,...,a_n{)}^T$和B(b1,b2,b3,...,bn)T$B(b_1,b_2,b_3,...,b_n{)}^T$，那么A⊙B=(a1b1,a2b2,a3b3,...,anbn)T$A\odot B=(a_1{b}_1,a_2{b}_2,a_3{b}_3,...,a_n{b}_n{)}^T$。之所以使用Hadamard积，是因为我们不了解**函数的形式，所以用Hadamard积来乘**函数的导数。另外补充矩阵求导的知识点，其中∂AB∂B=AT$\frac{\mathrm{\partial }AB}{\mathrm{\partial }B}=A^T$。

∂J(W,b,x,y)∂WL=∂J(W,b,x,y)∂zL∂zL∂WL=(aL−y)⊙σ′(zL)(aL−1)T$$\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{W}^L}=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^L}\frac{\mathrm{\partial }{z}^L}{\mathrm{\partial }{W}^L}=(a^L-y)\odot {\sigma }^{\prime }(z^L)(a^{L-1}{)}^T$$

∂J(W,b,x,y)∂bL=∂J(W,b,x,y)∂zL∂zL)∂bL=(aL−y)⊙σ′(zL)$$\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{b}^L}=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^L}\frac{\mathrm{\partial }{z}^L)}{\mathrm{\partial }{b}^L}=(a^L-y)\odot {\sigma }^{\prime }(z^L)$$

注意到在求解输出层W，b的时候，有公共部分∂J(W,b,x,y)∂zL$\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^L}$，因此我们可以把公共部分先算出来，记为

δL=∂J(W,b,x,y)∂zL=(aL−y)⊙σ′(zL)$${\delta }^L=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^L}=(a^L-y)\odot {\sigma }^{\prime }(z^L)$$

现在我们已经把输出层的梯度算出来了，那么如何求解L-1、L-2…层的梯度呢？这里我们需要进一步递推，对于第l$l$层的δl${\delta }^l$可以表示为

δl=∂J(W,b,x,y)∂zl=∂J(W,b,x,y)∂zL∂zL∂zL−1∂zL−1∂zL−2...∂zl+1∂zl$${\delta }^l=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^l}=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^L}\frac{\mathrm{\partial }{z}^L}{\mathrm{\partial }{z}^{L-1}}\frac{\mathrm{\partial }{z}^{L-1}}{\mathrm{\partial }{z}^{L-2}}...\frac{\mathrm{\partial }{z}^{l+1}}{\mathrm{\partial }{z}^l}$$

如果我们能够计算出第l$l$层的δl${\delta }^l$，那么对于该层的Wl,bl$W^l,b^l$也会很容易计算。为什么呢?注意到前向传播算法，我们有

zl=Wlal−1+bl$$z^l=W^l{a}^{l-1}+b^l$$

所以根据上式我们可以很方便的计算第l$l$层的Wl,bl$W^l,b^l$

∂J(W,b,x,y)∂Wl=∂J(W,b,x,y)∂zl∂zl∂Wl=δl(al−1)T$$\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{W}^l}=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^l}\frac{\mathrm{\partial }{z}^l}{\mathrm{\partial }{W}^l}={\delta }^l(a^{l-1}{)}^T$$

∂J(W,b,x,y)∂bl=∂J(W,b,x,y)∂zl∂zl)∂bl=δl$$\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{b}^l}=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^l}\frac{\mathrm{\partial }{z}^l)}{\mathrm{\partial }{b}^l}={\delta }^l$$

现在问题关键便是如何求解δl${\delta }^l$。假设我们已经得到第l+1$l+1$层的δl+1${\delta }^{l+1}$，那么如何得到第l$l$层的δl${\delta }^l$呢？我们注意到

δl=∂J(W,b,x,y)∂zl=∂J(W,b,x,y)∂zl+1∂zl+1∂zl=δl+1∂zl+1∂zl=$${\delta }^l=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^l}=\frac{\mathrm{\partial }J(W,b,x,y)}{\mathrm{\partial }{z}^{l+1}}\frac{\mathrm{\partial }{z}^{l+1}}{\mathrm{\partial }{z}^l}={\delta }^{l+1}\frac{\mathrm{\partial }{z}^{l+1}}{\mathrm{\partial }{z}^l}=$$

∂(δl+1)Tzl+1∂zl=∂(δl+1)T(Wl+1σ(zl)+bl+1)∂zl=∂(δl+1)TWl+1σ(zl)∂zl=$$\frac{\mathrm{\partial }({\delta }^{l+1}{)}^T{z}^{l+1}}{\mathrm{\partial }{z}^l}=\frac{\mathrm{\partial }({\delta }^{l+1}{)}^T(W^{l+1}\sigma (z^l)+b^{l+1})}{\mathrm{\partial }{z}^l}=\frac{\mathrm{\partial }({\delta }^{l+1}{)}^T{W}^{l+1}\sigma (z^l)}{\mathrm{\partial }{z}^l}=$$

((δl+1)TWl+1)T⊙σ′(zl)=(Wl+1)Tδl+1⊙σ′(zl)$$(({\delta }^{l+1}{)}^T{W}^{l+1}{)}^T\odot {\sigma }^{\prime }(z^l)=(W^{l+1}{)}^T{\delta }^{l+1}\odot {\sigma }^{\prime }(z^l)$$

现在我们已经得到δl${\delta }^l$的递推式，只要我们求出当前隐含层的δl${\delta }^l$，便能够得到Wl,bl$W^l,b^l$。

3.DNN反向传播算法过程

梯度下降算法有批量(Batch)，小批量(Mini-Batch)，随机三种方式，采用哪种方式取决于我们的问题而定。为简化描述，这里采用最基本的批量梯度下降法来描述反向传播算法。

输入：总层数L、各隐含层与输出层的神经元个数、**函数、损失函数、迭代步长α、最大迭代次数Max、停止迭代阈值ϵ、输入的m个训练样本(x1,y1),(x2,y2),...,(xm,ym)$(x_1,y_1),(x_2,y_2),...,(x_m,y_m)$。

输出：各隐含层与输出层的线性关系系数W和偏倚变量b。

• 初始化各隐藏层与输出层的线性关系系数矩阵W和偏倚向量b为随机值。
• for iter=1 to Max$foriter=1toMax$
  
  • for i =1 to m$fori=1tom$
    
    • 将a1$a^1$输入值设置为xi$x_i$
    • for l=2 to L$forl=2toL$，进行前向传播算法，计算ai,l=σ(zi,l)=σ(Wlai,l−1+bl)$a^{i,l}=\sigma (z^{i,l})=\sigma (W^l{a}^{i,l-1}+b^l)$
    • 通过损失函数计算输出层δi,L${\delta }^{i,L}$
    • for l=L to 2$forl=Lto2$，进行反向传播算法，计算δi,l=(Wl+1)Tδi,l+1⊙σ′(zi,l)${\delta }^{i,l}=(W^{l+1}{)}^T{\delta }^{i,l+1}\odot {\sigma }^{\prime }(z^{i,l})$
  • for l=2 to L$forl=2toL$，更新第l$l$层的Wl,bl$W^l,b^l$
    
    • Wl=Wl−α∑mi=1δi,l(ai,l−1)T$W^l=W^l-\alpha \sum _{i=1}^m{\delta }^{i,l}(a^{i,l-1}{)}^T$
    • bl=bl−α∑mi=1δi,l$b^l=b^l-\alpha \sum _{i=1}^m{\delta }^{i,l}$
  • 如果所有的W,b的变化值都小于停止迭代阈值ϵ，跳出循环。
• 输出各隐含层和输出层的线形关系系数矩阵W和偏倚向量b。

通过深度神经网络之中的前向传播算法和反向传播算法的结合，我们能够利用DNN模型去解决各种分类或回归问题，但对于不同问题，效果如何呢？是否会过拟合呢？我们将在下次文章中详细介绍损失函数和**函数的选择、正则化方面的知识点，来让深度神经网络能更精确的解决我们的问题。

参考

刘建平Pinard_深度神经网络(DNN)反向传播算法(BP)

4.推广

更多内容请关注公众号谓之小一，若有疑问可在公众号后台提问，随时回答，欢迎关注，内容转载请注明出处。