deep_feedforward_networks.tex

% !Mode:: "TeX:UTF-8"
% Translator: Kai Li
\chapter{\glsentrytext{deep_feedforward_network}}
\label{chap:deep_feedforward_networks}

\firstgls{deep_feedforward_network}，也叫作\firstgls{feedforward_neural_network}或者\firstall{MLP}，是典型的深度学习模型。\gls{feedforward_network}的目标是近似某个函数$f^*$。
例如，对于分类器，$y = f^*(\Vx)$将输入$\Vx$映射到一个类别$y$。
\gls{feedforward_network}定义了一个映射$\Vy = f(\Vx; \Vtheta)$，并且学习参数$\Vtheta$的值，使它能够得到最佳的函数近似。

这种模型被称为\firstgls{feedforward}的，是因为信息流过$\Vx$的函数，流经用于定义$f$的中间计算过程，最终到达输出$\Vy$。
%这里的feedforward是翻译成前向的，还是前馈的？
在模型的输出和模型本身之间没有\firstgls{feedback}连接。
当前馈神经网络被扩展成包含\gls{feedback}连接时，它们被称为\firstgls{RNN}，在\chapref{chap:sequence_modeling_recurrent_and_recursive_nets}介绍。

\gls{feedforward_network}对于机器学习的从业者是极其重要的。
它们是许多重要商业应用的基础。
例如，用于对照片中的对象进行识别的\gls{CNN}就是一种专门的\gls{feedforward_network}。
\gls{feedforward_network}是通往\gls{recurrent_network}之路的概念基石，后者在自然语言的许多应用中发挥着巨大作用。

\gls{feedforward_neural_network}被称作\firstgls{network}是因为它们通常用许多不同函数复合在一起来表示。
该模型与一个有向无环图相关联，而图描述了函数是如何复合在一起的。
例如，我们有三个函数$f^{(1)}, f^{(2)}$和$f^{(3)}$连接在一个链上以形成$f(\Vx) = f^{(3)}(f^{(2)}(f^{(1)}(\Vx)) )$。
这些链式结构是神经网络中最常用的结构。
在这种情况下，$f^{(1)}$被称为网络的\firstgls{first_layer}，$f^{(2)}$被称为\firstgls{second_layer}，以此类推。
链的全长称为模型的\firstgls{depth}。
正是因为这个术语才出现了``深度学习''这个名字。
\gls{feedforward_network}的最后一层被称为\firstgls{output_layer}。
在神经网络训练的过程中，我们让$f(\Vx)$去匹配$f^*(\Vx)$的值。
训练数据为我们提供了在不同训练点上取值的、含有噪声的$f^*(\Vx)$的近似实例。
每个样本$\Vx$都伴随着一个标签$y\approx f^*(\Vx)$。
训练样本直接指明了输出层在每一点$\Vx$上必须做什么；它必须产生一个接近$y$的值。
但是训练数据并没有直接指明其他层应该怎么做。
学习算法必须决定如何使用这些层来产生想要的输出，但是训练数据并没有说每个单独的层应该做什么。
相反，学习算法必须决定如何使用这些层来最好地实现$f^*$的近似。
因为训练数据并没有给出这些层中的每一层所需的输出，所以这些层被称为\firstgls{hidden_layer}。

% -- 163 --

最后，这些网络被称为\emph{神经网络}是因为它们或多或少地受到神经科学的启发。
网络中的每个\gls{hidden_layer}通常都是向量值的。这些\gls{hidden_layer}的维数决定了模型的\firstgls{width}。
向量的每个元素都可以被视为起到类似一个神经元的作用。
除了将层想象成向量到向量的单个函数，我们也可以把层想象成由许多并行操作的\firstgls{unit}组成，每个单元表示一个向量到标量的函数。
每个单元在某种意义上类似一个神经元，它接收的输入来源于许多其他的单元，并计算它自己的激活值。
使用多层向量值表示的想法来源于神经科学。
用于计算这些表示的函数$f^{(i)}(\Vx)$的选择，也或多或少地受到神经科学观测的指引，这些观测是关于生物神经元计算功能的。
然而，现代的神经网络研究受到更多的是来自许多数学和工程学科的指引，并且神经网络的目标并不是完美地给大脑建模。
我们最好将\gls{feedforward_neural_network}想成是为了实现统计\gls{generalization}而设计出的函数近似机，它偶尔从我们了解的大脑中提取灵感，但并不是大脑功能的模型。

一种理解\gls{feedforward_network}的方式是从线性模型开始，并考虑如何克服它的局限性。
线性模型，例如\gls{logistic_regression}和线性回归，是非常吸引人的，因为无论是通过闭解形式还是使用凸优化，它们都能高效且可靠地拟合。
线性模型也有明显的缺陷，那就是该模型的能力被局限在线性函数里，所以它无法理解任何两个输入变量间的相互作用。

为了扩展线性模型来表示$\Vx$的非线性函数，我们可以不把线性模型用于$\Vx$本身，而是用在一个变换后的输入$\phi(\Vx)$上，这里$\phi$是一个非线性变换。
同样，我们可以使用\secref{sec:support_vector_machines}中描述的\gls{kernel_trick}，来得到一个基于隐含地使用$\phi$映射的非线性学习算法。
我们可以认为$\phi$提供了一组描述$\Vx$的特征，或者认为它提供了$\Vx$的一个新的表示。

% -- 164 --

剩下的问题就是如何选择映射$\phi$。
\begin{enumerate}
\item 其中一种选择是使用一个通用的$\phi$，例如无限维的$\phi$，它隐含地用在基于~\glssymbol{RBF}~核的\gls{kernel_machines}上。
如果$\phi(\Vx)$具有足够高的维数，我们总是有足够的能力来拟合训练集，但是对于测试集的\gls{generalization}往往不佳。
非常通用的特征映射通常只基于局部光滑的原则，并且没有将足够的先验信息进行编码来解决高级问题。

\item 另一种选择是手动地设计$\phi$。
在深度学习出现以前，这一直是主流的方法。
这种方法对于每个单独的任务都需要人们数十年的努力，从业者各自擅长特定的领域（如语音识别或计算机视觉），并且不同领域之间很难迁移(transfer)。

\item 深度学习的策略是去学习$\phi$。
在这种方法中，我们有一个模型$y = f(\Vx;\theta, \Vw) = \phi(\Vx; \theta)^\top \Vw$。
我们现在有两种参数：用于从一大类函数中学习$\phi$的参数$\Vtheta$，以及用于将$\phi(\Vx)$映射到所需的输出的参数$\Vw$。
这是\gls{deep_feedforward_network}的一个例子，其中$\phi$定义了一个\gls{hidden_layer}。
这是三种方法中唯一一种放弃了训练问题的凸性的，但是利大于弊。
在这种方法中，我们将表示参数化为$\phi(\Vx; \Vtheta)$，并且使用优化算法来寻找$\Vtheta$，使它能够得到一个好的表示。
如果我们想要的话，这种方法也可以通过使它变得高度通用以获得第一种方法的优点——我们只需使用一个非常广泛的函数族$\phi(\Vx; \Vtheta)$。
这种方法也可以获得第二种方法的优点。
人类专家可以将他们的知识编码进网络来帮助\gls{generalization}，他们只需要设计那些他们期望能够表现优异的函数族$\phi(\Vx; \Vtheta)$即可。
这种方法的优点是人类设计者只需要寻找正确的函数族即可，而不需要去寻找精确的函数。
\end{enumerate}

这种通过学习特征来改善模型的一般化原则不仅仅适用于本章描述的\gls{feedforward_neural_network}。
它是深度学习中反复出现的主题，适用于全书描述的所有种类的模型。
\gls{feedforward_neural_network}是这个原则的应用，它学习从$\Vx$到$\Vy$的确定性映射并且没有反馈连接。
后面出现的其他模型会把这些原则应用到学习随机映射、学习带有反馈的函数以及学习单个向量的概率分布。

% -- 165 --

本章我们先从\gls{feedforward_network}的一个简单例子说起。
接着，我们讨论部署一个\gls{feedforward_network}所需的每个设计决策。
首先，训练一个\gls{feedforward_network}至少需要做和线性模型同样多的设计决策：选择一个优化模型、代价函数以及输出单元的形式。
我们先回顾这些基于梯度学习的基本知识，然后去面对那些只出现在\gls{feedforward_network}中的设计决策。
\gls{feedforward_network}已经引入了\gls{hidden_layer}的概念，这需要我们去选择用于计算\gls{hidden_layer}值的\firstgls{activation_function}。
我们还必须设计网络的结构，包括网络应该包含多少层、这些层应该如何连接，以及每一层包含多少单元。
在深度神经网络的学习中需要计算复杂函数的梯度。
我们给出\firstgls{back_propagation}算法和它的现代推广，它们可以用来高效地计算这些梯度。
最后，我们以某些历史观点来结束这一章。

\section{实例：学习XOR}
\label{sec:example_learning_xor}

为了使\gls{feedforward_network}的想法更加具体，我们首先从一个可以完整工作的\gls{feedforward_network}说起。
这个例子解决一个非常简单的任务：学习XOR函数。

XOR函数（``异或''逻辑）是两个二进制值$x_1$和$x_2$的运算。
当这些二进制值中恰好有一个为1时，XOR函数返回值为1。
其余情况下返回值为0。
XOR函数提供了我们想要学习的目标函数$y = f^*(\Vx)$。
我们的模型给出了一个函数$y=f(\Vx; \Vtheta)$并且我们的学习算法会不断调整参数$\Vtheta$来使得$f$尽可能接近$f^*$。

在这个简单的例子中，我们不会关心统计\gls{generalization}。
我们希望网络在这四个点$\SetX =\{[0, 0]^\top, [0, 1]^\top, [1, 0]^\top, [1, 1]^\top\}$上表现正确。
我们会用全部这四个点来训练我们的网络，唯一的挑战是拟合训练集。

我们可以把这个问题当作是回归问题，并使用\gls{mean_squared_error}损失函数。
我们选择这个损失函数是为了尽可能简化本例中用到的数学。
在应用领域，对于二进制数据建模时，\glssymbol{mean_squared_error}通常并不是一个合适的损失函数。
更加合适的方法将在\secref{sec:sigmoid_units_for_bernoulli_output_distributions}中讨论。

% -- 166 --

评估整个训练集上表现的~\glssymbol{mean_squared_error}~损失函数为
\begin{equation}
J(\Vtheta) = \frac{1}{4} \sum_{\Vx\in \SetX} (f^*(\Vx) - f(\Vx; \Vtheta))^2.
\end{equation}

我们现在必须要选择我们模型$f(\Vx; \Vtheta)$的形式。
假设我们选择一个线性模型，$\Vtheta$包含$\Vw$和$b$，那么我们的模型被定义成
\begin{equation}
f(\Vx; \Vw, b) = \Vx^\top \Vw + b.
\end{equation}
我们可以使用\gls{normal_equations}关于$\Vw$和$b$最小化$J(\Vtheta)$，来得到一个闭式解。

解\gls{normal_equations}以后，我们得到$\Vw = 0$以及$b = \frac{1}{2}$。
线性模型仅仅是在任意一点都输出0.5。
为什么会发生这种事？
\figref{fig:chap6_xor_space_gray}演示了线性模型为什么不能用来表示XOR函数。
解决这个问题的其中一种方法是使用一个模型来学习一个不同的特征空间，在这个空间上线性模型能够表示这个解。
% fig 6.1
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/xor_space_gray}}
\fi
\captionsetup{singlelinecheck=off}
\caption{通过学习一个表示来解决XOR问题。
图上的粗体数字标明了学得的函数必须在每个点输出的值。
\emph{(左)}直接应用于原始输入的线性模型不能实现XOR函数。
当$x_1 = 0$时，模型的输出必须随着$x_2$的增大而增大。
当$x_1 = 1$时，模型的输出必须随着$x_2$的增大而减小。
线性模型必须对$x_2$使用固定的系数$w_2$。
因此，线性模型不能使用$x_1$的值来改变$x_2$的系数，从而不能解决这个问题。 
\emph{(右)}在由神经网络提取的特征表示的变换空间中，线性模型现在可以解决这个问题了。
在我们的示例解决方案中，输出必须为1的两个点折叠到了特征空间中的单个点。
换句话说，非线性特征将$\Vx = [1, 0]^\top$和$\Vx = [0, 1]^\top$都映射到了特征空间中的单个点$\Vh = [1, 0]^\top$。
线性模型现在可以将函数描述为$h_1$增大和$h_2$减小。
在该示例中，学习特征空间的动机仅仅是使得模型的能力更大，使得它可以拟合训练集。
在更现实的应用中，学习的表示也可以帮助模型\gls{generalization}。}
\label{fig:chap6_xor_space_gray}
\end{figure}


具体来说，我们这里引入一个非常简单的\gls{feedforward_neural_network}，它有一层\gls{hidden_layer}并且\gls{hidden_layer}中包含两个单元。
见\figref{fig:chap6_example_mlp}中对该模型的解释。
这个\gls{feedforward_network}有一个通过函数$f^{(1)}(\Vx;\MW, \Vc)$计算得到的\gls{hidden_unit}的向量$\Vh$。
这些\gls{hidden_unit}的值随后被用作第二层的输入。
第二层就是这个网络的输出层。
输出层仍然只是一个线性回归模型，只不过现在它作用于$\Vh$而不是$\Vx$。
网络现在包含链接在一起的两个函数：$\Vh=f^{(1)}(\Vx; \MW, \Vc)$和$y = f^{(2)}(\Vh; \Vw, b)$，完整的模型是$f(\Vx; \MW, \Vc, \Vw, b) = f^{(2)}(f^{(1)}(\Vx))$。
% fig 6.2
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/example_mlp}}
\fi
\captionsetup{singlelinecheck=off}
\caption{使用两种不同样式绘制的\gls{feedforward_network}的示例。
具体来说，这是我们用来解决XOR问题的\gls{feedforward_network}。
它有单个\gls{hidden_layer}，包含两个单元。
\emph{(左)}在这种样式中，我们将每个单元绘制为图中的一个节点。
这种风格是清楚而明确的，但对于比这个例子更大的网络，它可能会消耗太多的空间。
\emph{(右)}在这种样式中，我们将表示每一层激活的整个向量绘制为图中的一个节点。
这种样式更加紧凑。
有时，我们对图中的边使用参数名进行注释，这些参数是用来描述两层之间的关系的。
这里，我们用矩阵$\MW$描述从$\Vx$到$\Vh$的映射，用向量$\Vw$描述从$\Vh$到$y$的映射。
当标记这种图时，我们通常省略与每个层相关联的截距参数。}
\label{fig:chap6_example_mlp}
\end{figure}


$f^{(1)}$应该是哪种函数？线性模型到目前为止都表现不错，让$f^{(1)}$也是线性的似乎很有诱惑力。
可惜的是，如果$f^{(1)}$是线性的，那么\gls{feedforward_network}作为一个整体对于输入仍然是线性的。
暂时忽略截距项，假设$f^{(1)}(\Vx)= \MW^\top \Vx$并且$f^{(2)}(\Vh)=\Vh^\top \Vw$，那么$f(\Vx) = \Vw^\top\MW^\top \Vx$。
我们可以将这个函数重新表示成$f(\Vx) = \Vx^\top\Vw'$其中$\Vw' = \MW\Vw$。

% -- 167 --

显然，我们必须用非线性函数来描述这些特征。
大多数神经网络通过仿射变换之后紧跟着一个被称为激活函数的固定非线性函数来实现这个目标，其中仿射变换由学得的参数控制。
我们这里使用这种策略，定义$\Vh=g(\MW^\top \Vx+\Vc)$，其中$\MW$是线性变换的权重矩阵，$\Vc$是偏置。
此前，为了描述线性回归模型，我们使用权重向量和一个标量的偏置参数来描述从输入向量到输出标量的仿射变换。
现在，因为我们描述的是向量$\Vx$到向量$\Vh$的仿射变换，所以我们需要一整个向量的偏置参数。
激活函数$g$通常选择对每个元素分别起作用的函数，有$h_i =g(\Vx^\top \MW_{:, i} + c_i)$。
在现代神经网络中，默认的推荐是使用由激活函数$g(z)=\max\{0, z\}$定义的\firstgls{ReLU}或者称为~\glssymbol{ReLU}~\citep{Jarrett-ICCV2009-small,Nair-2010-small,Glorot+al-AI-2011-small}，如\figref{fig:chap6_relu_color}所示。
% fig 6.3
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/relu_color}}
\fi
\caption{\gls{rectified_linear}激活函数。
该激活函数是被推荐用于大多数\gls{feedforward_neural_network}的默认激活函数。
将此函数用于线性变换的输出将产生非线性变换。 然而，函数仍然非常接近线性，在这种意义上它是具有两个线性部分的分段线性函数。
由于\gls{ReLU}几乎是线性的，因此它们保留了许多使得线性模型易于使用基于梯度的方法进行优化的属性。 它们还保留了许多使得线性模型能够\gls{generalization}良好的属性。 计算机科学的一个通用原则是，我们可以从最小的组件构建复杂的系统。
就像图灵机的内存只需要能够存储0或1的状态，我们可以从\gls{rectified_linear}函数构建一个\gls{universal_function_approximator}。}
\label{fig:chap6_relu_color}
\end{figure}


我们现在可以指明我们的整个网络是
\begin{equation}
f(\Vx; \MW, \Vc, \Vw, b) = \Vw^\top \max\{ 0, \MW^\top \Vx + \Vc \} +b.
\end{equation}

