CN112491424A - 低密度校验码的改进admm惩罚译码方法 - Google Patents
低密度校验码的改进admm惩罚译码方法 Download PDFInfo
- Publication number
- CN112491424A CN112491424A CN202011370520.3A CN202011370520A CN112491424A CN 112491424 A CN112491424 A CN 112491424A CN 202011370520 A CN202011370520 A CN 202011370520A CN 112491424 A CN112491424 A CN 112491424A
- Authority
- CN
- China
- Prior art keywords
- vector
- decoding
- iteration
- penalty
- admm
- 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.)
- Withdrawn
Links
Images
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
Landscapes
- Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Error Detection And Correction (AREA)
Abstract
本发明公开了一种低密度校验码的改进ADMM惩罚译码方法,主要解决现有译码方法复杂度高、译码速度慢的问题,其实现步骤为:(1)利用接收端接收到的消息向量计算对数似然比向量;(2)设置ADMM惩罚译码参数,并对参数进行初始化;(3);依次计算第k+1次迭代的中间向量、解向量、辅助向量以及拉格朗日乘子向量;(4)进行译码终止条件判断;(5)终止译码,并输出译码结果。本发明针对LDPC码译码问题采用改进罚函数,结合并替代了原有罚函数实现ADMM惩罚译码,从而有效提升了译码性能,且明显降低了译码复杂度,缩减了平均译码时间。
Description
技术领域
本发明属于通信技术领域,涉及编译码技术,具体为一种低密度校验码的改进交替方向乘子法(Alternating Direction Method ofMultipliers,ADMM)惩罚译码方法,可用于光纤通信、磁存储、卫星数字视频及声频广播。
背景技术
Gallager博士于1962年最早提出了低密度校验(Low-Density Parity-Check,LDPC)码的概念,这是一种线性分组码,其校验矩阵稀疏,且具有逼近香农限的良好性能。因此,许多编译码领域的工作人员对LDPC码产生了极大的兴趣,对该码进行了广泛的关注。目前,LDPC码在现代通信等领域得到了比较广泛的应用。
使用凸优化理论解决线性规划(LinearProgramming,LP)问题的LP译码方法是LDPC码的一类常见的译码方法,其优点是具有最大似然认证特性,但是译码复杂度高。基于交替方向乘子法,Barman等于2011年提出了一种迭代LP译码方法(Barman S,Liu X S,Draper S C,et al.Decomposition methods for large scale LP decoding.AnnualAllerton Conference on Communication,Control,and Computing(Allerton),253-260,2011[C].Monticello:IEEE,2011.),该方法能够降低原有LP译码的复杂度,但其不足之处是在低信噪比区域译码性能较差。
为了解决LDPC码ADMM译码方法在低信噪比区域译码性能差的问题,通常在ADMM译码问题的目标函数中加入罚函数,Liu X S等人提出了一种基于ADMM的惩罚译码方法(LiuX S,Draper S C.The ADMM penalized decoder for LDPC codes[J].IEEE Transactionson Information Theory,2016,62(6):2966-2984.),主要介绍了l1罚函数g1(x)=-|x-0.5|和l2罚函数g2(x)=-(x-0.5)2两种常用的罚函数,该方法在未增加译码复杂度的基础上改善了低信噪比区域的译码性能;然而,采用这两种罚函数进行ADMM惩罚译码的不足是译码复杂度仍然较高,译码速度有待进一步提高。
发明内容
针对上述已有技术的不足之处,本发明提供了一种低密度校验码的改进ADMM惩罚译码方法,用于解决现有译码方法复杂度高、译码速度慢的问题;在原有ADMM惩罚译码技术的基础上,通过设计一种改进的罚函数来提高译码性能,并且降低平均译码时间。
本发明实现上述目的具体步骤如下:
(1)发送端发送一个码长为n的LDPC码,通过加性高斯白噪声信道传输至接收端;
(2)接收端接收来自发送端的LDPC码,并将其记为消息向量r:
r={ri|i∈I},
其中,i表示变量节点;I表示LDPC码所有变量节点的集合,且I={1,2,…,n};
(3)利用消息向量r计算对数似然比向量γ:
γ={γi|i∈I}T,
其中,γi表示LDPC码ADMM惩罚译码问题模型的目标函数系数;{·}T表示转置操作;
(4)设计改进罚函数g(x),得到如下表达式:
(5)设置ADMM惩罚译码参数:
设定最大迭代次数N,且N为大于或等于500的整数;设置容差值ε为10-5、增广拉格朗日函数的罚参数μ为4.0,设置改进罚函数第一、三段的罚参数α为3.9、第二段的罚参数β为0.9;
(6)初始化:
(8)根据具体LDPC码的变量节点i得到其对应的度数di,利用第k+1次迭代的中间向量tk+1和变量节点i的度数di计算第k+1次迭代的解向量xk+1;具体计算公式如下:
其中,μ为增广拉格朗日函数的罚参数,为解向量xk+1的第i个分量,di表示变量节点i的度数,为中间向量tk+1的第i个分量,g′(x)为改进罚函数g(x)的一阶导数;∏[0,1](a)表示实数a在区间[0,1]上作投影运算;
(11)译码终止条件判断:
若满足终止条件,则进入步骤(12);否则,对k加1后返回步骤(7);
(12)将第k+1次迭代解向量xk+1作为译码结果输出,结束译码。
与现有技术相比,本发明具有以下优点:
第一、由于本发明在ADMM译码问题的目标函数中引入改进罚函数替代原有罚函数,从而能够更加灵活地惩罚伪码字实现ADMM惩罚译码,明显降低了译码复杂度,改善了译码性能并提高了译码速度;
第二、由于本发明采用的改进罚函数,其设计结合了原有罚函数,因此能够有效改善低信噪比区域的译码性能。
附图说明
图1是本发明的实现流程图;
图2是本发明与现有方法的译码性能仿真结果对比图;
图3是本发明与现有方法的译码平均迭代次数仿真结果对比图;
图4是本发明与现有方法的平均译码时间仿真结果对比图。
具体实施方式
一、技术原理
本发明是基于ADMM惩罚译码器求解LDPC码的线性规划译码问题进行的。
假设发送端发送一个码长为n的LDPC码,经过加性高斯白噪声信道传输后,接收端接收到的消息向量记为r={ri|i∈I},其中集合I={1,2,…,n}是LDPC码所有变量节点集合。利用向量r,根据如下公式计算对数似然比向量γ={γi|i∈I}T:
其中,Pr(·)表示括号内代表的事件发生的概率。将对数似然比γi作为线性规划目标函数的系数,引入罚函数的LDPC码LP译码问题模型:
其中,γi为目标函数系数;向量x={x1,x2,…,xn}表示ADMM惩罚译码的解向量;α为罚函数的罚参数;J={1,2,…,m}是LDPC码所有校验节点集合;dj是校验节点j所校验变量节点的个数;zj是维数为dj的辅助向量;是由维数为dj且所有含偶数个1的0-1向量所构成的校验多胞体;转换矩阵Tj是由校验节点j生成的,其定义为:假设校验节点j的邻接变量节点集合为其中则Tj中对于所有k∈dj的(k,ik)处的值为1,其余位置的值为0;g(x)是本发明设计的改进罚函数,且
上述问题对应的增广拉格朗日函数为:
其中,μ为增广拉格朗日函数的罚参数,λj为维数为dj的拉格朗子乘子向量,符号||·||2表示2-范数运算。
二、实施例
以IEEE802.16E标准中码率为0.75的(1152,288)非规则LDPC码为例进行译码,结合附图对本发明进行详细描述。
参照图1,本发明提供一种低密度校验码的改进ADMM惩罚译码方法,其实现步骤如下:
步骤1:发送端发送一个码长为n的LDPC码,通过加性高斯白噪声信道传输至接收端;
步骤2:接收端接收来自发送端的LDPC码,并将其记为消息向量r:
r={ri|i∈I},
其中,i表示变量节点;I表示LDPC码所有变量节点的集合,且I={1,2,…,n};
步骤3:计算对数似然比向量。
利用消息向量r计算对数似然比向量γ:
γ={γi|i∈I}T,
其中,γi作为LDPC码ADMM惩罚译码问题模型的目标函数系数;{·}T表示转置操作;
目标函数系数γi,按照下式得到:
其中,Pr(·)表示括号内代表的事件·的发生概率。
步骤4:基于现有罚函数,设计得到改进罚函数g(x):
步骤5:设置ADMM惩罚译码参数。
ADMM惩罚译码参数包括解向量x,拉格朗日乘子向量yj(j表示校验节点),辅助向量zj,中间向量t,容差值ε,迭代次数k,最大迭代次数N,罚函数的罚参数α、β,以及增广拉格朗日函数的罚参数μ。
译码的最大迭代次数设置得越大,译码性能会越好,但是译码时间也会越长。综合考虑译码的时间和性能,本发明的实施例以码率为0.75的(1152,288)非规则LDPC码为例,采用本发明所设计的改进罚函数进行ADMM惩罚译码。将译码的最大迭代次数N设置为1000;将容差值ε设置为10-5;由于改进罚函数为分段函数,为得到更好的译码效果,将第一、三段的罚参数α设置为3.9,第二段的罚参数β设置为0.9;增广拉格朗日函数的罚参数μ设置为4.0。
步骤6:初始化:
步骤7:计算第k+1次迭代的中间向量tk+1。
其中,为中间向量tk+1的第i个分量,j∈Nv(i)表示变量节点i所关联的校验节点j,和分别表示与第k次迭代辅助向量和第k次迭代拉格朗日乘子向量中第i个变量节点有关联的值,μ为增广拉格朗日函数的罚参数。
步骤8:计算第k+1次迭代解向量xk+1。
利用第k+1次迭代中间向量tk+1和变量节点i的度数di,通过如下公式计算第k+1次迭代解向量xk+1:
其中,度数di是根据具体LDPC码的变量节点i得到;为解向量xk+1的第i个分量,di表示变量节点i的度数,为中间向量tk+1的第i个分量,g′(x)为本发明所设计的改进罚函数g(x)的一阶导数,Π[0,1](a)表示实数a在区间[0,1]上作投影运算;
实数a在区间[0,1]上作投影运算的运算规则如下:
其中,Tj是由校验节点j生成的转换矩阵,是由维数为dj且所有含偶数个1的0-1向量所构成的校验多胞体,表示向量到校验多胞体的欧几里德投影运算,称为欧几里德投影算子,dj是校验节点j所校验变量节点的个数。
步骤11:译码终止条件判断。
判断是否满足ADMM译码终止条件:若满足终止条件,则进入步骤12;否则,对k加1后返回步骤7。
步骤12:将第k+1次迭代解向量xk+1作为译码结果输出,结束译码。
下面结合仿真实验对本发明的效果作进一步的说明。
1.仿真条件:
本发明的仿真实验在64位Win7操作系统、Intel(R)Core(TM)i53470处理器、3.2GHz主频、8GB内存的硬件环境和VC6.0的软件环境下进行的。
2.仿真内容:
在加性高斯白噪声信道下,分别用本发明设计的改进ADMM惩罚译码方法和原有ADMM惩罚译码方法对IEEE802.16E标准中码率为0.75的(1152,288)非规则LDPC码进行译码,其译码性能对比图如图2所示、译码平均迭代次数对比如图3所示、平均译码时间对比图4所示。
3.仿真结果分析:
参照图2,其中给出的3条曲线具体为:
带圆环形的曲线表示在加性高斯白噪声信道下,采用l1罚函数的原有ADMM惩罚译码方法的纠错性能仿真曲线;
带方框形的曲线表示在加性高斯白噪声信道下,采用l2罚函数的原有ADMM惩罚译码方法的纠错性能仿真曲线;
带星形的曲线表示在加性高斯白噪声信道下,采用本发明设计的改进ADMM惩罚译码方法的纠错性能仿真曲线。
由图2可以看出,本发明设计的改进ADMM惩罚译码方法的纠错性能明显优于原有ADMM惩罚译码方法的纠错性能。
参照图3,其中给出的3条曲线具体为:
带圆环形的曲线表示在加性高斯白噪声信道下,采用l1罚函数的原有ADMM惩罚译码方法译码的平均迭代次数仿真曲线;
带方框形的曲线表示在加性高斯白噪声信道下,采用l2罚函数的原有ADMM惩罚译码方法译码的平均迭代次数仿真曲线;
带星形的曲线表示在加性高斯白噪声信道下,采用本发明设计的改进ADMM惩罚译码方法译码的平均迭代次数仿真曲线。
由图3可以看出,本发明设计的改进ADMM惩罚译码方法译码的平均迭代次数明显少于原有ADMM惩罚译码方法译码的平均迭代次数。
参照图4,其中给出的3条曲线具体为:
带圆环形的曲线表示在加性高斯白噪声信道下,采用l1罚函数的原有ADMM惩罚译码方法的平均译码时间仿真曲线;
带方框形的曲线表示在加性高斯白噪声信道下,采用l2罚函数的原有ADMM惩罚译码方法译码的译码时间仿真曲线;
带星形的曲线表示在加性高斯白噪声信道下,采用本发明设计的改进ADMM惩罚译码方法的平均译码时间仿真曲线。
由图4可以看出,本发明设计的改进ADMM惩罚译码方法的平均译码时间明显少于原有ADMM惩罚译码方法的平均译码时间。
综上,本发明设计了一种改进ADMM惩罚译码方法,与原有LDPC码ADMM惩罚译码方法相比,该方法不仅取得了更好的译码性能,并且降低了平均译码时间。
上述仿真分析证明了本发明所提方法的正确性与有效性。
本发明未详细说明部分属于本领域技术人员公知常识。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,显然对于本领域的专业人员来说,在了解了本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变,但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。
Claims (7)
1.一种低密度校验码的改进ADMM惩罚译码方法,其特征在于,包括如下步骤:
(1)发送端发送一个码长为n的LDPC码,通过加性高斯白噪声信道传输至接收端;
(2)接收端接收来自发送端的LDPC码,并将其记为消息向量r:
r={ri|i∈I},
其中,i表示变量节点;I表示LDPC码所有变量节点的集合,且I={1,2,…,n};
(3)利用消息向量r计算对数似然比向量γ:
γ={γi|i∈I}T,
其中,γi表示LDPC码ADMM惩罚译码问题模型的目标函数系数;{·}T表示转置操作;
(4)设计改进罚函数g(x),得到如下表达式:
(5)设置ADMM惩罚译码参数:
设定最大迭代次数N,且N为大于或等于500的整数;设置容差值ε为10-5、增广拉格朗日函数的罚参数μ为4.0,设置改进罚函数第一、三段的罚参数α为3.9、第二段的罚参数β为0.9;
(6)初始化:
(8)根据具体LDPC码的变量节点i得到其对应的度数di,利用第k+1次迭代的中间向量tk +1和变量节点i的度数di计算第k+1次迭代的解向量xk+1;具体计算公式如下:
其中,μ为增广拉格朗日函数的罚参数,为解向量xk+1的第i个分量,di表示变量节点i的度数,为中间向量tk+1的第i个分量,g′(x)为改进罚函数g(x)的一阶导数;Π[0,1](a)表示实数a在区间[0,1]上作投影运算;
(11)译码终止条件判断:
若满足终止条件,则进入步骤(12);否则,对k加1后返回步骤(7);
(12)将第k+1次迭代解向量xk+1作为译码结果输出,结束译码。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011370520.3A CN112491424A (zh) | 2020-11-30 | 2020-11-30 | 低密度校验码的改进admm惩罚译码方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011370520.3A CN112491424A (zh) | 2020-11-30 | 2020-11-30 | 低密度校验码的改进admm惩罚译码方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112491424A true CN112491424A (zh) | 2021-03-12 |
Family
ID=74937036
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011370520.3A Withdrawn CN112491424A (zh) | 2020-11-30 | 2020-11-30 | 低密度校验码的改进admm惩罚译码方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112491424A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113055069A (zh) * | 2021-03-17 | 2021-06-29 | 西安电子科技大学 | 基于admm的单组多播单流系统预编码矩阵设计的方法 |
-
2020
- 2020-11-30 CN CN202011370520.3A patent/CN112491424A/zh not_active Withdrawn
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113055069A (zh) * | 2021-03-17 | 2021-06-29 | 西安电子科技大学 | 基于admm的单组多播单流系统预编码矩阵设计的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7395487B2 (en) | Common circuitry supporting both bit node and check node processing in LDPC (Low Density Parity Check) decoder | |
KR100651036B1 (ko) | 공간다중화 시스템에서 부분적인 구복호 기반 복잡도를줄이는 반복 복호 수신 장치 및 그 방법 | |
CN108847848B (zh) | 一种基于信息后处理的极化码的bp译码算法 | |
KR100804793B1 (ko) | 저밀도 패러티 검사 복호기에서의 검사 노드 갱신 방법 | |
CN108964669B (zh) | 基于度分解和交替乘子法的ldpc码二次规划译码方法 | |
US7971131B1 (en) | System and method for iterative denoising and error correction decoding | |
US20040153942A1 (en) | Soft input soft output decoder for turbo codes | |
KR102231278B1 (ko) | 이진 로그를 사용하는 저밀도 패리티 검사 복호기 및 그것의 복호 방법 | |
CN110784232B (zh) | 一种空间耦合ldpc码滑窗译码方法 | |
CN108988872B (zh) | 基于分层最小和算法的ldpc译码方法 | |
CN109586732B (zh) | 中短码ldpc编解码系统和方法 | |
CN110417512B (zh) | 一种用于cpm通信系统的联合迭代译码方法 | |
CN110830049B (zh) | 一种基于密度进化改进偏移最小和的ldpc译码方法 | |
Jin et al. | Reliability-based soft-decision decoding with multiple biases | |
WO2007068554A1 (en) | Serial concatenation scheme and its iterative decoding using an inner ldpc and an outer bch code | |
CN112491424A (zh) | 低密度校验码的改进admm惩罚译码方法 | |
Muller et al. | Parallelism efficiency in convolutional turbo decoding | |
Dou et al. | Soft-decision based sliding-window decoding of staircase codes | |
CN114759930B (zh) | 基于交替方向乘子法ldpc码的lp早停译码方法 | |
CN108365918B (zh) | 一种基于有效集中准则的多元ldpc码译码方法 | |
CN111835363A (zh) | 基于交替方向乘子法的ldpc码译码方法 | |
CN109921877B (zh) | 一种awgn信道下基于cuda架构的ldpc码并行译码方法 | |
CN112332868A (zh) | 一种基于DVB-RCS2的turbo并行译码方法 | |
CN108092672B (zh) | 一种基于折叠调度的bp译码方法 | |
CN111555759A (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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20210312 |