一种基于有限域的二进制和多进制空间耦合LDPC码构造方法
An Approach for Constructing Binary and Nonbinary Spatially-Coupled LDPC Codes Based on Finite Fields
DOI: 10.12677/AAM.2021.102063, PDF, HTML, XML, 下载: 338  浏览: 477 
作者: 梁 宇, 杨卫华*, 李玉瑛:太原理工大学数学学院,山西 晋中
关键词: 空间耦合LDPC码准循环LDPC码有限域围长Spatially-Coupled LDPC Codes Quasi-Cyclic LDPC Codes Finite Field Girth
摘要: 空间耦合LDPC (Spatially-Coupled LDPC, SC-LDPC)码由于阈值饱和特性,被证明是未来无线通信系统的有力候选码型。SC-LDPC码是一种卷积LDPC码,在二元无记忆对称信道下采用置信传播算法时具有逼近香农限的性能。本文提出了一种构造二进制和多进制空间耦合LDPC码的统一方法。基于有限域GF(q)上的本原元、加法子群、乘法子群,构造了6类二进制和多进制空间耦合准循环LDPC码(Spatially-Coupled Quasi­Cyclic LDPC, SC-QC-LDPC),用此构造的LDPC码的围长至少为6。
Abstract: Spatially-coupled LDPC codes are demonstrated to strong candidates for future optical communications systems due to the threshold saturating property. Spatially-coupled LDPC code is a kind of convolutional LDPC code. It has the performance of approaching Shannon limit when using belief propagation decoding algorithm in binary memoryless symmetric channel. A unified approach for constructing binary and nonbinary spatially-coupled LDPC codes is presented in this paper. Six classes of binary and nonbinary spatially-coupled Quasi-Cyclic LDPC codes are constructed based on primitive elements, additive subgroups, multiplicative subgroups of finite fields. Moreover, the girth of these codes is greater than or equal to 6.
文章引用:梁宇, 杨卫华, 李玉瑛. 一种基于有限域的二进制和多进制空间耦合LDPC码构造方法[J]. 应用数学进展, 2021, 10(2): 575-586. https://doi.org/10.12677/AAM.2021.102063

1. 引言

空间耦合LDPC码源自于LDPC卷积码 [1],由于其优良的阈值而受到人们的广泛关注。对于二元无记忆对称(Binary memoryless symmetric, BMS)信道,已经证明了规则LDPC码集 [2] 的最大后验概率(Maximum aposteriori probability, MAP)阈值可以通过空间耦合原始LDPC码集的生成码集的置信传播(Belief propagation, BP)阈值来逼近 [3] [4] [5] [6],这就是所谓的阈值饱和现象。除了渐近性能分析外,SC­LDPC码的有限长性能特别是利用窗口译码方法进行译码,也成为研究的重点 [7] [8]。最近,多进制SC­LDPC码也开始受到研究人员的关注。多进制SC­LDPC码的阈值饱和现象在二元删除信道(Binary erasure channel, BEC)被研究 [9]。经典的SC­LDPC码构造,被分类为图覆盖结构,是基于展开LDPC码的校验矩阵来形成SC­LDPC码的校验矩阵 [1] [10]。

准循环LDPC码 [11] 由于其低复杂度和高度并行的编码 [12] 和译码 [13],已被各种通信系统标准化,最近又应用于数据存储产品中。QC­LDPC码通常是由在有限域上相同大小的稀疏循环矩阵的阵列的零空间给出的 [15]。QC­LDPC码具有优良有限长性能和高效硬件实现的优点 [12] [13] [14]。广泛的仿真结果 [15] [16] [17] 已经表明设计良好的QC­LDPC码优于非结构化随机LDPC码。仿真结果还表明,非二进制原模图LDPC码的性能优于与其对应的二进制码 [18]。

由于QC­LDPC码和SC­LDPC码的优点,本文从理论上研究了在任意有限域上具有准循环结构的SC­LDPC码的构造,此类码称之为SC­QC­LDPC码。最近,二进制SC­QC­LDPC码的构造 [19] 使用了一种改进的渐进边缘增长(Progressive edge growth, PEG)算法 [20]。相反,本文构造SC­QC­LDPC码是基于有限域的确定性构造。这避免了耗时的计算机搜索,同时确保了校验矩阵的Tanner图的围长至少为6,这也通常保证了码的优良性能。

2. 预备知识

2.1. 有限域中元素的位置向量

考虑有限域 G F ( q ) ,q为一个素数的幂。令 α G F ( q ) 的本原元,则 α 的q个幂次 α = 0 α 0 = 1 α , , α q 2 表示出了 G F ( q ) 的所有元素,且 α q 1 = 1 G F ( q ) 包含两个群:加法群和乘法群。 G F ( q ) 的q个元素在加法运算下构成加法群, G F ( q ) q 1 个非零元素在乘法运算下构成乘法群。

对于任意一个非零元素 α i 0 i < q 1 ,我们定义一个 G F ( 2 ) 上的 q 1 维向量: z b ( α i ) = ( z 0 , z 1 , , z q 2 ) ,其中第i个分量 z i = 1 ,而其他 q 2 个分量等于零。将这样只有一个1分量的 G F ( 2 ) 上的 q 1 维向量称为元素 α i 的二进制位置向量(binary location­vector)。显然,域 G F ( q ) 中两个不同非零元素的二进制位置向量的1分量在不同位置上。元素0的二进制位置向量定义为 q 1 维全零向量 z b ( 0 ) = ( 0 , 0 , , 0 ) 。类似地,对于任意一个非零元素 α i 0 i < q 1 ,我们定义一个 G F ( q ) 上的 q 1 维向量 z q ( α i ) = ( z 0 , z 1 , , z q 2 ) ,其中第i个分量 z i = α i ,而其他 q 2 个分量等于零。将这样只有一个 α i 分量的 G F ( q ) 上的 q 1 维向量称为元素 α i 的q进制位置向量(q­ary location­vector)。元素0的q进制位置向量定义为 q 1 维全零向量 z q ( 0 ) = ( 0 , 0 , , 0 )

