二维码笔记系列:

本篇是学习编码理论时写的杂记,未做整理,内容可能枯燥且冗余。

唠唠闲话

二维码是一种存储二进制数据的矩阵条码,其使用 RS(Reed-Solomon) 编码进行纠错。RS 码是一种线性码(Linear code),其在现实生活中应用广泛,且存在高效的编解码算法。

本篇介绍编码原理,纠错算法以及拓展知识,跳转目录:


分组码

本节介绍分组码的定义例子,以及汉明距离等重要概念。

概述

在信息传输中,因传输通道噪音的干扰,接收端收到的信息与发送端发出的信息可能有误差。纠错码的引入,通过增加额外数据位,使传输的数据抵抗干扰的能力更加鲁棒(robust)。

编码理论中,分组码(block code)是一大类重要的纠错码,它以块的形式对数据进行编码。分组码的抽象定义在概念上很有用,因为它允许编码理论家、数学家和计算机科学家以统一的方式研究所有分组码的局限性。这种局限性通常对编码界限进行研究,研究其与分组码不同参数的联系,例如信息率与检测纠正错误能力。

定义及例子

纠错码用于在受信道噪声影响的不可靠通信信道上可靠地传输数字数据。当发送者想要使用分组码传输长数据流时,发送者将流分成一些固定大小的片段,每个这样的片段称为消息(message)。分组码过程中,将每 message 单独编码成块(block)。然后,发送者将所有 block 传输给接收者,接收者又可以使用某种解码机制(希望)从可能损坏的接收 block 中恢复原始消息。整体传输的性能和成功取决于信道和分组码的参数。

  1. Σ\Sigma 为长度为 qq 的字符集,Σn\Sigma^n 为 n 元字符构成的集合

  2. 将子集 CΣnC\subseteq\Sigma^n 称为块长(block length)为 nn分组码

  3. CC消息长度或维数定义为 k=logq(C)k=log_q(|C|)

  4. CC编码率或信息率(code rate/information rate) 定义为

    R=kn=logq(C)nR=\frac{k}{n}=\frac{log_q(|C|)}{n}

    注:RR = 消息长度/编码块长度

  5. CC 的元素 xx 称为码字(codeword)

  6. 设分组码 CC 的距离为 dd,将 CC 简称为 (n,k,d)Σ(n, k, d)_{\Sigma} 分组码,距离定义在下边给出

  7. 分组码可以理解为嵌入映射 CΣnC\hookrightarrow\Sigma^n,映射左侧为原始信息,右侧为加密后的信息

  8. 例1:{hi, am, ok}\{hi,\ am,\ ok\} 为字符集 Σ={az}\Sigma=\{a-z\} 的块长为 2 维数为 log263log_{26}3 的分组码

  9. 例2:{(1,1,1),(0,0,0)}\{(1,1,1), (0,0,0)\}Σ={0,1}\Sigma=\{0,1\} 的块长 3 维数 1 的分组码,写成映射形式

    Σ1Σ30(0,0,0)1(1,1,1)\begin{align*} \Sigma^1&\rightarrow\Sigma^3\\ 0&\rightarrow(0, 0, 0)\\ 1&\rightarrow(1, 1, 1) \end{align*}

  10. 用大空间 Σ3\Sigma^3 表示小空间 Σ1\Sigma^1 的数据,能更好地处理错误信息:

  • 检测错误(detect error):假设传输 (0,0,0)(0,0,0) 时,发生了至多 2 个错误,比如

    (0,0,0)(0,1,1)(0,0,0)\rightarrow(0,1,1)

    接收端判断发现 (0,1,1)Σ3(0,1,1)\notin\Sigma^3,立刻检测到错误
  • 填补缺失(correct erasure):假设传输 (0,0,0)(0,0,0) 时,丢失了至多 2 个位置的数据(丢失位置已知),比如

    (0,0,0)(,0,)(0,0,0)\rightarrow(*,0,*)

    接收端通过余下数据仍能还原初始数据
  • 纠错能力(correct error):假设传输 (0,0,0)(0,0,0) 时,有至多 1 个位置的数据错误,比如

    (0,0,0)(0,1,0)(0,0,0)\rightarrow(0,1,0)

    由于至多一个错误位置,初始数据只可能是 (0,0,0)(0,0,0)
  1. 显然,编码的纠错能力是有限的,假设数据 (0,0,0)(0,0,0) 传送时发生 2 个位置的错误,比如变成 (0,1,1)(0,1,1),接收端并不能确定原数据是 (0,0,0)(0,0,0) 还是 (1,1,1)(1,1,1)。类似地,检测错误和填补缺失的能力也是有限的,假设传送时发生 3 个位置的错误,这种错误显然不能被检测出来。

汉明距离

下边引入距离定义,用来衡量处理错误的能力,不妨设字符集 Σ\Sigma 包含数字 0

  1. 设码字 c=(c1,,cn)Σnc=(c_1,\cdots,c_n)\in\Sigma^n,定义权重(weight) wt(c)=#{ici0}wt(c) = \#\{i|c_i\neq 0\}cc 的非零成分的数目。

  2. c,cΣnc,c'\in\Sigma^n,定义 cccc'汉明距离(Hamming distance)

    Δ(c,c)=#{icici}\Delta(c,c') = \#\{i|c_i\neq c_i'\}

    注:容易验证 Σn\Sigma^n 关于 Δ\Delta 作成度量空间

  3. 编码 CC 的(最小)距离定义为

    δ(C)=min{Δ(c,c)  ccC}\delta(C) = \min\{\Delta(c,c')\ |\ c\neq c'\in C\}

    距离 d=δ(C)d=\delta(C) 描述了码字 cc 转变为另一码字 cc' 的最少修改数

  4. 定理)容易证明,分组码 (n,k,d)Σ(n,k,d)_{\Sigma} 的错误处理能力为:

    • 编码 CC 能检测至多 ed1e\leq d-1 个错误
    • 编码 CC 能填补至多 ed1e\leq d-1 个缺失
    • 编码 CC 能纠正至多 ed12e\leq \lfloor\frac{d-1}{2}\rfloor 个错误,即 2ed12e\leq d-1
  5. 从图像理解汉明距离,以码字为圆心,半径 ed12e\leq \lfloor\frac{d-1}{2}\rfloor 的“小球”互不相交
    20220504115713

衡量指标

编码在设计上有几个重要目标,或者说衡量指标:

  • 处理错误的能力,比如检测错误,纠正错误
  • 尽量小的冗余信息(minimize overhead)
  • 计算效率,比如编码,解码和纠错等操作的执行效率

