深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)

  • 作者:韩信子@ShowMeAI
  • 教程地址:https://www.showmeai.tech/tutorials/35
  • 本文地址:https://www.showmeai.tech/article-detail/217
  • 声明:版权所有,转载请联系平台与作者并注明出处
  • 收藏ShowMeAI查看更多精彩内容

第2门课 改善深层神经网络:超参数调试、正则化以及优化,第2周: 优化算法

本系列为吴恩达老师《深度学习专项课程(Deep Learning Specialization)》学习与总结整理所得,对应的课程视频可以在这里查看。

引言

在ShowMeAI前一篇文章 深度学习的实用层面 中我们对以下内容进行了介绍:

  • Train / Dev / Test sets的切分和比例选择
  • Bias和Variance的相关知识
  • 防止过拟合的方法:L2正则化和Dropout
  • 规范化输入以加快梯度下降速度和精度
  • 梯度消失和梯度爆炸的原因及处理方法
  • 梯度检查

本篇内容展开介绍深度神经网络中的一些优化算法,通过使用这些技巧和方法来提高神经网络的训练速度和精度。

1.Batch梯度下降法

Mini-batch梯度下降 Mini-batch Gradient Descent

Batch梯度下降法(批梯度下降法)是最常用的梯度下降形式,它是基于整个训练集的梯度下降算法,在更新参数时使用所有的样本来进行更新。

对整个训练集进行梯度下降法的时候,我们必须处理整个训练数据集,然后才能进行一步梯度下降,即每一步梯度下降法需要对整个训练集进行一次处理,如果训练数据集很大的时候,处理速度就会比较慢。

但是如果每次处理训练数据的一部分,基于这个子集进行梯度下降法,算法迭代速度会更快。而处理的这些一小部分训练子集即称为Mini-Batch,这个算法也就是我们说的Mini-Batch梯度下降法

2.Mini-Batch梯度下降法

理解 Mini-batch 梯度下降 Understanding Mini-batch Gradient Descent

Mini-Batch梯度下降法(小批量梯度下降法)每次同时处理单个的Mini-Batch,其他与Batch梯度下降法一致。

使用Batch梯度下降法,对整个训练集的一次遍历只能做一个梯度下降;而使用Mini-Batch梯度下降法,对整个训练集的一次遍历(称为一个epoch)能做Mini-Batch个数个梯度下降。之后,可以一直遍历训练集,直到最后收敛到一个合适的精度。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第1张图片

Batch梯度下降法和Mini-Batch梯度下降法代价函数的变化趋势如上图所示:

  • 使用Batch gradient descent,随着迭代次数增加,cost是不断减小的。
  • 使用Mini-batch gradient descent,随着在不同的mini-batch上迭代训练,cost并不是单调下降,而是振荡下降的,最终也能得到较低的cost值。出现细微振荡的原因是不同的mini-batch之间是有差异的。例如可能第一个子集 ( X { 1 } , Y { 1 } ) (X^{\{1\}},Y^{\{1\}}) (X{1},Y{1})是好的子集,而第二个子集 ( X { 2 } , Y { 2 } ) (X^{\{2\}},Y^{\{2\}}) (X{2},Y{2})包含了一些噪声noise。出现细微振荡是正常的。

2.1 Batch大小及影响

我们在训练神经网络的时候,使用mini-batch gradient descent,经常要指定一个batch批次的样本数量。而不同的batch大小会影响训练的过程,其中有2个特例,mini-batch gradient descent会退化为不同的算法:

  • Mini-Batch的大小为1,即是随机梯度下降法(stochastic gradient descent),每个样本都是独立的Mini-Batch。
  • Mini-Batch的大小为 m m m(数据集大小),即是Batch梯度下降法。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第2张图片

如上图,我们对比一下Batch gradient descentStachastic gradient descent的梯度下降曲线。

  • 图中蓝色的线代表Batch gradient descent。Batch gradient descent会比较平稳地接近全局最小值,但是因为使用了所有m个样本,每次前进的速度有些慢。
  • 图中紫色的线代表Stochastic gradient descent。Stochastic gradient descent每次前进速度很快,但是路线曲折,有较大的振荡,最终会在最小值附近来回波动,难以真正达到最小值处。而且在数值处理上就不能使用向量化的方法来提高运算速度。

