Transformer详解（附代码）

引言

$Transformer\mathrm{Transformer}$ 模型是 $Google\mathrm{Google}$ 团队在 $2017$ 年 $6$ 月由 $AshishVaswani\mathrm{Ashish\text{ }Vaswani}$ 等人在论文《 $AttentionIsAllYouNeed\mathrm{Attention\text{ }Is\text{ }All \text{ }You \text{ } Need}$ 》所提出，当前它已经成为 $NLP\mathrm{NLP}$ 领域中的首选模型。 $Transformer\mathrm{Transformer}$ 抛弃了 $RNN\mathrm{RNN}$ 的顺序结构，采用了 $Self\mathrm{Self}$ - $Attention\mathrm{Attention}$ 机制，使得模型可以并行化训练，而且能够充分利用训练资料的全局信息，加入 $Transformer\mathrm{Transformer}$ 的 $Seq2seq\mathrm{Seq2seq}$ 模型在 $NLP\mathrm{NLP}$ 的各个任务上都有了显著的提升。本文做了大量的图示目的是能够更加清晰地讲解 $Transformer\mathrm{Transformer}$ 的运行原理，以及相关组件的操作细节，文末还有完整可运行的代码示例。

注意力机制

$Transformer\mathrm{Transformer}$ 中的核心机制就是 $Self\mathrm{Self}$ - $Attention\mathrm{Attention}$ 。 $Self\mathrm{Self}$ - $Attention\mathrm{Attention}$ 机制的本质来自于人类视觉注意力机制。当人视觉在感知东西时候往往会更加关注某个场景中显著性的物体，为了合理利用有限的视觉信息处理资源，人需要选择视觉区域中的特定部分，然后集中关注它。注意力机制主要目的就是对输入进行注意力权重的分配，即决定需要关注输入的哪部分，并对其分配有限的信息处理资源给重要的部分。

Self-Attention

