深度学习概述:从基础概念、计算步骤到调优方法

25f783c52fdc6550c090d59f012bb4b8.png

自2012年深度学习崛起以来,这项技术已被应用到计算机视觉、NLP、推荐等各领域,推动了技术发展,创造了巨大的商业价值。同时,深度学习的算法和网络结构在这期间也在不断发展,但不管怎样,它们都遵循深度学习的基础理论。

温故而知新,如今回顾深度学习的基本原理,依然对我们理解它从何发展而来,为什么能发挥作用至关重要,也对我们在下个十年推进深度学习有指导性意义。

因此,本文将重点回顾深度学习相关的基础理论知识。首先介绍深度学习与人工智能的相关概念,然后从感知机等初等神经网络结构出发,阐述深度学习模型的前向传播和反向传播计算理论。同时系统性地介绍常见的前馈神经网络和卷积神经网络,以及深度学习实践中模型常用的调优方法,包括多种权重学习的最优化算法和模型过拟合的调整方法。

1

深度学习与人工智能

深度学习(Deep Learning, DL)作为人工智能的一个重要子分支,说到深度学习,不得不提人工智能(Artificial Intelligence, AI)和与之密切相关的机器学习(Machine Learning, ML)。刚接触人工智能、机器学习和深度学习的读者可能比较容易混淆三者之间的关系,它们之间的相互包含关系如图1-1所示。机器学习是一种实现人工智能的方法,深度学习是一种实现机器学习的技术。下面简要介绍相关概念。

fd64b32ceeda0f630eeee49321f1463a.png

图1-1 人工智能、机器学习、深度学习的包含关系

人工智能是在1956年由约翰·麦卡锡提出的概念,研究用于模拟、延伸和扩展人类智能的科学技术,属于现代计算机科学的一个重要分支。人工智能试图让计算机拥有人类的智慧,即具备理解语言、学习、记忆、推理、决策等诸多能力。人工智能延伸出了很多子领域,包括机器人、语音识别、图像识别、自然语言处理和专家系统等。

机器学习是实现人工智能的重要技术,采用算法解析观测到的大量数据,从中学习出更具一般性的规律,然后对真实世界中的事件作出预测。典型的机器学习算法包括决策树、随机森林、逻辑回归、支持向量机、朴素贝叶斯等。机器学习领域有一个经典共识,即数据和特征决定了机器学习性能的上限,而模型和算法只是不断朝着这个上限逼近。

而在模型和算法设计过程中,传统机器学习需要投入大量的人力在特征工程上,而理想的状态是让机器帮助我们自动找出应该使用的特征空间,无须人参与。为此,人们希望设计的机器学习算法能够自动学习特征和任务之间的关联,还能从简单特征中提取复杂特征,深度学习就是满足这个特点的机器学习算法。

深度学习的概念源于人工神经网络的研究,主要通过组合和抽取低层特征,形成更加抽象的高层表示属性类别或特征,以发现数据的分布特征表示。

2

感知机与神经网络

2.1 单层感知机

神经网络的概念源于生命科学中的神经系统。在生命科学中,神经元是动物脑神经系统中最基本的单元,数百亿的神经元相互连接,组成复杂的神经系统,用来完成学习、认知和体内对生理功能活动的调节。

840f817709827866014797e6d85b811c.png

图1-2 生物神经单元与感知机

如图1-2所示,神经细胞按照功能大致可分为树突、细胞体和轴突。按照对逻辑电路的理解,每个神经细胞可被视为一个只有兴奋或者抑制两种状态的器件,当某个神经元从其他神经细胞接收到的信号强度超过某个阈值时,细胞体就会兴奋,产生电脉冲并传递到其他神经元。

受生物神经细胞的启发,计算机科学家提出感知机(Perceptron)来模拟动物神经细胞,对人工神经网络的发展具有里程碑式的意义。感知机可以认为是一种人工神经元,其本质是一种线性模型,它接收多个信号输入,产生一个信号输出,如图1-2(b)所示。其数学含义可以表达为

8af371799f10d72f8366518ba9d22023.png

(1.1)

其中,b 是实数,称为偏置,包含神经细胞信息的阈值,04d511ac626fcc1547eeb5a6b42cd0d9.png可以视为各个树突传入信息的求和。f(x)为阶跃函数,满足

62e332388ebb43d083dd77b3053e0972.png

(1.2)

其函数图像如图1-3所示。

ba0843f63161ebfc9c4d90052fb60c6e.png

图1-3 单位阶跃函数

感知机输出两个数值结果,因此可以用于解决一些二分类问题。对于输入信息x1x2,…,xn,处理方式为

2b0eec77149bd592e9f985301fad82b3.png

(1.3)

在几何数学中,式(1.3)可以认为是一个超平面方程。因此,感知机可以认为是采用一个超平面将n 维空间中的数据切分为两部分。当然,现实中大部分的数据并不能够恰好被一个超平面分割。如果一个数据集的正负样本能够被一个超平面区分开,那么称这个数据集是线性可分的。总而言之,感知机已经可以用于解决部分线性问题,但是其结构相对简单,存在不能处理线性不可分问题的缺陷。

2.2 多层感知机

为了进一步挖掘感知机的能力,20世纪80年代,多层感知机(Multilayer Perception, MLP)被提了出来。多层感知机是单个感知机的推广,用来克服感知机不能对线性不可分数据进行识别的弱点。多层感知机在单层感知机的基础上引入了一到多个隐藏层,基本结构由三层组成。第一层为输入层(Input Layer),第二层为隐藏层(Hidden Layer),第三层为输出层(Output Layer),如图1-4所示。

b1832e7ed89368a6a2ef2a2dc48a53b5.png

图1-4 多层感知机示例

多层感知机相比单个感知机,除层数多之外,还增加了激活函数(Activation Function)。类似于图1-3中的阶跃函数的作用,神经元节点对上层输入权重求和之后,经历一个函数变换后作为最终输出,这个变换称为激活函数。阶跃函数就是一种激活函数,除此之外,神经网络中常见的激活函数还有 Sigmoid和线性整流(Rectified Linear Unit,ReLU)函数,如图1-5所示。下面分别介绍。

3a0e72d48d9e1df3a58e376d6868dee7.png

图1-5 激活函数

Sigmoid函数的表达式为

e6e8dd653b02d29b91ee66693c4d212b.png

(1.4)

Sigmoid函数的输出值映射在(0,1)之间,单调连续,可求导数,如图1-5(a)所示。但是包含指数计算和除法计算,导致计算复杂度高,同时求导后容易产生梯度趋近于0的问题,不利于神经网络的训练。

ReLU 函数可以有效解决梯度消失问题,其函数形式为

c91265da6d1755928d5c811e8b3646c0.png

(1.5)

如图1-5(b)所示,当输入为正时,输出为正,当输入为负时,输出为0。

对于激活函数,需要注意的一点是,若每个神经元的激活函数都是线性函数,那么任意层数多层感知机都可被简化成一个等价的单层感知器。下面做简单证明,假设线性激活函数为f1(x)=k1x+c1,f2(x)=k2x+c2,那么隐藏层为

e82f49d70e08da443d37fdf9ea098cc1.png

(1.6)

输出层为

b0f2a3047405726eaf4c622bd0e29c09.png

(1.7)

下面用矩阵形式简化表达,令

5b42182b5c7d7c365f35f1a6aaa9dd08.png

则式(1.6)、式(1.7)可以简化为

6483182236982fd8927b7c7b4cb19855.png

(1.8)

                      ac9a20c0ef72bcd00f23805cad678067.png

                           6cc4e2a290ccd209efedd93bacf71209.png                          297bcf1f900b168e0ecbd54ebf811999.png

(1.9)

其中,C 为偏置参数矩阵。由于2641eeb380ad6ad6279cd5607e3296eb.png,则图1-4 中的三层结构可以压缩为两层感知机,采用数学归纳法,可以证明对每个神经元的激活函数都是线性函数的多层神经网络,故压缩为单层感知机。

3

前馈神经网络

前馈神经网络(Feedforward Neural Network, FNN)通常由一个输入层、多个隐藏层和一个输出层构成,典型结构如图1-6所示。每层的神经元可以接收前一层神经元的信号,并产生输出到下一层,信号从输入层向输出层单向传播。其中隐藏层可代替人工特征工程进行自动特征提取,其数量以及每层的神经元个数可根据具体任务设定。前馈神经网络结构简单,能够以任意精度逼近任何连续函数,应用广泛。

083af8eb00855a2941b746058f438adb.png

图1-6 典型的深度神经网络结构

3.1 前馈神经网络的模型

如图1-6所示,若第(l-1)层神经元为75e6cbb85ec6d10fe8f25da1a161c3d1.png,l层神经元为c2af038431ed62a7cb0713b2b290fa15.png,其层间权重参数矩阵为5dfea41a99a703b92e32e6c5bf93d4e7.png,偏置向量为6815063294428db670fb55d0663cee9c.png,激活函数为f1185782a39648430a9f902407f06c46.png,则两层之间的传递计算关系为

3f0f3165b5362dcbab010cfc348368c4.png

(1.10)

以此迭代关系,可完成隐藏层前向计算。

输出层的激活函数往往需要与具体的任务相结合,深度学习中常见的任务有回归(Regression)任务和分类(Classification)任务。如果函数的输出是一个标量(即一个单独的数),我们定义其为回归任务;如果函数的输出为有限的几种可能(例如图片分类),我们定义其为分类任务;回归任务的激活函数为恒等函数,即直接输出输入信号线性变换后的值。分类任务的激活函数常常采用Softmax函数来实现多分类。

Softmax分类模型会有多个输出,且输出个数与类别个数相等,Softmax 函数定义如下:

70aa3bd711bb944c684ee3645692d2a0.png

(1.11)

其中,928bc4bd007f2537c0b911c6af504382.png为样本属于第 k 类的概率,4671dcc3349e66757ee7bd10c14f3482.png表示最后的隐藏层神经元值。Softmax函数计算示意图如图1-7所示。

52c94cc3a4026c04f34aedc3de4f9d3f.png

图1-7 Softmax激活函数计算示意图

3.2 前馈神经网络的学习

一般而言,对于已经搭建好的多层神经网络,输入的信息依次通过所有神经网络层,最终得到输出结果的过程称为前向传播(Forward Propagation)。因为神经网络的权重参数往往是随机初始化的,导致预测得到的结果与真实结果会存在误差。

现在我们希望根据误差来调节权重参数和偏置,使得预测值更加接近真实值,这种采用“误差反向传播算法”的方法简称反向传播(Back Propagation, BP)。反向传播是深度神经网络中更新权重参数最常用的方式。反向传播是利用预测值与真实值的误差来反向逐层计算权重梯度并进行权重参数的迭代,最终达到收敛。

1. 前向传播

一个完整的神经网络包括前向传播过程和后向传播过程,如图1-8所示。

ef699ae6485e76def68198a51e08a305.png

图1-8 前向、后向传播过程

其中4c40458813f2f78bccac5a7f9730d584.png表示第 l 层的第 i 的神经元的输出值。2a9926a3f571fd8e73658d9f82174ee1.png表示参与计算的第 l 层神经元的参数,f07071747176462c643c9806466bbabe.png对应第l层的第 i 个元素与上一层的第 j 个元素的权重系数。例如,第一层隐藏层的计算公式为

871af5a6c59421098ebb487e6eb24c20.png

(1.12)

第二层隐藏层的计算公式为

fcf8bad6b16132e5ac475ce6377aaf52.png

(1.13)

第三层隐藏层的计算公式为

5de2288a3382953780ca949c0d70c812.png

(1.14)

输出层公式为

eb0c30523d35f86924413d5348f1510a.png

(1.15)

