(转)零基础入门深度学习(6) - 长短时记忆网络(LSTM)

不管即将到来的是大数据时代仍是人工智能时代，亦或是传统行业使用人工智能在云上处理大数据的时代，做为一个有理想有追求的程序员，不懂深度学习（Deep Learning）这个超热的技术，会不会感受立刻就out了？如今救命稻草来了，《零基础入门深度学习》系列文章旨在讲帮助爱编程的你从零基础达到入门级水平。零基础意味着你不须要太多的数学知识，只要会写程序就好了，没错，这是专门为程序员写的文章。虽然文中会有不少公式你也许看不懂，但同时也会有更多的代码，程序员的你必定能看懂的（我周围是一群狂热的Clean Code程序员，因此我写的代码也不会不好）。

 

文章列表

零基础入门深度学习(1) - 感知器 
零基础入门深度学习(2) - 线性单元和梯度降低 
零基础入门深度学习(3) - 神经网络和反向传播算法 
零基础入门深度学习(4) - 卷积神经网络 
零基础入门深度学习(5) - 循环神经网络 
零基础入门深度学习(6) - 长短时记忆网络(LSTM) 
零基础入门深度学习(7) - 递归神经网络

 

往期回顾

在上一篇文章中，咱们介绍了循环神经网络以及它的训练算法。咱们也介绍了循环神经网络很难训练的缘由，这致使了它在实际应用中，很难处理长距离的依赖。在本文中，咱们将介绍一种改进以后的循环神经网络：长短时记忆网络(Long Short Term Memory Network, LSTM)，它成功的解决了原始循环神经网络的缺陷，成为当前最流行的RNN，在语音识别、图片描述、天然语言处理等许多领域中成功应用。但不幸的一面是，LSTM的结构很复杂，所以，咱们须要花上一些力气，才能把LSTM以及它的训练算法弄明白。在搞清楚LSTM以后，咱们再介绍一种LSTM的变体：GRU (Gated Recurrent Unit)。 它的结构比LSTM简单，而效果却和LSTM同样好，所以，它正在逐渐流行起来。最后，咱们仍然会动手实现一个LSTM。

 

长短时记忆网络是啥

咱们首先了解一下长短时记忆网络产生的背景。回顾一下零基础入门深度学习(5) - 循环神经网络中推导的，偏差项沿时间反向传播的公式：

 

 

 

 

咱们能够根据下面的不等式，来获取的模的上界（模能够看作对中每一项值的大小的度量）：

 

 

 

 

咱们能够看到，偏差项从t时刻传递到k时刻，其值的上界是的指数函数。分别是对角矩阵和矩阵W模的上界。显然，除非乘积的值位于1附近，不然，当t-k很大时（也就是偏差传递不少个时刻时），整个式子的值就会变得极小（当乘积小于1）或者极大（当乘积大于1），前者就是梯度消失，后者就是梯度爆炸。虽然科学家们搞出了不少技巧（好比怎样初始化权重），让的值尽量贴近于1，终究仍是难以抵挡指数函数的威力。

梯度消失到底意味着什么？在零基础入门深度学习(5) - 循环神经网络中咱们已证实，权重数组W最终的梯度是各个时刻的梯度之和，即：

 

 

 

 

假设某轮训练中，各时刻的梯度以及最终的梯度之和以下图：

咱们就能够看到，从上图的t-3时刻开始，梯度已经几乎减小到0了。那么，从这个时刻开始再往以前走，获得的梯度（几乎为零）就不会对最终的梯度值有任何贡献，这就至关于不管t-3时刻以前的网络状态h是什么，在训练中都不会对权重数组W的更新产生影响，也就是网络事实上已经忽略了t-3时刻以前的状态。这就是原始RNN没法处理长距离依赖的缘由。

既然找到了问题的缘由，那么咱们就能解决它。从问题的定位到解决，科学家们大概花了七、8年时间。终于有一天，Hochreiter和Schmidhuber两位科学家发明出长短时记忆网络，一举解决这个问题。

其实，长短时记忆网络的思路比较简单。原始RNN的隐藏层只有一个状态，即h，它对于短时间的输入很是敏感。那么，假如咱们再增长一个状态，即c，让它来保存长期的状态，那么问题不就解决了么？以下图所示：

