CN107968657B - 一种适用于低密度奇偶校验码的混合译码方法 - Google Patents

一种适用于低密度奇偶校验码的混合译码方法 Download PDF

Info

Publication number
CN107968657B
CN107968657B CN201711216506.6A CN201711216506A CN107968657B CN 107968657 B CN107968657 B CN 107968657B CN 201711216506 A CN201711216506 A CN 201711216506A CN 107968657 B CN107968657 B CN 107968657B
Authority
CN
China
Prior art keywords
node
variable node
decoding method
check
information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711216506.6A
Other languages
English (en)
Other versions
CN107968657A (zh
Inventor
姜明
孙凯飞
赵春明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201711216506.6A priority Critical patent/CN107968657B/zh
Publication of CN107968657A publication Critical patent/CN107968657A/zh
Application granted granted Critical
Publication of CN107968657B publication Critical patent/CN107968657B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/03Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits
    • H03M13/11Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits using multiple parity bits
    • H03M13/1102Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1105Decoding
    • H03M13/1111Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms
    • H03M13/1117Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms using approximations for check node processing, e.g. an outgoing message is depending on the signs and the minimum over the magnitudes of all incoming messages according to the min-sum rule
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/03Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits
    • H03M13/11Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits using multiple parity bits
    • H03M13/1102Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1105Decoding
    • H03M13/1111Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms
    • H03M13/1125Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms using different domains for check node and bit node processing, wherein the different domains include probabilities, likelihood ratios, likelihood differences, log-likelihood ratios or log-likelihood difference pairs

Abstract

本发明公开了一种适用于低密度奇偶校验码的混合译码方法,通过设定相应阈值来根据校验节点度数对校验节点进行分类,度数小于阈值的节点采用基于线性插值的APPROXIMATE‑MIN*方法进行更新,度数大于阈值的节点采用归一化修正最小和译码方法进行计算更新,适用于码率兼容结构,较好的均衡实现复杂度和译码性能,既降低了复杂度,同时也减少了解码性能损失。同时本发明用固定因子乘性修正解码输入信息,消除了译码方法对高斯信道方差估计的依赖,较好的均衡实现复杂度和译码性能。本发明缩小了BP译码方法和MS译码方法性能的差距,消除了解码方法对高斯信道方差估计的依赖。

Description

一种适用于低密度奇偶校验码的混合译码方法
技术领域
本发明属于译码方法技术领域,涉及一种适用于低密度奇偶校验码的混合译码方法。
背景技术
LDPC码即低密度奇偶校验(Low-Density Parity-Check,LDPC)码,它是由RobertG.Gallager博士于1963年首次提出的一类具有稀疏校验矩阵的线性分组码,并在1995被D.J.C.Mackay和R.W.Neal重新提出。LDPC码具有逼近Shannon限的良好性能,而且译码复杂度较低、结构灵活、低延时,可实现高速并行译码,近二十年来一直是信道编码领域的研究热点。目前LDPC码已广泛应用于无线通信(WiFi,WiMAX)、深空通信、光纤通信、卫星数字视频和音频广播等领域。现在LDPC码已经成为5G信道编码中数据信道的编码方案,eMMB场景已经确定使用LDPC码的编码方案。
除了关键的eMMB场景,5G uRLLC场景对编码的吞吐率和处理延时提出较高要求。5G系统对LDPC码编码方案要求支持的码率范围为1/6-8/9,对于码率小于1/3的低码率LDPC码要具有更低的解码延时。
G.Gallager在提出LDPC码时给出了两种迭代译码方法。其中比特翻转(bitflipping,BF)译码方法属于硬判决译码方法,该译码方法译码复杂度低,易于硬件实现,但译码性能较差;软判决译码是基于后验概率的置信传播(belief propagation,BP)译码方法,该方法在性能上可以逼近最大似然译码,但其译码复杂度很高。目前LDPC码的主要解码方案包括BP类译码方法及BP类简化方法,如MS(Min-Sum,MS)类方法。由于BP译码方法的实现复杂度较高,在实际应用中使用较多的为BP类简化方法。BP类简化方法主要分为两类,对BP方法的简化和MS类简化方法。
MS类译码方法仅通过选择输入变量节点信息的最小值和次小值来更新检验节点信息,降低了译码方法的计算复杂度,但会带来性能上的损失。为了降低MS译码方法的性能损失,一些针对MS译码方法的修正方案被提了出来。MS译码方法的性能损失来源于对校验节点单元信息的过高估计。因此归一化修正最小和译码方法(Normalized Min-Sum,NMS)和偏移最小和译码方法(Offset Min-Sum,OMS)通过对MS方法校验节点输出进行乘性和加性修正来提升解码性能。同时,一些译码方法针对NMS和OMS方法的修正因子来进一步NMS和OMS方法性能。
APPROXIMATE-MIN*译码方法是一类对BP简化的方法,该方法是通过选取输入的变量节点信息的最小值,然后对校验节点更新进行更精确的修正,性能更接近BP译码方法。但该方法对解码输入似然比计算精度要求较高,依赖信道方差估计的,具体实现复杂度太高。
因此,现有的BP译码方法及其简化方法复杂度太高,而MS译码方法虽然降低了复杂度,但其性能低,无法满足需求。
发明内容
为解决上述问题,本发明公开了一种适用于低密度奇偶校验码的混合译码方法,缩小了BP译码方法和MS译码方法性能的差距,消除了解码方法对高斯信道方差估计的依赖,同时适用于码率兼容结构,较好地均衡了实现复杂度和译码性能,既降低了复杂度,同时也减少了解码性能的损失。
为了达到上述目的,本发明提供如下技术方案:
一种适用于低密度奇偶校验码的混合译码方法,包括如下步骤:
步骤一:对每个变量节点进行迭代初始化,并在变量节点中引入常数固定乘性修正因子;
步骤二:迭代次数加一;
步骤三:对于各个校验节点,利用第k-1次迭代中产生的由相应变量节点传递而来的信息更新该校验节点向相应变量节点传递的信息:度数小于阈值的节点采用基于线性插值的APPROXIMATE-MIN*方法进行更新,度数大于阈值的节点采用归一化修正最小和译码方法进行计算更新;
步骤四:对于各个变量节点,根据由相应的校验节点传递而来的信息和该变量节点相对应的初始输入似然比信息,计算该变量节点的输出似然比信息;
步骤五:对每个变量节点进行硬判决并计算校正子;
步骤六:如果满足条件则结束译码,输出码字;否则顺序执行下一步骤;
步骤七:对于各个变量节点,利用由相应校验节点传递而来的信息和步骤四中计算出来的该变量节点的输出似然比信息,计算该变量节点向相连的各校验节点传递的信息;
顺序执行步骤二,直至迭代完毕。
进一步的,所述步骤一具体包括:
初始化迭代次数k=0;对每个变量节点vn,n∈[1,N]初始化:
Figure BDA0001485603320000021
其中β为常数固定乘性修正因子;Lnm为变量节点vn传递给校验节点cm的信息;yn为接收序列中的元素。
进一步的,所述步骤三具体包括:
对于各个校验节点cm(m∈[1,M]),利用第k-1次迭代中产生的由相应变量节点传递而来的信息
Figure BDA0001485603320000031
更新该校验节点向相应变量节点vn(n∈B(m))传递的信息
Figure BDA0001485603320000032
Figure BDA0001485603320000033
Figure BDA0001485603320000034
其中T为设定的阈值;dc为矩阵每行中非零元素的个数,Lmn为校验节点cm传递给变量节点vn的信息;
根据第k-1次迭代中产生的由相应变量节点传递而来的信息
Figure BDA0001485603320000035
找到
Figure BDA0001485603320000036
绝对值的最小值的下标索引jmin,初始化局部变量δ0=∞;对每个h∈{1...dc}\jmin,计算
Figure BDA0001485603320000037
Figure BDA0001485603320000038
Figure BDA0001485603320000039
其中,B(m)\n为除变量节点vn的B(m)集合,
Figure BDA00014856033200000310
Figure BDA00014856033200000311
归一化修正最小和方法,
Figure BDA00014856033200000312
α为修正因子。
进一步的,所述步骤四具体包括:
对于各个变量节点vn(n∈[1,N]),根据由相应的校验节点传递而来的信息
Figure BDA00014856033200000313
和该变量节点相对应的初始输入似然比信息
Figure BDA00014856033200000314
计算该变量节点vn(n∈[1,N])的输出似然比信息
Figure BDA00014856033200000315
Figure BDA00014856033200000316
其中,A(n)为变量节点vn参与的校验节点集合。
进一步的,所述步骤五具体包括:
对每个变量节点vn(n∈[1,N])进行硬判决并计算校正子:
Figure BDA0001485603320000041
Sk=WkHT
进一步的,所述步骤六中条件为
Figure BDA0001485603320000042
或k=K,输出的码字
Figure BDA0001485603320000043
进一步的,所述步骤七中计算该变量节点向相连的各校验节点传递的信息的过程具体包括:
对于各个变量节点vn(n∈[1,N]),利用由相应校验节点cm(m∈A(n))传递而来的信息
Figure BDA0001485603320000044
和步骤四中计算出来的该变量节点vn(n∈[1,N])的输出似然比信息
Figure BDA0001485603320000045
计算该变量节点向相连的各校验节点cm(m∈A(n))传递的信息
Figure BDA0001485603320000046
Figure BDA0001485603320000047
与现有技术相比,本发明具有如下优点和有益效果:
1.本发明通过设定相应阈值来根据校验节点度数对校验节点进行分类,度数小于阈值的节点采用基于线性插值的APPROXIMATE-MIN*方法进行更新,度数大于阈值的节点采用归一化修正最小和译码方法进行计算更新,适用于码率兼容结构,较好的均衡实现复杂度和译码性能,既降低了复杂度,同时也减少了解码性能损失。
2.本发明用固定因子乘性修正解码输入信息,消除了译码方法对高斯信道方差估计的依赖,较好的均衡实现复杂度和译码性能。
3.本译码方法缩小了最小和译码方法与BP译码方法的解码性能差距。
附图说明
图1为本发明提供的适用于低密度奇偶校验码的混合译码方法流程图。
图2为一个1/4码率的LDPC码(Code1-R14-N8448)校验矩阵的分块示意图。
图3为Code1-R14-N8448码分别采用BP译码方法、NMS译码方法、以及本发明所述的一种适用于低密度奇偶校验码的混合译码方法进行仿真的误帧率性能曲线。
图4为一个1/5码率的LDPC码(Code2-R15-N10496)校验矩阵的分块示意图。
图5为Code2-R15-N10496码分别采用BP译码方法、NMS译码方法、以及本发明所述的一种适用于低密度奇偶校验码的混合译码方法进行仿真的误帧率性能曲线。
具体实施方式
以下将结合具体实施例对本发明提供的技术方案进行详细说明,应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
本发明提出的适用于低密度奇偶校验码的混合译码方法,针对置信传播类解码方法与最小和类简化解码方法性能存在明显差距,提出一种混合型最小和解码方法。本译码方法根据校验节点度数的不同,设定对应的阈值,对度数不同的校验节点更新采用不同的更新计算方法。对于度数大于阈值的校验节点采用归一化修正最小和译码方法更新,度数小于阈值的校验节点采用基于线性插值的APPROXIMATE-MIN*译码更新方法,适用于码率兼容结构,降低运算复杂度。同时由于APPROXIMATE-MIN*译码方法对解码输入似然比计算精度要求较高,依赖对信道方差的估计,具体实现复杂度太高。本发明改用固定因子乘性修正解码输入信息,消除解码方法对高斯信道方差估计的依赖,较好的均衡实现复杂度和译码性能,既降低了复杂度,同时也减少了解码性能损失。
具体来说,本发明提出的适用于低密度奇偶校验码的混合译码方法中各主要参数定义如下:
LDPC码使用一个M×N的校验矩阵H来进行定义,M为矩阵行数,N为矩阵列数;矩阵每行中非零元素的个数dc代表行重即校验节点的度数,矩阵中每列非零元素的个数ds代表列重即变量节点的度数;vn代表第n个变量节点,cm代表第m个校验节点。W=[w1,w2,…,wN]为LDPC码码字;A(n)为变量节点vn参与的校验节点集合,B(m)校验节点cm包含的变量节点的集合,B(m)\n为除变量节点vn的B(m)集合;Lnm、Lmn分别为变量节点vn传递给校验节点cm的信息和校验节点cm传递给变量节点vn的信息。采用BPSK调制,X=[x1,x2,…,xN]其中xn=1-2wn,n∈[1,N],通过AWGN信道;接收序列Y=[y1,y2,…,yN]其中yn=xn+zn,zn为均值为0,方差为σ2=N0/2的独立高斯随机变量,最大迭代次数为K,sign(x)表示对x取符号位,
Figure BDA0001485603320000051
本发明方法如图1所示,包括如下步骤:
S101:迭代初始化:初始化迭代次数k=0;对每个变量节点vn,n∈[1,N]初始化:
Figure BDA0001485603320000052
其中β为常数固定乘性修正因子;完成初始化后,顺序执行步骤S102;
S102:k=k+1,即迭代次数加一,顺序执行步骤S103;
S103:更新校验节点:对于各个校验节点cm(m∈[1,M]),利用第k-1次迭代中产生的由相应变量节点传递而来的信息
Figure BDA0001485603320000061
更新该校验节点向相应变量节点vn(n∈B(m))传递的信息
Figure BDA0001485603320000062
Figure BDA0001485603320000063
Figure BDA0001485603320000064
其中T为设定的阈值;根据第k-1次迭代中产生的由相应变量节点传递而来的信息
Figure BDA0001485603320000065
找到
Figure BDA0001485603320000066
绝对值的最小值的下标索引jmin,初始化局部变量δ0=∞;对每个h∈{1...dc}\jmin,计算
Figure BDA0001485603320000067
Figure BDA0001485603320000068
Figure BDA0001485603320000069
其中,
Figure BDA00014856033200000610
Figure BDA00014856033200000611
基于归一化修正最小和方法,
Figure BDA00014856033200000612
α为修正因子;完成后,顺序执行S104;
S104:更新输出似然比信息:对于各个变量节点vn(n∈[1,N]),根据由相应的校验节点传递而来的信息
Figure BDA00014856033200000613
和该变量节点相对应的初始输入似然比信息
Figure BDA00014856033200000614
计算该变量节点vn(n∈[1,N])的输出似然比信息
Figure BDA00014856033200000615
Figure BDA00014856033200000616
完成计算后,顺序执行步骤S105;
S105:变量节点硬判决:对每个变量节点vn(n∈[1,N])进行硬判决并计算校正子:
Figure BDA00014856033200000617
Sk=WkHT
完成上述计算后,顺序执行步骤S106;
S106:若
Figure BDA0001485603320000071
或k=K,则结束译码,输出码字
Figure BDA0001485603320000072
否则顺序执行步骤S107;
S107:更新变量节点:对于各个变量节点vn(n∈[1,N]),利用由相应校验节点cm(m∈A(n))传递而来的信息
Figure BDA0001485603320000073
和步骤S104中计算出来的该变量节点vn(n∈[1,N])的输出似然比信息
Figure BDA0001485603320000074
计算该变量节点向相连的各校验节点cm(m∈A(n))传递的信息
Figure BDA0001485603320000075
Figure BDA0001485603320000076
完成计算后,顺序执行步骤S102。
为了进一步说明本发明效果,我们提供了实施例1和实施例2。
实施例1:
图2为一个1/4码率的LDPC码(Code1-R14-N8448)校验矩阵的分块示意图,其校验矩阵N=8448,M=6400,信息比特K=2048,编码比特打孔前256比特;该编码的分块矩阵形式具体描述如下,C矩阵为3行47列,每个元素值都为-1;D矩阵为47行47列,除对角元素为0外,其他元素都为-1;矩阵A和B为25行19列,具体元素值见表1表2。
表1:Code1-R14-N8448码分块式校验矩阵A分块(25*19)
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1
14 83 111 68 66 30 28 29 9 42 20 11 55 22 12 117 1 0 0
43 -1 101 62 121 108 103 114 91 70 33 116 3 120 61 72 0 -1 0
56 65 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
54 56 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 23 105 96 -1 -1
0 -1 -1 -1 -1 -1 -1 -1 65 119 8 38 48 68 -1 -1 -1 -1 -1
121 36 -1 -1 -1 86 107 -1 -1 -1 -1 -1 -1 -1 0 81 88 -1 -1
71 117 -1 -1 -1 -1 -1 -1 4 0 72 108 79 -1 -1 -1 -1 -1 -1
79 -1 -1 -1 -1 10 4 -1 -1 -1 -1 -1 -1 59 117 31 64 -1 -1
121 40 -1 2 -1 -1 -1 -1 44 17 -1 3 -1 -1 -1 -1 -1 -1 -1
121 -1 -1 -1 -1 -1 -1 67 -1 -1 29 -1 -1 -1 124 72 22 -1 -1
-1 38 -1 -1 -1 -1 -1 -1 11 114 -1 73 -1 15 -1 -1 -1 -1 -1
24 -1 99 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 12 96 26 -1 -1
36 48 -1 -1 14 -1 -1 -1 27 101 -1 -1 -1 -1 -1 -1 -1 -1 -1
81 -1 -1 -1 -1 -1 -1 70 -1 -1 -1 101 -1 -1 80 -1 31 -1 -1
-1 127 -1 -1 40 -1 -1 -1 -1 115 41 -1 -1 -1 -1 -1 -1 -1 -1
95 -1 116 102 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 75 -1 81 -1 -1
117 -1 -1 -1 -1 -1 120 -1 -1 25 -1 70 -1 -1 -1 -1 -1 -1 -1
12 -1 -1 8 5 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 38 -1 -1
79 -1 -1 -1 -1 105 -1 18 -1 -1 -1 -1 -1 -1 24 -1 -1 -1 -1
91 -1 33 -1 -1 -1 -1 -1 -1 -1 -1 -1 112 -1 -1 -1 21 -1 -1
38 -1 -1 -1 123 -1 80 -1 -1 -1 -1 -1 -1 -1 2 -1 -1 -1 -1
12 96 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 99 -1 12
-1 75 -1 -1 -1 -1 -1 -1 -1 73 -1 -1 -1 64 -1 -1 -1 -1 -1
50 41 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 60 -1 88 -1 42
表2:Code1-R14-N8448码分块式校验矩阵B分块(25*19)
14 -1 -1 -1 -1 -1 -1 -1 -1 84 -1 8 70 -1 -1 -1 -1 -1 -1
109 -1 -1 -1 -1 -1 -1 41 -1 -1 -1 -1 -1 -1 -1 -1 72 -1 14
119 -1 -1 -1 -1 22 -1 -1 16 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
-1 32 41 4 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
101 -1 -1 -1 -1 -1 60 -1 -1 -1 103 -1 -1 -1 -1 -1 -1 -1 -1
115 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 82 -1 -1 -1 49 -1 76
58 0 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 66 -1 -1 7 -1
-1 115 -1 -1 -1 -1 -1 -1 84 -1 -1 -1 -1 -1 -1 -1 95 -1 -1
60 -1 -1 -1 11 -1 -1 -1 -1 -1 -1 -1 -1 112 -1 -1 -1 -1 -1
30 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 35 -1 64 10 -1
99 -1 -1 -1 -1 -1 -1 -1 -1 -1 94 11 -1 -1 -1 -1 -1 -1 -1
79 -1 6 -1 -1 88 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 103 -1 -1
87 -1 -1 44 -1 -1 -1 123 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
81 -1 -1 -1 -1 -1 114 -1 -1 -1 -1 -1 71 -1 -1 -1 -1 -1 -1
-1 109 -1 -1 -1 -1 -1 -1 -1 16 -1 -1 -1 -1 -1 -1 -1 99 -1
54 -1 -1 -1 -1 60 -1 -1 -1 -1 -1 -1 -1 31 -1 -1 -1 -1 -1
91 -1 -1 -1 -1 -1 -1 -1 -1 -1 32 -1 -1 -1 -1 6 -1 -1 -1
22 -1 -1 -1 97 -1 -1 119 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
-1 89 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 18 -1 -1 -1 -1 4 -1
121 -1 77 33 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
-1 101 -1 -1 -1 -1 19 -1 -1 -1 -1 -1 -1 64 -1 -1 -1 -1 -1
81 -1 -1 -1 -1 120 -1 -1 -1 -1 -1 127 -1 -1 -1 -1 -1 -1 -1
-1 60 -1 -1 -1 -1 -1 -1 -1 -1 88 -1 -1 -1 -1 -1 -1 101 -1
111 -1 -1 -1 -1 -1 -1 -1 29 -1 -1 -1 -1 -1 -1 -1 -1 -1 69
67 -1 -1 -1 0 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 34 -1 -1 -1
Code1-R14-N8448码分块式校验矩阵中的每一个元素x代表一个小方阵:若该元素为-1,则为全零阵;否则为偏移量等于x的置换单位阵;对于该LDPC码,小方阵的维数为Z=128,仿真时前2Z比特打孔传输。表3中各个元素表示校验矩阵中不同行重的分布。
表3 Code1-R14-N8448码分块式校验矩阵行重分布
行重 19 18 17 8 7 6 5 4 3
Code1-R14-N8448码 行数 384 384 384 1536 768 2688 4992 7680 384
图3是针对Code1-R14-N8448码分别采用BP译码方法、NMS译码方法以及本发明所述的一种适用于低密度奇偶校验码的混合译码方法进行仿真的误帧率性能曲线。其中,译码器的最大迭代次数均设置为50次;NMS译码方法的修正因子为α=0.7;本发明所述的一种适用于低密度奇偶校验码的混合译码方法阈值T=10,初始似然比固定乘性修正因子β=1.047,NMS译码方法的修正因子为α=0.7;采用AWGN信道进行仿真;前256比特打孔传输。
从图3中可以看出对于Code1-R14-N8448码BP译码方法与NMS译码方法达到10-2误帧率有0.5dB左右性能损失。而本发明所述的一种低密度奇偶校验码译码的译码方法与BP译码方法性能损失在0.2dB以内。
实施例2:
图4为一个1/5码率的LDPC码(Code2-R15-N10496)校验矩阵的分块示意图,其校验矩阵N=10496,M=8448,信息比特K=2048,编码比特打孔前256比特;该编码的分块矩阵形式具体描述如下,C矩阵为3行63列,每个元素值都为-1;D矩阵为63行63列,除对角元素为0外,其他元素都为-1;矩阵A和B为33行19列,具体元素值见表4表5。
表4:Code2-R15-N10496码分块式校验矩阵A分块(33*19)
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1
14 83 111 68 66 30 28 29 9 42 20 11 55 22 12 117 1 0 0
43 -1 101 62 121 108 103 114 91 70 33 116 3 120 61 72 0 -1 0
56 65 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
54 56 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 23 105 96 -1 -1
0 -1 -1 -1 -1 -1 -1 -1 65 119 8 38 48 68 -1 -1 -1 -1 -1
121 36 -1 -1 -1 86 107 -1 -1 -1 -1 -1 -1 -1 0 81 88 -1 -1
71 117 -1 -1 -1 -1 -1 -1 4 0 72 108 79 -1 -1 -1 -1 -1 -1
79 -1 -1 -1 -1 10 4 -1 -1 -1 -1 -1 -1 59 117 31 64 -1 -1
121 40 -1 2 -1 -1 -1 -1 44 17 -1 3 -1 -1 -1 -1 -1 -1 -1
121 -1 -1 -1 -1 -1 -1 67 -1 -1 29 -1 -1 -1 124 72 22 -1 -1
-1 38 -1 -1 -1 -1 -1 -1 11 114 -1 73 -1 15 -1 -1 -1 -1 -1
24 -1 99 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 12 96 26 -1 -1
36 48 -1 -1 14 -1 -1 -1 27 101 -1 -1 -1 -1 -1 -1 -1 -1 -1
81 -1 -1 -1 -1 -1 -1 70 -1 -1 -1 101 -1 -1 80 -1 31 -1 -1
-1 127 -1 -1 40 -1 -1 -1 -1 115 41 -1 -1 -1 -1 -1 -1 -1 -1
95 -1 116 102 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 75 -1 81 -1 -1
117 -1 -1 -1 -1 -1 120 -1 -1 25 -1 70 -1 -1 -1 -1 -1 -1 -1
12 -1 -1 8 5 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 38 -1 -1
79 -1 -1 -1 -1 105 -1 18 -1 -1 -1 -1 -1 -1 24 -1 -1 -1 -1
91 -1 33 -1 -1 -1 -1 -1 -1 -1 -1 -1 112 -1 -1 -1 21 -1 -1
38 -1 -1 -1 123 -1 80 -1 -1 -1 -1 -1 -1 -1 2 -1 -1 -1 -1
12 96 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 99 -1 12
-1 75 -1 -1 -1 -1 -1 -1 -1 73 -1 -1 -1 64 -1 -1 -1 -1 -1
50 41 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 60 -1 88 -1 42
14 -1 -1 -1 -1 -1 -1 -1 -1 84 -1 8 70 -1 -1 -1 -1 -1 -1
109 -1 -1 -1 -1 -1 -1 41 -1 -1 -1 -1 -1 -1 -1 -1 72 -1 14
119 -1 -1 -1 -1 22 -1 -1 16 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
-1 32 41 4 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
101 -1 -1 -1 -1 -1 60 -1 -1 -1 103 -1 -1 -1 -1 -1 -1 -1 -1
115 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 82 -1 -1 -1 49 -1 76
58 0 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 66 -1 -1 7 -1
-1 115 -1 -1 -1 -1 -1 -1 84 -1 -1 -1 -1 -1 -1 -1 95 -1 -1
表5:Code2-R15-N10496码分块式校验矩阵B分块(33*19)
60 -1 -1 -1 11 -1 -1 -1 -1 -1 -1 -1 -1 112 -1 -1 -1 -1 -1
30 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 35 -1 64 10 -1
99 -1 -1 -1 -1 -1 -1 -1 -1 -1 94 11 -1 -1 -1 -1 -1 -1 -1
79 -1 6 -1 -1 88 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 103 -1 -1
87 -1 -1 44 -1 -1 -1 123 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
81 -1 -1 -1 -1 -1 114 -1 -1 -1 -1 -1 71 -1 -1 -1 -1 -1 -1
-1 109 -1 -1 -1 -1 -1 -1 -1 16 -1 -1 -1 -1 -1 -1 -1 99 -1
54 -1 -1 -1 -1 60 -1 -1 -1 -1 -1 -1 -1 31 -1 -1 -1 -1 -1
91 -1 -1 -1 -1 -1 -1 -1 -1 -1 32 -1 -1 -1 -1 6 -1 -1 -1
22 -1 -1 -1 97 -1 -1 119 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
-1 89 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 18 -1 -1 -1 -1 4 -1
121 -1 77 33 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
-1 101 -1 -1 -1 -1 19 -1 -1 -1 -1 -1 -1 64 -1 -1 -1 -1 -1
81 -1 -1 -1 -1 120 -1 -1 -1 -1 -1 127 -1 -1 -1 -1 -1 -1 -1
-1 60 -1 -1 -1 -1 -1 -1 -1 -1 88 -1 -1 -1 -1 -1 -1 101 -1
111 -1 -1 -1 -1 -1 -1 -1 29 -1 -1 -1 -1 -1 -1 -1 -1 -1 69
67 -1 -1 -1 0 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 34 -1 -1 -1
-1 123 -1 10 -1 -1 -1 -1 -1 -1 -1 -1 64 -1 -1 -1 -1 -1 -1
40 -1 0 -1 -1 -1 -1 51 -1 -1 -1 -1 -1 50 -1 -1 -1 -1 -1
85 -1 -1 -1 -1 -1 0 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 3 -1
103 -1 -1 -1 -1 -1 -1 45 108 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
36 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 3 -1 16 -1 -1 -1
56 -1 -1 -1 -1 84 -1 -1 -1 -1 31 -1 -1 -1 -1 -1 -1 -1 39
119 -1 80 -1 -1 -1 -1 -1 -1 -1 -1 -1 121 -1 -1 -1 -1 -1 -1
-1 -1 -1 26 68 -1 -1 -1 -1 -1 59 -1 -1 -1 -1 -1 -1 -1 -1
97 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 88 -1 -1 -1 -1 -1 82 -1
6 -1 -1 -1 -1 104 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 41
-1 -1 -1 -1 32 -1 43 -1 -1 -1 -1 -1 -1 9 -1 -1 -1 -1 -1
99 -1 -1 29 -1 -1 -1 -1 -1 2 -1 -1 -1 -1 -1 -1 -1 -1 -1
-1 -1 36 -1 -1 -1 -1 24 -1 -1 -1 -1 -1 -1 -1 51 -1 -1 -1
-1 81 -1 -1 -1 -1 -1 -1 83 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
-1 -1 -1 -1 -1 -1 40 -1 -1 -1 105 -1 102 -1 -1 -1 -1 -1 -1
-1 -1 -1 -1 2 83 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 44 -1
Code2-R15-N10496码分块式校验矩阵中的每一个元素x代表一个小方阵:若该元素为-1,则为全零阵;否则为偏移量等于x的置换单位阵;对于该LDPC码,小方阵的维数为Z=128,仿真时前2Z比特打孔传输。表6中各个元素表示校验矩阵中不同行重的分布。
表6 Code2-R15-N10496码分块式校验矩阵行重分布
行重 19 18 17 8 7 6 5 4 3
Code2-R15-N10496码 行数 384 384 384 1536 768 2688 5760 12672 768
图5是针对Code2-R15-N10496码分别采用BP译码方法、NMS译码方法以及本发明所述的一种适用于低密度奇偶校验码的混合译码方法进行仿真的误帧率性能曲线。其中,译码器的最大迭代次数均设置为50次;NMS译码方法的修正因子为α=0.7;本发明所述的一种适用于低密度奇偶校验码的混合译码方法阈值T=10,初始似然比固定乘性修正因子β=0.8,NMS译码方法的修正因子为α=0.7;采用AWGN信道进行仿真;前256比特打孔传输。
从图5中可以看出对于Code2-R15-N10496码BP译码方法与NMS译码方法存在0.6dB左右性能损失,本发明所述的一种低密度奇偶校验码译码的译码方法与BP译码方法损失缩小到0.2dB以内。
显然,本译码方法缩小了最小和译码方法与BP译码方法的解码性能差距。
本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段,还包括由以上技术特征任意组合所组成的技术方案。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。

Claims (7)

1.一种适用于低密度奇偶校验码的混合译码方法,其特征在于,包括如下步骤:
步骤一:对每个变量节点进行迭代初始化,并在变量节点中引入常数固定乘性修正因子;
步骤二:迭代次数加一;
步骤三:对于各个校验节点,利用第k-1次迭代中产生的由相应变量节点传递而来的信息更新该校验节点向相应变量节点传递的信息:校验节点度数小于阈值的节点采用基于线性插值的APPROXIMATE-MIN*方法进行更新,校验节点度数大于阈值的节点采用归一化修正最小和译码方法进行计算更新;
步骤四:对于各个变量节点,根据由相应的校验节点传递而来的信息和该变量节点相对应的初始输入似然比信息,计算该变量节点的输出似然比信息;
步骤五:对每个变量节点进行硬判决并计算校正子;
步骤六:如果满足条件则结束译码,输出码字;否则顺序执行下一步骤;
步骤七:对于各个变量节点,利用由相应校验节点传递而来的信息和步骤四中计算出来的该变量节点的输出似然比信息,计算该变量节点向相连的各校验节点传递的信息;
顺序执行步骤二,直至迭代完毕。
2.根据权利要求1所述的适用于低密度奇偶校验码的混合译码方法,其特征在于,所述步骤一具体包括:
初始化迭代次数k=0;对每个变量节点vn,n∈[1,N]初始化:
Figure FDA0002956403000000011
其中N为矩阵列数,β为常数固定乘性修正因子;Lnm为变量节点vn传递给校验节点cm的信息;yn为接收序列中的元素。
3.根据权利要求2所述的适用于低密度奇偶校验码的混合译码方法,其特征在于,所述步骤三具体包括:
对于各个校验节点cm,m∈[1,M],M为矩阵行数,利用第k-1次迭代中产生的由相应变量节点传递而来的信息
Figure FDA0002956403000000012
更新该校验节点向相应变量节点vn传递的信息
Figure FDA0002956403000000013
Figure FDA0002956403000000014
Figure FDA0002956403000000015
其中n∈B(m),B(m)为校验节点cm包含的变量节点的集合,T为设定的阈值;dc为矩阵每行中非零元素的个数,Lmn为校验节点cm传递给变量节点vn的信息;
根据第k-1次迭代中产生的由相应变量节点传递而来的信息
Figure FDA0002956403000000021
找到
Figure FDA0002956403000000022
绝对值的最小值的下标索引jmin,初始化局部变量δ0=∞;对每个h∈{1...dc}\jmin,计算
Figure FDA0002956403000000023
Figure FDA00029564030000000213
Figure FDA0002956403000000024
其中,B(m)\n为除变量节点vn的B(m)集合,
Figure FDA0002956403000000025
Figure FDA0002956403000000026
归一化修正最小和方法,
Figure FDA0002956403000000027
α为修正因子。
4.根据权利要求3所述的适用于低密度奇偶校验码的混合译码方法,其特征在于,所述步骤四具体包括:
对于各个变量节点vn,根据由相应的校验节点传递而来的信息
Figure FDA0002956403000000028
和该变量节点相对应的初始输入似然比信息
Figure FDA0002956403000000029
计算该变量节点vn的输出似然比信息
Figure FDA00029564030000000210
Figure FDA00029564030000000211
其中,m∈A(n),n∈[1,N],A(n)为变量节点vn参与的校验节点集合。
5.根据权利要求4所述的适用于低密度奇偶校验码的混合译码方法,其特征在于,所述步骤五具体包括:
对每个变量节点vn进行硬判决并计算校正子:
Figure FDA00029564030000000212
Sk=WkHT
6.根据权利要求5所述的适用于低密度奇偶校验码的混合译码方法,其特征在于,所述步骤六中条件为
Figure FDA0002956403000000031
或k=K,输出的码字
Figure FDA0002956403000000032
7.根据权利要求6所述的适用于低密度奇偶校验码的混合译码方法,其特征在于,所述步骤七中计算该变量节点向相连的各校验节点传递的信息的过程具体包括:
对于各个变量节点vn,利用由相应校验节点cm传递而来的信息
Figure FDA0002956403000000033
和步骤四中计算出来的该变量节点vn的输出似然比信息
Figure FDA0002956403000000034
m∈A(n),n∈[1,N],计算该变量节点向相连的各校验节点cm传递的信息
Figure FDA0002956403000000035
Figure FDA0002956403000000036
CN201711216506.6A 2017-11-28 2017-11-28 一种适用于低密度奇偶校验码的混合译码方法 Active CN107968657B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711216506.6A CN107968657B (zh) 2017-11-28 2017-11-28 一种适用于低密度奇偶校验码的混合译码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711216506.6A CN107968657B (zh) 2017-11-28 2017-11-28 一种适用于低密度奇偶校验码的混合译码方法

Publications (2)

Publication Number Publication Date
CN107968657A CN107968657A (zh) 2018-04-27
CN107968657B true CN107968657B (zh) 2021-05-18

Family

ID=61998850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711216506.6A Active CN107968657B (zh) 2017-11-28 2017-11-28 一种适用于低密度奇偶校验码的混合译码方法

Country Status (1)

Country Link
CN (1) CN107968657B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102475805B1 (ko) * 2018-08-20 2022-12-09 에스케이하이닉스 주식회사 디코더, 그것의 동작 방법 및 메모리 시스템
CN109541444B (zh) * 2018-10-18 2021-11-02 天津大学 基于混合粒度奇偶校验的集成电路故障注入检测方法
CN109889208A (zh) * 2019-01-04 2019-06-14 浙江大学 Navdat中基于归一化最小和算法的ldpc信道译码方法
CN111313911A (zh) * 2019-11-22 2020-06-19 南京大学 一种适用于5g ldpc码的度自适应定点化译码算法
CN111641416B (zh) * 2020-06-19 2023-04-07 重庆邮电大学 一种多归一化因子的低密度奇偶校验码译码方法
CN111865334A (zh) * 2020-08-17 2020-10-30 长安大学 一种低密度奇偶校验码的高速滑窗置信解码方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101107782A (zh) * 2005-05-20 2008-01-16 三菱电机株式会社 用于解码纠错码的方法
CN101277291A (zh) * 2006-10-17 2008-10-01 北京凌讯华业科技有限公司 低比特分辨率的混合最小和ldpc解码方法和装置
CN101465652A (zh) * 2007-12-20 2009-06-24 中兴通讯股份有限公司 一种低密度奇偶校验的译码方法
CN103199874A (zh) * 2012-01-05 2013-07-10 国民技术股份有限公司 一种低密度奇偶校验码译码方法
CN103607208A (zh) * 2013-11-25 2014-02-26 上海数字电视国家工程研究中心有限公司 基于归一化修正因子序列的ldpc最小和译码方法
CN106452455A (zh) * 2016-09-23 2017-02-22 华南理工大学 基于OpenCL 移动设备QC‑LDPC 的动态译码方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080109698A1 (en) * 2006-07-25 2008-05-08 Legend Silicon Hybrid min-sum decoding apparatus with low bit resolution for ldpc code

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101107782A (zh) * 2005-05-20 2008-01-16 三菱电机株式会社 用于解码纠错码的方法
CN101277291A (zh) * 2006-10-17 2008-10-01 北京凌讯华业科技有限公司 低比特分辨率的混合最小和ldpc解码方法和装置
CN101465652A (zh) * 2007-12-20 2009-06-24 中兴通讯股份有限公司 一种低密度奇偶校验的译码方法
CN103199874A (zh) * 2012-01-05 2013-07-10 国民技术股份有限公司 一种低密度奇偶校验码译码方法
CN103607208A (zh) * 2013-11-25 2014-02-26 上海数字电视国家工程研究中心有限公司 基于归一化修正因子序列的ldpc最小和译码方法
CN106452455A (zh) * 2016-09-23 2017-02-22 华南理工大学 基于OpenCL 移动设备QC‑LDPC 的动态译码方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"APPROXIMATE-MIN CONSTRAINT NODE UPDATING FOR LDPC CODE DECODING";Christopher Jones等;《IEEE Military Communications Conference》;20031016;157-162 *
"Near Optimum Universal Belief Propagation Based Decoding of Low-Density Parity Check Codes";Jinghu Chen等;《IEEE TRANSACTIONS ON COMMUNICATIONS》;20020331;第50卷(第3期);406-414 *
"Nonbinary Hybrid LDPC Codes";Lucile Sassatelli等;《IEEE TRANSACTIONS ON INFORMATION THEORY》;20101031;第56卷(第10期);5314-5334 *
"一种基于改进线性规划的LDPC码混合译码算法";陈紫强等;《电路与系统学报》;20150228;第18卷(第1期);107-112 *

Also Published As

Publication number Publication date
CN107968657A (zh) 2018-04-27

Similar Documents

Publication Publication Date Title
CN107968657B (zh) 一种适用于低密度奇偶校验码的混合译码方法
EP2093887B1 (en) Apparatus and method for channel encoding and decoding in a communication system using low-density parity-check codes
CN1830149B (zh) 用于对低密度奇偶校验码编码的方法
US8321745B2 (en) Apparatus and method for channel encoding and decoding in communication system using low-density parity-check codes
JP5506878B2 (ja) 低密度パリティ検査符号のパリティ検査行列生成方法
WO2017080249A1 (zh) 生成用于在信道中传输的低密度奇偶校验码的方法及设备
CN109286405B (zh) 一种低复杂度的极化码递进式比特翻转sc译码方法
JP5461728B2 (ja) 低密度パリティ検査符号を使用するデジタルブロードキャスト通信システムにおけるチャネル復号化装置及びその方法
JP2008541496A (ja) 誤り訂正符号を復号化する方法
EP3649737B1 (en) Generalized low-density parity check codes (gldpc)
CN107404321B (zh) 用于纠错码解码的方法和设备
CN110830050B (zh) 一种ldpc译码方法、系统、电子设备及存储介质
EP2472725A1 (en) Method and apparatus for encoding and decoding a channel in a communication system using a linear block code
CN104158550A (zh) 一种基于深空通信环境的码率兼容原模图ldpc码构造方法
Sridharan et al. Convergence analysis for a class of LDPC convolutional codes on the erasure channel
CN111034055A (zh) 在非二进制ldpc解码器中简化的校验节点处理
CN107615666A (zh) Ldpc截短码的译码方法和译码设备
CN101867449B (zh) 基于地面数字电视的高效ldpc译码器
CN108566211B (zh) 基于H矩阵层处理顺序动态变化的layered LDPC译码方法
KR20090012189A (ko) Ldpc 부호의 성능 개선을 위한 스케일링 기반의 개선된min-sum 반복복호알고리즘을 이용한 복호 장치 및그 방법
CN113055024B (zh) 用于5g-nr系统短块长低码率ldpc码的修正译码方法
CN109639394B (zh) 多边类型低密度奇偶校验码的分边类中继译码方法
Kumar et al. A review of channel coding schemes in the 5G standard
CN113364470B (zh) 用于qkd的基于ge模型的改进型ldpc码译码方法
CN110708139B (zh) 一种基于并行矢量消息传递算法的ldpc码字扩展方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant