CN106877884B - 一种减少译码路径分裂的极化码译码方法 - Google Patents

一种减少译码路径分裂的极化码译码方法 Download PDF

Info

Publication number
CN106877884B
CN106877884B CN201710062953.4A CN201710062953A CN106877884B CN 106877884 B CN106877884 B CN 106877884B CN 201710062953 A CN201710062953 A CN 201710062953A CN 106877884 B CN106877884 B CN 106877884B
Authority
CN
China
Prior art keywords
decoding
path
node
bit
splitting
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
CN201710062953.4A
Other languages
English (en)
Other versions
CN106877884A (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.)
White Box Shanghai Microelectronics Technology Co ltd
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 CN201710062953.4A priority Critical patent/CN106877884B/zh
Publication of CN106877884A publication Critical patent/CN106877884A/zh
Application granted granted Critical
Publication of CN106877884B publication Critical patent/CN106877884B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/13—Linear codes

Landscapes

  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Error Detection And Correction (AREA)

Abstract

本发明公开了一种减少译码路径分裂的极化码译码方法，本方法在进行比特判决时，以该比特的对数似然比为置信度，当该比特位的置信度达到一定的门限时，不需要分裂出新的译码路径，以减少列表中待选译码路径的数量，从而降低译码复杂度。译码时根据译码置信度是否超过门限值、以及译码节点的类型确定译码路径的分裂策略。由于置信度门限值计算过一次之后可以存到表格中待查，本方法在译码器实时译码时不会增加额外的计算复杂度。本发明在比特位置信度很高的情况下可以不分裂新的译码路径，又不会损失误码率性能。该方法能够自适应信道噪声条件，当信噪比越高时，译码器性能增益越明。

Description