新增长的状态c，称为单元状态(cell state)。咱们把上图按照时间维度展开：

上图仅仅是一个示意图，咱们能够看出，在t时刻，LSTM的输入有三个：当前时刻网络的输入值、上一时刻LSTM的输出值、以及上一时刻的单元状态；LSTM的输出有两个：当前时刻LSTM输出值、和当前时刻的单元状态。注意、、都是向量。

LSTM的关键，就是怎样控制长期状态c。在这里，LSTM的思路是使用三个控制开关。第一个开关，负责控制继续保存长期状态c；第二个开关，负责控制把即时状态输入到长期状态c；第三个开关，负责控制是否把长期状态c做为当前的LSTM的输出。三个开关的做用以下图所示：

接下来，咱们要描述一下，输出h和单元状态c的具体计算方法。

 

长短时记忆网络的前向计算

前面描述的开关是怎样在算法中实现的呢？这就用到了门（gate）的概念。门实际上就是一层全链接层，它的输入是一个向量，输出是一个0到1之间的实数向量。假设W是门的权重向量，是偏置项，那么门能够表示为：

 

 

 

 

门的使用，就是用门的输出向量按元素乘以咱们须要控制的那个向量。由于门的输出是0到1之间的实数向量，那么，当门输出为0时，任何向量与之相乘都会获得0向量，这就至关于啥都不能经过；输出为1时，任何向量与之相乘都不会有任何改变，这就至关于啥均可以经过。由于（也就是sigmoid函数）的值域是(0,1)，因此门的状态都是半开半闭的。

LSTM用两个门来控制单元状态c的内容，一个是遗忘门（forget gate），它决定了上一时刻的单元状态有多少保留到当前时刻；另外一个是输入门（input gate），它决定了当前时刻网络的输入有多少保存到单元状态。LSTM用输出门（output gate）来控制单元状态有多少输出到LSTM的当前输出值。

咱们先来看一下遗忘门：

 

 

式

 

上式中，是遗忘门的权重矩阵，表示把两个向量链接成一个更长的向量，是遗忘门的偏置项，是sigmoid函数。若是输入的维度是，隐藏层的维度是，单元状态的维度是（一般），则遗忘门的权重矩阵维度是。事实上，权重矩阵都是两个矩阵拼接而成的：一个是，它对应着输入项，其维度为；一个是，它对应着输入项，其维度为。能够写为：

 

 

 

 

下图显示了遗忘门的计算：

接下来看看输入门：

 

 

式

 

上式中，是输入门的权重矩阵，是输入门的偏置项。下图表示了输入门的计算：

接下来，咱们计算用于描述当前输入的单元状态，它是根据上一次的输出和本次输入来计算的：

 

 

式

 

下图是的计算：

如今，咱们计算当前时刻的单元状态。它是由上一次的单元状态按元素乘以遗忘门，再用当前输入的单元状态按元素乘以输入门，再将两个积加和产生的：

 

 

式

 

符号表示按元素乘。下图是的计算：

这样，咱们就把LSTM关于当前的记忆和长期的记忆组合在一块儿，造成了新的单元状态。因为遗忘门的控制，它能够保存好久好久以前的信息，因为输入门的控制，它又能够避免当前可有可无的内容进入记忆。下面，咱们要看看输出门，它控制了长期记忆对当前输出的影响：

 

 

式

 

下图表示输出门的计算：

LSTM最终的输出，是由输出门和单元状态共同肯定的：

 

 

式

 

下图表示LSTM最终输出的计算：

式1到式6就是LSTM前向计算的所有公式。至此，咱们就把LSTM前向计算讲完了。

 

长短时记忆网络的训练

熟悉咱们这个系列文章的同窗都清楚，训练部分每每比前向计算部分复杂多了。LSTM的前向计算都这么复杂，那么，可想而知，它的训练算法必定是很是很是复杂的。如今只有作几回深呼吸，再一头扎进公式海洋吧。

 

LSTM训练算法框架

LSTM的训练算法仍然是反向传播算法，对于这个算法，咱们已经很是熟悉了。主要有下面三个步骤：

1. 前向计算每一个神经元的输出值，对于LSTM来讲，即、、、、五个向量的值。计算方法已经在上一节中描述过了。
2. 反向计算每一个神经元的偏差项值。与循环神经网络同样，LSTM偏差项的反向传播也是包括两个方向：一个是沿时间的反向传播，即从当前t时刻开始，计算每一个时刻的偏差项；一个是将偏差项向上一层传播。
3. 根据相应的偏差项，计算每一个权重的梯度。

 