在分组码 (n,k,d)Σ(n,k,d)_{\Sigma} 中,错误处理能力用距离 dd相对距离 dd' 描述,冗余信息用编码率/信息率 RR 表述。

一般地,编码距离 dd 越大,错误处理能力越强,但相应地,冗余信息将越多,这导致了信息率 RR 的减小。因此,协调错误处理能力 dd 与信息率 RR 是编码领域的重要问题,然而,如何制定最佳权衡策略(best trade-off)还是个开问题。

尽管如此,有不少研究在某些方面给出了解答,例如下边要介绍的汉明界。

有时也用相对距离 d=dnd'=\frac{d}{n} 来衡量编码的纠错能力。

汉明界

在编码理论领域,汉明界(Hamming Bound) 是对分组码参数的限制,它对编码空间的利用率提出了重要限制。汉明界也被称为球包装界体积界,这可以从推导过程看出来。

  1. xΣnx\in\Sigma^n,如下定义 BΣn(x,e)B_{\Sigma^n}(x, e) ,并称为 xx 的以 ee 为半径的 Hamming Ball。

    BΣn(x,e)={yΣn  Δ(x,y)e}B_{\Sigma^n}(x, e) = \{y\in\Sigma^n\ |\ \Delta(x,y)\leq e\}

  2. 将汉明球 BΣn(x,e)B_{\Sigma^n}(x, e) 包含的 Σn\Sigma^n 元素个数称为球的体积(Volume) ,记为

    Volq(e,n)=BΣn(x,e), where q=ΣVol_q(e,n)=|B_{\Sigma^n}(x, e)|,\ where\ q=|\Sigma|

    由于对称性,球体积与 xx 的选取无关,因而体积定义中不含参数 xx

  3. 推导易得下边公式

    Volq(e,n)=i=0e(ni)(q1)iVol_q(e,n) = \sum_{i=0}^e\begin{pmatrix}n\\ i\end{pmatrix}(q-1)^i

    其中加项 (ni)(q1)i\begin{pmatrix}n\\ i\end{pmatrix}(q-1)^i 代表与 xx 的汉明距离为 iiΣn\Sigma^n 中的元素数目。换言之,这些元素可通过 xx 修改 i 个位置数据得到。

  4. ddCC 的距离,同前边示意图易知,汉明球 BΣn(x,d12),xCB_{\Sigma^n}(x,\lfloor\frac{d-1}{2}\rfloor),x\in C 之间两两不交,继而有

    CVolq(d12,n)Σnqn|C|\cdot Vol_q(\lfloor\frac{d-1}{2}\rfloor,n) \leq |\Sigma^n| \leq q^n

    取等当且仅当这些汉明球覆盖了整个空间;整理可得

    logq(C)+logq(Volq(d12,n))n k+logq(Volq(d12,n))n kn+logq(Volq(d12,n))n1 R1logq(Volq(d12,n))n\begin{align*} &\log_q(|C|) + \log_q(Vol_q(\lfloor\frac{d-1}{2}\rfloor,n))\leq n\\ \Rightarrow &\ k+\log_q(Vol_q(\lfloor\frac{d-1}{2}\rfloor,n))\leq n\\ \Rightarrow &\ \frac{k}{n}+\frac{\log_q(Vol_q(\lfloor\frac{d-1}{2}\rfloor,n))}{n}\leq 1\\ \Rightarrow &\ R \leq 1-\frac{\log_q(Vol_q(\lfloor\frac{d-1}{2}\rfloor,n))}{n} \end{align*}

    上式表明,当编码块长 nn,字符集长 qq 以及编码距离 dd 确定时,信息率 RR 存在不等式右侧给出的上界,这一上界便称为汉明界

  5. 汉明界在某些角度上回答了编码距离与信息率的关系,假设编码块长 nn 和字符集 Σ\Sigma 确定,汉明界表明:当编码距离 dd 确定时,信息率 RR 不会超过某个数值。而取到边界当且仅当所构造编码的半径为 d12\lfloor\frac{d-1}{2}\rfloor 的汉明球覆盖了整个空间。达到汉明界的编码被称为完美编码(perfect code)汉明码(Hamming code)

  6. 举个例子,考虑二进制字符集 Σ={0,1}\Sigma=\{0,1\},设编码的块长为 n=7n=7,距离为 d=3d=3,维数为 kk,计算汉明界

    R=kn1logq(Volq(d12,n))n=1log2(Vol2(1,7))7, Vol2(1,7)=1+7=1log2(8)7=47\begin{align*} R = \frac{k}{n}&\leq 1-\frac{\log_q(Vol_q(\lfloor\frac{d-1}{2}\rfloor,n))}{n}\\ &= 1-\frac{\log_2(Vol_2(1,7))}{7},\ Vol_2(1,7) = 1+7\\ &= 1-\frac{\log_2(8)}{7} = \frac{4}{7} \end{align*}

    根据汉明界,编码维数 k4k\leq 4,如下编码可以取到等号

    C:Σ4Σ7(x1x2x3x4)(1111011110111101)(x1x2x3x4)=(x1x2x3x4x2+x3+x4x1+x3+x4x1+x2+x4)\begin{align*} C:\Sigma^4&\rightarrow \Sigma^7\\ \begin{pmatrix}x_1\\ x_2\\ x_3\\ x_4\end{pmatrix}& \mapsto\begin{pmatrix} 1& & & \\ & 1& & \\ & &1 & \\ & & &1 \\ 0& 1& 1 &1 \\ 1& 0 &1 &1 \\ 1& 1 & 0 &1 \\ \end{pmatrix} \begin{pmatrix}x_1\\ x_2\\ x_3\\ x_4\end{pmatrix} = \begin{pmatrix}x_1\\ x_2\\ x_3\\ x_4\\ x_2+x_3+x_4 \\ x_1+x_3+x_4 \\ x_1+x_2+x_4\end{pmatrix} \end{align*}

例中,编码用线性变换定义,这种编码方式称为线性码,Reed Solomon 编码便是一类特殊的线性码。


线性码

线性码是一种纠错码,任何码字的线性组合也是一个码字。线性码传统上分为分组码和卷积码,以下仅讨论分组码。

下设 F=Fq\mathbb{F} = \mathbb{F}_qqq 个元素构成的有限域。