(1) Batch梯度下降法(Batch gradient descent)

  • 对所有 m 个训练样本执行一次梯度下降,每一次迭代时间较长,训练过程慢
  • 相对噪声低一些,幅度也大一些。
  • 成本函数总是向减小的方向下降。

(2) 随机梯度下降法(Stochastic gradient descent)

  • 对每一个训练样本执行一次梯度下降,训练速度快,但丢失了向量化带来的计算加速
  • 有很多噪声,减小学习率可以适当。
  • 成本函数总体趋势向全局最小值靠近,但永远不会收敛,而是一直在最小值附近波动。

(3) Mini-Batch gradient descent

实际使用中,batch size不能设置得太大(会倾向于Batch gradient descent),也不能设置得太小(倾向于Stochastic gradient descent)。

选择一个1的合适的大小进行Mini-Batch梯度下降,可以实现快速学习,也应用了向量化带来的好处,且成本函数的下降处于前两者之间。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第3张图片

mini-batch gradient descent的梯度下降曲线如图绿色曲线所示,每次前进速度较快,且振荡较小,基本能接近全局最小值

2.2 Batch大小的选择

吴恩达老师也给出了一些关于batch大小选择的经验:

  • 训练样本量小(如 m ≤ 2000 m \le 2000 m2000),选择Batch梯度下降法。
  • 训练样本量大,选择Mini-Batch梯度下降法。
  • 与计算机的信息存储方式相适应,代码在Batch大小为2的幂次时运行要快一些,典型的大小为 2 6 2^6 26 2 7 2^7 27、…、 2 9 2^9 29
  • Batch的大小要匹配CPU/GPU内存。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第4张图片

Batch的大小是重要的超参数,需要根据经验快速尝试,找到能够最有效地减少成本函数的值。

2.3 获得Mini-Batch的步骤

前面提到了batch大小的选择方法,当我们确定batch大小后,在应用mini-batch梯度下降算法时,可以通过以下方式获得1个Batch的数据:

  • 将数据集打乱
  • 按照既定的大小分割数据集

其中打乱数据集的代码:

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第5张图片

# 获得样本数量
m = X.shape[1] 
# 对m个样本进行乱序
permutation = list(np.random.permutation(m))
# 取出洗牌之后的样本特征和标签
shuffled_X = X[:, permutation]
shuffled_Y = Y[:, permutation].reshape((1,m))

(上述python代码使用到numpy工具库,想了解更多的同学可以查看ShowMeAI的 图解数据分析 系列中的numpy教程,也可以通过ShowMeAI制作的 numpy速查手册 快速了解其使用方法)

代码解读

np.random.permutationnp.random.shuffle有两处不同:

  • 如果传给permutation一个矩阵,它会返回一个洗牌后的矩阵副本;而shuffle只是对一个矩阵进行洗牌,没有返回值。
  • 如果传入一个整数,它会返回一个洗牌后的arange

2.4 符号表示

在进一步讲解优化算法之前,我们来对数学标记做一个统一和说明:

  • 我们使用小括号上标 i i i表示训练集里的值, x ( i ) x^{(i)} x(i)是第 i i i个训练样本。
  • 我们使用中括号上标 l l l表示神经网络的层数, z [ l ] z^{[l]} z[l]表示神经网络中第 l l l层的 z z z值。
  • 我们使用上标 t t t来代表不同的Batch数据,即 X t X^{t} Xt Y t Y^{t} Yt

3.指数加权平均

指数加权平均 Exponentially Weighted Averages

下面我们将介绍指数加权平均(Exponentially weighted averages)的概念。

举个例子,记录半年内伦敦市的气温变化,并在二维平面上绘制出来,如下图所示:

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第6张图片

看上去,温度数据似乎有noise,而且抖动较大。如果我们希望看到半年内气温的整体变化趋势,可以通过「移动平均」(moving average)的方法来对每天气温进行平滑处理。

例如我们可以设 V 0 = 0 V_0=0 V0=0,当成第0天的气温值。

第一天的气温与第0天的气温有关:

V 1 = 0.9 V 0 + 0.1 θ 1 V_1=0.9V_0+0.1\theta_1 V1=0.9V0+0.1θ1

第二天的气温与第一天的气温有关:

