深度学习中的身份很重要

莫里茨·哈特 谷歌大脑。 m@mrtz.org    马腾宇 斯坦福大学。 tengyu@stanford.edu 在 Google 进行的工作。
摘要

深度学习中的一个新兴设计原则是深度人工神经网络的每一层都应该能够轻松地表达恒等变换。 这个想法不仅激发了各种标准化技术,例如批量标准化,而且也是残差网络取得巨大成功的关键。

在这项工作中,我们将身份参数化的原理置于更坚实的理论基础上,同时取得了进一步的实证进展。 我们首先给出一个非常简单的证明,证明任意深度的线性残差网络没有虚假的局部最优值。 线性前馈网络在其标准参数化中的相同结果要微妙得多。 其次,我们表明,具有 ReLu 激活的残差网络具有通用的有限样本表达能力,只要模型的参数多于样本大小,网络就可以表示其样本的任何函数。

直接受到我们理论的启发,我们尝试了一种极其简单的残差架构,仅包含残差卷积层和 ReLu 激活,但没有批量归一化、dropout 或最大池。 我们的模型在 CIFAR10、CIFAR100 和 ImageNet 分类基准上比之前的全卷积网络有了显着改进。

1简介

用于图像分类的传统卷积神经网络,例如 AlexNet ([13]),其参数化方式为:当所有可训练权重均为 0 时,卷积层表示 0-映射。 此外,权重围绕 0. 对称初始化 这种标准参数化使得使用随机梯度方法训练的卷积层保留已经很好的特征变得非常重要。 换句话说,这样的卷积层在训练时不能轻易收敛到恒等变换。

[9] 通过批量归一化观察到并部分解决了这个缺点,即使用学习的均值和协方差对输入进行逐层白化。 但这个想法仍然有些隐含,直到残差网络([6];[7])明确引入了卷积层的重新参数化,这样当所有可训练的权重都是0,层代表恒等函数。 形式上,对于输入 x,,每个残差层的形式为 x+h(x), 而不是 h(x). 这种简单的重新参数化允许更深层次的架构,很大程度上避免了梯度消失(或爆炸)的问题。 从那时起,残差网络以及使用相同参数化的后续架构在各种计算机视觉基准(例如 CIFAR10 和 ImageNet)上始终取得了最先进的结果。

1.1 我们的贡献

在这项工作中,我们从理论角度考虑恒等参数化,同时将我们的一些理论见解转化为实验。 宽松地说,我们的第一个结果强调了身份参数化如何使优化变得更容易,而我们的第二个结果表明对于表示来说也是如此。

线性残差网络。

由于一般的非线性神经网络超出了当前优化理论方法的范围,因此我们将深度线性网络的情况视为简化模型。 线性网络将任意线性映射表示为矩阵序列 AA2A1. 目标函数为 𝔼yAA1x2,其中 y=Rx 用于某些未知的线性变换,Rx 是从分布中得出的。 近年来,人们对这种线性网络进行了积极的研究,作为通向一般非线性情况的垫脚石(参见1.2节)。 尽管 AA1 只是一个线性映射,但因子变量 (A,,A1) 的优化问题是非凸的。

与残差网络类似,我们将目标函数参数化为

minA1,,A𝔼y(I+A)(I+A1)x2. (1.1)

直观地说,当深度 足够大时,我们可以希望目标函数 R 具有因子表示,其中每个矩阵 Ai 具有较小的矩阵规范。 例如,任何对称正半inite 矩阵 O 都可以写成乘积 O=OO1,,其中每个 Oi=O1/ 对于大的 都非常接近同一性,因此 Ai=OiI 具有较小的谱规范。 我们首先证明类似的主张对于所有具有正行列式的线性变换 R 都是正确的111正如下面将要讨论的定理2.1,不失一般性地假设R的行列式是积极的。 . 具体来说,我们证明对于每个具有 det(R)>0 的线性变换 R,都存在一个 (1.1) 的全局优化器 (A1,,A),使得足够大的深度,

max1iAiO(1/). (1.2)

这里,A表示A.的谱范数 常数因子取决于R.的条件。我们在定理2.1中给出了正式的陈述。 该定理有一个有趣的结果,即随着深度的增加,存在较小的范数解,因此正则化可能会抵消参数的增加。

确定了小范数解的存在性后,我们对线性残差网络的主要结果表明,当所有矩阵都具有足够小的范数时,目标函数 (1.1) 实际上很容易优化。 更正式地说,让 A=(A1,,A)f(A) 表示 (1.1) 中的目标函数,我们可以证明仅当 f(A)=0 时梯度才消失> 前提是maxiAiO(1/). 参见定理2.2 这个结果意味着线性残差网络除了全局最优之外没有临界点 相比之下,对于标准线性神经网络,通过[12]的工作,我们只知道这些网络除了全局最优之外没有局部最优,但并不排除其他临界点。 事实上,设置 Ai=0 总是会导致标准参数化中出现错误的临界点。

通用有限样本表达能力。

回到具有 ReLU 激活的非线性残差网络,我们可以问:仅基于具有 ReLU 激活的残差层的深度神经网络的表达能力如何? 为了回答这个问题,我们给出了一个非常简单的结构,表明这种残差网络具有完美的有限样本表达能力。 换句话说,具有 ReLU 激活的残差网络可以轻松表达大小为 n, 的样本的任何函数,只要它具有足够多于 n 的参数。 请注意,这个要求在实践中很容易满足。 例如,在 CIFAR 10 (n=50000) 上,成功的残差网络通常具有多个 106 个参数。 更正式地说,对于大小为 n 且具有 r 类的数据集,我们的构造需要 O(nlogn+r2) 参数。 定理3.2给出了正式的陈述。

我们构造中的每个残差层的形式为 x+VReLU(Ux),,其中 UV 是线性变换。 这些层比标准残差层简单得多,标准残差层通常具有两个 ReLU 激活以及两个批量归一化实例。