基本概念

  1. (定义) 编码空间 Fn\mathbb{F}^n 的真子空间 C=FkC=\mathbb{F}^k 称为 F\mathbb{F} 上的维数为 kk,块长为 nn线性码(Linear Code)。将矩阵 GFn×kG\in\mathbb{F}^{n\times k} 称为 CC生成元矩阵(Generator Matrix),若

    C=column span(G)={Gx  xFk}C=column\ span(G)=\{G\cdot x\ |\ x\in\mathbb{F}^k\}

    借助生成矩阵 GG,线性码 CC 可写为映射形式

    C:FkFnXGX, where X=(x1x2xk)\begin{align*} C:\mathbb{F}^k&\rightarrow\mathbb{F}^n\\ X&\mapsto GX,\ where\ X=\begin{pmatrix} x_1\\ x_2\\\vdots\\ x_k \end{pmatrix} \end{align*}

  2. 生成元矩阵不唯一:任意可逆矩阵左乘生成元矩阵,仍得到生成元矩阵。特别地,通过调整 Fn\mathbb{F}^n 的索引,总能左乘得到上分块为单位阵的生成元矩阵:

    C:FkFnX(IG)X=(XGX)\begin{align*} C:\mathbb{F}^k&\rightarrow\mathbb{F}^n\\ X &\mapsto \begin{pmatrix}I\\ G'\end{pmatrix} X=\begin{pmatrix} X\\G'X \end{pmatrix} \end{align*}

    容易发现,初始信息 XX 嵌在编码后的信息中。我们将满足这一性质的编码称为系统码(Systematic Code)

  3. HF(nk)×nH\in\mathbb{F}^{(n-k)\times n} 为线性码 CC奇偶检验矩阵(Parity Check Matrix),若

    ker(H)={xFnHx=0}={GyFnyFk}=im(G)\begin{align*} \ker(H) &= \{x\in\mathbb{F}^n|Hx=0\}\\ &= \{Gy\in\mathbb{F}^n| y\in\mathbb{F}^k\}=im(G) \end{align*}

    由线性代数知识,易见 rank(H)=nk,rank(G)=krank(H)=n-k, rank(G)=k,即矩阵 GG 为列满秩,矩阵HH 为行满秩

  4. (C,G,H)(C,G,H) 为一组线性码,则 (C,Ht,Gt)(C^\perp, H^t, G^t) 也为一组线性码,称为 CC对偶码(Dual Code),其中

    C={xFnx,y=xty=0, yC}C^{\perp} = \{x\in\mathbb{F}^n|\lang x,y\rang=x^ty=0,\forall\ y\in C\}

注:编码可以理解为信息空间到表示空间的单射,线性码则是通过列满秩矩阵给出这个单射。

编码距离

