CN107612560B - 基于部分信息比特似然比的极化码早期迭代停止方法 - Google Patents

基于部分信息比特似然比的极化码早期迭代停止方法 Download PDF

Info

Publication number
CN107612560B
CN107612560B CN201710827229.6A CN201710827229A CN107612560B CN 107612560 B CN107612560 B CN 107612560B CN 201710827229 A CN201710827229 A CN 201710827229A CN 107612560 B CN107612560 B CN 107612560B
Authority
CN
China
Prior art keywords
iteration
decoding
likelihood ratio
information bit
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
CN201710827229.6A
Other languages
English (en)
Other versions
CN107612560A (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.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
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 Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN201710827229.6A priority Critical patent/CN107612560B/zh
Publication of CN107612560A publication Critical patent/CN107612560A/zh
Application granted granted Critical
Publication of CN107612560B publication Critical patent/CN107612560B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Error Detection And Correction (AREA)

Abstract

本发明公开了一种基于部分信息比特似然比的极化码早期迭代停止方法,包括如下步骤:S1)预设BP译码的最大迭代次数;S2)利用BP译码算法对极化编码信息进行译码;S3)一次迭代完成后,对相邻两次迭代译码输出部分信息比特似然比进行比较;如果其中相同的信息比特似然比在预设的比较空间中的比例达到预设阈值,则停止迭代并输出当前迭代所得到的译码结果,否则,继续进行迭代,直至达到预设的最大迭代次数。本发明在每一次迭代过程中,均对前后两次迭代译码器输出的部分信息比特似然比进行判断;如果相同的信息比特似然比在比较空间中的比例达到预设阈值则停止,从而大大降低译码的计算复杂度和译码延时,有效降低硬件资源消耗。

Description

