Layer Normalization

tags: Normalization
前言
在上一篇的文章中我们介绍了BN的计算方法并且讲解了BN如何应用在MLP以及CNN中如何使用BN。在文章的最后，我们指出BN并不适用于RNN等动态网络和batchsize较小的时候效果不好。Layer Normalization（LN）的提出有效的解决BN的这两个问题。LN和BN不同点是归一化的维度是互相垂直的，如图1所示。在图1中NNN
表示样本轴，CCC
表示通道轴，FFF
是每个通道的特征数量。BN如右侧所示，它是取不同样本的同一个通道的特征做归一化；LN则是如左侧所示，它取的是同一个样本的不同通道做归一化。
 
图1：LN(左)和BN(右)对比示意图
1. BN的问题
1.1 BN与Batch Size
如图1右侧部分，BN是按照样本数计算归一化统计量的，当样本数很少时，比如说只有4个。这四个样本的均值和方差便不能反映全局的统计分布息，所以基于少量样本的BN的效果会变得很差。在一些场景中，比如说硬件资源受限，在线学习等场景，BN是非常不适用的。
1.2 BN与RNN
RNN可以展开成一个隐藏层共享参数的MLP，随着时间片的增多，展开后的MLP的层数也在增多，最终层数由输入数据的时间片的数量决定，所以RNN是一个动态的网络。
在一个batch中，通常各个样本的长度都是不同的，当统计到比较靠后的时间片时，例如图2中t>4t>4t>4
时，这时只有一个样本还有数据，基于这个样本的统计信息不能反映全局分布，所以这时BN的效果并不好。
另外如果在测试时我们遇到了长度大于任何一个训练样本的测试样本，我们无法找到保存的归一化统计量，所以BN无法运行。
 
图2：RNN中使用BN会导致batchsize过小的问题
2. LN详解
2.1 MLP中的LN
通过第一节的分析，我们知道BN的两个缺点的产生原因均是因为计算归一化统计量时计算的样本数太少。LN是一个独立于batch size的算法，所以无论样本数多少都不会影响参与LN计算的数据量，从而解决BN的两个问题。LN的做法如图1左侧所示：根据样本的特征数做归一化。
先看MLP中的LN。设HHH
是一层中隐层节点的数量，lll
是MLP的层数，我们可以计算LN的归一化统计量μ\muμ
和σ\sigmaσ
：
μl=1H∑i=1Hailσl=1H∑i=1H(ail−μl)2\mu^l = \frac{1}{H}\sum_{i=1}^H a_i^l
\qquad
\sigma^l = \sqrt{\frac{1}{H} \sum_{i=1}^H(a_i^l
 - \mu^l)^2}μl=H1​i=1∑H​ail​σl=H1​i=1∑H​(ail​−μl)2​
注意上面统计量的计算是和样本数量没有关系的，它的数量只取决于隐层节点的数量，所以只要隐层节点的数量足够多，我们就能保证LN的归一化统计量足够具有代表性。通过μl\mu^lμl
和σl\sigma^lσl
可以得到归一化后的值a^l\hat{a}^la^l
：
a^l=al−μl(σl)2+ϵ\hat{\mathbf{a}}^l = \frac{\mathbf{a}^l - \mu^l}{(\sigma^l)^2 + \sqrt{\epsilon}}a^l=(σl)2+ϵ​al−μl​
其中ϵ\epsilonϵ
是一个很小的小数，防止除0（论文中忽略了这个参数）。
在LN中我们也需要一组参数来保证归一化操作不会破坏之前的信息，在LN中这组参数叫做增益（gain）ggg
和偏置（bias）bbb
（等同于BN中的γ\gammaγ
和β\betaβ
）。假设激活函数为fff
，最终LN的输出为：
hl=f(gl⊙a^l+bl)\mathbf{h}^l = f(\mathbf{g}^l \odot \hat{\mathbf{a}}^l + \mathbf{b}^l)hl=f(gl⊙a^l+bl)
合并公式(2)，(3)并忽略参数lll
，我们有：
h=f(gσ2+ϵ⊙(a−μ)+b)\mathbf{h} = f(\frac{\mathbf{g}}{\sqrt{\sigma^2 + \epsilon}} \odot (\mathbf{a} - \mu)+ \mathbf{b})h=f(σ2+ϵ​g​⊙(a−μ)+b)
2.2 RNN中的LN
在RNN中，我们可以非常简单的在每个时间片中使用LN，而且在任何时间片我们都能保证归一化统计量统计的是HHH
个节点的信息。对于RNN时刻ttt
时的节点，其输入是t−1t-1t−1
时刻的隐层状态ht−1h^{t-1}ht−1
和ttt
时刻的输入数据xt\mathbf{x}_txt​
，可以表示为：
at=Whhht−1+Wxhxt\mathbf{a}^t = W_{hh}h^{t-1} + W_{xh}\mathbf{x}^tat=Whh​ht−1+Wxh​xt
接着我们便可以在at\mathbf{a}^tat
上采取和1.1节中完全相同的归一化过程：
ht=f(g(σt)2+ϵ⊙(at−μt)+b)μt=1H∑i=1Haitσt=1H∑i=1H(ait−μt)2\mathbf{h}^t = f(\frac{\mathbf{g}}{\sqrt{(\sigma^t)^2 + \epsilon}} \odot (\mathbf{a}^t - \mu^t)+ \mathbf{b})
\qquad
\mu^t = \frac{1}{H}\sum_{i=1}^H a^t_i
\qquad
\sigma^t = \sqrt{\frac{1}{H} \sum_{i=1}^H(a_i^t
 - \mu^t)^2}ht=f((σt)2+ϵ​g​⊙(at−μt)+b)μt=H1​i=1∑H​ait​σt=H1​i=1∑H​(ait​−μt)2​