一种减少译码路径分裂的极化码译码方法
技术领域
本发明属于通信系统中的信道编译码技术领域，具体涉及通信系统信道译码器中一种基于串行抵消列表的极化码改进译码方法。
背景技术
极化码是一种新型编码方式，是目前唯一可理论证明在二进制离散无记忆信道下达到香农极限，并且具有可实用的线性复杂度编译码能力的信道编码技术，成为下一代通信系统(5G)中信道编码方案的强有力候选者。极化码在码长趋于无穷时，采用低复杂度的串行抵消(SC，Successive Cancellation)译码方法就可以无限逼近信道容量。极化码在长码时性能优异，但在中短码时需要采用串行抵消列表(SCL，Successive CancellationList)译码才能实现更好的性能。串行抵消译码是从译码树的根节点开始，通过递归计算的方式深入到叶节点，在叶节点进行判别。通过在译码树上定义特殊节点，当译码树深入到特殊节点时，直接判别而不再深入到叶节点，从而降低译码时延和运算复杂度。然而串行抵消列表译码方法的复杂度会随着列表长度增大而增加。现有的串行抵消列表译码方法在判决每个比特时都要分裂出新的译码路径，放入列表中作为待选路径。当信道条件较好(如信噪比较高或者子信道极化程度较高)时，仍需在判决每个比特时都要分裂出新的译码路径，串行抵消列表译码方法复杂度依然很高。
发明内容
发明目的：为了克服现有技术中存在的不足，本发明提供一种减少译码路径分裂的极化码译码方法，本发明在保证极化码译码器误码率性能的情况下，能够使串行抵消列表译码方法有效减少译码路径的平均数量，提高译码器吞吐量的极化码译码。该方法在进行比特判决时，以该比特的对数似然比为置信度，当该比特位的置信度达到一定的门限时，不需要分裂出新的译码路径，以减少列表中待选译码路径的数量，从而降低译码复杂度。译码时根据译码置信度是否超过门限值、以及译码节点的类型确定译码路径的分裂策略。由于置信度门限值计算过一次之后可以存到表格中待查，本方法在译码器实时译码时不会增加额外的计算复杂度；得益于该方法的自适应特性，该译码方法能够降低译码复杂度，提高译码器的吞吐量，在信噪比较高时性能增益尤为明显。
技术方案：为实现上述目的，本发明采用的技术方案为：
一种减少译码路径分裂的极化码译码方法，包括以下步骤：
步骤1，递归遍历译码树：对于一个码长为N，信息位长度为K的极化码，其译码树为深度为n＝log2N的二叉树，采用基于串行抵消的译码方法从根节点开始递归遍历整个译码树。其中，左子节点和右子节点以同样的方式向其子节点遍历，如果遇到以下4种类型的节点，则不再向下递归而直接译码后返回。这4种节点类型分别是：若其包含的叶节点全是休眠位比特，则记为Rate-0节点。若其包含的叶节点全是信息位比特，则记为Rate-1节点。若其包含的叶节点只有最后一位是信息位比特，则记为REP节点。若其包含的叶节点只有第一位是休眠位比特，则记为SPC节点。
步骤2，对步骤1中提到的4种类型的节点根据节点类型分别给出相应的译码策略，得到各节点相应的译码结果。
步骤3，根据译码时节点之间传递待译符号的对数似然比计算译码路径分裂门限，根据译码路径分裂门限决定是否分裂新的译码路径。译码时先估计信噪比Eb/No，当前Eb/No对应的译码路径分裂门限计算过一次之后就将译码路径分裂门限存入一张列表中。下一个码块译码时，如果Eb/No对应的译码路径分裂门限已存在则直接找到对应的列表待查，如果不存在则重新计算存到对应的列表中。
步骤4，译码路径分裂策略：译码开始时，列表里只有一条译码路径，每条译码路径都有一个评价路径优劣的度量值PM，初始路径的PM值设为0。第2步在译码时，根据当前信噪比估计值Eb/No，结合步骤3列表中的分裂门限，根据步骤1所述的节点类型确定不同的译码路径分裂策略，同时更新路径的PM值。分裂后产生的新路径添加到同一列表中，如果列表中路径数量多于设定的最大值Lmax，则只保留PM值最大的Lmax个路径。
步骤5，根据循环冗余校验选择最终译码结果：
如果译码树从根节点返回，对于列表中的第l个译码路径Pl，其返回的比特序列βv记为βv (l)，
Figure BDA0001219872570000021
为译码结果，则路径Pl的译码结果为：
Figure BDA0001219872570000031
其中，Gn为码长为n的极化码生成矩阵。如果
Figure BDA0001219872570000032
满足循环冗余校验，则将
Figure BDA0001219872570000033
选为最终译码结果。如果所有存活路径的译码结果
Figure BDA0001219872570000034
均不满足循环冗余校验，则选择路径度量值PM最小的路径所得到的
Figure BDA0001219872570000035
作为最终译码结果，整个译码过程结束。
优选的：所述采用基于串行抵消的译码方法如下：
译码时需要在节点之间传递待译符号的对数似然比LLR，对数似然比LLR初始值通过如下方式计算得到：
Figure BDA0001219872570000036
其中，L(ui)表示ui的对数似然比，ui表示发送比特，yi为译码器接收到的待译码的符号，Pr(yi|ui)表示发送符号为ui时接收到符号yi的概率。
对于译码树上的一个节点v，其包含的叶节点的个数为Nv，该节点从父节点接收一个长度为Nv的LLR值序列，记为av。译码开始时，由公式(1)计算接收符号的对数似然比并传递给根节点，然后由av计算要传递给左子节点的LLR值序列al：
al[i]＝sgn(av[i])sgn(av[i+Nv/2])min(|av[i]|，|av[i+Nv/2]|)，0≤i＜Nv/2(2)
其中，i是向量av的位置索引，av[i]表示av的第i个元素，
Figure BDA0001219872570000037
以递归的方式遍历左子节点，返回后从左子节点接收到长度为Nv/2的比特序列βl，
然后计算要传递给右子节点的LLR值序列ar：
ar[i]＝av[i+Nv/2]-(2βl[i]-1)av[i]，0≤i＜Nv/2 (3)
以递归的方式遍历右子节点，返回后从右字节点接收到长度为Nv的比特序列βr，然后计算要返回给父节点的比特序列βv：
Figure BDA0001219872570000041
其中，
Figure BDA0001219872570000042
为二进制异或运算。
优选的：所述步骤2中的译码策略如下：
如果是Rate-0节点，则：
βv[i]＝0，0≤i＜Nv (5)
如果是Rate-1节点，则：
βv[i]＝h(αv[i])，0≤i＜Nv (6)
其中，
Figure BDA0001219872570000043
如果是REP节点，则：
Figure BDA0001219872570000044
如果是SPC节点，则：
βv[i]＝h(αv[i])，0≤i＜Nv (9)
如果SPC节点的译码结果βv不满足奇偶校验，则找到序列αv中绝对值最小的，记为αv[min]，min是其在序列αv中的索引，然后将βv[min]中的比特位反转。
优选的：译码路径分裂门限的计算方法：
对于一个码长为N的极化码，第i个极化子信道记为
Figure BDA0001219872570000045
子信道
Figure BDA0001219872570000046
的截止速率为
Figure BDA0001219872570000047
通过如下方式递推计算得到：
Figure BDA0001219872570000048
Figure BDA0001219872570000049
其中，
Figure BDA0001219872570000051
K为信息位长度，N为码长，Eb/No为接收端信噪比的估计值。
码长为N的第i个比特，决定是否进行路径分裂的判决门限为
Figure BDA0001219872570000052
Figure BDA0001219872570000053
如果
Figure BDA0001219872570000054
则令
Figure BDA0001219872570000055
T为门限常数。
优选的：所述步骤4中根据节点类型确定不同的译码路径分裂策略，同时更新路径的PM值的方法：
当前节点深度为d，令N＝2d。当前信噪比估计值为Eb/No，按第3步所述方法找到Eb/No对应的列表待查。
如果是Rate-0节点，不产生新的译码路径，更新当前路径的度量值PM：
Figure BDA0001219872570000056
其中，h(·)为式(7)。
如果是Rate-1节点，原有路径的度量值PM不变。找到序列av中绝对值最小的两个，记为av[min 1]，av[min 2]，min 1、min 2为其分别在序列av中的位置索引，在该列表中查询
Figure BDA0001219872570000057
判断式(14)(15)是否成立：
Figure BDA0001219872570000058
Figure BDA0001219872570000059
若式(14)成立将当前路径复制一份新路径P1，复制时路径中的av、βv、PM都要复制，P1的度量值为PM1，将其βv[min 1]比特位反转，更新PM1：
PM1＝PM1-|αv[min 1]| (16)
若式(15)成立将当前路径复制一份新路径P2，将其βv[min 2]比特位反转，更新路径P2的度量值PM2：
PM2＝PM2-|αy[min 2]| (17)
若式(14)、式(15)同时成立，将当前路径复制一份新路径P3，将其βv[min 1]，βv[min 2]两个比特位反转，更新路径P3的度量值PM3：
PM3＝PM3-|av[min 1]|-|αv[min 2]| (18)
如果是REP节点，在该列表中查询门限值
Figure BDA0001219872570000061
如果
Figure BDA0001219872570000062
则不产生新路径，且路径的度量值PM不变。否则将当前路径复制一份新路径P1，将新路径中βv所有比特反转，新路径与原路径都按如下方式更新PM：
Figure BDA0001219872570000063
如果是SPC节点，若满足奇偶校验令q＝1，否则q＝0。更新原路径PM值：
PM＝PM-(1-q)|αv[min 1]| (20)
找到av中绝对值最小的4个，记为av[min 1]，av[min 2]，av[min 3]，av[min 4]，min1～min4分别为其在av中的位置索引，在列表中查询
Figure BDA0001219872570000064
判断式(21)-(23)是否成立：
Figure BDA0001219872570000065
Figure BDA0001219872570000066
Figure BDA0001219872570000067
若式(21)成立将当前路径复制一份新路径P1，将βv[min 1]，βv[min 2]比特位反转，更新路径P1的度量值PM1：
PM1＝PM1-q|αv[min 1]|-|αv[min 2]| (24)
若式(22)成立将当前路径复制一份新路径P2，将βv[min 1]，βv[min 3]比特位反转，更新路径P2的度量值PM2：
PM2＝PM2-q|αv[min 1]|-|αv[min 3]| (25)
若式(23)成立将当前路径复制一份新路径P3，将βv[min 1]，βv[min 4]比特位反转，更新路径P3的度量值PM3：
PM3＝PM3-q|αv[min 1]|-|αv[min 4]| (26)
对于列表中的每条译码路径，都需要独立进行译码和路径分裂，分裂后产生的新路径添加到同一列表中，如果列表中路径数量多于设定的最大值Lmax，则只保留PM值最大的Lmax个路径。
优选的：所述门限常数T大于等于10。
优选的：所述信噪比估计值Eb/No以dB为单位并四舍五入近似为整数。
本发明相比现有技术，具有以下有益效果：
本发明以每个比特判决时的置信度作为译码路径是否分裂的依据，设计了一整套适合各种类型译码树节点的译码方案。与现有技术相比，本发明具有如下增益效果：当置信度很高时，既可以保证误码率性能几乎不损失，又能够减少译码路径分裂次数，从而减少列表中译码路径的数量，降低了极化码译码器的计算复杂度，从而提高译码器的吞吐量；当信噪比越高时，本方法的性能增益越明显，可以极大提高译码器的吞吐量。
附图说明
图1为码长N＝8，信息位长度K＝4的一个极化码译码树实例；
图2为本发明设计的极化码译码方法完整流程。
具体实施方式
下面结合附图和具体实施例，进一步阐明本发明，应理解这些实例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明提出的基于串行抵消列表的减少译码路径分裂的极化码译码方法，以终端接收设备为例来给出一种译码器实施例：
如图1所述，包括如下步骤：
第一步：递归遍历译码树。
本发明采用树形译码结构，对于一个码长为N，信息位长度为K的极化码，其译码树为深度为n＝log2N的二叉树。图1给出了N＝8、K＝4一个译码树实例，图中白色叶节点表示休眠位比特，黑色叶节点表示信息位比特。基于串行抵消的译码方法从根节点开始递归遍历整个译码树。
译码时需要在节点之间传递待译符号的对数似然比(LLR，log-likelihoodratio)，LLR初始值通过如下方式计算得到：
Figure BDA0001219872570000081
其中，ui表示发送比特，L(ui)表示ui的对数似然比，yi为译码器接收到的待译码的符号，Pr(yi|ui)表示发送符号为ui时，接收到符号yi的概率。Pr(yi|ui)可从解调器的软判决信息中得到。
如图1所示，对于译码树上的一个节点v，其包含的叶节点的个数为Nv，该节点从父节点接收一个长度为Nv的LLR值序列，记为av。译码开始时，由公式(1)计算接收符号的对数似然比并传递给根节点，然后由av计算要传递给左子节点的LLR值序列al：
al[i]＝sgn(av[i])sgn(av[i+Nv/2])min(|av[i]|，|av[i+Nv/2]|)，0≤i＜Nv/2(2)
其中，i是向量av的位置索引，av[i]表示av的第i个元素，
Figure BDA0001219872570000082
以递归的方式遍历左子节点，返回后从左子节点接收到长度为Nv/2的比特序列βl，然后计算要传递给右子节点的LLR值序列ar：
ar[i]＝av[i+Nv/2]-(2βl[i]-1)av[i]，0≤i＜Nv/2 (3)
以递归的方式遍历右子节点，返回后从右字节点接收到长度为Nv的比特序列βr，然后计算要返回给父节点的比特序列βv：
Figure BDA0001219872570000091
其中，
Figure BDA0001219872570000092
为二进制异或运算。
在本步骤中，左子节点和右子节点以同样的方式向其子节点遍历。如果遇到以下几种类型的节点，则不再向下递归而直接译码后返回。这几种节点类型分别是：若其包含的叶节点全是休眠位比特，则记为Rate-0节点；若其包含的叶节点全是信息位比特，则记为Rate-1节点；若其包含的叶节点只有最后一位是信息位比特，则记为REP(RepetitionNodes)节点；若其包含的叶节点只有第一位是休眠位比特，则记为SPC(Single ParityCheck Nodes)节点。
第二步：针对第一步中提到的4种类型的节点分别给出相应的译码策略。
本步骤中Nv、av、βv的定义同第一步。
如果是Rate-0节点，则：
βv[i]＝0，0≤i＜Nv (5)
如果是Rate-1节点，则：
βv[i]＝h(αv[i])，0≤i＜Nv (6)
其中，
Figure BDA0001219872570000093
如果是REP节点，则：
Figure BDA0001219872570000094
如果是SPC节点，则：
βv[i]＝h(αv[i])，0≤i＜Nv (9)
如果SPC节点的译码结果βv不满足奇偶校验，则找到序列αv中绝对值最小的，记为αv[min]，min是其在序列αv中的索引，然后将βv[min]中的比特位反转。
第三步：计算译码路径分裂门限。
在上一步中每次进行译码时，接收到的LLR值的绝对值的大小可以看作该比特译码置信度的高低，因此对该节点从父节点接收到LLR值序列av进行门限判别，以决定是否分裂新的译码路径。首先对决定是否进行路径分裂的LLR门限值进行计算：
对于一个码长为N的极化码，第i个极化子信道记为
Figure BDA0001219872570000101
于信道
Figure BDA0001219872570000102
的截止速率为
Figure BDA0001219872570000103
可以看做是该子信道极化程度的度量标准，因此我们用
Figure BDA0001219872570000104
计算出一个与信道优劣有关的动态门限。
Figure BDA0001219872570000105
可以通过如下方式递推计算得到：
Figure BDA0001219872570000106
Figure BDA0001219872570000107
其中，
Figure BDA0001219872570000108
K为信息位长度，N为码块长度，Eb/No为接收端信噪比的估计值，可由常用的信道估计算法得到。
码长为N的第i个比特，决定是否进行路径分裂的判决门限为
Figure BDA0001219872570000109
Figure BDA00012198725700001010
如果
Figure BDA00012198725700001011
则令
Figure BDA00012198725700001012
(T为门限常数，运营商可以根据网络运行情况自行确定，推荐取值10以上)。当某个子信道极化程度越高时，
Figure BDA00012198725700001013
就越小，门限值
Figure BDA00012198725700001014
越大，因此
Figure BDA00012198725700001015
是一个与子信道优劣有关的动态门限。
译码时先估计信噪比Eb/No(以dB为单位并四舍五入近似为整数)，当前Eb/No对应的门限值计算过一次之后就将
Figure BDA00012198725700001016
存入一张表格中。下一个码块译码时，如果Eb/No对应的门限值已存在则直接找到对应的表格待查，如果不存在则重新计算存到对应的表格中。当有足够多的码块被译码后，系统工作范围内的所有Eb/No对应的门限都会存在于表格中，因而不会增加译码的实时复杂度。
第四步：译码路径分裂策略。
译码开始时，列表里只有一条译码路径，每条译码路径都有一个评价路径优劣的度量值PM(Path Metric)，初始路径的PM值设为0。第二步在译码时根据节点类型选择不同的译码路径分裂策略，要同时更新路径的PM值。
本步骤中Nv、av、βv的定义同第一步，当前节点深度为d，令N＝2d。当前信噪比估计值为Eb/No，按第三步所述方法找到Eb/No对应的表格待查。
如果是Rate-0节点，不产生新的译码路径，更新当前路径的度量值PM：
Figure BDA0001219872570000111
其中，h(·)为式(7)。
如果是Rate-1节点，原有路径的度量值PM不变；找到序列av中绝对值最小的两个，记为av[min 1]，av[min 2]，min 1、min 2为其分别在序列av中的位置索引，在该表格中查询
Figure BDA0001219872570000112
判断式(14)(15)是否成立：
Figure BDA0001219872570000113
Figure BDA0001219872570000114
若(14)成立将当前路径复制一份新路径P1(复制时路径中的av、βv、PM都要复制)，P1的度量值为PM1，将其βv[min 1]比特位反转，更新PM1：
PM1＝PM1-|αv[min 1]| (16)
若(15)成立将当前路径复制一份新路径P2，将其βv[min 2]比特位反转，更新路径P2的度量值PM2：
PM2＝PM2-|αv[min 2]| (17)
若(14)(15)同时成立，将当前路径复制一份新路径P3，将其βv[min 1]，βv[min 2]两个比特位反转，更新路径P3的度量值PM3：
PM3＝PM3-|αv[min 1]|-|αv[min 2]| (18)
如果是REP节点，在该表格中查询门限值
Figure BDA0001219872570000121
如果
Figure BDA0001219872570000122
则不产生新路径，且路径的度量值PM不变；否则将当前路径复制一份新路径P1，将新路径中βv所有比特反转，新路径与原路径都按如下方式更新PM：
Figure BDA0001219872570000123
如果是SPC节点，若满足奇偶校验令q＝1，否则q＝0；更新原路径PM值：
PM＝PM-(1-q)|αv[min 1]| (20)
找到av中绝对值最小的4个，记为av[min 1]，av[min 2]，av[min 3]，av[min 4]，min1～min4分别为其在av中的位置索引，在表格中查询
Figure BDA0001219872570000124
判断式(21)-(23)是否成立：
Figure BDA0001219872570000125
Figure BDA0001219872570000126
Figure BDA0001219872570000127
若(21)成立将当前路径复制一份新路径P1，将βv[min 1]，βv[min 2]比特位反转，更新路径P1的度量值PM1：
PM1＝PM1-q|αv[min 1]|-|αv[min 2]| (24)
若(22)成立将当前路径复制一份新路径P2，将βv[min 1]，βv[min 3]比特位反转，更新路径P2的度量值PM2：
PM2＝PM2-q|αv[min 1]|-|αv[min 3]| (25)
若(23)成立将当前路径复制一份新路径P3，将βv[min 1]，βv[min 4]比特位反转，更新路径P3的度量值PM3：
PM3＝PM3-q|αv[min 1]|-|αv[min 4]| (26)
对于列表中的每条译码路径，都需要独立进行译码和路径分裂，分裂后产生的新路径添加到同一列表中，如果列表中路径数量多于设定的最大值Lmax(Lmax取值为常整数，该取值可由运营商根据误码率要求情况自行确定，Lmax越大，误码率越低)，则只保留PM值最大的Lmax个路径。
第五步：根据循环冗余校验选择最终译码结果。
如果译码树从根节点返回，对于列表中的第l个译码路径Pl，其返回的比特序列βv记为βv (l)，
Figure BDA0001219872570000131
为译码结果，则路径Pl的译码结果为：
Figure BDA0001219872570000132
其中，Gn为码长为n的极化码生成矩阵。如果
Figure BDA0001219872570000133
满足循环冗余校验，则将
Figure BDA0001219872570000134
选为最终译码结果。如果所有存活路径的译码结果
Figure BDA0001219872570000135
均不满足循环冗余校验，则选择路径度量值PM最小的路径所得到的
Figure BDA0001219872570000136
作为最终译码结果，整个译码过程结束。
如图2所示，为本发明的流程图，包括以下步骤：
S1，计算LLR初始值并赋值给根节点，进入步骤S2。
S2，判断节点是否为需要译码的节点，若不是需要译码的节点，进入步骤S3；若是需要译码的节点，进入步骤S4。
S3，计算左子节点LLR值，计算出的左子节点LLR值赋值给当前节点的左子节点，进入步骤S2。
S4，根据节点类型进行译码和路径分裂，进入步骤S5。
S5，返回上一层节点，并判断上一层节点是否为根节点，若不是根节点，进入步骤S6；若是根节点，进入步骤S7。
S6，判断节点从哪个子节点返回，若是从右子节点返回，合并左右节点返回的比特序列，进入步骤S5；若是从左子节点返回，计算右子节点LLR值，计算出的右子节点LLR值赋值给当前节点的右子节点，进入步骤S2。
S7，根据循环冗余校验从待选路径选择最终译码结果。
本发明以串行抵消列表方法为基础，提出了一种提高译码器吞吐量并且保证误码率性能的极化码译码方法。该方法考虑了译码时的置信度水平，在比特位置信度很高的情况下可以不分裂新的译码路径，又不会损失误码率性能。该方法能够自适应信道噪声条件，当信噪比越高时，译码器性能增益越明。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (3)