基于部分信息比特似然比的极化码早期迭代停止方法
技术领域
本发明涉及一种极化码处理方法,尤其涉及一种基于部分信息比特似然比的极化码早期迭代停止方法。
背景技术
基于信道极化现象,2008年Arikan在ISIT会议上提出了一种容量“可达”的码字,称为Polar码,其在论文中严格证明了在二进制输入离散无记忆信道中当码长趋于无穷时信道容量可以达到香农界。在译码端,Arikan同时提出了一种串行译码方法,称串行抵消算法(Successive cancellation,SC)。由于其串行译码结构,SC算法译码延迟较高。为了减少译码延迟,迭代的置信度传播(Belief propagation,BP)算法也被用于Polar码译码。
BP算法是并行,其译码延迟为O(IlogN),其中I为迭代次数。由此可知,减少迭代次数对于减少BP算法译码延迟非常重要。传统BP算法到达预设最大迭代次数才会停止,实际中正确译码结果在迭代早期就会得到,所以相关研究者提出了多种早期迭代停止方法以避免多余的迭代,例如minLLR、G-Martix、LMA和CA等方法。早期迭代停止方法在译码过程中对译码结果进行检测,如果满足停止标准则停止迭代输出译码结果,能够有效减少平均迭代次数。已有的早期迭代停止方法一般复杂度较高,提出一种复杂度更低的早期迭代停止方法是重要的研究方向。
发明内容
本发明所要解决的技术问题是提供一种基于部分信息比特似然比的极化码早期迭代停止方法,能够大大降低译码的计算复杂度和译码延时,同时便于硬件实现。
本发明为解决上述技术问题而采用的技术方案是提供一种基于部分信息比特似然比的极化码早期迭代停止方法,包括如下步骤:S1)预设BP译码的最大迭代次数;S2)利用BP译码算法对极化码编码信息进行译码;S3)一次迭代完成后,对相邻两次迭代BP译码器输出的部分信息比特似然比进行比较;如果其中相同的信息比特似然比在预设的比较空间中的比例达到预设阈值,则停止迭代并输出当前迭代步所得到的译码结果,否则,继续进行迭代,直至达到预设的最大迭代次数。
上述的基于部分信息比特似然比的极化码早期迭代停止方法,其中,所述极化码的码长为N,所述极化码包含K位信息位,设集合A为信息位的集合,集合
Figure BDA0001407781700000021
称为比较空间,所述集合
Figure BDA0001407781700000022
包含A中错误概率最小的M位信息位,M称为比较空间容量,K≥M>0,每个极化信道的错误概率通过高斯近似方法仿真得到,所述步骤S3)若满足如下不等式,则停止迭代:
Figure BDA0001407781700000023
其中,R为比例阈值,其取值范围为{R|0<R≤1},
Figure BDA0001407781700000024
表示BP算法因子图中第1列第i行节点的L信息。
上述的基于部分信息比特似然比的极化码早期迭代停止方法,其中,对于参数为(N,K)的极化码,其对应因子图由n=log2N阶计算单元和n+1列节点构成,每阶由N/2个处理单元构成,(i,j)表示从左起第i行,第j列的节点;每个节点从右到左传递通过节点(i,j)的软信息记为Li,j,从左到右传递通过节点(i,j)的软信息记为Ri,j,对因子图最左端的1列节点中的软信息进行硬判决可得到信息比特序列u的估计值
Figure BDA0001407781700000025
所述步骤S2)先向左传播更新节点中的Li,j,到达最左侧后开始向右传播更新节点中的Ri,j;迭代终止后,如果不是信息位则该位译码为0,否则依照最左端节点中的Li,1的符号判断信息位是0还是1。
上述的基于部分信息比特似然比的极化码早期迭代停止方法,其中,所述步骤S3)通过组合使用比较器、加法器与阈值比较器实现早期迭代停止机制;每次迭代后使用M个比较器比较当前迭代所得到的
Figure BDA0001407781700000026
与上次迭代所得到的
Figure BDA0001407781700000027
t-1次迭代的
Figure BDA0001407781700000028
从BP译码器的存储器中读取,t次迭代的
Figure BDA0001407781700000029
直接从译码器处理单元中获得;比较器的比较结果为{q1,q2,...,qM},若
Figure BDA00014077817000000210
则qi=1,否则qi=0;加法器用于计算
Figure BDA00014077817000000211
结果为
Figure BDA00014077817000000212
阈值比较器判断Q是否大于等于M*R,若大于等于则输出D=1,否则输出0;D=1则BP译码器停止迭代输出译码结果,D=0则BP译码器继续迭代直到达到预设的最大迭代次数。
上述的基于部分信息比特似然比的极化码早期迭代停止方法,其中,所述最大迭代次数预设为15~80。
本发明对比现有技术有如下的有益效果:本发明提供的基于部分信息比特似然比的极化码早期迭代停止方法,一次迭代完成后,对相邻两次迭代输出的属于比较空间的信息比特似然比进行比较;如果其中相同的信息比特似然比在比较空间中的比例达到预设阈值,则停止迭代并输出当前迭代步所得到的译码结果。当最大迭代次数为40次且Eb/N0=3.5dB时,与固定迭代40次的原始BP译码器相比,本发明能使平均迭代次数减少83.16%,有效降低了计算复杂度与译码延迟。minLLR标准的加法运算复杂度为N,比较运算复杂度为N;LMA标准的加法运算复杂度为2N,比较运算复杂度为N。与之相比本发明的加法运算复杂度仅为N/32,比较运算的复杂度为N/32+1,可有效降低早期迭代停止标准的硬件复杂度。
附图说明
图1为本发明参数为(8,4)的极化码因子图;
图2为本发明极化码因子图的基本单元示意图;
图3是本发明早期迭代停止流程示意图;
图4是本发明早期迭代停止模块硬件结构;
图5是本发明参数为(1024,512)的极化码,最大迭代次数为40的BP译码方法与原始BP译码器在不同信噪比信道下的平均迭代次数比较示意图;
图6本发明参数为(1024,512)的极化码,最大迭代次数为40的BP译码方法与原始BP译码器的译码性能比较示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的描述。
为了减小BP译码算法的译码复杂度,早期迭代停止算法十分重要。早期迭代停止算法是指在译码迭代过程中自适应地检测是否已经得到可靠的译码输出,如果条件满足可以立刻结束译码。早期迭代停止算法可以线性地降低译码复杂度和译码延时。本发明提供的基于部分信息比特似然比的极化码早期迭代停止方法,包括如下步骤:
S1)预设BP译码的最大迭代次数;
S2)利用BP译码算法对极化编码信息进行译码;
S3)一次迭代完成后,对相邻两次迭代输出的属于比较空间的信息比特似然比进行比较;如果其中相同的信息比特似然比在比较空间中的比例达到预设阈值,则停止迭代并输出当前迭代步所得到的译码结果,否则,继续进行迭代,直至达到预设的最大迭代次数。
本发明的具有早期迭代停止机制的极化码BP译码方法,利用BP译码算法对信道接收值进行译码,信息更新公式如下:
Figure BDA0001407781700000041
f(x,y)≈α*sign(x)sign(y)min(|x|,|y|) (2)
码长为N的Polar码包含K位信息位,设集合A为信息位的集合。集合
Figure BDA0001407781700000042
称为比较空间,包含A中错误概率最小的M(K≥M>0)位信息位,M称为比较空间容量,每个极化信道的错误概率可通过高斯近似方法仿真得到。在利用BP译码器对极化码进行译码的每一次迭代过程中,若满足不等式(3)则停止迭代,输出译码结果。
Figure BDA0001407781700000043
其中0≤R≤1,
Figure BDA0001407781700000044
表示BP算法因子图中第1列第i行节点的L信息。若不满足不等式(3),则继续进行迭代,直至达到预设的最大迭代次数。
相比现有技术方案,本发明能够在不造成译码性能损失的情况下显著减小译码迭代次数。对于(1024,512)极化码,当最大迭代次数为40次且Eb/E0=3.5dB时,本发明能使平均迭代次数降低83.16%,有效降低了计算复杂度与译码延迟。本发明的加法运算复杂度为N/32,比较运算复杂度为N/32+1。
本实施例采用参数为(1024,512)的Polar码进行测试,码长为N=1024,K=512,码率为0.5。使用高斯近似方法在信噪比为1.5dB下仿真得到1024个子信道的错误概率,错误概率最小的512个子信道位置构成集合A,用于传输信息,称为信息位;剩下的512个信道用于传输固定信息,称为冻结位。调制方式为二进制相移键控(Binary Phase Shift Keying,BPSK),信道为加性高斯白噪声信道(Additive White Gaussian Noise,AWGN)。码字
Figure BDA0001407781700000045
由长度为1024的
Figure BDA0001407781700000046
与生成矩阵G相乘得到。生成矩阵
Figure BDA0001407781700000047
Figure BDA0001407781700000048
表示矩阵
Figure BDA0001407781700000049
的log21024=10次克罗内克积。信道接收值Y1 N使用对数似然比(Log-likelihood ratio,LLR)形式表示。本实施例中采用1位符号位,4位整数位,2位小数位的7位量化方案,便于硬件实现。
参数为(1024,512)的极化码的因子图由log21024=10阶构成(码长为8的因子图如图1所示),其中每阶由N/2=512个处理单元构成(图2为单个处理单元示意图)。因子图最左端的一列对应信息比特u。(i,j)表示从左起第i行,第j列的节点。每个节点都有两种信息,本发明将从右到左传递通过节点(i,j)的信息记为Li,j,将从左到右传递通过节点(i,j)的信息记为Ri,j,这些信息以LLR形式相互传递更新。
在译码过程中,预设迭代最大迭代次数为40次,先对Ri,1和Li,n+1进行初始化。Li,n+1初始化为信道接受值Yi,Ri,1根据位置信息分别初始化为0和7位量化方案能表示的最大值15.75。具体公式如下:
Figure BDA0001407781700000051
Li,n+1=Yi (5)
本实施例中设α=0.9,依据式(1),(2)先向左传播更新节点中的Li,j,到达最左侧后开始向右传播开始更新Ri,j
本发明方法的基本流程如图3所示。本实施例中参数M设为32,β值设为1/4。为了降低译码的计算复杂度和译码延时,若Li,1,i∈S中在前后两次迭代中保持不变的比例大于1/4,此时的译码输出就可以认为是可靠的译码输出。具体公式如下:
Figure BDA0001407781700000052
如果不满足式(6),则判断译码次数是否达到40次,若达到最大迭代次数则终止迭代;如果没有,则继续更新信息,进行下一次迭代。
迭代终止后,如果是冻结位则该位译码为0,否则依照最左端节点中的Li,1的符号判断信息位是0还是1,得到译码结果。BP译码器硬判决依据的公式如下:
Figure BDA0001407781700000053
图4为本实施例的早期迭代停止模块硬件架构,每次迭代后使用M个比较器比较当前迭代所得到的
Figure BDA0001407781700000054
与上次迭代所得到的
Figure BDA0001407781700000055
t-1次迭代的
Figure BDA0001407781700000056
从BP译码器的存储器中读取,t次迭代的
Figure BDA0001407781700000057
直接从译码器处理单元中获得。比较器的比较结果为{q1,q2,...,qM},若
Figure BDA0001407781700000058
则qi=1,否则qi=0。加法器用于计算
Figure BDA0001407781700000061
结果为
Figure BDA0001407781700000062
本实施例中阈值R设为1/4,阈值比较器判断Q是否大于等于M*R,若大于等于则输出D=1,否则输出0。BP译码器根据早期迭代停止模块的输出结果判断是否停止迭代,即D=1时译码器停止迭代输出译码结果,D=0时继续迭代直到达到最大迭代次数。
下表为本发明的方法与另外两种早期迭代停止标准的计算复杂度对比。本实施例中M=N/32,则本发明的加法运算复杂度为N/32,比较运算的复杂度为N/32+1,有效降低了早期迭代停止标准的计算复杂度。
早期停止标准 minLLR LMA 本发明
加法运算 N 2N N/32
比较运算 N N N/32+1
图5显示了本实施与传统极化码BP译码方法在不同信噪比信道中的平均迭代次数。图中Eb/N0为信噪比,Average number of iterations表示平均迭代次数。当最大迭代次数为40次且Eb/N0=3.5dB时,与迭代次数固定为40次的原始BP译码器相比,本发明能使平均迭代次数减少83.16%,与minLLR(β=9.5)和LMA相比,迭代减少性能分别上升6.48%和8.14%。
图6显示了本实施及传统极化码BP译码器在高斯加性白噪声信道中的测试结果。图中横坐标Eb/N0为信噪比,图例中fer为误帧率,ber为误比特率,40Fixed Iteration表示迭代次数固定为40次。根据图6可看出本发明能在迭代次数显著少于传统BP译码器的情况下,达到和传统BP译码器一样的译码性能,并未造成译码性能损失。
虽然本发明已以较佳实施例揭示如上,然其并非用以限定本发明,任何本领域技术人员,在不脱离本发明的精神和范围内,当可作些许的修改和完善,因此本发明的保护范围当以权利要求书所界定的为准。