关于公式和符号的说明

首先，咱们对推导中用到的一些公式、符号作一下必要的说明。

接下来的推导中，咱们设定gate的激活函数为sigmoid函数，输出的激活函数为tanh函数。他们的导数分别为：

 

 

 

 

从上面能够看出，sigmoid和tanh函数的导数都是原函数的函数。这样，咱们一旦计算原函数的值，就能够用它来计算出导数的值。

LSTM须要学习的参数共有8组，分别是：遗忘门的权重矩阵和偏置项、输入门的权重矩阵和偏置项、输出门的权重矩阵和偏置项，以及计算单元状态的权重矩阵和偏置项。由于权重矩阵的两部分在反向传播中使用不一样的公式，所以在后续的推导中，权重矩阵、、、都将被写为分开的两个矩阵：、、、、、、、。

咱们解释一下按元素乘符号。当做用于两个向量时，运算以下：

 

 

 

 

当做用于一个向量和一个矩阵时，运算以下：

 

 

 

 

当做用于两个矩阵时，两个矩阵对应位置的元素相乘。按元素乘能够在某些状况下简化矩阵和向量运算。例如，当一个对角矩阵右乘一个矩阵时，至关于用对角矩阵的对角线组成的向量按元素乘那个矩阵：

 

 

 

 

当一个行向量右乘一个对角矩阵时，至关于这个行向量按元素乘那个矩阵对角线组成的向量：

 

 

 

 

上面这两点，在咱们后续推导中会屡次用到。

在t时刻，LSTM的输出值为。咱们定义t时刻的偏差项为：

 

 

 

 

注意，和前面几篇文章不一样，咱们这里假设偏差项是损失函数对输出值的导数，而不是对加权输入的导数。由于LSTM有四个加权输入，分别对应、、、，咱们但愿往上一层传递一个偏差项而不是四个。但咱们仍然须要定义出这四个加权输入，以及他们对应的偏差项。

 

 

 

 

 

偏差项沿时间的反向传递

沿时间反向传递偏差项，就是要计算出t-1时刻的偏差项。

 

 

 

 

咱们知道，是一个Jacobian矩阵。若是隐藏层h的维度是N的话，那么它就是一个矩阵。为了求出它，咱们列出的计算公式，即前面的式6和式4：

 

 

 

 

显然，、、、都是的函数，那么，利用全导数公式可得：

 

 

式

 

下面，咱们要把式7中的每一个偏导数都求出来。根据式6，咱们能够求出：

 

 

 

 

根据式4，咱们能够求出：

 

 

 

 

由于：

 

 

 

 

咱们很容易得出：

 

 

 

 

将上述偏导数带入到式7，咱们获得：

 

 

式

 

根据、、、的定义，可知：

 

 

式式式式

 

式8到式12就是将偏差沿时间反向传播一个时刻的公式。有了它，咱们能够写出将偏差项向前传递到任意k时刻的公式：

 

 

式

 

 

将偏差项传递到上一层

咱们假设当前为第l层，定义l-1层的偏差项是偏差函数对l-1层加权输入的导数，即：

 

 

 

 

本次LSTM的输入由下面的公式计算：

 

 

 

 

上式中，表示第l-1层的激活函数。

由于、、、都是的函数，又是的函数，所以，要求出E对的导数，就须要使用全导数公式：

 

 

式

 

式14就是将偏差传递到上一层的公式。

 

权重梯度的计算

对于、、、的权重梯度，咱们知道它的梯度是各个时刻梯度之和（证实过程请参考文章零基础入门深度学习(5) - 循环神经网络），咱们首先求出它们在t时刻的梯度，而后再求出他们最终的梯度。

咱们已经求得了偏差项、、、，很容易求出t时刻的、的、的、的：

 

 

 

 

将各个时刻的梯度加在一块儿，就能获得最终的梯度：

 

 

 

 

对于偏置项、、、的梯度，也是将各个时刻的梯度加在一块儿。下面是各个时刻的偏置项梯度：

 

 

 

 

下面是最终的偏置项梯度，即将各个时刻的偏置项梯度加在一块儿：

 

 

 

 