其中,a9e3aa6872818abe72a0dc07f6298af2.png表示Sigmoid激活函数,f009146f17fdc8951adff86871403c87.png指经过激活函数的神经元数值。用矩阵的形式可以表示为

72d6868fb0e2d6fc895a0ba542198eae.png

(1.16)

96df2ab2fc7740f347689ec2e127aa33.png

(1.17)

其中,l=1, 2, 3, 4表示网络层数

aac0030e97b917ec3d049db59549ab67.png

2. 损失函数

通过前向过程进行结果预测得到d9be00ba2416619eb6a82110137848ca.png,损失函数(Loss Function)用来估量模型的输出49b8977b27c00e0b7507d389b4f96a78.png与真实值 y 之间的差距。损失函数越好,通常模型的预测效果越好,深度学习过程就是对损失函数求最小化的过程。

所谓真实值,通常指经过人工处理标注数据的值(LabeledData),往往对应为监督学习(Supervised Learning)的任务。损失函数的确定也需要依据具体问题而定,例如,回归问题一般采用欧氏距离(Euclidean Distance),如平方损失(Mean Squared Loss)函数,分类问题常用的损失函数如交叉熵损失(Cross Entropy Loss)函数。

平方损失函数采用均方差(Mean Squared Error, MSE)来衡量预测值与真实值的误差,其基本形式如下:

d87c016b7cbf6cfcb7ae5c5a0514847f.png

(1.18)

平方损失函数的导数连续,计算简单。假设模型预测与真实值之间的误差服从标准高斯分布

7f4f76644c60cd16abe3edbf580ff1da.png

其中aa3233360587a810d557c16e60478df3.png为方差。但是平方损失函数容易受异常点的影响,异常值往往带来极大的平方误差,不利于参数的更新。对于分类问题,均方差分布难以满足高斯分布,因此一般不使用平方损失函数。

交叉熵是信息论中的一个重要概念,主要用于度量两个概率分布间的差异性,往往用于分类问题。对于二分类问题,损失函数

b3c7069e40911517a871e691eae2ae9b.png

(1.19)

对于多分类问题的损失函数,则可以定义为

d6a75b403b25b0ee730feebe567a590e.png

(1.20)

其中,C 为类别的数量,194c61774568d63d35d7eb40bd3cca71.png等于0或者1,如果预测出的类别和样本标记相同则为1,否则为0,05f1c8e609fd5106965f8394d3158383.png为样本属于类别 c 的概率。

在实际深度学习训练任务中,并不一定能保证所有训练样本都是标注好的,若训练样本只有部分数据被标注,这种训练任务一般称为半监督学习任务(Semi-Supervised Learning),半监督学习只对有标注的训练样本求损失函数;若训练样本都没有标注,这种训练任务称为无监督学习(Unsupervised Learning),用于学习样本数据的统计规律或样本数据的内在结构。

3. 反向传播

损失函数给出了预测值与真实值之间的差异,为模型的优化指引方向,即需要调整模型参数505e75bfd4f1b0cda0917ae1abcb1af8.pnga34d1aceaa82e238cb188f28093d322d.png。著名数学家高斯提出的最小二乘法(Least Squares Method)是一个经典的方法,然而对于待定系数太多的问题,最小二乘法需要面对大型矩阵求逆计算量较大的问题,为此前人提出了基于迭代算法的优化方法。

具体而言,为了弥合预测值与真实值y 的差异,将误差进行后向传播,并使用梯度下降算法进行模型参数的调整,即采用反向传播算法,调整权重参数的优化算法可采用随机梯度下降(Stochastic Gradient Descent, SGD)算法。

随机梯度下降算法更新权重参数的具体形式如下

8a6fa425e36e04b62145ebaa3366f56d.png

(1.21)

383922173fcf48467cff1d73eb8bb3a3.png

(1.22)

其中,权重矩阵的偏导b6897f0ee48dcb7ac54127e6bbc7a6c5.png和偏置导数6f7e466b2f1259e27a794c3c48091f7f.png是计算的核心,b9cb4f3195cd02febdb628c125190fbe.png是学习率(Learning Rate),表示梯度更新步伐的大小。由于神经网络是逐层搭建的,f9a521d294bfb975b63a2206783ee952.png9dd79c0db1b1de9d1c771d5331e90940.png一般需要反向逐层计算得到。下面举一个简单例子来说明多元复合函数求导法则(又称为多元函数的链式法则)。假设可微函数z=f(u,v),且u=h(s,t),v=g(s,t),则z 是s、t 的复合函数。

3311f6772419cdb152ebb717688fe6ab.png

(1.23)

11c5fd3efa294e7730128fbe880262f5.png

(1.24)

依照链式法则可计算各参数梯度,如表1-1所示。

4dc4e93f8470bc8e9a1ba56bf5ad77be.png

表1-1 反向传播计算权重和偏置

对深度神经网络的计算方式,能够总结出什么规律呢? 若以图的角度来审视权重参数的传递路径,可能会更直观。

以图1-8中神经节点fbd45434ae35592e0891dd5bb0e416a2.png相应的权重参数为例,如果将图视为一个有向图,从损失函数 L 出发,有

cd417ae0807bb55f405c488aae548bda.png

dec06c16786ff18b05e8be74cbb6282c.png

f051dc5712b32bb07d07624b5be5c1e5.png

7a2be982bb38296f1c6f98832452b815.png

等四条路径。那么对10aa7ec5569bd5f4f40de34519d4adf6.png产生贡献的6d1fb1de2d0a43ed759afd4c5b15ac14.pngabb7b728e59c6e847538f0182f4fd67d.png1d9e853af9978eda28b63214262f1e65.png可以计算为

31bdda904b0b8ad2349f114d4f5e7376.png

(1.25)

1774c3825f31232b1b103395c0b54722.png

(1.26)

fd8fabbe18d1462a64e1b01222d719af.png

(1.27)