% -- 168 --

% -- 169 --

我们现在可以给出XOR问题的一个解。
令
\begin{align}
\MW &= \begin{bmatrix}
1 & 1\\
1 & 1
\end{bmatrix},\\
\Vc &= \begin{bmatrix}
0\\
-1
\end{bmatrix},\\
\Vw &= \begin{bmatrix}
1\\
-2
\end{bmatrix},
\end{align}
以及$b=0$。

我们现在可以了解这个模型如何处理一批输入。
令$\MX$表示设计矩阵，它包含二进制输入空间中全部的四个点，每个样本占一行，那么矩阵表示为：
\begin{equation} 
\MX = \begin{bmatrix}
0 & 0\\
0 & 1\\ 
1 & 0\\ 
1 & 1
\end{bmatrix}. 
\end{equation} 
神经网络的第一步是将输入矩阵乘以第一层的权重矩阵： 
\begin{equation}
\MX\MW = \begin{bmatrix} 
0 & 0\\ 
1 & 1\\ 
1 & 1\\ 
2 & 2 
\end{bmatrix}.
\end{equation} 
然后，我们加上偏置向量$\Vc$，得到 
\begin{equation} 
\begin{bmatrix} 
0 & -1\\
1 & 0\\ 
1 & 0\\ 
2 & 1 
\end{bmatrix}. 
\end{equation}
在这个空间中，所有的样本都处在一条斜率为1的直线上。
当我们沿着这条直线移动时，输出需要从0升到1，然后再降回0。
线性模型不能实现这样一种函数。
为了用$\Vh$对每个样本求值，我们使用\gls{rectified_linear_transformation}： 
\begin{equation}
\begin{bmatrix} 
0 & 0\\ 
1 & 0\\ 
1 & 0\\ 
2 & 1 
\end{bmatrix}. 
\end{equation}

% -- 170 --

这个变换改变了样本间的关系。它们不再处于同一条直线上了。
如\figref{fig:chap6_xor_space_gray}所示，它们现在处在一个可以用线性模型解决的空间上。

我们最后乘以一个权重向量$\Vw$:
\begin{equation}
\begin{bmatrix}
0\\
1\\
1\\
0
\end{bmatrix}.
\end{equation}
神经网络对这一批次中的每个样本都给出了正确的结果。

在这个例子中，我们简单地指定了解决方案，然后说明它得到的误差为零。
在实际情况中，可能会有数十亿的模型参数以及数十亿的训练样本，所以不能像我们这里做的那样进行简单地猜解。
与之相对的，基于梯度的优化算法可以找到一些参数使得产生的误差非常小。
我们这里给出的XOR问题的解处在损失函数的全局最小点，所以梯度下降算法可以收敛到这一点。
梯度下降算法还可以找到XOR问题一些其他的等价解。
梯度下降算法的收敛点取决于参数的初始值。
在实践中，梯度下降通常不会找到像我们这里给出的那种干净的、容易理解的、整数值的解。

\section{基于梯度的学习}
\label{sec:gradient_based_learning}

设计和训练神经网络与使用梯度下降训练其他任何机器学习模型并没有太大不同。
在\secref{sec:building_a_machine_learning_algorithm}中，我们描述了如何通过指定一个优化过程、代价函数和一个模型族来构建一个机器学习算法。

我们到目前为止看到的线性模型和神经网络的最大区别，在于神经网络的非线性导致大多数我们感兴趣的\gls{cost_function}都变得非凸。
这意味着神经网络的训练通常使用迭代的、基于梯度的优化，仅仅使得代价函数达到一个非常小的值；而不是像用于训练线性回归模型的线性方程求解器，或者用于训练\gls{logistic_regression}或~\glssymbol{SVM}~的凸优化算法那样保证全局收敛。
凸优化从任何一种初始参数出发都会收敛（理论上如此——在实践中也很鲁棒但可能会遇到数值问题）。
用于非凸损失函数的随机梯度下降没有这种收敛性保证，并且对参数的初始值很敏感。
对于\gls{feedforward_neural_network}，将所有的权重值初始化为小随机数是很重要的。
偏置可以初始化为零或者小的正值。
这种用于训练\gls{feedforward_neural_network}以及几乎所有深度模型的迭代的基于梯度的优化算法会在\chapref{chap:optimization_for_training_deep_models}详细介绍，参数初始化会在\secref{sec:parameter_initialization_strategies}中具体说明。
就目前而言，只需要懂得，训练算法几乎总是基于使用梯度来使得代价函数下降的各种方法即可。
一些特别的算法是对梯度下降思想的改进和提纯（在\secref{sec:gradient_based_optimization}中介绍）还有一些更特别的，大多数是对随机梯度下降算法的改进（在\secref{sec:stochastic_gradient_descent_chap5}中介绍）。

% -- 171 --

我们当然也可以用梯度下降来训练诸如线性回归和\gls{SVM}之类的模型，并且事实上当训练集相当大时这是很常用的。
从这点来看，训练神经网络和训练其他任何模型并没有太大区别。
计算梯度对于神经网络会略微复杂一些，但仍然可以很高效而精确地实现。
\secref{sec:back_propagation_and_other_differentiation_algorithms}将会介绍如何用\gls{BP}算法以及它的现代扩展算法来求得梯度。

和其他的机器学习模型一样，为了使用基于梯度的学习方法我们必须选择一个代价函数，并且我们必须选择如何表示模型的输出。
现在，我们重温这些设计上的考虑，并且特别强调神经网络的情景。

\subsection{代价函数}
\label{sec:cost_functions}

深度神经网络设计中的一个重要方面是代价函数的选择。
幸运的是，神经网络的代价函数或多或少是和其他的参数模型例如线性模型的代价函数相同的。

在大多数情况下， 我们的参数模型定义了一个分布$p(\Vy\mid\Vx;\Vtheta)$并且我们简单地使用最大似然原理。
这意味着我们使用训练数据和模型预测间的\gls{cross_entropy}作为代价函数。

有时，我们使用一个更简单的方法，不是预测$\Vy$的完整概率分布，而是仅仅预测在给定$\Vx$的条件下$\Vy$的某种统计量。
某些专门的损失函数允许我们来训练这些估计量的预测器。

用于训练神经网络的完整的代价函数，通常在我们这里描述的基本代价函数的基础上结合一个正则项。
我们已经在\secref{sec:regularization}中看到正则化应用到线性模型中的一些简单的例子。
用于线性模型的权重衰减方法也直接适用于深度神经网络，而且是最流行的正则化策略之一。
用于神经网络的更高级的正则化策略将在\chapref{chap:regularization_for_deep_learning}中讨论。

% -- 172 --

\subsubsection{使用最大似然学习条件分布}
\label{sec:learning_conditional_distributions_with_maximum_likelihood}

大多数现代的神经网络使用最大似然来训练。
这意味着代价函数就是负的对数似然，它与训练数据和模型分布间的\gls{cross_entropy}等价。
这个代价函数表示为
\begin{equation}
J(\Vtheta) = -\SetE_{\RVx, \RVy \sim \hat{p}_\text{data}} \log p_\text{model} (\Vy \mid \Vx).
\end{equation}

代价函数的具体形式随着模型而改变，取决于$\log p_\text{model}$的具体形式。
上述方程的展开形式通常会有一些项不依赖于模型的参数，我们可以舍去。
例如，正如我们在\secref{sec:the_task_t}中看到的，如果$p_\text{model}(\Vy\mid\Vx) = \CalN(\Vy;f(\Vx;\Vtheta), \MI)$，那么我们就重新得到了\gls{mean_squared_error}代价，
\begin{equation}
J(\theta) = \frac{1}{2} \SetE_{\RVx, \RVy \sim  \hat{p}_\text{data}} || \Vy - f(\Vx; \Vtheta) ||^2 + \text{const},
\end{equation}
至少系数$\frac{1}{2}$和常数项不依赖于$\Vtheta$。
舍弃的常数是基于\gls{gaussian_distribution}的方差，在这种情况下我们选择不把它参数化。
之前，我们看到了对输出分布的最大似然估计和对线性模型\gls{mean_squared_error}的最小化之间的等价性，但事实上，这种等价性并不要求$f(\Vx; \Vtheta)$用于预测\gls{gaussian_distribution}的均值。

使用最大似然来导出代价函数的方法的一个优势是，它减轻了为每个模型设计代价函数的负担。
明确一个模型$p(\Vy\mid\Vx)$则自动地确定了一个代价函数$\log p(\Vy\mid\Vx)$。

贯穿神经网络设计的一个反复出现的主题是代价函数的梯度必须足够的大和具有足够的预测性，来为学习算法提供一个好的指引。
饱和（变得非常平）的函数破坏了这一目标，因为它们把梯度变得非常小。
这在很多情况下都会发生，因为用于产生\gls{hidden_unit}或者输出单元的输出的激活函数会饱和。
负的对数似然帮助我们在很多模型中避免这个问题。
很多输出单元都会包含一个指数函数，这在它的变量取绝对值非常大的负值时会造成饱和。
负对数似然代价函数中的对数函数消除了某些输出单元中的指数效果。
我们将会在\secref{sec:output_units}中讨论代价函数和输出单元的选择间的相互作用。

% -- 173 --

用于实现最大似然估计的\gls{cross_entropy}代价函数有一个不同寻常的特性，那就是当它被应用于实践中经常遇到的模型时，它通常没有最小值。
对于离散型输出变量，大多数模型以一种特殊的形式来参数化，即它们不能表示概率零和一，但是可以无限接近。
\gls{logistic_regression}是其中一个例子。
对于实值的输出变量，如果模型可以控制输出分布的密度（例如，通过学习\gls{gaussian_output_distribution}的方差参数），那么它可能对正确的训练集输出赋予极其高的密度，这将导致\gls{cross_entropy}趋向负无穷。
\chapref{chap:regularization_for_deep_learning}中描述的正则化技术提供了一些不同的方法来修正学习问题，使得模型不会通过这种方式来获得无限制的收益。

\subsubsection{学习条件统计量}
\label{sec:learning_conditional_statistics}

有时我们并不是想学习一个完整的概率分布$p(\Vy\mid\Vx; \Vtheta)$，而仅仅是想学习在给定$\Vx$时$\Vy$的某个条件统计量。

例如，我们可能有一个预测器$f(\Vx; \Vtheta)$，我们想用它来预测$\Vy$的均值。
如果我们使用一个足够强大的神经网络，我们可以认为这个神经网络能够表示一大类函数中的任何一个函数$f$，这个类仅仅被一些特征所限制，例如连续性和有界，而不是具有特殊的参数形式。
从这个角度来看，我们可以把代价函数看作是一个\firstgls{functional}而不仅仅是一个函数。
泛函是函数到实数的映射。我们因此可以将学习看作是选择一个函数而不仅仅是选择一组参数。
我们可以设计代价泛函在我们想要的某些特殊函数处取得最小值。
例如，我们可以设计一个代价泛函，使它的最小值处于一个特殊的函数上，这个函数将$\Vx$映射到给定$\Vx$时$\Vy$的期望值。
对函数求解优化问题需要用到\firstgls{calculus_of_variations}这个数学工具，我们将在\secref{sec:calculus_of_variations}中讨论。
理解\gls{calculus_of_variations}对于理解本章的内容不是必要的。
目前，只需要知道\gls{calculus_of_variations}可以被用来导出下面的两个结果。

% -- 174 --

我们使用\gls{calculus_of_variations}导出的第一个结果是解优化问题
\begin{equation}
f^* = \underset{f}{\argmin}  \ \SetE_{\RVx, \RVy \sim  p_\text{data}} ||\Vy-f(\Vx)||^2
\end{equation}
得到
\begin{equation}
f^*(\Vx) = \SetE_{\RVy\sim p_\text{data}(\Vy|\Vx)} [\Vy],
\end{equation}
要求这个函数处在我们要优化的类里。
换句话说，如果我们能够用无穷多的、来源于真实的数据生成分布的样本进行训练，最小化\gls{mean_squared_error}代价函数将得到一个函数，它可以用来对每个$\Vx$的值预测出$\Vy$的均值。

不同的代价函数给出不同的统计量。
第二个使用\gls{calculus_of_variations}得到的结果是
\begin{equation}
f^* = \underset{f}{\argmin} \ \SetE_{\RVx, \RVy \sim  p_\text{data}} ||\Vy - f(\Vx)||_1
\end{equation}
将得到一个函数可以对每个$\Vx$预测$\Vy$取值的\emph{中位数}，只要这个函数在我们要优化的函数族里。
这个代价函数通常被称为\firstgls{mean_absolute_error}。

可惜的是，\gls{mean_squared_error}和\gls{mean_absolute_error}在使用基于梯度的优化方法时往往成效不佳。
一些饱和的输出单元当结合这些代价函数时会产生非常小的梯度。
这就是为什么\gls{cross_entropy}代价函数比\gls{mean_squared_error}或者\gls{mean_absolute_error}更受欢迎的原因之一了，即使是在没必要估计整个$p(\Vy\mid\Vx)$分布时。


\subsection{输出单元}
\label{sec:output_units}

代价函数的选择与输出单元的选择紧密相关。
大多数时候，我们简单地使用数据分布和模型分布间的\gls{cross_entropy}。
选择如何表示输出决定了\gls{cross_entropy}函数的形式。

任何可用作输出的神经网络单元，也可以被用作\gls{hidden_unit}。
这里，我们着重讨论将这些单元用作模型输出时的情况，不过原则上它们也可以在内部使用。
我们将在\secref{sec:hidden_units}中重温这些单元，并且给出当它们被用作\gls{hidden_unit}时一些额外的细节。

在本节中，我们假设\gls{feedforward_network}提供了一组定义为$\Vh=f(\Vx;\Vtheta)$的隐藏特征。
输出层的作用是随后对这些特征进行一些额外的变换来完成整个网络必须完成的任务。

% -- 175 --

\subsubsection{用于\glsentrytext{gaussian_output_distribution}的线性单元}
\label{sec:linear_units_for_gaussian_output_distributions}

一种简单的输出单元是基于仿射变换的输出单元，仿射变换不具有非线性。
这些单元往往被直接称为线性单元。

给定特征$\Vh$，线性输出单元层产生一个向量$\hat{\Vy} = \MW^\top \Vh+\Vb$。

线性输出层经常被用来产生条件\gls{gaussian_distribution}的均值：
\begin{equation}
p(\Vy\mid\Vx) = \CalN(\Vy; \hat{\Vy}, \MI ).
\end{equation}
最大化其对数似然此时等价于最小化\gls{mean_squared_error}。

最大似然框架也使得学习\gls{gaussian_distribution}的协方差矩阵更加容易，或更容易地使\gls{gaussian_distribution}的协方差矩阵作为输入的函数。
然而，对于所有输入，协方差矩阵都必须被限定成一个正定矩阵。
线性输出层很难满足这种限定，所以通常使用其他的输出单元来对协方差参数化。
对协方差建模的方法将在\secref{sec:other_output_types}中简要介绍。

因为线性模型不会饱和，所以它们易于采用基于梯度的优化算法，甚至可以使用其他多种优化算法。

\subsubsection{用于\glsentrytext{bernoulli_output_distribution}的\glsentrytext{sigmoid}单元}
\label{sec:sigmoid_units_for_bernoulli_output_distributions}

许多任务需要预测二值型变量$y$的值。
具有两个类的分类问题可以归结为这种形式。

此时最大似然的方法是定义$y$在$\Vx$条件下的~\gls{bernoulli_distribution}。

\gls{bernoulli_distribution}仅需单个参数来定义。
神经网络只需要预测$P(y =1\mid\Vx)$即可。
为了使这个数是有效的概率，它必须处在区间$[0, 1]$中。

为满足该约束条件需要一些细致的设计工作。
假设我们打算使用线性单元，并且通过阈值来限制它成为一个有效的概率：
\begin{equation}
P(y=1 \mid \Vx) = \max \left \{ 0, \min \{1, \Vw^\top \Vh+b \} \right \}.
\end{equation}
这的确定义了一个有效的条件概率分布，但我们无法使用梯度下降来高效地训练它。
当$\Vw^\top \Vh+b$处于单位区间外时，模型的输出对其参数的梯度都将为$\bm{0}$。
梯度为$\bm{0}$通常是有问题的，因为学习算法对于如何改善相应的参数不再具有指导意义。

% -- 176 --

相反，最好是使用一种新的方法来保证无论何时模型给出了错误的答案时，总能有一个较大的梯度。
这种方法是基于使用~\gls{sigmoid}~输出单元结合最大似然来实现的。

\gls{sigmoid}~输出单元定义为
\begin{equation}
\hat{y} = \sigma \left (\Vw^\top \Vh + b \right ),
\end{equation}
这里$\sigma$是\secref{sec:useful_properties_of_common_functions}中介绍的~\gls{logistic_sigmoid}~函数。

我们可以认为~\gls{sigmoid}~输出单元具有两个部分。
首先，它使用一个线性层来计算$z=\Vw^\top \Vh+b$。
接着，它使用~\gls{sigmoid}~激活函数将$z$转化成概率。

