旋转位置编码（ROPE）详解：从Transformer到现代前沿

旋转位置编码（ROPE）详解：从Transformer到现代前沿

标签：NLP, Transformer, 位置编码, ROPE, 深度学习, 机器学习

摘要：本文详细介绍了旋转位置编码（ROPE）在Transformer模型中的应用，包括其数学原理、代码实现以及与传统位置编码的对比。ROPE通过旋转嵌入向量来编码位置信息，能够自然捕捉相对位置关系，特别适用于长序列任务。文中提供了ROPE的PyTorch实现和使用示例，旨在帮助读者深入理解并应用这一先进技术。

引言

在自然语言处理（NLP）领域，Transformer模型因其强大的并行计算能力和对长距离依赖的捕捉而成为主流。然而，Transformer中的自注意力机制（Self-Attention）天生无序，无法感知输入序列中词语的相对位置。为了解决这一问题，研究者们提出了多种位置编码（Positional Encoding）方案，其中旋转位置编码（Rotary Position Embedding, ROPE）凭借其独特的旋转机制和在长序列任务中的优异表现，成为近年来备受关注的技术。

本文将从基础概念出发，逐步深入ROPE的数学原理、实现细节及其在Transformer模型中的应用，旨在为读者提供一篇最详细的中文博客。无论您是NLP初学者还是资深研究者，本文都将帮助您全面理解ROPE的精髓。我们将涵盖从简单概念到复杂结论的介绍，包括数学推导和代码实现。

---

1. 为什么需要位置编码？

在介绍ROPE之前，我们首先需要理解为什么Transformer模型需要位置编码。

1.1 Transformer的自注意力机制

Transformer模型的核心是自注意力机制，它通过计算输入序列中每个词与其他词的注意力权重，来捕捉词与词之间的关系。自注意力机制的计算公式如下：

$$\text{Attention}(Q,K,V)=\text{softmax}\left(\frac{Q{K}^T}{\sqrt{d_k}}\right)V$$

其中：

• ($Q$)：查询（Query）矩阵
• ($K$)：键（Key）矩阵
• ($V$)：值（Value）矩阵
• ($d_k$)：键向量的维度

关键点：自注意力机制是无序的，即对输入序列的位置信息不敏感。这意味着，如果我们将输入序列中的词语打乱顺序，自注意力机制的输出不会改变。这在NLP任务中是致命的，因为词语的顺序对语义至关重要。例如，“我喜欢你”和“你喜欢我”在自注意力机制中无法区分。

1.2 传统的位置编码方法

为了弥补这一缺陷，Vaswani等人在2017年的原始Transformer论文《Attention is All You Need》中提出了正弦/余弦位置编码（Sinusoidal Positional Encoding）。这种方法为每个位置生成一个固定的编码向量，并将其加到词嵌入向量上，以提供位置信息。

正弦/余弦位置编码的公式：

$$PE(pos,2i)=\sin \left(\frac{pos}{10000^{2i/d}}\right),PE(pos,2i+1)=\cos \left(\frac{pos}{10000^{2i/d}}\right)$$

其中：

• ($pos$)：词在序列中的位置
• ($i$)：嵌入向量的维度索引
• ($d$)：嵌入向量的总维度

工作原理：

• 通过不同频率的正弦和余弦函数，为每个位置生成独特的编码。
• 将 ( $PE$ ) 加到词嵌入向量 ( $\mathbf{x}$ ) 上：( ${\mathbf{x}}^{\prime }=\mathbf{x}+PE$ )。

局限性：

1. 绝对位置依赖：正弦/余弦位置编码主要捕捉绝对位置信息，难以直接反映词与词之间的相对位置关系。
2. 长序列表现欠佳：在处理超长序列时，位置编码的频率设计可能不适应，导致性能下降。
3. 静态编码：位置编码是预先计算好的，无法动态适应不同序列长度。