2.2. 有限域中元素的矩阵散列

δ G F ( q ) 的一个非零元素。当 0 i < q 1 时, α i δ 的二进制位置向量 z b ( α i δ ) α i 1 δ 的二进制位置向量 z b ( α i 1 δ ) 的循环移位(右移一位)。因为 α q 1 = 1 ,则 z b ( α q 1 δ ) = z b ( δ ) 。因而, z b ( δ ) 就是 z b ( α q 2 δ ) 的循环移位。以元素 δ , α δ , , α q 2 δ 的二进制位置向量作为行,可以得到一个 G F ( 2 ) 上的 ( q 1 ) × ( q 1 ) 矩阵 A b ( δ ) 。该矩阵是一个循环置换矩阵(Circulant permutation matrix, CPM),其每一行都是上一行的循环移位,而第一行是最后一行的循环移位。将矩阵 A b ( δ ) 称为元素 δ G F ( 2 ) 上的 ( q 1 ) 重二进制矩阵散列((q-1)­fold binary matrix dispersion)。类似地,我们以元素 δ , α δ , , α q 2 δ 的q进制位置向量作为行,可以得到一个 G F ( q ) 上的 ( q 1 ) × ( q 1 ) 矩阵 A q ( δ ) 。该矩阵是一个q进制 α 倍循环置换矩阵(q­ary α­multiplied circulant permutation matrix),其每一行都是上一行乘以 α 的循环移位,而第一行是最后一行乘以 α 的循环移位。将矩阵 A q ( δ ) 称为元素 δ G F ( q ) 上的 ( q 1 ) 重q进制矩阵散列((q-1)­fold q­ary matrix dispersion)。

3. 基于有限域构造QC­LDPC码的一般方法

α G F ( q ) 的本原元。构造过程始于一个 G F ( q ) 上的 m × n 矩阵

B = [ b 0 b 1 b m 1 ] = [ b 0 , 0 b 0 , 1 b 0 , n 1 b 1 , 0 b 1 , 1 b 1 , n 1 b m 1 , 0 b m 1 , 1 b m 1 , n 1 ]

其行满足如下两个约束条件:1) 若 0 i < m 0 k , l < q 1 ,且 k l ,则 α k b i α l b i 至多有一个对应位置上的值相同(即它们至少有 n 1 个对应位置上的值不同);2) 若 0 i , j < m i j ,且 0 k , l < q 1 ,则 α k b i α l b j 至少有 n 1 个对应位置上的值不同。矩阵B的两个约束条件分别称为 α 乘行约束1和2 (α­multiplied row constraint 1 and 2)。

定理1:如果矩阵B的每一行至多有一个0元素且矩阵B的任意一个 2 × 2 子矩阵是非奇异的,则矩阵B满足 α 乘行约束1和2。

证明:因为当且仅当 b i , s = 0 时, α k b i , s = α l b i , s ,其中 0 i < m 0 k , l < q 1 ,且 k l ,所以当且仅当矩阵B的每一行至多有一个0元素,矩阵B满足 α 乘行约束1。假设矩阵B不满足 α 乘行约束2,则 α k b i , s = α l b j , s α k b i , t = α l b j , t ,其中 0 i , j < m i j 0 k , l < q 1 0 s , t < n s t 。因此 α k b i , s α l b j , t = α l b j , s α k b i , t ,由此推得 b i , s b j , t b i , t b j , s = 0 即矩阵B存在 2 × 2 奇异子矩阵,矛盾。

定理2:如果矩阵B满足 α 乘行约束1和2,则矩阵B满足 2 × 2 子矩阵约束,即矩阵B中任意 2 × 2 子矩阵至少有一个0元素或是非奇异的。

证明:如果B满足 α 乘行约束1和2,则 α k b i , s = α l b i , s α k b i , s α l b i , s α k b i , s = α l b j , s α k b i , s α l b j , s ,其中 0 i , j < m i j 0 k , l < q 1 0 s , t < n s t

情况1:当 α k b i , s = α l b i , s α k b i , s = α l b j , s 时,因为 α k b i , s = α l b i , s ,所以 b i , s = 0 ;又因为 α k b i , s = α l b j , s ,所以 b j , s = b i , s = 0 。因此 B 2 × 2 = [ b i , s b i , t b j , s b j , t ] = [ 0 b i , t 0 b j , t ] ,故B满足 2 × 2 子矩阵约束。

情况2:当 α k b i , s = α l b i , s α k b i , s α l b j , s 时,因为 α k b i , s = α l b i , s ,所以 b i , s = 0 ;又因为 α k b i , s α l b j , s ,所以 b j , s 0 。因此 B 2 × 2 = [ b i , s b i , t b j , s b j , t ] = [ 0 b i , t b j , s b j , t ] det ( B 2 × 2 ) = b i , t b j , s 0 ,故B满足 2 × 2 子矩阵约束。