对于、、、的权重梯度，只须要根据相应的偏差项直接计算便可：

 

 

 

 

以上就是LSTM的训练算法的所有公式。由于这里面存在不少重复的模式，仔细看看，会发觉并非太复杂。

固然，LSTM存在着至关多的变体，读者能够在互联网上找到不少资料。由于你们已经熟悉了基本LSTM的算法，所以理解这些变体比较容易，所以本文就再也不赘述了。

 

长短时记忆网络的实现

在下面的实现中，LSTMLayer的参数包括输入维度、输出维度、隐藏层维度，单元状态维度等于隐藏层维度。gate的激活函数为sigmoid函数，输出的激活函数为tanh。

 

激活函数的实现

咱们先实现两个激活函数：sigmoid和tanh。

 

1. class SigmoidActivator(object):
2. def forward(self, weighted_input):
3. return 1.0 / (1.0 + np.exp(-weighted_input))
5. def backward(self, output):
6. return output * (1 - output)
9. class TanhActivator(object):
10. def forward(self, weighted_input):
11. return 2.0 / (1.0 + np.exp(-2 * weighted_input)) - 1.0
13. def backward(self, output):
14. return 1 - output * output

 

LSTM初始化

和前两篇文章代码架构同样，咱们把LSTM的实现放在LstmLayer类中。

根据LSTM前向计算和方向传播算法，咱们须要初始化一系列矩阵和向量。这些矩阵和向量有两类用途，一类是用于保存模型参数，例如、、、、、、、；另外一类是保存各类中间计算结果，以便于反向传播算法使用，它们包括、、、、、、、、、、，以及各个权重对应的梯度。

在构造函数的初始化中，只初始化了与forward计算相关的变量，与backward相关的变量没有初始化。这是由于构造LSTM对象的时候，咱们还不知道它将来是用于训练（既有forward又有backward）仍是推理（只有forward）。

 

1. class LstmLayer(object):
2. def __init__(self, input_width, state_width,
3. learning_rate):
4. self.input_width = input_width
5. self.state_width = state_width
6. self.learning_rate = learning_rate
7. # 门的激活函数
8. self.gate_activator = SigmoidActivator()
9. # 输出的激活函数
10. self.output_activator = TanhActivator()
11. # 当前时刻初始化为t0
12. self.times = 0
13. # 各个时刻的单元状态向量c
14. self.c_list = self.init_state_vec()
15. # 各个时刻的输出向量h
16. self.h_list = self.init_state_vec()
17. # 各个时刻的遗忘门f
18. self.f_list = self.init_state_vec()
19. # 各个时刻的输入门i
20. self.i_list = self.init_state_vec()
21. # 各个时刻的输出门o
22. self.o_list = self.init_state_vec()
23. # 各个时刻的即时状态c~
24. self.ct_list = self.init_state_vec()
25. # 遗忘门权重矩阵Wfh, Wfx, 偏置项bf
26. self.Wfh, self.Wfx, self.bf = (
27. self.init_weight_mat())
28. # 输入门权重矩阵Wfh, Wfx, 偏置项bf
29. self.Wih, self.Wix, self.bi = (
30. self.init_weight_mat())
31. # 输出门权重矩阵Wfh, Wfx, 偏置项bf
32. self.Woh, self.Wox, self.bo = (
33. self.init_weight_mat())
34. # 单元状态权重矩阵Wfh, Wfx, 偏置项bf
35. self.Wch, self.Wcx, self.bc = (
36. self.init_weight_mat())
38. def init_state_vec(self):
39. '''
40. 初始化保存状态的向量
41. '''
42. state_vec_list = []
43. state_vec_list.append(np.zeros(
44. (self.state_width, 1)))
45. return state_vec_list
47. def init_weight_mat(self):
48. '''
49. 初始化权重矩阵
50. '''
51. Wh = np.random.uniform(-1e-4, 1e-4,
52. (self.state_width, self.state_width))
53. Wx = np.random.uniform(-1e-4, 1e-4,
54. (self.state_width, self.input_width))
55. b = np.zeros((self.state_width, 1))
56. return Wh, Wx, b

 

前向计算的实现

forward方法实现了LSTM的前向计算：

 

