三种梯度下降算法的比较和几种优化算法

博客已迁至知乎，本文链接：https://zhuanlan.zhihu.com/p/70910873

前言

这篇文章介绍了三种梯度下降方法的原理与优缺点，详细地讲解了Momentum、RMSprop和Adam优化算法，给出了使用建议。

---

三种梯度下降方法

1.Batch Gradient Descent，全部样本梯度下降一次，训练样本很大时，单次迭代需要时间太长。

2.Stochastic Gradient Descent，单个样本梯度下降一次，没有了向量化加速，效率比Batch Gradient Descent低，到达loss最低区域后还可能会跳出来，当然这也可以使它从局部最小值区域跳出来，可以使用学习率衰减来缓解这个问题。

3.Mini-batch Gradient Descent，部分样本梯度下降一次，上两个方法的折中，它可能不会收敛也可能不会在很小的范围内波动(同样可以用学习率衰减的方法来缓解这个问题)。

下面是loss的梯度图，三条线是三种梯度下降方法每下降一次的路线，蓝色是Batch Gradient Descent，紫色是Stochastic Gradient Descent，绿色是Mini-batch Gradient Descent。

进阶理解：
相对于GD，SGD能更有效的利用信息，特别是信息比较冗余的时候。举个例子，比如所有样本都需要向一个方向优化，GD优化一次需要对整个样本集迭代一次，而SGD只对一个样本优化就可以达到同样的效果。SGD相对于GD的另外一个优点是可以跳出局部最小值区域。
而mini-batch GD综合了两者的优点，既有了GD的向量化加速，还能像SGD更有效利用样本信息、可以跳出局部最小值区域的优点。另外，使用mini-batch，你还会发现不需要等待整个训练集被处理完就可以开始进行后续工作。

下面总结一下mini-batch的优点：
1.有向量化加速，加快了训练速度。
2.能有效利用样本信息，特别是信息比较冗余的时候。
3.有随机性，可以跳出局部最小值区域。
4.不需要等待整个训练集被处理完就可以开始进行后续工作。