情况3:当 α k b i , s α l b i , s α k b i , s = α l b j , s 时,因为 α k b i , s = α l b i , s α k b i , s = α l b j , s ,所以 α k b i , s α l b j , t α k b i , t α l b j , s ,即 b i , s b j , t b i , t b j , s 。因此 B 2 × 2 = [ b i , s b i , t b j , s b j , t ] det ( B 2 × 2 ) = b i , s b j , t b i , t b j , s 0 ,故B满足 2 × 2 子矩阵约束。

情况4:当 α k b i , s α l b i , s α k b i , s α l b j , s 时,因为 α k b i , s α l b i , s α k b i , s α l b j , s ,所以 α k b i , s α l b j , t α k b i , t α l b j , s ,即 b i , s b j , t b i , t b j , s 。因此 B 2 × 2 = [ b i , s b i , t b j , s b j , t ] det ( B 2 × 2 ) = b i , s b j , t b i , t b j , s 0 ,故B满足 2 × 2 子矩阵约束。

因此如果矩阵B满足 α 乘行约束1和2,则矩阵B满足 2 × 2 子矩阵约束。

将矩阵B的每一项 b i , j 用其相应的 ( q 1 ) 重二进制矩阵散列 A b ( b i , j ) 替换后,可以得到一个 m × n 阵列 H b

H b = [ A b ( b 0 , 0 ) A b ( b 0 , 1 ) A b ( b 0 , n 1 ) A b ( b 1 , 0 ) A b ( b 1 , 1 ) A b ( b 1 , n 1 ) A b ( b m 1 , 0 ) A b ( b m 1 , 1 ) A b ( b m 1 , n 1 ) ]

其中,若 b i , j 0 ,则 A b ( b i , j ) G F ( 2 ) 上的一个 ( q 1 ) × ( q 1 ) CPM;若 b i , j = 0 ,则 A b ( b i , j ) 是一个全零矩阵。 H b G F ( 2 ) 上的一个 m ( q 1 ) × n ( q 1 ) 矩阵。

类似地,将矩阵B的每一项 b i , j 用其相应的 ( q 1 ) 重q进制矩阵散列 A q ( b i , j ) 替换后,可以得到一个 m × n 阵列 H q

H q = [ A q ( b 0 , 0 ) A q ( b 0 , 1 ) A q ( b 0 , n 1 ) A q ( b 1 , 0 ) A q ( b 1 , 1 ) A q ( b 1 , n 1 ) A q ( b m 1 , 0 ) A q ( b m 1 , 1 ) A q ( b m 1 , n 1 ) ]

其中,若 b i , j 0 ,则 A q ( b i , j ) G F ( q ) 上的一个 α 倍循环置换矩阵;若 b i , j = 0 ,则 A b ( b i , j ) 是一个全零矩阵。

H b H q 分别是 G F ( 2 ) G F ( q ) 上的一个 m ( q 1 ) × n ( q 1 ) 矩阵,将其分别称矩阵B的 ( q 1 ) 重二进制和q进制阵列散列((q-1)­fold binary and q­ary array dispersions)。

对于任意整数对 γ ρ 1 γ m 1 ρ n ,令 H b ( γ , ρ ) H q ( γ , ρ ) 分别表示 H b H q 的一个 γ × ρ 子阵列,则 H b ( γ , ρ ) H q ( γ , ρ ) 分别是 G F ( 2 ) G F ( q ) 上的一个 γ ( q 1 ) × ρ ( q 1 ) 矩阵。 H b ( γ , ρ ) G F ( 2 ) 上的零空间给出了一个二进制QC­LDPC码 C b ,类似地, H q ( γ , ρ ) G F ( q ) 上的零空间给出了一个q进制QC­LDPC码 C q ,其码长均为 ( q 1 ) ρ 。若 H b ( γ , ρ ) H q ( γ , ρ ) 均不包含全零矩阵,则其具有固定的行重 γ 和列重 ρ ,码 C b C q 均是规则的,反之,它们均为非规则的。

定理3:如果矩阵B满足 2 × 2 子矩阵约束,则 H b H q 满足RC­约束(Row­column constraint),即 H b H q 中任意两行(或两列)至多有一个对应位置都为非零项。

证明:假设 H b 不满足RC­约束。因为 H b 是一个循环置换矩阵的阵列,所以存在 i , j , s , t , k , l ,其中 0 i , j < m i j 0 s , t < n s t 0 k , l < q 1 ,使得 A b ( b i , s ) 的第k行与 A b ( b j , s ) 的第l行相等,即 z b ( α k b i , s ) = z b ( α l b j , s ) ;且 A b ( b i , t ) 的第k行与 A b ( b j , t ) 的第l行相等,即 z b ( α k b i , t ) = z b ( α l b j , t ) 。因为在 G F ( q ) 上不同的元素有不同的位置向量,所以 α k b i , s = α l b j , s α k b i , t = α l b j , t 。因此 α k b i , s α l b j , t = α k b i , t α l b j , s ,即 b i , s b j , t = b i , t b j , s ,与矩阵B满足 2 × 2 子矩阵约束矛盾,故 H b 满足RC­约束。同理,可证 H q 也满足RC­约束。

4. 基于有限域构造SC­QC­LDPC码的一般方法

考虑 G F ( q ) 上一个满足 2 × 2 子矩阵约束 m × n 矩阵 B = [ b i , j ] 0 i < m , 0 j < n ,复制矩阵B形成B的一个二维半无限阵列 B r e p

B r e p = [ b r e p , i , j ] 0 i , j < = [ B B B B B B B B B ] ,

