CN107689801A - Ldpc码admm迭代译码的早停止方法 - Google Patents
Ldpc码admm迭代译码的早停止方法 Download PDFInfo
- Publication number
- CN107689801A CN107689801A CN201710798802.5A CN201710798802A CN107689801A CN 107689801 A CN107689801 A CN 107689801A CN 201710798802 A CN201710798802 A CN 201710798802A CN 107689801 A CN107689801 A CN 107689801A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- code word
- admm
- represent
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M13/00—Coding, 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/03—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
- H03M13/05—Error 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/11—Error 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/1102—Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
- H03M13/1105—Decoding
- H03M13/1108—Hard decision decoding, e.g. bit flipping, modified or weighted bit flipping
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M13/00—Coding, 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/03—Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
- H03M13/05—Error 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/11—Error 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/1102—Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
- H03M13/1148—Structural properties of the code parity-check or generator matrix
Landscapes
- Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Error Detection And Correction (AREA)
Abstract
一种低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，主要解决现有技术的中低信噪比区域ADMM迭代译码收敛速度慢的问题。本发明实现方法的步骤：初始化译码参数；计算变量节点初始信息；更新变量节点信息；更新辅助变量；更新拉格朗日乘子向量；获取码字比特的硬判决值；计算ADMM迭代译码中每两次连续迭代之间发生变化的码字比特总数；计算码字比特的硬判决变化率；译码停止判决；停止译码。本发明根据ADMM迭代译码过程中每两次连续迭代之间码字比特的硬判决变化率CRHD的大小来判断是否能够及早停止迭代，从而减少译码的平均迭代次数，提高译码速度。
Description
技术领域
本发明属于通信技术领域，更进一步涉及信道编码技术领域中低密度奇偶检验码(Low-Density Parity-Check Codes，LDPC码)交替方向乘子法(Alternating DirectionMethod of Multipliers，ADMM)迭代译码的一种有效的早停止(Early Termination，ET)方法。本发明可实现及早停止LDPC码的ADMM迭代译码。
背景技术
具有低译码复杂度和逼近香农限良好性能的LDPC码已经被广泛应用于现代通信的深空通信、无线通信等领域中，并被802.11n、802.16e、10GBASE-T等各种现代通信标准采纳。因此，LDPC码及其译码方法已经成为近年来信道编码领域普遍关注的研究热点。
基于ADMM的LDPC码线性规划译码方法具有最大似然认证特性，但是收敛速度较慢，译码复杂度较高。利用早停止方法可以在ADMM译码过程的早期阶段结束译码，有效降低译码延迟并且提高收敛速度。而目前LDPC码ADMM译码只有两种迭代停止方法。
西安电子科技大学在其技术的专利文献“基于加速交替方向乘子法的LDPC码线性规划译码方法”(申请公布日：2014年10月8日，申请公布号：CN 104092468A，申请号：2014103209428)中公开了一种基于加速交替方向乘子法的LDPC码线性规划译码方法。该专利申请中的译码方法通过引入加速模块和改变增广拉格朗日分解式的更新次序，减少译码迭代次数，从而提高译码速度。该方法可以在高信噪比区域降低ADMM译码的平均迭代次数，其不足之处是在中低信噪比区域的平均迭代次数仍然较高。
Barman等人在其发表的论文“Decomposition methods for large scale LPdecoding”(IEEE Trans.Inf.Theory，2013，59，(12)，pp.7870-7886.)中提出了一种基于大规模低密度奇偶检验LDPC码线性规划译码的分解方法。该方法通过重新设计适用于LDPC码译码的线性规划问题求解方法，针对ADMM迭代译码的标准ε规则，将消息主残差和对偶残差分别与ε比较，若两个值都比ε小，则迭代停止，从而使译码效率大大提高。该方法存在的不足之处是：LDPC码ADMM译码迭代过程中更新校验节点信息需要执行耗时的欧几里德投影运算，译码收敛速度慢。
发明内容
本发明的目的在于针对上述已有技术的不足，提出LDPC码ADMM迭代译码的一种有效的早停止方法，可实现及早停止LDPC码的ADMM迭代译码，降低译码的平均迭代次数，提高译码器的收敛速度。
为了实现上述目的，本发明方法的思路是：通过跟踪ADMM迭代译码过程中每两个连续迭代之间码字比特的硬判决变化率(Change Rate of Hard Decision，CRHD)，并根据CRHD值的大小来判断是否及早停止迭代，从而降低译码的平均迭代次数，提高译码器的收敛速度。
本发明方法的实现包括如下步骤：
(1)初始化译码参数：
(1a)计算低密度奇偶检验LDPC码的码字比特对应的代价参数；
(1b)将交替方向乘子法ADMM译码中的容差值设置为10-5、超松弛参数设置为1.9、初始迭代次数设置为0；
(1c)在[1，20]范围内，设置交替方向乘子法ADMM迭代译码停止时的迭代次数，在[0.05，0.10]范围内，设置交替方向乘子法ADMM迭代译码中码字比特的硬判决变化率CRHD的阈值；
(1d)将拉格朗日乘子向量λ初始值设置为全0向量、辅助向量z初始值设置为全0向量；
(2)计算每个变量节点的初始信息：
xq＝γi
其中，xq表示第q个变量节点的初始信息，q的取值范围为[1，N]，N表示低密度奇偶校验LDPC码的码字长度，γi表示初始化中计算的低密度奇偶检验LDPC码的码字比特i对应的代价参数，i的取值范围为[1，N]，q与i的取值对应相等；
(3)更新每个变量节点的当前信息；
(4)更新辅助变量；
(5)更新拉格朗日乘子向量；
(6)获取每个码字比特的硬判决值：
其中，k表示第k次迭代次数，表示第i个码字比特在第k次迭代结束后对应的硬判决HD值，xi表示第i个变量节点的当前信息；
(7)计算交替方向乘子法ADMM译码中每两次连续迭代之间发生变化的码字比特总数；
(8)计算码字比特的硬判决变化率；
(9)判断是否满足早停止译码的条件，若是，则执行步骤(10)，否则，将当前迭代次数加1后执行步骤(3)；
(10)停止交替方向乘子法ADMM译码。
本发明与现有技术相比具有以下优点：
第一，由于本发明使用交替方向乘子法ADMM迭代译码过程中码字比特的硬判决HD信息设计译码及早停止的条件，从而克服了现有技术的置信传播(Belief Propagation，BP)译码中使用信号变化率(Sign-Changing Rate，SCR)设计早停止方法时需要额外计算对数似然比(Log Likelihood Ratios，LLRs)的值，译码延迟较大和译码效率不高的缺点，从而使得本发明能够提高译码速度。
第二，由于本发明通过跟踪ADMM迭代译码过程中每两次连续迭代之间码字比特的硬判决变化率CRHD，并根据CRHD值的大小来判断是否能够及早停止迭代，从而避免了现有技术的ADMM译码迭代过程中需要执行大量欧几里德投影运算的缺点，使得本发明能够降低译码复杂度。
附图说明
图1是本发明的流程图；
图2是本发明取不同Thr值时与现有方法收敛速度和译码性能对比图；
图3是本发明取不同Ite值时与现有方法收敛速度和译码性能对比图；
图4是本发明取最佳参数组合时与现有方法收敛速度和译码性能对比图。
具体实施方式
下面结合附图对本发明做进一步描述。
下面结合附图1，对本发明的具体步骤做进一步描述。
步骤1，初始化译码参数。
计算低密度奇偶检验LDPC码的码字比特对应的代价参数：
其中，γi表示低密度奇偶检验LDPC码的第i个码字比特对应的代价参数，i的取值范围为[1，N]，N表示低密度奇偶校验LDPC码的码字长度，log表示以2为底的对数操作，Pr(·)表示事件发生概率，rw表示加性高斯白噪声信道输出的第w位信息，w的取值范围为[1，N]，Pr(A|B)表示当事件B发生时事件A发生的概率，ct表示发送的第t位信息，t的取值范围为[1，N]，i、w与t的取值对应相等。
将交替方向乘子法ADMM译码中的容差值设置为10-5、超松弛参数设置为1.9、初始迭代次数设置为0。
在[1，20]范围内，设置交替方向乘子法ADMM迭代译码停止时的迭代次数，在[0.05，0.10]范围内，设置交替方向乘子法ADMM迭代译码中码字比特的硬判决变化率CRHD的阈值。
将拉格朗日乘子向量λ初始值设置为全0向量、辅助向量z初始值设置为全0向量。
步骤2，计算每个变量节点的初始信息：
xq＝γi
其中，xq表示第q个变量节点的初始信息，q的取值范围为[1，N]，N表示低密度奇偶校验LDPC码的码字长度，γi表示初始化中计算的低密度奇偶检验LDPC码的码字比特i对应的代价参数，i的取值范围为[1，N]，q与i的取值对应相等；
步骤3，更新每个变量节点的当前信息：
其中，xq表示第q个变量节点的当前信息，Π[0,1]表示在区间[0，1]内做欧几里德投影操作，dw表示与第w个变量节点相邻的校验节点的总数，∑表示求和操作，j表示第j个校验节点，∈表示属于符号，Nv(w)表示与第w个变量节点相邻的校验节点的索引集，表示第j个校验节点对应的辅助向量zj中与第t个变量节点对应的向量位，μ表示针对具体信道、具体码字，根据误码率大小优化得到的增广拉格朗日参数，表示第j个校验节点对应的拉格朗日乘子向量λj中与第t个变量节点对应的向量位，α表示针对具体信道、具体码字，根据误码率大小优化得到的罚函数的惩罚参数，sgn(·)表示符号函数，xi表示第i个变量节点更新前的信息，q、i、w与t的取值对应相等。
步骤4，更新辅助变量：
其中，zj表示第j个校验节点对应的辅助向量，表示在第j个校验节点对应的校验多胞体Pj上做欧几里德投影操作，Tj表示第j个校验节点对应的转换矩阵，x表示变量节点消息向量，λj表示第j个校验节点对应的拉格朗日乘子向量。
步骤5，更新拉格朗日乘子向量：
λj＝λj+μ(Tjx-zj)
其中，λr表示第r个校验节点对应的当前拉格朗日乘子向量，r与j的取值对应相等。
步骤6，获取每个码字比特的硬判决值：
其中，k表示第k次迭代次数，表示第i个码字比特在第k次迭代结束后对应的硬判决HD值，xi表示第i个变量节点的当前信息。
步骤7，计算交替方向乘子法ADMM译码中每两次连续迭代之间发生变化的码字比特总数：
其中，Num表示交替方向乘子法ADMM译码中每两次连续迭代之间发生变化的码字比特总数，i表示低密度奇偶校验LDPC码中的第i个码字比特，|·|表示取绝对值操作，表示对应第i个码字比特在第k-1次迭代结束后对应的硬判决HD值。
步骤8，计算码字比特的硬判决变化率：
其中，CRHD表示交替方向乘子法ADMM译码中每两次连续迭代之间码字比特的硬判决变化率。
步骤9，判断是否满足早停止译码的条件，若是，则执行步骤10，否则，将当前迭代次数加1后执行步骤3。
所述的译码早停止条件是指满足以下三个条件之一的情形：
条件1：交替方向乘子法ADMM译码中，每两次连续迭代之间发生变化的码字比特总数为0。
条件2：不等式k≥Ite与CRHD≥Thr同时成立，其中，Ite表示初始化时所选取的交替方向乘子法ADMM迭代译码停止时的迭代次数，Thr表示初始化时所选取的交替方向乘子法ADMM迭代译码中码字比特的硬判决变化率CRHD的阈值。
条件3：交替方向乘子法ADMM迭代次数达到20。
步骤10，停止交替方向乘子法ADMM译码。
下面通过本发明与现有ADMM迭代译码的早停止方法比较，对本发明LDPC码ADMM迭代译码的早停止方法的效果作进一步描述。
在加性高斯白噪声信道下，利用ADMM惩罚译码对IEEE802.16E标准中码率为1/2长度为576的LDPC码C1和码率为3/4长度为1152的LDPC码C2进行仿真，译码分别采用现有的两种停止方法和本发明提出的早停止方法。译码过程中，使用Liu,X.-S.等人在“The ADMMpenalized decoder for LDPC codes”(IEEE Trans.Inf.Theory，2016，62，(6)，pp.2966-2984)中提到的l1罚函数，令α和μ分别为l1罚函数的惩罚参数和增广拉格朗日函数的惩罚参数。针对码C1，l1罚函数的惩罚参数设置为1.5，增广拉格朗日函数的惩罚参数设置为4.1。针对码C2，l1罚函数的惩罚参数设置为2.1，增广拉格朗日函数的惩罚参数设置为3.6。
现有的ADMM迭代译码停止方法有两种。一种是Barman等人在其发表的论文“Decomposition methods for large scale LP decoding”(IEEE Trans.Inf.Theory，2013，59，(12)，pp.7870-7886.)中使用的针对ADMM迭代译码的标准ε规则停止方法，该方法将消息主残差和对偶残差分别与ε比较，若两个值都比ε小，则迭代停止。另一种是Siegel等人在其发表的论文“Efficient iterative LP decoding of LDPC codes withalternating direction method of multipliers”(Proc.IEEE Int.Symp.Inf.Theory，Istanbul，Turkey,July 2013，pp.1501-1505.)中提出的早停止方法，该方法在每次迭代结束时验证校验方程HxT＝0是否成立，若方程成立，则迭代停止，其中H为LDPC码的校验矩阵，x为码字比特的硬判决HD。
图2为本发明取不同Thr值时与现有技术的ADMM迭代译码停止方法的收敛速度和译码性能对比图。图2(a)表示现有技术ADMM迭代译码停止方法与本发明取不同Thr值时译码的平均迭代次数ANIs(Average Number of Iterations)的对比图。图2(a)中的横坐标表示信噪比，纵坐标表示平均迭代次数ANIs。图2(a)中以圆圈标识的曲线表示本发明提出的ADMM迭代译码的早停止方法取Ite＝8和Thr＝5％时码C1的平均迭代次数ANIs。图2(a)中以叉号标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝6％时码C1的平均迭代次数ANIs。图2(a)中以加号标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝7％时码C1的平均迭代次数ANIs。图2(a)中以倒三角形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝8％时码C1的平均迭代次数ANIs。图2(a)中以米形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝9％时码C1的平均迭代次数ANIs。图2(a)中以五角星形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝10％时码C1的平均迭代次数ANIs。图2(a)中以菱形标识的曲线表示利用HxT＝0早停止方法进行ADMM迭代译码时码C1的平均迭代次数ANIs。图2(a)中以方框标识的曲线表示利用标准ε规则停止方法进行ADMM迭代译码时码C1的平均迭代次数ANIs。
图2(b)表示现有技术的ADMM迭代译码停止方法与本发明取不同Thr值时译码的误比特率BER(Bit Error Rate)性能对比图。图2(b)中的横坐标表示信噪比，纵坐标表示误比特率BER。图2(b)中以圆圈标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝5％时码C1的误比特率BER性能。图2(b)中以叉号标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝6％时码C1的误比特率BER性能。图2(b)中以加号标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝7％时码C1的误比特率BER性能。图2(b)中以倒三角形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝8％时码C1的误比特率BER性能。图2(b)中以米形标识的曲线表示取本发明提出的ADMM迭代译码早停止方法Ite＝8和Thr＝9％时码C1的误比特率BER性能。图2(b)中以五角星形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Ite＝8和Thr＝10％时码C1的误比特率BER性能。图2(b)中以菱形标识的曲线表示利用HxT＝0早停止方法进行ADMM迭代译码时码C1的误比特率BER性能。图2(b)中以方框标识的曲线表示利用标准ε规则停止方法进行ADMM迭代译码时码C1的误比特率BER性能。
由图2可以看出，与现有技术的两种ADMM迭代译码的停止方法比较，在性能损失不大的情况下，本发明提出的早停止方法在低信噪比区域可以显著地降低平均迭代次数ANIs。而且，平均迭代次数ANIs的减少在低Thr下更显著，然而低Thr带来了较多的性能损失。因此，有必要在Thr和误比特率BER性能之间进行权衡。
图3为本发明取不同Ite值时与现有技术的ADMM迭代译码停止方法的收敛速度和译码性能对比图。图3(a)表示现有ADMM迭代译码停止方法与本发明取不同Ite值时译码的平均迭代次数ANIs对比图。图3(a)中的横坐标表示信噪比，纵坐标表示平均迭代次数ANIs。图3(a)中以圆圈标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝6时码C1的平均迭代次数ANIs。图3(a)中以叉号标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝7时码C1的平均迭代次数ANIs。图3(a)中以加号标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝8时码C1的平均迭代次数ANIs。图3(a)中以倒三角形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝9时码C1的平均迭代次数ANIs。图3(a)中以米形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝10时码C1的平均迭代次数ANIs。图3(a)中以五角星形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝11时码C1的平均迭代次数ANIs。图3(a)中以菱形标识的曲线表示利用HxT＝0早停止方法进行ADMM迭代译码时码C1的平均迭代次数ANIs。图3(a)中以方框标识的曲线表示利用标准ε规则停止方法进行ADMM迭代译码时码C1的平均迭代次数ANIs。
图3(b)表示现有技术的ADMM迭代译码停止译码方法与本发明取不同Ite值时译码的误比特率BER性能对比图。图3(b)中的横坐标表示信噪比，纵坐标表示误比特率BER。图3(b)中以圆圈标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝6时码C1的误比特率BER性能。图3(b)中以叉号标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝7时码C1的误比特率BER性能。图3(b)中以加号标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝8时码C1的误比特率BER性能。图3(b)中以倒三角形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝9时码C1的误比特率BER性能。图3(b)中以米形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝10时码C1的误比特率BER性能。图3(b)中以五角星形标识的曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝11时码C1的误比特率BER性能。图3(b)中以菱形标识的曲线表示利用HxT＝0早停止方法进行ADMM迭代译码时码C1的误比特率BER性能。图3(b)中以方框标识的曲线表示利用标准ε规则停止方法进行ADMM迭代译码时码C1的误比特率BER性能。
由图3可以看出，与现有技术的两种ADMM迭代译码的停止方法比较，在性能损失不大的情况下，本发明提出的早停止方法在低信噪比区域可以显著地降低平均迭代次数ANIs。而且，平均迭代次数ANIs的减少在低Ite下更显著，然而低Ite带来了较多的性能损失。因此，有必要在Ite和误比特率BER性能之间进行权衡。
通过本发明的仿真，得到码C1的最佳组合参数为(Ite，Thr)＝(10，7％)，同样地，得到码C2的最佳组合参数为(Ite，Thr)＝(11，7％)。
图4为本发明取最佳参数组合(Ite，Thr)时与现有ADMM迭代译码停止方法的收敛速度和译码性能对比图。图4(a)表示现有ADMM迭代译码停止方法与本发明取最佳参数组合(Ite，Thr)时译码的平均迭代次数ANIs对比图。图4(a)中的横坐标表示信噪比，纵坐标表示平均迭代次数ANIs。图4(a)中以米形标识的实曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝10时码C1的平均迭代次数ANIs。图4(a)中以菱形标识的实曲线表示利用HxT＝0早停止方法进行ADMM迭代译码时码C1的平均迭代次数ANIs。图4(a)中以方框标识的实曲线表示利用标准ε规则停止方法进行ADMM迭代译码时码C1的平均迭代次数ANIs。图4(a)中以米形标识的虚曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝11时码C2的平均迭代次数ANIs。图4(a)中以菱形标识的虚曲线表示利用HxT＝0早停止方法进行ADMM迭代译码时码C2的平均迭代次数ANIs。图4(a)中以方框标识的虚曲线表示利用标准ε规则停止方法进行ADMM迭代译码时码C2的平均迭代次数ANIs。
图4(b)表示现有ADMM迭代译码停止方法与本发明取最佳参数组合(Ite，Thr)时译码的误比特率BER性能对比图。图4(b)中的横坐标表示信噪比，纵坐标表示误比特率BER。图4(b)中以米形标识的实曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝10时码C1的误比特率BER性能。图4(b)中以菱形标识的实曲线表示利用HxT＝0早停止方法进行ADMM迭代译码时码C1的误比特率BER性能。图4(b)中以方框标识的实曲线表示利用标准ε规则停止方法进行ADMM迭代译码时码C1的误比特率BER性能。图4(b)中以米形标识的虚曲线表示本发明提出的ADMM迭代译码早停止方法取Thr＝7％和Ite＝11时码C2的误比特率BER性能。图4(b)中以菱形标识的虚曲线表示利用HxT＝0早停止方法进行ADMM迭代译码时码C2的误比特率BER性能。图4(b)中以方框标识的虚曲线表示利用标准ε规则停止方法进行ADMM迭代译码时码C2的误比特率BER性能。
由图4可以看出，与现有技术的两种停止方法相比，本发明所提出的ADMM迭代译码早停止方法在低信噪比区域可显著降低译码的平均迭代次数ANIs，且误比特率BER性能损失可以忽略不计。
Claims (8)
1.一种低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，其特征在于，包括如下步骤：
(1)初始化译码参数：
(1a)计算低密度奇偶检验LDPC码的码字比特对应的代价参数；
(1b)将交替方向乘子法ADMM译码中的容差值设置为10-5、超松弛参数设置为1.9、初始迭代次数设置为0；
(1c)在[1，20]范围内，设置交替方向乘子法ADMM迭代译码停止时的迭代次数，在[0.05，0.10]范围内，设置交替方向乘子法ADMM迭代译码中码字比特的硬判决变化率CRHD的阈值；
(1d)将拉格朗日乘子向量λ初始值设置为全0向量、辅助向量z初始值设置为全0向量；
(2)计算每个变量节点的初始信息：
xq＝γi
其中，xq表示第q个变量节点的初始信息，q的取值范围为[1，N]，N表示低密度奇偶校验LDPC码的码字长度，γi表示初始化中计算的低密度奇偶检验LDPC码的码字比特i对应的代价参数，i的取值范围为[1，N]，q与i的取值对应相等；
(3)更新每个变量节点的当前信息；
(4)更新辅助变量；
(5)更新拉格朗日乘子向量；
(6)获取每个码字比特的硬判决值：
<mrow>
<msubsup>
<mi>D</mi>
<mi>i</mi>
<mi>k</mi>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0.5</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo><</mo>
<mn>0.5</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，k表示第k次迭代次数，表示第i个码字比特在第k次迭代结束后对应的硬判决HD值，xi表示第i个变量节点的当前信息；
(7)计算交替方向乘子法ADMM译码中每两次连续迭代之间发生变化的码字比特总数；
(8)计算码字比特的硬判决变化率；
(9)判断是否满足早停止译码的条件，若是，则执行步骤(10)，否则，将当前迭代次数加1后执行步骤(3)；
(10)停止交替方向乘子法ADMM译码。
2.根据权利要求1所述的低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，其特征在于，步骤(1a)中所述的计算低密度奇偶检验LDPC码的码字比特对应的代价参数公式如下：
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>Pr</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>w</mi>
</msub>
<mo>|</mo>
<msub>
<mi>c</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Pr</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>w</mi>
</msub>
<mo>|</mo>
<msub>
<mi>c</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中，γi表示低密度奇偶检验LDPC码的第i个码字比特对应的代价参数，i的取值范围为[1，N]，log表示以2为底的对数操作，Pr(·)表示事件发生概率，rw表示加性高斯白噪声信道输出的第w位信息，w的取值范围为[1，N]，Pr(A|B)表示当事件B发生时事件A发生的概率，ct表示发送的第t位信息，t的取值范围为[1，N]，i、w与t的取值对应相等。
3.根据权利要求1所述的低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，其特征在于，步骤(3)中所述的更新每个变量节点的当前信息是按照下述公式实现的：
<mrow>
<msub>
<mi>x</mi>
<mi>q</mi>
</msub>
<mo>=</mo>
<msub>
<mo>&Pi;</mo>
<mrow>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</msub>
<mfrac>
<mn>1</mn>
<msub>
<mi>d</mi>
<mi>w</mi>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>v</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<msubsup>
<mi>z</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>&mu;</mi>
</mfrac>
<msubsup>
<mi>&lambda;</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>&mu;</mi>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>&alpha;</mi>
<mi>sgn</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
其中，xq表示第q个变量节点的当前信息，Π[0,1]表示在区间[0，1]内做欧几里德投影操作，dw表示与第w个变量节点相邻的校验节点的总数，∑表示求和操作，j表示第j个校验节点，∈表示属于符号，Nv(w)表示与第w个变量节点相邻的校验节点的索引集，表示第j个校验节点对应的辅助向量zj中与第t个变量节点对应的向量位，μ表示针对具体信道、具体码字，根据误码率大小优化得到的增广拉格朗日参数，表示第j个校验节点对应的拉格朗日乘子向量λj中与第t个变量节点对应的向量位，α表示针对具体信道、具体码字，根据误码率大小优化得到的罚函数的惩罚参数，sgn(·)表示符号函数，xi表示第i个变量节点更新前的信息，q、i、w与t的取值对应相等。
4.根据权利要求1所述的低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，其特征在于，步骤(4)中所述的更新辅助变量是按照下述公式实现的：
<mrow>
<msub>
<mi>z</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<msub>
<mo>&Pi;</mo>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mi>x</mi>
<mo>+</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mi>j</mi>
</msub>
<mi>&mu;</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中，zj表示第j个校验节点对应的辅助向量，表示在第j个校验节点对应的校验多胞体Pj上做欧几里德投影操作，Tj表示第j个校验节点对应的转换矩阵，x表示变量节点消息向量，λj表示第j个校验节点对应的拉格朗日乘子向量。
5.根据权利要求1所述的低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，其特征在于，步骤(5)中所述的更新拉格朗日乘子向量是按照下述公式实现的：
λr＝λj+μ(Tjx-zj)
其中，λr表示第r个校验节点对应的当前拉格朗日乘子向量，r与j的取值对应相等。
6.根据权利要求1所述的低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，其特征在于，步骤(7)中所述的计算交替方向乘子法ADMM译码中每两次连续迭代之间发生变化的码字比特总数是按照下述公式实现的：
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<mo>|</mo>
<msubsup>
<mi>D</mi>
<mi>i</mi>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>D</mi>
<mi>i</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>|</mo>
</mrow>
其中，Num表示交替方向乘子法ADMM译码中每两次连续迭代之间发生变化的码字比特总数，i表示低密度奇偶校验LDPC码中的第i个码字比特，|·|表示取绝对值操作，表示对应第i个码字比特在第k-1次迭代结束后对应的硬判决HD值。
7.根据权利要求1所述的低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，其特征在于，步骤(8)中所述的计算码字比特的硬判决变化率是按照下述公式实现的：
<mrow>
<mi>C</mi>
<mi>R</mi>
<mi>H</mi>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
</mrow>
<mi>N</mi>
</mfrac>
</mrow>
其中，CRHD表示交替方向乘子法ADMM译码中每两次连续迭代之间码字比特的硬判决变化率。
8.根据权利要求1所述的低密度奇偶校验LDPC码交替方向乘子法ADMM迭代译码的早停止方法，其特征在于，步骤(9)中所述的译码早停止条件是指满足以下三个条件之一的情形：
条件1：交替方向乘子法ADMM译码中，每两次连续迭代之间发生变化的码字比特总数为0；
条件2：不等式k≥Ite与CRHD≥Thr同时成立，其中，Ite表示初始化时所选取的交替方向乘子法ADMM迭代译码停止时的迭代次数，Thr表示初始化时所选取的交替方向乘子法ADMM迭代译码中码字比特的硬判决变化率CRHD的阈值；
条件3：交替方向乘子法ADMM迭代次数达到20。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710798802.5A CN107689801B (zh) | 2017-09-07 | 2017-09-07 | Ldpc码admm迭代译码的早停止方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710798802.5A CN107689801B (zh) | 2017-09-07 | 2017-09-07 | Ldpc码admm迭代译码的早停止方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107689801A true CN107689801A (zh) | 2018-02-13 |
| CN107689801B CN107689801B (zh) | 2019-10-25 |
Family
ID=61155169
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710798802.5A Active CN107689801B (zh) | 2017-09-07 | 2017-09-07 | Ldpc码admm迭代译码的早停止方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107689801B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108964669A (zh) * | 2018-07-06 | 2018-12-07 | 西安电子科技大学 | 基于度分解和交替乘子法的ldpc码二次规划译码方法 |
| CN109743062A (zh) * | 2019-01-03 | 2019-05-10 | 宝鸡文理学院 | 一种低密度奇偶检验码admm-lp译码的早停止方法 |
| CN110661532A (zh) * | 2019-11-12 | 2020-01-07 | 西安电子科技大学 | 基于多元ldpc码噪声增强的符号翻转译码方法 |
| CN110995277A (zh) * | 2019-12-06 | 2020-04-10 | 浙江大学 | 一种多层神经网络辅助的罚对偶分解信道译码方法 |
| CN112910472A (zh) * | 2021-01-21 | 2021-06-04 | 西安电子科技大学 | 基于2范数箱型admm的ldpc码惩罚译码方法 |
| CN113098528A (zh) * | 2021-03-16 | 2021-07-09 | 上海微波技术研究所（中国电子科技集团公司第五十研究所） | 基于ldpc译码的早停方法和系统 |
| CN114759930A (zh) * | 2022-04-11 | 2022-07-15 | 宝鸡文理学院 | 基于交替方向乘子法ldpc码的lp早停译码方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104092468A (zh) * | 2014-07-07 | 2014-10-08 | 西安电子科技大学 | 基于加速交替方向乘子法的ldpc码线性规划译码方法 |
| CN104753542A (zh) * | 2015-04-23 | 2015-07-01 | 西安工业大学 | 用于ldpc码的比特翻转和线性规划组合译码方法 |
| CN105530014A (zh) * | 2015-12-30 | 2016-04-27 | 西安电子科技大学 | 基于简化投影算子的ldpc码交替方向乘子译码方法 |
-
2017
- 2017-09-07 CN CN201710798802.5A patent/CN107689801B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104092468A (zh) * | 2014-07-07 | 2014-10-08 | 西安电子科技大学 | 基于加速交替方向乘子法的ldpc码线性规划译码方法 |
| CN104753542A (zh) * | 2015-04-23 | 2015-07-01 | 西安工业大学 | 用于ldpc码的比特翻转和线性规划组合译码方法 |
| CN105530014A (zh) * | 2015-12-30 | 2016-04-27 | 西安电子科技大学 | 基于简化投影算子的ldpc码交替方向乘子译码方法 |
Non-Patent Citations (5)
| Title |
|---|
| C.-H. LIN等: ""Efficient layer stopping technique for layered LDPC decoding"", 《ELECTRONICS LETTERS》 * |
| JIANGPENG LI等: ""Memory efficient layered decoder design with early termination for LDPC codes"", 《2011 IEEE INTERNATIONAL SYMPOSIUM OF CIRCUITS AND SYSTEMS (ISCAS)》 * |
| JIN LI等: ""Early stopping for LDPC decoding: convergence of mean magnitude (CMM)"", 《IEEE COMMUNICATIONS LETTERS》 * |
| M. FERRARI等: ""Safe Early Stopping for Layered LDPC Decoding"", 《IEEE COMMUNICATIONS LETTERS》 * |
| YU-HSIN CHEN等: ""A Channel-Adaptive Early Termination strategy for LDPC decoders"", 《2009 IEEE WORKSHOP ON SIGNAL PROCESSING SYSTEMS》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108964669A (zh) * | 2018-07-06 | 2018-12-07 | 西安电子科技大学 | 基于度分解和交替乘子法的ldpc码二次规划译码方法 |
| CN108964669B (zh) * | 2018-07-06 | 2021-07-06 | 西安电子科技大学 | 基于度分解和交替乘子法的ldpc码二次规划译码方法 |
| CN109743062A (zh) * | 2019-01-03 | 2019-05-10 | 宝鸡文理学院 | 一种低密度奇偶检验码admm-lp译码的早停止方法 |
| CN110661532A (zh) * | 2019-11-12 | 2020-01-07 | 西安电子科技大学 | 基于多元ldpc码噪声增强的符号翻转译码方法 |
| CN110661532B (zh) * | 2019-11-12 | 2023-02-10 | 西安电子科技大学 | 基于多元ldpc码噪声增强的符号翻转译码方法 |
| CN110995277A (zh) * | 2019-12-06 | 2020-04-10 | 浙江大学 | 一种多层神经网络辅助的罚对偶分解信道译码方法 |
| CN110995277B (zh) * | 2019-12-06 | 2021-06-01 | 浙江大学 | 一种多层神经网络辅助的罚对偶分解信道译码方法 |
| CN112910472A (zh) * | 2021-01-21 | 2021-06-04 | 西安电子科技大学 | 基于2范数箱型admm的ldpc码惩罚译码方法 |
| CN113098528A (zh) * | 2021-03-16 | 2021-07-09 | 上海微波技术研究所（中国电子科技集团公司第五十研究所） | 基于ldpc译码的早停方法和系统 |
| CN113098528B (zh) * | 2021-03-16 | 2023-01-24 | 上海微波技术研究所（中国电子科技集团公司第五十研究所） | 基于ldpc译码的早停方法和系统 |
| CN114759930A (zh) * | 2022-04-11 | 2022-07-15 | 宝鸡文理学院 | 基于交替方向乘子法ldpc码的lp早停译码方法 |
| CN114759930B (zh) * | 2022-04-11 | 2024-05-10 | 宝鸡文理学院 | 基于交替方向乘子法ldpc码的lp早停译码方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107689801B (zh) | 2019-10-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107689801A (zh) | Ldpc码admm迭代译码的早停止方法 | |
| CN103208995B (zh) | 一种低密度奇偶校验码译码的提前终止方法 | |
| CN109257148A (zh) | 一种基于高斯近似阈值判断的极化码bp译码方法 | |
| CN104025459A (zh) | 译码处理方法及译码器 | |
| CN108847850A (zh) | 一种基于crc-sscl的分段极化码编译码方法 | |
| CN107248866A (zh) | 一种降低极化码译码时延的方法 | |
| CN106888026A (zh) | 基于lsc‑crc译码的分段极化码编译码方法及系统 | |
| CN104092468B (zh) | 基于加速交替方向乘子法的ldpc码线性规划译码方法 | |
| CN106100795A (zh) | 一种基于倒位码字Plotkin 构造和信息位重休眠的Polar 码编码协作 | |
| CN101656541A (zh) | Rs码的译码方法和装置 | |
| Fathollahi et al. | Sparse multi-decoder recursive projection aggregation for Reed-Muller codes | |
| CN104218955A (zh) | 基于比特翻转的ldpc码局部搜索译码方法 | |
| CN107612560A (zh) | 基于部分信息比特似然比的极化码早期迭代停止方法 | |
| CN105530014A (zh) | 基于简化投影算子的ldpc码交替方向乘子译码方法 | |
| CN104393877A (zh) | 基于加权的非规则ldpc码线性规划译码方法 | |
| Hu et al. | A comparative study of polar code decoding algorithms | |
| Deng et al. | Reduced-complexity deep neural network-aided channel code decoder: A case study for BCH decoder | |
| CN103152057B (zh) | 一种基于双归一化修正因子的ldpc译码器及译码方法 | |
| Yan et al. | Simplified early stopping criterion for belief-propagation polar code decoder based on frozen bits | |
| CN104184480B (zh) | 一种降低复杂度的ldpc改进译码方法 | |
| CN102832954B (zh) | 一种基于软信息平均最小值的Turbo码译码迭代停止方法 | |
| US8984377B2 (en) | Stopping methods for iterative signal processing | |
| CN103595424A (zh) | 分量译码方法、译码器及Turbo译码方法、装置 | |
| CN104753542B (zh) | 用于ldpc码的比特翻转和线性规划组合译码方法 | |
| CN104242956A (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 |