我们暂时忽略对于$\Vx$的依赖性，只讨论如何用$z$的值来定义$y$的概率分布。
\gls{sigmoid}~可以通过构造一个非归一化（和不为1）的概率分布$\tilde{P}(y)$来得到。
我们可以随后除以一个合适的常数来得到有效的概率分布。
如果我们假定非归一化的对数概率对$y$ 和$z$是线性的，可以对它取指数来得到非归一化的概率。
我们然后对它归一化，可以发现这服从~\gls{bernoulli_distribution}，该分布受$z$的~\gls{sigmoid}~变换控制：
\begin{align}
\log \tilde{P}(y) &= yz,\\
\tilde{P}(y) &= \exp(yz),\\
P(y) &= \frac{\exp(yz)}{\sum_{y' = 0}^1 \exp(y' z)},\\
P(y) &= \sigma((2y-1)z).
\end{align}
基于指数和归一化的概率分布在统计建模的文献中很常见。
用于定义这种二值型变量分布的变量$z$被称为\firstgls{logit}。

% -- 177 --

这种在对数空间里预测概率的方法可以很自然地使用最大似然学习。
因为用于最大似然的代价函数是$-\log P(y\mid\Vx)$，代价函数中的log抵消了~\gls{sigmoid}~中的exp。
如果没有这个效果，\gls{sigmoid}~的饱和性会阻止基于梯度的学习做出好的改进。
我们使用最大似然来学习一个由~\gls{sigmoid}~参数化的~\gls{bernoulli_distribution}，它的损失函数为
\begin{align}
J(\Vtheta) &= -\log P(y\mid\Vx)\\
&= -\log \sigma ((2y-1)z)\\
&= \zeta((1-2y)z).
\end{align}

这个推导使用了\secref{sec:useful_properties_of_common_functions}中的一些性质。
通过将损失函数写成~\gls{softplus_function}的形式，我们可以看到它仅仅在$(1-2y)z$取绝对值非常大的负值时才会饱和。
因此饱和只会出现在模型已经得到正确答案时——当$y=1$且$z$取非常大的正值时，或者$y=0$且$z$取非常小的负值时。
当$z$的符号错误时，\gls{softplus_function}的变量$(1-2y)z$可以简化为$|z|$。
当$|z|$变得很大并且$z$的符号错误时，\gls{softplus_function}渐近地趋向于它的变量$|z|$。
对$z$求导则渐近地趋向于$\text{sign}(z)$，所以，对于极限情况下极度不正确的$z$，\gls{softplus_function}完全不会收缩梯度。
这个性质很有用，因为它意味着基于梯度的学习可以很快地改正错误的$z$。

当我们使用其他的损失函数，例如\gls{mean_squared_error}之类的，损失函数会在$\sigma(z)$饱和时饱和。
\gls{sigmoid}~激活函数在$z$取非常小的负值时会饱和到0，当$z$取非常大的正值时会饱和到1。
这种情况一旦发生，梯度会变得非常小以至于不能用来学习，无论此时模型给出的是正确还是错误的答案。
因此，最大似然几乎总是训练~\gls{sigmoid}~输出单元的优选方法。

理论上，\gls{sigmoid}~的对数总是确定和有限的，因为~\gls{sigmoid}~的返回值总是被限制在开区间$(0, 1)$上，而不是使用整个闭区间$[0, 1]$的有效概率。
在软件实现时，为了避免数值问题，最好将负的对数似然写作$z$的函数，而不是$\hat{y}=\sigma(z)$的函数。
如果~\gls{sigmoid}~函数下溢到零，那么之后对$\hat{y}$取对数会得到负无穷。

\subsubsection{用于
\glsentrytext{multinoulli_output_distribution}的softmax单元}
\label{sec:softmax_units_for_multinoulli_output_distributions}

任何时候当我们想要表示一个具有$n$个可能取值的离散型随机变量的分布时，我们都可以使用softmax函数。
它可以看作是~\gls{sigmoid}~函数的扩展，其中~\gls{sigmoid}~函数用来表示二值型变量的分布。

% -- 178 --

softmax函数最常用作分类器的输出，来表示$n$个不同类上的概率分布。
比较少见的是，softmax函数可以在模型内部使用，例如如果我们想要在某个内部变量的$n$个不同选项中进行选择。

在二值型变量的情况下，我们希望计算一个单独的数
\begin{equation}
\hat{y} = P(y=1\mid\Vx).
\end{equation}
因为这个数需要处在0和1之间，并且我们想要让这个数的对数可以很好地用于对数似然的基于梯度的优化，我们选择去预测另外一个数$z=\log \hat{P}(y=1\mid\Vx)$。
对其指数化和归一化，我们就得到了一个由~\gls{sigmoid}~函数控制的~\gls{bernoulli_distribution}。

为了推广到具有$n$个值的离散型变量的情况，我们现在需要创造一个向量$\hat{\Vy}$，它的每个元素是$\hat{y}_i = P(y=i\mid\Vx)$。
我们不仅要求每个$\hat{y}_i$元素介于0和1之间，还要使得整个向量的和为1，使得它表示一个有效的概率分布。
用于~\gls{bernoulli_distribution}的方法同样可以推广到~\gls{multinoulli_distribution}。
首先，线性层预测了未归一化的对数概率：
\begin{equation}
\Vz = \MW^\top \Vh+\Vb,
\end{equation}
其中$z_i=\log \hat{P}(y=i\mid\Vx)$。
softmax函数然后可以对$z$指数化和归一化来获得需要的$\hat{\Vy}$。
最终，softmax函数的形式为
\begin{equation}
\text{softmax}(\Vz)_i = \frac{\exp(z_i)}{\sum_j \exp(z_j)}.
\end{equation}

和~\gls{logistic_sigmoid}一样，当使用最大化对数似然训练softmax来输出目标值$\RSy$时，使用指数函数工作地非常好。
这种情况下，我们想要最大化$\log P(\RSy =i; \Vz)=\log \text{softmax}(\Vz)_i$。
将softmax定义成指数的形式是很自然的因为对数似然中的log可以抵消softmax中的exp：
\begin{equation}
\log \text{softmax}(\Vz)_i = z_i - \log \sum_j \exp(z_j).
\label{eq:6.30}
\end{equation}

% -- 179 --

\eqnref{eq:6.30}中的第一项表示输入$z_i$总是对代价函数有直接的贡献。
因为这一项不会饱和，所以即使$z_i$对\eqnref{eq:6.30}的第二项的贡献很小，学习依然可以进行。
当最大化对数似然时，第一项鼓励$z_i$被推高，而第二项则鼓励所有的$\Vz$被压低。
为了对第二项$\log \sum_j \exp(z_j)$有一个直观的理解，注意到这一项可以大致近似为$\max_j z_j$。
这种近似是基于对任何明显小于$\max_j z_j$的$z_k$，$\exp(z_k)$都是不重要的。
我们能从这种近似中得到的直觉是，负对数似然代价函数总是强烈地惩罚最活跃的不正确预测。
如果正确答案已经具有了softmax的最大输入，那么$-z_i$项和$\log\sum_j \exp(z_j) \approx \max_j z_j = z_i$项将大致抵消。
这个样本对于整体训练代价贡献很小，这个代价主要由其他未被正确分类的样本产生。

到目前为止我们只讨论了一个例子。
总体来说，未正则化的最大似然会驱动模型去学习一些参数，而这些参数会驱动softmax函数来预测在训练集中观察到的每个结果的比率：
\begin{equation}
\text{softmax}(\Vz(\Vx; \Vtheta))_i \approx \frac{\sum_{j=1}^m \bm{1}_{y^{(j)}=i, \Vx^{(j)} = \Vx}  }{ \sum_{j=1}^{m} \bm{1}_{\Vx^{(j)} = \Vx} }.
\end{equation}
因为最大似然是一致的估计量，所以只要模型族能够表示训练的分布，这就能保证发生。
在实践中，有限的模型能力和不完美的优化将意味着模型只能近似这些比率。

除了对数似然之外的许多目标函数对softmax函数不起作用。
具体来说，那些不使用对数来抵消softmax中的指数的目标函数，当指数函数的变量取非常小的负值时会造成梯度消失，从而无法学习。
特别是，平方误差对于softmax单元来说是一个很差的损失函数，即使模型做出高度可信的不正确预测，也不能训练模型改变其输出\citep{Bridle90}。
要理解为什么这些损失函数可能失败，我们需要检查softmax函数本身。

像~\gls{sigmoid}~一样，softmax激活函数可能会饱和。
\gls{sigmoid}~函数具有单个输出，当它的输入极端负或者极端正时会饱和。
对于softmax的情况，它有多个输出值。
当输入值之间的差异变得极端时，这些输出值可能饱和。
当softmax饱和时，基于softmax的许多代价函数也饱和，除非它们能够转化饱和的激活函数。

% -- 180 --

为了说明softmax函数对于输入之间差异的响应，观察到当对所有的输入都加上一个相同常数时softmax的输出不变：
\begin{equation}
\text{softmax}(\Vz) = \text{softmax}(\Vz+c).
\end{equation}
使用这个性质，我们可以导出一个数值方法稳定的softmax函数的变体：
\begin{equation}
\text{softmax}(\Vz) = \text{softmax}(\Vz- \max_i z_i).
\end{equation}
变换后的形式允许我们在对softmax函数求值时只有很小的数值误差，即使是当$z$包含极正或者极负的数时。
观察softmax数值稳定的变体，可以看到softmax函数由它的变量偏离$\max_i z_i$的量来驱动。

当其中一个输入是最大（$z_i = \max_i z_i$）并且$z_i$远大于其他的输入时，相应的输出$\text{softmax}(\Vz)_i$会饱和到1。
当$z_i$不是最大值并且最大值非常大时，相应的输出$\text{softmax}(\Vz)_i$也会饱和到0。
这是~\gls{sigmoid}~单元饱和方式的一般化，并且如果损失函数不被设计成对其进行补偿，那么也会造成类似的学习困难。

softmax函数的变量$\Vz$可以通过两种方式产生。
最常见的是简单地使神经网络中softmax之前的层输出$\Vz$的每个元素，就像先前描述的使用线性层$\Vz={W}^\top\Vh+\Vb$。
虽然很直观，但这种方法是对分布的过度参数化。
$n$个输出总和必须为1的约束意味着只有$n-1$个参数是必要的；第$n$个概率值可以通过1减去前面$n-1$个概率来获得。
因此，我们可以强制要求$\Vz$的一个元素是固定的。
例如，我们可以要求$z_n=0$。
事实上，这正是~\gls{sigmoid}~单元所做的。
定义$P(y=1\mid\Vx)=\sigma(z)$等价于用二维的$\Vz$以及$z_1=0$来定义$P(y=1\mid\Vx)=\text{softmax}(\Vz)_1$。
无论是$n-1$个变量还是$n$个变量的方法，都描述了相同的概率分布，但会产生不同的学习机制。
在实践中，无论是过度参数化的版本还是限制的版本都很少有差别，并且实现过度参数化的版本更为简单。

从神经科学的角度看，有趣的是认为softmax是一种在参与其中的单元之间形成竞争的方式：softmax输出总是和为1，所以一个单元的值增加必然对应着其他单元值的减少。
这与被认为存在于皮质中相邻神经元间的侧抑制类似。
在极端情况下（当最大的$a_i$和其他的在幅度上差异很大时），它变成了\firstgls{winner_take_all}的形式（其中一个输出接近1，其他的接近0）。

``softmax''的名称可能会让人产生困惑。
这个函数更接近于argmax函数而不是max函数。
``soft''这个术语来源于softmax函数是连续可微的。
``argmax''函数的结果表示为一个\gls{one_hot}向量（只有一个元素为1，其余元素都为0的向量），不是连续和可微的。
softmax函数因此提供了argmax的``软化''版本。max函数相应的软化版本是$\text{softmax}(\Vz)^\top \Vz$。
可能最好是把softmax函数称为``softargmax''，但当前名称已经是一个根深蒂固的习惯了。

% -- 181 --

\subsubsection{其他的输出类型}
\label{sec:other_output_types}

之前描述的线性、\gls{sigmoid}~和softmax输出单元是最常见的。
神经网络可以推广到我们希望的几乎任何种类的输出层。
最大似然原则给如何为几乎任何种类的输出层设计一个好的代价函数提供了指导。

一般的，如果我们定义了一个条件分布$p(\Vy\mid\Vx; \Vtheta)$，最大似然原则建议我们使用$-\log p(\Vy\mid \Vx;\Vtheta)$作为代价函数。

一般来说，我们可以认为神经网络表示函数$f(\Vx;\Vtheta)$。
这个函数的输出不是对$\Vy$值的直接预测。
相反，$f(\Vx;\Vtheta)=\Vomega$提供了$y$分布的参数。
我们的损失函数就可以表示成$-\log p(\RVy; \Vomega(\Vx))$。

例如，我们想要学习在给定$\RVx$时，$\RVy$的条件\gls{gaussian_distribution}的方差。
简单情况下，方差$\sigma^2$是一个常数，此时有一个解析表达式，这是因为方差的最大似然估计量仅仅是观测值$\RVy$与它们的期望值的差值的平方平均。
一种计算上代价更加高但是不需要写特殊情况代码的方法是简单地将方差作为分布$p(\RVy\mid\Vx)$的其中一个属性，这个分布由$\Vomega=f(\Vx;\Vtheta)$控制。
负对数似然$-\log p(\Vy;\Vomega(\Vx))$将为代价函数提供一个必要的合适项来使我们的优化过程可以逐渐地学到方差。
在标准差不依赖于输入的简单情况下，我们可以在网络中创建一个直接复制到$\Vomega$中的新参数。
这个新参数可以是$\sigma$本身，或者可以是表示$\sigma^2$的参数$v$，或者可以是表示$\frac{1}{\sigma^2}$的参数$\beta$，取决于我们怎样对分布参数化。
我们可能希望模型对不同的$\RVx$值预测出$\RVy$不同的方差。
这被称为\firstgls{heteroscedastic}模型。
在异方差情况下，我们简单地把方差指定为$f(\RVx;\Vtheta)$其中一个输出值。
实现它的典型方法是使用精度而不是方差来表示\gls{gaussian_distribution}，就像\eqnref{eq:3.22}所描述的。
在多维变量的情况下，最常见的是使用一个对角精度矩阵
\begin{equation}
\text{diag}(\Vbeta).
\end{equation}
这个公式适用于梯度下降，因为由$\Vbeta$参数化的\gls{gaussian_distribution}的对数似然的公式仅涉及$\beta_i$的乘法和$\log \beta_i$的加法。
乘法、加法和对数运算的梯度表现良好。
相比之下，如果我们用方差来参数化输出，我们需要用到除法。
除法函数在零附近会变得任意陡峭。
虽然大梯度可以帮助学习，但任意大的梯度通常导致不稳定。
如果我们用标准差来参数化输出，对数似然仍然会涉及除法，并且还将涉及平方。
通过平方运算的梯度可能在零附近消失，这使得学习被平方的参数变得困难。
无论我们使用的是标准差，方差还是精度，我们必须确保\gls{gaussian_distribution}的协方差矩阵是正定的。
因为精度矩阵的特征值是协方差矩阵特征值的倒数，所以这等价于确保精度矩阵是正定的。
如果我们使用对角矩阵，或者是一个常数乘以单位矩阵\footnote{译者注：这里原文是``If we use a diagonal matrix, or a scalar times the diagonal matrix...''即``如果我们使用对角矩阵，或者是一个标量乘以对角矩阵...''，但一个标量乘以对角矩阵和对角矩阵没区别，结合上下文可以看出，这里原作者误把``identity''写成了``diagonal matrix''，因此这里采用``常数乘以单位矩阵''的译法。}，那么我们需要对模型输出强加的唯一条件是它的元素都为正。
如果我们假设$\Va$是用于确定对角精度的模型的原始激活，那么可以用softplus 函数来获得正的精度向量：$\Vbeta=\Vzeta(\Va)$。
这种相同的策略对于方差或标准差同样适用，也适用于常数乘以单位阵的情况。

% -- 182 --

学习一个比对角矩阵具有更丰富结构的协方差或者精度矩阵是很少见的。
如果协方差矩阵是满的和有条件的，那么参数化的选择就必须要保证预测的协方差矩阵是正定的。
这可以通过写成$\VSigma(\Vx)=\MB(\Vx)\MB^\top (\Vx)$来实现，这里$\MB$是一个无约束的方阵。
如果矩阵是满秩的，那么一个实际问题是计算似然的代价是很高的，计算一个$d\times d$的矩阵的行列式或者$\VSigma(\Vx)$的逆（或者等价地并且更常用地，对它特征值分解或者$\MB(\Vx)$的特征值分解）需要$O(d^3)$的计算量。

% -- 183 --

我们经常想要执行多峰回归(multimodal regression)，即预测条件分布$p(\Vy\mid\Vx)$的实值，该条件分布对于相同的$\Vx$值在$\Vy$空间中有多个不同的峰值。
在这种情况下，高斯混合是输出的自然表示\citep{Jacobs-nc91,bishop1994mixture}。
将高斯混合作为其输出的神经网络通常被称为\firstgls{mixture_density_network}。
具有$n$个分量的高斯混合输出由下面的条件分布定义：
\begin{equation}
p(\Vy\mid\Vx) = \sum_{i=1}^n p(\RSc = i \mid \Vx) \CalN(\Vy; \Vmu^{(i)}(\Vx), \VSigma^{(i)}(\Vx)).
\end{equation}
神经网络必须有三个输出：定义$p(\RSc=i\mid\Vx)$的向量，对所有的$i$给出$\Vmu^{(i)}(\Vx)$的矩阵，以及对所有的$i$给出$\VSigma^{(i)}(\Vx)$的张量。
这些输出必须满足不同的约束：
\begin{enumerate}
\item 混合组件$p(\RSc=i\mid\Vx)$：它们由\gls{latent_variable}\footnote{我们之所以认为$\RSc$是\gls{latent}的，是因为我们不能直接在数据中观测到它：给定输入$\RVx$和目标$\RVy$，不可能确切地知道是哪个高斯组件产生$\RVy$，但我们可以想象$\RVy$是通过选择其中一个来产生的，并且将那个未被观测到的选择作为随机变量。}~$\RSc$关联着，在$n$个不同组件上形成~\gls{multinoulli_distribution}。
这个分布通常可以由$n$维向量的softmax来获得，以确保这些输出是正的并且和为1。

