一种非规则QC-LDPC码的低错误平层构造方法与流程

一种非规则qc-ldpc码的低错误平层构造方法
技术领域
1.本发明属于信道编码技术领域，涉及一种基于基本陷阱集(elementary trapping sets,ets) 消除与围长约束的非规则准循环低密度奇偶检验(quasi-cyclic low-density parity-check,qc
‑ꢀ
ldpc)码构造方法。


背景技术：

2.近年来，低密度奇偶检验(low-density parity-check,ldpc)码的优越性得到国内外科研工作者关注，并且已成为现代通信系统不可或缺的部分，被用来检测和修正由信道效应如噪声、衰减和干扰等引起的信息传输错误。然而，其性能提高的同时，编码复杂度也同样提高了，进而导致实际应用中成本增加和资源浪费。为了解决该问题，国内外学者提出了qc-ldpc码， qc-ldpc码因其高效并行硬件实现、较低的编译码复杂度和出色的纠错性能成为ldpc码的一个重要类别。通过在qc-ldpc码中引入非规则的度分布，提高了qc-ldpc码的瀑布区性能，但这种瀑布区的性能提升往往是以较明显的错误平层为代价的。这是因为非规则度分布的引入增加了tanner图中导致译码失败的陷阱集(trapping sets,ts)等拓扑结构，这类陷阱集等拓扑结构会引起非规则qc-ldpc码出现错误平层现象。非规则qc-ldpc码在光链路，储存设备，空间通信等中具有广泛应用。因此如何消除陷阱集，构造低错误平层的非规则qc
‑ꢀ
ldpc码是一项具有挑战性的任务。
3.近年来，许多学者针对错误平层现象进行了大量研究，包括通过测量围长、短环多样性以及短环连接方式等方式间接构造低错误平层ldpc码。也有通过消除陷阱集直接构造低错误平层ldpc码的方法，如naseri通过约束tanner图中8环的类型而避免规则qc-ldpc码中的小陷阱集。但此方法不适用于非规则qc-ldpc码，且对扩展系数提出了要求。hashemi提出了非规则qc-ldpc码tanner图中基本陷阱集的dot,path,lollipop(dpl)扩展方法用于搜索无叶基本陷阱集(leafless elementary trapping sets,lets)和有叶基本陷阱集(elementary trapping setswith leaf,etsl)。此方法能对基本陷阱集进行穷尽搜索但搜索复杂度很高。karimi在全连通基矩阵上构造不含无叶基本陷阱集(leafless elementary trapping sets,lets)的规则qc-ldpc 码。之后又针对非规则qc-ldpc提出基于dot,path,lollipop(dpl)扩展和lets结构中的父子关系的无叶基本陷阱集的最小父结构搜索方法。此搜索方法相较于dpl扩展复杂度有所降低，但对于tanner图中lets的父结构识别筛选过程依然复杂。
4.本发明针对非规则qc-ldpc码中的搜索消除ets的过程复杂度较高的问题，通过度分布的巧妙选择和围长约束，不需对目标范围内的ets进行dpl扩展搜索，只对由简单环构成的基本陷阱进行简单的搜索和消除，便能消除设置范围内的绝大多数基本陷阱集。此方法计算复杂度相对较低，且仿真结果表明，利用本发明提出的方法构造的qc-ldpc码型，具有良好的纠错性能，并且未出现明显错误平层现象。


技术实现要素：