全卷积残差网络的强大功能。

直接受到我们表达性结果的简单性的启发,我们在 CIFAR10、CIFAR100 和 ImageNet 数据集上尝试了非常相似的架构。 我们的架构只是一系列卷积残差层,每个层都有一个 ReLU 激活,但没有标准架构中常见的批量归一化、丢失或最大池化。 最后一层是未经训练的固定随机投影。 根据我们的理论,卷积权重在 0, 附近初始化,主要使用高斯噪声作为对称破坏器。 唯一的正则化器是标准权重衰减(2-正则化),并且不需要 dropout。 尽管很简单,我们的架构在 CIFAR10 基准(使用标准数据增强)上达到了 6.38% top-1 分类错误。 这与[6]中报告的最佳残差网络具有竞争力,后者实现了6.43% 此外,它还改进了之前由 [15] 实现的最佳全卷积网络 7.25% 的性能。 与我们的不同,之前的全卷积架构还需要对整个数据集进行 dropout 和非标准预处理(ZCA)。 我们的架构还显着改进了 Cifar100 和 ImageNet 上的 [15]

1.2相关工作

自从残差网络([6];[7])出现以来,大多数最先进的图像分类网络都采用了卷积网络的残差参数化层。 [8] 报告了残差网络变体(称为密集网络)的进一步令人印象深刻的改进。 这些网络不是将原始输入添加到卷积层的输出,而是直接通过串联保留原始特征。 在此过程中,密集网络还能够轻松地在更高维空间中对身份嵌入进行编码。 看看我们的理论结果是否也适用于这种残差网络变体将会很有趣。

尽管全面的答案仍然难以捉摸,但在理解神经网络的优化领域方面最近取得了进展。 [5][4]中的实验表明,训练目标具有有限数量的具有较大函数值的不良局部最小值。 [3] 的工作在神经网络的优化景观与物理学中的自旋玻璃模型的优化景观之间进行了类比 ([1])。 [14]表明2层神经网络没有不好的可微局部最小值,但他们没有证明一个好的可微局部最小值确实存在。 [2][12]表明线性神经网络没有坏的局部极小值。 相比之下,我们表明深度线性残差网络的优化景观没有坏的临界点,这是一个更强、更理想的特性。 我们的证明也明显更简单,说明了重新参数化优化的力量。 我们的结果还表明,与较浅的网络相比,更深的网络可能具有更理想的优化景观。

2 线性残差网络的优化景观

考虑从噪声测量 y=Rx+ξ, 学习线性变换 R:dd 的问题,其中 ξ𝒩(0,Id)d 维球形高斯向量。 𝒟表示输入数据x,的分布,令Σ=𝔼x𝒟[xx]为其协方差矩阵。

当然,有很多方法可以解决这个经典问题,但我们的目标是深入了解神经网络,特别是残差网络的优化领域。 因此,我们通过一系列权重矩阵 A1,,Ad×d 来参数化我们的学习模型,

h0 =x,hj=hj1+Ajhj1,y^=h. (2.1)

这里 h1,,h11 隐藏层,y^=h 是学习模型对输入 x. 的预测 更简洁地说,我们有

y^=(I+A)(I+A1)x.

很容易看出这个模型可以表达任何线性变换R.我们将使用A作为所有权重矩阵的简写,即×d×d 维张量,包含 A1,,A 作为切片。 我们的目标函数是最大似然估计,

f(A,(x,y))=y^y2=(I+A)(I+A1)xRxξ2. (2.2)

我们将分析人口风险的情况,定义为:

f(A):=𝔼[f(A,(x,y))].

回想一下,AiAi 的谱范数。 我们将张量 A 的范数 |||||| 定义为其切片的谱范数的最大值,

|A|:=max1iAi.

本节第一个定理指出,目标函数f具有较小||||||范数的最优解,其与层数成反比 因此,当架构很深时,我们可以寻找相当小的范数解决方案。 我们定义γ:=max{|logσmax(R)|,|logσmin(R)|} 这里σmin(),σmax()分别表示R的最小和最大奇异值。

Theorem 2.1.

假设3γdet(R)>0 那么,存在一个具有范数的群体风险f()的全局最优解A

|A|(4π+3γ)/.

我们首先注意到条件det(R)>0在以下意义上不失一般性。 给定任何具有负行列式的线性变换 R ,我们可以通过使用附加维度增强数据和标签来有效地翻转行列式:让 x=[x,b]y=[y,b] ,其中b是一个独立的随机变量(例如,来自标准正态分布),并让R=[R001] 然后,我们有 y=Rx+ξdet(R)=det(R)>0222当维度为奇数时,有一个更简单的方法可以看到这一点 - 翻转标签对应于翻转 R,我们有 det(R)=det(R)

其次,我们注意到这里 γ 应该被认为是一个常量,因为如果 R 太大(或太小),我们可以适当缩放数据,以便 σmin(R)1σmax(R). 具体来说,如果σmax(R)/σmin(R)=κ,那么我们可以正确缩放输出,以便σmin(R)=1/κσmax(R)=κ 在本例中,我们有 γ=logκ,对于相当大的条件数 κ,它将保持一个小常数。 我们还指出,我们在分析中没有尝试优化常量因素。 定理2.1的证明比较复杂,推迟到A节。

鉴于对定理 2.1 的观察,我们将注意力限制在分析 A|||||| 集合中的 f() 的景观 -范数小于τ

τ={A×d×d:|A|τ}.

这里使用定理2.1,半径τ应该被认为是1/的量级。 我们在本节中的主要定理声称,对于任何 τ<1,域 τ 中不存在坏临界点。 回想一下,临界点具有消失梯度。

Theorem 2.2.

对于任何τ<1,我们有域τ内目标函数f()的任何临界点A也必须是全局最小值。

