大家好,又见面了,我是你们的朋友全栈君。
BP算法是适用于多层神经网络的一种算法,它是建立在梯度下降法的基础上的。本文着重推导怎样利用梯度下降法来minimise Loss Function。
给出多层神经网络的示意图:
1.定义Loss Function
假设有一组数据样本 x 1 x^{1} x1, x 2 x^{2} x2,… ,每一个x都有很多个特征,输入x,会得到一个输出y,每一个输出都对应一个损失函数L,将所有L加起来就是total loss。
那么每一个L该如何定义呢?这里还是采用了交叉熵,如下所示:
这里的 y i y^{i} yi是真实输出, y ^ y\hat{} y^是y target,是人为定义的。最终Total Loss的表达式如下:
2.Gradient Descent
L对应了一个参数,即Network parameters θ(w1,w2…b1,b2…),那么Gradient Descent就是求出参数 θ ∗ \theta^{*} θ∗来minimise Loss Function,即:
梯度下降的具体步骤为:
3.求偏微分
从上图可以看出,这里难点主要是求偏微分,由于L是所有损失之和,因此我们只需要对其中一个损失求偏微分,最后再求和即可。
先抽取一个简单的神经元来解释:
先理一理各个变量之间的关系:我们要求的是Total Loss对参数w的偏导,而Total Loss是一个个小的l累加得到的,因此,我们只需要求得 ∂ l ∂ w \frac{\partial l}{\partial w} ∂w∂l,而 l = − y ^ l n y l=-\hat{y}lny l=−y^lny,其中 y ^ \hat{y} y^是人为定义的,跟w没有关系,因此我们只需要知道 ∂ y ∂ w \frac{\partial y}{\partial w} ∂w∂y。l跟z有关系,根据链式求导法则,我们需要求 ∂ l ∂ z \frac{\partial l}{\partial z} ∂z∂l和 ∂ z ∂ w \frac{\partial z}{\partial w} ∂w∂z,其中 ∂ z ∂ w \frac{\partial z}{\partial w} ∂w∂z的求解较为容易,如下图所示:
求 ∂ l ∂ z \frac{\partial l}{\partial z} ∂z∂l是一个难点,因为我们并不知道后面到底有多少层,也不知道情况到底有多复杂,我们不妨先取一种最简单的情况,如下所示:
4.反向传播
在第一张图里面,我们经过正向传播很容易求出了 ∂ z ∂ w \frac{\partial z}{\partial w} ∂w∂z,而对于 ∂ l ∂ z \frac{\partial l}{\partial z} ∂z∂l,则并不是那么好求。上图其实就是运用了反向传播的思想, 对于上图中 ∂ l ∂ z \frac{\partial l}{\partial z} ∂z∂l最后的表达式,我们可以换一种结构,如下所示:
l对两个z的偏导我们假设是已知的,并且在这里是作为输入,三角形结构可以理解为一个乘法运算电路,其放大系数为 σ ′ ( z ) \sigma {}'(z) σ′(z)。但是在实际情况中,l对两个z的偏导是未知的。假设神经网络最终的结构就是如上图所示,那么我们的问题已经解决了:
其中:
但是假如该神经元不是最后一层,我们又该如何呢?比如又多了一层,如下所示:
那么我们只要知道 ∂ l ∂ z a \frac{\partial l}{\partial z_{a}} ∂za∂l和 ∂ l ∂ z b \frac{\partial l}{\partial z_{b}} ∂zb∂l,我们同样可以算出 ∂ l ∂ z ′ \frac{\partial l}{\partial z{}’} ∂z′∂l以及 ∂ l ∂ z ′ ′ \frac{\partial l}{\partial z{}”} ∂z′′∂l,原理跟上面类似,如下所示:
∂ l ∂ z b \frac{\partial l}{\partial z_{b}} ∂zb∂l同样是l先对y求导,y再对 z b z_{b} zb求导。
那假设我们再加一层呢?再加两层呢?再加三层呢?。。。,情况还是一样的,还是先求l对最后一层z的导数,乘以权重相加后最后再乘上 σ ′ ( z ′ ′ , z ′ ′ ′ , . . . ) \sigma {}'(z{}”,z{}”’,…) σ′(z′′,z′′′,...)即可。
最后给一个实例:
它的反向传播图长这样:
我们可以很轻松的算出 ∂ l ∂ z 5 \frac{\partial l}{\partial z_{5}} ∂z5∂l和 ∂ l ∂ z 6 \frac{\partial l}{\partial z_{6}} ∂z6∂l,算出这两个之后,根据上面我们找到的关系式,我们也可以轻易算出 ∂ l ∂ z 3 \frac{\partial l}{\partial z_{3}} ∂z3∂l和 ∂ l ∂ z 4 \frac{\partial l}{\partial z_{4}} ∂z4∂l,最后再算出 ∂ l ∂ z 1 \frac{\partial l}{\partial z_{1}} ∂z1∂l和 ∂ l ∂ z 2 \frac{\partial l}{\partial z_{2}} ∂z2∂l。然后 ∂ l ∂ z 1 \frac{\partial l}{\partial z_{1}} ∂z1∂l和 ∂ l ∂ z 2 \frac{\partial l}{\partial z_{2}} ∂z2∂l再分别乘上x1和x2,就是我们最终要找的 ∂ l ∂ w 1 \frac{\partial l}{\partial w_{1}} ∂w1∂l和 ∂ l ∂ w 2 \frac{\partial l}{\partial w_{2}} ∂w2∂l。
我们不难发现,这种计算方式很清楚明了地体现了“反向传播”四个字。
好了,目标达成!!
5.总结
通过Forward Pass我们求得 ∂ z ∂ w = a \frac{\partial z}{\partial w}=a ∂w∂z=a,然后通过Backward Pass我们求得 ∂ l ∂ z \frac{\partial l}{\partial z} ∂z∂l,二者相乘,就是 ∂ l ∂ w \frac{\partial l}{\partial w} ∂w∂l。利用上述方法求得所有参数的值之后,我们就可以用梯度下降法来更新参数,直至找到最优解。
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/132274.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...