机器学习（十三）卷积神经网络(CNN)详细公式推导

摘要

CNN 是为了识别二维形状而特殊设计的多层感知器, 对二维形状的缩放, 倾斜或其它形式的变形具有高度不变性. 每一个神经元从上一层的局部区域得到输入, 这迫使神经元提取局部特征. 一旦一个特征被提取出来, 它相对于其它特征的位置被近似保留下来, 而忽略掉精确的位置. 每个卷积层后面跟着一个池化, 使得特征图的分辨率降低, 而有利于降低对二维图形的平移或其他形式的敏感度.

权值共享是 CNN 的一大特点, 每一个卷积核滑遍整个输入矩阵, 而参数的数目等于卷积核矩阵元素的数目加上一个偏置, 相对于 BP 网络大大减少了训练的参数. 例如对于一张 300 × 300 的彩色图片, 如果第一层网络有25个神经元, 使用 BP网络, 参数的数目可达 675万个( 300 × 300 × 3 × 25 + 1), 而卷积网络只有76个( 25 × 3 + 1 )

Abstract

CNN is a multilayer perceptron specially designed to recognize 2D shapes and is highly invariant to scaling, tilting or other forms of deformation. Each neuron gets its input from the local region of the previous layer, which forces the neuron to extract local features. Once a feature is extracted, its position relative to other features is approximately preserved, while the exact position is ignored. Each convolutional layer is followed by a pooling, which reduces the resolution of the feature map in favor of reducing sensitivity to translation or other forms of 2D graphics.

Weight sharing is a major feature of CNN. Each convolution kernel slides through the entire input matrix, and the number of parameters is equal to the number of convolution kernel matrix elements plus a bias, which greatly reduces the training parameters compared with BP network. For example, for a 300×300 300\times300 300×300 color image, if the first layer network has 25 neurons, using BP network, the number of parameters can reach 6.75 million (300×300 ×3 × 25 + 1). Convolutional networks have only 76 (25 × 3 + 1)

1、 CNN 结构介绍

上面是一个简单的 CNN 结构图, 第一层输入图片, 进行卷积(Convolution)操作, 得到第二层深度为 3 的特征图(Feature Map). 对第二层的特征图进行池化(Pooling)操作, 得到第三层深度为 3 的特征图. 重复上述操作得到第五层深度为 5 的特征图, 最后将这 5 个特征图, 也就是 5 个矩阵, 按行展开连接成向量, 传入全连接(Fully Connected)层, 全连接层就是一个 BP 神经网络. 图中的每个特征图都可以看成是排列成矩阵形式的神经元, 与 BP神经网络中的神经元大同小异. 下面是卷积和池化的计算过程.

卷积

对于一张输入图片, 将其转化为矩阵, 矩阵的元素为对应的像素值. 假设有一个 5 × 5 的图像，使用一个 3 × 3 的卷积核进行卷积，可得到一个 3 × 3 的特征图. 卷积核也称为滤波器(Filter).

具体的操作过程如下图所示:
黄色的区域表示卷积核在输入矩阵中滑动, 每滑动到一个位置, 将对应数字相乘并求和, 得到一个特征图矩阵的元素. 注意到, 动图中卷积核每次滑动了一个单位, 实际上滑动的幅度可以根据需要进行调整. 如果滑动步幅大于 1, 则卷积核有可能无法恰好滑到边缘, 针对这种情况, 可在矩阵最外层补零, 补一层零后的矩阵如下图所示:

可根据需要设定补零的层数. 补零层称为 Zero Padding, 是一个可以设置的超参数, 但要根据卷积核的大小, 步幅, 输入矩阵的大小进行调整, 以使得卷积核恰好滑动到边缘.

一般情况下, 输入的图片矩阵以及后面的卷积核, 特征图矩阵都是方阵, 这里设输入矩阵大小为 w, 卷积核大小为 k, 步幅为 s, 补零层数为 p , 则卷积后产生的特征图大小计算公式为:
$w^{\prime }=\frac{(w+2p-k)}{s}+1$
上图是对一个特征图采用一个卷积核卷积的过程, 为了提取更多的特征, 可以采用多个卷积核分别进行卷积, 这样便可以得到多个特征图. 有时, 对于一张三通道彩色图片, 或者如第三层特征图所示, 输入的是一组矩阵, 这时卷积核也不再是一层的, 而要变成相应的深度.
上图中, 最左边是输入的特征图矩阵, 深度为 3, 补零(Zero Padding)层数为 1, 每次滑动的步幅为 2. 中间两列粉色的矩阵分别是两组卷积核, 一组有三个, 三个矩阵分别对应着卷积左侧三个输入矩阵, 每一次滑动卷积会得到三个数, 这三个数的和作为卷积的输出. 最右侧两个绿色的矩阵分别是两组卷积核得到的特征图.