\item 均值$\Vmu^{(i)}(\Vx)$：它们指明了与第$i$个高斯组件相关联的中心或者均值，并且是无约束的（通常对于这些输出单元完全没有非线性）。
如果$\RVy$是个$d$维向量，那么网络必须输出一个由$n$个这种$d$维向量组成的$n\times d$ 的矩阵。
用最大似然来学习这些均值要比学习只有一个输出模式的分布的均值稍稍复杂一些。我们只想更新那个真正产生观测数据的组件的均值。
在实践中，我们并不知道是哪个组件产生了观测数据。
负对数似然表达式将每个样本对每个组件的贡献进行赋权，权重的大小由相应的组件产生这个样本的概率来决定。

\item 协方差$\VSigma^{(i)}(\Vx)$：它们指明了每个组件$i$的协方差矩阵。
和学习单个高斯组件时一样，我们通常使用对角矩阵来避免计算行列式。
和学习混合均值时一样，最大似然是很复杂的，它需要将每个点的部分责任分配给每个混合组件。
如果给定了混合模型的正确的负对数似然，梯度下降将自动地遵循正确的过程。
\end{enumerate}
有报告说基于梯度的优化方法对于混合条件高斯（作为神经网络的输出）可能是不可靠的，部分是因为涉及到除法（除以方差）可能是数值不稳定的（当某个方差对于特定的实例变得非常小时，会导致非常大的梯度）。
一种解决方法是\firstgls{clip_gradients}（见\secref{sec:clipping_gradients}），另外一种是启发式缩放梯度\citep{Uria+al-ICML2014}。

% -- 184 --

高斯混合输出在语音生成模型\citep{schuster1999supervised}和物理运动\citep{Graves-arxiv2013}中特别有效。
混合密度策略为网络提供了一种方法来表示多种输出模式，并且控制输出的方差，这对于在这些实数域中获得高质量的结果是至关重要的。
混合密度网络的一个实例如\figref{fig:chap6_mdn_color}所示。
% fig 6.4
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/mdn_color}}
\fi
\captionsetup{singlelinecheck=off}
\caption{从具有混合密度输出层的神经网络中抽取的样本。
输入$x$从\gls{uniform_distribution}中采样，输出$y$从$p_{\text{model}}(y\mid x)$中采样。 神经网络能够学习从输入到输出分布的参数的非线性映射。 这些参数包括控制三个组件中的哪一个将产生输出的概率，以及每个组件各自的参数。
每个混合组件都是\gls{gaussian_distribution}，具有预测的均值和方差。 输出分布的这些方面都能够相对输入$x$变化，并且以非线性的方式改变。}
\label{fig:chap6_mdn_color}
\end{figure}


一般的，我们可能希望继续对包含更多变量的、更大的向量$\Vy$来建模，并在这些输出变量上施加更多更丰富的结构。
例如，我们可能希望神经网络输出字符序列形成一个句子。
在这些情况下，我们可以继续使用最大似然原理应用到我们的模型$p(\Vy;\Vomega(\Vx))$上，但我们用来描述$\Vy$的模型会变得非常复杂，超出了本章的范畴。
\chapref{chap:sequence_modeling_recurrent_and_recursive_nets}描述了如何使用循环神经网络来定义这种序列上的模型，第\ref{part:deep_learning_research}部分描述了对任意概率分布进行建模的高级技术。

\section{\glsentrytext{hidden_unit}}
\label{sec:hidden_units}

到目前为止，我们集中讨论了神经网络的设计选择，这对于使用基于梯度的优化方法来训练的大多数参数化机器学习模型都是通用的。
现在我们转向一个\gls{feedforward_neural_network}独有的问题：该如何选择\gls{hidden_unit}的类型，这些\gls{hidden_unit}用在模型的\gls{hidden_layer}中。

% -- 185 --

\gls{hidden_unit}的设计是一个非常活跃的研究领域，并且还没有许多明确的指导性理论原则。

\gls{ReLU}是\gls{hidden_unit}极好的默认选择。
许多其他类型的\gls{hidden_unit}也是可用的。
决定何时使用哪种类型的\gls{hidden_unit}是困难的事（尽管\gls{ReLU}通常是一个可接受的选择）。
我们这里描述对于每种\gls{hidden_unit}的一些基本直觉。
这些直觉可以用来建议我们何时来尝试一些单元。
通常不可能预先预测出哪种\gls{hidden_unit}工作得最好。
设计过程充满了试验和错误，先直觉认为某种\gls{hidden_unit}可能表现良好，然后用它组成神经网络进行训练，最后用验证集来评估它的性能。

这里列出的一些\gls{hidden_unit}可能并不是在所有的输入点上都是可微的。
例如，\gls{ReLU} $g(z)=\max\{0, z\}$在$z=0$处不可微。
这似乎使得$g$对于基于梯度的学习算法无效。
在实践中，梯度下降对这些机器学习模型仍然表现得足够好。
部分原因是神经网络训练算法通常不会达到代价函数的局部最小值，而是仅仅显著地减小它的值，如\figref{fig:chap4_approx_opt_color}所示。
这些想法会在\chapref{chap:optimization_for_training_deep_models}中进一步描述。
因为我们不再期望训练能够实际到达梯度为$\bm{0}$的点，所以代价函数的最小值对应于梯度未定义的点是可以接受的。
不可微的\gls{hidden_unit}通常只在少数点上不可微。
一般来说，函数$g(z)$具有左导数和右导数，左导数定义为紧邻在$z$左边的函数的斜率，右导数定义为紧邻在$z$右边的函数的斜率。
只有当函数在$z$处的左导数和右导数都有定义并且相等时，函数在$z$点处才是可微的。
神经网络中用到的函数通常对左导数和右导数都有定义。
在$g(z)=\max\{0,z\}$的情况下，在$z=0$处的左导数是0，右导数是1。
神经网络训练的软件实现通常返回左导数或右导数的其中一个，而不是报告导数未定义或产生一个错误。
这可以通过观察到在数字计算机上基于梯度的优化总是会受到数值误差的影响来启发式地给出理由。
当一个函数被要求计算$g(0)$时，底层值真正为0 是不太可能的。
相对的，它可能是被舍入为0的一个小量$\epsilon$。
在某些情况下，理论上有更好的理由，但这些通常对神经网络训练并不适用。
重要的是，在实践中，我们可以放心地忽略下面描述的\gls{hidden_unit}激活函数的不可微性。

% -- 186 --

除非另有说明，大多数的\gls{hidden_unit}都可以描述为接受输入向量$\Vx$，计算仿射变换$\Vz=\MW^\top \Vx+\Vb$，然后使用一个逐元素的非线性函数$g(\Vz)$。
大多数\gls{hidden_unit}的区别仅仅在于激活函数$g(\Vz)$的形式。

\subsection{\glsentrytext{ReLU}及其扩展}
\label{sec:rectified_linear_units_and_their_generalizations}

\gls{ReLU}使用激活函数$g(z)=\max\{0, z\}$。

\gls{ReLU}易于优化，因为它们和线性单元非常类似。
线性单元和\gls{ReLU}的唯一区别在于\gls{ReLU}在其一半的定义域上输出为零。
这使得只要\gls{ReLU}处于激活状态，它的导数都能保持较大。
它的梯度不仅大而且一致。
整流操作的二阶导数几乎处处为0，并且在\gls{ReLU}处于激活状态时，它的一阶导数处处为1。
这意味着相比于引入二阶效应的激活函数来说，它的梯度方向对于学习来说更加有用。

\gls{ReLU}通常作用于仿射变换之上：
\begin{equation}
\Vh = g(\MW^\top \Vx + \Vb).
\end{equation}
当初始化仿射变换的参数时，可以将$\Vb$的所有元素设置成一个小的正值，例如0.1。
这使得\gls{ReLU}很可能初始时就对训练集中的大多数输入呈现激活状态，并且允许导数通过。

有很多\gls{ReLU}的扩展存在。
大多数这些扩展的表现比得上\gls{ReLU}，并且偶尔表现得更好。

\gls{ReLU}的一个缺陷是它们不能通过基于梯度的方法学习那些使它们激活为零的样本。
\gls{ReLU}的各种扩展保证了它们能在各个位置都接收到梯度。

\gls{ReLU}的三个扩展基于当$z_i<0$时使用一个非零的斜率$\alpha_i$：$h_i =g(\Vz, \Valpha)_i = \max(0, z_i) + \alpha_i \min(0, z_i)$。
\firstgls{absolute_value_rectification}固定$\alpha_i=-1$来得到$g(z)=|z|$。
它用于图像中的对象识别\citep{Jarrett-ICCV2009}，其中寻找在输入照明极性反转下不变的特征是有意义的。
\gls{ReLU}的其他扩展比这应用地更广泛。
\firstgls{leaky_ReLU}\citep{Maas-et-al-ICML2013}将$\alpha_i$固定成一个类似0.01的小值，\firstgls{PReLU}或者~\textbf{\glssymbol{PReLU}}~将$\alpha_i$作为学习的参数\citep{He-et-al-arxiv2015}。

% -- 187 --

\firstgls{maxout_unit}\citep{Goodfellow_maxout_2013}进一步扩展了\gls{ReLU}。
\gls{maxout_unit}将$\Vz$划分为每组具有$k$个值的组，而不是使用作用于每个元素的函数$g(z)$。
每个\gls{maxout_unit}则输出每组中的最大元素：
\begin{equation}
g(\Vz)_i = \underset{j\in \SetG^{(i)}}{\max} z_j
\end{equation}
这里$\SetG^{(i)}$是组$i$的输入索引集$\{(i-1)k+1, \ldots, ik\}$。
这提供了一种方法来学习对输入$\Vx$空间中多个方向响应的分段线性函数。

\gls{maxout_unit}可以学习具有多达$k$段的分段线性的凸函数。
\gls{maxout_unit}因此可以视为\emph{学习激活函数}本身而不仅仅是单元之间的关系。
使用足够大的$k$，\gls{maxout_unit}可以以任意的精确度来近似任何凸函数。
特别地，具有两块的maxout层可以学习实现和传统层相同的输入$\Vx$的函数，这些传统层可以使用\gls{rectified_linear}激活函数、\gls{absolute_value_rectification}、\gls{leaky_ReLU} 或\gls{PReLU}，或者可以学习实现与这些都不同的函数。
maxout层的参数化当然也将与这些层不同，所以即使是maxout学习去实现和其他种类的层相同的$\Vx$的函数这种情况下，学习的机理也是不一样的。

每个~\gls{maxout_unit}现在由$k$个权重向量来参数化，而不仅仅是一个，所以~\gls{maxout_unit}通常比\gls{ReLU}需要更多的正则化。
如果训练集很大并且每个单元的块数保持很低的话，它们可以在没有正则化的情况下工作得不错\citep{cai2013deep}。

\gls{maxout_unit}还有一些其他的优点。
在某些情况下，要求更少的参数可以获得一些统计和计算上的优点。
具体来说，如果由$n$个不同的线性过滤器描述的特征可以在不损失信息的情况下，用每一组$k$个特征的最大值来概括的话，那么下一层可以获得$k$倍更少的权重数 。

因为每个单元由多个过滤器驱动，\gls{maxout_unit}具有一些冗余来帮助它们抵抗一种被称为\firstgls{catastrophic_forgetting}的现象，这个现象是说神经网络忘记了如何执行它们过去训练的任务\citep{Goodfellow+al-ICLR2014-small}。

\gls{ReLU}和它们的这些扩展都是基于一个原则，那就是如果它们的行为更接近线性，那么模型更容易优化。
使用线性行为更容易优化的一般性原则同样也适用于除深度线性网络以外的情景。
循环网络可以从序列中学习并产生状态和输出的序列。
当训练它们时，需要通过一些\gls{time_step}来传播信息，当其中包含一些线性计算（具有大小接近1的某些方向导数）时，这会更容易。
作为性能最好的循环网络结构之一，LSTM通过求和在时间上传播信息，这是一种特别直观的线性激活。
它将在\secref{sec:the_long_short_term_memory_and_other_gated_rnns}中进一步讨论。

% -- 188 --

\subsection{\gls{logistic_sigmoid}与双曲正切函数}
\label{sec:logistic_sigmoid_and_hyperbolic_tangent}

在引入\gls{ReLU}之前，大多数神经网络使用~\gls{logistic_sigmoid}~激活函数
\begin{equation}
g(z) = \sigma(z)
\end{equation}
或者是双曲正切激活函数
\begin{equation}
g(z) = \text{tanh}(z).
\end{equation}
这些激活函数紧密相关，因为$\text{tanh}(z)=2\sigma(2z)-1$。

我们已经看过~\gls{sigmoid}~单元作为输出单元用来预测二值型变量取值为1的概率。
与分段线性单元不同，\gls{sigmoid}~单元在其大部分定义域内都饱和——当$z$取绝对值很大的正值时，它们饱和到一个高值，当$z$取绝对值很大的负值时，它们饱和到一个低值，并且仅仅当$z$接近0时它们才对输入强烈敏感。
\gls{sigmoid}~单元的广泛饱和性会使得基于梯度的学习变得非常困难。
因为这个原因，现在不鼓励将它们用作前馈网络中的\gls{hidden_unit}。
当使用一个合适的代价函数来抵消~\gls{sigmoid}~的饱和性时，它们作为输出单元可以与基于梯度的学习相兼容。

当必须要使用~\gls{sigmoid}~激活函数时，双曲正切激活函数通常要比~\gls{logistic_sigmoid}~函数表现更好。
在$\text{tanh}(0)=0$而$\sigma(0)=\frac{1}{2}$的意义上，它更像是单位函数。
因为$\text{tanh}$在0附近与单位函数类似，训练深层神经网络$\hat{y}=\Vw^\top \text{tanh}(\MU^\top \text{tanh}(\MV^\top \Vx))$类似于训练一个线性模型$\hat{y}= \Vw^\top \MU^\top \MV^\top \Vx$，只要网络的激活能够被保持地很小。
这使得训练$\text{tanh}$网络更加容易。

% -- 189 --

\gls{sigmoid}~激活函数在除了\gls{feedforward_network}以外的情景中更为常见。
循环网络、许多概率模型以及一些\gls{AE}有一些额外的要求使得它们不能使用分段线性激活函数，并且使得~\gls{sigmoid}~单元更具有吸引力，尽管它存在饱和性的问题。

\subsection{其他\gls{hidden_unit}}
\label{sec:other_hidden_units}

也存在许多其他种类的\gls{hidden_unit}，但它们并不常用。

一般来说，很多种类的可微函数都表现得很好。
许多未发布的激活函数与流行的激活函数表现得一样好。
为了提供一个具体的例子，作者在MNIST数据集上使用$\Vh=\cos(\MW\Vx+\Vb)$测试了一个\gls{feedforward_network}，并获得了小于1\%的误差率，这可以与更为传统的激活函数获得的结果相媲美。
在新技术的研究和开发期间，通常会测试许多不同的激活函数，并且会发现许多标准方法的变体表现非常好。
这意味着，通常新的\gls{hidden_unit}类型只有在被明确证明能够提供显著改进时才会被发布。
新的\gls{hidden_unit}类型如果与已有的\gls{hidden_unit}表现大致相当的话，那么它们是非常常见的，不会引起别人的兴趣。

列出文献中出现的所有\gls{hidden_unit}类型是不切实际的。
我们只对一些特别有用和独特的类型进行强调。

其中一种是完全没有激活函数$g(z)$。
也可以认为这是使用单位函数作为激活函数的情况。
我们已经看过线性单元可以用作神经网络的输出。
它也可以用作\gls{hidden_unit}。
如果神经网络的每一层都仅由线性变换组成，那么网络作为一个整体也将是线性的。
然而，神经网络的一些层是纯线性也是可以接受的。
考虑具有$n$个输入和$p$个输出的神经网络层$\Vh=g(\MW^\top \Vx+\Vb)$。
我们可以用两层来代替它，一层使用权重矩阵$\MU$，另一层使用权重矩阵$\MV$。
如果第一层没有激活函数，那么我们（用$\MU$和$\MV$替代$\MW$）实际上是对基于$\MW$的原始层的权重矩阵进行了因式分解。
分解方法是计算$\Vh=g(\MV^\top \MU^\top \Vx+\Vb)$。
如果$U$产生了$q$个输出，那么$\MU$和$V$一起仅包含$(n+p)q$个参数，而$\MW$包含$np$个参数。
如果$q$很小，这可以在很大程度上节省参数。
这是以将线性变换约束为低秩的代价来实现的，但这些低秩关系往往是足够的。
线性\gls{hidden_unit}因此提供了一种减少网络中参数数量的有效方法。

% -- 190 --


softmax单元是另外一种经常用作输出的单元（如\secref{sec:softmax_units_for_multinoulli_output_distributions}中所描述的），但有时也可以用作\gls{hidden_unit}。
softmax单元很自然地表示具有$k$个可能值的离散型随机变量的概率分布，所以它们可以用作一种开关。
这些类型的\gls{hidden_unit}通常仅用于明确地学习操作内存的高级结构中，将在\secref{sec:explicit_memory}中描述。