$Self\mathrm{Self}$ - $Attention\mathrm{Attention}$ 工作原理如上图所示，给定输入 $wordembedding\mathrm{word\text{ }embedding}$ 向量 $a1,a2,a3∈Rdl×1a^1,a^2,a^3 \in \mathbb{R}^{d_l \times 1}$ ，然后对于输入向量 $ai,i∈{1,2,3}a^i,i\in \{1,2,3\}$ 通过矩阵 $Wq∈Rdk×dl,Wk∈Rdk×dl,Wv∈Rdl×dlW^q\in \mathbb{R}^{d_k \times d_l},W^k\in \mathbb{R}^{d_k \times d_l},W^v\in \mathbb{R}^{d_l\times d_l}$ 进行线性变换得到 $Query\mathrm{Query}$ 向量 $qi∈Rdk×1q^i\in\mathbb{R}^{d_k \times 1}$ ， $Key\mathrm{Key}$ 向量 $ki∈Rdk×1k^i\in \mathbb{R}^{d_k \times 1}$ ，以及 $Value\mathrm{Value}$ 向量 $vi∈Rdl×1v^i\in \mathbb{R}^{d_l \times 1}$ ，即 ${qi=Wq⋅aiki=Wk⋅ai,i∈{1,2,3}vi=Wv⋅ai\left\{\begin{aligned}q^i&=W^q \cdot a^i\\k^i&=W^k \cdot a^i,\quad i\in\{1,2,3\}\\v^i&=W^v \cdot a^i\end{aligned}\right.$

⎨

⎧qikivi=Wq⋅ai=Wk⋅ai,i∈{1,2,3}=Wv⋅ai如果令矩阵

A=(a1,a2,a3)∈Rdl×3A=(a^1,a^2,a^3)\in\mathbb{R}^{d_l \times 3}

，

Q=(q1,q2,q3)∈Rdk×3Q=(q^1,q^2,q^3)\in\mathbb{R}^{d_k \times 3}

，

K=(k1,k2,k3)∈Rdk×3K=(k^1,k^2,k^3)\in\mathbb{R}^{d_k \times 3}

，

V=(v1,v2,v3)∈Rdl×3V=(v^1,v^2,v^3)\in\mathbb{R}^{d_l \times 3}

，则此时则有

{Q=Wq⋅AK=Wk⋅AV=Wv⋅A\left\{\begin{aligned}Q&=W^q \cdot A\\K&=W^k \cdot A\\V&=W^v \cdot A\end{aligned}\right.

接着再利用得到的

Query\mathrm{Query}

向量和

Key\mathrm{Key}

向量计算注意力得分，论文中采用的注意力计算公式为点积缩放公式

αli=(qi)⊤⋅kldk=dkdk∑n=1dkknl⋅qni,i,l∈{1,2,3}\alpha^{i}_l=\frac{(q^i)^{\top}\cdot k^l}{\sqrt{d^k}}=\frac{\sqrt{d^k}}{d^k}\sum\limits_{n=1}^{d^k}k^l_n\cdot q^i_n,\quad i,l \in \{1,2,3\}

论文中假定

Key\mathrm{Key}

向量

kl=(k1l,k2l,k3l)k^l=(k^l_1,k^l_2,k^l_3)

的元素和

Query\mathrm{Query}

向量

qi=(q1i,q2i,q3i)q^i=(q^i_1,q^i_2,q^i_3)

的元素独立同分布，且令均值为

0

，方差为

1

，则此时注意力向量

ai∈R3×1a^{i}\in \mathbb{R}^{3 \times 1}

的第

l

个分量

αli\alpha^{i}_l

的均值为

0

，方差

1

具体的计算公式如下

E[αli]=dkdk∑n=1dkE[knl]⋅E[qni]=0,i,l∈{1,2,3}Var[αli]=1dk∑n=1dkVar[knl]⋅Var[qni]=1,i,l∈{1,2,3}\begin{aligned}\mathbb{E}\left[\alpha^i_l\right]&=\frac{\sqrt{d^k}}{d^k}\sum\limits_{n=1}^{d^k}\mathbb{E}\left[k^l_n\right]\cdot \mathbb{E}\left[q^i_n\right]=0,\quad i,l \in \{1,2,3\}\\ \mathrm{Var}\left[\alpha^i_l\right]&=\frac{1}{d^k}\sum\limits_{n=1}^{d^k}\mathrm{Var}\left[k^l_n\right]\cdot \mathrm{Var}\left[q^i_n\right]=1,\quad i,l \in \{1,2,3\}\end{aligned}

令注意力分数矩阵

Λ=(α1,α2,α3)∈R3×3\Lambda=(\alpha^1,\alpha^2,\alpha^3)\in \mathbb{R}^{3 \times 3}

，则有

Λ=K⊤⋅Qdk\Lambda=\frac{K^{\top}\cdot Q}{\sqrt{d^k}}

注意分数向量

αi\alpha^i

经过

softmax\mathrm{softmax}

层得到归一化后的注意力分布

βi\beta^i

，即为

βji=eαji∑n=13eαni,i,j={1,2,3}\beta^i_j = \frac{e^{\alpha^{i}_j}}{\sum\limits_{n=1}^3e^{\alpha^{i}_n}},\quad i,j=\{1,2,3\}

最后利用得到的注意力分布向量

βi\beta^i

和

Value\mathrm{Value}

矩阵

V

获得最后的输出

bi∈Rdl×1b^i\in \mathbb{R}^{d_l \times 1}

，则有

bi=∑l=13βli⋅vl,i∈{1,2,3}b^i=\sum\limits^{3}_{l=1}\beta^{i}_l \cdot v^{l},\quad i \in \{1,2,3\}

令输出矩阵

B=(b1,b2,b3)∈Rdl×3B=(b^1,b^2,b^3)\in\mathbb{R}^{d_l\times 3}

，则有

B=Attention(Q,K,V)=V⋅softmax(K⊤⋅Qdk)B=\mathrm{Attention}(Q,K,V)=V\cdot\mathrm{softmax}\left(\frac{K^{\top}\cdot Q}{\sqrt{d^k}}\right)

Multi-Head Attention

$Multi\mathrm{Multi}$ - $HeadAttention\mathrm{Head\text{ }Attention}$ 的工作原理与 $Self\mathrm{Self}$ - $Attention\mathrm{Attention}$ 的工作原理非常类似。为了方便图解可视化将 $Multi\mathrm{Multi}$ - $Head\mathrm{Head}$ 设置为 $2$ - $Head\mathrm{Head}$ ，如果 $Multi\mathrm{Multi}$ - $Head\mathrm{Head}$ 设置为 $8$ - $Head\mathrm{Head}$ ，则上图的 $qi,ki,vi,i∈{1,2,3}q^i,k^i,v^i,i\in\{1,2,3\}$ 的下一步的分支数为 $8$ 。给定输入 $wordembedding\mathrm{word\text{ }embedding}$ 向量 $a1,a2,a3∈Rdl×1a^1,a^2,a^3 \in \mathbb{R}^{d_l \times 1}$ ，然后对于输入向量 $ai,i∈{1,2,3}a^i,i\in \{1,2,3\}$ 通过矩阵 $Wq∈Rdk×dl,Wk∈Rdk×dl,Wv∈Rdl×dlW^q\in \mathbb{R}^{d_k \times d_l},W^k\in \mathbb{R}^{d_k \times d_l},W^v\in \mathbb{R}^{d_l\times d_l}$ 进行第一次线性变换得到 $Query\mathrm{Query}$ 向量 $qi∈Rdk×1q^i\in\mathbb{R}^{d_k \times 1}$ ， $Key\mathrm{Key}$ 向量 $ki∈Rdk×1k^i \in\mathbb{R}^{d_k \times 1}$ ，以及 $Value\mathrm{Value}$ 向量 $vi∈Rdl×1v^i \in\mathbb{R}^{d_l \times 1}$ 。然后再对 $Query\mathrm{Query}$ 向量 $q^i$ 通过矩阵 $Wq1∈Rdm×dkW^{q1}\in \mathbb{R}^{d_m \times d_k}$ 和 $Wq2∈Rdm×dkW^{q2}\in \mathbb{R}^{d_m\times d_k}$ 进行第二次线性变换得到 $qi1∈Rdm×1q^{i1}\in \mathbb{R}^{d_m \times 1}$ 和 $qi2∈Rdm×1q^{i2}\in \mathbb{R}^{d_m\times 1}$ ，同理对 $Key\mathrm{Key}$ 向量 $k^i$ 通过矩阵 $Wk1∈Rdm×dkW^{k1}\in \mathbb{R}^{d_m \times d_k}$ 和 $Wk2∈Rdm×dkW^{k2}\in \mathbb{R}^{d_m\times d_k}$ 进行第二次线性变换得到 $ki1∈Rdm×1k^{i1}\in \mathbb{R}^{d_m\times 1}$ 和 $ki2∈Rdm×1k^{i2}\in \mathbb{R}^{d_m\times 1}$ ，对 $Value\mathrm{Value}$ 向量 $v^i$ 通过矩阵 $Wv1∈Rdl2×dlW^{v1}\in \mathbb{R}^{\frac{d_l}{2}\times d_l}$ 和 $Wv2∈Rdl2×dlW^{v2}\in \mathbb{R}^{\frac{d_l}{2}\times d_l}$ 进行第二次线性变换得到 $vi1∈Rdl2×1v^{i1}\in \mathbb{R}^{\frac{d_l}{2}\times 1}$ 和 $vi2∈Rdl2×1v^{i2}\in \mathbb{R}^{\frac{d_l}{2}\times 1}$ ,具体的计算公式如下所示： ${qih=Wqh⋅Wq⋅aikih=Wkh⋅Wk⋅ai,i={1,2,3},h={1,2}vih=Wvh⋅Wv⋅ai\left\{\begin{aligned}q^{ih}&=W^{qh}\cdot W^{q} \cdot a^i\\ k^{ih}&=W^{kh}\cdot W^{k} \cdot a^i,\quad i=\{1,2,3\},\quad h=\{1,2\}\\v^{ih}&=W^{vh}\cdot W^{v} \cdot a^i\end{aligned}\right.$

⎨

⎧qihkihvih=Wqh⋅Wq⋅ai=Wkh⋅Wk⋅ai,i={1,2,3},h={1,2}=Wvh⋅Wv⋅ai令矩阵

Q1=(q11,q21,q31)∈Rdm×3Q2=(q12,q22,q32)∈Rdm×3K1=(k11,k21,k31)∈Rdm×3K2=(k12,k22,k32)∈Rdm×3V1=(v11,v21,v31)∈Rdl2×3V2=(v12,v22,v32)∈Rdl2×3\begin{array}{ll}Q^{1}=(q^{11},q^{21},q^{31})\in \mathbb{R}^{d_m\times 3}&\quad Q^2=(q^{12},q^{22},q^{32})\in\mathbb{R}^{d_m\times 3}\\K^{1}=(k^{11},k^{21},k^{31})\in \mathbb{R}^{d_m\times 3}&\quad K^2=(k^{12},k^{22},k^{32})\in\mathbb{R}^{d_m\times 3}\\V^{1}=(v^{11},v^{21},v^{31})\in \mathbb{R}^{\frac{d_l}{2}\times 3}&\quad V^2=(v^{12},v^{22},v^{32})\in\mathbb{R}^{\frac{d_l}{2}\times 3}\end{array}

此时则有

Q1=Wq1⋅Wq⋅AQ2=Wq2⋅Wq⋅AK1=Wk1⋅Wk⋅AK2=Wk2⋅Wk⋅AV1=Wv1⋅Wv⋅AV2=Wv2⋅Wv⋅A\begin{array}{ll}Q^{1}=W^{q1}\cdot W^{q} \cdot A &\quad Q^2=W^{q2}\cdot W^{q} \cdot A\\K^{1}=W^{k1}\cdot W^{k} \cdot A&\quad K^2=W^{k2}\cdot W^{k} \cdot A\\V^{1}=W^{v1}\cdot W^{v} \cdot A&\quad V^2=W^{v2}\cdot W^{v} \cdot A\end{array}

对于每个

Head\mathrm{Head}

利用得到对于

Query\mathrm{Query}

向量和

Key\mathrm{Key}

向量计算对应的注意力得分，其中注意力向量

αih\alpha^{ih}

的第

l

个分量的计算公式为

αlih=(qih)⊤⋅klh,i∈{1,2,3},h∈{1,2},l∈{1,2,3}\alpha^{ih}_l=(q^{ih})^{\top}\cdot k^{lh},\quad i\in\{1,2,3\},h\in\{1,2\},l\in\{1,2,3\}

令注意力分数矩阵

Λ1=(α11,α21,α31)\Lambda^1=(\alpha^{11},\alpha^{21},\alpha^{31})

，

Λ2=(α12,α22,α32)\Lambda^2=(\alpha^{12},\alpha^{22},\alpha^{32})

，则有

Λ1=(K1)⊤⋅Q1dm,Λ2=(K2)⊤⋅Q2dm\Lambda^{1}=\frac{(K^1)^{\top}\cdot Q^1}{\sqrt{d_m}},\quad\Lambda^{2}=\frac{(K^2)^{\top}\cdot Q^2}{\sqrt{d_m}}

注意分数向量

αih\alpha^{ih}

经过

softmax\mathrm{softmax}

层得到归一化后的注意力分布

βih\beta^{ih}

，即为

βjih=eαjih∑n=13eαnih,i,j={1,2,3},h={1,2}\beta^{ih}_j = \frac{e^{\alpha^{ih}_j}}{\sum\limits_{n=1}^3e^{\alpha^{ih}_n}},\quad i,j=\{1,2,3\}, h=\{1,2\}

对于每一个

Head\mathrm{Head}

利用得到的注意力分布向量

βih\beta^{ih}

和

Value\mathrm{Value}

矩阵

V^h

获得最后的输出

bih∈Rdl2×1b^{ih}\in \mathbb{R}^{\frac{d_l}{2} \times 1}

，则有

bih=∑l=13βlih⋅vlh,i∈{1,2,3},h∈{1,2}b^{ih}=\sum\limits^{3}_{l=1}\beta^{ih}_l \cdot v^{lh},\quad i \in \{1,2,3\}, h\in\{1,2\}

两个

Head\mathrm{Head}

的

b^{ih}

的向量按照如下方式拼接在一起，则有

B=(b11b21b31b12b22b32)∈Rdl×3B=\left(\begin{array}{lll}b^{11}&b^{21}&b^{31}\\b^{12}&b^{22}&b^{32}\end{array}\right)\in \mathbb{R}^{d_l \times 3}

给定参数矩阵

WO∈Rdl×dlW^{O}\in \mathbb{R}^{d_l\times d_l}

，则输出矩阵为

O=WO⋅B∈Rdl×3O=W^{O}\cdot B\in \mathbb{R}^{d_l \times 3}

综上所述则有

O=MultiHead(Q,K,V)=WO⋅Concat(V1⋅softmax((K1)⊤⋅Q1dm)V2⋅softmax((K2)⊤⋅Q2dm))O=\mathrm{MultiHead}(Q,K,V)=W^O\cdot\mathrm{Concat}\left(\begin{array}{l}V^1\cdot\mathrm{softmax}\left(\frac{(K^1)^{\top}\cdot Q^1}{\sqrt{d_m}}\right)\\ \\V^2\cdot\mathrm{softmax}\left(\frac{(K^2)^{\top}\cdot Q^2}{\sqrt{d_m}}\right)\end{array}\right)

Mask Self-Attention

如下图左半部分所示， $Self\mathrm{Self}$ - $Attention\mathrm{Attention}$ 的输出向量 $bi,i∈{1,2,3,4}b^i, i \in \{1,2,3,4\}$ 综合了输入向量 $ai,i∈{1,2,3,4}a^i, i \in \{1,2,3,4\}$ 的全部信息，由此可见， $Self\mathrm{Self}$ - $Attention\mathrm{Attention}$ 在实际编程中支持并行运算。如下图右半部分所示， $MaskSelf\mathrm{Mask \text{ } Self}$ - $Attention\mathrm{Attention}$ 的输出向量 $b^i$ 只利用了已知部分输入的向量 $a^i$ 的信息。例如， $b 1$ 只是与 $a^1$ 有关； $b^2$ 与 $a^1$ 和 $a^2$ 有关； $b^3$ 与 $a^1$ ， $a^2$ 和 $a^3$ 有关； $b^4$ 与 $a^1$ ， $a^2$ ， $a^3$ 和 $a^4$ 有关。 $MaskSelf\mathrm{Mask \text{ } Self}$ - $Attention\mathrm{Attention}$ 在 $Transformer\mathrm{Transformer}$ 中被用到过两次。

$Transformer\mathrm{Transformer}$ 的 $Encoder\mathrm{Encoder}$ 中如果输入一句话的 $word\mathrm{word}$ 长度小于指定的长度，为了能够让长度一致往往会用 $0$ 进行填充，此时则需要用 $MaskSelf\mathrm{Mask \text{ } Self}$ - $Attention\mathrm{Attention}$ 来计算注意力分布。
$Transformer\mathrm{Transformer}$ 的 $Decoder\mathrm{Decoder}$ 的输出是有时序关系的，当前的输出只与之前的输入有关，所以此时算注意力分布时需要用到 $MaskSelf\mathrm{Mask \text{ } Self}$ - $Attention\mathrm{Attention}$ 。

Transformer模型

以上对 $Transformer\mathrm{Transformer}$ 中的核心内容即自注意力机制进行了详细解剖，接下来会对 $Transformer\mathrm{Transformer}$ 模型架构进行介绍。 $Transformer\mathrm{Transformer}$ 模型是由 $Encoder\mathrm{Encoder}$ 和 $Decoder\mathrm{Decoder}$ 两个模块组成，具体的示意图如下所示，为了能够对 $Transformer\mathrm{Transformer}$ 内部的操作细节进行更清晰的展示，下图以矩阵运算的视角对 $Transformer\mathrm{Transformer}$ 的原理进行讲解。
$Encoder\mathrm{Encoder}$ 模块操作的具体流程如下所示：

$Encoder\mathrm{Encoder}$ 的输入由两部分组成分别是词编码矩阵 $\in \mathbb{R}^{n \times l \times d}$ 和位置编码矩阵 $\in \mathbb{R}^{n \times l \times d}$ ，其中 $n$ 表示句子数目， $l$ 表示一句话单词的最大数目， $d$ 表示的是词向量的维度。位置编码矩阵 $P$ 表示的是每个单词在一句里的所有位置信息，因为 $Self\mathrm{Self}$ - $Attention\mathrm{Attention}$ 计算注意力分布的时候只能给出输出向量和输入向量之间的权重关系，但是不能给出词在一句话里的位置信息，所以需要在输入里引入位置编码矩阵 $P$ 。位置编码向量生成方法有很多。一种比较简单粗暴的方式就是根据单词在句子中的位置生成一个 $one\mathrm{one}$ - $hot\mathrm{hot}$ 的位置编码；还有的方法是将位置编码当成参数进行训练学习；在该论文里是利用三角函数对位置进行编码，具体的公式如下所示 $PE(pos,2i)=sin⁡(pos10002i/d),PE(pos,2i+1)=cos⁡(pos10002i/d)\mathrm{PE}(pos,2i)=\sin(\frac{pos}{1000^{2i/d}}),\quad \mathrm{PE}(pos,2i+1)=\cos(\frac{pos}{1000^{2i/d}})$ 其中 $PE\mathrm{PE}$ 表示的是位置编码向量， $p os$ 表示词在句子中的位置， $i$ 表示编码向量的位置索引。
输入矩阵 $I + P$ 通过线性变换生成矩阵 $Q$ ， $K$ ， $V$ 。在实际编程中是将输入 $I + P$ 直接赋值给 $Q$ ， $K$ ， $V$ 。如果输入单词长度小于最大长度并 $0$ 来填充的时候，还要相应引入 $Mask\mathrm{Mask}$ 矩阵。
将矩阵 $Q$ ， $K$ ， $V$ 输入到 $Multi\mathrm{Multi}$ - $HeadAttention\mathrm{Head\text{ }Attention}$ 模块中进行注意分布的计算得到矩阵 $I′∈Rn×l×dI^{\prime}\in \mathbb{R}^{n \times l \times d}$ ，计算公式为 $I′=MultiHead(Q,K,V)I^{\prime}=\mathrm{MultiHead}(Q,K,V)$ 具体的计算细节参考上文关于 $Multi\mathrm{Multi}$ - $HeadAttention\mathrm{Head\text{ }Attention}$ 原理的讲解不在这里赘述。然后将原始输入 $I + P$ 与注意力分布 $I′I^{\prime}$ 进行残差计算得到输出矩阵 $I+P+I′∈Rn×l×dI+P+I^{\prime}\in \mathbb{R}^{n \times l \times d}$ 。
对矩阵 $I+P+I′={xijk}nldI+P+I^{\prime}=\{x_{ijk}\}^{nld}$ 进行层归一化操作得到 $I′′∈Rn×l×dI^{\prime\prime}\in\mathbb{R}^{n \times l \times d}$ ，具体的计算公式为 ${μij=∑k=1dxijkσij=∑k=1d(xijk−μij)2⟹x^ijk=xijk−uijσij,i∈{1,⋯,n},j∈{1,⋯,l},k∈{1,⋯,d}\left\{\begin{aligned}\mu^{ij}&=\sum\limits_{k=1}^d x_{ijk}\\\sigma^{ij}&=\sqrt{\sum\limits_{k=1}^d\left(x_{ijk}-\mu^{ij}\right)^2}\end{aligned}\right. \Longrightarrow \hat{x}_{ijk}=\frac{x_{ijk}-u^{ij}}{\sigma^{ij}},\quad i\in\{1,\cdots,n\},j\in\{1,\cdots,l\},k\in\{1,\cdots,d\}$
将 $I′′I^{\prime\prime}$ 输入到全连接神经网络中得到 $I′′′∈Rn×l×dI^{\prime\prime\prime}\in \mathbb{R}^{n \times l \times d}$ ，然后再让全连接神经网络的输入 $I′′I^{\prime\prime}$ 与输出 $I′′′I^{\prime\prime\prime}$ 进行残差计算得到 $I′′+I′′′I^{\prime\prime}+I^{\prime\prime\prime}$ ，接着对 $I′′+I′′′I^{\prime\prime}+I^{\prime\prime\prime}$ 进行层归一化操作。
以上是一个 $Block\mathrm{Block}$ 的操作原理，将 $N$ 个 $Block\mathrm{Block}$ 进行堆叠就组成了 $Encoder\mathrm{Encoder}$ 的模块，得到的最后输出为 $IN∈Rn×l×dI^N \in \mathbb{R}^{n \times l \times d}$ 。这里需要注意的是 $Encoder\mathrm{Encoder}$ 模块中的各个组件的操作顺序并不是固定的，也可以先进行归一化操作，然后再计算注意力分布，再归一化，再预测等。

$Decoder\mathrm{Decoder}$ 模块操作的具体流程如下所示：

$Decoder\mathrm{Decoder}$ 的输入也由两部分组成分别是词编码矩阵 $\in \mathbb{R}^{n_1 \times l_1 \times d}$ 和位置编码矩阵 $PO∈Rn1×l1×dP^O \in \mathbb{R}^{n_1 \times l_1 \times d}$ 。因为 $Decoder\mathrm{Decoder}$ 的输入是具有时顺序关系的（即上一步的输出为当前步输入）所以还需要输入 $Mask\mathrm{Mask}$ 矩阵 $M$ 以便计算注意力分布。
输入矩阵 $O+P^O$ 通过线性变换生成矩阵 $Q^\hat{Q}$ ， $K^\hat{K}$ ， $V^\hat{V}$ 。在实际编程中是将输入 $O+P^O$ 直接赋值给 $Q^\hat{Q}$ ， $K^\hat{K}$ ， $V^\hat{V}$ 。如果输入单词长度小于最大长度并 $0$ 来填充的时候，还要相应引入 $Mask\mathrm{Mask}$ 矩阵。
将矩阵 $Q^\hat{Q}$ ， $K^\hat{K}$ ， $V^\hat{V}$ 以及 $Mask\mathrm{Mask}$ 矩阵 $M$ 输入到 $MaskMulti\mathrm{Mask\text{ }Multi}$ - $HeadAttention\mathrm{Head\text{ }Attention}$ 模块中进行注意分布的计算得到矩阵 $O′∈Rn1×l1×dO^{\prime}\in \mathbb{R}^{n_1 \times l_1 \times d}$ ，计算公式为 $O′=MaskMultiHead(Q^,K^,V^,M)O^{\prime}=\mathrm{MaskMultiHead}(\hat{Q},\hat{K},\hat{V},M)$ 具体的计算细节参考上文关于 $MaskSelf\mathrm{Mask \text{ }Self}$ - $Attention\mathrm{Attention}$ 的讲解不在这里赘述。然后将原始输入 $O+P^O$ 与注意力分布 $O′O^{\prime}$ 进行残差计算得到输出矩阵 $O+PO+O′∈Rn1×l1×dO+P^O+O^{\prime}\in \mathbb{R}^{n_1 \times l_1 \times d}$ 。接着再对矩阵 $O+PO+O′O+P^O+O^{\prime}$ 进行层归一化操作得到 $O′′∈Rn1×l1×dO^{\prime\prime}\in\mathbb{R}^{n_1 \times l_1 \times d}$ 。
$Encoder\mathrm{Encoder}$ 的输出 $I^N$ 通过线性变换得到 $Q^N$ 和 $K^N$ ， $O′O^{\prime}$ 进行线性变换得到 $V^′\hat{V}^{\prime}$ ，利用矩阵 $Q^N$ 和 $K^N$ 和 $V^′\hat{V}^{\prime}$ 进行交叉注意力分布的计算得到 $O′′′O^{\prime\prime\prime}$ ，计算公式为 $O′′′=MultiHead(QN,KN,V^′)O^{\prime\prime\prime}=\mathrm{MultiHead}(Q^N,K^N,\hat{V}^{\prime})$ 这里的交叉注意力分布综合 $Encoder\mathrm{Encoder}$ 输出结果和 $Decoder\mathrm{Decoder}$ 中间结果的信息。实际编程编程中将 $I^N$ 直接赋值给 $Q^\hat{Q}$ 和 $K^\hat{K}$ ， $O′O^{\prime}$ 直接赋值给 $V^′\hat{V}^{\prime}$ 。然后将 $O′′O^{\prime\prime}$ 与注意力分布 $O′′′O^{\prime\prime\prime}$ 进行残差计算得到输出矩阵 $O′′+O′′′O^{\prime\prime}+O^{\prime\prime\prime}$ 。
接着对 $O′′+O′′′O^{\prime\prime}+O^{\prime\prime\prime}$ 进行层归一操作得到 $O′′′′O^{\prime\prime\prime\prime}$ ，再将 $O′′′′O^{\prime\prime\prime\prime}$ 输入到全连接神经网络中得到 $O′′′′′O^{\prime\prime\prime\prime\prime}$ ，接着再做一步残差操作得到 $O′′′′+O′′′′′O^{\prime\prime\prime\prime}+O^{\prime\prime\prime\prime\prime}$ ，最后再进行一层归一化操作。
以上是一个 $Block\mathrm{Block}$ 的操作原理，将 $N$ 个 $Block\mathrm{Block}$ 进行堆叠就组成了 $Decoder\mathrm{Decoder}$ 的模块，得到的输出为 $ON∈Rn1×l1×dO^N \in \mathbb{R}^{n_1 \times l_1 \times d}$ 。然后在词汇字典中找到当前预测最大概率的单词，并将该单词词向量作为下一阶段的输入，重复以上步骤，直到输出“ $end\mathrm{end}$ ”字符为止。

代码示例

$Transformer\mathrm{Transformer}$ 具体的代码示例如下所示为一个国外博主视频里的代码，并根据上文对代码的一些细节进行了探讨。根据上文中 $Multi\mathrm{Multi}$ - $HeadAttention\mathrm{Head\text{ }Attention}$ 原理示例图可知，严格来看 $Multi\mathrm{Multi}$ - $HeadAttention\mathrm{Head\text{ }Attention}$ 在求注意分布的时候中间其实是有两步线性变换。给定输入向量 $x∈R256×1x\in \mathbb{R}^{256\times 1}$ 第一步线性变换直接让向量 $x$ 赋值给 $q$ ， $k$ ， $v$ ，这一过程以下程序中有所体现，在这里并不会产生歧义。第二步线性变换产生多 $Head\mathrm{Head}$ ，假设 $Head=8\mathrm{Head}=8$ 的时候，按理说 $q$ 要与 $8$ 个矩阵 $Wq1,⋯,Wq8W^{q1},\cdots,W^{q8}$ 进行线性变换得到 $8$ 个 $q1,⋯,q8q^{1},\cdots,q^{8}$ ，同理 $k$ 要与 $8$ 个矩阵 $Wk1,⋯,Wk8W^{k1},\cdots,W^{k8}$ 进行线性变换得到 $8$ 个 $k1,⋯,k8k^{1},\cdots,k^{8}$ ， $v$ 要与 $8$ 个矩阵 $Wv1,⋯,Wv8W^{v1},\cdots,W^{v8}$ 进行线性变换得到 $8$ 个 $v1,⋯,v8v^{1},\cdots,v^{8}$ ，如果按照这个方式在程序实现则需要定义24个权重矩阵，非常的麻烦。以下程序中有一个简单的权重定义方法，通过该方法也可以实现以上多 $Head\mathrm{Head}$ 的线性变换，以向量 $(q_1,\cdots, q_{256})^{\top}\in \mathbb{R}^{256 \times 1}$ 为例：

首先将向量 $q$ 进行截断分成 $Head=8\mathrm{Head}=8$ 个向量，即为 ${q(1)=(E,0,0,0,0,0,0,0)⋅qq(2)=(0,E,0,0,0,0,0,0)⋅qq(3)=(0,0,E,0,0,0,0,0)⋅qq(4)=(0,0,0,E,0,0,0,0)⋅qq(5)=(0,0,0,0,E,0,0,0)⋅qq(6)=(0,0,0,0,0,E,0,0)⋅qq(7)=(0,0,0,0,0,0,E,0)⋅qq(8)=(0,0,0,0,0,0,0,E)⋅q\left\{\begin{aligned}q^{(1)}&=({\bf{E},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})\cdot q\\q^{(2)}&=({\bf{0},\bf{E},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})\cdot q\\q^{(3)}&=({\bf{0},\bf{0},\bf{E},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})\cdot q\\q^{(4)}&=({\bf{0},\bf{0},\bf{0},\bf{E},\bf{0},\bf{0},\bf{0},\bf{0}})\cdot q\\q^{(5)}&=({\bf{0},\bf{0},\bf{0},\bf{0},\bf{E},\bf{0},\bf{0},\bf{0}})\cdot q\\q^{(6)}&=({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{E},\bf{0},\bf{0}})\cdot q\\q^{(7)}&=({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{E},\bf{0}})\cdot q\\q^{(8)}&=({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{E}})\cdot q \end{aligned}\right.$ 其中 $q(i)∈R32×1q^{(i)}\in \mathbb{R}^{32\times 1}$ 是 $q$ 的第 $i$ 个截断向量， $E∈R32×32{\bf{E}}\in \mathbb{R}^{32 \times 32}$ 是单位矩阵， $0∈R32×32{\bf{0}}\in \mathbb{R}^{32 \times 32}$ 是零矩阵。
然后对 $q(i),i∈{1,⋯,8}q^{(i)},i\in \{1,\cdots,8\}$ 用相同的权重矩阵 $\in \mathbb{R}^{32 \times 32}$ 进行线性变换，此时可以发现，训练过程的时候只需要更新权重矩阵 $W$ 即可，而且可以进行多 $Head\mathrm{Head}$ 线性变换， $8$ 个权重矩阵可以表示为： ${Wq1=W⋅(E,0,0,0,0,0,0,0)=(W,0,0,0,0,0,0,0)Wq2=W⋅(0,E,0,0,0,0,0,0)=(0,W,0,0,0,0,0,0)Wq3=W⋅(0,0,E,0,0,0,0,0)=(0,0,W,0,0,0,0,0)Wq4=W⋅(0,0,0,E,0,0,0,0)=(0,0,0,W,0,0,0,0)Wq5=W⋅(0,0,0,0,E,0,0,0)=(0,0,0,0,W,0,0,0)Wq6=W⋅(0,0,0,0,0,E,0,0)=(0,0,0,0,0,W,0,0)Wq7=W⋅(0,0,0,0,0,0,E,0)=(0,0,0,0,0,0,W,0)Wq8=W⋅(0,0,0,0,0,0,0,E)=(0,0,0,0,0,0,0,W)\left\{\begin{aligned}W^{q1}&=W\cdot ({\bf{E},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})=(W,{\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})\\W^{q2}&=W\cdot ({\bf{0},\bf{E},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})=({\bf{0},}W{,\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})\\W^{q3}&=W\cdot ({\bf{0},\bf{0},\bf{E},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})=({\bf{0},\bf{0},}W{,\bf{0},\bf{0},\bf{0},\bf{0},\bf{0}})\\W^{q4}&=W\cdot ({\bf{0},\bf{0},\bf{0},\bf{E},\bf{0},\bf{0},\bf{0},\bf{0}})=({\bf{0},\bf{0},\bf{0},}W{,\bf{0},\bf{0},\bf{0},\bf{0}})\\W^{q5}&=W\cdot ({\bf{0},\bf{0},\bf{0},\bf{0},\bf{E},\bf{0},\bf{0},\bf{0}})=({\bf{0},\bf{0},\bf{0},\bf{0},}W{,\bf{0},\bf{0},\bf{0}})\\W^{q6}&=W\cdot ({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{E},\bf{0},\bf{0}})=({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},}W{,\bf{0},\bf{0}})\\W^{q7}&=W\cdot ({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{E},\bf{0}})=({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},}W{,\bf{0}})\\W^{q8}&=W\cdot ({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{E}})=({\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},\bf{0},}W{})\end{aligned}\right.$ 其中权重矩阵 $Wqi∈R32×256,i∈{1,⋯,8}W^{qi}\in\mathbb{R}^{32 \times 256},i\in\{1,\cdots,8\}$ 。

import torch
import torch.nn as nn
import osclass SelfAttention(nn.Module):def __init__(self, embed_size, heads):super(SelfAttention, self).__init__()self.embed_size = embed_sizeself.heads = headsself.head_dim = embed_size // headsassert (self.head_dim * heads == embed_size), "Embed size needs to be div by heads"self.values = nn.Linear(self.head_dim, self.head_dim, bias=False)self.keys = nn.Linear(self.head_dim, self.head_dim, bias=False)self.queries = nn.Linear(self.head_dim, self.head_dim, bias=False)self.fc_out = nn.Linear(heads * self.head_dim, embed_size)def forward(self, values, keys, query, mask):N =query.shape[0]value_len , key_len , query_len = values.shape[1], keys.shape[1], query.shape[1]# split embedding into self.heads piecesvalues = values.reshape(N, value_len, self.heads, self.head_dim)keys = keys.reshape(N, key_len, self.heads, self.head_dim)queries = query.reshape(N, query_len, self.heads, self.head_dim)values = self.values(values)keys = self.keys(keys)queries = self.queries(queries)energy = torch.einsum("nqhd,nkhd->nhqk", queries, keys)# queries shape: (N, query_len, heads, heads_dim)# keys shape : (N, key_len, heads, heads_dim)# energy shape: (N, heads, query_len, key_len)if mask is not None:energy = energy.masked_fill(mask == 0, float("-1e20"))attention = torch.softmax(energy/ (self.embed_size ** (1/2)), dim=3)out = torch.einsum("nhql, nlhd->nqhd", [attention, values]).reshape(N, query_len, self.heads*self.head_dim)# attention shape: (N, heads, query_len, key_len)# values shape: (N, value_len, heads, heads_dim)# (N, query_len, heads, head_dim)out = self.fc_out(out)return outclass TransformerBlock(nn.Module):def __init__(self, embed_size, heads, dropout, forward_expansion):super(TransformerBlock, self).__init__()self.attention = SelfAttention(embed_size, heads)self.norm1 = nn.LayerNorm(embed_size)self.norm2 = nn.LayerNorm(embed_size)self.feed_forward = nn.Sequential(nn.Linear(embed_size, forward_expansion*embed_size),nn.ReLU(),nn.Linear(forward_expansion*embed_size, embed_size))self.dropout = nn.Dropout(dropout)def forward(self, value, key, query, mask):attention = self.attention(value, key, query, mask)x = self.dropout(self.norm1(attention + query))forward = self.feed_forward(x)out = self.dropout(self.norm2(forward + x))return outclass Encoder(nn.Module):def __init__(self,src_vocab_size,embed_size,num_layers,heads,device,forward_expansion,dropout,max_length,):super(Encoder, self).__init__()self.embed_size = embed_sizeself.device = deviceself.word_embedding = nn.Embedding(src_vocab_size, embed_size)self.position_embedding = nn.Embedding(max_length, embed_size)self.layers = nn.ModuleList([TransformerBlock(embed_size,heads,dropout=dropout,forward_expansion=forward_expansion,)for _ in range(num_layers)])self.dropout = nn.Dropout(dropout)def forward(self, x, mask):N, seq_length = x.shapepositions = torch.arange(0, seq_length).expand(N, seq_length).to(self.device)out = self.dropout(self.word_embedding(x) + self.position_embedding(positions))for layer in self.layers:out = layer(out, out, out, mask)return outclass DecoderBlock(nn.Module):def __init__(self, embed_size, heads, forward_expansion, dropout, device):super(DecoderBlock, self).__init__()self.attention = SelfAttention(embed_size, heads)self.norm = nn.LayerNorm(embed_size)self.transformer_block = TransformerBlock(embed_size, heads, dropout, forward_expansion)self.dropout = nn.Dropout(dropout)def forward(self, x, value, key, src_mask, trg_mask):attention = self.attention(x, x, x, trg_mask)query = self.dropout(self.norm(attention + x))out = self.transformer_block(value, key, query, src_mask)return outclass Decoder(nn.Module):def __init__(self,trg_vocab_size,embed_size,num_layers,heads,forward_expansion,dropout,device,max_length,):super(Decoder, self).__init__()self.device = deviceself.word_embedding = nn.Embedding(trg_vocab_size, embed_size)self.position_embedding = nn.Embedding(max_length, embed_size)self.layers = nn.ModuleList([DecoderBlock(embed_size, heads, forward_expansion, dropout, device)for _ in range(num_layers)])self.fc_out = nn.Linear(embed_size, trg_vocab_size)self.dropout = nn.Dropout(dropout)def forward(self, x ,enc_out , src_mask, trg_mask):N, seq_length = x.shapepositions = torch.arange(0, seq_length).expand(N, seq_length).to(self.device)x = self.dropout((self.word_embedding(x) + self.position_embedding(positions)))for layer in self.layers:x = layer(x, enc_out, enc_out, src_mask, trg_mask)out =self.fc_out(x)return outclass Transformer(nn.Module):def __init__(self,src_vocab_size,trg_vocab_size,src_pad_idx,trg_pad_idx,embed_size = 256,num_layers = 6,forward_expansion = 4,heads = 8,dropout = 0,device="cuda",max_length=100):super(Transformer, self).__init__()self.encoder = Encoder(src_vocab_size,embed_size,num_layers,heads,device,forward_expansion,dropout,max_length)self.decoder = Decoder(trg_vocab_size,embed_size,num_layers,heads,forward_expansion,dropout,device,max_length)self.src_pad_idx = src_pad_idxself.trg_pad_idx = trg_pad_idxself.device = devicedef make_src_mask(self, src):src_mask = (src != self.src_pad_idx).unsqueeze(1).unsqueeze(2)# (N, 1, 1, src_len)return src_mask.to(self.device)def make_trg_mask(self, trg):N, trg_len = trg.shapetrg_mask = torch.tril(torch.ones((trg_len, trg_len))).expand(N, 1, trg_len, trg_len)return trg_mask.to(self.device)def forward(self, src, trg):src_mask = self.make_src_mask(src)trg_mask = self.make_trg_mask(trg)enc_src = self.encoder(src, src_mask)out = self.decoder(trg, enc_src, src_mask, trg_mask)return outif __name__ == '__main__':device = torch.device("cuda" if torch.cuda.is_available() else "cpu")print(device)x = torch.tensor([[1,5,6,4,3,9,5,2,0],[1,8,7,3,4,5,6,7,2]]).to(device)trg = torch.tensor([[1,7,4,3,5,9,2,0],[1,5,6,2,4,7,6,2]]).to(device)src_pad_idx = 0trg_pad_idx = 0src_vocab_size = 10trg_vocab_size = 10model = Transformer(src_vocab_size, trg_vocab_size, src_pad_idx, trg_pad_idx, device=device).to(device)out = model(x, trg[:, : -1])print(out.shape)