池化

池化又叫下采样(Dwon sampling), 与之相对的是上采样(Up sampling). 卷积得到的特征图一般需要一个池化层以降低数据量. 池化的操作如下图所示:
和卷积一样, 池化也有一个滑动的核, 可以称之为滑动窗口, 上图中滑动窗口的大小为 2 × 2, 步幅为 2 , 每滑动到一个区域, 则取最大值作为输出, 这样的操作称为 Max Pooling. 还可以采用输出均值的方式, 称为 Mean Pooling.

全连接

经过若干层的卷积, 池化操作后, 将得到的特征图依次按行展开, 连接成向量, 输入全连接网络

2、公式

输入
$V=\text{conv2}(W,X,\text{"valid"})+b$
输出
$Y=\phi (V)$
上面的输入输出公式是对每一个卷积层而言的, 每一个卷积层都有一个不同的权重矩阵 W , 并且 W , X , Y 是矩阵形式.对于最后一层全连接层, 设为第 L层, 输出是向量形式的$y^L$, 期望输出是 d , 则有总误差公式
总误差
$E=\frac{1}{2}||d-y^L|{|}^2$
conv2() 是 卷积运算的卷积函数, 第三个参数 valid 指明卷积运算的类型, 前面介绍的卷积方式就是 valid 型. W是卷积核矩阵, X是输入矩阵, b是偏置, φ ( x ) 是激活函数. 总误差中的 d ， y分别是期望输出和网络输出的向量. 全连接层神经元的输入输出计算公式与 BP网络完全相同.

3、 梯度下降和反向传播算法

与 BP神经网络一样, CNN 也是通过梯度下降和反向传播算法进行训练的, 则全连接层的梯度公式与 BP网络完全一样, 这里就不重复了. 下面介绍卷积层和池化层的梯度公式.
$\frac{\partial E}{\partial w_{ij}}=\frac{\partial E}{\partial v_{ij}}\frac{\partial v_{ij}}{\partial w_{ij}}={\delta }_{ij}\frac{\partial v_{ij}}{\partial }$
现在的目标是求出${\delta }_{ij}和\frac{\partial v_{ij}}{\partial w_{ij}}$​​的表达式. 首先求${\delta }_{ij}$​, 从一个简单情况开始.
设输入是 3 × 3的矩阵, 卷积核是 2 × 2 的矩阵.
(1) 当步长为 1 时, 则可得到 2 × 2 的特征图矩阵. 如下图所示:

绿色矩阵是第 l-1层神经元, 红色的是第 l 层, 黄色的是第 l层的卷积核, 图中省略掉了第 l − 1层的卷积核. 每个矩阵的元素可以看成一个神经元, 于是和 BP网络一样, 可以对每个神经元定义局部梯度, 为了看着清楚, 把矩阵的元素换成$\delta$. 神经网络中数值的传递也可以看成导数的传递. 假设已经求出了 ${\delta }^l,则{\delta }^{l-1}$ 的计算公式为
${\delta }^{l-1}=\text{conv2}(\text{rot180}(W^l),{\delta }^l,\text{’full’}){\phi }^{\prime }(v^{l-1})$
这是一个递归公式. $\delta$, w , v 是矩阵形式, 公式中 rot180() 函数表示将矩阵逆时针旋转 180 度, 参数 full 表示进行的完全的卷积运算, 这种运算形象的解释如下:

首先将特征图矩阵补一圈零, 然后用权重矩阵对补零后的特征图矩阵进行 vaild 型卷积运算. 这里补零的层数等于卷积核矩阵的大小减去 1, full 型卷积的内涵在于使卷积从边缘开始, 即先计算 $w_{11}^l{\delta }_{11}^l$​, 然后是 $w_{12}^l{\delta }_{11}^l+w_{11}^l{\delta }_{12}^l$, 以此类推.
(2) 当滑动步幅大于 1 时