在训练神经网络时,正向传播和后向传播相互依赖。一方面,在正向传播期间认为a1565341939729f7514ba0488d588b58.pngc507f3cdfef98a29b7bace20b65ff54e.png是已知模型参数,用来逐层计算,得出神经元f758f2cc5f4502b730b2680aff8c4a8b.png预测结果62291b1588f435d6a757c3c499d886a3.png。另一方面,反向传播期间模型参数的偏导数37877c536e39e20f896ba7004c426b8a.pngcb634c79a3d7fdc5077afe758606407e.png的计算取决于由正向传播给出的隐藏变量ee312ef6cbb7359cfdbd8041b73f1f7d.png的当前值。更新的模型参数则是由优化算法根据最近迭代的反向传播给出的。

总的来说,在初始化模型参数后,交替使用正向传播和反向传播,利用反向传播给出的梯度更新模型参数,直到模型参数收敛。

4

卷积神经网络

传统方法中,读取图像的像素矩阵后,根据具体任务的需要设计相应的滤波器(Filter)。常用的滤波器是由具有专业背景的人设计出来的。同样是边缘检测的任务,由于滤波器的不同,其结果也具有显著差异。

现实生活中,图像或者视频应用的场景各种各样,为每一种场景设计相应的滤波器成本较高。如何根据任务来自适应求解相应的过滤器,显得十分必要。深度学习的设计理念是从大数据中学习,得到过滤器参数,以解决人工依赖问题。

4.1 图像数据的存储

当给出一幅图像时,首先需要思考图像是如何被计算机存储和理解的。图像由像素构成,每一像素由颜色构成。在计算机中,生成电子图像的最小单元为像素,整个图像是由规则排列的像素点阵构成。

675fd834e0b78ddd87aa2b021e883f79.png

图1-9 彩色图像示意图

在如图1-9所示的彩色图像中,每一个像素点由RGB颜色空间表示。换句话说,每一像素通过三个颜色通道表示:红色、绿色和蓝色。在一个8位图像中,每个颜色通道能取到的颜色值为28=256种,即0~255。越亮的区域,像素值越高,对应着较高的强度。相反,较暗的区域像素值较低,对应较低的强度。

目前智能手机拍摄的照片大都是 1000×1000像素以上的彩色照片,按照像素给出的数据规模可以达到1000×1000×3=3000000的量级。如果直接采用深度神经网络处理图像信息,会存在以下问题:①在全连接神经网络中,第一个隐藏层会接收百万级的输入参数,导致参数规模非常巨大,神经网络训练效果非常差,且容易产生过拟合现象。②直接使用深度神经网络,难以保证图像的局部不变性,即缩放、平移、旋转不变性。因此,设计适合图像数据的神经网络是十分必要的。

4.2 传统图像处理算子

对于图像信息的处理,人们首先想到的是借鉴动物大脑处理视觉信息的方法。Hubel和Wiesel等在20世纪50年代研究猫与猴子的视觉大脑皮层时,发现一些神经元能分别对某一小块视觉区域进行回应。

受生物意义上的图像识别的启发,人们在设计计算机识别图像的算法上,采用先通过“感受野”识别局部信息,然后合并成整个视觉图像的策略。为了模拟视神经细胞对图像的感知,图感知算法通常会采用某个算子f对图像进行处理。对于区块信息的处理,假设算子f一次性处理的区块面积为3×3的像素区域,计算模式如图1-10所示。

2eeb4cc39d48aad24ab3b5936ac72cba.png

图1-10 按“感受野”分片处理图像数据

边缘检测(Edge Detection)是图像处理和计算机视觉的基本问题,其目的是标识数字图像中亮度变化明显的点,图像属性中的显著变化通常反映了属性的重要事件和变化。简单来说,边缘检测就是寻找边缘,即区域或对象之间的边界,也即像素变化最快的位置。在实际的边缘检测中,往往只用到一阶和二阶差分。二阶导数还可以说明灰度突变的类型。

在某些情况下,如灰度变化均匀的图像,只利用一阶导数可能找不到边界,此时二阶导数就能提供很有用的信息。二阶导数对噪声比较敏感,解决的方法是先对图像进行平滑滤波,消除部分噪声,再进行边缘检测。利用二阶导数信息的算法是基于过零检测的,因此得到的边缘点数比较少,有利于后续的处理和识别工作。拉普拉斯边缘检测算子是典型的二阶边缘检测算子,该算子采用二阶导数的差分格式。下面讲解一下该式的推导过程。

aacb8d14655cd2ad7142a74de2e8096f.png

图1-11为拉普拉斯边缘检测算子的实际效果图。

27965994fa9c1e9210d76dd0b6e78e99.png

图1-11 拉普拉斯边缘检测算子过滤效果图

4.3 卷积

在传统方法中,我们所提到的处理边缘的Sobel滤波器和拉普拉斯滤波器均为一种卷积核。通过卷积运算提取图上的特征,如图像像素变化的梯度特征等。同时,一般将卷积运算前后的图数据统称为特征映射(Feature Map),输入图数据称作输入特征图(InputFeature Map),输出图数据称作输出特征图(Output Feature Map)。

图像中的卷积运算又是如何定义的呢? 如图1-12所示为以3×3的两个矩阵计算的形式展示卷积的计算过程。

b6d5f4ce8ae66b4f6e2039042c88175e.png

图1-12 卷积运算示意图

原始图片与特定的神经元做卷积运算,两个3×3的矩阵相乘后再相加,以图1-12为例,有

1×1+6×0+7×(-1)+20×2+5×0+8×(-2)+20×1+4×0+9×(-1)=29

值得注意的是,卷积运算不一定能得到正数,在后续特征图像生成时,会取绝对值或者采用激活函数(如ReLU)处理负值。图1-13更加形象化地展示了传统的图像处理中Sobel算子和拉普拉斯算子等进行卷积计算后的特点。经过Gx运算后,可以观察到在水平方向上有像素梯度的位置会很亮。

同样地,经过Gy运算后,可以观察到在垂直方向上有像素梯度的位置会很亮。拉普拉斯算子没有具体的分量算子,是在矩阵上做二阶导数来提取轮廓。从数学运算上看,卷积作为一种运算方式,用数学语言可以归纳为

