极市导读
本文介绍了Batch Normalization的过程及其实现,它能帮助解决深度神经网络中输入不稳定的问题,并使得神经网络的梯度大小相对固定。
统计学习中的一个很重要的假设就是输入的分布是相对稳定的。如果这个假设不满足,则模型的收敛会很慢,甚至无法收敛。所以,对于一般的统计学习问题,在训练前将数据进行归一化或者白化(whitening)是一个很常用的trick。
但这个问题在深度神经网络中变得更加难以解决。在神经网络中,网络是分层的,可以把每一层视为一个单独的分类器,将一个网络看成分类器的串联。这就意味着,在训练过程中,随着某一层分类器的参数的改变,其输出的分布也会改变,这就导致下一层的输入的分布不稳定。分类器需要不断适应新的分布,这就使得模型难以收敛。
对数据的预处理可以解决第一层的输入分布问题,而对于隐藏层的问题无能为力,这个问题就是Internal Covariate Shift。而Batch Normalization其实主要就是在解决这个问题。
除此之外,一般的神经网络的梯度大小往往会与参数的大小相关(仿射变换),且随着训练的过程,会产生较大的波动,这就导致学习率不宜设置的太大。Batch Normalization使得梯度大小相对固定,一定程度上允许我们使用更高的学习率。
Batch Normalization的过程很简单。我们假定我们的输入是一个大小为 的mini-batch ,通过下面的四个式子计算得到的 就是Batch Normalization(BN)的值。
首先,由(2.1)和(2.2)得到mini-batch的均值和方差,之后进行(2.3)的归一化操作,在分母加上一个小的常数是为了避免出现除0操作。最后的(2.4)再对
BN层通常添加在隐藏层的激活函数之前,线性变换之后。如果我们把(2.4)和之后的激活函数放在一起看,可以将他们视为一层完整的神经网络(线性+激活)。(注意BN的线性变换和一般隐藏层的线性变换仍有区别,前者是element-wise的,后者是矩阵乘法。)
此时,
另外,
在训练过程中,还需要维护一个移动平均的均值和方差,这两个移动平均会用于推断过程。
再来用反向传播求梯度。我们先画出BN的计算图,由于节点的文本不支持Tex,其中x_即为
Batch Normalization的计算图
假定损失函数为L,已知L相对于
前两个比较直观,求
从计算图和(2.3)中可以看出,要求
我们还需要求
(3.6)可以由式(2.1)直接得到,(3.7)是运用链式法则的结果。
算
(3.8)可以由式(2.2)直接得到。
最后我们还需要方差的导数。
至此,(3.4)所需要的偏导都已求出,组合起来,我们得到最终的式子。
这个式子仍然具有进一步化简的空间。在展开之后,由式(2.3),通过拼凑,可以将
参考cs231n的assignment2,我们将Batch Normalization分成正向(只包括训练)和反向两个过程。
正向过程的参数x是一个mini-batch的数据,gamma和beta是BN层的参数,bn_param是一个字典,包括
反向过程的参数是来自上一层的误差信号dout,以及正向过程中存储的中间变量cache,最后返回
实现与推导的不同在于,实现是对整个batch的操作。
import numpy as np
def batchnorm_forward(x, gamma, beta, bn_param):
# read some useful parameter
N, D = x.shape
eps = bn_param.get('eps', 1e-5)
momentum = bn_param.get('momentum', 0.9)
running_mean = bn_param.get('running_mean', np.zeros(D, dtype=x.dtype))
running_var = bn_param.get('running_var', np.zeros(D, dtype=x.dtype))
# BN forward pass
sample_mean = x.mean(axis=0)
sample_var = x.var(axis=0)
x_ = (x - sample_mean) / np.sqrt(sample_var + eps)
out = gamma * x_ + beta
# update moving average
running_mean = momentum * running_mean + (1-momentum) * sample_mean
running_var = momentum * running_var + (1-momentum) * sample_var
bn_param['running_mean'] = running_mean
bn_param['running_var'] = running_var
# storage variables for backward pass
cache = (x_, gamma, x - sample_mean, sample_var + eps)
return out, cache
def batchnorm_backward(dout, cache):
# extract variables
N, D = dout.shape
x_, gamma, x_minus_mean, var_plus_eps = cache
# calculate gradients
dgamma = np.sum(x_ * dout, axis=0)
dbeta = np.sum(dout, axis=0)
dx_ = np.matmul(np.ones((N,1)), gamma.reshape((1, -1))) * dout
dx = N * dx_ - np.sum(dx_, axis=0) - x_ * np.sum(dx_ * x_, axis=0)
dx *= (1.0/N) / np.sqrt(var_plus_eps)
return dx, dgamma, dbeta
推荐阅读