1. 摘要
在 R e L U ReLU ReLU 的基础上作者提出了 P R e L U PReLU PReLU,在几乎没有增加额外参数的前提下既可以提升模型的拟合能力,又能减小过拟合风险。
针对 R e L U / P R e L U ReLU/PReLU ReLU/PReLU 的矫正非线性,作者设计了一个鲁棒的的参数初始化方法。
2. 介绍
在过去几年,随着更强大网络模型的构建和有效防止过拟合策略的设计,我们在视觉识别任务上取得了长足的进步。一方面,由于网络复杂性的增加(比如增加深度、扩大宽度)、小步长的使用、新的激活函数和更成熟的网络层设计,神经网络变得更加能够拟合训练数据。另一方面,有效的正则化技术、数据增广和大规模的数据让网络获得了更好的泛化能力。
其中,激活函数 R e L U ReLU ReLU 是其中一个非常关键的因素,本文在此基础上做了两点主要改进。首先,我们提出了一个新的激活函数 P R e L U PReLU PReLU (Parametric Rectified Linear Unit),该激活函数可以自适应地学习矫正线性单元的参数,并且能够在增加可忽略的额外计算成本下提高准确率。其次,我们研究了模型训练的难度,得出了一种理论上合理的初始化方法,有助于深层网络模型的收敛。
3. PReLU
我们设计的激活函数定义如下:
其中, y i y_i yi 是非线性激活函数 f f f 在第 i i i 个通道的输入, a i a_i ai 负责控制负半轴的斜率。在这里,我们允许不同通道的激活函数不一样。当 a i = 0 a_i=0 ai=0 时, P R e L U PReLU PReLU 就变成了 R e L U ReLU ReLU,而且 a i a_i ai 是一个可以学习的参数。
如果 a i a_i ai 是一个小的固定值,那么 P R e L U PReLU PReLU 就变成了 L R e L U LReLU LReLU(Leaky ReLU)。 L R e L U LReLU LReLU 的动机是为了避免零梯度,实验表明 L R e L U LReLU LReLU 相较 R e L U ReLU ReLU 对准确率几乎没有提高。但是,作者希望模型能够自适应地学习 P R e L U PReLU PReLU 的参数,从而能够得到更专门的激活函数。
另外,作者也设计了一个通道共享的变体,也就是同一层网络所有通道共享同一个参数,这样,每一层仅仅引入了一个额外的参数。
在反向传播过程中, a i {a_i} ai 的更新公式可以通过链式法则求出, ε \varepsilon ε 代表目标函数。
然后,采用动量的方法来更新 a i {a_i} ai :
作者没有使用权重衰减(L2 正则化),因为这样会让 a i {a_i} ai 趋向于 0。在实验中,即使没有正则化,学习到的系数幅度也很少会大于 1。而且,对 a i {a_i} ai 的范围也未做限制,因此激活函数也可能是非单调的。
作者设计了一个简单的网络来验证 P R e L U PReLU PReLU 的有效性,学习到的系数如下表所示。
这其中,有两个有意思的现象。一,第一个卷积层的系数远远比 0 要大(0.681, 0.596),因为这一层的卷积核大部分是类 Gabor 过滤器,比如说边缘、纹理检测器,学习到的结果表明卷积核的正负响应都被接受。这在卷积核数量有限的情况下,可以被认为是一个更加经济地利用低层次信息的方式。二,对于通道独立的激活函数,较深的卷积层通常具有较小的系数。 这意味着激活函数在深度增加时逐渐变得“更加非线性”。 换句话说,学习模型倾向于在较早阶段保留更多信息,并在更深层次阶段变得更具辨别力。
针对通道共享的版本, P R e L U PReLU PReLU 仅仅引入了 13 个额外的参数,却取得了相对于基线情况 1.1% 的提升,这也表明了自适应学习激活函数形状的重要性。
4. 卷积核权重的初始化
近来,深层卷积神经网络的权重都是用正态分布来初始化的。若采用固定的标准差,非常深的模型很难收敛。对此,我们可以采取预训练的方式来初始化模型,但这样就需要更长的训练时间。先前的学者也提出了一个 X a v i e r Xavier Xavier 初始化方法,但它针对的是激活函数满足线性的情况, R e L U / P R e L U ReLU/PReLU ReLU/PReLU 却不满足这个特性。
4.1. 前向传播情况
针对一个卷积层,其前向传播为:
x x x 是一个 k 2 c ∗ 1 k^2c*1 k2c∗1 的向量,代表 c c c 个通道 k × k k×k k×k 区域的像素, k k k 是卷积核的大小, n = k 2 c n=k^2c n=k2c 代表一个激活值局部连接的像素个数。 W W W 是一个 d × n d×n d×n 的矩阵, d d d 是卷积核的个数, W W W 的每一行代表一个卷积核的权重。我们用 l l l 来标识某一层,则有 x l = f ( y l − 1 ) x_l=f(y_{l-1}) xl=f(yl−1), f f f 是激活函数,且有 c l = d l − 1 c_l=d_{l-1} cl=dl−1。
我们让 W l W_l Wl 初始化的元素互相独立且同分布,然后假设 x l x_l xl 中的元素也是互相独立并且同分布的,而且 x l x_l xl 与 W l W_l Wl 互相独立,那么有:
V a r [ X 1 + X 2 + ⋯ + X n ] = V a r [ X 1 ] + V a r [ X 2 ] + ⋯ + V a r [ X n ] Var[X_1+X_2+\cdots+X_n] = Var[X_1] + Var[X_2] + \cdots +Var[X_n] Var[X1+X2+⋯+Xn]=Var[X1]+Var[X2]+⋯+Var[Xn],如果 X 1 , X 2 , ⋯   , X n X_1,X_2,\cdots,X_n X1,X2,⋯,Xn 互相独立,这里每个 w l i ∗ x l i w_l^i*x_l^i wli∗xli 是互相独立的,总共对 n l n_l nl 项求和。
让 w l w_l wl 具有零均值,然后有:
V a r [ X ] = E [ X 2 ] − ( E [ X ] ) 2 Var[X]=E[X^2]-(E[X])^2 Var[X]=E[X2]−(E[X])2
V a r [ w l ] = E [ w l 2 ] − ( E [ w l ] ) 2 = E [ w l 2 ] Var[w_l] = E[w_l^2]-(E[w_l])^2=E[w_l^2] Var[wl]=E[wl2]−(E[wl])2=E[wl2]
E [ w l ] E [ x l ] = 0 ∗ E [ x l ] = 0 E[w_l]E[x_l] = 0*E[x_l]=0 E[wl]E[xl]=0∗E[xl]=0
V a r [ w l x l ] = E [ ( w l x l ) 2 ] − ( E [ w l x l ] ) 2 = E [ w l 2 ] E [ x l 2 ] − ( E [ w l ] E [ x l ] ) 2 = V a r [ w l ] E [ x l 2 ] Var[w_lx_l]=E[(w_lx_l)^2]-(E[w_lx_l])^2=E[w_l^2]E[x_l^2]-(E[w_l]E[x_l])^2=Var[w_l]E[x_l^2] Var[wlxl]=E[(wlxl)2]−(E[wlxl])2=E[wl2]E[xl2]−(E[wl]E[xl])2=Var[wl]E[xl2]
这里, E [ x l 2 ] ̸ = V a r [ x l ] E[x_l^2] \not = Var[x_l] E[xl2]̸=Var[xl],因为 R e L U ReLU ReLU 激活函数关于原点不对称, x l x_l xl 均值不为 0。
如果让 w l − 1 w_{l-1} wl−1 关于原点对称分布并且 b l − 1 = 0 b_{l-1}=0 bl−1=0,那么 y l − 1 y_{l-1} yl−1 也具有对称分布并且均值为 0。那么有:
V a r [ y l − 1 ] = E [ y l − 1 2 ] − ( E [ y l − 1 ] ) 2 = E [ y l − 1 2 ] Var[y_{l-1}] =E[y_{l-1}^2]-(E[y_{l-1}])^2=E[y_{l-1}^2] Var[yl−1]=E[yl−12]−(E[yl−1])2=E[yl−12]
经过 R e L U ReLU ReLU 激活函数后,一半的 y l − 1 y_{l-1} yl−1 变成了 0,所以有:
E [ x l 2 ] = 1 2 E [ y l − 1 2 ] = 1 2 V a r [ y l − 1 ] E[x_l^2]=\frac{1}{2}E[y_{l-1}^2]=\frac{1}{2}Var[y_{l-1}] E[xl2]=21E[yl−12]=21Var[yl−1]
代入前面的式子,可得:
把 L L L 层网络放在一起,则有:
这个乘积是初始化参数的关键。一个合适的初始化方式应该避免指数级减少或者放大输入信号的幅度,因此我们期望上面的乘积有一个合理的尺度,比如 1。
这就是一个零均值标准差为 2 / n l \sqrt{2/n_l} 2/nl 的正态分布,也就是作者提出的参数初始化方式,同时,初始化偏置 b = 0 b=0 b=0。
针对第一层网络,对输入我们没有应用激活函数,因此应该有 n 1 V a r [ w 1 ] = 1 n_1Var[w_1]=1 n1Var[w1]=1,但这只有一层,为了简化,所有层的参数我们都按照上面的式子来初始化。
4.2. 反向传播情况
同理,一个卷积层的反向传播可以表示为:
这里,为了简化, Δ x \Delta x Δx 和 Δ y \Delta y Δy 就代表相应变量的梯度。
Δ y \Delta y Δy 代表 d d d 个通道 k × k k×k k×k 区域的像素,被调整为了一个 k 2 d ∗ 1 k^2d*1 k2d∗1 的向量。 n ^ = k 2 d ̸ = n \hat n=k^2d\not =n n^=k2d̸=n, W ^ \hat W W^ 是一个 c × n ^ c×\hat n c×n^ 的矩阵, Δ x \Delta x Δx 是一个 c ∗ 1 c*1 c∗1 的向量,代表这一层每个像素的梯度值。
假设 w l w_l wl 和 Δ y l \Delta y_l Δyl 互相独立, w l w_l wl 关于原点对称且具有零均值,那么 Δ x l \Delta x_l Δxl 也具有零均值。同时,我们有:
Δ y l = f ′ ( y l ) Δ x l + 1 \Delta y_l=f'(y_l)\Delta x_{l+1} Δyl=f′(yl)Δxl+1
如果激活函数是 R e L U ReLU ReLU,那么 f ′ ( y l ) f'(y_l) f′(yl) 要么为 0 要么为 1,而且概率相等,因此有:
因子 1 2 \frac{1}{2} 21 的出现都是因为激活函数是 R e L U ReLU ReLU,把 L L L 层网络放在一起,则有:
如果梯度不会指数级增大或者缩小的话,需要有:
这是一个零均值标准差为 2 / n ^ l \sqrt{2/\hat n_l} 2/n^l 的正态分布,和前面唯一不同的是 n ^ l = k l 2 d l \hat n_l=k_l^2d_l n^l=kl2dl 而 n l = k l 2 c l = k l 2 d l − 1 n_l=k_l^2c_l=k_l^2d_{l-1} nl=kl2cl=kl2dl−1。针对第一层网络,对输入我们不需要求梯度,但同理为了简化,所有层的参数我们都采用一样的方式来初始化。
通过前向传播得到的式子和通过反向传播得到的式子都是可行的,因为二者只是相差了一个因子。也就是说如果初始化能够适当地缩放后向信号,那么前向信号也是如此; 反之亦然。对于作者实验的所有模型,两种形式都可以使它们收敛。
针对 P R e L U PReLU PReLU,我们则很容易得到:
5. 实验
5.1 和 X a v i e r Xavier Xavier 初始化方法的对比
针对一个 22 层的模型,用本文提出的初始化方式比 X a v i e r Xavier Xavier 能更快收敛。
针对一个 30 层的模型,用本文提出的初始化方式可以使模型收敛,但用 X a v i e r Xavier Xavier 则无法收敛。
5.2 R e L U ReLU ReLU 和 P R e L U PReLU PReLU 的对比
5.3 和其他模型的对比
获取更多精彩,请关注「seniusen」!