其他一些常见的\gls{hidden_unit}类型包括：
\begin{itemize}
\item \firstall{RBF}：$h_i = \exp \left (-\frac{1}{\sigma_i^2}|| \MW_{:,i}-\Vx||^2 \right )$。
这个函数在$\Vx$接近模板$\MW_{:,i}$时更加活跃。
因为它对大部分$\Vx$都饱和到0，因此很难优化。

\item \textbf{\gls{softplus}}函数：$g(a)=\zeta(a)=\log(1+e^a)$。
这是\gls{ReLU}的平滑版本，由~\cite{Dugas01}引入用于函数近似，由~\cite{Nair-2010-small}引入用于无向概率模型的条件分布。
\cite{Glorot+al-AI-2011-small}比较了softplus和\gls{ReLU}，发现后者的结果更好。
通常不鼓励使用~\gls{softplus_function}。
softplus表明\gls{hidden_unit}类型的性能可能是非常反直觉的——因为它处处可导或者因为它不完全饱和，人们可能希望它具有优于\gls{ReLU}的点，但根据经验来看，它并没有。

\item \firstgls{hard_tanh}：它的形状和$\text{tanh}$以及\gls{ReLU}类似，但是不同于后者，它是有界的，$g(a)=\max(-1, \min(1,a))$。
它由~\cite{Collobert04}引入。
\end{itemize}

\gls{hidden_unit}的设计仍然是一个活跃的研究领域，许多有用的\gls{hidden_unit}类型仍有待发现。

\section{架构设计}
\label{sec:architecture_design}

神经网络设计的另一个关键点是确定它的架构。
\firstgls{architecture}一词是指网络的整体结构：它应该具有多少单元，以及这些单元应该如何连接。

% -- 191 --

大多数神经网络被组织成称为层的单元组。
大多数神经网络架构将这些层布置成链式结构，其中每一层都是前一层的函数。
在这种结构中，第一层由下式给出：
\begin{equation}
\Vh^{(1)}= g^{(1)}\left ( \MW^{(1)\top} \Vx + \Vb^{(1)}\right );
\end{equation}
第二层由
\begin{equation}
\Vh^{(2)} = g^{(2)}\left ( \MW^{(2)\top}\Vh^{(1)}+\Vb^{(2)} \right );
\end{equation}
给出，以此类推。

在这些链式架构中，主要的架构考虑是选择网络的深度和每一层的宽度。
我们将会看到，即使只有一个\gls{hidden_layer}的网络也足够适应训练集。
更深层的网络通常能够对每一层使用更少的单元数和更少的参数，并且经常容易泛化到测试集，但是通常也更难以优化。
对于一个具体的任务，理想的网络架构必须通过实验，观测在验证集上的误差来找到。

\subsection{万能近似性质和深度}
\label{sec:universal_approximation_properties_and_depth}

线性模型，通过矩阵乘法将特征映射到输出，顾名思义，仅能表示线性函数。
它具有易于训练的优点，因为当使用线性模型时，许多损失函数会导出凸优化问题。
可惜的是，我们经常希望我们的系统学习非线性函数。

乍一看，我们可能认为学习非线性函数需要为我们想要学习的那种非线性专门设计一类模型族。
幸运的是，具有\gls{hidden_layer}的\gls{feedforward_network}提供了一种万能近似框架。
具体来说，\firstgls{universal_approximation_theorem}\citep{Hornik89,Cybenko89}表明，一个\gls{feedforward_neural_network}如果具有线性输出层和至少一层具有任何一种``挤压''性质的激活函数（例如\gls{logistic_sigmoid}激活函数）的\gls{hidden_layer}，只要给予网络足够数量的\gls{hidden_unit}，它可以以任意的精度来近似任何从一个有限维空间到另一个有限维空间的Borel可测函数。
\gls{feedforward_network}的导数也可以任意好地来近似函数的导数\citep{hornik1990universal}。
Borel可测的概念超出了本书的范畴；对于我们想要实现的目标，只需要知道定义在$\SetR^n$的有界闭集上的任意连续函数是Borel可测的，因此可以用神经网络来近似。
神经网络也可以近似从任何有限维离散空间映射到另一个的任意函数。
虽然原始定理最初以具有特殊激活函数的单元的形式来描述，这个激活函数当变量取绝对值非常大的正值和负值时都会饱和，\gls{universal_approximation_theorem}也已经被证明对于更广泛类别的激活函数也是适用的，其中就包括现在常用的\gls{ReLU}~\citep{Leshno-et-al-1993}。

% -- 192 --

\gls{universal_approximation_theorem}意味着无论我们试图学习什么函数，我们知道一个大的MLP一定能够\emph{表示}这个函数。
然而，我们不能保证训练算法能够\emph{学得}这个函数。
即使MLP能够表示该函数，学习也可能因两个不同的原因而失败。
首先，用于训练的优化算法可能找不到用于期望函数的参数值。
其次，训练算法可能由于过拟合而选择了错误的函数。
回忆\secref{sec:the_no_free_lunch_theorem}中的``没有免费的午餐''定理，说明了没有普遍优越的机器学习算法。
\gls{feedforward_network}提供了表示函数的万能系统，在这种意义上，给定一个函数，存在一个\gls{feedforward_network}能够近似该函数。
不存在万能的过程能够通过检查训练集上的样本，来找到一个能对训练集上没有的点也适用的函数。

\gls{universal_approximation_theorem}说明了，存在一个足够大的网络能够达到我们所希望的任意精度，但是定理并没有说这个网络有多大。
\cite{Barron93}提供了单层网络近似一大类函数所需大小的一些界。
不幸的是，在最坏情况下，可能需要指数数量的\gls{hidden_unit}（可能一个\gls{hidden_unit}对应着一个需要区分的输入配置）。
这在二进制值的情况下很容易看到：向量$\Vv \in \{0,1\}^n$上的可能的二值型函数的数量是$2^{2^n}$，并且选择一个这样的函数需要$2^n$位，这通常需要$O(2^n)$的自由度。

总之，具有单层的\gls{feedforward_network}足以表示任何函数，但是网络层可能大得不可实现，并且可能无法正确地学习和\gls{generalization}。
在很多情况下，使用更深的模型能够减少表示期望函数所需的单元的数量，并且可以减少\gls{generalization}误差。

% -- 193 --

存在一些函数族能够被深度大于某个值$d$的网络高效地近似，而当深度被限制到小于或等于$d$时需要一个远远大于之前的模型。
在很多情况下，浅层模型所需的\gls{hidden_unit}的数量是$n$的指数级。
这个结果最初是在那些不与连续可微的神经网络类似的机器学习模型中被证明，但现在已经扩展到了这些模型。
第一个结果是关于逻辑门电路的\citep{Hastad86}。
后来的工作将这些结果扩展到了具有非负权重的线性阈值单元\citep{Hastad91,Hajnal-et-al-1993}，然后扩展到了具有连续值激活的网络\citep{Maass-1992,Maass-et-al-1994}。
许多现代神经网络使用\gls{ReLU}。
\cite{Leshno-et-al-1993}证明带有一大类非多项式激活函数族的浅层网络，包括\gls{ReLU}，具有万能的近似性质，但是这些结果并没有强调深度或效率的问题——它们仅指出足够宽的\gls{rectifier_network}能够表示任意函数。
\cite{Montufar-et-al-NIPS2014}指出一些用深度\gls{rectifier_network}表示的函数可能需要浅层网络（一个\gls{hidden_layer}）指数级的\gls{hidden_unit}才能表示。
更确切的说，他们说明分段线性网络（可以通过整流非线性或~\gls{maxout_unit}获得）可以表示区域的数量是网络深度的指数级的函数。%这里为什么是非线性？
\figref{fig:chap6_space_folding}解释了带有\gls{absolute_value_rectification}的网络是如何创建函数的镜像图像的，这些函数在某些\gls{hidden_unit}的顶部计算，作用于\gls{hidden_unit}的输入。
每个\gls{hidden_unit}指定在哪里折叠输入空间，来创造镜像响应（在绝对值非线性的两侧）。
通过组合这些折叠操作，我们获得指数级的分段线性区域，他们可以概括所有种类的规则模式（例如，重复）。
% fig 6.5
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics[width=\textwidth]{Chapter6/figures/space_folding.jpg}}
\fi
\caption{ % This is a screenshot.
关于更深的整流网络具有指数优势的一个直观的几何解释，来自~\cite{Montufar-et-al-NIPS2014}。 
\emph{(左)}\gls{absolute_value_rectification}单元对其输入中的每对镜像点有相同的输出。
镜像的对称轴由单元的权重和偏置定义的超平面给出。 在该单元顶部计算的函数（绿色决策面）将是横跨该对称轴的更简单模式的一个镜像。 
\emph{(中)}该函数可以通过折叠对称轴周围的空间来得到。
\emph{(右)}另一个重复模式可以在第一个的顶部折叠（由另一个下游单元）以获得另外的对称性（现在重复四次，使用了两个\gls{hidden_layer}）。 
经~\cite{Montufar-et-al-NIPS2014}许可改编此图。}
\label{fig:chap6_space_folding}
\end{figure}


\cite{Montufar-et-al-NIPS2014}的主要定理指出，具有$d$个输入、深度为$l$、每个\gls{hidden_layer}具有$n$个单元的深度\gls{rectifier_network}可以描述的线性区域的数量是
\begin{equation}
O \left ( {n \choose d}^{d(l-1)} n^d \right ),
\end{equation}
意味着，这是深度$l$的指数级。
在每个单元具有$k$个过滤器的maxout网络中，线性区域的数量是
\begin{equation}
O \left ( k^{(l-1)+d} \right ).
\end{equation}

% -- 194 --

当然，我们不能保证在机器学习（特别是AI）的应用中我们想要学得的函数类型享有这样的属性。

我们还可能出于统计原因来选择深度模型。
任何时候，当我们选择一个特定的机器学习算法时，我们隐含地陈述了一些先验，这些先验是关于算法应该学得什么样的函数的。
选择深度模型默许了一个非常普遍的信念，那就是我们想要学得的函数应该涉及几个更加简单的函数的组合。
这可以从表示学习的观点来解释， 我们相信学习的问题包含发现一组潜在的\gls{factors_of_variation}，它们可以根据其他更简单的潜在的\gls{factors_of_variation}来描述。
或者，我们可以将深度结构的使用解释为另一种信念，那就是我们想要学得的函数是包含多个步骤的计算机程序，其中每个步骤使用前一步骤的输出。
这些中间输出不一定是\gls{factors_of_variation}，而是可以类似于网络用来组织其内部处理的计数器或指针。
根据经验，更深的模型似乎确实在广泛的任务中泛化得更好\citep{Bengio-NIPS2007,Erhan2009,Bengio-2009-book,UTLC+LISA-2011-small,Ciresan-et-al-2012,Krizhevsky-2012-small,sermanet-cvpr-13,Farabet-et-al-2013,couprie-iclr-13,LISA-EmotiW2013,Goodfellow+et+al-ICLR2014a,Szegedy-et-al-arxiv2014}。
\figref{fig:chap6_depth_color}和\figref{fig:chap6_model_size_color}展示了一些实验结果的例子。
这表明使用深层架构确实在模型学习的函数空间上表示了一个有用的先验。
% fig 6.6
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/depth_color}}
\fi
\caption{深度的影响。
实验结果表明，当从地址照片转录多位数字时，更深层的网络能够更好地\gls{generalization}。
数据来自~\cite{Goodfellow+et+al-ICLR2014a}。
测试集上的准确率随着深度的增加而不断增加。
\figref{fig:chap6_model_size_color}给出了一个对照实验，它说明了对模型尺寸其他方面的增加并不能产生相同的效果。}
\label{fig:chap6_depth_color}
\end{figure}

% fig 6.7
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/model_size_color}}
\fi
\caption{参数数量的影响。
更深的模型往往表现更好。
这不仅仅是因为模型更大。
\cite{Goodfellow+et+al-ICLR2014a}的这项实验表明，增加\gls{convolutional_network}层中参数的数量，但是不增加它们的深度，在提升测试集性能方面没有同样显著的效果，如此图所示。
图例标明了每条曲线对应的网络深度，以及曲线表示的是卷积层还是全连接层的大小变化。
我们可以观察到，在这种情况下，浅层模型在参数数量达到2000万时就过拟合，而深层模型在参数数量超过6000万时仍然表现良好。
这表明，使用深层模型表达出了对模型可以学习的函数空间的有用偏好。
具体来说，它表达了一种信念，即该函数应该由许多更简单的函数复合在一起而得到。
这可能导致学习由更简单的表示所组成的表示（例如，由边所定义的角）或者学习具有顺序依赖步骤的程序（例如，首先定位一组对象，然后分割它们，之后识别它们）。}
\label{fig:chap6_model_size_color}
\end{figure}

\subsection{其他架构上的考虑}
\label{sec:other_architectural_considerations}

目前为止，我们都将神经网络描述成层的简单链式结构，主要的考虑因素是网络的深度和每层的宽度。
在实践中，神经网络显示出相当的多样性。

许多神经网络架构已经被开发用于特定的任务。
用于计算机视觉的卷积神经网络的特殊架构将在\chapref{chap:convolutional_networks}中介绍。
\gls{feedforward_network}也可以推广到用于序列处理的循环神经网络，但有它们自己的架构考虑，将在\chapref{chap:sequence_modeling_recurrent_and_recursive_nets}中介绍。

% -- 195 --

一般的，层并不一定都需要连接在一条链中，尽管这是最常见的做法。
许多架构构建了一个主链，但随后又添加了额外的架构特性，例如从层$i$到层$i+2$或者更高层的跳跃连接。
这些跳跃连接使得梯度更容易从输出层流向更接近输入的层。

架构设计考虑的另外一个关键点是如何将层与层之间连接起来。
默认的神经网络层采用矩阵$\MW$描述的线性变换，每个输入单元连接到每个输出单元。
在之后章节中的许多专用网络具有较少的连接，使得输入层中的每个单元仅连接到输出层单元的一个小子集。
这些用于减少连接数量的策略减少了参数的数量以及用于评估网络的计算量，但通常高度依赖于问题。
例如，\chapref{chap:convolutional_networks}描述的卷积神经网络使用对于计算机视觉问题非常有效的稀疏连接的专用模式。
在这一章中，很难对通用神经网络的架构给出更多具体的建议。
我们在随后的章节中介绍一些特殊的架构策略，可以在不同的领域工作良好。

% -- 196 --

\section{\glsentrytext{BP}和其他的微分算法}
\label{sec:back_propagation_and_other_differentiation_algorithms}

当我们使用\gls{feedforward_neural_network}接收输入$\Vx$并产生输出$\hat{\Vy}$时，信息通过网络向前流动。
输入$\Vx$提供初始信息，然后传播到每一层的\gls{hidden_unit}，最终产生输出$\hat{\Vy}$。
这称之为~\firstgls{forward_propagation}。
在训练过程中，前向传播可以持续向前直到它产生一个标量代价函数$J(\Vtheta)$。
\firstgls{back_propagation}算法\citep{RHW}，经常简称为\textbf{\glsdesc{BP}}，允许来自代价函数的信息通过网络向后流动，以便计算梯度。

% -- 197 --

计算梯度的解析表达式是很直观的，但是数值化地求解这样的表达式在计算上的代价可能很大。
反向传播算法使用简单和廉价的程序来实现这个目标。

反向传播这个术语经常被误解为用于多层神经网络的整个学习算法。
实际上，反向传播仅指用于计算梯度的方法，而另一种算法，例如随机梯度下降，使用该梯度来进行学习。
此外，反向传播经常被误解为仅适用于多层神经网络，但是原则上它可以计算任何函数的导数（对于一些函数，正确的响应是报告函数的导数是未定义的）。
特别地，我们会描述如何计算一个任意函数$f$的梯度$\nabla_{\Vx} f(\Vx, \Vy)$，其中$\Vx$是一组变量，我们需要它们的导数，而$\Vy$是函数的另外一组输入变量，但我们并不需要它们的导数。
在学习算法中，我们最常需要的梯度是代价函数关于参数的梯度，即$\nabla_{\Vtheta} J(\Vtheta)$。
许多机器学习任务需要计算其他导数，来作为学习过程的一部分，或者用来分析学得的模型。
反向传播算法也适用于这些任务，不局限于计算代价函数关于参数的梯度。
通过在网络中传播信息来计算导数的想法非常普遍，它还可以用于计算诸如多输出函数$f$的~\gls{jacobian}~的值。
我们这里描述的是最常用的情况，其中$f$只有单个输出。

\subsection{计算图}
\label{sec:computational_graphs}

目前为止，我们已经用相对非正式的图形语言讨论了神经网络。
为了更精确地描述反向传播算法，使用更精确的\firstgls{computational_graph}语言是很有帮助的。

将计算形式化为图形的方法有很多。

这里，我们使用图中的每一个节点来表示一个变量。
变量可以是标量、向量、矩阵、张量、或者甚至是另一类型的变量。

为了形式化我们的图形，我们还需引入\firstgls{operation}这一概念。
操作是指一个或多个变量的简单函数。
我们的图形语言伴随着一组被允许的操作。
我们可以通过将多个操作复合在一起来描述更为复杂的函数。

% -- 198 --