1. def forward(self, x):
2. '''
3. 根据式1-式6进行前向计算
4. '''
5. self.times += 1
6. # 遗忘门
7. fg = self.calc_gate(x, self.Wfx, self.Wfh,
8. self.bf, self.gate_activator)
9. self.f_list.append(fg)
10. # 输入门
11. ig = self.calc_gate(x, self.Wix, self.Wih,
12. self.bi, self.gate_activator)
13. self.i_list.append(ig)
14. # 输出门
15. og = self.calc_gate(x, self.Wox, self.Woh,
16. self.bo, self.gate_activator)
17. self.o_list.append(og)
18. # 即时状态
19. ct = self.calc_gate(x, self.Wcx, self.Wch,
20. self.bc, self.output_activator)
21. self.ct_list.append(ct)
22. # 单元状态
23. c = fg * self.c_list[self.times - 1] + ig * ct
24. self.c_list.append(c)
25. # 输出
26. h = og * self.output_activator.forward(c)
27. self.h_list.append(h)
29. def calc_gate(self, x, Wx, Wh, b, activator):
30. '''
31. 计算门
32. '''
33. h = self.h_list[self.times - 1] # 上次的LSTM输出
34. net = np.dot(Wh, h) + np.dot(Wx, x) + b
35. gate = activator.forward(net)
36. return gate

从上面的代码咱们能够看到，门的计算都是相同的算法，而门和的计算仅仅是激活函数不一样。所以咱们提出了calc_gate方法，这样减小了不少重复代码。

 

反向传播算法的实现

backward方法实现了LSTM的反向传播算法。须要注意的是，与backword相关的内部状态变量是在调用backward方法以后才初始化的。这种延迟初始化的一个好处是，若是LSTM只是用来推理，那么就不须要初始化这些变量，节省了不少内存。

 

1. def backward(self, x, delta_h, activator):
2. '''
3. 实现LSTM训练算法
4. '''
5. self.calc_delta(delta_h, activator)
6. self.calc_gradient(x)

算法主要分红两个部分，一部分使计算偏差项：

 

1. def calc_delta(self, delta_h, activator):
2. # 初始化各个时刻的偏差项
3. self.delta_h_list = self.init_delta() # 输出偏差项
4. self.delta_o_list = self.init_delta() # 输出门偏差项
5. self.delta_i_list = self.init_delta() # 输入门偏差项
6. self.delta_f_list = self.init_delta() # 遗忘门偏差项
7. self.delta_ct_list = self.init_delta() # 即时输出偏差项
9. # 保存从上一层传递下来的当前时刻的偏差项
10. self.delta_h_list[-1] = delta_h
12. # 迭代计算每一个时刻的偏差项
13. for k in range(self.times, 0, -1):
14. self.calc_delta_k(k)
16. def init_delta(self):
17. '''
18. 初始化偏差项
19. '''
20. delta_list = []
21. for i in range(self.times + 1):
22. delta_list.append(np.zeros(
23. (self.state_width, 1)))
24. return delta_list
26. def calc_delta_k(self, k):
27. '''
28. 根据k时刻的delta_h，计算k时刻的delta_f、
29. delta_i、delta_o、delta_ct，以及k-1时刻的delta_h
30. '''
31. # 得到k时刻前向计算的值
32. ig = self.i_list[k]
33. og = self.o_list[k]
34. fg = self.f_list[k]
35. ct = self.ct_list[k]
36. c = self.c_list[k]
37. c_prev = self.c_list[k-1]
38. tanh_c = self.output_activator.forward(c)
39. delta_k = self.delta_h_list[k]
41. # 根据式9计算delta_o
42. delta_o = (delta_k * tanh_c *
43. self.gate_activator.backward(og))
44. delta_f = (delta_k * og *
45. (1 - tanh_c * tanh_c) * c_prev *
46. self.gate_activator.backward(fg))
47. delta_i = (delta_k * og *
48. (1 - tanh_c * tanh_c) * ct *
49. self.gate_activator.backward(ig))
50. delta_ct = (delta_k * og *
51. (1 - tanh_c * tanh_c) * ig *
52. self.output_activator.backward(ct))
53. delta_h_prev = (
54. np.dot(delta_o.transpose(), self.Woh) +
55. np.dot(delta_i.transpose(), self.Wih) +
56. np.dot(delta_f.transpose(), self.Wfh) +
57. np.dot(delta_ct.transpose(), self.Wch)
58. ).transpose()
60. # 保存所有delta值
61. self.delta_h_list[k-1] = delta_h_prev
62. self.delta_f_list[k] = delta_f
63. self.delta_i_list[k] = delta_i
64. self.delta_o_list[k] = delta_o
65. self.delta_ct_list[k] = delta_ct