定理2.2表明优化器足以收敛到总体风险的临界点,因为所有临界点也是全局最小值。

此外,除了定理2.2之外,我们还有域τ内的任何A满足

f(A)F24(1τ)22σmin(Σ)(f(A)Copt). (2.3)

这里Coptf()的全局最小值,f(A)F表示欧氏范数333TF:=ijkTijk2 ×d×d维张量f(A)的。 注意,σmin(Σ)表示Σ的最小奇异值。

方程 (2.3) 表示,与误差相比,梯度具有相当大的范数,这保证了迭代时梯度下降收敛到全局最小值 ([11])留在域τ,内,这不能由定理2.2本身保证。

为了证明定理2.2,我们从一个简化群体风险的简单主张开始。 我们用F表示矩阵的Frobenius范数,A,B表示标准基中AB的内积(即 A,B=tr(AB),其中 tr() 表示矩阵的迹。)

断言 2.3

在本节的设置中,我们有,

f(A)=((I+A)(I+A1)R)Σ1/2F2+C. (2.4)

这里C是一个不依赖于A的常数,Σ1/2表示Σ的平方根,即满足B2=Σ的唯一对称矩阵B

断言 2.3 的证明。

tr(A)表示矩阵A的迹。 E=(I+A)(I+A1)R 回顾 f(A) 的定义并使用方程 (2.2),我们有

f(A) =𝔼[Exξ2] (by equation (2.2) )
=𝔼[Ex2+ξ22Ex,ξ]
=𝔼[tr(ExxE)]+𝔼[ξ2] (since 𝔼[Ex,ξ]=𝔼[Ex,𝔼[ξ|x]]=0)
=tr(E𝔼[xx]E)+C (where C=𝔼[ξ2])
=tr(EΣE)+C=EΣ1/2F2+C. (since 𝔼[xx]=Σ)

接下来,我们通过简单的矩阵演算计算目标函数 f() 的梯度。 我们将完整的证明推迟到 A 部分。

Lemma 2.4.

f() 的梯度可以写为,

fAi=2(I+Ai+1)(I+A)EΣ(I+A1)(I+Ai1), (2.5)

其中 E=(I+A)(I+A1)R

现在我们准备证明定理2.2 关键的观察是每个矩阵 Aj 具有较小的范数并且不能取消单位矩阵。 因此,方程(2.5)中的梯度是除误差矩阵E之外的非零矩阵的乘积。因此,如果梯度消失,那么唯一的可能是矩阵E消失,这又意味着A是最优解。

定理证明2.2

使用引理 2.4,我们有,

fAiF =2(I+Ai+1)(I+A)EΣ(I+A1)(I+Ai1)F (by Lemma 2.4)
2jiσmin(I+Ai)σmin(Σ)EF (by Claim C.2)
2(1τ)1σmin(Σ1/2)EΣ1/2F. (since σmin(I+A)1A)

它遵循

f(A)F2 =i=1fAiF24(1τ)2(1)σmin(Σ)EΣ1/2F2
=4(1τ)2(1)σmin(Σ)(f(A)C) (by the definition of E and Claim 2.3)
4(1τ)2(1)σmin(Σ)(f(A)Copt).

因此我们完成了方程的证明(2.3)。 最后,如果A是一个临界点,即f(A)=0,那么通过方程(2.3)我们可以得到f(A)=Copt 也就是说,A是全局最小值。

3 残差网络的表征能力

在本节中,我们描述残差网络的有限样本表达能力。 我们考虑具有单个 ReLU 激活且没有批量归一化的残差层。 基本残差构建块是一个函数 𝒯U,V,s():kk,它由两个权重矩阵 Uk×k,Vk×k 和一个偏差向量 sk 参数化,

𝒯U,V,s(h)=VReLu(Uh+s). (3.1)

残差网络由一系列此类残差块组成。 [7]中的完整预激活架构相比,我们在每个构建块中删除了两个批量归一化层和一个ReLU层。

我们假设数据具有 r 标签,编码为 r 中的 r 标准基向量,用 e1,,er 表示。 我们有 n 训练示例 (x(1),y(1)),,(x(n),y(n)),其中 x(i)d 表示第 i 数据,y(i){e1,,er} 表示 i-第一个标签。 不失一般性,我们假设数据已标准化,因此 x(i)=1. 我们还做出温和的假设,即没有两个数据点彼此非常接近。

Assumption 3.1.

我们假设对于每个 1i<jn,我们都有 x(i)x(j)2ρ 对于某个绝对常数 ρ>0.

例如,图像总是可以在像素空间中受到难以察觉的扰动,从而满足小但恒定的 ρ. 的假设

在这个温和的假设下,我们证明只要参数数量是大于 n 的对数因子,残差网络就有能力表达数据的任何可能的标签。

Theorem 3.2.

假设训练示例满足假设3.1 然后,存在一个残差网络N(如下指定),其参数O(nlogn+r2)完美地表达了训练数据,即对于所有i{1,,n},网络Nx(i) 映射到 y(i).

在实践中,n>r2, 很常见,例如 Imagenet 数据集的情况,其中 n>106r=1000.

我们使用等式 (3.1) 中定义的 𝒯U,V,s 形式的构建块构建以下残差网络。 该网络由+1隐藏层h0,,h组成,输出由y^r表示。 第一层权重矩阵A0d维输入映射到k维隐藏变量h0 然后我们应用带有权重矩阵Aj,Bjk×k的构建块𝒯层。 最后,我们应用另一个层将隐藏变量 h 映射到 k 中的标签 y^ 从数学上来说,我们有

h0 =A0x,
hj =hj1+𝒯Aj,Bj,bj(hj1),j{1,,}
y^ =𝒯A+1,B+1,s+1(h).