2.3 LN与ICS和损失平面平滑
LN能减轻ICS吗？当然可以，至少LN将每个训练样本都归一化到了相同的分布上。而在BN的文章中介绍过几乎所有的归一化方法都能起到平滑损失平面的作用。所以从原理上讲，LN能加速收敛速度的。
3. 对照实验
这里我们设置了一组对照试验来对比普通网络，BN以及LN在MLP和RNN上的表现。这里使用的框架是Keras：
3.1 MLP上的归一化
这里使用的是MNIST数据集，但是归一化操作只添加到了后面的MLP部分。Keras官方源码中没有LN的实现，我们可以通过pip install keras-layer-normalization进行安装，使用方法见下面代码
from keras_layer_normalization import LayerNormalization
​
# 构建LN CNN网络
model_ln = Sequential()
model_ln.add(Conv2D(input_shape = (28,28,1), filters=6, kernel_size=(5,5), padding='valid', activation='tanh'))
model_ln.add(MaxPool2D(pool_size=(2,2), strides=2))
model_ln.add(Conv2D(input_shape=(14,14,6), filters=16, kernel_size=(5,5), padding='valid', activation='tanh'))
model_ln.add(MaxPool2D(pool_size=(2,2), strides=2))
model_ln.add(Flatten())
model_ln.add(Dense(120, activation='tanh'))
model_ln.add(LayerNormalization()) # 添加LN运算
model_ln.add(Dense(84, activation='tanh'))
model_ln.add(LayerNormalization())
model_ln.add(Dense(10, activation='softmax'))
另外两个对照试验也使用了这个网络结构，不同点在于归一化部分。图3左侧是batchsize=128时得到的收敛曲线，从中我们可以看出BN和LN均能取得加速收敛的效果，且BN的效果要优于LN。图3右侧是batchsize=8是得到的收敛曲线，这时BN反而会减慢收敛速度，验证了我们上面的结论，对比之下LN要轻微的优于无归一化的网络，说明了LN在小尺度批量上的有效性。图3的完整代码见连接：https://github.com/senliuy/keras_layerNorm_mlp_lstm/blob/master/mnist.ipynb 。
 
图3：batchsize=128(左)和batchsize=8(右)损失收敛曲线示意图
3.2 LSTM上的归一化
另外一组对照实验是基于imdb的二分类任务，使用了glove作为词嵌入。这里设置了无LN的LSTM和带LN的LSTM的作为对照试验。LN_LSTM源码参考 https://github.com/cleemesser/keras-layer-norm-work 其网络结构如下面代码：
# https://github.com/cleemesser/keras-layer-norm-work
from lstm_ln import LSTM_LN
model_ln = Sequential()
​
model_ln.add(Embedding(max_features,100))
model_ln.add(LSTM_LN(128))
model_ln.add(Dense(1, activation='sigmoid'))
model_ln.summary()
从图4的实验结果我我们可以看出LN对于RNN系列动态网络的收敛加速上的效果是略有帮助的。LN的有点主要体现在两个方面：
1.LN得到的模型更稳定；
2.LN有正则化的作用，得到的模型更不容易过拟合。
至于论文中所说的加速收敛的效果，从我的实验上结果上看不到明显的加速。源码见：https://github.com/senliuy/keras_layerNorm_mlp_lstm/blob/master/imdb.ipynb 。
 
图4：训练集损失值(左)验证集准确率(右)示意图
3.3 CNN上的归一化
我们也尝试了将LN添加到CNN之后，实验结果发现LN破坏了卷积学习到的特征，模型无法收敛，所以在CNN之后使用BN是一个更好的选择。
4. 总结
LN是和BN非常近似的一种归一化方法，不同的是BN取的是不同样本的同一个特征，而LN取的是同一个样本的不同特征。在BN和LN都能使用的场景中，BN的效果一般优于LN，原因是基于不同数据，同一特征得到的归一化特征更不容易损失信息。
但是有些场景是不能使用BN的，例如batchsize较小或者在RNN中，这时候可以选择使用LN，LN得到的模型更稳定且起到正则化的作用。RNN能应用到小批量和RNN中是因为LN的归一化统计量的计算是和batchsize没有关系的。
最近更新 2yr ago