Claims (4)

1.一种基于部分信息比特似然比的极化码早期迭代停止方法,其特征在于,包括如下步骤:
S1)预设BP译码的最大迭代次数;
S2)利用BP译码算法对极化码编码信息进行译码;
S3)一次迭代完成后,对相邻两次迭代BP译码器输出的部分信息比特似然比进行比较;如果其中相同的信息比特似然比在预设的比较空间中的比例达到预设阈值,则停止迭代并输出当前迭代所得到的译码结果,否则,继续进行迭代,直至达到预设的最大迭代次数;
所述极化码的码长为N,所述极化码包含K位信息位,设集合A为信息位的集合,集合
Figure FDA0002443359430000011
称为比较空间,所述集合
Figure FDA0002443359430000012
包含A中错误概率最小的M位信息位,M称为比较空间容量,K≥M>0,每个极化信道的错误概率通过高斯近似方法仿真得到,所述步骤S3)若满足如下不等式,则停止迭代:
Figure FDA0002443359430000013
其中,R为比例阈值,其取值范围为{R|0<R≤1},
Figure FDA0002443359430000014
表示BP算法因子图中第1列第i行节点的L信息。
2.如权利要求1所述的基于部分信息比特似然比的极化码早期迭代停止方法,其特征在于,对于参数为(N,K)的极化码,其对应因子图由n=log2N阶计算单元和n+1列节点构成,每阶由N/2个处理单元构成,(i,j)表示从左起第i行,第j列的节点;每个节点从右到左传递通过节点(i,j)的软信息记为Li,j,从左到右传递通过节点(i,j)的软信息记为Ri,j,对因子图最左端的1列节点中的软信息进行硬判决可得到信息比特序列u的估计值
Figure FDA0002443359430000015
所述步骤S2)先向左传播更新节点中的Li,j,到达最左侧后开始向右传播更新节点中的Ri,j;迭代终止后,如果不是信息位则该位译码为0,否则依照最左端节点中的Li,1的符号判断信息位是0还是1。
3.如权利要求2所述的基于部分信息比特似然比的极化码早期迭代停止方法,其特征在于,所述步骤S3)通过组合使用比较器、加法器与阈值比较器实现早期迭代停止机制;每次迭代后使用M个比较器比较当前迭代所得到的
Figure FDA0002443359430000016
i∈S与上次迭代所得到的
Figure FDA0002443359430000017
i∈S,t-1次迭代的
Figure FDA0002443359430000018
从BP译码器的存储器中读取,t次迭代的
Figure FDA0002443359430000019
直接从译码器处理单元中获得;比较器的比较结果为{q1,q2,...,qM},若
Figure FDA00024433594300000110
则qi=1,否则qi=0;加法器用于计算
Figure FDA0002443359430000021
结果为
Figure FDA0002443359430000022
阈值比较器判断Q是否大于等于M*R,若大于等于则输出D=1,否则输出0;D=1则BP译码器停止迭代输出译码结果,D=0则BP译码器继续迭代直到达到预设的最大迭代次数。
4.如权利要求3所述的基于部分信息比特似然比的极化码早期迭代停止方法,其特征在于,所述最大迭代次数预设为15~80。
CN201710827229.6A 2017-09-14 2017-09-14 基于部分信息比特似然比的极化码早期迭代停止方法 Active CN107612560B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710827229.6A CN107612560B (zh) 2017-09-14 2017-09-14 基于部分信息比特似然比的极化码早期迭代停止方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710827229.6A CN107612560B (zh) 2017-09-14 2017-09-14 基于部分信息比特似然比的极化码早期迭代停止方法