我们注意这里的A+1k×rB+1r×r以便维度兼容。 我们假设标签数 r 和输入维度 d 都小于 n,这在实际应用中是安全的。444In computer vision, typically r is less than 103 and d is less than 105 while n is larger than 106 超参数 k 将选择为 O(logn),层数选择为 =n/k 因此,第一层具有 dk 参数,中间的每个 构建块包含 2k2 参数,最终构建块具有 kr+r2 参数。 因此,参数总数为O(kd+k2+rk+r2)=O(nlogn+r2)

为了构建适合数据的上述形式的网络 N,我们首先采用一个随机矩阵 A0k×d 将所有数据点 x(i) 映射到向量 h0(i):=A0x(i) 这里我们将使用 hj(i) 来表示第 i 示例的第 j 层隐藏变量。 根据 Johnson-Lindenstrauss 定理([10],或参见 [17]),很有可能,结果向量 h0(i) 仍然满足假设3.1(缩放比例略有不同,常数ρ较大),即任意两个向量h0(i)h0(j)不是很相关的。

然后我们构造 中间层,将每个 i{1,,n}h0(i) 映射到 h(i) 这些向量 h(i) 将根据标签聚集到 r 组中,尽管它们位于 k 中,而不是所需的 r 中。 具体来说,我们通过在k中挑选r随机单位向量q1,,qr来设计这个聚类中心。 我们将它们视为维度 k 中的代理标签向量(请注意,k 可能比 r 小得多)。 在高维度(技术上,如果k>4logr)随机单位向量q1,,qr与小于<0.5的内积成对不相关。 我们将第 i 个示例与定义如下的目标代理标签向量 v(i) 相关联,

if y(i)=ej, then v(i)=qj. (3.2)

然后,我们将构造矩阵 (A1,B1),,(A,B),以便网络的前 层将向量 h0(i) 映射到代理标签向量 v(i) 从数学上讲,我们将构造 (A1,B1),,(A,B) 使得

i{1,,n},h(i)=v(i). (3.3)

最后,我们将构造最后一层 𝒯A+1,B+1,b+1,以便它将向量 q1,,qrk 映射到 e1,,err

j{1,,r},𝒯A+1,B+1,b+1(qj)=ej. (3.4)

将它们放在一起,通过定义 (3.2) 和方程 (3.3),对于每个 i,如果标签为 y(i)ej,那么 h(i) 将是 qj 然后通过方程(3.4),我们得到y^(i)=𝒯A+1,B+1,b+1(qj)=ej 因此我们得到y^(i)=y(i) 该计划的关键部分是构建中间层权重矩阵,以便h(i)=v(i) 我们将其封装到以下非正式引理中。 正式声明和完整证明推迟到B节。

Lemma 3.3 (引理B.2的非正式版本)

在上面的设置中,对于(几乎)任意向量 h0(1),,h0(n)v(1),,v(n){q1,,qr},存在权重矩阵 (A1,B1),,(A,B),这样,

i{1,,n},h(i)=v(i).

我们简要地概述引理的证明以提供直觉,并将完整的证明推迟到 B 节。 每个残差块应用于隐藏变量的操作可以抽象地写为:

h^h+𝒯U,V,s(h). (3.5)

其中 h 对应于块之前的隐藏变量,h^ 对应于块之后的隐藏变量。 我们声称,对于(几乎)任意的向量序列h(1),,h(n),存在𝒯U,V,s(),使得操作(3.5)变换k h(i) 的向量到我们可以自由选择的任意一组其他 k 向量,并保持其余 nk 向量的值。 具体来说,对于任何大小为 k 的子集 S 和任何所需的向量 v(i)(iS),都存在 U,V,s 使得

v(i) =h(i)+𝒯U,V,s(h(i))iS
h(i) =h(i)+𝒯U,V,s(h(i))iS (3.6)

该主张在引理 B.1 中得到了形式化。 我们可以重复使用它来构建 层构建块,每个构建块将 {h0(1),,h0(n)} 中的 k 向量子集转换为 {v(1),,v(n)},并维护他人的价值观。 回想一下,我们有 =n/k 层,因此在 层之后,所有向量 h0(i) 都会转换为 v(i) 的,这完成校样草图。

4 全卷积残差网络的功效

受我们理论的启发,我们在标准图像分类基准上尝试了全卷积残差网络。

4.1 CIFAR10 和 CIFAR100

我们的 CIFAR10 和 CIFAR100 架构是相同的,除了最终维度分别对应于类 10100 的数量。 在表 1 中,我们概述了我们的架构。 每个残差块的形式为x+C2(ReLU(C1x)),,其中C1,C2是指定维度的卷积(内核宽度、内核高度、输入通道数、输出通道数) )。 每个块中的第二个卷积始终具有步幅 1,而第一个卷积可能具有指定的步幅 2 在变换不保持维度的情况下,原始输入 x 使用平均池化和填充进行调整,这是残差层中的标准。

我们使用 Tensorflow 框架训练模型,使用动量优化器,动量为 0.9,,批量大小为 128 所有卷积权重均通过权重衰减0.0001.进行训练 初始学习率为 0.05,,下降了 103000050000 步长。 该模型在大约 50k 步时达到峰值性能,在单个 NVIDIA Tesla K40 GPU 上大约需要 24h 步。 我们的代码可以轻松地从开源实现555https://github.com/tensorflow/models/tree/master/resnet 通过删除批量归一化,调整剩余组件和模型架构。 与代码的一个重要区别是,我们使用标准差 σ=1/k2c, 的随机正态初始化器来初始化内核大小 k×kc 输出通道的残差卷积层,而不是比用于标准卷积层的 1/kc 更好。 这种小得多的权重初始化有助于训练,同时不影响表示。

与标准模型的一个显着区别是最后一层没有经过训练,而只是一个固定的随机投影。 一方面,这稍微改善了测试误差(可能是由于正则化效应)。 另一方面,这意味着我们模型中唯一可训练的权重是卷积的权重,使我们的架构成为“全卷积”。