其中对于任意 i , j 0 b r e p , i , j = b ( i mod m ) , ( j mod n ) 。其次,构造一个具有带对角线结构的 s × t 修饰矩阵 W b a s e = [ w b a s e , i , j ] 0 i < s , 0 j < t ,使用 W b a s e 修饰 B r e p ( s , t ) 形成一个SC基矩阵 B s c = B r e p ( s , t ) W b a s e 。最后,利用 B s c 的二进制阵列散列形成一个 G F ( 2 ) 上的 s ( q 1 ) × t ( q 1 ) 矩阵 H s c , q c , b ,其零空间给出一个二进制SC­QC­LDPC码 C s c , q c , b ;类似地,利用 B s c 的q进制阵列散列形成一个 G F ( q ) 上的 s ( q 1 ) × t ( q 1 ) 矩阵 H s c , q c , q ,其零空间给出一个q进制SC­QC­LDPC码 C s c , q c , q

在原模图构造中,我们可以把修饰矩阵 W b a s e 看作是与SC原模图相对应的原模矩阵(每对变量节点(Variable node, VN)和校验节点(Check node, CN)之间的边数不超过1)。此外,上述结构等价于提升 W b a s e 相对应的SC原模图,其提升因子为 q 1 。然而,提升(边的置换)是使用代数方法进行的。更一般地,这种构造也可以看作是SC­LDPC码的叠加构造 [21],叠加基矩阵与修饰矩阵 W b a s e 完全相同。

对于给定的修饰矩阵 W b a s e ,我们可以得出最终的SC­QC­LDPC码的校验矩阵 H s c , q c 完全取决于B的选择。下面的定理是关于B的选择如何影响最终的SC­QC­LDPC码的围长性质。

定理4:假设B是 G F ( q ) 上一个满足 2 × 2 子矩阵约束的 m × n 矩阵,且 B s c = B r e p ( s , t ) W b a s e 。对于任意的 i 1 , i 2 , j 1 , j 2 使得 w b a s e , i 1 , j 1 = w b a s e , i 1 , j 2 = w b a s e , i 2 , j 1 = w b a s e , i 2 , j 2 = 1 ,若 i 1 i 2 j 1 j 2 分别不能被m、n整除,则 B s c 满足 2 × 2 子矩阵约束。

证明:假设 B s c 不满足 2 × 2 子矩阵约束,则存在 i 1 , i 2 , j 1 , j 2 ,使得 B s c , i 1 , j 1 B s c , i 2 , j 2 = B s c , i 2 , j 1 B s c , i 1 , j 2 ,且 w b a s e , i 1 , j 1 = w b a s e , i 1 , j 2 = w b a s e , i 2 , j 1 = w b a s e , i 2 , j 2 = 1 。再根据B和 B r e p ( s , t ) 的关系,可得出 i 1 i 2 j 1 j 2 分别能被m、n整除,矛盾。因此, B s c 满足 2 × 2 子矩阵约束。

上述二进制和多进制SC­QC­LDPC码的结构是基于 G F ( q ) 上基矩阵B的 ( q 1 ) 重二进制矩阵散列或q进制矩阵散列和修饰矩阵 W b a s e 构造,其中基矩阵B满足 2 × 2 子矩阵约束。现在的问题是如何构造基矩阵B和修饰矩阵 W b a s e 。接下来的三节将研究上述问题。

5. 基于有限域本原元的SC­QC­LDPC码构造

基矩阵B的构造:

对于 0 < u < q 1 ,当且仅当u和 q 1 互素, G F ( q ) 上的非零元素 α u 是一个本原元。 G F ( q ) 中本原元的个数K可由下式给出的欧拉函数确定,其中 p 1 , p 2 , , p k q 1 的不同素因子。

K = ( q 1 ) i = 1 k ( 1 1 p i )

{ α u 1 , α u 2 , , α u K } G F ( q ) 中K个本原元组成的集合,令 u 0 = 0 。构造 G F ( q ) 上的 ( K + 1 ) × ( K + 1 ) 矩阵:

B ( 1 ) = [ b 0 b 1 b m 1 ] = [ α u 0 u 0 1 α u 1 u 0 1 α u K u 0 1 α u 0 u 1 1 α u 1 u 1 1 α u K u 1 1 α u 0 u K 1 α u 1 u K 1 α u K u K 1 ] (5.1)

由上式,可以看出,矩阵 B ( 1 ) 具有以下结构性质:1) 每一行(列)的所有元素在 G F ( q ) 上是不同的元素;2) 每一行(列)包含一个有且仅有一个零元素;3) 任何两行(列)对应位置上的元素均不同;4) K + 1 个0元素均位于矩阵的主对角线上。

推论1:矩阵 B ( 1 ) 满足 2 × 2 子矩阵约束。

证明:基于定理1和定理2,对于矩阵 B ( 1 ) 的任意 2 × 2 子矩阵 B 2 × 2 ( 1 ) ,其元素属于第 i 1 行、第 i 2 行和第 j 1 列、第 j 2 列,其中 0 i 1 , i 2 , j 1 , j 2 K + 1 i 1 i 2 j 1 j 2 ,其行列式 det ( B 2 × 2 ( 1 ) ) = ( α u j 1 u i 1 1 ) ( α u j 2 u i 2 1 ) ( α u j 2 u i 1 1 ) ( α u j 1 u i 2 1 ) = ( α u i 1 α u i 2 ) ( α u j 1 α u j 2 ) α u i 1 + u i 2 0 ,因此,矩阵 B ( 1 ) 满足 2 × 2 子矩阵约束。

因此, B ( 1 ) 可作为进行 ( q 1 ) 重阵列散列的基矩阵来构造二进制和多进制SC­QC­LDPC码。

