|
| 1 | +# coding: utf-8 |
| 2 | +# 作者: Robert Guthrie |
| 3 | + |
| 4 | +import torch |
| 5 | +import torch.autograd as autograd |
| 6 | +import torch.nn as nn |
| 7 | +import torch.optim as optim |
| 8 | + |
| 9 | +torch.manual_seed(1) |
| 10 | + |
| 11 | + |
| 12 | +def to_scalar(var): |
| 13 | + # 返回 python 浮点数 (float) |
| 14 | + return var.view(-1).data.tolist()[0] |
| 15 | + |
| 16 | + |
| 17 | +def argmax(vec): |
| 18 | + # 以 python 整数的形式返回 argmax |
| 19 | + _, idx = torch.max(vec, 1) |
| 20 | + return to_scalar(idx) |
| 21 | + |
| 22 | + |
| 23 | +def prepare_sequence(seq, to_ix): |
| 24 | + idxs = [to_ix[w] for w in seq] |
| 25 | + tensor = torch.LongTensor(idxs) |
| 26 | + return autograd.Variable(tensor) |
| 27 | + |
| 28 | + |
| 29 | +# 使用数值上稳定的方法为前向算法计算指数和的对数 |
| 30 | +def log_sum_exp(vec): |
| 31 | + max_score = vec[0, argmax(vec)] |
| 32 | + max_score_broadcast = max_score.view(1, -1).expand(1, vec.size()[1]) |
| 33 | + return max_score + \ |
| 34 | + torch.log(torch.sum(torch.exp(vec - max_score_broadcast))) |
| 35 | + |
| 36 | + |
| 37 | +class BiLSTM_CRF(nn.Module): |
| 38 | + def __init__(self, vocab_size, tag_to_ix, embedding_dim, hidden_dim): |
| 39 | + super(BiLSTM_CRF, self).__init__() |
| 40 | + self.embedding_dim = embedding_dim |
| 41 | + self.hidden_dim = hidden_dim |
| 42 | + self.vocab_size = vocab_size |
| 43 | + self.tag_to_ix = tag_to_ix |
| 44 | + self.tagset_size = len(tag_to_ix) |
| 45 | + |
| 46 | + self.word_embeds = nn.Embedding(vocab_size, embedding_dim) |
| 47 | + self.lstm = nn.LSTM( |
| 48 | + embedding_dim, hidden_dim // 2, num_layers=1, bidirectional=True) |
| 49 | + |
| 50 | + # 将LSTM的输出映射到标记空间 |
| 51 | + self.hidden2tag = nn.Linear(hidden_dim, self.tagset_size) |
| 52 | + |
| 53 | + # 过渡参数矩阵. 条目 i,j 是 |
| 54 | + # *从 * j *到* i 的过渡的分数 |
| 55 | + self.transitions = nn.Parameter( |
| 56 | + torch.randn(self.tagset_size, self.tagset_size)) |
| 57 | + |
| 58 | + # 这两句声明强制约束了我们不能 |
| 59 | + # 向开始标记标注传递和从结束标注传递 |
| 60 | + self.transitions.data[tag_to_ix[START_TAG], :] = -10000 |
| 61 | + self.transitions.data[:, tag_to_ix[STOP_TAG]] = -10000 |
| 62 | + |
| 63 | + self.hidden = self.init_hidden() |
| 64 | + |
| 65 | + def init_hidden(self): |
| 66 | + return (autograd.Variable(torch.randn(2, 1, self.hidden_dim // 2)), |
| 67 | + autograd.Variable(torch.randn(2, 1, self.hidden_dim // 2))) |
| 68 | + |
| 69 | + def _forward_alg(self, feats): |
| 70 | + # 执行前向算法来计算分割函数 |
| 71 | + init_alphas = torch.Tensor(1, self.tagset_size).fill_(-10000.) |
| 72 | + # START_TAG 包含所有的分数 |
| 73 | + init_alphas[0][self.tag_to_ix[START_TAG]] = 0. |
| 74 | + |
| 75 | + # 将其包在一个变量类型中继而得到自动的反向传播 |
| 76 | + forward_var = autograd.Variable(init_alphas) |
| 77 | + |
| 78 | + # 在句子中迭代 |
| 79 | + for feat in feats: |
| 80 | + alphas_t = [] # 在这个时间步的前向变量 |
| 81 | + for next_tag in range(self.tagset_size): |
| 82 | + # 对 emission 得分执行广播机制: 它总是相同的, |
| 83 | + # 不论前一个标注如何 |
| 84 | + emit_score = feat[next_tag].view( |
| 85 | + 1, -1).expand(1, self.tagset_size) |
| 86 | + # trans_score 第 i 个条目是 |
| 87 | + # 从i过渡到 next_tag 的分数 |
| 88 | + trans_score = self.transitions[next_tag].view(1, -1) |
| 89 | + # next_tag_var 第 i 个条目是在我们执行 对数-求和-指数 前 |
| 90 | + # 边缘的值 (i -> next_tag) |
| 91 | + next_tag_var = forward_var + trans_score + emit_score |
| 92 | + # 这个标注的前向变量是 |
| 93 | + # 对所有的分数执行 对数-求和-指数 |
| 94 | + alphas_t.append(log_sum_exp(next_tag_var)) |
| 95 | + forward_var = torch.cat(alphas_t).view(1, -1) |
| 96 | + terminal_var = forward_var + self.transitions[self.tag_to_ix[STOP_TAG]] |
| 97 | + alpha = log_sum_exp(terminal_var) |
| 98 | + return alpha |
| 99 | + |
| 100 | + def _get_lstm_features(self, sentence): |
| 101 | + self.hidden = self.init_hidden() |
| 102 | + embeds = self.word_embeds(sentence).view(len(sentence), 1, -1) |
| 103 | + lstm_out, self.hidden = self.lstm(embeds, self.hidden) |
| 104 | + lstm_out = lstm_out.view(len(sentence), self.hidden_dim) |
| 105 | + lstm_feats = self.hidden2tag(lstm_out) |
| 106 | + return lstm_feats |
| 107 | + |
| 108 | + def _score_sentence(self, feats, tags): |
| 109 | + # 给出标记序列的分数 |
| 110 | + score = autograd.Variable(torch.Tensor([0])) |
| 111 | + tags = torch.cat([torch.LongTensor([self.tag_to_ix[START_TAG]]), tags]) |
| 112 | + for i, feat in enumerate(feats): |
| 113 | + score = score + \ |
| 114 | + self.transitions[tags[i + 1], tags[i]] + feat[tags[i + 1]] |
| 115 | + score = score + self.transitions[self.tag_to_ix[STOP_TAG], tags[-1]] |
| 116 | + return score |
| 117 | + |
| 118 | + def _viterbi_decode(self, feats): |
| 119 | + backpointers = [] |
| 120 | + |
| 121 | + # 在对数空间中初始化维特比变量 |
| 122 | + init_vvars = torch.Tensor(1, self.tagset_size).fill_(-10000.) |
| 123 | + init_vvars[0][self.tag_to_ix[START_TAG]] = 0 |
| 124 | + |
| 125 | + # 在第 i 步的 forward_var 存放第 i-1 步的维特比变量 |
| 126 | + forward_var = autograd.Variable(init_vvars) |
| 127 | + for feat in feats: |
| 128 | + bptrs_t = [] # 存放这一步的后指针 |
| 129 | + viterbivars_t = [] # 存放这一步的维特比变量 |
| 130 | + |
| 131 | + for next_tag in range(self.tagset_size): |
| 132 | + # next_tag_var[i] 存放先前一步标注i的 |
| 133 | + # 维特比变量, 加上了从标注 i 到 next_tag 的过渡 |
| 134 | + # 的分数 |
| 135 | + # 我们在这里并没有将 emission 分数包含进来, 因为 |
| 136 | + # 最大值并不依赖于它们(我们在下面对它们进行的是相加) |
| 137 | + next_tag_var = forward_var + self.transitions[next_tag] |
| 138 | + best_tag_id = argmax(next_tag_var) |
| 139 | + bptrs_t.append(best_tag_id) |
| 140 | + viterbivars_t.append(next_tag_var[0][best_tag_id]) |
| 141 | + # 现在将所有 emission 得分相加, 将 forward_var |
| 142 | + # 赋值到我们刚刚计算出来的维特比变量集合 |
| 143 | + forward_var = (torch.cat(viterbivars_t) + feat).view(1, -1) |
| 144 | + backpointers.append(bptrs_t) |
| 145 | + |
| 146 | + # 过渡到 STOP_TAG |
| 147 | + terminal_var = forward_var + self.transitions[self.tag_to_ix[STOP_TAG]] |
| 148 | + best_tag_id = argmax(terminal_var) |
| 149 | + path_score = terminal_var[0][best_tag_id] |
| 150 | + |
| 151 | + # 跟着后指针去解码最佳路径 |
| 152 | + best_path = [best_tag_id] |
| 153 | + for bptrs_t in reversed(backpointers): |
| 154 | + best_tag_id = bptrs_t[best_tag_id] |
| 155 | + best_path.append(best_tag_id) |
| 156 | + # 弹出开始的标签 (我们并不希望把这个返回到调用函数) |
| 157 | + start = best_path.pop() |
| 158 | + assert start == self.tag_to_ix[START_TAG] # 健全性检查 |
| 159 | + best_path.reverse() |
| 160 | + return path_score, best_path |
| 161 | + |
| 162 | + def neg_log_likelihood(self, sentence, tags): |
| 163 | + feats = self._get_lstm_features(sentence) |
| 164 | + forward_score = self._forward_alg(feats) |
| 165 | + gold_score = self._score_sentence(feats, tags) |
| 166 | + return forward_score - gold_score |
| 167 | + |
| 168 | + def forward(self, sentence): # 不要把这和上面的 _forward_alg 混淆 |
| 169 | + # 得到 BiLSTM 输出分数 |
| 170 | + lstm_feats = self._get_lstm_features(sentence) |
| 171 | + |
| 172 | + # 给定特征, 找到最好的路径 |
| 173 | + score, tag_seq = self._viterbi_decode(lstm_feats) |
| 174 | + return score, tag_seq |
| 175 | + |
| 176 | + |
| 177 | +START_TAG = "<START>" |
| 178 | +STOP_TAG = "<STOP>" |
| 179 | +EMBEDDING_DIM = 5 |
| 180 | +HIDDEN_DIM = 4 |
| 181 | + |
| 182 | +# 编造一些训练数据 |
| 183 | +training_data = [( |
| 184 | + "the wall street journal reported today that apple corporation made money". |
| 185 | + split(), "B I I I O O O B I O O".split()), |
| 186 | + ("georgia tech is a university in georgia".split(), |
| 187 | + "B I O O O O B".split())] |
| 188 | + |
| 189 | +word_to_ix = {} |
| 190 | +for sentence, tags in training_data: |
| 191 | + for word in sentence: |
| 192 | + if word not in word_to_ix: |
| 193 | + word_to_ix[word] = len(word_to_ix) |
| 194 | + |
| 195 | +tag_to_ix = {"B": 0, "I": 1, "O": 2, START_TAG: 3, STOP_TAG: 4} |
| 196 | + |
| 197 | +model = BiLSTM_CRF(len(word_to_ix), tag_to_ix, EMBEDDING_DIM, HIDDEN_DIM) |
| 198 | +optimizer = optim.SGD(model.parameters(), lr=0.01, weight_decay=1e-4) |
| 199 | + |
| 200 | +# 在训练之前检查预测结果 |
| 201 | +precheck_sent = prepare_sequence(training_data[0][0], word_to_ix) |
| 202 | +precheck_tags = torch.LongTensor([tag_to_ix[t] for t in training_data[0][1]]) |
| 203 | +print(model(precheck_sent)) |
| 204 | + |
| 205 | +# 确认从之前的 LSTM 部分的 prepare_sequence 被加载了 |
| 206 | +for epoch in range(300): # 又一次, 正常情况下你不会训练300个 epoch, 这只是示例数据 |
| 207 | + for sentence, tags in training_data: |
| 208 | + # 第一步: 需要记住的是Pytorch会累积梯度 |
| 209 | + # 我们需要在每次实例之前把它们清除 |
| 210 | + model.zero_grad() |
| 211 | + |
| 212 | + # 第二步: 为我们的网络准备好输入, 即 |
| 213 | + # 把它们转变成单词索引变量 (Variables) |
| 214 | + sentence_in = prepare_sequence(sentence, word_to_ix) |
| 215 | + targets = torch.LongTensor([tag_to_ix[t] for t in tags]) |
| 216 | + |
| 217 | + # 第三步: 运行前向传递. |
| 218 | + neg_log_likelihood = model.neg_log_likelihood(sentence_in, targets) |
| 219 | + |
| 220 | + # 第四步: 计算损失, 梯度以及 |
| 221 | + # 使用 optimizer.step() 来更新参数 |
| 222 | + neg_log_likelihood.backward() |
| 223 | + optimizer.step() |
| 224 | + |
| 225 | +# 在训练之后检查预测结果 |
| 226 | +precheck_sent = prepare_sequence(training_data[0][0], word_to_ix) |
| 227 | +print(model(precheck_sent)) |
| 228 | +# 我们完成了! |
0 commit comments