不失一般性，我们定义一个操作仅返回单个输出变量。
这并没有失去一般性，是因为输出变量可以有多个条目，例如向量。
反向传播的软件实现通常支持具有多个输出的操作，但是我们在描述中避免这种情况，因为它引入了对概念理解不重要的许多额外细节。

如果变量$y$是变量$x$通过一个操作计算得到的，那么我们画一条从$x$到$y$的有向边。
我们有时用操作的名称来注释输出的节点，当上下文很明确时，有时也会省略这个标注。

计算图的实例可以参考\figref{fig:chap6_computational_graph}。
% fig 6.8
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/computational_graph}}
\fi
\captionsetup{singlelinecheck=off}
\caption{一些计算图的示例。
\emph{(a)}使用$\times$操作计算$z = xy$的图。 
\emph{(b)}用于\gls{logistic_regression}预测$\hat{y} = \sigma(\Vx^\top \Vw + b)$的图。 一些中间表达式在代数表达式中没有名称，但在图形中却需要。
我们简单地将第$i$个这样的变量命名为$\Vu^{(i)}$。
\emph{(c)}表达式$\MH = \max \{ 0, \MX\MW+ \Vb \}$的计算图，在给定包含小批量输入数据的设计矩阵$\MX$时，它计算\gls{ReLU}激活的设计矩阵$\MH$。
\emph{(d)}示例a-c对每个变量最多只实施一个操作，但是对变量实施多个操作也是可能的。 这里我们展示一个计算图，它对线性回归模型的权重$\Vw$实施多个操作。
这个权重不仅用于预测$\hat{y}$，也用于权重衰减罚项$\lambda \sum_i w_i^2$。}
\label{fig:chap6_computational_graph}
\end{figure}

\subsection{微积分中的链式法则}
\label{sec:chain_rule_of_calculus}

微积分中的链式法则（为了不与概率中的链式法则相混淆）用于计算复合函数的导数。
反向传播是一种计算链式法则的算法，使用高效的特定运算顺序。

设$x$是实数，$f$和$g$是从实数映射到实数的函数。
假设$y=g(x)$并且$z=f(g(x))=f(y)$。
那么链式法则是说
\begin{equation}
\frac{dz}{dx}=\frac{dz}{dy} \frac{dy}{dx}.
\label{eq:6.44}
\end{equation}

我们可以将这种标量情况进行扩展。
假设$\Vx\in \SetR^m, \Vy\in \SetR^n$，$g$是从$\SetR^m$到$\SetR^n$的映射，$f$是从$\SetR^n$到$\SetR$的映射。
如果$\Vy=g(\Vx)$并且$z=f(\Vy)$，那么
\begin{equation}
\frac{\partial z}{\partial x_i} = \sum_j \frac{\partial z}{\partial y_j} \frac{\partial y_j}{\partial x_i}.
\end{equation}
使用向量记法，可以等价地写成
\begin{equation}
\nabla_{\Vx}z = \left ( \frac{\partial \Vy}{\partial \Vx} \right )^\top \nabla_{\Vy} z,
\end{equation}
这里$\frac{\partial \Vy}{\partial \Vx}$是$g$的$n\times m$的~\gls{jacobian_matrix}。

从这里我们看到，变量$\Vx$的梯度可以通过~\gls{jacobian_matrix}$\frac{\partial \Vy}{\partial \Vx}$和梯度$\nabla_{\Vy} z$相乘来得到。
反向传播算法由图中每一个这样的~\gls{jacobian}~梯度的乘积操作所组成。

% -- 199 --

通常我们将反向传播算法应用于任意维度的张量，而不仅仅用于向量。
从概念上讲，这与使用向量的反向传播完全相同。 
唯一的区别是如何将数字排列成网格以形成张量。 
我们可以想象，在我们运行反向传播之前，将每个张量变平为一个向量，计算一个向量值梯度，然后将该梯度重新构造成一个张量。
从这种重新排列的观点上看，反向传播仍然只是将~\gls{jacobian}~乘以梯度。

% -- 200 --

为了表示值$z$关于张量$\TSX$的梯度，我们记为$\nabla_\TSX z$，就像$\TSX$是向量一样。
$\TSX$的索引现在有多个坐标——例如，一个3维的张量由三个坐标索引。
我们可以通过使用单个变量$i$来表示完整的索引元组，从而完全抽象出来。
对所有可能的元组$i$，$(\nabla_\TSX z)_i$给出$\frac{\partial z}{\partial \TEX_i}$。
这与向量中索引的方式完全一致，$(\nabla_{\Vx} z)_i$给出$\frac{\partial z}{\partial x_i}$。
使用这种记法，我们可以写出适用于张量的链式法则。
如果$\TSY=g(\TSX)$并且$z=f(\TSY)$，那么
\begin{equation}
  \nabla_\TSX z = \sum_j (\nabla_\TSX \TEY_j)\frac{\partial z}{\partial \TEY_j}.
  \label{eq:6.47}
\end{equation}

\subsection{递归地使用链式法则来实现\glsentrytext{BP}}
\label{sec:recursively_applying_the_chain_rule_to_obtain_backprop}

使用链式规则，我们可以直接写出某个标量关于计算图中任何产生该标量的节点的梯度的代数表达式。
然而，实际在计算机中计算该表达式时会引入一些额外的考虑。

具体来说，许多子表达式可能在梯度的整个表达式中重复若干次。
任何计算梯度的程序都需要选择是存储这些子表达式还是重新计算它们几次。
\figref{fig:chap6_repeated_subexpression}给出了一个例子来说明这些重复的子表达式是如何出现的。
在某些情况下，计算两次相同的子表达式纯粹是浪费。
在复杂图中，可能存在指数多的这种计算上的浪费，使得简单的链式法则不可实现。
在其他情况下，计算两次相同的子表达式可能是以较高的运行时间为代价来减少内存开销的有效手段。

我们首先给出一个版本的反向传播算法，它指明了梯度的直接计算方式（\algref{alg:bprop}以及相关的正向计算的\algref{alg:fprop}），按照它实际完成的顺序并且递归地使用链式法则。
我们可以直接执行这些计算或者将算法的描述视为用于计算反向传播的计算图的符号表示。
然而，这些公式并没有明确地操作和构造用于计算梯度的符号图。
这些公式将在后面的\secref{sec:general_back_propagation}和\algref{alg:backprop}中给出，其中我们还推广到了包含任意张量的节点。

首先考虑描述如何计算单个标量$u^{(n)}$（例如训练样本上的损失函数）的计算图。
我们想要计算这个标量对$n_i$个输入节点$u^{(1)}$到$u^{(n_i)}$的梯度。
换句话说，我们希望对所有的$i\in\{1,2,\ldots,n_i\}$计算$\frac{\partial u^{(n)}}{\partial u^{(i)}}$。
在使用反向传播计算梯度来实现参数的梯度下降时，$u^{(n)}$将对应单个或者小批量实例的代价函数，而$u^{(1)}$到$u^{(n_i)}$则对应于模型的参数。

% -- 201 --

我们假设图的节点已经以一种特殊的方式被排序，使得我们可以一个接一个地计算他们的输出，从$u^{(n_i+1)}$开始，一直上升到$u^{(n)}$。
如\algref{alg:fprop}中所定义的，每个节点$u^{(i)}$与操作$f^{(i)}$相关联，并且通过对以下函数求值来得到
\begin{equation}
  u^{(i)} = f(\SetA^{(i)}),
\end{equation}
其中$\SetA^{(i)}$是$u^{(i)}$所有父节点的集合。
% alg 6.1
\begin{algorithm}[htbp]
\caption{计算将$n_i$个输入$u^{(1)}$到$u^{(n_i)}$映射到一个输出$u^{(n)}$的程序。
这定义了一个计算图，其中每个节点通过将函数$f^{(i)}$应用到变量集合$\SetA^{(i)}$上来计算$u^{(i)}$的值，$\SetA^{(i)}$包含先前节点$u^{(j)}$的值满足$j<i$且$j \in Pa(u^{(i)})$。
计算图的输入是向量$\Vx$，并且被分配给前$n_i$个节点$u^{(1)}$到$u^{(n_i)}$。计算图的输出可以从最后一个（输出）节点$u^{(n)}$读出。}
\label{alg:fprop}
\begin{algorithmic}
\FOR {$i=1, \ldots, n_i$}
 \STATE $u^{(i)} \leftarrow x_i$
\ENDFOR
\FOR {$i=n_i+1, \ldots, n$}
 \STATE $\SetA^{(i)} \leftarrow \{ u^{(j)} \mid j \in Pa(u^{(i)}) \}$
 \STATE $u^{(i)} \leftarrow f^{(i)}(\SetA^{(i)})$
\ENDFOR
\STATE {\bf return} $u^{(n)}$
\end{algorithmic}
\end{algorithm}

% -- 202 --

该算法详细说明了前向传播的计算，我们可以将其放入图$\CalG$中。
为了执行反向传播，我们可以构造一个依赖于$\CalG$并添加额外一组节点的计算图。
这形成了一个子图$\CalB$，它的每个节点都是$\CalG$的节点。
$\CalB$中的计算和$\CalG$中的计算顺序完全相反，而且$\CalB$中的每个节点计算导数$\frac{\partial u^{(n)}}{\partial u^{(i)}}$与前向图中的节点$u^{(i)}$相关联。
这通过对标量输出$u^{(n)}$使用链式法则来完成：
\begin{equation}
  \frac{\partial u^{(n)}}{\partial u^{(j)}} = \sum_{i:j \in Pa(u^{(i)})} \frac{\partial u^{(n)} }{ \partial u^{(i)} } \frac{ \partial u^{(i)} }{ \partial u^{(j)} }
  \label{eq:6.49}
\end{equation}
这在\algref{alg:bprop}中详细说明。
子图$\CalB$恰好包含每一条对应着$\CalG$中从节点$u^{(j)}$到节点$u^{(i)}$的边。
从$u^{(j)}$到$u^{(i)}$的边对应着计算$\frac{\partial u^{(i)}}{\partial u^{(j)}}$。
另外，对于每个节点都要执行一个内积，内积的一个因子是对于$u^{j}$子节点$u^{(i)}$的已经计算的梯度，另一个因子是对于相同子节点$u^{(i)}$ 的偏导数$\frac{\partial u^{(i)}}{\partial u^{(j)}}$组成的向量。
总而言之，执行反向传播所需的计算量与$\CalG$中的边的数量成比例，其中每条边的计算包括计算偏导数（节点关于它的一个父节点的偏导数）以及执行一次乘法和一次加法。
下面，我们将此分析推广到张量值节点，这只是在同一节点中对多个标量值进行分组并能够更高效地实现。
% fig 6.9
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/repeated_subexpression}}
\fi
\captionsetup{singlelinecheck=off}
\caption[.]{
计算梯度时导致重复子表达式的计算图。
令$w \in \SetR$为图的输入。
我们对链中的每一步使用相同的操作函数$f: \SetR \to \SetR$，这样$x=f(w), y=f(x), z=f(y)$。
为了计算$\frac{\partial z}{\partial w}$，我们应用\eqnref{eq:6.44}得到：
\begin{align}
& \frac{\partial z}{\partial w}\\
=& \frac{\partial z}{\partial y} \frac{\partial y}{\partial x} \frac{\partial x}{\partial w}\\
\label{eq:6.52}
=& f'(y)f'(x)f'(w)\\ 
\label{eq:6.53}
=& f'(f(f(w))) f'(f(w)) f'(w). 
\end{align}
\eqnref{eq:6.52}建议我们采用的实现方式是，仅计算$f(w)$的值一次并将它存储在变量$x$中。
这是\gls{BP}算法所采用的方法。
\eqnref{eq:6.53}提出了一种替代方法，其中子表达式$f(w)$出现了不止一次。
在替代方法中，每次只在需要时重新计算$f(w)$。
当存储这些表达式的值所需的存储较少时，\eqnref{eq:6.52}的\gls{BP}方法显然是较优的，因为它减少了运行时间。
然而，\eqnref{eq:6.53}也是链式法则的有效实现，并且当存储受限时它是有用的。}
\label{fig:chap6_repeated_subexpression}
\end{figure}

% -- 203 --

反向传播算法被设计为减少公共子表达式的数量而不考虑存储的开销。
具体来说，它大约对图中的每个节点执行一个~\gls{jacobian}~乘积。
这可以从\algref{alg:bprop}中看出，反向传播算法访问了图中的节点$u^{(j)}$到节点$u^{(i)}$的每条边一次，以获得相关的偏导数$\frac{\partial u^{(i)}}{\partial u^{(j)}}$。
反向传播因此避免了重复子表达式的指数爆炸。
然而，其他算法可能通过对计算图进行简化来避免更多的子表达式，或者也可能通过重新计算而不是存储这些子表达式来节省内存。
我们将在描述完反向传播算法本身后再重新审视这些想法。
% alg 6.2
\begin{algorithm}[htb!]
\caption{\gls{BP}算法的简化版本，用于计算$u^{(n)}$关于图中变量的导数。
这个示例旨在通过演示所有变量都是标量的简化情况来进一步理解\gls{BP}算法，这里我们希望计算关于$u^{(1)},\ldots,u^{(n_i)}$的导数。
这个简化版本计算了关于图中所有节点的导数。
假定与每条边相关联的偏导数计算需要恒定的时间的话，该算法的计算成本与图中边的数量成比例。
这与\gls{forward_propagation}的计算次数具有相同的阶。
每个$\frac{\partial u^{(i)}}{\partial u^{(j)}}$是$u^{(i)}$的父节点$u^{(j)}$的函数，从而将前向图的节点链接到\gls{BP}图中添加的节点。}
\label{alg:bprop}
\begin{algorithmic}
\STATE 运行前向传播 (对于此例是\algref{alg:fprop}) 获得网络的激活。
\STATE 初始化 {\tt grad\_table}，用于存储计算好的导数的数据结构。 ${\tt grad\_table}[u^{(i)}]$将存储$\frac{\partial u^{(n)}}{\partial u^{(i)}}$计算好的值。
\STATE ${\tt grad\_table}[u^{(n)}] \leftarrow 1$
\FOR {$j=n-1$ down to 1}
\STATE 下一行使用存储的值计算 $\frac{\partial u^{(n)}}{\partial u^{(j)}} =
  \sum_{i: j \in Pa(u^{(i)})} \frac{\partial u^{(n)}}{\partial u^{(i)}} \frac{\partial u^{(i)}}{\partial u^{(j)}}$：
\STATE ${\tt grad\_table}[ u^{(j)}] \leftarrow 
\sum_{i: j \in Pa(u^{(i)})} {\tt grad\_table}[u^{(i)}]
\frac{\partial u^{(i)}}{\partial u^{(j)}}$
\ENDFOR
\STATE {\bf return} $\{ {\tt grad\_table}[ u^{(i)}] \mid i=1, \dots, n_i \} $
\end{algorithmic}
\end{algorithm}


% -- 204 --

\subsection{全连接MLP中的\glsentrytext{BP}计算}
\label{sec:back_propagation_computation_in_fully_connected_mlp}

为了阐明反向传播的上述定义，让我们考虑一个与全连接的多层MLP相关联的特定图。

\algref{alg:mlp-fprop}首先给出了前向传播，它将参数映射到与单个训练样本（输入，目标）$(\Vx,\Vy)$相关联的监督损失函数$L(\hat{\Vy}, \Vy)$，其中$\hat{\Vy}$是当$\Vx$提供输入时神经网络的输出。
% alg 6.3
\begin{algorithm}[ht]
\caption{典型深度神经网络中的\gls{forward_propagation}和代价函数的计算。
损失函数$L(\hat{\Vy}, \Vy)$取决于输出$\hat{\Vy}$和目标$\Vy$（参考\secref{sec:learning_conditional_distributions_with_maximum_likelihood}中损失函数的示例）。
为了获得总代价$J$，损失函数可以加上正则项$\Omega(\theta)$，其中$\theta$包含所有参数（权重和偏置）。
\algref{alg:mlp-bprop}说明了如何计算$J$关于参数$\MW$和$\Vb$的梯度。 为简单起见，该演示仅使用单个输入样本$\Vx$。
实际应用应该使用\gls{minibatch}。
请参考\secref{sec:example_back_propagation_for_mlp_training}以获得更加真实的演示。}
\label{alg:mlp-fprop}
\begin{algorithmic}
\REQUIRE 网络深度， $l$
\REQUIRE $\MW^{(i)}, i \in \{ 1, \dots, l\}$， 模型的权重矩阵
\REQUIRE $\Vb^{(i)}, i \in \{ 1, \dots, l\}$， 模型的偏置参数
\REQUIRE $\Vx$，程序的输入
\REQUIRE $\Vy$，目标输出
\STATE $\Vh^{(0)}=\Vx$
\FOR {$k=1, \ldots, l$}
 \STATE $\Va^{(k)} = \Vb^{(k)} + \MW^{(k)} \Vh^{(k-1)}$
 \STATE $\Vh^{(k)} = f(\Va^{(k)})$
\ENDFOR
\STATE $\hat{\Vy} = \Vh^{(l)}$
\STATE $J = L(\hat{\Vy},\Vy) + \lambda \Omega(\theta)$
\end{algorithmic}
\end{algorithm}