下面是mini-batch的伪代码，中括号上标代表层数：
$$\begin{array}{rl}将样本分为n个minibatch& \\ fort=1,...n:& \\ & 前向传播:\\ & \{\begin{array}{l}{a}^{[1]}=g(W^{[1]}{X}_t+b^{[1]})\\ a^{[2]}=g(W^{[2]}{X}_t+b^{[2]})\\ ...\\ a^{[l]}=g(W^{[l]}{X}_t+b^{[l]})\end{array}\\ & 计算loss:L_{总}=\frac{1}{n}\sum _{i=1}^{l}L({\hat{y}}^{[i]},y^{[i]})\\ & 反向传播:\\ & \{\begin{array}{l}计算各层梯度dw和db\\ W^{[l]}=W^{[l]}-\alpha d{W}^{[l]}\\ b^{[l]}=b^{[l]}-\alpha d{b}^{[l]}\end{array}\end{array}$$

用法总结

首先，如果训练集较小，直接使用Batch Gradient Descent梯度下降法，样本集较小就没必要使用mini-batch梯度下降法，这里的少是说小于差不多2000个样本，这样比较适合使用Batch Gradient Descent梯度下降法。

样本数目较大的话，一般的mini-batch大小为64到512，考虑到电脑内存设置和使用的方式，如果mini-batch大小是2的次方，代码会运行地快一些。64到512的mini-batch比较常见。

---

下面讲几种常见的梯度下降优化算法：

动量梯度下降法(Momentum)

Gradient descent with Momentum，这个梯度下降方法，基本的想法就是计算梯度的指数加权平均数，并利用它更新权重。直观来讲，就是给普通的梯度下降加了个“惯性”，就像开车，你不能开着开着想往右拐就瞬间拐到右边，它有个向前再往右的过程，换言之，你想改变行驶方向，是需要从之前的行驶方向慢慢改变的，并不能瞬间改变。同理，Momentum梯度下降也一样，比如这次迭代算出来你需要向a方向优化，但你并不能直接将你的方向改成a，需要综合考虑之前的方向。
下图左边是普通随机梯度下降，右边是Momentum随机梯度下降，可以看出后者加快了优化速度，抑制了震荡。

因为mini-batch相比标准的梯度下降来说，更新参数更快，所以收敛过程会有浮动(loss下降曲线)，使用动量梯度下降法可以减小该浮动，还能加速训练。

看下mini-batch GD with Momentum的公式：

$$\begin{array}{rl}初始化每层的v_{dW}、v_{db}，& 形状和dW、db一致，元素全为0\\ 将样本分为n个minibatch& \\ fort=1,...n:& \\ & 前向传播：\\ & \{\begin{array}{l}{a}^{[1]}=g(W^{[1]}{X}_t+b^{[1]})\\ a^{[2]}=g(W^{[2]}{X}_t+b^{[2]})\\ ...\\ a^{[l]}=g(W^{[l]}{X}_t+b^{[l]})\end{array}\\ & 计算loss:L_{总}=\frac{1}{n}\sum _{i=1}^{l}L({\hat{y}}^{[i]},y^{[i]})\\ & 反向传播:\\ & \{\begin{array}{l}计算各层梯度dw和db\\ v_{d{W}^{[l]}}=\beta v_{d{W}^{[l-1]}}+(1-\beta )d{W}^{[l]}\\ \\ v_{d{b}^{[l]}}=\beta v_{d{b}^{[l-1]}}+(1-\beta )d{b}^{[l]}\\ \\ W^{[l]}=W^{[l]}-\alpha v_{d{W}^{[l]}}\\ \\ W^{[l]}=b^{[l]}-\alpha v_{d{b}^{[l]}}\end{array}\end{array}$$

$\beta$越大，收敛过程越平滑，一般取值为0.8~0.999，0.9会是一个不错的选择。

RMSprop

RMSprop的算法，全称是root mean square prop算法，它也可以加速收敛，我们来看看它是如何运作的。

$$\begin{array}{rl}初始化每层的s_{dW}、s_{db}，& 形状和dW、db一致，元素全为0\\ 将样本分为n个minibatch& \\ fort=1,...,n:& \\ & 前向传播：\\ & \{\begin{array}{l}{a}^{[1]}=g(W^{[1]}{X}_t+b^{[1]})\\ a^{[2]}=g(W^{[2]}{X}_t+b^{[2]})\\ ...\\ a^{[l]}=g(W^{[l]}{X}_t+b^{[l]})\end{array}\\ & 计算总loss:L_{总}=\frac{1}{n}\sum _{i=1}^{l}L({\hat{y}}^{[i]},y^{[i]})\\ & 反向传播:\\ & \{\begin{array}{l}计算各层梯度dw和db\\ s_{d{W}^{[l]}}=\beta s_{d{W}^{[l-1]}}+(1-\beta )(d{W}^{[l]}{)}^2\\ \\ s_{d{b}^{[l]}}=\beta s_{d{b}^{[l-1]}}+(1-\beta )(d{b}^{[l]}{)}^2\\ \\ W^{[l]}=W^{[l]}-\frac{\alpha }{\sqrt{s_{dW}+\epsilon }}d{W}^{[l]}\\ \\ b^{[l]}=b^{[l]}-\frac{\alpha }{\sqrt{s_{db}+\epsilon }}d{b}^{[l]}\end{array}\end{array}$$

$\epsilon$是一个很小的数，使后两个式子无论如何都不会除以一个接近于零的数，一般 $\epsilon =10^{-8}$

原理：参数更新时 $\frac{\alpha }{\sqrt{s_{dW}+\epsilon }}$ 的作用是使梯度大的参数不要更新地太猛。因为梯度$dW$越大，$s_{dW}$就越大，而$\frac{\alpha }{\sqrt{s_{dW}+\epsilon }}$就越小，参数更新的幅度就越小。至于为什么使用平方指数加权平均，可能是因为想要先平方再开方来取绝对值？有懂的大佬请不吝赐教。
总而言之，RMSprop使得梯度大的参数更新幅度不那么大，很大程度上缓解了梯度下降震荡的问题，如下图所示，蓝线代表普通梯度下降，绿线代表RMSprop。

为何不把Momentum和RMSprop结合在一起用呢？那就有了Adam。

Adam

Adam 优化算法(Adam optimization algorithm)，基本上就是将Momentum和RMSprop结合在一起，那么来看看如何使用Adam算法。

$$\begin{array}{rl}初始化每层的v_{dW}、v_{db}、s_{dW}、s_{db}，& 形状和dW、db一致，元素全为0\\ 将样本分为n个minibatch& \\ fort=1,...,n:& \\ & 前向传播：\\ & \{\begin{array}{l}{a}^{[1]}=g(W^{[1]}{X}_t+b^{[1]})\\ a^{[2]}=g(W^{[2]}{X}_t+b^{[2]})\\ ...\\ a^{[l]}=g(W^{[l]}{X}_t+b^{[l]})\end{array}\\ & 计算总loss:L_{总}=\frac{1}{n}\sum _{i=1}^{l}L({\hat{y}}^{[i]},y^{[i]})\\ & 反向传播:\\ & \{\begin{array}{l}计算各层梯度dw(db同理)\\ {\hat{v}}_{d{W}^{[l]}}={\beta }_1{\hat{v}}_{d{W}^{[l-1]}}+(1-{\beta }_1)d{W}^{[l]}\\ \\ {\hat{s}}_{d{W}^{[l]}}={\beta }_2{\hat{s}}_{d{W}^{[l-1]}}+(1-{\beta }_2)(d{W}^{[l]}{)}^2\\ \\ v_{d{W}^{[l]}}=\frac{{\hat{v}}_{d{W}^{[l]}}}{1-({\beta }_1{)}^t}\\ \\ s_{d{W}^{[l]}}=\frac{{\hat{s}}_{d{W}^{[l]}}}{1-({\beta }_2{)}^t}\\ \\ W^{[l]}=W^{[l]}-\alpha \frac{v_{d{W}^{[l]}}}{\sqrt{s_{d{W}^{[l]}}}+\epsilon }\end{array}\end{array}$$

• $t$ 表示Adam的步数
• $l$ 神经网络当前所在的层
• ${\beta }_1$ 一般为$0.9$；${\beta }_2$ 一般为$0.99$
• $\alpha$ 是 learning rate
• $\epsilon$ 是一个很小的数，防止除数为0，一般为$10^{-8}$

看一下反向传播中红色公式，前两个式子分别是Momentum和RMSprop。

至于第三和第四个式子，是为了使各个梯度的权值之和为1。以Momentum的式子为例，设定${\beta }_1=0.9$：

t=1和t=2时：
$$\begin{array}{rl}{\hat{v}}_{d{W}^{[1]}}& =0.1d{W}^{[1]}\\ v_{d{W}^{[1]}}& =\frac{0.1d{W}^{[1]}}{1-0.9^1}=d{W}^{[1]}\\ {\hat{v}}_{d{W}^{[2]}}& =0.09d{W}^{[1]}+0.1d{W}^{[2]}\\ v_{d{W}^{[2]}}& =\frac{0.09d{W}^{[1]}+0.1d{W}^{[2]}}{1-0.9^2}=0.47d{W}^{[1]}+0.53d{W}^{[2]}\end{array}$$

可以看到每一步中所有梯度的权值和都会变成1。

总结

梯度下降除非样本很少，一般来说都用mini-batch了，梯度下降的优化算法我推荐Adam，不需要花费太多时间来选择优化算法，效果都差不多，不如花时间去优化模型里别的超参。
下面贴了两张比较各个优化算法的动图，很直观。动图来自这篇博客：An overview of gradient descent optimization algorithms

---

References：

[1] Ruder (2017) An overview of gradient descent optimization algorithms
[2] Bottou et.al(2018) Optimization Methods for Large-Scale Machine Learning
[3] D. Lee et.al (2016) Gradient Descent Only Converges to Minimizers
[4] http://zh.gluon.ai/chapter_optimization/adam.html