修饰矩阵 W b a s e 的构造:

基于定理4,选取一个具有SC结构原模图的原模矩阵 [22] 作为修饰矩阵来修饰基矩阵B。

B = [ 3 3 ] B 0 = B 1 = B 2 = [ 1 1 ]

下面以几个例子来说明基于有限域本原元的SC­QC­LDPC码的构造。

5.1. 一类二进制SC­QC­LDPC码

例5.1.

基矩阵B的选取:将 G F ( 2 8 ) 作为构造码的域。注意到 2 8 1 = 255 可分解为 3 × 5 × 17 。根据欧拉方程可知 G F ( 2 8 ) K = 255 × ( 1 1 / 3 ) ( 1 1 / 5 ) ( 1 1 / 17 ) = 128 个本原元。利用式(5.1),我们可以构造一个128 × 128基矩阵B。

修饰矩阵 W b a s e 的选取:选择 L = 40 a = 3 b = 2 c = 1 ,且选择40个分量规则矩阵 { Δ k } ( 0 k 4 ) 为全一矩阵。

二进制SC­QC­LDPC码的构造:利用算法1给出的构造方法,可形成一个 G F ( 2 ) 上的5334 × 10,160矩阵 H s c , q c , b ,其零空间给出一个二进制SC­QC­LDPC码 C s c , q c , b

5.2. 一类多进制SC­QC­LDPC码

例5.2.

基矩阵B的选取:将 G F ( 2 6 ) 作为构造码的域。注意到 2 6 1 = 63 可分解为 3 × 3 × 7 。根据欧拉方程可知 G F ( 2 8 ) K = 63 × ( 1 1 / 3 ) ( 1 1 / 7 ) = 36 个本原元。利用式(5.1),我们可以构造一个36 × 36基矩阵B。

修饰矩阵 W b a s e 的选取:选择 L = 40 a = 3 b = 2 c = 1 ,且选择40个分量规则矩阵 { Δ k } ( 0 k 4 ) 为全一矩阵。

多进制SC­QC­LDPC码的构造:利用算法2给出的构造方法,可形成一个 G F ( q ) 上2646 × 5040矩阵 H s c , q c , q ,其零空间给出一个q进制SC­QC­LDPC码 C s c , q c , q

6. 基于有限域的加法子群构造SC­QC­LDPC码

基矩阵B的构造:令 q = p m ,其中p为素数且m为正整数。设 G F ( q ) 是素域 G F ( q ) 扩域, α G F ( q ) 的本原元,则 α 0 , α 1 , , α m 1 G F ( q ) 上线性无关,它们构成了 G F ( q ) 的一组基。 G F ( q ) 的任意元素 α i 都可以表示为 α 0 , α 1 , , α m 1 的线性组合,即 α i = c i , 0 α 0 + c i , 1 α + + c i , m 1 α m 1 ,其中, c i , j G F ( p ) 。当 1 t < m 时,设 G 1 = { β 0 = 0 , β 1 , , β P t 1 } 是由 α 0 , α 1 , , α t 的线性组合生成的 G F ( q ) 上的加法子群,即 β i = c i , 0 α 0 + c i , 1 α + + c i , t α t ,其中, c i , j G F ( p ) 。设 G 2 = { γ 0 = 0 , γ 1 , , γ P m t 1 } 是由 α t , α t + 1 , , α m 1 的线性组合生成的 G F ( q ) 上的加法子群,即 γ i = c i , t α t + c i , t + 1 α t + 1 + + c i , m 1 α m 1 ,其中, c i , j G F ( p ) 。显然, G 1 G 2 = { 0 } 。当 0 i < p m t 时,定义如下元素集合:

γ i + G 1 = { γ i , γ i + β 1 , , γ i + β P t 1 }

该集合其实就是以 γ i 为陪集首的 G 1 的陪集,包括其自身。这 P m t 个陪集构成了有限域 G F ( q ) 中q个元素的一个划分。 G 1 的两个陪集是不相交的。

构造 G F ( q ) 上大小为 p m t × p t 的矩阵:

B ( 2 ) = [ b 0 b 1 b m 1 ] = [ 0 β 1 β p t 1 γ 1 γ 1 + β 1 γ 1 + β p t 1 γ p m t 1 γ p m t + β 1 γ p m t + β p t 1 ] (6.1)

这里对于 0 i < P m t ,第i行由 G 1 的第i个陪集的元素组成。 G F ( q ) 的每一个元素在 B ( 2 ) 中出现且仅出现一次。除第一行, B ( 2 ) 中的每一行 b i 都由 G F ( q ) 的非零元素组成。 B ( 2 ) 的第一行第一个位置为 G F ( q ) 上的0元素。除第一列, B ( 2 ) 中的每一行都由 G F ( q ) 的非零元素组成。 B ( 2 ) 的第一列第一个位置为 G F ( q ) 上的0元素。由于 G 1 的任意两个陪集互不相交,故 B ( 2 ) 任意不同两行对应位置的元素也互不相同。

推论2:矩阵 B ( 2 ) 满足 2 × 2 子矩阵约束。

证明:基于定理1和定理2,对于矩阵 B ( 2 ) 的任意 2 × 2 子矩阵 B 2 × 2 ( 2 ) ,其元素属于第 i 1 行、第 i 2 行和第 j 1 列、第 j 2 列,其中 0 i 1 , i 2 P m t 0 j 1 , j 2 P t i 1 i 2 j 1 j 2 ,其行列式 det ( B 2 × 2 ( 2 ) ) = ( γ i 1 + β j 1 ) ( γ i 2 + β j 2 ) ( γ i 1 + β j 2 ) ( γ i 2 + β j 1 ) = ( γ i 1 γ i 2 ) ( β j 1 β j 2 ) 0 ,因此,矩阵 B ( 2 ) 满足 2 × 2 子矩阵约束。