为了克服这些局限，研究者们探索了多种改进方案，其中ROPE作为一种创新的位置编码方法，凭借其旋转机制和对相对位置的自然捕捉，逐渐成为Transformer模型中的首选技术。

---

2. 旋转位置编码（ROPE）简介

ROPE（Rotary Position Embedding）最早由Su等人在2021年的论文《RoFormer: Enhanced Transformer with Rotary Position Embedding》中提出。它与传统的加性位置编码不同，通过旋转嵌入向量来编码位置信息，具体是将位置信息融入到自注意力机制的查询（Query）和键（Key）向量中。

2.1 ROPE的核心思想

ROPE的核心思想是：通过旋转矩阵对嵌入向量进行旋转，使得旋转后的向量在计算点积时，能够自然地反映出词与词之间的相对位置关系。

具体来说，ROPE将嵌入向量视为复数空间中的向量，并通过旋转矩阵对这些向量进行旋转。旋转的角度与词在序列中的位置成正比，从而将位置信息编码到向量的方向中。

独特优势：ROPE不仅能编码绝对位置，还能自然捕捉相对位置关系。例如，两个词之间的注意力得分会随着它们在序列中的相对距离而变化，这在处理长距离依赖时尤为重要。

---

3. ROPE的数学原理

要深入理解ROPE，我们需要从数学上推导其工作机制。

3.1 旋转矩阵与复数表示

ROPE的关键在于使用旋转矩阵来编码位置信息。首先，我们将嵌入向量的每一对维度（例如，第0和1维、第2和3维等）视为一个复数。假设嵌入维度 ($d$) 为偶数，我们可以将嵌入向量 ( $\mathbf{x}\in {\mathbb{R}}^d$ ) 分成 ( $d/2$ ) 个复数对。

例如，对于 ( $d=4$ )，嵌入向量 ( $\mathbf{x}=[x_0,x_1,x_2,x_3]$ ) 可以表示为两个复数：

• ( $x_0+x_{1}i$ )
• ( $x_2+x_{3}i$ )

在复数空间中，旋转一个复数 ( $z=a+bi$ ) 可以通过乘以 ( $e^{i\theta }$ ) 实现：

$$z^{\prime }=z\cdot e^{i\theta }=(a+bi)(\cos \theta +i\sin \theta )$$

展开后：

$$z^{\prime }=(a\cos \theta -b\sin \theta )+i(a\sin \theta +b\cos \theta )$$

这等价于对向量 ( $[a,b]$ ) 应用旋转矩阵：

$$\left[\begin{array}{cc}\cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{array}\right]\left[\begin{array}{c}a\\ b\end{array}\right]=\left[\begin{array}{c}a\cos \theta -b\sin \theta \\ a\sin \theta +b\cos \theta \end{array}\right]$$

因此，对于每个复数对（即每两个维度），我们可以通过旋转矩阵来编码位置信息。

3.2 ROPE的旋转角度设计

在ROPE中，旋转角度 ( $\theta$ ) 由词的位置 ( $p$ ) 和预定义的频率决定。具体地，对于第 ( $k$ ) 个复数对（对应第 ( $2k$ ) 和 ( $2k+1$ ) 维），旋转角度为：

$${\theta }_k=p\cdot {\omega }_k$$

其中，频率 ( ${\omega }_k$ ) 定义为：

$${\omega }_k=\frac{1}{10000^{2k/d}},k=0,1,\dots ,d/2-1$$

注意：这里的频率设计与正弦/余弦位置编码中的频率类似，但ROPE将其用于旋转角度，而不是直接加到嵌入向量上。

3.3 ROPE的旋转公式

对于位置 ( $p$ ) 处的嵌入向量 ( $\mathbf{x}\in {\mathbb{R}}^d$ )，ROPE通过以下方式旋转它：