448ab784dc59ffb3c2be4d4563ab72d1.png

从计算效果上看,经过上述几种卷积运算后,可以得到不同的处理结果,实现不同的功能。卷积参数深刻影响着卷积结果,将传统图像处理中的卷积核泛化为一般形式,即卷积核为待定系数,通过深度学习的方式来最终确定。

c1d197f87da23335e7057a60a2d8619b.png

图1-13 二维Sobel卷积核与拉普拉斯卷积核运算和可视化

4.4 池化

对于分类任务而言,即使物体局部比较模糊,依然能够对图片进行分类,而这种局部模糊化的处理可以用于压缩数据和参数的量。保留主要的特征同时减少参数(或者说降维,效果类似于主成分分析)和计算量。经过局部模糊化后的图片仍然具有平移不变性和特征不变性。

平移不变性指经过平移、旋转、缩放后,卷积仍能检测到图片的特征。特征不变性,也就是我们在图像处理中经常提到的特征的尺度不变性,例如,一只猫的图片被缩小了数倍,我们还能认出这是一张猫的图片,说明这张图片中仍保留着猫最重要的特征,图像压缩时去掉的信息不会影响整图判断,留下的信息具有尺度不变性的特征,可以表达图像的特征。同时,压缩过的图片像素矩阵大大减小,可以提高运算速度。

那么如何进行这种局部模糊化处理呢? 我们采用一种称为池化的算法。池化函数在计算某一位置的输出时,会计算该位置相邻区域的输出的某种总体统计特征,作为该位置的输出。池化算子又分为最大池化(MaxPooling)和平均池化(Average Pooling),计算过程如图1-14所示。

81914f34e0249bd46ee75e9ef5498363.png

图1-14 二维平均池化与最大池化

平均池化的输出为区域内像素的平均值,而最大池化的输出为区域内像素的最大值。在图像处理中,一般采用最大池化算子。池化算子与卷积算子存在几点差异: ①移动方式有所差异,池化算子每次移动的步幅与其本身窗格大小一致; ②池化过程没有需要学习的参数; ③不改变通道数(高维池化过程)。

一般来说,池化算子包括以下几方面的作用: ①逐渐降低数据体的空间尺寸; ②减少网络中参数的数量; ③减少耗费的计算资源; ④有效控制过拟合。

4.5 填充

填充(Padding)是指在输入高和宽的两侧填充元素。图1-15中,在原输入高和宽的两侧分别添加了值为0的元素,将输入高和宽从4调整为5,使得输出的高与宽由2调整为4。

当没有填充时,原来4×4的矩阵经过卷积后,变为2×2的矩阵,卷积运算对输入矩阵周边的运算次数少于中间的元素,调整后可使周边元素的计算与中间元素的计算次数相当,因此填充可以起到调节输出图像特征大小和输入图像边缘像素计算次数的作用。

c6b25b5bb108809d0efb9971b1190576.png

图1-15 填充示意图

4.6 步幅

从计算过程上看,从左边的输入数据到右侧的输出数据,是通过以卷积核窗口大小为滑动间隔计算完成的。需要注意的是,更为泛化的卷积运算过程,滑动间隔也是可以变化的,我们将每次滑动的行数和列数定义为步幅(Stride),图1-16是步幅为1和2的卷积计算。

8918987387f166c3f774abd4e1b77920.png

图1-16 步幅为1和2的卷积计算实例

填充和步幅都能改变输出特征图的大小,那么具体如何计算二者对于输出图的影响呢?具体而言,填充可以调大输出特征图的尺寸,而增大步幅则会减小输出图的尺寸。假设输入特征图矩阵的尺寸为90fd1bcf13430da1fc6a3b209967722a.png,卷积核大小为8ef84680b82ff2112a1d152c78836f44.png,步幅为S,则计算输出的特征图尺寸6a82191b26d88951dc83fa72647cc2b8.png

c79313f0a71a34db5b3ff3ad4548c58e.png

4.7 典型的卷积神经网络结构

卷积神经网络(Convolutional Neural Network,CNN)的网络结构经过长期发展,目前包含多种网络结构。卷积神经网络主要由三种模块构成: 卷积层、采样层和全连接层。这里介绍图神经网络的先驱工作LeNet-5,1988年由Yann LeCun提出,其中5表示五层结构的意思,用于实现手写字识别功能,其网络结构如图1-17所示,其符号说明如表1-2所示。

2b423201f252f7a988ec717c5fb56ca2.png

图1-17 LetNet-5示意图

e7cb943f9c96e905f248b0b48ae8a37d.png

表1-2 LetNet-5符号说明

LeNet5是早期非常经典的卷积神经网络,也是网络成功商业化的代表,但它的输入图像太小,加上数据不足,在早期并没有在除手写数字识别之外的其他计算机视觉任务上取得大的突破。

4.8 卷积神经网络与多层感知机的差别

多层感知机(MLP)其实是多个全连接层叠加组合而成的神经网络模型。与多层感知机相比,卷积神经网络把全连接层改成了卷积层和池化,也就是把传统的由一个个多层感知机组成的网络层变成由卷积和池化组成的神经网络层,如图1-18所示。

94e1a70e13656a9ba74832f01a139c66.png

图1-18 多层感知机与卷积神经网络的比较

对比一下多层感知机的层和由卷积与池化层组合成的卷积神经网络层: 假设存在一个8×8像素的图像,也就是64像素,同时假设有一个9个单元的全连接层,如图1-19所示。

240a6c4521874c630963b5a95e276f10.png

图1-19 卷积神经网络与多层感知机的数据处理模式对比

对于多层感知机模型,这一层总共需要64×9×2=1152个学习参数。因为每一个连接都对应一个权重W 和一个偏置b。对于同样有9个单元的卷积核则是有几个单元就有几个参数,即共有9个学习参数。因为对于不同的区域,共享同一个卷积核,因此就共用同一组参数。由于卷积核的参数共享,即使图片进行了一定的平移操作,图片的特征同样可以被识别出来,确保图片信息的“平移不变性”。