Publications (2)

Publication Number Publication Date
CN107612560A CN107612560A (zh) 2018-01-19
CN107612560B true CN107612560B (zh) 2020-07-24

Family

ID=61063909

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710827229.6A Active CN107612560B (zh) 2017-09-14 2017-09-14 基于部分信息比特似然比的极化码早期迭代停止方法

Country Status (1)

Country Link
CN (1) CN107612560B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110752891B (zh) * 2018-07-24 2022-04-29 中兴通讯股份有限公司 极化码译码方法及装置、存储介质、电子装置
CN109450455B (zh) * 2018-10-26 2023-09-08 中国计量大学 高性能极化码信息位选取方案
CN109257148B (zh) * 2018-11-26 2020-05-22 北京理工大学 一种基于高斯近似阈值判断的极化码bp译码方法
CN110061747A (zh) * 2019-04-28 2019-07-26 中国石油大学(华东) 一种极化码的基于阈值的比特反转译码方法
CN110752852B (zh) * 2019-09-26 2023-10-03 浙江科睿微电子技术有限公司 极化码的bp译码方法、装置、***、设备及存储介质
CN110855298B (zh) * 2019-12-02 2023-03-31 重庆邮电大学 基于子信道冻结条件的低迭代次数极化码bp译码方法
CN111726202B (zh) * 2020-06-16 2022-05-31 杭州电子科技大学 一种极化码置信传播译码的提前终止迭代方法
CN113938227A (zh) * 2021-10-08 2022-01-14 天津津航计算技术研究所 一种基于迭代译码的信噪比强弱动态判定方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8810633B2 (en) * 2009-11-30 2014-08-19 Canon Kabushiki Kaisha Robust image alignment for distributed multi-view imaging systems
CN104539296A (zh) * 2015-01-21 2015-04-22 西安电子科技大学 一种基于提前终止迭代策略的极化码改进bp译码方法
CN105187073A (zh) * 2015-10-13 2015-12-23 东南大学 一种极化码的bp译码方法及装置
CN105262494A (zh) * 2015-10-13 2016-01-20 东南大学 一种具有迭代早停止机制的极化码bp译码方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8810633B2 (en) * 2009-11-30 2014-08-19 Canon Kabushiki Kaisha Robust image alignment for distributed multi-view imaging systems
CN104539296A (zh) * 2015-01-21 2015-04-22 西安电子科技大学 一种基于提前终止迭代策略的极化码改进bp译码方法
CN105187073A (zh) * 2015-10-13 2015-12-23 东南大学 一种极化码的bp译码方法及装置
CN105262494A (zh) * 2015-10-13 2016-01-20 东南大学 一种具有迭代早停止机制的极化码bp译码方法