1.一种减少译码路径分裂的极化码译码方法，其特征在于，包括以下步骤：
步骤1，递归遍历译码树：对于一个码长为N,信息位长度为K的极化码，其译码树为深度为n＝log2N的二叉树，采用基于串行抵消的译码方法从根节点开始递归遍历整个译码树；其中，左子节点和右子节点以同样的方式向其子节点遍历，如果遇到以下4种类型的节点，则不再向下递归而直接译码后返回；这4种节点类型分别是：若其包含的叶节点全是休眠位比特，则记为Rate-0节点；若其包含的叶节点全是信息位比特，则记为Rate-1节点；若其包含的叶节点只有最后一位是信息位比特，则记为REP节点；若其包含的叶节点只有第一位是休眠位比特，则记为SPC节点；
采用基于串行抵消的译码方法如下：
译码时需要在节点之间传递待译符号的对数似然比LLR，对数似然比LLR初始值通过如下方式计算得到：
Figure FDA0002326117030000011
其中，L(ui)表示ui的对数似然比，ui表示发送比特，yi为译码器接收到的待译码的符号，Pr(yi|ui)表示发送符号为ui时接收到符号yi的概率；
对于译码树上的一个节点v，其包含的叶节点的个数为Nv，该节点从父节点接收一个长度为Nv的LLR值序列，记为av；译码开始时，由公式(1)计算接收符号的对数似然比并传递给根节点，然后由av计算要传递给左子节点的LLR值序列al：
al[i]＝sgn(av[i])sgn(av[i+Nv/2])min(|av[i]|,|av[i+Nv/2]|)，0≤i＜Nv/2 (2)
其中，i是向量av的位置索引，av[i]表示av的第i个元素，
Figure FDA0002326117030000012
以递归的方式遍历左子节点，返回后从左子节点接收到长度为Nv/2的比特序列βl，然后计算要传递给右子节点的LLR值序列ar：
ar[i]＝av[i+Nv/2]-(2βl[i]-1)av[i]，0≤i＜Nv/2 (3)
以递归的方式遍历右子节点，返回后从右字节点接收到长度为Nv的比特序列βr，然后计算要返回给父节点的比特序列βv：
Figure FDA0002326117030000021
其中，
Figure FDA0002326117030000022
为二进制异或运算；
步骤2，对步骤1中提到的4种类型的节点根据节点类型分别给出相应的译码策略，得到各节点相应的译码结果；
译码策略如下：
如果是Rate-0节点，则：
βv[i]＝0,0≤i＜Nv (5)
如果是Rate-1节点，则：
βv[i]＝h(αv[i])，0≤i＜Nv (6)
其中，
Figure FDA0002326117030000023
如果是REP节点，则：
Figure FDA0002326117030000024
如果是SPC节点，则：
βv[i]＝h(αv[i])，0≤i＜Nv (9)
如果SPC节点的译码结果βv不满足奇偶校验，则找到序列αv中绝对值最小的，记为αv[min]，min是其在序列αv中的索引，然后将βv[min]中的比特位反转；
步骤3，根据译码时节点之间传递待译符号的对数似然比计算译码路径分裂门限，根据译码路径分裂门限决定是否分裂新的译码路径；译码时先估计信噪比Eb/No，当前Eb/No对应的译码路径分裂门限计算过一次之后就将译码路径分裂门限存入一张列表中；下一个码块译码时，如果Eb/No对应的译码路径分裂门限已存在则直接找到对应的列表待查，如果不存在则重新计算存到对应的列表中；
译码路径分裂门限的计算方法：
对于一个码长为N的极化码，第i个极化子信道记为
Figure FDA0002326117030000031
子信道
Figure FDA0002326117030000032
的截止速率为
Figure FDA0002326117030000033
通过如下方式递推计算得到：
Figure FDA0002326117030000034
Figure FDA0002326117030000035
其中，
Figure FDA0002326117030000036
K为信息位长度，N为码长，Eb/No为接收端信噪比的估计值；
码长为N的第i个比特，决定是否进行路径分裂的判决门限为
Figure FDA0002326117030000037
Figure FDA0002326117030000038
如果
Figure FDA0002326117030000039
则令
Figure FDA00023261170300000310
T为门限常数；
步骤4，译码路径分裂策略：译码开始时，列表里只有一条译码路径，每条译码路径都有一个评价路径优劣的度量值PM，初始路径的PM值设为0；第2步在译码时，根据当前信噪比估计值Eb/No，结合步骤3列表中的分裂门限，根据步骤1所述的节点类型确定不同的译码路径分裂策略，同时更新路径的PM值；分裂后产生的新路径添加到同一列表中，如果列表中路径数量多于设定的最大值Lmax，则只保留PM值最大的Lmax个路径；
根据节点类型确定不同的译码路径分裂策略，同时更新路径的PM值的方法：
当前节点深度为d，令N＝2d；当前信噪比估计值为Eb/No，按第3步所述方法找到Eb/No对应的列表待查；
如果是Rate-0节点，不产生新的译码路径，更新当前路径的度量值PM：
Figure FDA00023261170300000311
其中，h(·)为式(7)；
如果是Rate-1节点，原有路径的度量值PM不变；找到序列av中绝对值最小的两个，记为av[min1],av[min2]，min1、min2为其分别在序列av中的位置索引，在该列表中查询
Figure FDA0002326117030000041
判断式(14)(15)是否成立：
Figure FDA0002326117030000042
Figure FDA0002326117030000043
若式(14)成立将当前路径复制一份新路径P1，复制时路径中的av、βv、PM都要复制，P1的度量值为PM1，将其βv[min1]比特位反转，更新PM1：
PM1＝PM1-|αv[min1]| (16)
若式(15)成立将当前路径复制一份新路径P2，将其βv[min2]比特位反转，更新路径P2的度量值PM2：
PM2＝PM2-|αv[min2]| (17)
若式(14)、式(15)同时成立，将当前路径复制一份新路径P3，将其βv[min1]，βv[min2]两个比特位反转，更新路径P3的度量值PM3：
PM3＝PM3-|αv[min1]|-|αv[min2]| (18)
如果是REP节点，在该列表中查询门限值
Figure FDA0002326117030000044
如果
Figure FDA0002326117030000045
则不产生新路径，且路径的度量值PM不变；否则将当前路径复制一份新路径P1，将新路径中βv所有比特反转，新路径与原路径都按如下方式更新PM：
Figure FDA0002326117030000046
如果是SPC节点，若满足奇偶校验令q＝1，否则q＝0；更新原路径PM值：
PM＝PM-(1-q)|αv[min1]| (20)
找到av中绝对值最小的4个，记为av[min1],av[min2],av[min3],av[min4]，min1～min4分别为其在av中的位置索引，在列表中查询
Figure FDA0002326117030000047
判断式(21)-(23)是否成立：
Figure FDA0002326117030000048
Figure FDA0002326117030000049
Figure FDA0002326117030000051
若式(21)成立将当前路径复制一份新路径P1，将βv[min1]，βv[min2]比特位反转，更新路径P1的度量值PM1：
PM1＝PM1-q|αv[min1]|-|αv[min2]| (24)
若式(22)成立将当前路径复制一份新路径P2，将βv[min1]，βv[min3]比特位反转，更新路径P2的度量值PM2：
PM2＝PM2-q|αv[min1]|-|αv[min3]| (25)
若式(23)成立将当前路径复制一份新路径P3，将βv[min1]，βv[min4]比特位反转，更新路径P3的度量值PM3：
PM3＝PM3-q|αv[min1]|-|αv[min4]| (26)
对于列表中的每条译码路径，都需要独立进行译码和路径分裂，分裂后产生的新路径添加到同一列表中，如果列表中路径数量多于设定的最大值Lmax，则只保留PM值最大的Lmax个路径；
步骤5，根据循环冗余校验选择最终译码结果：
如果译码树从根节点返回，对于列表中的第l个译码路径Pl，其返回的比特序列βv记为βv (l)，
Figure FDA0002326117030000056
为译码结果，则路径Pl的译码结果为：
Figure FDA0002326117030000052
其中，Gn为码长为n的极化码生成矩阵；如果
Figure FDA0002326117030000053
满足循环冗余校验，则将
Figure FDA0002326117030000054
选为最终译码结果；如果所有存活路径的译码结果
Figure FDA0002326117030000055
均不满足循环冗余校验，则选择路径度量值PM最小的路径所得到的
Figure FDA0002326117030000057
作为最终译码结果，整个译码过程结束。
2.根据权利要求1所述减少译码路径分裂的极化码译码方法，其特征在于：所述门限常数T大于等于10。
3.根据权利要求1所述减少译码路径分裂的极化码译码方法，其特征在于：所述信噪比估计值Eb/No以dB为单位并四舍五入近似为整数。
CN201710062953.4A 2017-02-01 2017-02-01 一种减少译码路径分裂的极化码译码方法 Active CN106877884B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710062953.4A CN106877884B (zh) 2017-02-01 2017-02-01 一种减少译码路径分裂的极化码译码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710062953.4A CN106877884B (zh) 2017-02-01 2017-02-01 一种减少译码路径分裂的极化码译码方法