另外一部分是计算梯度：

 

1. def calc_gradient(self, x):
2. # 初始化遗忘门权重梯度矩阵和偏置项
3. self.Wfh_grad, self.Wfx_grad, self.bf_grad = (
4. self.init_weight_gradient_mat())
5. # 初始化输入门权重梯度矩阵和偏置项
6. self.Wih_grad, self.Wix_grad, self.bi_grad = (
7. self.init_weight_gradient_mat())
8. # 初始化输出门权重梯度矩阵和偏置项
9. self.Woh_grad, self.Wox_grad, self.bo_grad = (
10. self.init_weight_gradient_mat())
11. # 初始化单元状态权重梯度矩阵和偏置项
12. self.Wch_grad, self.Wcx_grad, self.bc_grad = (
13. self.init_weight_gradient_mat())
15. # 计算对上一次输出h的权重梯度
16. for t in range(self.times, 0, -1):
17. # 计算各个时刻的梯度
18. (Wfh_grad, bf_grad,
19. Wih_grad, bi_grad,
20. Woh_grad, bo_grad,
21. Wch_grad, bc_grad) = (
22. self.calc_gradient_t(t))
23. # 实际梯度是各时刻梯度之和
24. self.Wfh_grad += Wfh_grad
25. self.bf_grad += bf_grad
26. self.Wih_grad += Wih_grad
27. self.bi_grad += bi_grad
28. self.Woh_grad += Woh_grad
29. self.bo_grad += bo_grad
30. self.Wch_grad += Wch_grad
31. self.bc_grad += bc_grad
32. print '-----%d-----' % t
33. print Wfh_grad
34. print self.Wfh_grad
36. # 计算对本次输入x的权重梯度
37. xt = x.transpose()
38. self.Wfx_grad = np.dot(self.delta_f_list[-1], xt)
39. self.Wix_grad = np.dot(self.delta_i_list[-1], xt)
40. self.Wox_grad = np.dot(self.delta_o_list[-1], xt)
41. self.Wcx_grad = np.dot(self.delta_ct_list[-1], xt)
43. def init_weight_gradient_mat(self):
44. '''
45. 初始化权重矩阵
46. '''
47. Wh_grad = np.zeros((self.state_width,
48. self.state_width))
49. Wx_grad = np.zeros((self.state_width,
50. self.input_width))
51. b_grad = np.zeros((self.state_width, 1))
52. return Wh_grad, Wx_grad, b_grad
54. def calc_gradient_t(self, t):
55. '''
56. 计算每一个时刻t权重的梯度
57. '''
58. h_prev = self.h_list[t-1].transpose()
59. Wfh_grad = np.dot(self.delta_f_list[t], h_prev)
60. bf_grad = self.delta_f_list[t]
61. Wih_grad = np.dot(self.delta_i_list[t], h_prev)
62. bi_grad = self.delta_f_list[t]
63. Woh_grad = np.dot(self.delta_o_list[t], h_prev)
64. bo_grad = self.delta_f_list[t]
65. Wch_grad = np.dot(self.delta_ct_list[t], h_prev)
66. bc_grad = self.delta_ct_list[t]
67. return Wfh_grad, bf_grad, Wih_grad, bi_grad, \
68. Woh_grad, bo_grad, Wch_grad, bc_grad

 

梯度降低算法的实现

下面是用梯度降低算法来更新权重：

 

1. def update(self):
2. '''
3. 按照梯度降低，更新权重
4. '''
5. self.Wfh -= self.learning_rate * self.Whf_grad
6. self.Wfx -= self.learning_rate * self.Whx_grad
7. self.bf -= self.learning_rate * self.bf_grad
8. self.Wih -= self.learning_rate * self.Whi_grad
9. self.Wix -= self.learning_rate * self.Whi_grad
10. self.bi -= self.learning_rate * self.bi_grad
11. self.Woh -= self.learning_rate * self.Wof_grad
12. self.Wox -= self.learning_rate * self.Wox_grad
13. self.bo -= self.learning_rate * self.bo_grad
14. self.Wch -= self.learning_rate * self.Wcf_grad
15. self.Wcx -= self.learning_rate * self.Wcx_grad
16. self.bc -= self.learning_rate * self.bc_grad

 