$${\mathbf{x}}_{\text{rot}}=\left[\begin{array}{c}{x}_0\cos {\theta }_0-x_1\sin {\theta }_0\\ x_0\sin {\theta }_0+x_1\cos {\theta }_0\\ x_2\cos {\theta }_1-x_3\sin {\theta }_1\\ x_2\sin {\theta }_1+x_3\cos {\theta }_1\\ ⋮\\ x_{d-2}\cos {\theta }_{d/2-1}-x_{d-1}\sin {\theta }_{d/2-1}\\ x_{d-2}\sin {\theta }_{d/2-1}+x_{d-1}\cos {\theta }_{d/2-1}\end{array}\right]$$

简言之，ROPE对嵌入向量的每一对维度应用一个旋转矩阵，旋转角度由位置 ( $p$ ) 和频率 ( ${\omega }_k$ ) 决定。

3.4 ROPE在自注意力中的应用

在Transformer的自注意力机制中，ROPE被应用于查询（Query）和键（Key）向量。具体地，对于序列中的每个位置 ( $i$ )，我们旋转其查询向量 ( ${\mathbf{q}}_i$ ) 和键向量 ( ${\mathbf{k}}_i$ )：

$${\mathbf{q}}_i^{\text{rot}}=R(p_i){\mathbf{q}}_i,{\mathbf{k}}_i^{\text{rot}}=R(p_i){\mathbf{k}}_i$$

其中，( $R(p_i)$ ) 是位置 ( $p_i$ ) 对应的旋转矩阵。

然后，注意力得分计算为：

$${\text{score}}_{ij}=\frac{{\mathbf{q}}_i^{\text{rot}}\cdot {\mathbf{k}}_j^{\text{rot}}}{\sqrt{d_k}}$$

关键洞察：由于旋转矩阵的性质，( ${\mathbf{q}}_i^{\text{rot}}\cdot {\mathbf{k}}_j^{\text{rot}}$ ) 可以表示为原始向量 ( ${\mathbf{q}}_i$ ) 和 ( ${\mathbf{k}}_j$ ) 的点积，乘以一个与相对位置 ( $p_i-p_j$ ) 相关的旋转项。这使得注意力得分自然地依赖于词与词之间的相对位置。

---

4. ROPE的相对位置编码特性

ROPE的一个显著优势是它能够自然地捕捉相对位置信息。让我们通过数学推导来理解这一点。

4.1 旋转后的点积

考虑两个位置 ( $m$ ) 和 ( $n$ )，对应的旋转查询和键向量为：

$${\mathbf{q}}_m^{\text{rot}}=R(m){\mathbf{q}}_m,{\mathbf{k}}_n^{\text{rot}}=R(n){\mathbf{k}}_n$$

它们的点积为：

$${\mathbf{q}}_m^{\text{rot}}\cdot {\mathbf{k}}_n^{\text{rot}}=(R(m){\mathbf{q}}_m{)}^T(R(n){\mathbf{k}}_n)={\mathbf{q}}_m^{T}R(m{)}^{T}R(n){\mathbf{k}}_n$$

由于旋转矩阵 ( $R(\theta )$ ) 是正交矩阵，满足 ( $R(\theta {)}^T=R(-\theta )$ )，因此：

$$R(m{)}^{T}R(n)=R(-m)R(n)=R(n-m)$$

（后文有推导过程。）

所以：

$${\mathbf{q}}_m^{\text{rot}}\cdot {\mathbf{k}}_n^{\text{rot}}={\mathbf{q}}_m^{T}R(n-m){\mathbf{k}}_n$$

这表明，旋转后的点积等价于将键向量 ( ${\mathbf{k}}_n$ ) 旋转一个与相对位置 ( $n-m$ ) 相关的角度后，再与查询向量 ( ${\mathbf{q}}_m$ ) 做点积。

4.2 相对位置的自然捕捉