\algref{alg:mlp-bprop}随后说明了将反向传播应用于该图所需的相关计算。
% alg 6.4
\begin{algorithm}[htbp]
\caption{深度神经网络中\algref{alg:mlp-fprop}的反向计算，它不止使用了输入$\Vx $和目标$\Vy$。
该计算对于每一层$k$都产生了对激活$\Va^{(k)}$的梯度，从输出层开始向后计算一直到第一个\gls{hidden_layer}。
这些梯度可以看作是对每层的输出应如何调整以减小误差的指导，根据这些梯度可以获得对每层参数的梯度。
权重和偏置上的梯度可以立即用作随机梯度更新的一部分（梯度算出后即可执行更新），或者与其他基于梯度的优化方法一起使用。
}
\label{alg:mlp-bprop}
\begin{algorithmic}
\STATE 在前向计算完成后，计算顶层的梯度：
\STATE $\Vg \leftarrow \nabla_{\hat{\Vy}} J = \nabla_{\hat{\Vy}} L(\hat{\Vy},\Vy)$
\FOR {$k=l, l-1, \dots, 1$}
\STATE 将关于层输出的梯度转换为非线性激活输入前的梯度（如果 $f$ 是逐元素的，则逐元素地相乘）：
  \STATE $\Vg \leftarrow \nabla_{\Va^{(k)}} J = \Vg \odot f'(\Va^{(k)})$
  \STATE 计算关于权重和偏置的梯度（如果需要的话，还要包括正则项）：
  \STATE $\nabla_{\Vb^{(k)}} J = \Vg + \lambda \nabla_{\Vb^{(k)}} \Omega(\theta)$
  \STATE $\nabla_{\MW^{(k)}} J = \Vg \; \Vh^{(k-1)\top} + \lambda \nabla_{\MW^{(k)}} \Omega(\theta)$
  \STATE 关于下一更低层的\gls{hidden_layer}传播梯度：
  \STATE $\Vg \leftarrow \nabla_{\Vh^{(k-1)}} J = \MW^{(k)\top} \; \Vg$
\ENDFOR
\end{algorithmic}
\end{algorithm}

\algref{alg:mlp-fprop}和\algref{alg:mlp-bprop}是简单而直观的演示。
然而，它们专门针对特定的问题。

现在的软件实现基于之后\secref{sec:general_back_propagation}中描述的一般形式的反向传播，它可以通过显式地操作表示符号计算的数据结构，来适应任何计算图。


% -- 205 --

\subsection{符号到符号的导数}
\label{sec:symbol_to_symbol_derivatives}

代数表达式和计算图都对\firstgls{symbol}或不具有特定值的变量进行操作。
这些代数或者基于图的表达式被称为\firstgls{symbolic_representation}。
当我们实际使用或者训练神经网络时，我们必须给这些符号赋特定的值。
我们用一个特定的\firstgls{numeric_value}来替代网络的符号输入$\Vx$，例如$[1.2, 3,765, -1.8]^\top$。

一些反向传播的方法采用计算图和一组用于图的输入的数值，然后返回在这些输入值处梯度的一组数值。
我们将这种方法称为\,\textbf{符号到数值}的微分。
这种方法用在诸如Torch~\citep{Torch-2011}和Caffe~\citep{Jia13caffe}之类的库中。

% -- 206 --
 
另一种方法是采用计算图以及添加一些额外的节点到计算图中，这些额外的节点提供了我们所需导数的符号描述。
这是Theano~\citep{bergstra+al:2010-scipy-small,Bastien-Theano-2012}和TensorFlow~\citep{tensorflow}所采用的方法。
\figref{fig:chap6_symbol_to_symbol}给出了该方法如何工作的一个例子。
这种方法的主要优点是导数可以使用与原始表达式相同的语言来描述。
因为导数只是另外一张计算图，我们可以再次运行反向传播，对导数再进行求导就能得到更高阶的导数。
高阶导数的计算在\secref{sec:higher_order_derivatives}中描述。
% fig 6.10
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/symbol_to_symbol}}
\fi
\captionsetup{singlelinecheck=off}
\caption{使用符号到符号的方法计算导数的示例。
在这种方法中，\gls{BP}算法不需要访问任何实际的特定数值。
相反，它将节点添加到计算图中来描述如何计算这些导数。
通用图形求值引擎可以在随后计算任何特定数值的导数。
\emph{(左)}在这个例子中，我们从表示$z=f(f(f(w)))$的图开始。
\emph{(右)}我们运行\gls{BP}算法，指导它构造表达式$\frac{dz}{dw}$对应的图。 在这个例子中，我们不解释\gls{BP}算法如何工作。
我们的目的只是说明想要的结果是什么：符号描述的导数的计算图。}
\label{fig:chap6_symbol_to_symbol}
\end{figure}

我们将使用后一种方法，并且使用构造导数的计算图的方法来描述反向传播算法。
图的任意子集之后都可以使用特定的数值来求值。
这允许我们避免精确地指明每个操作应该在何时计算。
相反，通用的图计算引擎只要当一个节点的父节点的值都可用时就可以进行求值。

% -- 207 --
  
基于符号到符号的方法的描述包含了符号到数值的方法。
符号到数值的方法可以理解为执行了与符号到符号的方法中构建图的过程中完全相同的计算。
关键的区别是符号到数值的方法不向外暴露出计算图。

\subsection{一般化的\glsentrytext{BP}}
\label{sec:general_back_propagation}

反向传播算法非常简单。
为了计算某个标量$z$关于图中它的一个祖先$\Vx$的梯度，我们首先观察到它关于$z$的梯度由$\frac{dz}{dz}=1$给出。
然后，我们可以计算对图中$z$的每个父节点的梯度，通过现有的梯度乘以产生$z$的操作的~\gls{jacobian}。
我们继续乘以~\gls{jacobian}，以这种方式向后穿过图，直到我们到达$\Vx$。
对于从$z$出发可以经过两个或更多路径向后行进而到达的任意节点，我们简单地对该节点来自不同路径上的梯度进行求和。

更正式地，图$\CalG$中的每个节点对应着一个变量。
为了实现最大的一般化，我们将这个变量描述为一个张量$\TSV$。
张量通常可以具有任意维度，并且包含标量、向量和矩阵。

我们假设每个变量$\TSV$与下列子程序相关联：
\begin{itemize}
    \item \verb|get_operation|($\TSV$)：它返回用于计算$\TSV$的操作，代表了在计算图中流入$\TSV$的边。
    例如，可能有一个Python或者C++的类表示矩阵乘法操作，以及~\verb|get_operation|函数。
    假设我们的一个变量是由矩阵乘法产生的，$\MC=\MA\MB$。
    那么，\verb|get_operation|($\TSV$)返回一个指向相应C++类的实例的指针。

    \item \verb|get_consumers|($\TSV, \CalG$)：它返回一组变量，是计算图$\CalG$中$\TSV$的子节点。

    \item \verb|get_inputs|($\TSV, \CalG$)：它返回一组变量，是计算图$\CalG$中$\TSV$的父节点。
\end{itemize}

% -- 208 --
  
每个操作~\verb|op|也与~\verb|bprop|操作相关联。
该~\verb|bprop|操作可以计算如\eqnref{eq:6.47}所描述的~\gls{jacobian}~向量积。
这是反向传播算法能够实现很大通用性的原因。
每个操作负责了解如何通过它参与的图中的边来反向传播。
例如，我们可以使用矩阵乘法操作来产生变量$\MC=\MA\MB$。
假设标量$z$关于$\MC$的梯度是$\MG$。
矩阵乘法操作负责定义两个反向传播规则，每个规则对应于一个输入变量。
如果我们调用~\verb|bprop|方法来请求关于$\MA$的梯度，那么在给定输出的梯度为$\MG$的情况下，矩阵乘法操作的~\verb|bprop|方法必须说明关于$\MA$的梯度是$\MG\MB^\top$。
类似的，如果我们调用~\verb|bprop|方法来请求关于$\MB$的梯度，那么矩阵操作负责实现~\verb|bprop|方法并指定希望的梯度是$\MA^\top\MG$。
反向传播算法本身并不需要知道任何微分法则。
它只需要使用正确的参数调用每个操作的~\verb|bprop|方法即可。
正式地，~\verb|op.bprop(inputs|$, \TSX, \TSG)$必须返回
\begin{equation}
  \sum_i (\nabla_{\TSX} \verb|op.f(inputs|)_i) \textsf{G}_i,
\end{equation}
这只是如\eqnref{eq:6.47}所表达的链式法则的实现。
这里，\verb|inputs|是提供给操作的一组输入，\verb|op.f|是操作实现的数学函数，$\TSX$是输入，我们想要计算关于它的梯度，$\TSG$是操作对于输出的梯度。

\verb|op.bprop|方法应该总是假装它的所有输入彼此不同，即使它们不是。
例如，如果~\verb|mul|操作传递两个$x$来计算$x^2$，\verb|op.bprop|方法应该仍然返回$x$作为对于两个输入的导数。
反向传播算法后面会将这些变量加起来获得$2x$，这是$x$上总的正确的导数。

反向传播算法的软件实现通常提供操作和其~\verb|bprop|方法，所以深度学习软件库的用户能够对使用诸如矩阵乘法、指数运算、对数运算等等常用操作构建的图进行反向传播。
构建反向传播新实现的软件工程师或者需要向现有库添加自己的操作的高级用户通常必须手动为新操作推导~\verb|op.bprop|方法。

反向传播算法的正式描述参考\algref{alg:backprop}。