卷积神经网络相对于多层感知机具有以下三个特点。

(1) 局部连接。卷积计算过程中,卷积核心只与特征图矩阵中卷积核心大小的矩阵块进行计算,分片处理,让每个神经元只与输入数据的一个局部区域连接,该连接的空间大小叫作神经元的感受野,会大大减少网络的参数。而一般的深度神经网络,在处理图像这样的高维输入时,让每个神经元都与前一层中的所有神经元进行全连接,这对于像素较多的图是不现实的。

(2) 权值共享。在卷积计算中使用相同的卷积核,用来实现共享参数。每个卷积核与上一层局部连接,同时每个卷积核的所有局部连接都使用同样的参数,这样会大大减少网络的参数,且各个小区块做卷积运算时,并无相互依赖,故可以并行计算。

(3) 池化。它的作用是逐渐降低数据的空间尺寸,这样就能减少网络中参数的数量,使计算资源耗费变少,也能有效控制过拟合。

5

深度学习训练的最优化算法

深度学习模型训练的目的是寻找模型参数来最小化损失函数,属于典型的最优化问题(Optimization Problem,OP)。优化算法的选择是深度学习的重要环节,即使在相同模型结构和数据集中,不同的优化算法也可能导致不同的训练效果。

当神经网络模型结构较为简单且参数较少时,尚能采用最小二乘法求得参数。然而,模型参数维度往往较高,无法直接给出数学上的解析解。最小二乘法需要计算逆矩阵,但其逆矩阵不存在,则无法直接用最小二乘法。梯度下降为代表的迭代法则不受矩阵可逆的约束,仍然可以使用,因此在神经网络中通常采用这类算法。

假设损失函数19ee6462e57a3f7dbbffbd54e10a7f56.png与权重参数b764611b31eb4f34265b22d5e7e118c9.png的关系可用图1-20所示的函数关系描述。下面我们利用梯度下降法来寻找使损失函数最小的参数50a8dddd51ce06942ba0b9f3ba9d8435.png。随机选取一个初始点4581d498d6fe8c977d561696eef63a7b.png和学习率e52d4f42dded5d44bff5e5c31d39c661.png来训练模型。计算损失函数8fefcf4f34d992c2ff9690d3a75dbe26.pngcc3cdcba7a0ab68e2a3a36a1b312be3d.png点的微分2efd78573c2606905433a7bb672e4ae1.png,如果得到的微分值为负,则增加0fcfbde4102ef974c9582561028ac17f.png的值(向右移动),如果得到的微分值为正,则减小04eb98e8ada82c672e05896ecd1b23ae.png的值(向左移动)。移动的步长(即学习率)用f008039a9375b6079a64de6afdf33b19.png表示,第一步训练结束后,参数更新为caa9ee81be9b9bdcdd6911e92d6e395b.png,接着进一步计算损失函数d9cf41f92c568e76bb0ff7014d4d4df5.png1ec6246480602b4e373e6f28ee53d5c4.png点的微分 6a4d2299668d6309bceaac0026d8c310.png,参数更新为2597918225662f807d6cc2af5e8da2cf.png,如此不断循环,每次参数的更新都向使损失函数更小的方向移动,经过多轮迭代后,4386282ee2adc7b1f97065c6d9837ecc.png落在一个局部最小点2b39b3f2daf179b4e58085113141c88c.png处,损失函数在此处保持稳定不变(微分为0),3bfb76bd19cbd09238ac015272e9bc30.png便是通过梯度下降法求得的相对最优解,如图1-20所示。

439f1acde0b63a6ccbcd9144fe6f3993.png

图1-20 梯度降求解过程

高维函数的损失函数在某点的梯度,梯度下降过程类似于在等高线地图上从某一点沿着等高线法线的方向不断移动,直至找到最低点,如图1-21所示。

996f2d9384f65c150324444c32faa1fc.png

图1-21 梯度曲面

梯度下降算法自身也存在一些问题,如学习速率0b6f4fda6aaf9f7461392599cd0712df.png的选择,鞍点(Saddle Point)问题以及局域最优化的问题。学习速率设置得太小,收敛速度会非常慢,训练时间过长,学习率设置得过大,则会越过最低点,无法达到最低点。若以一阶导数是否为零来作为判定最小值的标准,则会出现鞍点和局域最小值的问题。

在深度学习中,一般假设损失函数为凸函数,而实际中的数据难以保证预测函数满足凸函数条件。而在多元连续函数中,导数为零的点只能称为驻点,是极值点的必要不充分条件,需要结合黑塞矩阵(Hessian Matrix)做进一步判断。鞍点则是在该点导数为零的点,而在该点的函数值在一个方向是函数的极大值点,在另一个方向是函数的极小值点。梯度下降法会沿着梯度下降的方向搜寻极小值,沿着梯度上升的方向搜寻极大值。

若损失函数为非凸,则梯度下降找到的不一定是全局最优解,可能只是局部最优解,需要改变不同初始值来改变优化过程中的搜索路径,以求得最优解。对于鞍点问题,如果模型曲线是凹凸不平的,从不同的初始点出发,可能会落至不同的局部最小点,导致模型性能不一致,更糟糕的情况是模型最终不是陷入某个局部最优点,而是来回振荡或者停留在鞍点。

如图1-22所示,鞍点处的梯度同样为0,但很明显这个位置并不是最优解,此外,鞍点通常被相同误差值的平面所包围,在高维的情形下,这个鞍点附近的平坦区域范围可能非常大,使得梯度下降法很难脱离区域,可能会长时间滞留在该点附近。

14c90fdad60fc4cc37c38b32308b4bcf.png

图1-22 鞍点示意图

梯度下降法有三种常见的变体: 批量梯度下降法(Batch Gradient Descent,BGD)、随机梯度下降法(Stochastic Gradient Descent,SGD)和小批量梯度下降法(Mini-Batch GradientDescent,MBGD)。