上图分别是步长为 1, 2 时的卷积过程, 可以发现步长为 1 是一个比较完全的卷积, 当步长增大时, 会跳过某些区域, 因此可以通过对特征图补零的方式将步幅为 2 时生成的特征图转换为步幅 1 的特征图, 然后直接利用前面推导的公式即可求出${\delta }^{l-1}$. 补零如下图所示

(3) 当卷积核数目为 N时, 与 BP 网络隐藏层神经元局部梯度计算类似, 每一个权重的改变, 都会影响到下一层每一个神经元输入的改变, 所以误差的反响传递需要把误差之和传到前一层, 则公式可简单地修改为
${\delta }^{l-1}={\sum }_{i=0}^N\text{conv2}(\text{rot180}(W_i^l),{\delta }_i^l,\text{’full’}){\phi }^{\prime }(v^{l-1})$
(4) 对于池化层的梯度传递, 情况就简单了一些, 下面以 Max Pooling 为例
图中滑动窗口大小为 2 × 2, 步幅为 2, 则得到一个 2 × 2的特征图, Max Pooling 方法是选择最大的数值作为输出, 由于过程中没有任何运算, 不存在梯度的变化, 所以误差可以直接由第 l层传递到第 l − 1层, 如下图所示

图中假设四个角上的元素值最大, 则对应位置的局部梯度可直接反向传递, 其它位置由于没有连接, 则梯度设为 0.

如果池化类型是 Mean Pooling, 即取滑动窗口内的均值作为输出, 则将 l层的梯度的$\frac{1}{4}$ 倍反向传递到 l − 1层, 如下图所示

到此为止已经推导出了不同情况下的$\delta$的递归公式, 下面推导 $\frac{\partial v_{ij}}{\partial w_{ij}}$的表达式. 下图所示, 为了表达清楚, 将符号做一些替换, 相对于第 l层, 第 l − 1层的输出变量符号可由 y 变为 x.

根据卷积的原理, 可以得到
$\begin{array}{rl}{y}_{11}^l& =w_{11}^l{x}_{11}^{l-1}+w_{12}^l{x}_{12}^{l-1}+w_{21}^l{x}_{21}^{l-1}+w_{22}^l{x}_{22}^{l-1}\\ y_{12}^l& =w_{11}^l{x}_{12}^{l-1}+w_{12}^l{x}_{13}^{l-1}+w_{21}^l{x}_{22}^{l-1}+w_{22}^l{x}_{23}^{l-1}\\ y_{21}^l& =w_{11}^l{x}_{21}^{l-1}+w_{12}^l{x}_{22}^{l-1}+w_{21}^l{x}_{31}^{l-1}+w_{22}^l{x}_{32}^{l-1}\\ y_{22}^l& =w_{11}^l{x}_{22}^{l-1}+w_{12}^l{x}_{23}^{l-1}+w_{21}^l{x}_{32}^{l-1}+w_{22}^l{x}_{33}^{l-1}\end{array}$
可以发现, 由于共享权重, 一个权重$w_{ij}$的改变, 会影响所有的 y. 因而对某一个权重求梯度时, 需要运用全微分公式. 现在以 $\frac{\partial E}{\partial w_{11}^l}$ 为例, 求得表达式
$\begin{array}{l}\frac{\partial E}{\partial w_{11}^l}=\frac{\partial E}{\partial y_{11}^l}\frac{\partial y_{11}^l}{\partial w_{11}^l}+\frac{\partial E}{\partial y_{12}^l}\frac{\partial y_{12}^l}{\partial w_{12}^l}+\frac{\partial E}{\partial y_{21}^l}\frac{\partial y_{21}^l}{\partial w_{21}^l}+\frac{\partial E}{\partial y_{22}^l}\frac{\partial y_{22}^l}{\partial w_{22}^l}\\ ={\delta }_{11}^l{x}_{11}^{l-1}+{\delta }_{12}^l{x}_{12}^{l-1}+{\delta }_{21}^l{x}_{21}^{l-1}+{\delta }_{22}^l{x}_{22}^{l-1}\end{array}$
同理可得
$\begin{array}{l}\frac{\partial E}{\partial w_{12}^{\prime }}={\delta }_{11}^l{x}_{12}^{l-1}+{\delta }_{12}^l{x}_{13}^{l-1}+{\delta }_{21}^l{x}_{22}^{l-1}+{\delta }_{22}^l{x}_{23}^{l-1}\\ \frac{\partial E}{\partial w_{21}^l}={\delta }_{11}^l{x}_{21}^{l-1}+{\delta }_{12}^l{x}_{22}^{l-1}+{\delta }_{21}^l{x}_{31}^{l-1}+{\delta }_{22}^l{x}_{32}^{l-1}\\ \frac{\partial E}{\partial u_2^l}={\delta }_{11}^l{x}_{22}^{l-1}+{\delta }_{12}^l{x}_{23}^{l-1}+{\delta }_{21}^l{x}_{32}^{l-1}+{\delta }_{22}^l{x}_{33}^{l-1}\end{array}$
根据得到的结果, 可以发现计算 $\frac{\partial E}{\partial w_{ij}^l}$ 的规律:
首先将前面的图中卷积核和特征图调换一下位置, 将特征图矩阵的元素替换为$\delta$, 将卷积核矩阵的元素替换为 E对 w的偏导. 根据前面得到的偏导公式, 可以发现特征图对输入矩阵卷积可得到右边黄色的偏导矩阵. 即

