无锡专业网站推广,网络营销组合策略,单片机程序员开发网站,万网站长工具往期回顾 在前面的文章系列文章中#xff0c;我们介绍了全连接神经网络和卷积神经网络#xff0c;以及它们的训练和使用。他们都只能单独的取处理一个个的输入#xff0c;前一个输入和后一个输入是完全没有关系的。但是#xff0c;某些任务需要能够更好的处理序列的信息我们介绍了全连接神经网络和卷积神经网络以及它们的训练和使用。他们都只能单独的取处理一个个的输入前一个输入和后一个输入是完全没有关系的。但是某些任务需要能够更好的处理序列的信息即前面的输入和后面的输入是有关系的。比如当我们在理解一句话意思时孤立的理解这句话的每个词是不够的我们需要处理这些词连接起来的整个序列当我们处理视频的时候我们也不能只单独的去分析每一帧而要分析这些帧连接起来的整个序列。这时就需要用到深度学习领域中另一类非常重要神经网络循环神经网络(Recurrent Neural Network)。RNN种类很多也比较绕脑子。不过读者不用担心本文将一如既往的对复杂的东西剥茧抽丝帮助您理解RNNs以及它的训练算法并动手实现一个循环神经网络。 语言模型 RNN是在自然语言处理领域中最先被用起来的比如RNN可以为语言模型来建模。那么什么是语言模型呢 我们可以和电脑玩一个游戏我们写出一个句子前面的一些词然后让电脑帮我们写下接下来的一个词。比如下面这句 我昨天上学迟到了老师批评了____。 我们给电脑展示了这句话前面这些词然后让电脑写下接下来的一个词。在这个例子中接下来的这个词最有可能是『我』而不太可能是『小明』甚至是『吃饭』。 语言模型就是这样的东西给定一个一句话前面的部分预测接下来最有可能的一个词是什么。 语言模型是对一种语言的特征进行建模它有很多很多用处。比如在语音转文本(STT)的应用中声学模型输出的结果往往是若干个可能的候选词这时候就需要语言模型来从这些候选词中选择一个最可能的。当然它同样也可以用在图像到文本的识别中(OCR)。 使用RNN之前语言模型主要是采用N-Gram。N可以是一个自然数比如2或者3。它的含义是假设一个词出现的概率只与前面N个词相关。我们以2-Gram为例。首先对前面的一句话进行切词 我 昨天 上学 迟到 了 老师 批评 了 ____。 如果用2-Gram进行建模那么电脑在预测的时候只会看到前面的『了』然后电脑会在语料库中搜索『了』后面最可能的一个词。不管最后电脑选的是不是『我』我们都知道这个模型是不靠谱的因为『了』前面说了那么一大堆实际上是没有用到的。如果是3-Gram模型呢会搜索『批评了』后面最可能的词感觉上比2-Gram靠谱了不少但还是远远不够的。因为这句话最关键的信息『我』远在9个词之前 现在读者可能会想可以提升继续提升N的值呀比如4-Gram、5-Gram.......。实际上这个想法是没有实用性的。因为我们想处理任意长度的句子N设为多少都不合适另外模型的大小和N的关系是指数级的4-Gram模型就会占用海量的存储空间。 所以该轮到RNN出场了RNN理论上可以往前看(往后看)任意多个词。 循环神经网络是啥 循环神经网络种类繁多我们先从最简单的基本循环神经网络开始吧。 基本循环神经网络 下图是一个简单的循环神经网络如它由输入层、一个隐藏层和一个输出层组成 纳尼相信第一次看到这个玩意的读者内心和我一样是崩溃的。因为循环神经网络实在是太难画出来了网上所有大神们都不得不用了这种抽象艺术手法。不过静下心来仔细看看的话其实也是很好理解的。如果把上面有W的那个带箭头的圈去掉它就变成了最普通的全连接神经网络。x是一个向量它表示输入层的值这里面没有画出来表示神经元节点的圆圈s是一个向量它表示隐藏层的值这里隐藏层面画了一个节点你也可以想象这一层其实是多个节点节点数与向量s的维度相同U是输入层到隐藏层的权重矩阵读者可以回到第三篇文章零基础入门深度学习(3) - 神经网络和反向传播算法看看我们是怎样用矩阵来表示全连接神经网络的计算的o也是一个向量它表示输出层的值V是隐藏层到输出层的权重矩阵。那么现在我们来看看W是什么。循环神经网络的隐藏层的值s不仅仅取决于当前这次的输入x还取决于上一次隐藏层的值s。权重矩阵 W就是隐藏层上一次的值作为这一次的输入的权重。 如果我们把上面的图展开循环神经网络也可以画成下面这个样子 现在看上去就比较清楚了这个网络在t时刻接收到输入之后隐藏层的值是输出值是。关键一点是的值不仅仅取决于还取决于。我们可以用下面的公式来表示循环神经网络的计算方法 式式 式1是输出层的计算公式输出层是一个全连接层也就是它的每个节点都和隐藏层的每个节点相连。V是输出层的权重矩阵g是激活函数。式2是隐藏层的计算公式它是循环层。U是输入x的权重矩阵W是上一次的值作为这一次的输入的权重矩阵f是激活函数。 从上面的公式我们可以看出循环层和全连接层的区别就是循环层多了一个权重矩阵 W。 如果反复把式2带入到式1我们将得到 从上面可以看出循环神经网络的输出值是受前面历次输入值、、、、...影响的这就是为什么循环神经网络可以往前看任意多个输入值的原因。 双向循环神经网络 对于语言模型来说很多时候光看前面的词是不够的比如下面这句话 我的手机坏了我打算____一部新手机。 可以想象如果我们只看横线前面的词手机坏了那么我是打算修一修换一部新的还是大哭一场这些都是无法确定的。但如果我们也看到了横线后面的词是『一部新手机』那么横线上的词填『买』的概率就大得多了。 在上一小节中的基本循环神经网络是无法对此进行建模的因此我们需要双向循环神经网络如下图所示 当遇到这种从未来穿越回来的场景时难免处于懵逼的状态。不过我们还是可以用屡试不爽的老办法先分析一个特殊场景然后再总结一般规律。我们先考虑上图中的计算。 从上图可以看出双向卷积神经网络的隐藏层要保存两个值一个A参与正向计算另一个值A参与反向计算。最终的输出值取决于和。其计算方法为 和则分别计算 现在我们已经可以看出一般的规律正向计算时隐藏层的值与有关反向计算时隐藏层的值与有关最终的输出取决于正向和反向计算的加和。现在我们仿照式1和式2写出双向循环神经网络的计算方法 从上面三个公式我们可以看到正向计算和反向计算不共享权重也就是说U和U、W和W、V和V都是不同的权重矩阵。 深度循环神经网络 前面我们介绍的循环神经网络只有一个隐藏层我们当然也可以堆叠两个以上的隐藏层这样就得到了深度循环神经网络。如下图所示 我们把第i个隐藏层的值表示为、则深度循环神经网络的计算方式可以表示为 循环神经网络的训练 循环神经网络的训练算法BPTT BPTT算法是针对循环层的训练算法它的基本原理和BP算法是一样的也包含同样的三个步骤 前向计算每个神经元的输出值反向计算每个神经元的误差项值它是误差函数E对神经元j的加权输入的偏导数计算每个权重的梯度。最后再用随机梯度下降算法更新权重。 循环层如下图所示 前向计算 使用前面的式2对循环层进行前向计算 注意上面的、、都是向量用黑体字母表示而U、V是矩阵用大写字母表示。向量的下标表示时刻例如表示在t时刻向量s的值。 我们假设输入向量x的维度是m输出向量s的维度是n则矩阵U的维度是矩阵W的维度是。下面是上式展开成矩阵的样子看起来更直观一些 在这里我们用手写体字母表示向量的一个元素它的下标表示它是这个向量的第几个元素它的上标表示第几个时刻。例如表示向量s的第j个元素在t时刻的值。表示输入层第i个神经元到循环层第j个神经元的权重。表示循环层第t-1时刻的第i个神经元到循环层第t个时刻的第j个神经元的权重。 误差项的计算 BTPP算法将第l层t时刻的误差项值沿两个方向传播一个方向是其传递到上一层网络得到这部分只和权重矩阵U有关另一个是方向是将其沿时间线传递到初始时刻得到这部分只和权重矩阵W有关。 我们用向量表示神经元在t时刻的加权输入因为 因此 我们用a表示列向量用表示行向量。上式的第一项是向量函数对向量求导其结果为Jacobian矩阵 同理上式第二项也是一个Jacobian矩阵 其中diag[a]表示根据向量a创建一个对角矩阵即 最后将两项合在一起可得 上式描述了将沿时间往前传递一个时刻的规律有了这个规律我们就可以求得任意时刻k的误差项 式 式3就是将误差项沿时间反向传播的算法。 循环层将误差项反向传递到上一层网络与普通的全连接层是完全一样的这在前面的文章零基础入门深度学习(3) - 神经网络和反向传播算法中已经详细讲过了在此仅简要描述一下。 循环层的加权输入与上一层的加权输入关系如下 上式中是第l层神经元的加权输入(假设第l层是循环层)是第l-1层神经元的加权输入是第l-1层神经元的输出是第l-1层的激活函数。 所以 式 式4就是将误差项传递到上一层算法。 权重梯度的计算 现在我们终于来到了BPTT算法的最后一步计算每个权重的梯度。 首先我们计算误差函数E对权重矩阵W的梯度。 上图展示了我们到目前为止在前两步中已经计算得到的量包括每个时刻t 循环层的输出值以及误差项。 回忆一下我们在文章零基础入门深度学习(3) - 神经网络和反向传播算法介绍的全连接网络的权重梯度计算算法只要知道了任意一个时刻的误差项以及上一个时刻循环层的输出值就可以按照下面的公式求出权重矩阵在t时刻的梯度 式 在式5中表示t时刻误差项向量的第i个分量表示t-1时刻循环层第i个神经元的输出值。 我们下面可以简单推导一下式5。 我们知道 因为对W求导与无关我们不再考虑。现在我们考虑对权重项求导。通过观察上式我们可以看到只与有关所以 按照上面的规律就可以生成式5里面的矩阵。 我们已经求得了权重矩阵W在t时刻的梯度最终的梯度是各个时刻的梯度之和 式 式6就是计算循环层权重矩阵W的梯度的公式。 ----------数学公式超高能预警---------- 前面已经介绍了的计算方法看上去还是比较直观的。然而读者也许会困惑为什么最终的梯度是各个时刻的梯度之和呢我们前面只是直接用了这个结论实际上这里面是有道理的只是这个数学推导比较绕脑子。感兴趣的同学可以仔细阅读接下来这一段它用到了矩阵对矩阵求导、张量与向量相乘运算的一些法则。 我们还是从这个式子开始 因为与W完全无关我们把它看做常量。现在考虑第一个式子加号右边的部分因为W和都是W的函数因此我们要用到大学里面都学过的导数乘法运算 因此上面第一个式子写成 我们最终需要计算的是 式 我们先计算式7加号左边的部分。是矩阵对矩阵求导其结果是一个四维张量(tensor)如下所示 接下来我们知道它是一个列向量。我们让上面的四维张量与这个向量相乘得到了一个三维张量再左乘行向量最终得到一个矩阵 接下来我们计算式7加号右边的部分 于是我们得到了如下递推公式 这样我们就证明了最终的梯度是各个时刻的梯度之和。 ----------数学公式超高能预警解除---------- 同权重矩阵W类似我们可以得到权重矩阵U的计算方法。 式 式8是误差函数在t时刻对权重矩阵U的梯度。和权重矩阵W一样最终的梯度也是各个时刻的梯度之和 具体的证明这里就不再赘述了感兴趣的读者可以练习推导一下。 RNN的梯度爆炸和消失问题 不幸的是实践中前面介绍的几种RNNs并不能很好的处理较长的序列。一个主要的原因是RNN在训练中很容易发生梯度爆炸和梯度消失这导致训练时梯度不能在较长序列中一直传递下去从而使RNN无法捕捉到长距离的影响。 为什么RNN会产生梯度爆炸和消失问题呢我们接下来将详细分析一下原因。我们根据式3可得 上式的定义为矩阵的模的上界。因为上式是一个指数函数如果t-k很大的话也就是向前看很远的时候会导致对应的误差项的值增长或缩小的非常快这样就会导致相应的梯度爆炸和梯度消失问题取决于大于1还是小于1。 通常来说梯度爆炸更容易处理一些。因为梯度爆炸的时候我们的程序会收到NaN错误。我们也可以设置一个梯度阈值当梯度超过这个阈值的时候可以直接截取。 梯度消失更难检测而且也更难处理一些。总的来说我们有三种方法应对梯度消失问题 合理的初始化权重值。初始化权重使每个神经元尽可能不要取极大或极小值以躲开梯度消失的区域。使用relu代替sigmoid和tanh作为激活函数。原理请参考上一篇文章零基础入门深度学习(4) - 卷积神经网络的激活函数一节。使用其他结构的RNNs比如长短时记忆网络LTSM和Gated Recurrent UnitGRU这是最流行的做法。我们将在以后的文章中介绍这两种网络。 RNN的应用举例——基于RNN的语言模型 现在我们介绍一下基于RNN语言模型。我们首先把词依次输入到循环神经网络中每输入一个词循环神经网络就输出截止到目前为止下一个最可能的词。例如当我们依次输入 我 昨天 上学 迟到 了 神经网络的输出如下图所示 其中s和e是两个特殊的词分别表示一个序列的开始和结束。 向量化 我们知道神经网络的输入和输出都是向量为了让语言模型能够被神经网络处理我们必须把词表达为向量的形式这样神经网络才能处理它。 神经网络的输入是词我们可以用下面的步骤对输入进行向量化 建立一个包含所有词的词典每个词在词典里面有一个唯一的编号。任意一个词都可以用一个N维的one-hot向量来表示。其中N是词典中包含的词的个数。假设一个词在词典中的编号是iv是表示这个词的向量是向量的第j个元素则 上面这个公式的含义可以用下面的图来直观的表示 使用这种向量化方法我们就得到了一个高维、稀疏的向量稀疏是指绝大部分元素的值都是0。处理这样的向量会导致我们的神经网络有很多的参数带来庞大的计算量。因此往往会需要使用一些降维方法将高维的稀疏向量转变为低维的稠密向量。不过这个话题我们就不再这篇文章中讨论了。 语言模型要求的输出是下一个最可能的词我们可以让循环神经网络计算计算词典中每个词是下一个词的概率这样概率最大的词就是下一个最可能的词。因此神经网络的输出向量也是一个N维向量向量中的每个元素对应着词典中相应的词是下一个词的概率。如下图所示 Softmax层 前面提到语言模型是对下一个词出现的概率进行建模。那么怎样让神经网络输出概率呢方法就是用softmax层作为神经网络的输出层。 我们先来看一下softmax函数的定义 这个公式看起来可能很晕我们举一个例子。Softmax层如下图所示 从上图我们可以看到softmax layer的输入是一个向量输出也是一个向量两个向量的维度是一样的在这个例子里面是4。输入向量x[1 2 3 4]经过softmax层之后经过上面的softmax函数计算转变为输出向量y[0.03 0.09 0.24 0.64]。计算过程为 我们来看看输出向量y的特征 每一项为取值为0-1之间的正数所有项的总和是1。我们不难发现这些特征和概率的特征是一样的因此我们可以把它们看做是概率。对于语言模型来说我们可以认为模型预测下一个词是词典中第一个词的概率是0.03是词典中第二个词的概率是0.09以此类推。 语言模型的训练 可以使用监督学习的方法对语言模型进行训练首先需要准备训练数据集。接下来我们介绍怎样把语料 我 昨天 上学 迟到 了 转换成语言模型的训练数据集。 首先我们获取输入-标签对 输入标签s我我昨天昨天上学上学迟到迟到了了e然后使用前面介绍过的向量化方法对输入x和标签y进行向量化。这里面有意思的是对标签y进行向量化其结果也是一个one-hot向量。例如我们对标签『我』进行向量化得到的向量中只有第2019个元素的值是1其他位置的元素的值都是0。它的含义就是下一个词是『我』的概率是1是其它词的概率都是0。 最后我们使用交叉熵误差函数作为优化目标对模型进行优化。 在实际工程中我们可以使用大量的语料来对模型进行训练获取训练数据和训练的方法都是相同的。 交叉熵误差 一般来说当神经网络的输出层是softmax层时对应的误差函数E通常选择交叉熵误差函数其定义如下 在上式中N是训练样本的个数向量是样本的标记向量是网络的输出。标记是一个one-hot向量例如如果网络的输出那么交叉熵误差是假设只有一个训练样本即N1 我们当然可以选择其他函数作为我们的误差函数比如最小平方误差函数(MSE)。不过对概率进行建模时选择交叉熵误差函数更make sense。具体原因感兴趣的读者请阅读参考文献7。 RNN的实现 完整代码请参考GitHub: https://github.com/hanbt/learn_dl/blob/master/rnn.py (python2.7) 为了加深我们对前面介绍的知识的理解我们来动手实现一个RNN层。我们复用了上一篇文章零基础入门深度学习(4) - 卷积神经网络中的一些代码所以先把它们导入进来。 import numpy as npfrom cnn import ReluActivator, IdentityActivator, element_wise_op我们用RecurrentLayer类来实现一个循环层。下面的代码是初始化一个循环层可以在构造函数中设置卷积层的超参数。我们注意到循环层有两个权重数组U和W。 class RecurrentLayer(object): def __init__(self, input_width, state_width, activator, learning_rate): self.input_width input_width self.state_width state_width self.activator activator self.learning_rate learning_rate self.times 0 # 当前时刻初始化为t0 self.state_list [] # 保存各个时刻的state self.state_list.append(np.zeros( (state_width, 1))) # 初始化s0 self.U np.random.uniform(-1e-4, 1e-4, (state_width, input_width)) # 初始化U self.W np.random.uniform(-1e-4, 1e-4, (state_width, state_width)) # 初始化W在forward方法中实现循环层的前向计算这部分比较简单。 def forward(self, input_array): 根据『式2』进行前向计算 self.times 1 state (np.dot(self.U, input_array) np.dot(self.W, self.state_list[-1])) element_wise_op(state, self.activator.forward) self.state_list.append(state)在backword方法中实现BPTT算法。 def backward(self, sensitivity_array, activator): 实现BPTT算法 self.calc_delta(sensitivity_array, activator) self.calc_gradient() def calc_delta(self, sensitivity_array, activator): self.delta_list [] # 用来保存各个时刻的误差项 for i in range(self.times): self.delta_list.append(np.zeros( (self.state_width, 1))) self.delta_list.append(sensitivity_array) # 迭代计算每个时刻的误差项 for k in range(self.times - 1, 0, -1): self.calc_delta_k(k, activator) def calc_delta_k(self, k, activator): 根据k1时刻的delta计算k时刻的delta state self.state_list[k1].copy() element_wise_op(self.state_list[k1], activator.backward) self.delta_list[k] np.dot( np.dot(self.delta_list[k1].T, self.W), np.diag(state[:,0])).T def calc_gradient(self): self.gradient_list [] # 保存各个时刻的权重梯度 for t in range(self.times 1): self.gradient_list.append(np.zeros( (self.state_width, self.state_width))) for t in range(self.times, 0, -1): self.calc_gradient_t(t) # 实际的梯度是各个时刻梯度之和 self.gradient reduce( lambda a, b: a b, self.gradient_list, self.gradient_list[0]) # [0]被初始化为0且没有被修改过 def calc_gradient_t(self, t): 计算每个时刻t权重的梯度 gradient np.dot(self.delta_list[t], self.state_list[t-1].T) self.gradient_list[t] gradient有意思的是BPTT算法虽然数学推导的过程很麻烦但是写成代码却并不复杂。 在update方法中实现梯度下降算法。 def update(self): 按照梯度下降更新权重 self.W - self.learning_rate * self.gradient上面的代码不包含权重U的更新。这部分实际上和全连接神经网络是一样的留给感兴趣的读者自己来完成吧。 循环层是一个带状态的层每次forword都会改变循环层的内部状态这给梯度检查带来了麻烦。因此我们需要一个reset_state方法来重置循环层的内部状态。 def reset_state(self): self.times 0 # 当前时刻初始化为t0 self.state_list [] # 保存各个时刻的state self.state_list.append(np.zeros( (self.state_width, 1))) # 初始化s0最后是梯度检查的代码。 def gradient_check(): 梯度检查 # 设计一个误差函数取所有节点输出项之和 error_function lambda o: o.sum() rl RecurrentLayer(3, 2, IdentityActivator(), 1e-3) # 计算forward值 x, d data_set() rl.forward(x[0]) rl.forward(x[1]) # 求取sensitivity map sensitivity_array np.ones(rl.state_list[-1].shape, dtypenp.float64) # 计算梯度 rl.backward(sensitivity_array, IdentityActivator()) # 检查梯度 epsilon 10e-4 for i in range(rl.W.shape[0]): for j in range(rl.W.shape[1]): rl.W[i,j] epsilon rl.reset_state() rl.forward(x[0]) rl.forward(x[1]) err1 error_function(rl.state_list[-1]) rl.W[i,j] - 2*epsilon rl.reset_state() rl.forward(x[0]) rl.forward(x[1]) err2 error_function(rl.state_list[-1]) expect_grad (err1 - err2) / (2 * epsilon) rl.W[i,j] epsilon print weights(%d,%d): expected - actural %f - %f % ( i, j, expect_grad, rl.gradient[i,j])需要注意每次计算error之前都要调用reset_state方法重置循环层的内部状态。下面是梯度检查的结果没问题 小节 至此我们讲完了基本的循环神经网络、它的训练算法BPTT以及在语言模型上的应用。RNN比较烧脑相信拿下前几篇文章的读者们搞定这篇文章也不在话下吧然而循环神经网络这个话题并没有完结。我们在前面说到过基本的循环神经网络存在梯度爆炸和梯度消失问题并不能真正的处理好长距离的依赖虽然有一些技巧可以减轻这些问题。事实上真正得到广泛的应用的是循环神经网络的一个变体长短时记忆网络。它内部有一些特殊的结构可以很好的处理长距离的依赖我们将在下一篇文章中详细的介绍它。现在让我们稍事休息准备挑战更为烧脑的长短时记忆网络吧。 参考资料 RECURRENT NEURAL NETWORKS TUTORIALUnderstanding LSTM NetworksThe Unreasonable Effectiveness of Recurrent Neural NetworksAttention and Augmented Recurrent Neural NetworksOn the difficulty of training recurrent neural networks, Bengio et al.Recurrent neural network based language model, Mikolov et al.Neural Network Classification, Categorical Data, Softmax Activation, and Cross Entropy Error, McCaffrey转载自https://zybuluo.com/hanbingtao/note/541458转载于:https://www.cnblogs.com/tu6ge/p/9492338.html