表格1: CIFAR10/100 的架构(55 卷积,13.5M 个参数)
variable dimensions initial stride description
3×3×3×16 1 1 standard conv
3×3×16×64 1 9 residual blocks
3×3×64×128 2 9 residual blocks
3×3×128×256 2 9 residual blocks
8×8 global average pool
256×𝚗𝚞𝚖_𝚌𝚕𝚊𝚜𝚜𝚎𝚜 random projection (not trained)

我们模型的一个有趣的方面是,尽管其大小为 13.59 万个可训练参数,但即使数据集大小为 50000.,模型似乎也不会过快过拟合。 相比之下,我们发现很难在 CIFAR10 上训练具有这种大小的批量归一化的模型而不会出现明显的过度拟合。

Refer to caption
Refer to caption
图1: CIFAR10(左)和 CIFAR (100)(右)最佳模型的收敛图。 第一步是批量大小 128 的梯度更新。

2 总结了我们的模型与之前工作的非详尽列表相比的 top-1 分类误差,仅限于 之前的最佳全卷积结果[ 15],第一个残差结果[6],以及[8]在CIFAR上的最新结果。 所有结果均采用标准数据增强。

表2: 不同基准上top-1分类误差的比较
Method CIFAR10 CIFAR100 ImageNet

remarks

All-CNN 7.25 32.39 41.2

all-convolutional, dropout
extra data processing

Ours 6.38 24.64 35.29

all-convolutional

ResNet 6.43 25.16 19.38
DenseNet 3.74 19.25 N/A

4.2ImageNet

ImageNet ILSVRC 2012 数据集具有 1,281,167 个数据点和 1000 类。 每个图像的大小都调整为具有 3 通道的 224×224 像素。 我们在 [6] 中尝试了 34 层网络的全卷积变体。 原始模型实现了25.03%分类错误。 我们的派生模型具有 35.7M 可训练参数。 我们使用动量优化器(动量为 0.9)和学习率时间表对模型进行训练,学习率时间表从初始学习率 0.1. 开始,每两个历时衰减 0.94 倍。 训练分布在 6 机器上,异步更新。 每台机器都配备了 8 GPU (NVIDIA Tesla K40),并使用批量大小 256 分割在 8 GPU 上,以便每个 GPU 更新的批量大小 32.

与 CIFAR10 和 CIFAR100 的情况相比,在 ImageNet 上,我们的全卷积模型的表现明显比其原始模型差。 具体来说,我们经历了大量的欠拟合,这表明较大的模型可能会表现更好。

尽管存在这个问题,我们的模型在测试集(50000 数据点)上仍然达到了 35.29% top-1 分类错误,并且 14.17% top-5700,000 步骤后测试训练错误(大约一周)。 虽然不再是最先进的,但这种性能明显优于 [13] 报告的 40.7% 以及 [15] 我们相信,我们模型的更好的学习率计划和超参数设置很可能可以大大提高此处报告的初步性能。

5结论

我们的理论强调了训练深度人工神经网络时身份参数化的重要性。 一个突出的开放问题是将我们的优化结果扩展到非线性情况,其中每个残差都有一个 ReLU 激活,就像我们的表达结果一样。 我们推测类似于定理2.2的结果对于一般非线性情况是正确的。 与标准参数化不同,我们认为获得这样的结果没有根本障碍。

我们希望我们的理论和实验共同帮助简化深度学习的状态,旨在用一些基本原理来解释深度学习的成功,而不是需要巧妙组合的大量技巧。 我们相信,图像识别方面的大部分进步可以仅通过残差卷积层和 ReLU 激活来实现。 这可能会导致极其简单(尽管很深)的架构,与所有图像分类基准的最新技术相匹配。

致谢:我们感谢 Jason D. Lee、Qixing Huang 和 Jonathan Shewchuk 进行了有益的讨论,并善意地指出了本文早期版本中的错误。 我们还感谢 Jonathan Shewchuk 建议改进方程 (2.3),并将其纳入当前版本。 马腾宇感谢Dodds Fellowship和Siebel奖学金的支持。

参考

  • [1] Antonio Auffinger, Gérard Ben Arous, and Jiří Černỳ. Random matrices and complexity of spin glasses. Communications on Pure and Applied Mathematics, 66(2):165–201, 2013.
  • [2] P. Baldi and K. Hornik. Neural networks and principal component analysis: Learning from examples without local minima. Neural Netw., 2(1):53–58, January 1989.
  • [3] Anna Choromanska, Mikael Henaff, Michael Mathieu, Gérard Ben Arous, and Yann LeCun. The loss surfaces of multilayer networks. In AISTATS, 2015.
  • [4] Yann N Dauphin, Razvan Pascanu, Caglar Gulcehre, Kyunghyun Cho, Surya Ganguli, and Yoshua Bengio. Identifying and attacking the saddle point problem in high-dimensional non-convex optimization. In Advances in neural information processing systems, pages 2933–2941, 2014.
  • [5] I. J. Goodfellow, O. Vinyals, and A. M. Saxe. Qualitatively characterizing neural network optimization problems. ArXiv e-prints, December 2014.
  • [6] Kaiming He, Xiangyu Zhang, Shaoqing Ren, and Jian Sun. Deep residual learning for image recognition. In arXiv prepring arXiv:1506.01497, 2015.
  • [7] Kaiming He, Xiangyu Zhang, Shaoqing Ren, and Jian Sun. Identity mappings in deep residual networks. In Computer Vision - ECCV 2016 - 14th European Conference, Amsterdam, The Netherlands, October 11-14, 2016, Proceedings, Part IV, pages 630–645, 2016.
  • [8] Gao Huang, Zhuang Liu, and Kilian Q. Weinberger. Densely connected convolutional networks. CoRR, abs/1608.06993, 2016.
  • [9] Sergey Ioffe and Christian Szegedy. Batch normalization: Accelerating deep network training by reducing internal covariate shift. In Proceedings of the 32nd International Conference on Machine Learning, ICML 2015, Lille, France, 6-11 July 2015, pages 448–456, 2015.
  • [10] William B Johnson and Joram Lindenstrauss. Extensions of lipschitz mappings into a hilbert space. Contemporary mathematics, 26(189-206):1, 1984.
  • [11] H. Karimi, J. Nutini, and M. Schmidt. Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-\L{}ojasiewicz Condition. ArXiv e-prints, August 2016.
  • [12] K. Kawaguchi. Deep Learning without Poor Local Minima. ArXiv e-prints, May 2016.
  • [13] Alex Krizhevsky, Ilya Sutskever, and Geoffrey E Hinton. Imagenet classification with deep convolutional neural networks. In Advances in neural information processing systems, pages 1097–1105, 2012.
  • [14] D. Soudry and Y. Carmon. No bad local minima: Data independent training error guarantees for multilayer neural networks. ArXiv e-prints, May 2016.
  • [15] J. T. Springenberg, A. Dosovitskiy, T. Brox, and M. Riedmiller. Striving for Simplicity: The All Convolutional Net. ArXiv e-prints, December 2014.
  • [16] Eric W. Weisstein. Normal matrix, from mathworld–a wolfram web resource., 2016.
  • [17] Wikipedia. Johnson–lindenstrauss lemma — wikipedia, the free encyclopedia, 2016.