上述推导显示，注意力得分 ( ${\mathbf{q}}_m^{\text{rot}}\cdot {\mathbf{k}}_n^{\text{rot}}$ ) 直接依赖于相对位置 ( $n-m$ )，而不是绝对位置 ( $m$ ) 和 ( $n$ )。这使得ROPE在处理序列任务时，能够更自然地捕捉词与词之间的相对位置关系，尤其是在长序列中，相对位置信息比绝对位置信息更为重要。

简单示例：
假设 ( $d=2$ )，位置 ( $m=0$ )，( $n=1$ )，嵌入向量 ( ${\mathbf{q}}_0=[1,0]$ )，( ${\mathbf{k}}_1=[0,1]$ )，频率 ( ${\omega }_0=1$ )。那么：

• ( ${\theta }_0=0\cdot 1=0$ )，( $R(0)=\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]$ )
• (${\theta }_0=1\cdot 1=1$ )，( $R(1)=\left[\begin{array}{cc}\cos 1& -\sin 1\\ \sin 1& \cos 1\end{array}\right]$ )
• ( ${\mathbf{q}}_0^{\text{rot}}=[1,0]$ )
• ( ${\mathbf{k}}_1^{\text{rot}}=[\cos 1,\sin 1]$ )
• 点积：(${\mathbf{q}}_0^{\text{rot}}\cdot {\mathbf{k}}_1^{\text{rot}}=1\cdot \cos 1+0\cdot \sin 1=\cos 1$ )

如果位置变为 ( $m=1$ )，( $n=2$ )，相对位置仍为1，点积结果依然依赖 ( $\cos 1$ ) 和 ( $\sin 1$ )，体现出相对位置特性。

---

5. ROPE的代码实现

为了帮助读者更好地理解ROPE的实现细节，下面提供一个使用PyTorch的ROPE实现。

5.1 ROPE类定义

import torch
import torch.nn as nn