相比分组码,线性码的编码距离可以用矩阵 GGHH 计算。

  1. 以下命题等价

    1. 编码 CC 的汉明距离为 dd
    2. GG 列空间非零向量最小权重,即

      min{wt(x)xcolumn span of(G)}\min\{wt(x) | x\in column\ span\ of(G)\}

    3. HH 的任意 d1d-1 个列向量线性无关,且存在长为 dd 的线性相关列向量组;换言之,HH 的线性相关列向量组的最小长度为 dd
    4. GG 的任意 nd+1n-d+1 个行向量构成的向量组中,存在 kk 个线性无关向量;换言之,GG 的任意 nd+1n-d+1 行构成的矩阵列满秩
  2. 我们依次证明这几个命题

    • (1)(2)(1)\Leftrightarrow (2)
      根据线性性:

    d=δ(c,c)=δ(cc,0)=wt(cc)=wt(c),cCd=\delta(c,c')=\delta(c-c', 0) = wt(c-c')=wt(c''),c''\in C

    • (2)(3)(2)\Leftrightarrow(3)
      先将 HH 列分块化,并由 ker(H)=C\ker(H)=C

      H=(H1Hn)x=(x1,,xn)tCHx=Oi=1nxiHi=O\begin{align*} H = \begin{pmatrix} H_1&\cdots&H_n \end{pmatrix}\\ x=(x_1,\cdots,x_n)^t\in C &\Leftrightarrow Hx = O\\ &\Leftrightarrow \sum_{i=1}^nx_iH_i=O \end{align*}

    xCx\in Crr 个非零元数目对应得到 HHrr 个线性相关的列向量

    • (1)(4)(1)\Leftrightarrow(4)
      根据定义:CC 的编码距离为 dd 等价于 yC\forall y\in Cyy 的任意 d1d-1 项出错不会得到另一个码字;等价于方程组 Gx=yGx=y 中,将 G,yG,y 的任意 d1d-1 个行向量丢弃得到的方程 Gx=yG'x=y',其只有唯一解 xx;等价于 GG 去掉任意 d1d-1 行,得到的矩阵列满秩
  3. 一般地,选取列向量 xCx\in C 得距离上界 dwt(x)d\leq wt(x),选取 HHrr 个线性相关向量得到距离下界 drd\geq r,比如前边例子

    G=(1111011110111101),H=(011110010110101101001)G = \begin{pmatrix} 1& & & \\ & 1& & \\ & &1 & \\ & & &1 \\ 0& 1& 1 &1 \\ 1& 0 &1 &1 \\ 1& 1 & 0 &1 \\ \end{pmatrix},H=\begin{pmatrix} 0&1&1&1&1&0&0\\ 1&0&1&1&0&1&0\\ 1&1&0&1&0&0&1 \end{pmatrix}

    GG 的第一列非零数得 d3d\leq 3,由 HH 前三列向量线性相关性 d3d\geq 3,于是 d=3d=3


信息率

前边提到衡量编码的三个重要指标:处理错误能力尽量小的冗余信息,以及计算效率。关于错误处理能力,我们用编码距离 dd 来衡量,对于线性码还可以借助生成元矩阵 GG 和奇偶检验矩阵 HH 做相关计算;对于冗余信息,Hamming Bound 给出了某些方面的回答,本节将对信息率的边界做更深入的探讨。

本节侧重介绍辛格尔顿界,其与 Reed Solomon 纠错直接关联,其他作为拓展兴趣只做科普性地介绍,不讨论证明。

辛格尔顿界

在编码理论中,辛格尔顿界(Singleton Bound) 以 Richard Collom Singleton 命名,其表明 (n,k,d)q(n,k,d)_q 分组码中,信息长满足 knd+1k\leq n-d+1,继而

R1d1nR\leq 1-\frac{d-1}{n}

这个定理的证明非常简单,考虑复合映射 φ\varphi 为单射

φ:ΣkΣnΣnd+1c(c1,,cn)(cd,,cn)\begin{align*} \varphi:\Sigma^k&\rightarrow\Sigma^n\rightarrow\Sigma^{n-d+1}\\ c&\mapsto(c_1,\cdots, c_{n})\mapsto(c_{d},\cdots,c_n) \end{align*}

φ\varphi 必为单射,由于 cccc' 至多前 d-1 个位置不同,即

φ(c)=φ(c)(cd,,cn)=(cd,,cn)Δ(c,c)d1c=c\begin{align*} \varphi(c)=\varphi(c')&\Rightarrow(c_{d},\cdots,c_n)=(c_{d}',\cdots,c_n')\\ &\Rightarrow \Delta(c,c')\leq d-1\\ &\Rightarrow c=c' \end{align*}

φ\varphi 为单射,比较集合立得 knd+1k\leq n-d+1

对于线性码,从上一节编码距离的刻画(4) 立知,nd+1kn-d+1\geq k

GV 界

在编码理论中,GV Bound(Gilbert–Varshamov Bound) 是由 Edgar Gilbert 和 Rom Varshamov 分别独立提出的,关于编码(未必线性)参数的边界,Varshamov 通过使用线性码的概率方法证明了这一界限。有关该证明的更多信息请参阅这里

定理 (Gilbert–Varshamov) 设 qq 为素数幂次,dnd\leq n,则存在线性码 CC 使得

  • CC 的块长为 nn,字符集数目为 qq, 汉明边界为 dd
  • CC 的信息率 R=knR=\frac{k}{n} 满足下边性质

R1logq(Volq(d1,n))1n\begin{align*} R\geq 1- \frac{log_q(Vol_q(d-1, n))-1}{n} \end{align*}

定理证明过程并不复杂,但使用了分析概率的方法,属于非构造性证明,感性趣看这个视频

渐进分析(Asymptotics)

GV Bound 给出了 RR 的下界,结合 Hamming Bound 得到

1logq(Volq(d1,n))1n R1logq(Volq(d12,n))n1- \frac{log_q(Vol_q(d-1, n))-1}{n} \leq \ R \leq 1-\frac{\log_q(Vol_q(\lfloor\frac{d-1}{2}\rfloor,n))}{n}

这个不等式中,信息率 RR 与参数 d,n,qd,n,q 均有关,且形式复杂,不利于理解 RR 的性质。为此,我们使用渐进分析的方式,通过取极限去除某些参数,以研究 RR 最本质的性质。

信息率与距离

  1. C={Ci}i=1\mathfrak{C}=\{C_i\}_{i=1}^\infin 为编码构成的序列,其中 CiC_i(ni,ki,di)qi(n_i,k_i,d_i)_{q_i} 分组码,定义序列 C\mathfrak{C} 的信息率

    R(C):=limikiniR(\mathfrak{C}) :=\lim_{i\rightarrow\infty}\frac{k_i}{n_i}

    定义 C\mathfrak{C}相对距离(简称距离) δ(C)\delta(C)

    δ(C):=limidini\delta(\mathfrak{C}) :=\lim_{i\rightarrow\infty}\frac{d_i}{n_i}

  2. 前边讨论汉明界时,举了 (7,4,3)2(7,4,3)_2-汉明码的例子,类似地,容易构造编码序列 C={Ci}i=1\mathfrak{C}=\{C_i\}_{i=1}^\infty,其中 CiC_i(2i1,2ii1,3)2(2^i-1, 2^i-i-1, 3)_2 汉明码,此时

    R(C)=limi2ii12i1=1δ(C)=limi3ni=0\begin{align*} R(\mathfrak{C}) &= \lim_{i\rightarrow\infty}\frac{2^i-i-1}{2^i-1}=1\\ \delta(\mathfrak{C}) &=\lim_{i\rightarrow\infty}\frac{3}{n_i}=0 \end{align*}

    CC 的信息率趋近于 11,但距离趋近于 0,即纠错能力趋于没有。

  3. 我们希望在 δ(C)\delta(\mathfrak{C})R(C)R(\mathfrak{C}) 中找到最佳的协调方式,特别地,将满足 δ(C)0R(C)\delta(\mathfrak{C})\neq 0\neq R(\mathfrak{C}) 的编码序列称为渐进好的(Asymptotically Good)

边界分析

  1. GV Bound 和 Hamming Bound 公式中的体积球 VolVol 不利于讨论,为此先定义 q-ary 熵函数(q-ary entropy function),用于化简公式

    Hq:[0,1][0,1]xxlogq(q1)xlogq(x)(1x)logq(1x)\begin{align*} H_q:[0,1]&\rightarrow[0,1]\\ x&\mapsto xlog_q(q-1)-xlog_q(x)-(1-x)log_q(1-x) \end{align*}

    • x0x\rightarrow 0 时,上式趋 0
    • x1x\rightarrow 1 时,上式趋 logq(q1)log_q(q-1)
    • qq\rightarrow \infty 时,Hq(x)xH_q(x)\rightarrow x
    • 图像整体开口朝下
  2. q2Z, 0p11qq\geq 2 \in\mathbb{Z},\ 0\leq p\leq 1-\frac{1}{q},则(略去证明)

    qnHq(p)o(n)Volq(pn,n)qnHq(p)\begin{align*} q^{n\cdot H_q(p)-o(n)}\leq Vol_q(pn, n)\leq q^{n\cdot H_q(p)} \end{align*}

  3. 根据 Hamming Bound,对任意编码序列 C\mathfrak{C},有

    R11nlogq(Volq(d12,n))n1Hq(δ(C)2)\begin{align*} R &\leq 1-\frac{1}{n}\log_q(Vol_q(\lfloor\frac{d-1}{2}\rfloor,n))\\ &\underset{n\rightarrow \infty}{\leq}1-H_q(\frac{\delta(\mathfrak{C})}2) \end{align*}

    根据 GV Bound,存在编码序列 C\mathfrak{C},使得

    R11n(logq(Volq(d1,n))1)11n(nHq(d1n)o(n)1)n1Hq(δ(C))\begin{align*} R&\geq 1- \frac{1}{n}(log_q(Vol_q(d-1, n))-1)\\ &\geq 1-\frac{1}{n}(n\cdot H_q(\frac{d-1}{n})-o(n)-1)\\ &\underset{n\rightarrow \infty}{\geq}1-H_q(\delta(\mathfrak{C}))\\ \end{align*}

  4. 综上,有

    1Hq(δ(C))R(C)1Hq(δ(C)2)1-H_q(\delta(\mathfrak{C}))\leq R(\mathfrak{C})\leq 1-H_q(\frac{\delta(\mathfrak{C})}2)

    我们将这两个边界分别称为编码序列 C\mathfrak{C}渐进 GV Bound渐进 Hamming Bound。当 q=2q=2 时,图像表现为
    20220504112629

  5. 延伸

    • Q1:是否存在编码序列,满足渐进 GV Bound,即落在黄色区域
      • q49q\geq 49 时,问题成立
      • q=2q=2 时,这还是个开问题
    • 能否构造编码序列,满足渐进 GV Bound,同样地,这个问题也仅在 qq 足够大时存在
    • 渐进 GV Bound 给出蓝线以下的存在性,继而说明了存在“渐近好的”编码序列
    • 关于边界分析,还有 Plotkin Bound 等,将图像进一步细化,这些都描述了比率 R(C)R(\mathfrak{C}) 与距离 δ(C)\delta(\mathfrak{C}) 之间的关系
      20220504121301
    • 特别注意,本节讨论对编码序列的渐进分析,从图像上看 Singleton Bound 始终在 Plotkin Bound 上方,但这是在 qq 固定下取极限得到的图像;落实到具体编码时,Singleton Bound 是可以取到的,比如 RS 码。

编码效率

概述

本节讨论编码效率问题,特别地,博客只讨论线性分组码的情形。

编码过程的这几个步骤涉及计算:

  • 信息编码,即 FqkFqn\mathbb{F}_q^k\rightarrow\mathbb{F}_q^n
  • 信息解码,即 Im(C)FqkIm(C)\rightarrow\mathbb{F}_q^k
  • 填补缺失,比如 (,0,0)(0,0,0)(*,0,0) \rightarrow (0,0,0)
  • 检测错误,比如 (1,1,0)(1,1,0)
  • 计算编码距离
  • 纠正错误,比如 (1,0,0)(0,0,0)(1,0,0)\rightarrow(0,0,0)

前四步存在高效算法(多项式时间),但计算编码距离纠错不存在处理一般情形的高效算法。

这种纠错上的困难性,既是坏处也是好处:

  • 坏处在于:我们不能用任意线性码来制定编码规则,否则纠错带来的时间开销过大,效率过低
  • 好处在于,这种困难性可以应用于非对称加密系统

坏处是可以避免的,比如构造特殊编码,使编码具备更多性质,允许更多操作,继而存在更高效的算法,应用上我们并不需要设计能处理一般编码的算法。

好处是这类数论难题通常导致加密算法的产生,这部分理论与二维码关联很低,放在公钥系统中单独介绍。

编码算法

下边介绍线性码中编码,解码,补缺,检错,求距离以及纠错这些步骤的一般计算方法。

下设 CC(n,k,d)q(n, k,d)_q 线性分组码,生成元矩阵为 GFn×kG\in \mathbb{F}^{n\times k},奇偶检验矩阵为 HFk×nH\in \mathbb{F}^{k\times n},设 x=(x1,,xk)tCx=(x_1,\cdots,x_k)^t\in C 为信息。

  1. 编码:用矩阵 GG 左乘信息 xx

    C:FqkFqnxGx=x\begin{align*} C:\mathbb{F}_q^k&\rightarrow\mathbb{F}_q^n\\ x&\mapsto Gx=x' \end{align*}

  2. 解码:由线性代数知识易得矩阵 GFk×nG'\in \mathbb{F}^{k\times n} 使得 GGG'Gkk 阶单位矩阵。左乘矩阵 GG' 即从 xx' 中解码信息 xx

  3. 检错:根据 ker(H)=Im(C)\ker(H)=Im(C),只需判断 HxHx' 是否为 0。

  4. 补缺:不妨设 xx' 缺失了前 e<de<d 个数据,将 GG 相应分块化,得

    Gx=(G0G1)x=(x1xe1xexn)=(xexn)\begin{align*} Gx =\begin{pmatrix} G_0 \\ G_1 \end{pmatrix} x=\begin{pmatrix} x_1'\\\vdots\\ x_{e-1}'\\ x_e'\\ \vdots\\ x_n' \end{pmatrix}=\begin{pmatrix} *\\\vdots\\ *\\ x_e'\\ \vdots\\ x_n' \end{pmatrix} \end{align*}

    反解得到原始数据 xx

    G1x=(xexn),x=G11(xexn)\begin{align*} G_1x=\begin{pmatrix} x_e'\\ \vdots\\ x_n' \end{pmatrix}, x=G_1^{-1}\begin{pmatrix} x_e'\\ \vdots\\ x_n' \end{pmatrix} \end{align*}

    e<de<d 时,左侧方程组必有解,而 RS 码中 G1G_1 还是可逆的,此时直接左成逆矩阵还原数据。

  5. 求距离:回顾介绍了编码距离 dd 的两个计算性质:

    • dd 为生成元矩阵 GG 列空间中非零向量的最小权
    • dd 为奇偶检验矩阵 HH 的线性相关列向量组的最小长度

    这两条性质给出了 dd 的上下界,因此随机选取列向量来计算 dd,但文献论证了计算一般线性码的编码距离不是 RP 问题(nondeterministic polynomial time),也即不存在针对一般线性码的高效算法

  6. 纠错:假设码字 x=Gx+zx'=Gx+z,其中 zz 为纠错范围内的错误扰动,左乘矩阵 HH

    b=Hx=H(Gx+z)=HGx+Hz=Hzb = Hx' = H(Gx+z) = HGx + Hz = Hz

    求解线性方程组 Hz=bHz=b 计算错误扰动 zz,由于 zz 处在纠错范围内,存在唯一向量 zz 使得

    Hz=b, wt(z)d12Hz=b,\ wt(z) \leq \lfloor\frac{d-1}{2}\rfloor

    线性方程组的解构成了一个仿射空间(线性空间+平移向量),问题等同于求该空间中权重最小的向量。然而这个问题的求解仍是困难的,哪怕对 RS 码,参阅文献:Maximum-Likelihood Decoding of Reed-Solomon Codes is NP-hard

注:表面上看,纠错只是解线性方程,但由于这是对有限域上讨论,且问题不是求出方程的一个解,而是计算最小解。由于解空间是随着维数以指数形式增长,暴力破解自然行不通。这种表明简单,但应用到有限域就变得复杂的现象在数论中经常出现,比如计算离散对数 logq(a)log_q(a)


最后,回到文章重点,RS 码及纠错算法,由于 RS 码取到了 Singleton Bound,编码距离计算简单,难点在纠错上。

Reed-Solomon 纠错

Reed-Solomon 码是一组纠错码,由 Irving S. Reed 和Gustave Solomon 在 1960 年引入。它有许多应用,比如 MiniDiscs、CD、DVD、蓝光光盘、二维码、DSL、WiMAX 等等。

RS 码具有以下特点:

  • RS 码满足 Singleton Bound
  • RS 码存在高效的解码算法
  • RS 码在现实生活中被广泛应用

背景知识

具备高效算法的编码通常具有丰富的结构性质,RS 码便是如此,其为一类特殊线性码,以范德蒙德矩阵(Vandermonde matrix) 为生成元矩阵。

我们先介绍理论基础,再给出 RS 码的定义和性质(这么安排并不好,有空调整一下)。

下设 F\mathbb{F} 为数域,α=(α0,,αn1)Fn1\alpha=(\alpha_0,\cdots,\alpha_{n-1})\in\mathbb{F}^{n-1},且 αi\alpha_i 互不相等

  1. 数域可以分两大类,有理数域 Q\mathbb{Q},实数域 R\mathbb{R} 和复数域 C\mathbb{C} 等属于特征 0 域;有限域 FpF_p 及其扩张,属于特征 p 域,线性代数理论对两类数域均适用。

  2. 称数域 F\mathbb{F}有限域,若其仅有有限个元素,并将包含 qq 个元素的有限域记为 Fq\mathbb{F}_q,有限域也被称为伽罗瓦域(Galois Filed)。由近世代数可知

    • 有限域均形如 Fq\mathbb{F}_q,其中 q=prq=p^rpp 为素数
    • 同阶的有限域相互同构
    • Fq=Fq{0}\mathbb{F}_q^*=\mathbb{F}_q-\{0\},其乘法结构作成 q1q-1 阶循环群
    • 特别地,rr 取 1 时,Fp\mathbb{F}_p 称为素域,其结构简单,同构于 Z/pZ\mathbb{Z}/p\mathbb{Z}
    • 一般地,Fq\mathbb{F}_q 同构于素域 Fp\mathbb{F}_p 关于 xr1x^r-1 的分裂域。当 r1r\neq 1时, F\mathbb{F} 结构复杂,而不是简单的模 qq 结构。
  3. 性质F\mathbb{F} 上的 nn 次多项式 f(x)f(x) 至多存在 nn 个根。若干说明:

    • nn 个根”指 F\mathbb{F} 上的根,通过带余除法和归纳法容易证明
    • 代数基本定理表明 nn 次多项式恰有 nn 个根,这里 nn 个根指代数闭包 Fˉ\bar{\mathbb{F}} 上的根,比如实数域上 x2+1=0x^2+1=0 没有根,但在代数闭域 C\mathbb{C} 上有两个根
    • 零多项式的次数一般假定为无穷
  4. 如下定义范德蒙德矩阵 V(α,k)Fn×kV(\alpha,k)\in\mathbb{F}^{n\times k}

    V(α,k)=(1α0α02α0k11α1α12α1k11αn1αn12αn1k1)V(\alpha,k)= \begin{pmatrix} 1 & \alpha_0 & \alpha_0^2 & \cdots & \alpha_0^{k-1}\\ 1 & \alpha_1 & \alpha_1^2 & \cdots & \alpha_1^{k-1}\\ \vdots & \vdots & \vdots & \ddots & \vdots\\ 1 & \alpha_{n-1} & \alpha_{n-1}^2 & \cdots & \alpha_{n-1}^{k-1} \end{pmatrix}

    knk\leq n 时,由范德蒙德行列式易证 V(α,k)V(\alpha,k) 的秩为 kk

    1α0α02α0k11α1α12α1k11αk1αk12αk1k1=i>j(αiαj)0\begin{align*} \begin{vmatrix} 1 & \alpha_0 & \alpha_0^2 & \cdots & \alpha_0^{k-1}\\ 1 & \alpha_1 & \alpha_1^2 & \cdots & \alpha_1^{k-1}\\ \vdots & \vdots & \vdots & \ddots & \vdots\\ 1 & \alpha_{k-1} & \alpha_{k-1}^2 & \cdots & \alpha_{k-1}^{k-1} \end{vmatrix}=\prod_{i>j}(\alpha_i-\alpha_j)\neq 0 \end{align*}

  5. 范德蒙德矩阵与多项式有紧密联系:左乘运算可转化为多项式运算

    (1α0α02α0k11α1α12α1k11αn1αn12αn1k1)(a0a1ak1)=(i=0k1aiα0ii=0k1aiα1ii=0k1aiαn1i)=(f(α0)f(α1)f(αn1))\begin{align*} \begin{pmatrix} 1 & \alpha_0 & \alpha_0^2 & \cdots & \alpha_0^{k-1}\\ 1 & \alpha_1 & \alpha_1^2 & \cdots & \alpha_1^{k-1}\\ \vdots & \vdots & \vdots & \ddots & \vdots\\ 1 & \alpha_{n-1} & \alpha_{n-1}^2 & \cdots & \alpha_{n-1}^{k-1} \end{pmatrix} \begin{pmatrix} a_0\\ a_1\\ \vdots\\ a_{k-1} \end{pmatrix}= \begin{pmatrix} \sum\limits_{i=0}^{k-1}a_i\alpha_0^i\\ \sum\limits_{i=0}^{k-1}a_i\alpha_1^i\\ \vdots\\ \sum\limits_{i=0}^{k-1}a_i\alpha_{n-1}^i \end{pmatrix}= \begin{pmatrix} f(\alpha_0)\\ f(\alpha_1) \\ \vdots \\ f(\alpha_{n-1}) \end{pmatrix} \end{align*}

    其中 f(x)=i=0k1aixif(x)=\sum\limits_{i=0}^{k-1}a_ix^i,从这一性质也能看出 V(α,k)V(\alpha,k) 为列满秩矩阵

    V(α,k)X=Of(αi)=0,  0inf=0(deg(f)<=n)X=(a0,,ak)t=O\begin{align*} V(\alpha,k)X=O&\Leftrightarrow f(\alpha_i)=0,\ \forall\ 0\leq i\leq n\\ &\Rightarrow f=0\quad (\deg(f)<=n)\\ &\Rightarrow X=(a_0,\cdots, a_k)^t=O \end{align*}

    特别地,当 k=nk=n 时,V(α,k)V(\alpha,k) 为可逆矩阵

  6. 推论:对任意 {yi}i=0k1Fqn\{y_i\}_{i=0}^{k-1}\subseteq\mathbb{F}_q^n,存在唯一次数不大于 k1k-1 的多项式 ff,使得 f(αi)=yif(\alpha_i)=y_i。具体地,ff 的各项系数{ai}i=0k1\{a_i\}_{i=0}^{k-1} 由以下线性方程组得到

    V(α,n)(a0a1ak1)=(y0y1yn1)V(\alpha,n) \begin{pmatrix} a_0\\ a_1\\ \vdots \\ a_{k-1} \end{pmatrix}= \begin{pmatrix} y_0\\ y_1\\ \vdots \\ y_{n-1} \end{pmatrix}

  7. 推论2:有限域上,任意映射 FqF\mathbb{F}_q\rightarrow\mathbb{F} 可写为次数不超过 q1q-1 的多项式。此外,利用性质 xq=xx^q=x,任意多项式可以转化为次数不超过 q1q-1 的多项式。

这个现象很有趣:有限域上的函数没有可积可微的概念,所有映射都是多项式!且函数集是有限集。

基本概念

Reed Solomon 码是满足 Singleton Bound 且具有高效算法的一类编码。(以下是 RS 码的 Original View)

下设 qnkq\geq n\geq kα=(α0,,αn1)Fqn\alpha=(\alpha_0,\cdots,\alpha_{n-1})\in\mathbb{F}_q^n 且其中元素互不相同

  1. 如下定义 Fqn\mathbb{F}_q^n 的子集 RSq(α,n,k)RS_q(\alpha,n,k),称为信息长 kk,块长 nn 的关于向量 α\alpha 的 Reed-Solomon 码

    RSq(α,n,k)={(f(α0),,f(αn1))f(x)Fq[x], deg(f)k1}RS_q(\alpha,n,k)=\{(f(\alpha_0), \cdots, f(\alpha_{n-1}))|f(x)\in\mathbb{F}_q[x],\ \deg(f)\leq k-1\}

  2. 写成线性变换形式

    RSq(α,n,k):FqkFq[x]Fqn(a0an1)f(x)=i=0k1aixi(f(α0)f(αn1))\begin{align*} RS_q(\alpha,n,k):\mathbb{F}_q^k&\hookrightarrow \mathbb{F}_q[x]&\rightarrow \mathbb{F}_q^n\qquad\qquad\\ \begin{pmatrix}a_0\\\vdots\\a_{n-1}\end{pmatrix} &\mapsto f(x)=\sum_{i=0}^{k-1}a_ix^i &\mapsto \begin{pmatrix}f(\alpha_0)\\\vdots\\ f(\alpha_{n-1})\end{pmatrix} \end{align*}

    根据上一节讨论,其生成元矩阵为 V(α,k)V(\alpha, k)

  3. RS 码满足 Singleton Bound,根据以下等价关系可证

    • ddRSq(α,k,n)RS_q(\alpha,k,n) 的编码距离
    • dd 为生成元矩阵 V(α,k)V(\alpha,k) 列空间的最小权重
    • dd 为 n - (V(α,k)V(\alpha,k) 列空间的最大零数)

    一方面,由于 deg(f)k1\deg(f)\leq k-1f(αi)f(\alpha_i) 中至多有 k1k-1 项为 0,因而 dn(k1)d\geq n-(k-1);另一方面,取 f(x)=(xα0)(xαk1)f(x)=(x-\alpha_0)\cdots(x-\alpha_{k-1}) 得到权重为 nk+1n-k+1 的列向量,证毕。

  4. 将满足 Singleton Bound 的编码称为 MDS(Maximum Distance Separable) 码 。容易证明:(n,k,d)q(n,k,d)_q 线性分组码为 MDSMDS 码当且仅当其生成元矩阵 GG 的任意 k×kk\times k 阶子矩阵可逆。

特殊 RS 码

我们考虑在编码中更常用一种的 RS 码:取 n=q1n=q-1

  1. 我们先给出一个引理,用于求解该 RS 码的奇偶检验矩阵

    aFqad=0,  0<d<q1\begin{align*} \sum_{a\in\mathbb{F}_q}a^d=0,\ \forall\ 0<d<q-1 \end{align*}

    证明:设 γ\gammaFq\mathbb{F}_q^* 乘法群的生成元

    aFqad=aFad=i=1q1γid=i=1q1(γd)i=1(γd)q1γd1=0\begin{align*} \sum_{a\in\mathbb{F}_q}a^d &=\sum_{a\in\mathbb{F}^*}a^d =\sum_{i=1}^{q-1}\gamma^{i\cdot d}\\ &=\sum_{i=1}^{q-1}(\gamma^{d})^i =\frac{1-(\gamma^d)^q}{1-\gamma^d}-1=0 \end{align*}

    根据循环群性质, aq11,a0a^{q-1}\equiv 1,\forall a\neq 0,继而有 (γd)qγd(\gamma^d)^q\equiv \gamma^d

  2. α=(γ1,γ2,γp1)\alpha=(\gamma^1,\gamma^2\cdots,\gamma^{p-1})RSq(α,k,n)RS_q(\alpha,k,n) 的生成元矩阵 V(α,k)V(\alpha,k)

    V(α,k)=(1γ1γ2γk11(γ2)1(γ2)2(γ2)k11(γp2)1(γp2)2(γp2)k11111)\begin{align*} V(\alpha,k)&=\begin{pmatrix} 1 & \gamma^1 & \gamma^2 & \cdots & \gamma^{k-1}\\ 1 & (\gamma^2)^1 & (\gamma^2)^2 & \cdots & (\gamma^2)^{k-1}\\ \vdots & \vdots & \vdots & \ddots & \vdots\\ 1 & (\gamma^{p-2})^1 & (\gamma^{p-2})^2 & \cdots & (\gamma^{p-2})^{k-1}\\ 1 & 1 & 1 & \cdots & 1 \end{pmatrix} \end{align*}

    其奇偶检验矩阵由下式确定

    RSq(α,n,k)={(c1,c2,,cq1)c(γj)=0,1jnk}where c(x)=l=1q1clxl=l=1q1f(γl)xl\begin{align*} RS_q(\alpha,n,k)&=\{(c_1,c_2,\cdots,c_{q-1})|c(\gamma^j)=0,\forall 1\leq j\leq n-k\}\\ where\ c(x)&=\sum\limits_{l=1}^{q-1}c_lx^l =\sum\limits_{l=1}^{q-1}f(\gamma^l)x^l\\ \end{align*}

    利用前一引理易证

    c(γj)=l=1p1f(γl)(γj)l=l=1p1(i=0k1ai(γl)i)(γj)l=i=0k1ai(l=1p1(γl)(i+j))=0\begin{align*} c(\gamma^j) &= \sum_{l=1}^{p-1}f(\gamma^l)\cdot(\gamma^j)^l\\ &=\sum_{l=1}^{p-1}\left(\sum_{i=0}^{k-1}a_i(\gamma^l)^i\right)\cdot(\gamma^j)^l\\ &=\sum_{i=0}^{k-1}a_i\left(\sum_{l=1}^{p-1}(\gamma^l)^{(i+j)}\right) =0 \end{align*}

    最后一式,1i+jnk+in1q11\leq i+j\leq n-k+i\leq n-1\leq q-1
    c(γj)=0c(\gamma^j)=0 的必要性通过比较维数立证。

  3. 推论:RSq((γ1,,γq1),n,k)RS_q((\gamma^1,\cdots,\gamma^{q-1}),n,k) 的奇偶检验矩阵为

    H=V((γ1,,γq1),nk)t=(1γγ2γn11(γ2)1(γ2)2(γ2)n11(γnk)1(γnk)2(γnk)n1)\begin{align*} H &=V((\gamma^1,\cdots,\gamma^{q-1}),n-k)^t\\ &=\begin{pmatrix} 1&\gamma&\gamma^2&\cdots&\gamma^{n-1}\\ 1&(\gamma^2)^1&(\gamma^2)^2&\cdots&(\gamma^2)^{n-1}\\ \vdots&\vdots&\vdots&\ddots&\vdots\\ 1&(\gamma^{n-k})^1&(\gamma^{n-k})^2&\cdots&(\gamma^{n-k})^{n-1}\\ \end{pmatrix} \end{align*}

  4. 推广的RS码(Generalized RS code) 定义如下

    GRSq(α,n,k;λ)={(λ0f(α0),λ1f(α1),,λn1f(αn1)),fFq[x],deg(f)k1}GRS_q(\alpha,n,k;\overrightarrow\lambda)= \{(\lambda_0f(\alpha_0),\lambda_1f(\alpha_1),\cdots,\lambda_{n-1}f(\alpha_{n-1})),f\in\mathbb{F}_q[x],\deg(f)\leq k-1\}

    其中 α(Fq)n\alpha\in(\mathbb{F}_q^*)^n ,此外生成元矩阵为

    G=MV(α,k), diag(M)=(λ0,,λn1)\begin{align*} G &=M\cdot V(\alpha,k),\ diag(M)=(\lambda_0,\cdots,\lambda_{n-1}) \end{align*}

    一般地,可以推导 GRSq(α,n,k;λ)GRS_q(\alpha,n,k;\overrightarrow\lambda) 的对偶码

简单说,RS 码为多项式 fff.αf.\alpha 的映射;另一方面,当编码维数为 q1q-1 时,也可以将 RS 码视为多项式到多项式的映射,即 f(x)c(x)f(x) \rightarrow c(x)

纠错算法

本节介绍 Reed Solomon 的 Berlekamp–Welch 算法。

  1. 我们先表述纠错问题:

    • RSq(α,n,k)RS_q(\alpha,n,k) 为 RS 码,编码距离 d=nk+1d=n-k+1,纠错能力 ed12=nk2e\leq\lfloor\frac{d-1}{2}\rfloor=\lfloor\frac{n-k}{2}\rfloor
    • w=(w0,,wn1)Fqn\overrightarrow{w}=(w_0,\cdots,w_{n-1})\in\mathbb{F}_q^n 为接收到的编码,且至多存在 ee 个错误
    • 求多项式 f(x)Fq[x]f(x)\in\mathbb{F}_q[x] 使得:deg(f)k1\deg(f)\leq k-1 且至多存在 ee 个数值 αi\alpha_i,使得 f(αi)wif(\alpha_i)\neq w_i
  2. 我们表述算法的过程,再论证细节:

    • 计算次数为 ee 的首一多项式 E(x)E(x) 以及次数 e+k1e+k-1 的多项式 Q(x)Q(x) 使得

    wiE(αi)=Q(αi),  i=0,,n1w_iE(\alpha_i) = Q(\alpha_i),\ \forall\ i=0,\cdots,n-1

    这里 E(x)E(x) 称为错误检查多项式(Error locator polynomial)

    • 若找不到满足这一性质的多项式 E(x)E(x)Q(x)Q(x),则说明 w\overrightarrow w 超出了纠错范围
    • f~(x)=Q(x)/E(x)\tilde{f}(x)=Q(x)/E(x),若 Δ(f~,w)e\Delta(\tilde{f},w)\leq e,则 ff 为所求,否则说明 w\overrightarrow w 超出了纠错范围。这里 Δ(f~,w)\Delta(\tilde{f},w) 为满足 f~(αi)wi\tilde{f}(\alpha_i)\neq w_iii 的个数
  3. 我们先证明,若 ww 在纠错范围内,则满足性质的多项式 E(x)E(x)Q(x)Q(x)。设 ff 为原始信息,如下定义 E(x)E(x)Q(x)Q(x)

    E(x)=i:f(αi)wi(xαi)xeΔ(f,w)Q(x)=E(x)f(x)\begin{align*} E(x) &= \prod_{i:f(\alpha_i)\neq w_i}(x-\alpha_i)\cdot x^{e-\Delta(f,w)}\\ Q(x) &= E(x)\cdot f(x) \end{align*}

    f(αi)wif(\alpha_i)\neq w_i 时,

    wiE(αi)=Q(αi)=0w_iE(\alpha_i)=Q(\alpha_i)=0

    f(αi)=wif(\alpha_i)=w_i 时,亦有

    Q(αi)=E(αi)f(αi)=wiE(αi)Q(\alpha_i)=E(\alpha_i)f(\alpha_i)=w_iE(\alpha_i)

  4. 另一方面,若存在另一组 Q(x),E(x)Q'(x),E'(x) 满足同一性质,我们断言

    Q(x)/E(x)=Q(x)/E(x)Q'(x)/E'(x)=Q(x)/E(x)

    换言之,即证

    R(x)=Q(x)E(x)Q(x)E(x)=0R(x)=Q(x)\cdot E'(x)-Q'(x)\cdot E(x)=0

    注意到 R(x)R(x) 的次数为

    deg(R(x)=e+e+k1nk+k1<n\deg(R(x)= e+e+k-1\leq n-k+k-1<n

    然而 R(x)R(x) 存在 nn 个根,继而 R(x)=0R(x)=0

    R(αi)=Q(αi)E(αi)Q(αi)E(αi)=wiE(αi)E(αi)wiE(αi)E(αi)=0, i=0,,n1\begin{align*} R(\alpha_i)&=Q(\alpha_i)\cdot E'(\alpha_i)-Q'(\alpha_i)\cdot E(\alpha_i)\\ &=w_iE(\alpha_i)E'(\alpha_i)-w_iE'(\alpha_i)E(\alpha_i)\\ &=0,\quad \forall\ i=0,\cdots, n-1 \end{align*}

  5. 结合两部分讨论,若 ww 在纠错范围内,满足这一性质的 E(x)E(x)Q(x)Q(x) 必然存在,且计算 Q(x)/E(x)Q(x)/E(x) 立得原始信息 f(x)f(x)