% -- 209 --
% alg 6.5
\begin{algorithm}[ht]
\caption{\gls{BP}算法最外围的骨架。
这部分做简单的设置和清理工作。
大多数重要的工作发生在\algref{alg:build_grad}的子程序{\tt build\_grad}中。
}
\label{alg:backprop}
\begin{algorithmic}
\REQUIRE $\SetT$，需要计算梯度的目标变量集
\REQUIRE $\CalG$，计算图
\REQUIRE $z$， 要微分的变量
\STATE 令 $\CalG'$ 为$\CalG$剪枝后的计算图，其中仅包括$z$的祖先以及$\SetT$中节点的后代。
\STATE 初始化 {\tt grad\_table}，它是关联张量和对应导数的数据结构。
\STATE ${\tt grad\_table}[z] \leftarrow 1$
\FOR{$\TSV$ in $\SetT$}
\STATE ${\tt build\_grad}(\TSV, \CalG, \CalG', {\tt grad\_table})$
\ENDFOR
\STATE Return {\tt grad\_table} restricted to $\SetT$
\end{algorithmic}
\end{algorithm}

\begin{algorithm}[ht]
\caption{\gls{BP}算法的内循环子程序${\tt build\_grad}(\TSV, \CalG, \CalG', {\tt grad\_table})$，
由\algref{alg:backprop}中定义的\gls{BP}算法调用。
}
\label{alg:build_grad}
\begin{algorithmic}
\REQUIRE $\TSV$，应该被加到$\CalG$和{\tt grad\_table}的变量。
\REQUIRE $\CalG$，要修改的图。
\REQUIRE $\CalG'$，根据参与梯度的节点$\CalG$的受限图。
\REQUIRE {\tt grad\_table}，将节点映射到对应梯度的数据结构。
\IF{$\SetV$ is in {\tt grad\_table} }
 \STATE Return ${\tt grad\_table}[\TSV]$
\ENDIF
\STATE $i \leftarrow 1$
\FOR{$\TSC$ in ${\tt get\_consumers}(\TSV, \CalG')$ }
\STATE ${\tt op} \leftarrow {\tt get\_operation}(\TSC)$
\STATE $\TSD \leftarrow {\tt build\_grad}(\TSC, \CalG, \CalG', {\tt grad\_table})$
\STATE $\TSG^{(i)} \leftarrow {\tt op.bprop}({\tt get\_inputs}(\TSC, \CalG'), \TSV, \TSD)$ 
\STATE $i \leftarrow i + 1$
\ENDFOR
\STATE $\TSG \leftarrow \sum_i \TSG^{(i)}$
\STATE ${\tt grad\_table}[\TSV] = \TSG$
\STATE 插入 $\TSG$ 和将其生成到$\CalG$中的操作
\STATE Return $\TSG$
\end{algorithmic}
\end{algorithm}


% -- 210 --
  
在\secref{sec:chain_rule_of_calculus}中，我们使用反向传播作为一种策略来避免多次计算链式法则中的相同子表达式。
由于这些重复子表达式的存在，简单的算法可能具有指数运行时间。
现在我们已经详细说明了反向传播算法，我们可以去理解它的计算成本。
如果我们假设每个操作的执行都有大致相同的开销，那么我们可以依据执行操作的数量来分析计算成本。
注意这里我们将一个操作记为计算图的基本单位，它实际可能包含许多算术运算（例如，我们可能将矩阵乘法视为单个操作）。
在具有$n$个节点的图中计算梯度，将永远不会执行超过$O(n^2)$个操作，或者存储超过$O(n^2)$个操作的输出。
这里我们是对计算图中的操作进行计数，而不是由底层硬件执行的单独操作，所以重要的是要记住每个操作的运行时间可能是高度可变的。
例如，两个矩阵相乘可能对应着图中的一个单独的操作，但这两个矩阵可能每个都包含数百万个元素。
我们可以看到，计算梯度至多需要$O(n^2)$的操作，因为在最坏的情况下，前向传播的步骤将在原始图的全部$n$个节点上运行（取决于我们想要计算的值，我们可能不需要执行整个图）。
反向传播算法在原始图的每条边添加一个~\gls{jacobian}~向量积，可以用$O(1)$个节点来表达。
因为计算图是有向无环图，它至多有$O(n^2)$条边。对于实践中常用图的类型，情况会更好。
大多数神经网络的代价函数大致是链式结构的，使得反向传播只有$O(n)$的成本。
这远远胜过简单的方法，简单方法可能需要在指数级的节点上运算。
这种潜在的指数级代价可以通过非递归地扩展和重写递归链式法则（\eqnref{eq:6.49}）来看出：
\begin{equation}
  \frac{\partial u^{(n)}}{\partial u^{(j)}} =
  \sum_{\substack{\text{path}(u^{(\pi_1)}, u^{(\pi_2)}, \ldots, u^{(\pi_t)}  ),\\ \text{from } \pi_1=j \text{ to }\pi_t = n}}
  \prod_{k=2}^t \frac{\partial u^{(\pi_k)}}{\partial u^{(\pi_{k-1})}}.
\end{equation}
由于节点$j$到节点$n$的路径数目可以关于这些路径的长度上指数地增长，所以上述求和符号中的项数（这些路径的数目），可能以前向传播图的深度的指数级增长。
会产生如此大的成本是因为对于$\frac{\partial u^{(i)}}{\partial u^{(j)}}$，相同的计算会重复进行很多次。
为了避免这种重新计算，我们可以将反向传播看作一种表填充算法，利用存储的中间结果$\frac{\partial u^{(n)}}{\partial u^{(i)}}$来对表进行填充。
图中的每个节点对应着表中的一个位置，这个位置存储对该节点的梯度。
通过顺序填充这些表的条目，反向传播算法避免了重复计算许多公共子表达式。
这种表填充策略有时被称为\firstgls{dynamic_programming}。
  
% -- 211 --
  
\subsection{实例：用于MLP训练的\glsentrytext{BP}}
\label{sec:example_back_propagation_for_mlp_training}

作为一个例子，我们利用反向传播算法来训练\gls{MLP}。

这里，我们考虑一个具有单个\gls{hidden_layer}的非常简单的\gls{MLP}。
为了训练这个模型，我们将使用\gls{minibatch}随机梯度下降算法。
反向传播算法用于计算单个\gls{minibatch}上的代价的梯度。
具体来说，我们使用训练集上的一\gls{minibatch}实例，将其规范化为一个设计矩阵$\MX$以及相关联的类标签向量$\Vy$。
网络计算隐藏特征层$\MH=\max\{0, \MX\MW^{(1)}\}$。
为了简化表示，我们在这个模型中不使用偏置。
假设我们的图语言包含~\verb|relu|操作，该操作可以对$\max\{0,\MZ\}$表达式的每个元素分别进行计算。
类的非归一化对数概率的预测将随后由$\MH\MW^{(2)}$ 给出。
假设我们的图语言包含~\verb|cross_entropy|操作，用以计算目标$\Vy$和由这些未归一化对数概率定义的概率分布间的\gls{cross_entropy}。
所得到的\gls{cross_entropy}定义了代价函数$J_\text{MLE}$。最小化这个\gls{cross_entropy}将执行对分类器的最大似然估计。
然而，为了使得这个例子更加真实，我们也包含一个正则项。
总的代价函数为
\begin{equation}
  J = J_{\text{MLE}} + \lambda \left ( \sum_{i, j} \left (W_{i, j}^{(1)} \right )^2 + \sum_{i, j} \left (W_{i, j}^{(2)} \right)^2 \right )
\end{equation}
包含了\gls{cross_entropy}和系数为$\lambda$的权重衰减项。
它的计算图在\figref{fig:chap6_mlp_example}中给出。
% fig 6.11
\begin{figure}[!htb]
\ifOpenSource
\centerline{\includegraphics{figure.pdf}}
\else
\centerline{\includegraphics{Chapter6/figures/mlp_example}}
\fi
\caption{用于计算代价函数的计算图，这个代价函数是使用\gls{cross_entropy}损失以及权重衰减训练我们的单层~\glssymbol{MLP}~示例所产生的。}
\label{fig:chap6_mlp_example}
\end{figure}

这个示例的梯度计算图实在太大，以致绘制或者阅读都将是乏味的。
这显示出了反向传播算法的优点之一，即它可以自动生成梯度，而这种计算对于软件工程师来说需要进行直观但冗长的手动推导。

我们可以通过观察\figref{fig:chap6_mlp_example}中的正向传播图来粗略地描述反向传播算法的行为。
为了训练，我们希望计算$\nabla_{\MW^{(1)}} J$和$\nabla_{\MW^{(2)}} J$。
有两种不同的路径从$J$后退到权重：一条通过\gls{cross_entropy}代价，另一条通过权重衰减代价。
权重衰减代价相对简单，它总是对$\MW^{(i)}$上的梯度贡献$2\lambda \MW^{(i)}$。
  
% -- 212 --
  
另一条通过\gls{cross_entropy}代价的路径稍微复杂一些。
令$\MG$是由~\verb|cross_entropy|操作提供的对未归一化对数概率$\MU^{(2)}$的梯度。
反向传播算法现在需要探索两个不同的分支。
在较短的分支上，它使用对矩阵乘法的第二个变量的反向传播规则，将$\MH^\top \MG$加到$\MW^{(2)}$的梯度上。
另一条更长些的路径沿着网络逐步下降。
首先，反向传播算法使用对矩阵乘法的第一个变量的反向传播规则，计算$\nabla_{\MH} J = \MG\MW^{(2)\top}$。
接下来，\verb|relu|操作使用其反向传播规则对先前梯度的部分位置清零，这些位置对应着$\MU^{(1)}$中所有小于0的元素。记上述结果为$\MG'$。 
反向传播算法的最后一步是使用对~\verb|matmul|操作的第二个变量的反向传播规则，将$\MX^\top \MG'$加到$\MW^{(1)}$的梯度上。

在计算了这些梯度以后，梯度下降算法或者其他优化算法所要做的就是使用这些梯度来更新参数。

对于MLP，计算成本主要来源于矩阵乘法。
在前向传播阶段，我们乘以每个权重矩阵，得到了$O(w)$数量的乘-加，其中$w$是权重的数量。
在反向传播阶段，我们乘以每个权重矩阵的转置，这具有相同的计算成本。
算法主要的存储成本是我们需要将输入存储到\gls{hidden_layer}的非线性中去。
这些值从被计算时开始存储，直到反向过程回到了同一点。
因此存储成本是$O(mn_h)$，其中$m$是\gls{minibatch}中样本的数目，$n_h$是\gls{hidden_unit}的数量。
  
% -- 213 --
  
\subsection{复杂化}
\label{sec:complications}

我们这里描述的反向传播算法要比实践中实际使用的实现要简单。

正如前面提到的，我们将操作的定义限制为返回单个张量的函数。
大多数软件实现需要支持可以返回多个张量的操作。 
例如，如果我们希望计算张量中的最大值和该值的索引，则最好在单次运算中计算两者，因此将该过程实现为具有两个输出的操作效率更高。

我们还没有描述如何控制反向传播的内存消耗。
反向传播经常涉及将许多张量加在一起。
在朴素方法中，将分别计算这些张量中的每一个，然后在第二步中对所有这些张量求和。 
朴素方法具有过高的存储瓶颈，可以通过保持一个缓冲器，并且在计算时将每个值加到该缓冲器中来避免该瓶颈。

反向传播的现实实现还需要处理各种数据类型，例如32位浮点数、64位浮点数和整型。
处理这些类型的策略需要特别的设计考虑。

一些操作具有未定义的梯度，重要的是跟踪这种情况并且确定用户请求的梯度是否是未定义的。

各种其他技术的特性使现实世界的微分更加复杂。 这些技术性并不是不可逾越的，本章已经描述了计算微分所需的关键知识工具，但重要的是要知道还有许多的精妙之处存在。
  
% -- 214 --
  
\subsection{深度学习界以外的微分}
\label{sec:differentiation_outside_the_deep_learning_community}

深度学习界在某种程度上已经与更广泛的计算机科学界隔离开来，并且在很大程度上发展了自己关于如何进行微分的文化态度。
更一般地，\firstgls{automatic_differentiation}领域关心如何以算法方式计算导数。 
这里描述的反向传播算法只是自动微分的一种方法。 
它是一种称为\firstgls{reverse_mode_accumulation}的更广泛类型的技术的特殊情况。 
其他方法以不同的顺序来计算链式法则的子表达式。 
一般来说，确定一种计算的顺序使得计算开销最小，是困难的问题。 
找到计算梯度的最优操作序列是NP完全问题\citep{naumann2008optimal}，在这种意义上，它可能需要将代数表达式简化为它们最廉价的形式。

例如，假设我们有变量$p_1,p_2\ldots,p_n$表示概率，以及变量$z_1,z_2,\ldots,z_n$表示未归一化的对数概率。
假设我们定义
\begin{equation}
  q_i = \frac{\exp(z_i)}{\sum_i \exp(z_i)},
\end{equation}
其中我们通过指数化、求和与除法运算构建softmax函数，并构造\gls{cross_entropy}损失函数$J=-\sum_i p_i\log q_i$。
人类数学家可以观察到$J$对$z_i$的导数有一个非常简单的形式：$q_i-p_i$
%\footnote{译者注：$\frac{\partial J}{\partial z_i} = -\frac{p_i}{q_i}\frac{\partial q_i}{\partial z_i} -\sum_{j \ne i} \frac{p_j}{q_j} \frac{\partial q_j}{\partial z_i}  = -\frac{p_i}{q_i}(q_i - q_i^2) - \sum_{j \ne i} \frac{p_j}{q_j}(-q_jq_i) =p_i(q_i - 1) + (1-p_i)q_i = q_i - p_i.$}。
反向传播算法不能够以这种方式来简化梯度，而是会通过原始图中的所有对数和指数操作显式地传播梯度。
一些软件库如Theano~\citep{bergstra+al:2010-scipy-small,Bastien-Theano-2012}能够执行某些种类的代数替换来改进由纯反向传播算法提出的图。
  
% -- 215 --
  
当前向图$\CalG$具有单个输出节点，并且每个偏导数$\frac{\partial u^{(i)}}{\partial u^{(j)}}$都可以用恒定的计算量来计算时，反向传播保证梯度计算的计算数目和前向计算的计算数目是同一个量级：这可以在\algref{alg:bprop}中看出，因为每个局部偏导数$\frac{\partial u^{(i)}}{\partial u^{(j)}}$以及递归链式公式（\eqnref{eq:6.49}）中相关的乘和加都只需计算一次。
因此，总的计算量是$O(\#\text{edges})$。
然而，可能通过对反向传播算法构建的计算图进行简化来减少这些计算量，并且这是NP完全问题。
诸如Theano和TensorFlow的实现使用基于匹配已知简化模式的试探法，以便重复地尝试去简化图。
我们定义反向传播仅用于计算标量输出的梯度，但是反向传播可以扩展到计算~\gls{jacobian_matrix}（该~\gls{jacobian_matrix}或者来源于图中的$k$个不同标量节点，或者来源于包含$k$个值的张量值节点）。
朴素的实现可能需要$k$倍的计算：对于原始前向图中的每个内部标量节点，朴素的实现计算$k$个梯度而不是单个梯度。
当图的输出数目大于输入的数目时，有时更偏向于使用另外一种形式的自动微分，称为\firstgls{forward_mode_accumulation}。
前向模式计算已经被提出用于循环神经网络梯度的实时计算，例如\citep{Williams89b}。
这也避免了存储整个图的值和梯度的需要，是计算效率和内存使用的折中。
前向模式和后向模式的关系类似于左乘和右乘一系列矩阵之间的关系，例如
\begin{equation}
  \MA \MB \MC \MD,
\end{equation}
其中的矩阵可以认为是~\gls{jacobian_matrix}。
例如，如果$\MD$是列向量，而$\MA$有很多行，那么这对应于一幅具有单个输出和多个输入的图，并且从最后开始乘，反向进行，只需要矩阵-向量的乘积。
这对应着反向模式。
相反，从左边开始乘将涉及一系列的矩阵-矩阵乘积，这使得总的计算变得更加昂贵。
然而，如果$\MA$的行数小于$D$的列数，则从左到右乘更为便宜，这对应着前向模式。

在机器学习以外的许多社区中，更常见的是使用传统的编程语言来直接实现微分软件，例如用Python或者C来编程，并且自动生成使用这些语言编写的不同函数的程序。
在深度学习界中，计算图通常使用由专用库创建的明确的数据结构表示。
专用方法的缺点是需要库开发人员为每个操作定义~\verb|bprop|方法，并且限制了库的用户仅使用定义好的那些操作。
然而，专用方法也允许定制每个操作的反向传播规则，允许开发者以非显而易见的方式提高速度或稳定性，对于这种方式自动的过程可能不能复制。

因此，反向传播不是计算梯度的唯一方式或最佳方式，但它是一个非常实用的方法，继续为深度学习社区服务。 
在未来，深度网络的微分技术可能会提高，因为深度学习的从业者更加懂得了更广泛的自动微分领域的进步。
  
% -- 216 --
  
\subsection{高阶微分}
\label{sec:higher_order_derivatives}

一些软件框架支持使用高阶导数。 
在深度学习软件框架中，这至少包括Theano和TensorFlow。
这些库使用一种数据结构来描述要被微分的原始函数，它们使用相同类型的数据结构来描述这个函数的导数表达式。
这意味着符号微分机制可以应用于导数（从而产生高阶导数）。

在深度学习的相关领域，很少会计算标量函数的单个二阶导数。
相反，我们通常对Hessian矩阵的性质比较感兴趣。
如果我们有函数$f:\SetR^n \to \SetR$，那么Hessian矩阵的大小是$n\times n$。
在典型的深度学习应用中，$n$将是模型的参数数量，可能很容易达到数十亿。
因此，完整的Hessian矩阵甚至不能表示。

典型的深度学习方法是使用\firstgls{Krylov_methods}，而不是显式地计算Hessian矩阵。
Krylov方法是用于执行各种操作的一组迭代技术，这些操作包括像近似求解矩阵的逆、或者近似矩阵的特征值或特征向量等，而不使用矩阵-向量乘法以外的任何操作。

为了在Hesssian矩阵上使用Krylov方法，我们只需要能够计算Hessian矩阵$\MH$和一个任意向量$\Vv$间的乘积即可。
实现这一目标的一种直观方法\citep{christianson1992automatic}是
\begin{equation}
  \MH \Vv=\nabla_{\Vx} \left [ (\nabla_{\Vx} f(x))^\top \Vv\right ].
\end{equation}
该表达式中两个梯度的计算都可以由适当的软件库自动完成。
注意，外部梯度表达式是内部梯度表达式的函数的梯度。

如果$\Vv$本身是由计算图产生的一个向量，那么重要的是指定自动微分软件不要对产生$\Vv$的图进行微分。

虽然计算Hessian通常是不可取的，但是可以使用Hessian向量积。
可以对所有的$i=1,\ldots,n$简单地计算$\MH \Ve^{(i)}$，其中$\Ve^{(i)}$是$e_i^{(i)}=1$并且其他元素都为0的\gls{one_hot}向量。

\section{历史小记}
\label{sec:historical_notes}

\gls{feedforward_network}可以被视为一种高效的非线性函数近似器，它以使用梯度下降来最小化函数近似误差为基础。
从这个角度来看，现代\gls{feedforward_network}是一般函数近似任务的几个世纪进步的结晶。
  
% -- 217 --
  
处于反向传播算法底层的链式法则是17世纪发明的\citep{Leibniz-1676,LHopital-1696}。
微积分和代数长期以来被用于求解优化问题的封闭形式，但梯度下降直到19世纪才作为优化问题的一种迭代近似的求解方法被引入\citep{cauchy1847}。

从20世纪40年代开始，这些函数近似技术被用于导出诸如感知机的机器学习模型。 
然而，最早的模型都是基于线性模型。 
来自包括Marvin Minsky的批评指出了线性模型族的几个缺陷，例如它无法学习XOR函数，这导致了对整个神经网络方法的抵制。

学习非线性函数需要\gls{MLP}的发展和计算该模型梯度的方法。
基于动态规划的链式法则的高效应用开始出现在20世纪60年代和70年代，主要用于控制领域\citep{Kelley-1960,Bryson-et-al-1961,Dreyfus-1962,Bryson-1969,Dreyfus-1973}，也用于灵敏度分析\citep{Linnainmaa-1976}。 
\cite{Werbos-1981}提出应用这些技术来训练人工神经网络。
这个想法以不同的方式被独立地重新发现后\citep{LeCun85,Parker-1985,Rumelhart86b-small}，最终在实践中得以发展。
\firstgls{parallel_distributed_processing}一书在其中一章提供了第一次成功使用反向传播的一些实验的结果\citep{Rumelhart86c}，这对反向传播的普及做出了巨大的贡献，并且开启了一个研究多层神经网络非常活跃的时期。
然而，该书作者提出的想法，特别是Rumelhart 和Hinton提出的想法远远超过了反向传播。
它们包括一些关键思想，关于可能通过计算实现认知和学习的几个核心方面，后来被冠以`` 联结主义''的名称，因为它强调了神经元之间的连接作为学习和记忆的轨迹的重要性。
特别地，这些想法包括分布式表示的概念\citep{Hinton-et-al-PDP1986}。

在反向传播的成功之后，神经网络研究获得了普及，并在20世纪90年代初达到高峰。 
随后，其他机器学习技术变得更受欢迎，直到2006年开始的现代深度学习复兴。

现代\gls{feedforward_network}的核心思想自20世纪80年代以来没有发生重大变化。
仍然使用相同的反向传播算法和相同的梯度下降方法。
1986年至2015年神经网络性能的大部分改进可归因于两个因素。
首先，较大的数据集减少了统计\gls{generalization}对神经网络的挑战的程度。
第二，神经网络由于更强大的计算机和更好的软件基础设施已经变得更大。
然而，少量算法上的变化也显著改善了神经网络的性能。
  
% -- 218 --
  
其中一个算法上的变化是用\gls{cross_entropy}族损失函数替代\gls{mean_squared_error}损失函数。
\gls{mean_squared_error}在20世纪80年代和90年代流行，但逐渐被\gls{cross_entropy}损失替代，并且最大似然原理的想法在统计学界和机器学习界之间广泛传播。
使用\gls{cross_entropy}损失大大提高了具有~\gls{sigmoid}~和softmax输出的模型的性能，而当使用\gls{mean_squared_error}损失时会存在饱和和学习缓慢的问题。

另一个显著改善\gls{feedforward_network}性能的算法上的主要变化是使用分段线性\gls{hidden_unit}来替代~\gls{sigmoid}~\gls{hidden_unit}，例如用\gls{ReLU}。
使用$\max\{0, z\}$函数的整流在早期神经网络中已经被引入，并且至少可以追溯到认知机（Cognitron）和神经认知机(Neocognitron)\citep{Fukushima75,Fukushima80}。
这些早期的模型没有使用\gls{ReLU}，而是将整流用于非线性函数。
尽管整流在早期很普及，在20世纪80年代，整流很大程度上被~\gls{sigmoid}~所取代，也许是因为当神经网络非常小时，\gls{sigmoid}~表现更好。
到21世纪初，由于有些迷信的观念，认为必须避免具有不可导点的激活函数，所以避免了\gls{ReLU}。
这在2009年开始发生改变。
\cite{Jarrett-ICCV2009-small}观察到，在神经网络结构设计的几个不同因素中``使用整流非线性是提高识别系统性能的最重要的唯一因素''。

对于小的数据集，\cite{Jarrett-ICCV2009-small}观察到，使用整流非线性甚至比学习\gls{hidden_layer}的权重值更加重要。
随机的权重足以通过\gls{rectifier_network}传播有用的信息，允许在顶部的分类器层学习如何将不同的特征向量映射到类标识。

当有更多数据可用时，学习开始提取足够的有用知识来超越随机选择参数的性能。
\cite{Glorot+al-AI-2011-small}说明，在深度\gls{rectifier_network}中的学习比在激活函数具有曲率或两侧饱和的深度网络中的学习更容易。

\gls{ReLU}还具有历史意义，因为它们表明神经科学继续对深度学习算法的发展产生影响。
\cite{Glorot+al-AI-2011-small}从生物学考虑\gls{ReLU}的导出。
半整流非线性旨在描述生物神经元的这些性质：(1) 对于某些输入，生物神经元是完全不活跃的。
(2) 对于某些输入，生物神经元的输出和它的输入成比例。
(3) 大多数时间，生物神经元是在它们不活跃的状态下进行操作（即它们应该具有\firstgls{sparse_activation}）。
  
% -- 219 --
  
当2006年深度学习开始现代复兴时，\gls{feedforward_network}仍然有不良的声誉。
从2006年至2012年，人们普遍认为，\gls{feedforward_network}不会表现良好，除非它们得到其他模型的辅助，例如概率模型。
现在已经知道，只要具备适当的资源和工程实践，\gls{feedforward_network}表现得非常好。
今天，\gls{feedforward_network}中基于梯度的学习被用作发展概率模型的工具，例如\chapref{chap:deep_generative_models}中描述的\gls{VAE}和生成式对抗网络。
\gls{feedforward_network}中基于梯度的学习自2012年以来一直被视为一种强大的技术，并应用于许多其他机器学习任务，而不是被视为必须由其他技术支持的不可靠技术。
在2006 年，业内使用无监督学习来支持监督学习，现在更讽刺的是，更常见的是使用监督学习来支持无监督学习。

\gls{feedforward_network}还有许多未实现的潜力。
未来，我们期望它们用于更多的任务，优化算法和模型设计的进步将进一步提高它们的性能。
本章主要描述了神经网络模型族。
在接下来的章节中，我们将讨论如何使用这些模型——如何对它们进行正则化和训练。

  
% -- 220 --