$\frac{\partial E}{\partial w^l}=\text{conv2}({\delta }^l,x^{l-1},\text{’valid’})$
公式中 w, $\delta$, x是以矩阵形式表达的. 同时, 可得到
$\begin{array}{l}\frac{\partial E}{\partial b^l}=\frac{\partial E}{\partial v_{11}^l}\frac{\partial v_{11}^l}{\partial b^l}+\frac{\partial E}{\partial v_{12}^l}\frac{\partial v_{12}^l}{\partial b^l}+\frac{\partial E}{\partial v_{21}^l}\frac{\partial v_{21}^l}{\partial b^l}+\frac{\partial E}{\partial v_{22}^l}\frac{\partial v_{22}^l}{\partial b^l}\\ ={\delta }_{11}^l+{\delta }_{12}^l+{\delta }_{21}^l+{\delta }_{22}^l\\ ={\sum }_i{\sum }_j{\delta }_{ij}^l\end{array}$

$\delta$ 递归公式的推导

CNN 是一种前馈型神经网络, 每一个神经元只与前一层神经元相连, 接受前一层的输出, 经过运算后输出给下一层, 各层之间没有反馈. 因此对于计算局部梯度, 可以根据后一层神经元的梯度, 递归地向前传递. 首先根据前面演示过的一张图

局部梯度的定义公式
${\delta }_{ij}^{l-1}=\frac{\partial E}{\partial v_{ij}^{l-1}}=\frac{\partial E}{\partial y_{ij}^{l-1}}\frac{\partial y_{ij}^{l-1}}{\partial v_{ij}^{l-1}}=\frac{\partial E}{\partial y_{ij}^{l-1}}\phi \left(v_{ij}^{l-1}\right)$
现在我们求出几项 $\frac{\partial E}{\partial v_{ij}^{l-1}}$​ 来寻找规律, 首先假设图中第 l层中的 δ 已知, 根据卷积运算的原理可得
$\begin{array}{l}{v}_{11}^l=w_{11}^l{x}_{11}^{l-1}+w_{12}^l{x}_{12}^{l-1}+w_{21}^l{x}_{21}^{l-1}+w_{22}^l{x}_{22}^{l-1}+b\\ \frac{\partial E}{\partial y_{11}^l}=\frac{\partial E}{\partial v_{11}^l}\frac{\partial v_{11}^l}{\partial y_{11}^l}={\delta }_{11}^l{w}_{11}^l\end{array}$
同理可得
$\begin{array}{l}\frac{\partial E}{\partial y_{12}^1}={\delta }_{11}^l{w}_{12}^l+{\delta }_{12}^l{w}_{11}^l\\ \frac{\partial E}{\partial y_{22}^l}={\delta }_{11}^l{w}_{22}^l+{\delta }_{12}^l{w}_{21}^l+{\delta }_{21}^l{w}_{12}^l+{\delta }_{22}^l{w}_{11}^l\end{array}$
则可以发现 $\frac{\partial E}{\partial y_{ij}^{l-1}}$ 的计算公式可以表达为上图中黄色卷积核在补零后的特征图卷积, 步幅为 1.
$\frac{\partial E}{\partial y^{l-1}}=\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{v}2\left(\mathrm{rot}180\left(w^l\right),{\delta }^l,\text{ ’full’ }\right)$
代入局部梯度公式得

${\delta }^{l-1}=\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{v}2\left(\mathrm{rot}180\left(w^l\right),{\delta }^l,\text{ ’full’ }\right){\phi }^{\prime }(v^{l-1})$