5.有鉴于此，本发明的目的为提供一种基于基本陷阱集消除与围长约束的非规则qc-ldpc 码低错误平层构造方法，通过度分布的巧妙选择和围长约束，不需对目标范围内的ets进行 dpl扩展搜索，只对由简单环构成的基本陷阱进行简单的搜索和消除，便能消除设置范围内的绝大多数基本陷阱集。该方法除了能够改善高信噪比区域的错误平层，还具有码长、码率的任意可设性。因而该方法能满足通信系统对纠错码具有码率可灵活选择、低错误平层和计算复杂低的需求为达到上述目的，本发明提供如下技术方法：
6.首先，根据在构造基矩阵时可以通过设置不含奇数度数的度分布避免任意b值为奇数的 ets出现，以减少ets搜索计算量这一理论，以基矩阵大小是24
×
48前提将初始度分布为λ1(x)＝0.38345x 0.04237x2 0.57409x3变量节点集中的第27、28位变量节点度数由3改为 4。同时第25、26位变量节点度数由2改为4，得到新的校验节点度分布：λ2(x)＝0.33334x 0.66666x3。由此度分布构造的基矩阵中不存在任意b为奇数的(a,b)ets，由于算法1已经确保基矩阵中不存在a≤8,b≤4的(a,b)简单环lets。但a≤8,b＞4的简单环 lets依旧可能存在，并且这些简单环lets可能通过dpl扩展成为a≤8,b≤4范围内的lets。通过对扩展过程的分析可知，要使b值减小到能进入a≤8,b≤4中，必定会有4环和6环产生，因而本方法利用算法2中围长约束和ap序列扩展避免4环和6环的产生。
7.其次，选取了合适的校验节点度分布之后，结合基本陷阱集消除和围长约束对渐进边增长(progressive edge growth,peg)算法改进：peg算法在添加边时，为了最大化局部环长，每个变量节点的每条边都应该连接当前度最小的校验节点。通过添加基本陷阱集消除和围长约束到peg算法中，使每个变量节点连接最后一条边时进行简单环结构的基本陷阱集搜索，若搜索到a≤a
max
，b≤b
max
范围内的简单环基本陷阱集则在大小为a
max
×bmax
矩阵a的相应位置计数。对a中对应的各类基本陷阱集求和，并存放在b的对应行中。查询b中每行，找出基本陷阱集总数最小的行，该行对应的校验节点即为最优节点，计算当前最优节点的平均围长 g
av
，若平均围长值大于等于设置的平均围长阈值γ
av
，则连接当前节点。本专利中设置a
max
＝8， b
max
＝4，平均围长阈值γ
av
＝10。并使用改进的peg算法构造一个24
×
48的基矩阵m，构造出的基矩阵中没有a≤a
max
,b≤b
max
范围内的简单环基本陷阱集，且通过围长约束和ap序列扩展抑制基矩阵tanner图中4环和6环产生，从而抑制未在a≤a
max
,b≤b
max
范围内ets通过 dpl扩展进入a≤a
max
,b≤b
max
范围，就能确保构造的基矩阵中不存在设置范围内的绝大多数 ets，从而降低错误平层现象。
8.最后，基于等差(arithmetic progression,ap)序列构造一个24
×
48的循环移位矩阵p，从而通过循环移位矩阵p来扩展基本矩阵m，以此得到奇偶校验矩阵h，扩展的规则为：将m中的0用维数为q
×
q的零矩阵o替换，其中q在本专利里为50，1用维数为q
×
q的循环置换矩阵p
ij
替换，p
ij
由维数为q
×
q的单位矩阵e向右循环移位p
ij
次得到。p
ij
的取值来源于移位矩阵p中第i行第j列的元素，因为矩阵p的维数与基本矩阵m的维数相同，所以循环移位系数p
ij
对应于m矩阵中第i行第j列的元素。
9.最后，在相同的仿真环境下，将本专利所提出的非规则qc-ldpc码型构造方法与其它编码构造方法进行仿真对比分析。
10.本发明的有益效果在于：利用基于基本陷阱集消除和围长约束的非规则qc-ldpc码构造方法。通过度分布的巧妙选取，利用基本陷阱集搜索和围长约束以改进渐进边增长
peg算法构造基矩阵，然后通过ap序列扩展得到所需的校验矩阵。该方法仅需对简单环形式的ets 进行搜索和消除，就能确保构造的基矩阵中不存在设置范围内的绝大多数ets，从而降低错误平层现象，且该方法计算复杂度相对较低，可灵活设计码长码率。仿真结果表明，所构造的peg-trap sets-ave-ap(ptaa)-qc-ldpc(2400,1200)码纠错性能要优于与peg-trap sets
‑ꢀ
ap(ptap)-qc-ldpc(2400,1200)码、ieee802.16标准中qc-ldpc(2400,1200)码、采用peg与 ap序列所构造的peg-ap-qc-ldpc(2400,1200)码、非规则protograph-based(pb)qc
‑ꢀ
ldpc(2400,1200)码与和规则girth-8qc-ldpc(2400,1200)码，且没有明显的错误平层现象。
附图说明
11.为了使本发明的目的、技术方法和有益效果更加清楚，本发明提供如下附图进行说明：
12.图1为本发明算法的技术流程图；
13.图2为dv＝3的(4,4)ets
14.图3(a)为(4,2)lets(b)为(4,4)etsl1(c)为(4,4)etsl215.图4为本专利非规则qc-ldpc码与其他三种码型的性能对比图；
16.图5为本专利非规则qc-ldpc码与其他两种码型的性能对比图；
具体实施方式
17.下面将结合附图，对本发明的优选实施例进行详细的描述。
18.(1)结合附图2、3说明，为了方便理解基本陷阱集消除算法和围长约束，下面详细解释了校验矩阵tanner图中的陷阱集、陷阱集扩展、围长约束以及ap序列。
19.定义1.1：一个(a,b)陷阱集s是变量节点集合v的子集，s的大小即s中变量节点个数 a，与s相连的度数为奇数的校验节点个数为b。我们也称陷阱集s属于(a,b)类，a,b值越小，对纠错性能影响越大。特别的，如果与s相连的所有校验节点度数都不大于2，则s是一个 (a,b)基本陷阱集。图2是一个变量节点度数dv＝3的(4,4)基本陷阱集。白色矩形表示度数为偶数的校验节点，用γe(s)表示。黑色矩形表示度数为奇数的校验节点，用γo(s)表示。圆形表示变量节点。
20.定义1.2：将ets的tanner图g(s)中的γo(s)移除，将γe(s)用直线替代，同时去除γo(s) 的连接线，得到的图称为s的正态图(normal graph)，用表示。
21.定义1.3：若ets中的每个变量节点都至少连接两个γe(s)，则此基本陷阱集为无叶基本陷阱集。
22.定义1.4：若ets中至少有一个变量节点只连接一个γe(s)，则此ets为etsl。若etsl 中至少含有一个环则是etsl1，否则是etsl2。附图3(a),(b),(c)分别为(4,2)lets,(4,4)etsl1, (3,3)etsl2与各自的正态图。
23.定义1.5：tanner图中从任意变量节点出发，依次经过校验节点和变量节点，最终回到出发节点的路径称为环。若环中没有弦(chord)存在，此环称为简单环。记环中变量节点数量为 d。则环长为l＝2d。tanner图中最短的环长为围长。通过一个变量节点的最短环长度称为变量节点的局部围长。
24.定义1.6：不在环中，但与环中的任意两个点相邻的一条边称为弦。
25.定义1.7：校验矩阵的平均围长(girth average,ga)定义如式(1)所示：
[0026][0027]
其中，k为第i个变量节点局部围长的一半，l
max
为所有变量节点的最大围长，gi(l)为局部围长为l的变量节点占总变量节点的百分比。
[0028]
对校验矩阵进行平均围长ga约束可以有效降低矩阵中短环的比重，从而改善纠错性能。
[0029]
引理1.1：对非规则tanner图中的(a,b)lets，分别使用pam(m≥2)和扩展会得到(a 1,b k-2m)，和其中d
vi
,i＝1,...,m,是扩展的变量节点度数。
[0030]
引理1.2：在度分布为λ(x)的非规则tanner图中，对(a,b)lets使用 pam(m≥2)和得到(a',b')lets，b'的最小值分别等于 (b min{z-2a,-y})