附录 A 2 部分缺少证明

在本节中,我们给出了定理2.1和引理2.4的完整证明,在2节中省略了它们。

A.1 定理证明2.1

事实证明,如果假设R是一个对称正半定(PSD)矩阵,或者如果我们允许变量是复矩阵,那么证明会容易得多。 这里我们首先给出第一个特殊情况的证明草图。 读者可以跳过它并跳转到下面的完整证明。 对于特殊情况,我们还将证明更强的结果,即 |A|3γ/

R是PSD时,它可以通过正交矩阵U进行对角化,在R=UZU的意义上,其中Z=diag(z1,,zd)是对角矩阵非负对角线条目z1,,zd A1==A=Udiag(zi1/)UI,然后我们有

(I+A)(I+A1) =(Udiag(zi1/)U)=Udiag(zi1/)U (since UU=I)
=UZU=R.

我们看到由 A 定义的网络重建了变换 R,因此它是总体风险的全局最小值(正式参见下面的声明 2.3)。 接下来,我们验证每个 Aj 具有较小的谱范数:

Aj =IUdiag(zi1/)U)=U(Idiag(zi)1/)U=Idiag(zi)1/ (since U is orthonormal)
=maxi|zi1/1|. (A.1)

σmin(R)ziσmax(R)开始,我们有3γ|logzi| 它遵循

|zi1/1|=|e(logzi)/1|3|(logzi)/|3γ/. (since |ex1|3|x| for all |x|1)

然后使用方程 (A.1) 和上面的方程,我们得到 |A|maxjAj3γ/,它完成了特殊情况的证明。

为了充分证明定理2.1,我们从以下声明开始:

Claim A.1.

假设 Q2×2 是正交矩阵。 那么对于任意整数q,存在矩阵W1,,Wq2×2和对角矩阵Λ满足(a)Q=W1WqΛWjIπ/q,(b) Λ 是一个对角矩阵,对角线上有 ±1,并且 (c) 如果 Q 是旋转,则 Λ=I.

证明。

我们首先考虑Q是旋转的情况。 每个旋转矩阵可以写为T(θ):=[cosθsinθsinθcosθ] 假设Q=T(θ) 然后我们可以取W1==Wq=T(θ/q)Λ=I 我们可以验证

WjIθ/q.

接下来,我们考虑 Q 是反射的情况。 那么我们可以将 Q 写为 Q=T(θ)diag(1,1),其中 diag(1,1) 是相对于 y 轴的反射。 然后我们就可以取W1==Wq=T(θ/q)Λ=diag(1,1)来完成证明。

接下来给出定理2.1的正式完整证明。 主要思想是减少块对角线情况并应用上述权利要求。

定理证明2.1

R=UKVR 的奇异值分解,其中 UV 为两个正交矩阵,K是对角线上具有非负项的对角矩阵。 由于det(R)=det(U)det(K)det(V)>0det(K)>0,我们可以适当翻转U,V,以便det(U)=det(V)=1 由于 U 是一个正矩阵(即 U 满足 UU=UU),根据权利要求 C.1,我们可以将 U 通过正交矩阵 S 分块对角化为 U=SDS1,其中 D=diag(D1,,Dm) 是一个实数分块对角矩阵,每个分块 Di 的大小最多为 2×2 使用声明A.1,我们可以得出对于任何Di,都存在Wi,1,,Wi,q,Λi,使得

Di=Wi,1Wi,qΛi (A.2)

Wi,jIπ/q Λ=diag(Λ1,,Λm)Wj=diag(W1,j,Wm,j) 我们可以将方程 (A.2) 重写为

D=W1WqΛ. (A.3)

此外,我们有 Λ 是一个对角矩阵,对角线上有 ±1 由于 Wi,j 是行列式为 1 的正交矩阵,因此我们有 det(Λ)=det(D)=det(U)=1 也就是说,Λ 在对角线上有偶数个 1 然后我们可以将 1 分组为 2×2 块。 请注意,[1001] 是旋转矩阵T(π) 因此,我们可以将 Λ 写为对角线上的 +1 和块 T(π) 的串联。 然后应用声明A.1(在每个块T(π)上),我们得到W1,,Wq使得

Λ=W1Wq (A.4)

其中WjIπ/q 因此,使用方程 (A.3) 和 (2.3),我们得到

U=SDS1=SW1S1SWqS1SW1S1SWqS1.

