卷积神经网络——池化层学习——最大池化

池化层（Pooling layers）

除了卷积层，卷积网络也经常使用池化层来缩减模型的大小，提高计算速度，同时提高所提取特征的鲁棒性，我们来看一下。

先举一个池化层的例子，然后我们再讨论池化层的必要性。

假如输入是一个4×4矩阵，用到的池化类型是最大池化（max pooling）。

执行最大池化的树池是一个2×2矩阵。

执行过程非常简单，把4×4的输入拆分成不同的区域，我把这个区域用不同颜色来标记。

对于2×2的输出，输出的每个元素都是其对应颜色区域中的最大元素值。
左上区域的最大值是9，右上区域的最大元素值是2，左下区域的最大值是6，右下区域的最大值是3。

为了计算出右侧这4个元素值，我们需要对输入矩阵的2×2区域做最大值运算。

这就像是应用了一个规模为2的过滤器，因为我们选用的是2×2区域，步幅是2，这些就是最大池化的超参数。

因为我们使用的过滤器为2×2，最后输出是9。

然后向右移动2个步幅，计算出最大值2。

然后是第二行，向下移动2步得到最大值6。

最后向右移动3步，得到最大值3。

这是一个2×2矩阵，即f=2，步幅是2，即s=2。

这是对最大池化功能的直观理解，你可以把这个4×4区域看作是某些特征的集合，也就是神经网络中某一层的非激活值集合。

---

---

数字大意味着可能探测到了某些特定的特征，左上象限具有的特征可能是一个垂直边缘，一只眼睛，或是大家害怕遇到的CAP特征。

显然左上象限中存在这个特征，这个特征可能是一只猫眼探测器。

然而，右上象限并不存在这个特征。

最大化操作的功能就是只要在任何一个象限内提取到某个特征，它都会保留在最大化的池化输出里。

所以最大化运算的实际作用就是，如果在过滤器中提取到某个特征，那么保留其最大值。

如果没有提取到这个特征，可能在右上象限中不存在这个特征，那么其中的最大值也还是很小，这就是最大池化的直观理解。

必须承认，人们使用最大池化的主要原因是此方法在很多实验中效果都很好。

尽管刚刚描述的直观理解经常被引用，不知大家是否完全理解它的真正原因，不知大家是否理解最大池化效率很高的真正原因。

其中一个有意思的特点就是，它有一组超参数，但并没有参数需要学习。

实际上，梯度下降没有什么可学的，一旦确定了f和s，它就是一个固定运算，梯度下降无需改变任何值。

---

---

我们来看一个有若干个超级参数的示例，输入是一个5×5的矩阵。

我们采用最大池化法，它的过滤器参数为3×3，即f=3，步幅为1，s=1，输出矩阵是3×3。

之前讲的计算卷积层输出大小的公式同样适用于最大池化，
即(n+2p-f)/s+1，这个公式也可以计算最大池化的输出大小。

此例是计算3×3输出的每个元素，我们看左上角这些元素，注意这是一个3×3区域，因为有3个过滤器，取最大值9。

然后移动一个元素，因为步幅是1，蓝色区域的最大值是9。

继续向右移动，蓝色区域的最大值是5。

然后移到下一行，因为步幅是1，我们只向下移动一个格，所以该区域的最大值是9。

这个区域也是9。

这两个区域的最大值都是5。

最后这三个区域的最大值分别为8，6和9。

超参数f=3，s=1，最终输出如图所示。

以上就是一个二维输入的最大池化的演示，如果输入是三维的，那么输出也是三维的。

例如，输入是5×5×2，那么输出是3×3×2。

计算最大池化的方法就是分别对每个通道执行刚刚的计算过程。

如上图所示，第一个通道依然保持不变。

对于第二个通道，我刚才画在下面的，在这个层做同样的计算，得到第二个通道的输出。

一般来说，如果输入是5×5×n_c，输出就是3×3×n_c，n_c个通道中每个通道都单独执行最大池化计算，以上就是最大池化算法。

func MaxPooling(data [][]float64, pool_size int, step int)([][]float64){
    result := [][]float64{}
    for i := 0; i < len(data) - pool_size + 1; i += step{
        temp := []float64{}
        for j := 0; j < len(data[0]) - pool_size + 1; j += step{
            var max_num float64 = 0
            for cur_i := i; cur_i < i + pool_size; cur_i++{
                for cur_j :=j; cur_j < j + pool_size; cur_j++{
                    max_num = math.Max(max_num, data[cur_i][cur_j])
                }
            }
            temp = append(temp, max_num)
        }
        result = append(result, temp)
    }
    return result
}