class RotaryEmbedding(nn.Module):
    def __init__(self, dim, max_position_embeddings=2048, base=10000, device=None):
        super(RotaryEmbedding, self).__init__()
        self.dim = dim
        self.max_position_embeddings = max_position_embeddings
        self.base = base
        # 计算频率
        inv_freq = 1.0 / (self.base ** (torch.arange(0, self.dim, 2).float().to(device) / self.dim))
        self.register_buffer("inv_freq", inv_freq)

    def forward(self, x, seq_len=None):
        if seq_len is None:
            seq_len = x.size(-2)
        # 生成位置索引
        position_ids = torch.arange(0, seq_len, dtype=torch.float, device=x.device).unsqueeze(-1)
        # 计算旋转角度
        angles = position_ids * self.inv_freq
        # 计算sin和cos
        sin = torch.sin(angles)
        cos = torch.cos(angles)
        # 扩展到与x相同的维度
        sin = sin.unsqueeze(0).repeat(x.size(0), 1, 1)
        cos = cos.unsqueeze(0).repeat(x.size(0), 1, 1)
        # 应用旋转
        x_rot = x.clone()
        for i in range(0, self.dim, 2):
            x_rot[..., i] = x[..., i] * cos[..., i // 2] - x[..., i + 1] * sin[..., i // 2]
            x_rot[..., i + 1] = x[..., i] * sin[..., i // 2] + x[..., i + 1] * cos[..., i // 2]
        return x_rot

代码说明：

• dim：嵌入维度，必须为偶数。
• max_position_embeddings：支持的最大序列长度。
• base：频率计算中的基数，默认值为10000。
• inv_freq：预计算的频率，用于生成旋转角度。
• forward：对输入张量 ( $x$ ) 应用旋转，返回旋转后的张量。

5.2 使用示例

假设我们有一个批次大小为1、序列长度为3、嵌入维度为4的输入：

dim = 4
seq_len = 3
x = torch.randn(1, seq_len, dim)  # 形状: [batch_size, seq_len, dim]
rotary_emb = RotaryEmbedding(dim)
x_rot = rotary_emb(x)
print("Rotated embeddings:", x_rot)

在Transformer模型中，我们通常在自注意力层的查询和键向量上应用ROPE：

import math
# 假设 q 和 k 是查询和键向量，形状: [batch_size, num_heads, seq_len, head_dim]
batch_size, num_heads, seq_len, head_dim = 2, 8, 10, 64
q = torch.randn(batch_size, num_heads, seq_len, head_dim)
k = torch.randn(batch_size, num_heads, seq_len, head_dim)
rotary_emb = RotaryEmbedding(head_dim)
q_rot = rotary_emb(q)
k_rot = rotary_emb(k)
# 计算注意力得分
scores = torch.einsum("bhqd,bhkd->bhqk", q_rot, k_rot) / math.sqrt(head_dim)
print("Attention scores shape:", scores.shape)  # [batch_size, num_heads, seq_len, seq_len]

---

6. ROPE的优势与局限性

6.1 优势

1. 自然捕捉相对位置信息：通过旋转机制，注意力得分直接依赖于相对位置，适合处理长距离依赖。
2. 对序列长度鲁棒：ROPE不需要为特定序列长度预先计算位置编码，可以灵活适应不同长度的输入。
3. 计算效率高：旋转操作可以通过向量化实现，计算开销小。
4. 与Transformer架构兼容：ROPE可以无缝集成到现有的Transformer模型中，无需大幅修改架构。

6.2 局限性

1. 嵌入维度要求：ROPE要求嵌入维度为偶数，这在某些模型设计中可能需要调整。
2. 超长序列下的频率调整：在处理极长序列时，可能需要调整频率参数（如NTK-aware scaled ROPE），以避免旋转角度过大导致的信息丢失。
3. 实现复杂度：相较于简单的加性位置编码，ROPE的实现稍显复杂，需要仔细处理旋转矩阵的计算。

---

7. ROPE在现代大型语言模型中的应用

ROPE因其优异的性能，已被广泛应用于现代大型语言模型中，例如：

• Llama：Meta的Llama模型系列采用了ROPE作为位置编码方案。
• PaLM：Google的PaLM模型使用ROPE处理长序列任务。
• RoBERTa：Fairseq中的RoBERTa实现也集成了ROPE。

在实际应用中，ROPE在文本生成、机器翻译、长文本总结等领域展现了强大的能力，尤其是在处理长序列和捕捉复杂语义关系时，优势尤为明显。

---

8. 总结

旋转位置编码（ROPE）作为一种创新的位置编码方法，通过旋转嵌入向量来编码位置信息，成功解决了传统位置编码在捕捉相对位置信息和处理长序列方面的不足。其数学原理巧妙地利用了旋转矩阵的性质，使得注意力机制能够自然地融入相对位置信息。ROPE的高效实现和与Transformer架构的良好兼容性，使其成为现代NLP模型的首选位置编码方案。

通过本文的详细介绍，从基础概念到数学推导，再到代码实现和实际应用，相信读者已经对ROPE有了深入的理解。ROPE不仅是理论上的突破，也是实践中的利器，在推动NLP技术发展的道路上扮演着重要角色。

---

补充：$R(-m)R(n)=R(n-m)$推导

我们来一步步推导 ( $R(m{)}^{T}R(n)=R(-m)R(n)=R(n-m)$ ) 这个等式，深入理解旋转位置编码（ROPE）中旋转矩阵的性质。这部分推导涉及到旋转矩阵的正交性和旋转角度的加法性质，是ROPE能够捕捉相对位置信息的关键数学基础。

---

1. 旋转矩阵的定义

在ROPE中，旋转矩阵 ( $R(p)$ ) 是为序列中位置 ( $p$ ) 定义的，用于旋转嵌入向量。假设嵌入维度 ( $d$ ) 为偶数，( $R(p)$ ) 是分块对角矩阵，每个块是一个 2×2 的旋转矩阵，对应嵌入向量的每一对维度（例如，第0和1维、第2和3维）。对于第 ( $k$ ) 个维度对，旋转矩阵定义为：

$$R_k(p)=\left[\begin{array}{cc}\cos (p\cdot {\omega }_k)& -\sin (p\cdot {\omega }_k)\\ \sin (p\cdot {\omega }_k)& \cos (p\cdot {\omega }_k)\end{array}\right]$$

其中：

• ( $p$ ) 是位置索引（整数）。
• ( ${\omega }_k=\frac{1}{10000^{2k/d}}$ ) 是第 ( $k$ ) 个维度对的频率，( $k=0,1,\dots ,d/2-1$ )。

对于整个 ( $d$ ) 维嵌入向量，( $R(p)$ ) 是 ( $d\times d$ ) 的分块对角矩阵：

$$R(p)=\text{diag}(R_0(p),R_1(p),\dots ,R_{d/2-1}(p))$$

但在推导中，我们只需关注单个 2×2 旋转矩阵的性质，因为整体矩阵的运算可以分解为每个块的独立运算。

---

2. 旋转矩阵的正交性

旋转矩阵的一个重要性质是正交性，即：

$$R(p{)}^{T}R(p)=I$$

其中 ( $I$ ) 是单位矩阵。这意味着 ( $R(p{)}^T$ ) 是 ( $R(p)$ ) 的逆矩阵：

$$R(p{)}^T=R(p{)}^{-1}$$

对于 2×2 旋转矩阵：

$$R(p)=\left[\begin{array}{cc}\cos (p{\omega }_k)& -\sin (p{\omega }_k)\\ \sin (p{\omega }_k)& \cos (p{\omega }_k)\end{array}\right]$$

其转置为：

$$R(p{)}^T=\left[\begin{array}{cc}\cos (p{\omega }_k)& \sin (p{\omega }_k)\\ -\sin (p{\omega }_k)& \cos (p{\omega }_k)\end{array}\right]$$

我们可以验证：

$$R(p{)}^{T}R(p)=\left[\begin{array}{cc}\cos (p{\omega }_k)& \sin (p{\omega }_k)\\ -\sin (p{\omega }_k)& \cos (p{\omega }_k)\end{array}\right]\left[\begin{array}{cc}\cos (p{\omega }_k)& -\sin (p{\omega }_k)\\ \sin (p{\omega }_k)& \cos (p{\omega }_k)\end{array}\right]$$

计算第一行第一列：

$$\cos (p{\omega }_k)\cdot \cos (p{\omega }_k)+\sin (p{\omega }_k)\cdot \sin (p{\omega }_k)={\cos }^2(p{\omega }_k)+{\sin }^2(p{\omega }_k)=1$$

计算完整矩阵：

• 第一行第二列：( $\cos (p{\omega }_k)\cdot (-\sin (p{\omega }_k))+\sin (p{\omega }_k)\cdot \cos (p{\omega }_k)=-\cos (p{\omega }_k)\sin (p{\omega }_k)+\sin (p{\omega }_k)\cos (p{\omega }_k)=0$ )
• 第二行第一列：( $(-\sin (p{\omega }_k))\cdot \cos (p{\omega }_k)+\cos (p{\omega }_k)\cdot \sin (p{\omega }_k)=-\sin (p{\omega }_k)\cos (p{\omega }_k)+\cos (p{\omega }_k)\sin (p{\omega }_k)=0$ )
• 第二行第二列：( $(-\sin (p{\omega }_k))\cdot (-\sin (p{\omega }_k))+\cos (p{\omega }_k)\cdot \cos (p{\omega }_k)={\sin }^2(p{\omega }_k)+{\cos }^2(p{\omega }_k)=1$ )

结果为：

$$R(p{)}^{T}R(p)=\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]=I$$

这确认了 ( $R(p)$ ) 是正交矩阵。

---

3. 旋转矩阵的逆与负角度

旋转矩阵的逆可以通过旋转角度的负值表示。考虑：

$$R(-p)=\left[\begin{array}{cc}\cos (-p{\omega }_k)& -\sin (-p{\omega }_k)\\ \sin (-p{\omega }_k)& \cos (-p{\omega }_k)\end{array}\right]$$

由于 ( $\cos (-\theta )=\cos (\theta )$ ) 和 ( $\sin (-\theta )=-\sin (\theta )$ )：

$$R(-p)=\left[\begin{array}{cc}\cos (p{\omega }_k)& \sin (p{\omega }_k)\\ -\sin (p{\omega }_k)& \cos (p{\omega }_k)\end{array}\right]$$

这与 ( $R(p{)}^T$ ) 完全相同，因此：

$$R(p{)}^T=R(-p)$$

这表明，旋转矩阵的转置等价于反向旋转，即旋转角度取负。

---

4. 推导 ( $R(m{)}^{T}R(n)=R(-m)R(n)$ )

现在推导第一个等式 ( $R(m{)}^{T}R(n)=R(-m)R(n)$ )。

根据旋转矩阵的正交性：

$$R(m{)}^T=R(-m)$$

所以：

$$R(m{)}^{T}R(n)=R(-m)R(n)$$

这个等式是直接成立的，因为 ( $R(m{)}^T$ ) 定义为 ( $R(-m)$ )。我们可以进一步验证：

• ( $R(m)=\left[\begin{array}{cc}\cos (m{\omega }_k)& -\sin (m{\omega }_k)\\ \sin (m{\omega }_k)& \cos (m{\omega }_k)\end{array}\right]$ )
• ( $R(m{)}^T=\left[\begin{array}{cc}\cos (m{\omega }_k)& \sin (m{\omega }_k)\\ -\sin (m{\omega }_k)& \cos (m{\omega }_k)\end{array}\right]=R(-m)$ )
• ( $R(n)=\left[\begin{array}{cc}\cos (n{\omega }_k)& -\sin (n{\omega }_k)\\ \sin (n{\omega }_k)& \cos (n{\omega }_k)\end{array}\right]$ )

左侧 ( $R(m{)}^{T}R(n)$ ) 和右侧 ( $R(-m)R(n)$ ) 是同一个矩阵乘法，因此：

$$R(m{)}^{T}R(n)=R(-m)R(n)$$

---

5. 推导 ( $R(-m)R(n)=R(n-m)$ )

接下来推导第二个等式 ( $R(-m)R(n)=R(n-m)$ )，这涉及到旋转矩阵的角度加法性质。

计算 ( $R(-m)R(n)$ )：

$$R(-m)=\left[\begin{array}{cc}\cos (-m{\omega }_k)& -\sin (-m{\omega }_k)\\ \sin (-m{\omega }_k)& \cos (-m{\omega }_k)\end{array}\right]=\left[\begin{array}{cc}\cos (m{\omega }_k)& \sin (m{\omega }_k)\\ -\sin (m{\omega }_k)& \cos (m{\omega }_k)\end{array}\right]$$

$$R(n)=\left[\begin{array}{cc}\cos (n{\omega }_k)& -\sin (n{\omega }_k)\\ \sin (n{\omega }_k)& \cos (n{\omega }_k)\end{array}\right]$$

矩阵乘法 ( $R(-m)R(n)$ )：

• 第一行第一列：( $\cos (m{\omega }_k)\cdot \cos (n{\omega }_k)+\sin (m{\omega }_k)\cdot \sin (n{\omega }_k)$ )
• 第一行第二列：( $\cos (m{\omega }_k)\cdot (-\sin (n{\omega }_k))+\sin (m{\omega }_k)\cdot \cos (n{\omega }_k)$ )
• 第二行第一列：( $(-\sin (m{\omega }_k))\cdot \cos (n{\omega }_k)+\cos (m{\omega }_k)\cdot \sin (n{\omega }_k)$ )
• 第二行第二列：( $(-\sin (m{\omega }_k))\cdot (-\sin (n{\omega }_k))+\cos (m{\omega }_k)\cdot \cos (n{\omega }_k)$ )

使用三角恒等式：

• ( $\cos (a)\cos (b)+\sin (a)\sin (b)=\cos (a-b)$ )
• ( $\cos (a)(-\sin (b))+\sin (a)\cos (b)=-\sin (a)\cos (b)+\sin (a)\cos (b)=-\sin (a-b)$ )
• ( $-\sin (a)\cos (b)+\cos (a)\sin (b)=\sin (a-b)$ )

代入 ( $a=m{\omega }_k$ )，( $b=n{\omega }_k$ )：

• 第一行第一列：( $\cos (m{\omega }_k-n{\omega }_k)=\cos ((n-m){\omega }_k)$ )
• 第一行第二列：( $-\sin (m{\omega }_k-n{\omega }_k)=-\sin ((n-m){\omega }_k)$ )
• 第二行第一列：( $\sin (m{\omega }_k-n{\omega }_k)=\sin ((n-m){\omega }_k)$ )
• 第二行第二列：( $\cos (m{\omega }_k-n{\omega }_k)=\cos ((n-m){\omega }_k)$ )

因此：

$$R(-m)R(n)=\left[\begin{array}{cc}\cos ((n-m){\omega }_k)& -\sin ((n-m){\omega }_k)\\ \sin ((n-m){\omega }_k)& \cos ((n-m){\omega }_k)\end{array}\right]=R(n-m)$$

---

6. 完整推导结论

结合两步：

$$R(m{)}^{T}R(n)=R(-m)R(n)=R(n-m)$$

这个等式表明，两个位置 ( $m$ ) 和 ( $n$ ) 的旋转矩阵相乘后，结果仅依赖于它们的相对位置 ( $n-m$ )。这正是ROPE能够捕捉相对位置信息的核心数学依据。

---

7. 在ROPE中的意义

在自注意力机制中，查询向量 ( ${\mathbf{q}}_m$ ) 和键向量 ( ${\mathbf{k}}_n$ ) 被旋转为：

$${\mathbf{q}}_m^{\text{rot}}=R(m){\mathbf{q}}_m,{\mathbf{k}}_n^{\text{rot}}=R(n){\mathbf{k}}_n$$

点积为：

$${\mathbf{q}}_m^{\text{rot}}\cdot {\mathbf{k}}_n^{\text{rot}}=(R(m){\mathbf{q}}_m{)}^T(R(n){\mathbf{k}}_n)={\mathbf{q}}_m^{T}R(m{)}^{T}R(n){\mathbf{k}}_n={\mathbf{q}}_m^{T}R(n-m){\mathbf{k}}_n$$

这表明注意力得分依赖于相对位置 ( $n-m$ )，而不是绝对位置 ( $m$ ) 和 ( $n$ )，从而实现了ROPE的相对位置编码特性。

---

总结

通过旋转矩阵的正交性（( $R(m{)}^T=R(-m)$ )）和角度加法性质（( $R(-m)R(n)=R(n-m)$ )），我们推导出了 ( $R(m{)}^{T}R(n)=R(-m)R(n)=R(n-m)$ )。这个性质是ROPE在Transformer模型中高效编码相对位置的数学基础。

参考资料

1. Su, J., et al. (2021). “RoFormer: Enhanced Transformer with Rotary Position Embedding.” arXiv:2104.09864
2. Vaswani, A., et al. (2017). “Attention is All You Need.” arXiv:1706.03762
3. Rotary Embeddings Explained | Papers With Code
4. GitHub - lucidrains/rotary-embedding-torch
5. Understanding Rotary Positional Encoding | Medium

---

后记

2025年2月23日20点19分于上海，在Grok 3 大模型辅助下完成。