V 2 = 0.9 V 1 + 0.1 θ 2 = 0.9 ( 0.9 V 0 + 0.1 θ 1 ) + 0.1 θ 2 = 0. 9 2 V 0 + 0.9 ⋅ 0.1 θ 1 + 0.1 θ 2 \begin{aligned} V_2 &=0.9V_1+0.1\theta_2\\ &=0.9(0.9V_0+0.1\theta_1)+0.1\theta_2\\ &=0.9^2V_0+0.9\cdot0.1\theta_1+0.1\theta_2 \end{aligned} V2=0.9V1+0.1θ2=0.9(0.9V0+0.1θ1)+0.1θ2=0.92V0+0.90.1θ1+0.1θ2

第三天的气温与第二天的气温有关:

V 3 = 0.9 V 2 + 0.1 θ 3 = 0.9 ( 0. 9 2 V 0 + 0.9 ⋅ 0.1 θ 1 + 0.1 θ 2 ) + 0.1 θ 3 = 0. 9 3 V 0 + 0. 9 2 ⋅ 0.1 θ 1 + 0.9 ⋅ 0.1 θ 2 + 0.1 θ 3 \begin{aligned} V_3 &= 0.9V_2+0.1\theta_3\\ &= 0.9(0.9^2V_0+0.9\cdot0.1\theta_1+0.1\theta_2)+0.1\theta_3\\ &= 0.9^3V_0+0.9^2\cdot 0.1\theta_1+0.9\cdot 0.1\theta_2+0.1\theta_3 \end{aligned} V3=0.9V2+0.1θ3=0.9(0.92V0+0.90.1θ1+0.1θ2)+0.1θ3=0.93V0+0.920.1θ1+0.90.1θ2+0.1θ3

即第 t t t天与第 t − 1 t-1 t1天的气温迭代关系为:

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第7张图片

经过「移动平均」(moving average)处理得到的气温如下图红色曲线所示:

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第8张图片

这种滑动平均算法称为指数加权平均(exponentially weighted average)。根据前面的例子,我们可以看到它的推导公式一般形式为: V t = β V t − 1 + ( 1 − β ) θ t V_t=\beta V_{t-1}+(1-\beta)\theta_t Vt=βVt1+(1β)θt

其中指数加权平均的天数由 β \beta β值决定,近似表示为 1 1 − β \frac{1}{1-\beta} 1β1。上面的例子中:

  • β = 0.9 \beta=0.9 β=0.9,则 1 1 − β = 10 \frac{1}{1-\beta}=10 1β1=10,表示将前10天进行指数加权平均。
  • β = 0.98 \beta=0.98 β=0.98,则 1 1 − β = 50 \frac{1}{1-\beta}=50 1β1=50,表示将前50天进行指数加权平均。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第9张图片

β \beta β值越大,则指数加权平均的天数越多,平均后的趋势线就越平缓,但是同时也会向右平移。上图中绿色曲线和橙色曲线分别表示了 β = 0.98 \beta=0.98 β=0.98 β = 0.5 \beta=0.5 β=0.5时,指数加权平均的结果。

公式解释

这里的 1 1 − β \frac{1}{1-\beta} 1β1是怎么来的呢?就标准数学公式来说,指数加权平均算法跟之前所有天的数值都有关系。

但是指数是衰减的,一般认为衰减到 1 e \frac1e e1就可以忽略不计了。因此,根据之前的推导公式,我们只要证明 β 1 1 − β = 1 e \beta^{\frac{1}{1-\beta}}=\frac1e β1β1=e1就好了。

1 1 − β = N \frac{1}{1-\beta}=N 1β1=N N > 0 N>0 N>0,则 β = 1 − 1 N \beta=1-\frac{1}{N} β=1N1 1 N < 1 \frac1N<1 N1<1。即证明转化为 ( 1 − 1 N ) N = 1 e (1-\frac1N)^N=\frac1e (1N1)N=e1

显然,当 N > > 0 N>>0 N>>0时,上述等式是近似成立的。这就简单解释了为什么指数加权平均的天数的计算公式为 1 1 − β \frac{1}{1-\beta} 1β1

综上,**指数加权平均(Exponentially Weight Average)**是一种常用的序列数据处理方式,计算公式为:

