非监督特征学习与深度学习（十二）---- 卷积神经网络

注：本文转载自https://github.com/ysh329/Chinese-UFLDL-Tutorial
因为github上的makedown格式显示的不够完全，看的非常不方便，因此放到CSDN上比较好查阅学习。

卷积神经网络（Convolutional Nerual Network）

概述（Overview）

卷积神经网络（ CNN ）是有一个或多个卷积层（常伴有下采样步骤）并后面跟一个或多个全连接层的标准多层神经网络。卷积神经网络在体系结构的设计利用了输入图像的二维结构（其它的二维输入还有语音信号等）。

卷积神经网络的实现是借助局部连接和在其之后的绑定权重，其中池化操作有平移不变特征。卷积神经网络的另一个优点在于更容易训练，并且卷积神经网络的参数虽多但却比相同隐含单元的全连接网络要少。

在本节内容中，我们将会讨论卷积神经网络的架构和用来计算模型中不同参数梯度的反向传播算法。卷积和池化更详细的具体操作见本教程的各自章节。

网络架构（Architecture）

一个卷积神经网络是由很多个卷积层、后接可选的（一层或多层）下采样层以及后接的全连接层组成的。卷积层的输入是一个规模为 m×m×r 的图像，其中（第一个） m 是图像的高度，第二个 m 是宽度， r 是图像的通道个数。，例如一个 RGB 图像的通道数 r=3 。卷积层有 k 个规模为 n×n×q 的滤波器（或称为核），其中 n 小于图片的维度（译者注：图片的高度或宽度）， q 既可以是通道个数 r ，也可能对于不同的滤波器（或称为核）而不同。

过滤器的规模增加了局部连接的结构（译者注：不是很理解这句话的深层含义：The size of the filters gives rise to the locally connected structure.），原图被（这种结构）卷积成为规模为 m−n+1 的 k 个特征图（译者注：这里的“特征图”，即 feature map ，是二维的。其中 m−n+1 是一个特征图的宽度或者高度，原图是正方形，所以这里特征图的边长为 m−n+1 ）。

之后，每个（特征）图通常在 p×p 的相邻区域进行平均值或最大值下采样（译者注： CNN 中的下采样即池化），对于相邻区域 p 的值范围从 2 （对于小图片值为 2 ，例如 MNIST 手写图片数据集）并通常不超过 5 （对于大图片）。

在下采样层的之前或之后，（译者注：会对结果）应用一个附加的偏置和 S 型的非线性（译者注：的映射）。

下图描述了卷积神经网络中一个完整的卷积和下采样层。一样颜色的（译者注：神经）单元有着连接权重。

图1：一个卷积神经网络的第一层的池化过程。
相同颜色的神经单元有连接权重，不同颜色的神经元表示不同的滤波器（图）
（译者注：这里不同颜色的神经元代表着不同的滤波器或核。
对于“RF”，个人理解为感受野（ Receptive Field ）的缩写，
感受野是视觉系统信息处理的基本结构和功能单元。
注意：滤波器和核是一个意思，滤波器和通道不是一个意思）。

在卷积层后可能有任意个全连接层。在一个标准的多层神经网络中，被密集连接（译者注：即全连接）的这些层是一样的。

反向传播（Back Propagation）

对于网络中的成本函数 J(W,b;x,y) ， δ(l+1) 是第 (l+1) 层的误差项，其中 (W,b) 是（译者注：成本函数的）参数， (x,y) 是带有标签的训练数据。如果第 l 层与第 (l+1) 层是密集连接（译者注：即全连接）的，那么可计算第 l 层的误差项为：

δ(l)=((W(l))Tδ(l+1))∙f′(z(l))

计算第 l 层梯度为：

∇W(l)J(W,b;x,y)∇b(l)J(W,b;x,y)=δ(l+1)(a(l))T,=δ(l+1).

如果第 l 层是一个卷积和下采样（译者注：即池化）层，那么误差被传播（的过程，可通过如下公式）表示为：

δ(l)k=upsample((W(l)k)Tδ(l+1)k)∙f′(z(l)k)

其中， k 是滤波器（译者注：或称为核，同义）的索引编号， f′(z(l)k) 是激活函数的导数。上采样操作的误差传播是通过池化层（译者注：即下采样层）与进到池化层的各神经元（译者注：“进到池化层的各神经元”，即池化层前一层的神经元）计算误差。

举个例子，若做完了平均值池化，然后对（前一层的）池化（神经）单元进行简单的均匀分布（译者注：均匀分布，即 P(X=k)=1m,k=1,...,m ）上采样。在最大值池化中，即使（上层）输入（译者注：在池化的相邻区域内）的变化很小，（池化）单元会对结果产生扰动（译者注：本句翻译不确定，”In max pooling the unit which was chosen as the max receives all the error since very small changes in input would perturb the result only through that unit. “）。

Finally, to calculate the gradient w.r.t to the filter maps, we rely on the border handling convolution operation again and flip the error matrix δ(l)k the same way we flip the filters in the convolutional layer.

最后，为了计算特征图（即 feature map ）的梯度，我们再次借助边界来处理卷积运算，并翻转误差矩阵 δ(l)k （译者注：这里的“翻转”，即逆时针旋转 180° 。逆时针旋转 90° ，对应 MATLAB 中的函数 rot90 ， rot90(δ(l+1)k,2) 表示对二维矩阵 δ(l+1)k 逆时针旋转 90° 2 次），这个过程与在卷基层翻转（译者注：即逆时针旋转 180° ）滤波器是一样的。

∇W(l)kJ(W,b;x,y)∇b(l)kJ(W,b;x,y)=∑i=1m(a(l)i)∗rot90(δ(l+1)k,2),=∑a,b(δ(l+1)k)a,b.

其中， a(l) 是第 l 层的输入，并且有 a(1) 是输入的图像。操作 (a(l)i)∗δ(l+1)k 是第 l 层的第 i 个输入关于第 k 个滤波器（或称为核）的“有效的”卷积（操作，译者注： δ(l+1)k 是误差矩阵）。