Also Published As

Publication number Publication date
CN107612560A (zh) 2018-01-19

Similar Documents

Publication Publication Date Title
CN107612560B (zh) 基于部分信息比特似然比的极化码早期迭代停止方法
US7539920B2 (en) LDPC decoding apparatus and method with low computational complexity algorithm
CN109660264B (zh) 高性能极化码译码算法
CN108462558B (zh) 一种极化码scl译码方法、装置及电子设备
USRE44421E1 (en) Decoding apparatus for low-density parity-check codes using sequential decoding, and method thereof
CN108039891B (zh) 一种基于多级更新流程的极化码bp译码方法及装置
CN108847848B (zh) 一种基于信息后处理的极化码的bp译码算法
CN109286405B (zh) 一种低复杂度的极化码递进式比特翻转sc译码方法
US9178533B2 (en) Method of setting number of iteration counts of iterative decoding, and apparatus and methods of iterative decoding
CN110233628B (zh) 极化码的自适应置信传播列表译码方法
CN109951190B (zh) 一种自适应Polar码SCL译码方法及译码装置
CN111726202B (zh) 一种极化码置信传播译码的提前终止迭代方法
CN111835364B (zh) 一种极化码的低复杂度神经bp译码方法
CN114157309A (zh) 极化码译码方法、装置及***
WO2007044991A2 (en) Broadcast message passing decoding of low density parity check codes
CN109547035B (zh) 流水bp极化译码器硬件架构的建立方法及译码器硬件架构
US8930790B1 (en) Method and apparatus for identifying selected values from among a set of values
CN107707333B (zh) 一种基于码字估计值的极化码早期迭代停止方法及装置
CN111313913B (zh) 低时延的交叉调度极化码bp译码方法及装置
KR20090012189A (ko) Ldpc 부호의 성능 개선을 위한 스케일링 기반의 개선된min-sum 반복복호알고리즘을 이용한 복호 장치 및그 방법
KR20090064268A (ko) 가변 보정값을 이용한 복호화 장치 및 그 방법
CN111835363B (zh) 基于交替方向乘子法的ldpc码译码方法
Albayrak et al. Low-complexity early termination method for rateless soft decoder
KR102045438B1 (ko) Ldpc 부호의 복호를 위한 방법 및 장치
CN114422084A (zh) 基于高低llr占比的ad-scl译码方法

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