BP神经网络基础和数学原理
BP神经网络(Back Propagation, 反向传播), 是一种前馈型神经网络, 也是全连接层的结构. 感知器, BP神经网络, 卷积神经网络和径向基网络都是常见的前馈型神经网络, BP神经网络也是一种多层感知器.
基本概念
上图是一个简单的BP神经网络结构, 由输入层、隐藏层、输出层组成,其中每个绿色的圆圈是一个神经元, 蓝色的连接线表示数值的传递. 输入层有两个神经元, 表示可以传入一个二维向量, 最终输出一个二维向量. 这个神经网络可以应用于二维平面上样本点的分类, 输入的 2 维向量是点的坐标, 输出的 2 维向量表示分类信息, 例如第一维输出值为 1, 第二维输出 0, 表示这个点属于第一类, 第一维输出 0, 第二维输出 1, 表示这个点属于第二类. 以其中隐藏层1的第二个神经元为例, 其内部计算过程如下图:
从左向右看,
x
i
x_i
xi 表示由上一个神经元传入的数值. 权重
w
w
w 中,
i
i
i 表示前一层神经元的序号,
j
j
j 表示指向的神经元的序号, 上标
l
l
l 为层次. 偏置
b
j
b_j
bj 是属于这个神经元的,
j
j
j 表示这个神经元的序号. 经加权求和后得到
v
j
v_j
vj, 其中的
M
M
M 为前一层的神经元总数. 再将结果代入激活函数
φ
\varphi
φ, 得到
y
j
y_j
yj 作为此神经元的输出. (关于激活函数的内容, 请看后文)
公式
为了公式清晰,省略掉权重 w w w 的上标 l l l, 相对于输出层的神经元, 有以下公式:
输入:
v
j
=
∑
i
=
1
M
w
i
j
x
i
+
b
j
.
v_j=\sum_{i=1}^{M} w_{ij} x_i +b_j.
vj=i=1∑Mwijxi+bj.
输出:
y j = φ ( v j ) . y_j=\varphi(v_j). yj=φ(vj).
总误差(Total instantaneous error energy):
E = 1 2 ∑ j = 1 N ( d j − y j ) 2 . E=\frac{1}{2} \sum_{j=1}^{N} (d_j-y_j)^2. E=21j=1∑N(dj−yj)2.
N N N 表示输出层神经元总数. 误差函数 E E E 是优化的目标函数, 通过迭代地改变权重 w w w 和 偏置 b b b, 使得函数值最小化, 即使得输出结果与输出期望差距缩小.
梯度下降与反向传播算法
我们的目标是寻找合适的 w w w 和 b b b, 使得网络的输出与期望输出尽可能接近, 这样当我们输入一个未经过训练的样本点时, 便可以准确地做出分类. 为了求解, 在训练开始时随时生成 w w w 和 b b b, 随机生成的值可以满足均匀分布或者正态分布. 再应用反向传播和梯度下降算法, 不断修正 w w w 和 b b b , 使得总误差最小化. 现在, 我们设输出层神经元代号为 k k k, 隐藏层2 代号为 j j j, 隐藏层1 代号为 i i i. 下面分两种情形推导梯度公式.
情形一: 对于输出层神经元
由梯度下降算法, 求得总误差关于权值 w j k w_{jk} wjk 的梯度
∂
E
∂
w
j
k
=
∂
E
∂
y
k
∂
y
k
∂
v
k
∂
v
k
∂
w
j
k
=
∂
E
∂
v
k
x
j
.
\frac{\partial E}{\partial w_{j k}}=\frac{\partial E}{\partial y_{k}} \frac{\partial y_{k}}{\partial v_{k}} \frac{\partial v_{k}}{\partial w_{j k}}=\frac{\partial E}{\partial v_{k}} x_{j}.
∂wjk∂E=∂yk∂E∂vk∂yk∂wjk∂vk=∂vk∂Exj.
同理得到关于偏置
b
k
b_k
bk 的梯度
∂
E
∂
b
k
=
∂
E
∂
y
k
∂
y
k
∂
v
k
∂
v
k
∂
b
k
=
∂
E
∂
v
k
.
\frac{\partial E}{\partial b_{k}}=\frac{\partial E}{\partial y_{k}} \frac{\partial y_{k}}{\partial v_{k}} \frac{\partial v_{k}}{\partial b_{k}}=\frac{\partial E}{\partial v_{k}}.
∂bk∂E=∂yk∂E∂vk∂yk∂bk∂vk=∂vk∂E.
公式推导中我们保留了
∂
E
∂
v
k
\frac{\partial E}{\partial v_{k}}
∂vk∂E, 是为了导出局部梯度
δ
k
=
∂
E
∂
v
k
=
∂
E
∂
y
k
∂
y
k
∂
v
k
=
−
(
d
k
−
y
k
)
φ
′
(
v
k
)
.
\delta_{k}=\frac{\partial E}{\partial v_{k}}=\frac{\partial E}{\partial y_{k}} \frac{\partial y_{k}}{\partial v_{k}}=-\left(d_{k}-y_{k}\right) \varphi^{\prime}\left(v_{k}\right).
δk=∂vk∂E=∂yk∂E∂vk∂yk=−(dk−yk)φ′(vk).
关于局部梯度, 我的理解是将一个神经元所有输入的加权和看做一个整体, 权重
w
w
w 和
b
b
b 仅通过这个整体影响神经网络的其他部分.
则权值调整和偏置的调整分别为
Δ
w
j
k
=
−
η
∂
E
∂
w
j
k
=
−
η
δ
k
x
j
Δ
b
k
=
−
η
∂
E
∂
b
k
=
−
η
δ
k
.
\begin{array}{l} \Delta w_{j k}=-\eta \frac{\partial E}{\partial w_{j k}}=-\eta \delta_{k} x_{j} \\ \Delta b_{k}=-\eta \frac{\partial E}{\partial b_{k}}=-\eta \delta_{k}. \end{array}
Δwjk=−η∂wjk∂E=−ηδkxjΔbk=−η∂bk∂E=−ηδk.
其中,
η
\eta
η 是学习率, 或者称为步长, 根据实际情况, 一般取一个较小的数, 如 0.1, 0.01. 前面加个负号意味着沿着负梯度方向(即函数值下降的方向)搜索.
情形二: 对于隐藏层要分一种情况讨论的原因在于无法直接得出隐藏层的误差, 因为只有输出层才有期望的输出, 而隐藏层没有定义. 所以, 对总误差求偏导, 需要穿过输出层, 到达隐藏层, 或者可以理解为原先的复合函数内部又复合了一个函数.
(1) 按照情形一的推导有
∂
E
∂
w
i
j
=
∂
E
∂
v
k
∂
v
j
∂
w
i
j
=
∂
E
∂
v
j
x
i
.
\frac{\partial E}{\partial w_{i j}}=\frac{\partial E}{\partial v_{k}} \frac{\partial v_{j}}{\partial w_{i j}}=\frac{\partial E}{\partial v_{j}} x_{i}.
∂wij∂E=∂vk∂E∂wij∂vj=∂vj∂Exi.
(2) 根据网络结构图, 可以看出, 隐藏层2 中一个神经元输出值的改变会影响每个输出层神经元的误差, 则对于隐藏层2 的
d
e
l
t
a
j
delta_j
deltaj, 需要由每一个输出层神经元反向传播误差.
δ j = ∂ E ∂ v j = ∑ k = 1 N ∂ E ∂ v k ∂ v k ∂ v j = ∑ k = 1 N δ k ∂ v k ∂ v j = ∑ k = 1 N δ k ∂ v k ∂ y j ∂ y j ∂ v j = ∑ k = 1 N δ k w k j φ ′ ( v j ) . \begin{aligned} \delta_{j} &=\frac{\partial E}{\partial v_{j}} \\ &=\sum_{k=1}^{N} \frac{\partial E}{\partial v_{k}} \frac{\partial v_{k}}{\partial v_{j}} \\ &=\sum_{k=1}^{N} \delta_{k} \frac{\partial v_{k}}{\partial v_{j}} \\ &=\sum_{k=1}^{N} \delta_{k} \frac{\partial v_{k}}{\partial y_{j}} \frac{\partial y_{j}}{\partial v_{j}} \\ &=\sum_{k=1}^{N} \delta_{k} w_{k j} \varphi^{\prime}\left(v_{j}\right). \end{aligned} δj=∂vj∂E=k=1∑N∂vk∂E∂vj∂vk=k=1∑Nδk∂vj∂vk=k=1∑Nδk∂yj∂vk∂vj∂yj=k=1∑Nδkwkjφ′(vj).
(3) 得
Δ w i j = − η ∂ E ∂ w i j = − η δ j x i Δ b j = − η ∂ E ∂ b j = − η δ j . \begin{aligned} &\Delta w_{i j}=-\eta \frac{\partial E}{\partial w_{i j}}=-\eta \delta_{j} x_{i}\\ &\Delta b_{j}=-\eta \frac{\partial E}{\partial b_{j}}=-\eta \delta_{j} \end{aligned}. Δwij=−η∂wij∂E=−ηδjxiΔbj=−η∂bj∂E=−ηδj.
可以注意到, 对于不同层的神经元只要求出对应的 δ \delta δ 即可, 在后面要讲的卷积神经网络同样沿用这种推到思路.
激活函数
多层的网络结构, 是一种线性变换, 而使用非线性的激活函数, 可以加入非线性因素, 增强模型的表达能力.
Logistic函数:
φ
(
x
)
=
1
1
+
e
−
x
.
\varphi(x)=\frac{1}{1+e^{-x}}.
φ(x)=1+e−x1.
对函数求导得:
φ ′ ( x ) = − − e − x ( 1 + e − x ) 2 = 1 1 + e − x e − x 1 + e − x = 1 1 + e − x 1 + e − x − 1 1 + e − x = 1 1 + e − x ( 1 − 1 1 + e − x ) = φ ( x ) ( 1 − φ ( x ) ) . \begin{aligned} \varphi^{\prime}(x)=&-\frac{-e^{-x}}{\left(1+e^{-x}\right)^{2}} \\ &=\frac{1}{1+e^{-x}} \frac{e^{-x}}{1+e^{-x}} \\ &=\frac{1}{1+e^{-x}} \frac{1+e^{-x}-1}{1+e^{-x}} \\ &=\frac{1}{1+e^{-x}}\left(1-\frac{1}{1+e^{-x}}\right) \\ &=\varphi(x)(1-\varphi(x)). \end{aligned} φ′(x)=−(1+e−x)2−e−x=1+e−x11+e−xe−x=1+e−x11+e−x1+e−x−1=1+e−x1(1−1+e−x1)=φ(x)(1−φ(x)).
Logistic 函数可以将非常大的输入值域映射到 0 到 1 的小范围内. 另外, 可以看到 logistic函数的导数可以用它本身表示, 这便于计算导数. 函数中的 e − x e^{-x} e−x 有时被写为 e − k x e^{-kx} e−kx , 其中 k k k 为正的常数, 可改变函数曲线的陡峭程度.
双曲正切函数:
tanh
(
x
)
=
e
x
−
e
−
x
e
x
+
e
−
x
.
\tanh(x)=\frac{e^{x}-e^{-x}}{e^{x}+e^{-x}}.
tanh(x)=ex+e−xex−e−x.
求导得:
tanh
′
=
1
−
tanh
2
(
x
)
.
\tanh^{\prime}=1-\tanh^{2}(x).
tanh′=1−tanh2(x).
Rectifier 函数:
φ
(
x
)
=
max
(
0
,
x
)
.
\varphi(x)=\max(0,x).
φ(x)=max(0,x).
如果一个神经元部署了这个激活函数, 则这个神经元被称为 ReLU (rectified linear unit).
有时也使用一个平滑的 softplus 函数逼近 Rectifier 函数.
φ
(
x
)
=
ln
(
1
+
e
x
)
.
\varphi(x)=\ln(1+e^x).
φ(x)=ln(1+ex).
这两个函数的图像:
误差函数
从上图中可以发现, 当
x
x
x 变大或变小时, 曲线趋于平缓, 即梯度变小, 由导出的公式可以看出, 即使有较大的误差, 也会使得
D
e
l
t
a
w
Delta w
Deltaw,
D
e
l
t
a
b
Delta b
Deltab 很小, 这会导致权重
w
w
w 和偏置
b
b
b 更新缓慢. 这个问题源于我们选择的误差函数是方差代价函数. 即
E
j
=
1
2
(
d
j
−
y
j
)
2
.
E_j=\frac{1}{2}(d_j-y_j)^{2}.
Ej=21(dj−yj)2.
为了解决上面遇到的问题, 引入交叉熵代价函数( Cross-entropy cost function ), 应用于输出层, 原有的激活函数可保留下来.
交叉熵代价函数:
C
=
−
∑
j
=
1
N
(
d
j
ln
y
j
+
(
1
−
d
j
)
ln
(
1
−
y
j
)
)
.
C=-\sum_{j=1}^{N}\left(d_{j} \ln y_{j}+\left(1-d_{j}\right) \ln \left(1-y_{j}\right)\right).
C=−j=1∑N(djlnyj+(1−dj)ln(1−yj)).
其中
N
N
N 是维度。
按照前面的推导过程, 以情形一为例, 推导出代价函数 C 关于 w ( i j ) w_(ij) w(ij) 的偏导数:
∂ C ∂ w i j = ∂ C ∂ y j ∂ y j ∂ v j ∂ v j ∂ w i j = − ∑ j = 1 N ( d j y j − 1 − d j 1 − y j ) y j ′ y i = − ∑ j = 1 N d j ( 1 − y j ) − ( 1 − d j ) y j y j ( 1 − y j ) y j ( 1 − y j ) y i = − ∑ j = 1 N ( d j − y j ) y i . \begin{array}{l}{\frac{\partial C}{\partial w_{i j}}=\frac{\partial C}{\partial y_{j}} \frac{\partial y_{j}}{\partial v_{j}} \frac{\partial v_{j}}{\partial w_{i j}}} \\ {=-\sum_{j=1}^{N}\left(\frac{d_{j}}{y_{j}}-\frac{1-d_{j}}{1-y_{j}}\right) y_{j}^{\prime} y_{i}} \\ {=-\sum_{j=1}^{N} \frac{d_{j}\left(1-y_{j}\right)-\left(1-d_{j}\right) y_{j}}{y_{j}\left(1-y_{j}\right)} y_{j}\left(1-y_{j}\right) y_{i}} \\ {=-\sum_{j=1}^{N}\left(d_{j}-y_{j}\right) y_{i}}.\end{array} ∂wij∂C=∂yj∂C∂vj∂yj∂wij∂vj=−∑j=1N(yjdj−1−yj1−dj)yj′yi=−∑j=1Nyj(1−yj)dj(1−yj)−(1−dj)yjyj(1−yj)yi=−∑j=1N(dj−yj)yi.
结果中不包含激活函数 y j y_j yj 的导数, 而是期望值与激活函数值的差, 则当误差大时, 权重更新快, 误差小时, 权重更新慢, 这是我们期望达到的效果.
另外, 对于一个多分类的问题, 我们希望输出的结果可以被视为属于某一类的概率, 即将输出向量转换后, 使得各维度和为 1. 这里引入 Softmax函数, 或者叫归一化指数函数(Normalized exponential function):
σ ( z ) i = e z i ∑ j = 1 N e z j , i = 1 , … , N . \sigma(z)_{i}=\frac{e^{z_{i}}}{\sum_{j=1}^{N} e^{z_{j}}}, i=1, \ldots, N. σ(z)i=∑j=1Nezjezi,i=1,…,N.
其中, z z z 是一个 N N N 维的向量. 函数值是向量中每一维归一化后的值.
反向传播和梯度下降算法进一步讨论
(1) 反向传播算法实现了一种对可能的权重空间的梯度下降搜索, 它迭代地减小网络输出值和期望值的差距. 对于多层的网络, 误差函数可能有多个局部极小值, 梯度下降算法可能会陷入某个局部极小值, 而不能保证达到全局最小. 但是在实际应用中, 人们发现这些局部极小值仍然可以较好地符合期望. 可以这样直观的理解这个问题: 由于多层的网络含有很多的权重, 即函数的维度很高, 当陷入一个局部极小值时, 对于其它的权重, 这未必是局部极小值(因为关于某个权重的偏导数与其他权重无关), 而有可能离开这个局部极小值.
(2) 在实际应用中, 训练一个 BP 神经网络有两种方式: 一是在线学习, 即随机梯度下降算法(SGD), 方法是一次输入一个样本点, 计算出总误差, 运用反向传播算法, 调整权重和偏置. 二是批量学习, 即批梯度下降算法, 方法是先后输入多个样本点, 累计输出的总误差, 最后求出每个数据的总误差平均值, 以这个平均值, 运用反向传播算法, 调整权重和偏置. 每次输入的多个数据叫做一个 batch(一批), 数据的数目叫做 batch的大小, 每进行一次这样的过程叫做一个epoch(回合). 下面对比来看这两种方法:
- 随机梯度下降: 对每个样本点做一次反向传播, 由于对单个样本点来说, 它在空间中的分布具有不确定性, 所以梯度的搜索方向具有一定的随机性, 从而不容易陷入局部极小值. 其次, 它能够追踪到训练数据很小的改变, 尤其是数据产生于不稳定的环境下.
- 批梯度下降: 由于一次对一批样本点的学习, 在简单的情况下, 能够相对较快地到达局部极小值. 从统计的角度来看, 批梯度下降可以看成某种形式的统计推断(个人理解为中心极限定理), 适合于解非线性回归问题.
(3) 避免局部最小值的启发式方法:
对权重的更新增加一个冲量项. 即分别在 Δ w , Δ b \Delta w, \Delta b Δw,Δb 两式后面加一项 α Δ w i j ( n − 1 ) \alpha \Delta w_{ij}(n-1) αΔwij(n−1), 其中 n − 1 n-1 n−1 表示上一次迭代时的调整量, α \alpha α 称为冲量常数, 加上这一项后, 相当于增加下一次权重调整的步长, 从而可能一次越过某一局部极小值, 到达下一个更小的极小值.
使用同样的数据训练多个网络, 但用不同的初始化方法初始化权重. 选出最优的训练结果. 或者保留训练出来的多个网络模型, 在进行预测时, 计算多个网络的加权平均值做为最终的结果.
应用举例
首先, 由左边的神经网络结构图可以看出, 输入一个 2 维的向量, 经过一个隐藏层, 输出 10 维的向量. 这个网络学习一个语音识别的模型. 首先根据单词发音的声波, 进行频谱分析(快速傅里叶变换), 选出两个振幅最大的频率, 作为输入向量, 右图是将输入向量映射到二维平面的结果, 经过学习, 可以得到一个复杂的非线性分类平面.