修饰矩阵 W b a s e 的构造:

基于定理 ,选取一个码率为1⁄3的结构化不规则LDPC码的原模图,其相对应的原模矩阵为:

V = [ 1 1 1 0 1 1 0 1 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 1 1 0 1 1 1 0 1 0 0 0 1 1 1 1 1 ] ,

为了保持原结构化不规则分组码的度分布,使用展开技术(unwrapping technique)来构造修饰矩阵 W b a s e 。将V分解为两个矩阵 V 1 V 2 ,使得 V = V 1 + V 2 ,其中

V 1 = [ 1 0 1 0 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 1 1 0 0 1 0 0 0 0 0 0 1 0 1 0 1 ] ,

V 2 = [ 0 1 0 0 0 1 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 0 0 1 0 1 0 ] ,

W b a s e = [ V 1 V 2 V 1 V 2 V 1 V 2 V 1 V 2 ] (6.2)

下面以几个例子来说明基于有限域的加法子群SC­QC­LDPC码的构造。

6.1. 一类二进制SC­QC­LDPC码

例6.1.

基矩阵B的选取:将 G F ( 2 8 ) 作为构造码的域。令 α G F ( 2 8 ) 的本原元( G F ( 2 ) 的一个扩域)。设 t = 5 ,则 m t = 3 。设 G 1 G 2 G F ( 2 8 ) 上加法群的两个子群,分别由 { α 0 , α 1 , α 2 , α 3 , α 4 } { α 5 , α 6 , α 7 } 张成。因此 G 1 = { β 0 = 0 , β 1 , , β 31 } G 2 = { γ 0 = 0 , γ 1 , , γ 7 } 。利用式(6.1),我们可以构造一个8 × 32基矩阵B。

修饰矩阵 W b a s e 的选取:选择 L = 26 ,利用式(6.2),我们可以构造一个162 × 208修饰矩阵 W b a s e

二进制SC­QC­LDPC码的构造:利用算法1给出的构造方法,可形成一个 G F ( 2 ) 上20,574 × 26,416矩阵 H s c , q c , b ,其零空间给出一个二进制SC­QC­LDPC码 C s c , q c , b

6.2. 一类多进制SC­QC­LDPC码

例6.2.

基矩阵B的选取:将 G F ( 2 6 ) 作为构造码的域。令 α G F ( 2 6 ) 的本原元( G F ( 2 ) 的一个扩域)。设 t = 4 ,则 m t = 2 。设 G 1 G 2 G F ( 2 6 ) 上加法群的两个子群,分别由 { α 0 , α 1 , α 2 , α 3 } { α 4 , α 5 } 张成。因此 G 1 = { β 0 = 0 , β 1 , , β 15 } G 2 = { γ 0 = 0 , γ 1 , γ 2 , γ 3 } 。利用式(6.1),我们可以构造一个4 × 16基矩阵B。

修饰矩阵 W b a s e 的选取:选择 L = 26 ,利用式(6.2),我们可以构造一个162 × 208修饰矩阵 W b a s e

多进制SC­QC­LDPC码的构造:利用算法2给出的构造方法,可形成一个 G F ( q ) 上10,206 × 13,104矩阵 H s c , q c , q ,其零空间给出一个q进制SC­QC­LDPC码 C s c , q c , q

7. 基于有限域的乘法子群构造SC­QC­LDPC码

基矩阵B的构造:令 α G F ( q ) 的本原元。假设 q 1 不是素数。令 β = α k γ = α m ,其中k、m互素且 q 1 = k m ,则 M 1 = { β 0 = 1 , β , , β m 1 } M 2 = { γ 0 = 1 , γ , , γ k 1 } 构成了 G F ( q ) 上的两个循环乘法子群,且 M 1 M 2 = { 1 } 。当 0 i < k ,集合 γ i M 1 = { γ i , γ i β , , γ i β m 1 } M 1 的一个乘法陪集。循环子群 M 1 共有k个乘法陪集,包括其自身。当 0 j < m ,集合 β j M 2 = { β j , β j γ , , β j γ k 1 } M 2 的一个乘法陪集。循环子群 M 2 共有m个乘法陪集,包括其自身。构造 G F ( q ) 上大小为 k × m 的矩阵:

B ( 3 ) = [ b 0 b 1 b m 1 ] = [ 0 β 1 β m 1 1 γ 1 γ β 1 γ β m 1 1 γ k 1 1 γ k 1 β 1 γ k 1 β m 1 1 ] (7.1)

其中,第i行的分量可通过将 M 1 的第i个乘法陪集中每一个元素减去 G F ( q ) 中的元素1得到的;第j行的分量可通过将 M 2 的第j个乘法陪集中每一个元素减去 G F ( q ) 中的元素1得到的。

由上式,可以看出,矩阵 B ( 1 ) 具有以下结构性质:1) 有且仅有一个0元素,其位于左上角,其余 k m 1 个元素均是 G F ( q ) 上的非零元素;2) 任意两行的m个对应位置上的元素均不相同;3) 任意两列的k个对应位置上的元素均不相同;4) 所有元素均是 G F ( q ) 上的不同元素。

推论3:矩阵 B ( 3 ) 满足 2 × 2 子矩阵约束。