梯度检查的实现

和RecurrentLayer同样，为了支持梯度检查，咱们须要支持重置内部状态：

 

1. def reset_state(self):
2. # 当前时刻初始化为t0
3. self.times = 0
4. # 各个时刻的单元状态向量c
5. self.c_list = self.init_state_vec()
6. # 各个时刻的输出向量h
7. self.h_list = self.init_state_vec()
8. # 各个时刻的遗忘门f
9. self.f_list = self.init_state_vec()
10. # 各个时刻的输入门i
11. self.i_list = self.init_state_vec()
12. # 各个时刻的输出门o
13. self.o_list = self.init_state_vec()
14. # 各个时刻的即时状态c~
15. self.ct_list = self.init_state_vec()

最后，是梯度检查的代码：

 

1. def data_set():
2. x = [np.array([[1], [2], [3]]),
3. np.array([[2], [3], [4]])]
4. d = np.array([[1], [2]])
5. return x, d
7. def gradient_check():
8. '''
9. 梯度检查
10. '''
11. # 设计一个偏差函数，取全部节点输出项之和
12. error_function = lambda o: o.sum()
14. lstm = LstmLayer(3, 2, 1e-3)
16. # 计算forward值
17. x, d = data_set()
18. lstm.forward(x[0])
19. lstm.forward(x[1])
21. # 求取sensitivity map
22. sensitivity_array = np.ones(lstm.h_list[-1].shape,
23. dtype=np.float64)
24. # 计算梯度
25. lstm.backward(x[1], sensitivity_array, IdentityActivator())
27. # 检查梯度
28. epsilon = 10e-4
29. for i in range(lstm.Wfh.shape[0]):
30. for j in range(lstm.Wfh.shape[1]):
31. lstm.Wfh[i,j] += epsilon
32. lstm.reset_state()
33. lstm.forward(x[0])
34. lstm.forward(x[1])
35. err1 = error_function(lstm.h_list[-1])
36. lstm.Wfh[i,j] -= 2*epsilon
37. lstm.reset_state()
38. lstm.forward(x[0])
39. lstm.forward(x[1])
40. err2 = error_function(lstm.h_list[-1])
41. expect_grad = (err1 - err2) / (2 * epsilon)
42. lstm.Wfh[i,j] += epsilon
43. print 'weights(%d,%d): expected - actural %.4e - %.4e' % (
44. i, j, expect_grad, lstm.Wfh_grad[i,j])
45. return lstm

咱们只对作了检查，读者能够自行增长对其余梯度的检查。下面是某次梯度检查的结果：

 

GRU

前面咱们讲了一种普通的LSTM，事实上LSTM存在不少变体，许多论文中的LSTM都或多或少的不太同样。在众多的LSTM变体中，GRU (Gated Recurrent Unit)也许是最成功的一种。它对LSTM作了不少简化，同时却保持着和LSTM相同的效果。所以，GRU最近变得愈来愈流行。

GRU对LSTM作了两个大改动：

1. 将输入门、遗忘门、输出门变为两个门：更新门（Update Gate）和重置门（Reset Gate）。
2. 将单元状态与输出合并为一个状态：。

GRU的前向计算公式为：

 

 

 

 

下图是GRU的示意图：

GRU的训练算法比LSTM简单一些，留给读者自行推导，本文就再也不赘述了。

 

小结

至此，LSTM——也许是结构最复杂的一类神经网络——就讲完了，相信拿下前几篇文章的读者们搞定这篇文章也不在话下吧！如今咱们已经了解循环神经网络和它最流行的变体——LSTM，它们均可以用来处理序列。可是，有时候仅仅拥有处理序列的能力还不够，还须要处理比序列更为复杂的结构（好比树结构），这时候就须要用到另一类网络：递归神经网络(Recursive Neural Network)，巧合的是，它的缩写也是RNN。在下一篇文章中，咱们将介绍递归神经网络和它的训练算法。如今，漫长的烧脑暂告一段落，休息一下吧:)