这三种变体的区别体现在计算梯度所用的数据量大小上,目的是在参数更新准确率和训练时间之间实现平衡。批量梯度下降法每次迭代都使用整个训练集的数据来计算损失函数对参数的梯度,能准确反应模型的移动方向,其缺点是使用全量数据会导致训练速度慢;随机梯度下降法每次迭代只使用一个训练样本,可大大加快训练速度,其缺点是单个样本求得的梯度携带大量噪声,每次迭代并不是都朝着整体最优的方向移动,所以随机梯度下降法虽然训练速度快,但模型准确率会降低;小批量梯度下降法是对随机梯度下降法的一种优化,每次迭代使用一小批样本,这样可以降低随机梯度的方差,使模型收敛更稳定。

为了抑制随机梯度下降法的振荡,若在梯度下降过程加入惯性,也就是在随机梯度下降法的基础上引入一阶动量,则可以在一定程度上抑制随机梯度下降法的振荡,此改进方法称为动量随机梯度下降法(SGD with Momentum,SGD-M)。基于一阶动量的梯度下降法,其参数的更新不仅依赖于当前梯度,还依赖于历史梯度。可以理解为小球从山上滚下来时,不会在鞍点停留,而是由于惯性继续向前冲过鞍点,甚至有可能冲出局部最小点继续前进的过程。基于动量的梯度下降法的迭代更新规则可表示为

8f34f985e4944fb223c895e12301fb65.png

(1.28)

7d7c664a590fa1b2467c4d7a7c073c1d.png

(1.29)

其中,abc2d77d1797c47c43d5c80222961f9e.png是t-1轮迭代的梯度,b78da9efcf02a54f9e8b610e2115909f.png表示动量的累加,且初始值4c57ebc0cbd1632b218d6202df7f4d0b.png=0,2b03014847454fe3a1c456729a8b4bd5.png为阻尼系数,fae367f67659c16dba23c018fbcbaa30.png为学习速率,ca2be57d1f21e8721898ea605d36a7a3.png为权重矩阵历经第 t 轮迭代后的值。

从式(1.28)、式(1.29)可以看出,动量随机梯度下降法中存在两个额超参数: 27461829e645e94cbb3259e68aab2417.png(阻尼系数)和6efcb0b108c1e2a3773f9d9e86279499.png(学习速率)。

特别地,若f0d578f47ab5a2f214e8263f507a124f.png=0,则动量随机梯度下降法退化为随机梯度下降法,若a72358297a07cda5f35c6dc7de7b976b.png=1,则b4eb6dca733e978fb7625031ecb6caf7.png可直观地表达为af31cc0234eac1891f97184d2805a687.png的累计。

如图1-23所示,迭代之初,c7e172857f95f3566f035b40d1c61f87.png符号未变,故724a8a88734462e21bad48f03caf8a49.png的符号与d32ab1075b84bdcf3557cb1727565fcb.png一致,e328f6f58a30768af337260c460b7ce1.png能促进权重系数的学习; 随着迭代的进行,c361c2793ab1f78e4df387cbde304227.png跨过0变换符号后,b1f24de601aade422cdc9d3c64192ae3.png尚未来得及变换符号,597cac545896c0fc455242c97e04e846.png又能避免132095cf19910f53ff9655b0c10420f1.png偏离极值太远; 综合来看,动量随机梯度下降法相较于随机梯度下降法更能加速学习效率。

811121335e7adffee75cb2b322361db8.png

图1-23 考虑动量的梯度更新过程

梯度自适应优化算法(Adaptive Gradient Algorithm,AdaGrad)独立地适应所有模型参数的学习率,对于每一维的参数空间,采用的学习速率不同,其梯度更新方程如下:

92b63fd19fb296b014b08f718fc31636.png

(1.30)

ad4774b5388c6ed00ba0b86df5b6c2c0.png

(1.31)

fb54cbb70170b14efb606c695d6ccbde.png

(1.32)

其中,☉表示向量的元素积,4daef4499e4468cbb78ab71617c42e00.png为梯度累积变量,初始化为0,e8cc2301691698fa059dad125f8daf66.png是全局学习速率,5ce34a3aa9fe21c19dc7611908e36f1e.png是防止除零操作设定的正常实数值,一般取10-6。梯度自适应优化算法会记录过去所有梯度的平方和,因此,随着迭代次数的增加,更新量会趋近于0,直至完全不更新。即初始阶段激励收敛,后期则逐渐惩罚收敛。

梯度自适应优化算法的学习率不断衰退,会使得很多任务尚未达到最优解,其学习率已经过量减小。为了改善这个问题,有人提出了RMSProp算法。RMSProp算法将这些梯度按元素平方做指数加权移动平均,其他部分则与梯度自适应优化算法保持一致,其梯度更新规则为

4dba9ee2c59140d33cdc22e525740a22.png

(1.33)

c30c87c63891fe73f5134e8e7ad8a0e0.png

(1.34)

13da30d7ab3fa72cdda8c6e75f7855c8.png

(1.35)

其中,0≤61dab79e8d5c5bf152c855445bd1fa4b.png<1为超参数,自变量每个元素的学习速率在迭代过程中不再一直降低,改善了AdaGrad算法。

在实战中,更常被使用的是自适应动量随机优化算法(Adaptive Moment Estimation,Adam)。自适应动量随机优化算法可以理解为加了动量的RMSProp算法,同时使用动量和自适应学习率来加快模型收敛的速度。除了像RMSProp一样存储了过去梯度的平方848cd34407a3f8462fbe1937356aff1d.png的指数衰减平均值,也像动量一样保持了过去梯度cf9fa4116593b0145766127ece54c7db.png的指数衰减平均值。对于动量的累计采用指数移动平均,具体如下:

6e6cdf8dcc42951c687f9ef23082f719.png

(1.36)

69469d9f4cc3067e0a3e38fe61a4ffae.png

(1.37)