S t = { Y 1 , t = 1 β S t − 1 + ( 1 − β ) Y t , t > 1 S_t = \begin{cases} Y_1, &t = 1 \\ \beta S_{t-1} + (1-\beta)Y_t, &t > 1 \end{cases} St={Y1,βSt1+(1β)Yt,t=1t>1

其中 Y t Y_t Yt t t t下的实际值, S t S_t St t t t下加权平均后的值, β \beta β为权重值。

指数加权平均数在统计学中被称为“指数加权移动平均值”。

3.1 理解指数平均加权

理解指数加权平均 Understanding Exponentially Weighted Averages

我们将指数加权平均公式的一般形式写下来:

V t = β V t − 1 + ( 1 − β ) θ t = ( 1 − β ) θ t + ( 1 − β ) ⋅ β ⋅ θ t − 1 + ( 1 − β ) ⋅ β 2 ⋅ θ t − 2 + ⋯ + ( 1 − β ) ⋅ β t − 1 ⋅ θ 1 + β t ⋅ V 0 \begin{aligned} V_t &=\beta V_{t-1}+(1-\beta)\theta_t\\ & =(1-\beta)\theta_t+(1-\beta)\cdot\beta\cdot\theta_{t-1}+(1-\beta)\cdot \beta^2\cdot\theta_{t-2}+\cdots+(1-\beta)\cdot \beta^{t-1}\cdot \theta_1+\beta^t\cdot V_0 \end{aligned} Vt=βVt1+(1β)θt=(1β)θt+(1β)βθt1+(1β)β2θt2++(1β)βt1θ1+βtV0

观察上述推导得到的计算公式,其中:

  • θ t \theta_t θt, θ t − 1 \theta_{t-1} θt1, θ t − 2 \theta_{t-2} θt2,…, θ 1 \theta_1 θ1是原始数据值。
  • ( 1 − β ) (1-\beta) (1β), ( 1 − β ) β (1-\beta)\beta (1β)β, ( 1 − β ) β 2 (1-\beta)\beta^2 (1β)β2,…, ( 1 − β ) β t − 1 (1-\beta)\beta^{t-1} (1β)βt1是类似指数曲线,从右向左,呈指数下降的。

如果我们把每个时间点的 θ \theta θ和衰减指数写成向量形式,则最终指数加权平均结果 V t V_t Vt相当于两者的点乘。将原始数据值与衰减指数点乘,相当于做了指数衰减,随距离越远衰减越厉害(注意到 β \beta β小于1),有如下结论:

  • 离得越近的数据点,影响越大,离得越远的数据点,影响越小。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第10张图片

β = 0.9 \beta = 0.9 β=0.9时,

v 100 = 0.9 v 99 + 0.1 θ 100 v_{100} = 0.9v_{99} + 0.1 \theta_{100} v100=0.9v99+0.1θ100

v 99 = 0.9 v 98 + 0.1 θ 99 v_{99} = 0.9v_{98} + 0.1 \theta_{99} v99=0.9v98+0.1θ99

v 98 = 0.9 v 97 + 0.1 θ 98 v_{98} = 0.9v_{97} + 0.1 \theta_{98} v98=0.9v97+0.1θ98

展开:

v 100 = 0.1 θ 100 + 0.1 ∗ 0.9 θ 99 + 0.1 ∗ ( 0.9 ) 2 θ 98 + … v_{100} = 0.1 \theta_{100} + 0.1 * 0.9 \theta_{99} + 0.1 * {(0.9)}^2 \theta_{98} + \dots v100=0.1θ100+0.10.9θ99+0.1(0.9)2θ98+

其中, θ i \theta_i θi指第 i i i天的实际数据。所有 θ \theta θ前面的系数(不包括0.1)相加起来为1或者接近于1,这些系数被称作偏差修正(Bias Correction)

根据函数极限的一条定理:

lim ⁡ β → 0 ( 1 − β ) 1 β = 1 e ≈ 0.368 {\lim_{\beta\to 0}}(1 - \beta)^{\frac{1}{\beta}} = \frac{1}{e} \approx 0.368 β0lim(1β)β1=e10.368

β = 0.9 \beta = 0.9 β=0.9时,可以当作把过去10天的气温指数加权平均作为当日的气温,因为10天后权重已经下降到了当天的1/3左右。同理,当 β = 0.98 \beta = 0.98 β=0.98时,可以把过去50天的气温指数加权平均作为当日的气温。

因此,在计算当前时刻的平均值时,只需要前一天的平均值和当前时刻的值。

v t = β v t − 1 + ( 1 − β ) θ t v_t = \beta v_{t-1} + (1 - \beta)\theta_t vt=βvt1+(1β)θt

在实际代码中,只需要不断迭代赋值更新 v v v即可:

v : = β v + ( 1 − β ) θ t v := \beta v + (1 - \beta)\theta_t v:=βv+(1β)θt

指数平均加权并不是最精准的计算平均数的方法,你可以直接计算过去10天或50天的平均值来得到更好的估计,但缺点是保存数据需要占用更多内存,执行更加复杂,计算成本更加高昂。

指数加权平均数公式的好处之一在于它只需要一行代码,且占用极少内存,因此效率极高,且节省成本。

3.2 指数平均加权的偏差修正

指数加权平均的偏差修正 Bias Correction in Exponentially Weighted Averages

β = 0.98 \beta=0.98 β=0.98时,前面提到的气温示例,指数加权平均结果如绿色曲线。但实际上真实曲线如紫色曲线所示:

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第11张图片

紫色曲线与绿色曲线的区别是,紫色曲线开始的时候相对较低一些。因为开始时设置 v 0 = 0 v_0 = 0 v0=0,所以初始值会相对小一些,直到后面受前面的影响渐渐变小,趋于正常。

修正这种问题的方法是进行偏移校正(bias correction),即在每次计算完 v t v_t vt后,对 v t v_t vt进行下式处理:

V t = V t 1 − β t {V_t}=\frac{V_t}{1-\beta^t} Vt=1βtVt

换算到迭代公式中,即有 v t = β v t − 1 + ( 1 − β ) θ t 1 − β t v_t = \frac{\beta v_{t-1} + (1 - \beta)\theta_t}{{1-\beta^t}} vt=1βtβvt1+(1β)θt

观察上式:随着 t t t的增大, β \beta β t t t次方趋近于0。因此当 t t t很大的时候,偏差修正几乎没有作用,但是在前期学习可以帮助更好的预测数据。

4.动量梯度下降法

Momentum梯度下降 Gradient Descent with Momentum

4.1 从指数加权平均到动量梯度下降

大家已经了解了指数加权平均,现在我们回到神经网络优化算法,介绍一下动量梯度下降算法,其速度要比传统的梯度下降算法快很多。做法是在每次训练时,计算梯度的指数加权平均数,并利用该值来更新权重 W W W和常数项 b b b

具体过程为:for l = 1, .. , L

v d W [ l ] = β v d W [ l ] + ( 1 − β ) d W [ l ] v_{dW^{[l]}} = \beta v_{dW^{[l]}} + (1 - \beta) dW^{[l]} vdW[l]=βvdW[l]+(1β)dW[l]

v d b [ l ] = β v d b [ l ] + ( 1 − β ) d b [ l ] v_{db^{[l]}} = \beta v_{db^{[l]}} + (1 - \beta) db^{[l]} vdb[l]=βvdb[l]+(1β)db[l]

W [ l ] : = W [ l ] − α v d W [ l ] W^{[l]} := W^{[l]} - \alpha v_{dW^{[l]}} W[l]:=W[l]αvdW[l]

b [ l ] : = b [ l ] − α v d b [ l ] b^{[l]} := b^{[l]} - \alpha v_{db^{[l]}} b[l]:=b[l]αvdb[l]

其中,将动量衰减参数 β \beta β设置为0.9是超参数的一个常见且效果不错的选择。当 β \beta β被设置为0时,显然就成了Batch梯度下降法。

4.2 梯度下降 vs 动量梯度下降

我们用下图来对比一下优化算法的优化过程

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第12张图片

图中:

  • 蓝色曲线:使用一般的梯度下降的优化过程,由于存在上下波动,减缓了梯度下降的速度,因此只能使用一个较小的学习率进行迭代。
  • 紫色曲线:使用一般梯度下降+较大的学习率,结果可能偏离函数的范围。
  • 红色曲线:使用动量梯度下降,通过累加过去的梯度值来减少抵达最小值路径上的波动,加速了收敛,因此在横轴方向下降得更快。

当前后梯度方向一致时,动量梯度下降能够加速学习;而前后梯度方向不一致时,动量梯度下降能够抑制震荡。

另外,在10次迭代之后,移动平均已经不再是一个具有偏差的预测。因此实际在使用梯度下降法或者动量梯度下降法时,不会同时进行偏差修正。

补充:在其它文献资料中,动量梯度下降还有另外一种写法:

V d W = β V d W + d W V_{dW}=\beta V_{dW}+dW VdW=βVdW+dW

V d b = β V d b + d b V_{db}=\beta V_{db}+db Vdb=βVdb+db

即消去了 d W dW dW d b db db前的系数 ( 1 − β ) (1-\beta) (1β)。这样简化了表达式,但是学习因子 α \alpha α相当于变成了 α 1 − β \frac{\alpha}{1-\beta} 1βα,表示 α \alpha α也受 β \beta β的影响。从效果上来说,这种写法也是可以的,但是不够直观,且调参涉及到 α \alpha α,不够方便。所以,实际应用中,推荐第一种动量梯度下降的表达式。

动量梯度下降法的形象解释

将成本函数想象为一个碗状,从顶部开始运动的小球向下滚,其中 d w dw dw d b db db想象成球的加速度;而 v d w v_{dw} vdw v d b v_{db} vdb相当于速度。

小球在向下滚动的过程中,因为加速度的存在速度会变快,但是由于 β \beta β的存在,其值小于1,可以认为是摩擦力,所以球不会无限加速下去。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第13张图片

5.RMSProp 算法

RMSprop—— Root Mean Square Prop

RMSProp(Root Mean Square Propagation,均方根传播)是另外一种优化梯度下降速度的算法,它在对梯度进行指数加权平均的基础上,引入平方和平方根。具体过程为(省略了 l l l):

s d w = β s d w + ( 1 − β ) ( d w ) 2 s_{dw} = \beta s_{dw} + (1 - \beta)(dw)^2 sdw=βsdw+(1β)(dw)2

s d b = β s d b + ( 1 − β ) ( d b ) 2 s_{db} = \beta s_{db} + (1 - \beta)(db)^2 sdb=βsdb+(1β)(db)2

w : = w − α d w s d w + ϵ w := w - \alpha \frac{dw}{\sqrt{s_{dw} + \epsilon}} w:=wαsdw+ϵ dw

b : = b − α d b s d b + ϵ b := b - \alpha \frac{db}{\sqrt{s_{db} + \epsilon}} b:=bαsdb+ϵ db

其中, ε \varepsilon ε是一个实际操作时加上的较小数(例如 1 0 − 8 10^{-8} 108),为了防止分母太小而导致的数值不稳定。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第14张图片

如图所示,蓝色轨迹代表初始的移动,可以看到在 b b b方向上走得比较陡峭(即 d b db db较大),相比起来 d w dw dw较小,这影响了优化速度。

因此,在采用RMSProp算法后,由于 ( d w ) 2 (dw)^2 (dw)2较小、 ( d b ) 2 (db)^2 (db)2较大,进而 s d w s_{dw} sdw也会较小、 s d b s_{db} sdb也会较大,最终使得 d w s d w + ε \frac{dw}{\sqrt{s_{dw} + \varepsilon}} sdw+ε dw较大,而 d b s d b + ε \frac{db}{\sqrt{s_{db} + \varepsilon}} sdb+ε db较小。后面的更新就会像绿色轨迹一样,明显好于蓝色的更新曲线。RMSProp减小某些维度梯度更新波动较大的情况,使下降速度变得更快。

RMSProp有助于减少抵达最小值路径上的摆动,并允许使用一个更大的学习率 α \alpha α,从而加快算法学习速度。并且,它和Adam优化算法已被证明适用于不同的深度学习网络结构。

注意, β \beta β也是一个超参数。

对比原始梯度下降与RMSProp算法优化过程,如下图所示(上方为原始梯度下降,下方为RMSProp)

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第15张图片

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第16张图片

6.Adam 优化算法

Adam优化算法 Adam Optimization Algorithm

6.1 Adam算法介绍

**Adam (Adaptive Moment Estimation,自适应矩估计)**算法结合了动量梯度下降算法和RMSprop算法,通常有超越二者单独时的效果。具体过程如下(省略了 l l l):

首先进行初始化:

v d W = 0 , s d W = 0 , v d b = 0 , s d b = 0 v_{dW} = 0, s_{dW} = 0, v_{db} = 0, s_{db} = 0 vdW=0,sdW=0,vdb=0,sdb=0

用每一个Mini-Batch计算 d W dW dW d b db db,第 t t t次迭代时:

v d W = β 1 v d W + ( 1 − β 1 ) d W v_{dW} = \beta_1 v_{dW} + (1 - \beta_1) dW vdW=β1vdW+(1β1)dW

v d b = β 1 v d b + ( 1 − β 1 ) d b v_{db} = \beta_1 v_{db} + (1 - \beta_1) db vdb=β1vdb+(1β1)db

s d W = β 2 s d W + ( 1 − β 2 ) ( d W ) 2 s_{dW} = \beta_2 s_{dW} + (1 - \beta_2) {(dW)}^2 sdW=β2sdW+(1β2)(dW)2

s d b = β 2 s d b + ( 1 − β 2 ) ( d b ) 2 s_{db} = \beta_2 s_{db} + (1 - \beta_2) {(db)}^2 sdb=β2sdb+(1β2)(db)2

一般使用Adam算法时需要计算偏差修正:

v d W c o r r e c t e d = v d W 1 − β 1 t v^{corrected}_{dW} = \frac{v_{dW}}{1-{\beta_1}^t} vdWcorrected=1β1tvdW

v d b c o r r e c t e d = v d b 1 − β 1 t v^{corrected}_{db} = \frac{v_{db}}{1-{\beta_1}^t} vdbcorrected=1β1tvdb

s d W c o r r e c t e d = s d W 1 − β 2 t s^{corrected}_{dW} = \frac{s_{dW}}{1-{\beta_2}^t} sdWcorrected=1β2tsdW

s d b c o r r e c t e d = s d b 1 − β 2 t s^{corrected}_{db} = \frac{s_{db}}{1-{\beta_2}^t} sdbcorrected=1β2tsdb

所以,更新 W W W b b b时有:

W : = W − α v d W c o r r e c t e d s d W c o r r e c t e d + ε W := W - \alpha \frac{v^{corrected}_{dW}}{{\sqrt{s^{corrected}_{dW}} + \varepsilon}} W:=WαsdWcorrected +εvdWcorrected

b : = b − α v d b c o r r e c t e d s d b c o r r e c t e d + ε b := b - \alpha \frac{v^{corrected}_{db}}{{\sqrt{s^{corrected}_{db}} + \varepsilon}} b:=bαsdbcorrected +εvdbcorrected

6.2 Adam超参数的选择

Adam优化算法有很多的超参数,其中

  • 学习率 α \alpha α:需要尝试一系列的值,来寻找比较合适的
  • β 1 \beta_1 β1:常用的缺省值为0.9
  • β 2 \beta_2 β2:Adam算法的作者建议为0.999
  • ε \varepsilon ε:不重要,不会影响算法表现,Adam算法的作者建议为 1 0 − 8 10^{-8} 108

β 1 \beta_1 β1 β 2 \beta_2 β2 ε \varepsilon ε通常不需要调试。

对比原始梯度下降与RMSProp算法优化过程,如下图所示(上方为原始梯度下降,下方为Adam)

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第17张图片

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第18张图片

7.学习率衰减

学习率衰减 Learning Rate Decay

减小学习率 α \alpha α也能有效提高神经网络训练速度,这种方法被称为学习率衰减法(learning rate decay)。

学习率衰减就是随着迭代次数增加,学习率 α \alpha α逐渐减小。如下图示例。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第19张图片

蓝色折线表示设置一个固定的学习率 α \alpha α

  • 在最小值点附近,由于不同的Batch中存在一定的噪声,因此不会精确收敛,而是始终在最小值周围一个较大的范围内波动。

绿色折线表示随着时间慢慢减少学习率 α \alpha α的大小

  • 在初期 α \alpha α较大时,下降的步长较大,能以较快的速度进行梯度下降;
  • 后期逐步减小 α \alpha α的值,即减小步长,有助于算法的收敛,更容易接近最优解。

最常用的学习率衰减方法:

α = 1 1 + d e c a y _ r a t e ∗ e p o c h _ n u m ∗ α 0 \alpha = \frac{1}{1 + decay\_rate \ast epoch\_num} \ast \alpha_0 α=1+decay_rateepoch_num1α0

其中,decay_rate为衰减率(超参数),epoch_num为将所有的训练样本完整过一遍的次数。

  • 指数衰减: α = 0.9 5 e p o c h _ n u m ∗ α 0 \alpha = 0.95^{epoch\_num} \ast \alpha_0 α=0.95epoch_numα0
  • 其他: α = k e p o c h _ n u m ∗ α 0 \alpha = \frac{k}{\sqrt{epoch\_num}} \ast \alpha_0 α=epoch_num kα0
  • 离散下降

对于较小的模型,也有人会在训练时根据进度手动调小学习率。

8.局部最优问题

局部最优问题 the Problem of Local Optima

在使用梯度下降算法不断减小cost function时,可能会得到局部最优解(local optima)而不是全局最优解(global optima)

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第20张图片

之前我们对局部最优解的理解是形如碗状的凹槽,如图左边所示。但是在神经网络中,local optima的概念发生了变化。准确来说,大部分梯度为零的“最优点”并不是这些凹槽处,而是形如右边所示的马鞍状,称为saddle point。

所以在深度学习损失函数中,梯度为零并不能保证都是convex(极小值),也有可能是concave(极大值)。特别是在神经网络中参数很多的情况下,所有参数梯度为零的点很可能都是右边所示的马鞍状的saddle point,而不是左边那样的local optimum。

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第21张图片

类似马鞍状的plateaus会降低神经网络学习速度。Plateaus是梯度接近于零的平缓区域,如图所示。在plateaus上梯度很小,前进缓慢,到达saddle point需要很长时间。到达saddle point后,由于随机扰动,梯度一般能够沿着图中绿色箭头,离开saddle point,继续前进,只是在plateaus上花费了太多时间。

结论

  • 在训练较大的神经网络、存在大量参数,并且成本函数被定义在较高的维度空间时,困在极差的局部最优中是不大可能的;
  • 鞍点附近的平稳段会使得学习非常缓慢,而这也是动量梯度下降法、RMSProp以及Adam优化算法能够加速学习的原因,它们能帮助尽早走出平稳段。

参考资料

  • 图解数据分析
  • numpy速查手册

ShowMeAI系列教程推荐

  • 大厂技术实现:推荐与广告计算解决方案
  • 大厂技术实现:计算机视觉解决方案
  • 大厂技术实现:自然语言处理行业解决方案
  • 图解Python编程:从入门到精通系列教程
  • 图解数据分析:从入门到精通系列教程
  • 图解AI数学基础:从入门到精通系列教程
  • 图解大数据技术:从入门到精通系列教程
  • 图解机器学习算法:从入门到精通系列教程
  • 机器学习实战:手把手教你玩转机器学习系列
  • 深度学习教程:吴恩达专项课程 · 全套笔记解读
  • 自然语言处理教程:斯坦福CS224n课程 · 课程带学与全套笔记解读
  • 深度学习与计算机视觉教程:斯坦福CS231n · 全套笔记解读

推荐文章

  • ShowMeAI 深度学习教程(1) | 深度学习概论
  • ShowMeAI 深度学习教程(2) | 神经网络基础
  • ShowMeAI 深度学习教程(3) | 浅层神经网络
  • ShowMeAI 深度学习教程(4) | 深层神经网络
  • ShowMeAI 深度学习教程(5) | 深度学习的实用层面
  • ShowMeAI 深度学习教程(6) | 神经网络优化算法
  • ShowMeAI 深度学习教程(7) | 网络优化:超参数调优、正则化、批归一化和程序框架
  • ShowMeAI 深度学习教程(8) | AI应用实践策略(上)
  • ShowMeAI 深度学习教程(9) | AI应用实践策略(下)
  • ShowMeAI 深度学习教程(10) | 卷积神经网络解读
  • ShowMeAI 深度学习教程(11) | 经典CNN网络实例详解
  • ShowMeAI 深度学习教程(12) | CNN应用:目标检测
  • ShowMeAI 深度学习教程(13) | CNN应用:人脸识别和神经风格转换
  • ShowMeAI 深度学习教程(14) | 序列模型与RNN网络
  • ShowMeAI 深度学习教程(15) | 自然语言处理与词嵌入
  • ShowMeAI 深度学习教程(16) | Seq2seq序列模型和注意力机制

深度学习教程(6) | 神经网络优化算法(吴恩达·完整版)_第22张图片

你可能感兴趣的:(#,深度学习教程,◉,吴恩达专项课程最全笔记,深度学习,神经网络,优化算法,梯度下降,吴恩达)