此外,对于每个 jSWjS1I=S(WjI)S1=WjIπ/q,我们都有它,因为 S 是正交矩阵。 对于Wj也可以证明同样的道理。 因此让 Bj=SWjS1I 代表 jqBj+q=SWjS1I,我们可以重写,

U=(I+B1)(I+Bq).

我们可以通过将 V 分解为 2q 矩阵(与单位矩阵 π/q 接近)来类似地处理 V

V=(I+B1)(I+B2q).

最后,我们处理对角矩阵K。让K=diag(ki) 我们有minki=σmin(R),maxki=σmax(R) 然后,我们可以写K=(K)p,其中K=diag(ki1/p)p是稍后选择的整数。 我们有 KImax|ki1/p1|max|elogki1/p1| pγ=max{logmaxki,logminki}=max{logσmax(R),logσmin(R)} 时,我们有

KImax|elogki1/p1|3max|logki1/p|=3γ/p. (since |ex1|3|x| for |x|1)

B1′′==Bp′′=KI,然后我们有K=(I+Bp′′)(I+B1′′) 最后,我们选择p=3γ4π+3γq=π4π+3γ 666这里为了符号方便,不选择p,q为整数。 但是将它们四舍五入到最接近的整数将会以小的常数因子改变范数的最终界限。 并让Ap+4q=B2q,=Ap+2q+1=B1,Ap+2q=Bp′′,,A2q+1=B1′′,A2q=B2q,,A1=B1 我们有 4q+p=

R=UKV=(I+A)(I+A1).

此外,根据需要,我们还有 |||A|||max{Bj,Bj.Bj′′}max{π/q,3γ/p}4π+3γ

A.2 引理证明2.4

我们根据定义计算部分梯度。 Δjd×dAj 的无穷小变化。 使用声明 2.3,考虑 f(A1,,A+Δj,,A) 的泰勒展开式

f(A1,,A+Δj,,A)
=((I+A)(I+Aj+Δj)(I+A1)R)Σ1/2F2
=((I+A)(I+A1)R)Σ1/2+(I+A)Δj(I+A1)Σ1/2F2
=((I+A)(I+A1)R)Σ1/2F2+
2((I+A)(I+A1)R)Σ1/2,(I+A)Δj(I+A1)Σ1/2+O(ΔjF2)
=f(A)+2(I+Aj+1)(I+A)EΣ(I+A1)(I+Aj1),Δj+O(ΔjF2).

根据定义,这意味着 fAj=2(I+A)(I+Aj+1)EΣ(I+Aj1)(I+A1)。 ∎

附录B3部分缺少证明

在本节中,我们提供定理3.2的完整证明。 我们从以下引理开始,构建一个构建块 𝒯,将任意 n 向量序列的 k 向量转换为任意向量集,并且 main他人的价值。 为了更好的抽象,我们使用 α(i),β(i) 来表示向量序列。

Lemma B.1.

S[n] 的大小为 k。假设 α(1),,α(n)n 向量的序列,对于每个 1in 满足 a),我们有 1ραi21+ρ,并且 b) 如果 ijS 至少包含 i,j 之一,然后是 α(i)α(j)3ρ β(1),,β(n) 为任意向量序列。 那么,存在U,Vk×k,s,这样对于每个iS,我们有𝒯U,V,s(α(i))=β(i)α(i),此外,对于每个i[n]\S,我们有𝒯U,V,s(α(i))=0

我们可以看出结论意味着

β(i) =α(i)+𝒯U,V,s(α(i))iS
α(i) =α(i)+𝒯U,V,s(α(i))iS

这是方程 (3.6) 的另一种写法。

引理B.1的证明。

不失一般性,假设S={1,,k} 我们构造 U,V,s 如下。 U 的第 i 行为 i[k]α(i),并令 s=(12ρ)𝟏 其中 𝟏表示全1的向量。 Vi 列为 i[k]1α(i)2(12ρ)(β(i)α(i)) 接下来我们验证一下构造的正确性。 我们首先考虑1ik 我们知道 Uα(i) 是一个向量,其 i 坐标等于 α(i)21ρ Uα(i) 的第 j 坐标等于 α(j),α(i),可以使用引理假设将其上限设置为

α(j),α(i)=12(α(i)2+α(j)2)α(i)α(j)21+ρ3ρ12ρ. (B.1)

因此,这意味着 Uα(i)(12ρ)𝟏 包含单个正条目(其值至少为 α(i)2(12ρ)ρ),而所有其他条目均为非正条目。 这意味着 ReLu(Uα(i)+b)=(α(i)2(12ρ))ei 其中 ei 是第 i 个自然基向量。 由此可见VReLu(Uα(i)+b)=(α(i)2(12ρ))Vei=β(i)α(i) 最后,考虑ni>k。那么与式(B.1)中的计算类似,Uα(i)是所有坐标都小于12ρ的向量。 因此 Uα(i)+b 是一个具有负项的向量。 因此我们有 ReLu(Uα(i)+b)=0,这意味着 VReLu(Uα(i)+b)=0

现在我们准备好声明引理 3.3 的正式版本。

Lemma B.2.

假设 n 向量 z(1),,z(n) 的序列满足假设 3.1 的宽松版本:a) 对于每个 i,1ρz(i)21+ρ b) 对于每个 ij,我​​们有 z(i)z(j)2ρ; v(1),,v(n) 被定义在上面。 然后存在权重矩阵(A1,B1),,(A,B),这样给定i,h0(i)=z(i),我们有,

i{1,,n},h(i)=v(i).

我们将重复使用引理B.1来构造积木𝒯Aj,Bk,sj(),从而证明引理B.2 每个构建块 𝒯Aj,Bk,sj() 获取 {z(1),,z(n)}k 向量的子集,并将它们转换为 v(i) 的,同时将所有其他向量保持为固定的。 由于它们完全是 n/k 层,因此我们最终将所有 z(i) 映射到目标向量 v(i)