Publications (2)

Publication Number Publication Date
CN106877884A CN106877884A (zh) 2017-06-20
CN106877884B true CN106877884B (zh) 2020-04-28

Family

ID=59165866

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710062953.4A Active CN106877884B (zh) 2017-02-01 2017-02-01 一种减少译码路径分裂的极化码译码方法

Country Status (1)

Country Link
CN (1) CN106877884B (zh)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107666325B (zh) * 2017-09-20 2021-07-30 南京大学 基于列表连续删除算法的极化码译码路径选择方法
CN109547034B (zh) * 2017-09-22 2021-10-01 华为技术有限公司 译码方法及设备、译码器
CN107896137B (zh) * 2017-11-03 2020-10-27 东南大学 一种适用于极化码译码路径分裂的排序方法
CN109802690B (zh) * 2017-11-17 2022-04-22 中兴通讯股份有限公司 译码方法、装置和计算机可读存储介质
CN108063649B (zh) * 2017-12-14 2020-10-02 东南大学 一种低时延低复杂度的极化码译码方法
CN110324111B (zh) * 2018-03-31 2020-10-09 华为技术有限公司 一种译码方法及设备
CN110391817B (zh) * 2018-04-17 2021-02-09 华为技术有限公司 极化码的译码方法和装置
CN108566213B (zh) * 2018-04-25 2022-05-27 东南大学 极化码的串行抵消列表比特翻转译码方法
CN110830051B (zh) * 2018-08-07 2023-06-23 普天信息技术有限公司 一种译码方法和装置
WO2020042089A1 (zh) * 2018-08-30 2020-03-05 华为技术有限公司 Scl并行译码方法、装置及设备
CN111200439B (zh) 2018-11-16 2022-05-06 华为技术有限公司 译码方法、装置及设备
US11777533B2 (en) 2018-11-27 2023-10-03 Telefonaktiebolagget LM Ericsson (Publ) Method for polar decoding with dynamic successive cancellation list size and polar decoder
CN111404558B (zh) * 2019-01-03 2023-06-27 大唐移动通信设备有限公司 一种Polar码译码的方法、译码器及计算机存储介质
CN109951190B (zh) * 2019-03-15 2020-10-13 北京科技大学 一种自适应Polar码SCL译码方法及译码装置
CN109995381B (zh) * 2019-04-09 2022-09-13 合肥工业大学 一种基于预判机制的极化码译码系统及其方法
CN110190857B (zh) * 2019-04-28 2023-04-11 深圳大学 一种crc辅助校验极化码译码方法和智能终端
CN110535477A (zh) * 2019-08-30 2019-12-03 山东科技大学 并行的极化码译码方法
CN111181573B (zh) * 2020-03-09 2023-08-18 北京华力创通科技股份有限公司 数据译码方法、装置及电子设备
CN111786683B (zh) * 2020-05-29 2022-11-08 东南大学 一种低复杂度的极化码多码块译码器
CN111654292B (zh) * 2020-07-20 2023-06-02 中国计量大学 一种基于动态阈值的分裂简化极化码连续消除列表译码器
CN114268327A (zh) * 2021-11-26 2022-04-01 北京中科晶上科技股份有限公司 极化码译码方法、装置、存储介质和处理器
US20230177232A1 (en) * 2021-12-08 2023-06-08 International Business Machines Corporation Circuit serialization for parameterized circuit simulation
US11791844B2 (en) 2021-12-22 2023-10-17 Huawei Technologies Co., Ltd. Systems and methods for using special nodes for polar encoding in polar codes

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102694625A (zh) * 2012-06-15 2012-09-26 北京邮电大学 一种循环冗余校验辅助的极化码译码方法
CN105227189A (zh) * 2015-09-24 2016-01-06 电子科技大学 分段crc辅助的极化码编译码方法
CN105281785A (zh) * 2015-10-22 2016-01-27 东南大学 一种列表连续消除极化码译码方法、装置
CN105515590A (zh) * 2015-12-09 2016-04-20 东南大学 一种基于随机二进制数据流的有效低复杂度串行抵消列表极化码译码算法及其译码构架

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102694625A (zh) * 2012-06-15 2012-09-26 北京邮电大学 一种循环冗余校验辅助的极化码译码方法
CN105227189A (zh) * 2015-09-24 2016-01-06 电子科技大学 分段crc辅助的极化码编译码方法
CN105281785A (zh) * 2015-10-22 2016-01-27 东南大学 一种列表连续消除极化码译码方法、装置
CN105515590A (zh) * 2015-12-09 2016-04-20 东南大学 一种基于随机二进制数据流的有效低复杂度串行抵消列表极化码译码算法及其译码构架

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Efficient Soft Cancelation Decoder Architectures for Polar Codes;Jun Lin 等;《IEEE Transactions on Very Large Scale Integration (VLSI) Systems》;20170131;第25卷(第1期);第87-99页 *
Multi-Mode Unrolled Architectures for Polar Decoders;Pascal Giard 等;《IEEE Transactions on Circuits and Systems I: Regular Papers》;20160930;第63卷(第9期);第1443-1453页 *