，符号 (f)

＝max(f,0)，y是λ(x)的最大值且小于等于a，z是λ(x)最小值且大于a，时， d'是大于2的λ(x)最小值，时，
[0031]
引理1.2给出了一个(a,b)lets通过dpl扩展之后b'的下限值，不难看出通过 pam(m≥2)和扩展之后值分别最多减小2和1，而通过扩展后的b值能减少更多。在非规则ldpc码中除了lets，etsl的存在也是错误平层产生的原因。前人已经证明了非规则qc-ldpc码的错误平层中没有etsl2的影响，同时基于etsl1的扩展原理提出引理1.3说明不需要单独进行etsl1的搜索和消除。
[0032]
引理1.3：在最小度分布δ(g)≥2的非规则tanner图g中如果不存在a≤a
max
，b≤b
max
的 lets，那么g中也不存在a≤a
max
1，b≤b
max
的etsl1。
[0033]
由以上分析可知，想要消除非规则qc-ldpc码的错误平层现象只需要进行lets的搜索和消除即可。
[0034]
引理1.4：不含奇数度分布的非规则tanner图g中不存在任意b为奇数的(a,b)ets。
[0035]
由引理1.4可知，在构造校验矩阵时可以通过设置不含奇数度数的度分布避免任意b为奇数的(a,b)ets出现，以减少ets搜索计算量。
[0036]
引理1.5：如果基矩阵中各个元素差值d
i,j
满足式(2)，则对任意m
×
n的基矩阵，其对应的 tanner图的围长至少为8。
[0037]
[0038]
式中，d
i,j
＝p(i,j 1)-p(i.j)，1≤i≤m，1≤j≤n-1。d
i,j
表示基矩阵中第i行第j列与前一个元素之差，d是一个可任意设置的值，为第一行元素的公差。由引理1.5可知，通过等差 (arithmetic progression,ap)序列构造的基矩阵对应的tanner图的围长至少为8。
[0039]
(2)结合附图4、5说明，通过基本陷阱集搜索、度分布的巧妙选取和围长约束以改进peg 算法构造基矩阵，然后通过ap序列扩展得到所需的校验矩阵。该方法仅需对简单环形式的 ets进行搜索和消除，就能确保构造的基矩阵中不存在设置范围内的绝大多数ets，从而降低错误平层现象，且构造出的码型有较好的纠错性能。基本陷阱集消除算法和改进后的peg 算法伪代码如下：
[0040]
算法1：基本陷阱集消除算法
[0041][0042][0043]
算法2：结合基本陷阱集搜索与围长约束的改进peg算法
standard802.16.ieee standard for air interface for broadband wireless access systems[s].2018”中非规则qc-ldpc(2400,1200)码型和文献[2]“袁建国,王宏森,张希瑞,赖春红.一种基于ap与消除陷集的低错误平层qc-ldpc码构造方法[j].半导体光电,2020,41(01):146-150.”中非规则 ptap-qc-ldpc(2400,1200)码型进行仿真对比分析，其仿真结果如附图4所示。
[0049]
仿真表明：在误码率为10-6
时，本方法所构造的非规则ptaa-qc-ldpc(2400,1200)码相较于ptap-qc-ldpc(2400,1200)码的净编码增益(net coding gain，ncg)约提升0.15db，相较于ieee 801.16标准中qc-ldpc(2400,1200)码的ncg约提升0.25db，相较于peg-qc
‑ꢀ
ldpc(2400,1200)码的ncg约提升0.40db。且随着信噪比增加，ncg提升会进一步增大。
[0050]
为了进一步论证本方法所述构造方法的优越性，将本方法构造的非规则ptaa-qc
‑ꢀ
ldpc(2400,1200)码型与文献[3]“karimi b,banihashemi a h.construction of irregularprotograph-based qc-ldpc codes with low error floor[j].ieee transactions oncommunications,2021,69(1):3-18.”中基于原模图所构造的非规则protograph-based(pb)qc
‑ꢀ
ldpc(2400,1200)码和文献[4]“naseri s,banihashemi a h.construction of girth-8 qc
‑ꢀ
ldpc codes free of small trapping sets[j].ieee communications letters,2019,23(11):1904
‑ꢀ
1908.”中消除小陷阱集所构造的规则girth-8qc-ldpc(2400,1200)码进行仿真对比分析，其仿真结果如附图5所示。
[0051]
仿真表明：本方法构造的非规则ptaa-qc-ldpc(2400,1200)码型与非规则pb qc
‑ꢀ
ldpc(2400,1200)码相较于规则girth-8qc-ldpc(2400,1200)码均有更好的纠错性能和瀑布区表现。在误码率为10-7
时，本方法所构造的非规则ptaa-qc-ldpc(2400,1200)码相较于非规则pb qc-ldpc(2400,1200)码以及规则girth-8qc-ldpc(2400,1200)码的ncg改善了约 0.05db和1.68db。且随着信噪比增加，ncg提升进一步增大。还可由附图5对比分析可知，相较于非规则pb qc-ldpc(2400,1200)码，本方法所构造的非规则ptaa-qc-ldpc(2400,1200) 码在构造过程中消除了设置范围内的所有ets。此外，相较于文献[3]中基于原模图对lets 结构的父子关系进行dpl扩展消除的构造方法，本方法所提出的构造方法在构造过程没有进行dpl扩展，因而其计算复杂度更低，同时因为围长约束也具有更好的纠错性能。
[0052]
因而，利用本方法所提出的构造方法所构造的非规则ptaa-qc-ldpc码具有较好的纠错性能且没有明显的错误平层现象。
[0053]
最后说明的是，以上优选实施例仅用以说明本发明的技术方案而非限制，尽管通过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。