其中,80a88c83d0aef5bccf0e209f9e25318a.png=0,8c971fb40c80899926b8566d9f013eb6.png∈[0,1)为移动平均的超参。梯度按元素平方做指数加权移动平均:

188d858e547b27eb550d8dcd0c72e4a9.png

(1.38)

其中,e8a1d5e81af8b37e1a43c56be249c86b.png=0,f165cabea403651458e912920042920d.png∈[0,1)为移动平均的超参。因为f5461c29516f32f17fc978f196893bc3.png053c5c7331c200b7c086bc3a0934dc65.png被初始化为0向量,这样的结果会向0偏置,做了偏差校正,通过计算偏差校正后的99b491061a955368d2ec269e7a5d1024.png6bfe906d6188d1c6a4978c71ae255db1.png来抵消这些偏差。

bcd8985dac13575436702ffef4f6c3fa.png

(1.39)

3ee485a9d0acdebf3574649d0e036b17.png

(1.40)

其梯度更新规则为

6d1ed4cea95c5cf6a22b5c78d08cac78.png

(1.41)

6

深度学习中的过拟合和欠拟合

前面介绍了深度学习的计算步骤,训练完成后,通常需要对生成的神经网络模型进行效果评估。训练完的神经网络往往存在过拟合(Overfitting)或欠拟合(Underfitting)。

在训练神经网络模型前,通常会对样本集进行切分,分为训练集、测试集和验证集。深度学习的最终目标是使神经网络模型在训练数据集上学到的模型参数能很好地适用于“新样本”,也就是做预测。这种适用于新样本的能力,称为泛化能力。

泛化能力强的模型可以在新样本上表现出与训练样本一致的性能,而泛化能力差的模型通常会表现为两种情况:欠拟合和过拟合。如果一个模型在训练数据上表现非常好,但是在新数据集上却并不能取得一致的性能,称为过拟合。反之,如果在训练数据集和新数据集上表现都很差,就是欠拟合。

若发生欠拟合则需要进一步调整模型,例如,通过增加模型深度来增强模型的表现能力,或者通过增加训练次数。发生过拟合的原因与样本数据和神经网络的结构有关。对于数据集,存在一个假设就是训练集上的分布和真实数据集的分布是相同分布的。但是训练集只是真实数据集的一小部分,二者的样本分布可能有偏差。因为训练集可能包含较多噪声,若存在一些较强的噪声,那么模型很容易将它学习为一个强特征,造成过拟合。

在数据集上,可以增大训练集,让训练集和预测数据集的分布尽量接近。对于模型,当神经网络模型参数空间较大时,模型有较强的学习能力,能把包括噪声数据在内的每个数据拟合进来,导致模型过拟合。为了让模型抑制过拟合,典型的方法有早停法(Early Stopping)和正则化(Regularization)方法。

早停法是通过将原始训练数据集划分成训练集和验证集,在训练过程中以某一间隔来不断评估模型在验证集上的表现,当模型在验证集上的表现开始下降时,立即停止训练,这样就能避免继续训练导致的过拟合问题。正则化方法的本质是在损失函数中引入一个正则项,以加大对权重参数的惩罚,抑制过拟合。

正则化方法中常用的正则项为L1和L2正则,即在原来的损失函数上增加权重矩阵的L1范数和L2范数作为惩罚项,对权重系数进行约束。L1正则化对应的修正损失函数J(W, b)为

ced6a70dfe801ae3b4ad8de26c2695ab.png

(1.42)

其中,950c6423e3ef1a4334ed378a5ad80ef3.png为正则约束项,9649f00fb246dafe3673934058210a9b.png是一个超参数,用来控制惩罚项的相对贡献,cdd65b58adacd36e4c32c6006a940dc2.png是权重矩阵元的L1范数求和,也就是权值向量90fce5475a298b43a3f9055fc1821eb5.png中各元素的绝对值之和。每一层的权重矩阵3c320678b57cbfd496c3e823a1b1063e.png的平方和,假设第 l 层的权重矩阵维度553b0cc2ef9abe857d719b82680d73a5.png,计算形式为

32004eddab43c512fd81afeeee97e93f.png

(1.43)

L2正则损失部分则采用权重矩阵的L2范数,计算形式为

961ad63e48ad04a793db50fb5d7219f1.png

(1.44)

正则化模型相对于不加正则项的模型,在每一步更新模型参数优化损失函数时,使权重更靠近0。L1正则化采用绝对值形式的惩罚项,绝对值函数不完全可微,输出权重具有稀疏性,即产生一个稀疏模型,可以用于特征选择。L2范数采用二次函数,是完全可微的,能使权重平滑。

7

小结

深度学习是人工智能领域内最活跃的一个分支,是构建在深度神经网络上的一种模型。本文首先介绍了深度学习与人工智能的相关概念,然后从感知机等简单的神经网络结构出发,扩展出多种激活函数,以及包含多个隐藏层的深度神经网络,进一步介绍了深度学习模型的前向传播和反向传播计算理论。

前馈神经网络与卷积神经网络是深度学习中的两个基础而重要的模型。前馈神经网络并未对数据的形式有太多的要求,而卷积神经网络则是针对图像数据而设定的,卷积神经网络有权重共享、分片连接和池化等重要特征。

实践中,往往不是很容易得到模型的最优解,因此本文也介绍了包括多种权重学习的最优化算法,如随机梯度下降、基于动量的随机梯度下降算法、自适应动量随机优化算法、梯度自适应优化算法等,最后介绍了深度学习在实践中常遇到的过拟合和欠拟合问题。

其他人都在看

欢迎下载体验OneFlow v0.7.0:Stargazers · Oneflow-Inc/oneflow · GitHubOneFlow is a performance-centered and open-source deep learning framework. - Stargazers · Oneflow-Inc/oneflowhttps://github.com/Oneflow-Inc/oneflow/stargazers

猜你喜欢

转载自blog.csdn.net/OneFlow_Official/article/details/125342462
今日推荐