Also Published As

Publication number Publication date
CN106877884A (zh) 2017-06-20

Similar Documents

Publication Publication Date Title
CN106877884B (zh) 一种减少译码路径分裂的极化码译码方法
CN107248866B (zh) 一种降低极化码译码时延的方法
CN108462558B (zh) 一种极化码scl译码方法、装置及电子设备
CN109660264B (zh) 高性能极化码译码算法
TWI401892B (zh) 在使用低密度同位檢查碼的通訊系統中的通道編碼與解碼的方法及其裝置
CN109286405B (zh) 一种低复杂度的极化码递进式比特翻转sc译码方法
CN109921804B (zh) 一种自适应融合串行抵消列表极化码译码方法及系统
CN110995278B (zh) 一种改进极性码串行消除列表比特翻转译码方法及系统
US8000411B2 (en) Decoding scheme using multiple hypotheses about transmitted messages
CN111277277B (zh) 一种减少极化码连续对消表译码算法译码延迟的方法及装置
JP2006523985A (ja) ターボデコーダの反復の早期終了
JP2007515122A (ja) 移動通信システムで不均一誤り確率を有する符号化によるデータ送受信方法及び装置
CN111726202B (zh) 一种极化码置信传播译码的提前终止迭代方法
CN114285418A (zh) 一种改进的极化码epc-ms-scl译码方法
JP3756525B2 (ja) 固定長デシジョンウィンドウを使用するデータ信号の復号方法
CN110190857B (zh) 一种crc辅助校验极化码译码方法和智能终端
CN110943745A (zh) 一种提前终止迭代输出结果的极化码bp译码方法及系统
US7975212B2 (en) Sequential decoding method and apparatus thereof
US7500173B2 (en) Method of decoding a data word
CN113556135B (zh) 基于冻结翻转列表的极化码置信传播比特翻转译码方法
CN114421975A (zh) 一种基于翻转集的极化码sclf译码方法
CN114696953A (zh) 一种自由空间光通信的信道编码和译码方法
JP3285090B2 (ja) デジタルデータ復号器
CN110212924B (zh) 一种lt码编解码方法及系统
US20070283232A1 (en) Method for near maximum-likelihood sequential decoding

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
TR01 Transfer of patent right

Effective date of registration: 20210331

Address after: 201306 building C, No. 888, Huanhu West 2nd Road, Lingang New Area, Pudong New Area, Shanghai

Patentee after: Shanghai Hanxin Industrial Development Partnership (L.P.)

Address before: 211189 No. 2 Southeast University Road, Jiangning District, Nanjing, Jiangsu

Patentee before: SOUTHEAST University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230919

Address after: 201615 room 301-6, building 6, no.1158, Jiuting Central Road, Jiuting Town, Songjiang District, Shanghai

Patentee after: White box (Shanghai) Microelectronics Technology Co.,Ltd.

Address before: 201306 building C, No. 888, Huanhu West 2nd Road, Lingang New Area, Pudong New Area, Shanghai

Patentee before: Shanghai Hanxin Industrial Development Partnership (L.P.)

TR01 Transfer of patent right