证明:基于定理1和定理2,对于矩阵 B ( 3 ) 的任意 2 × 2 子矩阵 B 2 × 2 ( 3 ) ,其元素属于第 i 1 行、第 i 2 行和第 j 1 列、第 j 2 列,其中 0 i 1 , i 2 k 0 j 1 , j 2 m i 1 i 2 j 1 j 2 ,其行列式 det ( B 2 × 2 ( 3 ) ) = ( γ i 1 β j 1 1 ) ( γ i 2 β j 2 1 ) ( γ i 1 β j 2 1 ) ( γ i 2 β j 1 1 ) = ( γ i 1 γ i 2 ) ( β j 2 β j 1 ) 0 ,因此,矩阵 B ( 3 ) 满足 2 × 2 子矩阵约束。

修饰矩阵 W b a s e 的构造:

基于定理4,选取一个结构化不规则LDPC码的原模图,其相对应的原模矩阵为:

V * = [ 1 2 1 0 1 1 0 1 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 1 1 0 1 1 1 0 1 0 0 0 1 1 1 1 1 ]

将多边原模图分解成二进制矩阵之和,将V分解为 V 3 V 4 ,使得 V * = V 3 + V 4 ,其中

V 3 = [ 1 1 1 0 1 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 1 1 0 0 1 0 0 0 0 0 0 1 0 1 0 1 ] ,

V 4 = [ 0 1 0 0 0 1 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 0 0 1 0 1 0 ] ,

利用使用展开技术(unwrapping technique)来构造修饰矩阵 W b a s e

W b a s e = [ V 3 V 4 V 3 V 4 V 3 V 4 V 3 V 4 ] . (7.2)

下面以几个例子来说明基于有限域的乘法子群SC­QC­LDPC码的构造。

7.1. 一类二进制SC­QC­LDPC码

例7.1.

基矩阵B的选取:将 G F ( 2 8 ) 作为构造码的域。令 α G F ( 2 8 ) 的本原元( G F ( 2 ) 的一个扩域)。可以看到, 2 8 1 = 255 ,其可以分解为5和51的乘积,且它们互素。令 β = α 5 γ = α 51 。集合 M 1 = { β 0 = 1 , β , , β 50 } M 2 = { γ 0 = 1 , γ , , γ 4 } 构成了 G F ( q ) 上的两个循环乘法子群,且 M 1 M 2 = { 1 } 。利用式(7.1),我们可以构造一个5 × 51基矩阵B。

修饰矩阵 W b a s e 的选取:选择 L = 26 ,利用式(7.2),我们可以构造一个162 × 208修饰矩阵 W b a s e

二进制SC­QC­LDPC码的构造:利用算法1给出的构造方法,可形成一个 G F ( 2 ) 上20,574 × 26,416矩阵 H s c , q c , b ,其零空间给出一个二进制SC­QC­LDPC码 C s c , q c , b

7.2. 一类多进制SC­QC­LDPC码

例7.2.

基矩阵B的选取:将 G F ( 2 6 ) 作为构造码的域。令 α G F ( 2 6 ) 的本原元( G F ( 2 ) 的一个扩域)。可以看到, 2 6 1 = 63 ,其可以分解为7和9的乘积,且它们互素。令 β = α 7 γ = α 9 。集合 M 1 = { β 0 = 1 , β , , β 8 } M 2 = { γ 0 = 1 , γ , , γ 6 } 构成了 G F ( q ) 上的两个循环乘法子群,且 M 1 M 2 = { 1 } 。利用式 ,我们可以构造一个7 × 9基矩阵B。

修饰矩阵 W b a s e 的选取:选择 L = 26 ,利用式(7.2),我们可以构造一个162 × 208修饰矩阵 W b a s e

多进制SC-QC-LDPC码的构造:利用算法2给出的构造方法,可形成一个 G F ( q ) 上20,574 × 26,416矩阵 H s c , q c , q ,其零空间给出一个q进制SC-QC-LDPC码 C s c , q c , q

8. 结论

本文提出了一种构造二进制和多进制SC-QC-LDPC码的一般方法,所提出的方法是基于在有限域 G F ( q ) 上对满足 2 × 2 子矩阵约束的修饰基矩阵的阵列散列,由此导出的稀疏校验矩阵所对应的Tanner图的围长至少为6。基于一般方法,本文提出了3种构造二进制和多进制SC-QC-LDPC码的特定方法,由此构造出了6类二进制SC-QC-LDPC码和多进制SC-QC-LDPC码。

参考文献