引理B.2的证明。

我们重复使用引理B.1 S1=[1,,k] 然后使用引理B.1α(i)=z(i)β(i)=v(i)对于i[n],我们得到存在A1,B1,b1这样对于ik,它保持h1(i)=z(i)+𝒯A1,B1,b1(z(i))=v(i),对于ik,它保持h1(i)=z(i)+𝒯A1,B1,b1(z(i))=z(i) 现在我们归纳构建其他层。 我们将构建这样的图层:图层 j 的隐藏变量满足 hj(i)=v(i) 中的每一个 1ijk,以及 hj(i)=z(i) 中的每一个 ni>jk 假设我们已经构造了第一个 j 层,接下来我们使用引理 B.1 构造 j+1 层。 然后我们认为α(1)=v(1),,α(jk)=v(jk)α(jk+1)=z(jk+1),,α(n)=z(n)S={jk+1,,(j+1)k}的选择满足引理B.1的假设。 事实上,因为 qi 是均匀随机选择的,所以我们对每个 siqs,z(i)1ρ 都有 w.h.p。 因此,从 v(i){q1,,qr} 开始,我们就知道 v(i) 也与任何 z(i) 都不相关。 然后我们应用引理B.1并得出结论,存在Aj+1=U,Bj+1=V,bj+1=s使得𝒯Aj+1,bj+1,bj+1(v(i))=0对于ijk,𝒯Aj+1,bj+1,bj+1(z(i))=v(i)z(i) 代表 jk<i(j+1)k,𝒯Aj+1,bj+1,bj+1(z(i))=0 代表 ni>(j+1)k。这些意味着

hj+1(i) =hj(i)+𝒯Aj+1,bj+1,bj+1(v(i))=v(i)1ijk
hj+1(i) =hj(i)+𝒯Aj+1,bj+1,bj+1(z(i))=v(i)jk+1i(j+1)k
hj+1(i) =hj(i)+𝒯Aj+1,bj+1,bj+1(z(i))=z(i)(j+1)k<in

因此,我们构建了满足层 j+1 归纳假设的 j+1 层。 因此,通过归纳我们得到了所有层,最后一层对于每个示例i都满足h(i)=v(i)

现在我们准备按照3节中概述的总体计划证明定理3.2

定理证明3.2

我们使用形式化下面讨论的直觉定理3.2 First, take k=c(logn)/ρ2 for sufficiently large absolute constant c (for example, c=10 works), by Johnson-Lindenstrauss Theorem ([10], or see [17]) we have that when A0 is a random matrix with standard normal entires, with high probability, all the pairwise distance between the the set of vectors {0,x(1),,x(n)} are preserved up to 1±ρ/3 factor. 也就是说,我们对每个 i1ρ/3A0x(i)1+ρ/3 和每个 ijA0x(i)A0x(j)ρ(1ρ/3)2ρ/3 都有它。 z(i)=A0x(i)ρ=ρ/3 那么我们就有 z(i) 满足 Lemam B.2 的条件。 我们在k中选择r个随机向量q1,,qr v(1),,v(n)定义为等式(3.2)。 然后通过引理 B.2,我们可以构造矩阵 (A1,B1),,(A,B) 使得

h(i)=v(i). (B.2)

请注意,v(i){q1,,qr}qi 是随机单位向量。 因此,选择α(1)=q1,,α(r)=qrβ(1)=e1,,β(r)=er,并且满足引理B.1的条件,并使用引理B.1 我们得出结论,存在 A+1,B+1,s+1 使得

ej=𝒯A+1,B+1,b+1(vj), for every j{1,,r}.. (B.3)

根据方程(3.2)和方程(B.2)中v(i)的定义,我们得出y^(i)=h(i)+𝒯A+1,B+1,b+1(h(i))=y(i).,其中完成证明。

附录C工具箱

在本节中,我们将陈述两个民间传说的线性代数陈述。 下面的声明应该是已知的,但我们在文献中找不到它。 为了完整性,我们在这里提供证明。

Claim C.1

Ud×d为实数正规矩阵(即满足UU=UU)。 那么,存在一个正交矩阵 Sd×d 使得

U=SDS,

其中D是一个真正的块对角矩阵,由大小最多为2×2的块组成。

证明。

由于 U 是正规矩阵,因此它可酉对角化(有关背景,请参阅 [16])。 因此,d×d中存在酉矩阵V,d×d中存在对角矩阵,使得U具有特征分解U=VΛV* 由于 U 本身是一个实数矩阵,因此特征值(Λ 的对角线条目)以共轭对的形式出现,特征向量(V)。 也就是说,我们可以将V的列分组为(v1,v¯1),,(vs,vs¯),vs+1,,vt对,并令相应的特征值为λ1,λ¯1,,λλs,λ¯s,λs+1,,λt 这里λs+1,,λt 然后我们就得到了U=i=1s2(viλivi*)+i=s+1tviλivi Qi=(viλivi*),则我们有Qi是一个2阶实数矩阵。 Sid×2Qi 列跨度的正交基,然后我们可以将 Qi 写为 Qi=SiDiSi,其中 Di 是一个 2×2 矩阵。 最后让S=[S1,,Ss,vs+1,,vt]D=diag(D1,,Ds,λs+1,,λt)我们完成证明。

下面的权利要求用于证明定理2.2 我们在这里提供完整性证明。

Claim C.2 (民间传说)

对于任意两个矩阵 A,Bd×d,我们有

ABFσmin(A)BF.
证明。

由于 σmin(A)2AA 的最小特征值,因此我们有

BAABBσmin(A)2IB.

因此,可以得出

ABF2 =tr(BAAB)tr(Bσmin(A)2IB)
=σmin(A)2tr(BB)=σmin(A)2BF2.

两边取平方根完成证明。