[1] Felstrom, A.J. and Zigangirov, K. (1999) Time-Varying Periodic Convolutional Codes with Low-Density Parity-Check Matrix. IEEE Transactions on Information Theory, 45, 2181-2191.
https://doi.org/10.1109/18.782171
[2] Gallager, R.G. (1962) Low-Density Parity-Check Codes. IRE Transactions on Information Theory, 8, 21-28.
https://doi.org/10.1109/TIT.1962.1057683
[3] Kudekar, S., Richardson, T. and Urbanke, R. (2011) Threshold Saturation via Spatial Coupling: Why Convolutional LDPC Ensembles Perform So Well over the BEC. IEEE Transactions on Information Theory, 57, 803-834.
https://doi.org/10.1109/TIT.2010.2095072
[4] Kudekar, S., Measson, C., Richardson, T. and Urbankez, R. (2010) Threshold Saturation on BMS Channels via Spatial Coupling. 6th International Symposium on Turbo Codes and Iterative Information Processing, Brest, 6-10 September 2010, 309-313.
https://doi.org/10.1109/ISTC.2010.5613887
[5] Yedla, A., Jian, Y.-Y., Nguyen, P. and Pfifister, H. (2012) A Simple Proof of Threshold Saturation for Coupled Scalar Recursions. 2012 7th International Symposium on Turbo Codes and Iterative Information Processing, Gothenburg, 27-31 August 2012, 51-55.
https://doi.org/10.1109/ISTC.2012.6325197
[6] Lentmaier, M., Sridharan, A., Costello Jr., D.J. and Zigangirov, K. (2010) Iterative Decoding Threshold for LDPC Convolutional Codes. IEEE Transactions on Information Theory, 56, 5274-5289.
https://doi.org/10.1109/TIT.2010.2059490
[7] Huang, K., Mitchell, D.G.M., Wei, L., Ma, X. and Costello, D. (2015) Performance Comparison of LDPC Block and Spatially Coupled Codes over GF(q). IEEE Transactions on Communications, 63, 592-604.
https://doi.org/10.1109/TCOMM.2015.2397433
[8] Iyengar, A., Papaleo, M., Siegel, P., Wolf, J., Vanelli-Coralli, A. and Corazza, G. (2012) Windowed Decoding of Protograph-Based LDPC Convolutional Codes over Erasure Channels. IEEE Transactions on Information Theory, 58, 2303-2320.
https://doi.org/10.1109/TIT.2011.2177439
[9] Andriyanova, I. and Graelli Amat, A. (2013) Threshold Saturation for Nonbinary SC-LDPC Codes on the Binary Erasure Channel. IEEE Transactions on Information Theory, 62, 2622-2638.
https://doi.org/10.1109/TIT.2016.2540800
[10] Pusane, A., Smarandache, R., Vontobel, P. and Costello, D. (2011) Deriving Good LDPC Convolutional Codes from LDPC Block Codes. IEEE Transactions on Information Theory, 57, 835-857.
https://doi.org/10.1109/TIT.2010.2095211
[11] Fossorier, M. (2004) Quasi Cyclic Low-Density Parity-Check Codes from Circulant Permutation Matrices. IEEE Transactions on Information Theory, 50, 1788-1793.
https://doi.org/10.1109/TIT.2004.831841
[12] Li, Z., Chen, L., Zeng, L., Lin, S. and Fong, W. (2005) Efficient Encoding of Quasi-Cyclic Low-Density Parity-Check Codes. IEEE Transactions on Communications, 53, 1973-1973.
https://doi.org/10.1109/TCOMM.2005.858628
[13] Chen, Y. and Parhi, K. (2004) Overlapped Message Passing for Quasi-Cyclic Low-Density Parity Check Codes. IEEE Transactions on Circuits and Systems I: Regular Papers, 51, 1106-1113.
https://doi.org/10.1109/TCSI.2004.826194
[14] Uchikawa, H. Kasai, K. and Sakaniwa, K. (2011) Design and Performance of Rate-Compatible Non-Binary LDPC Convolutional Codes. IEICE Transactions on Fundamentals of Electronics, Communications and Computer Sciences, E94, 2135-2143.
[15] Ryan, W. and Lin, S. (2009) Channel Codes: Classical and Modern. Cambridge University Press, Cambridge.
https://doi.org/10.1017/CBO9780511803253
[16] Lan, L., Zeng, L., Tai, Y., Chen, L., Lin, S. and Abdel-Ghaffar, K. (2007) Construction of Quasi-Cyclic LDPC Codes for AWGN and Binary Erasure Channels: A Finite Field Approach. IEEE Transactions on Information Theory, 53, 2429-2458.
https://doi.org/10.1109/TIT.2007.899516
[17] Zhang, L., Huang, Q., Lin, S., Abdel-Ghaffar, K. and Blake, I.F. (2010) Quasi-Cyclic LDPC Codes: An Algebraic Construction, Rank Analysis, and Codes on Latin Squares. IEEE Transactions on Communications, 58, 3126-3139.
https://doi.org/10.1109/TCOMM.2010.091710.090721
[18] Dolecek, L., Divsalar, D., Sun, Y. and Amiri, B. (2014) Non-Binary Protograph-Based LDPC Codes: Enumerators, Analysis, and Designs, Information Theory. IEEE Transactions on Information Theory, 60, 3913-3941.
https://doi.org/10.1109/TIT.2014.2316215
[19] Chandrasetty, V.A., Johnson, S.J. and Lechner, G. (2013) Memory Efficient Decoders Using Spatially Coupled Quasi-Cyclic LDPC Codes. Clinical Orthopaedics and Related Research, arXiv:1305.5625
[20] Hu, X.-Y., Eleftheriou, E. and Arnold, D.-M. (2005) Regular and Irregular Progressive Edge-Growth Tanner Graphs. IEEE Transactions on Information Theory, 51, 386-398.
https://doi.org/10.1109/TIT.2004.839541
[21] Li, J., Liu, K., Lin, S., Abdel-Ghaffar, K. and Ryan, W.E. (2015) An Unnoticed Strong Connection between Algebraic-Based and Protograph-Based LDPC Codes. 2015 Information Theory and Applications Workshop, San Diego, 1-6 February 2015, 36-45.
https://doi.org/10.1109/ITA.2015.7308964
[22] Wei, L., Koike-Akino, T., Mitchell, D., Fuja, T. and Costello, D.J. (2014) Threshold Analysis of Non-Binary Spatially-Coupled LDPC Codes with Windowed Decoding, in Information Theory (ISIT). 2014 IEEE International Symposium on, Honolulu, 29 June-4 July 2014, 881-885.
https://doi.org/10.1109/ISIT.2014.6874959