CN105720992A - 一种极化码的简化译码方法 - Google Patents

一种极化码的简化译码方法 Download PDF

Info

Publication number
CN105720992A
CN105720992A CN201610045755.2A CN201610045755A CN105720992A CN 105720992 A CN105720992 A CN 105720992A CN 201610045755 A CN201610045755 A CN 201610045755A CN 105720992 A CN105720992 A CN 105720992A
Authority
CN
China
Prior art keywords
node
layer
information
decoding
fixed value
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
Application number
CN201610045755.2A
Other languages
English (en)
Other versions
CN105720992B (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.)
Shenzhen Graduate School Harbin Institute of Technology
Original Assignee
Shenzhen Graduate School Harbin Institute of 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 Shenzhen Graduate School Harbin Institute of Technology filed Critical Shenzhen Graduate School Harbin Institute of Technology
Priority to CN201610045755.2A priority Critical patent/CN105720992B/zh
Publication of CN105720992A publication Critical patent/CN105720992A/zh
Application granted granted Critical
Publication of CN105720992B publication Critical patent/CN105720992B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/13Linear codes
    • 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/1191Codes on graphs other than LDPC codes

Landscapes

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

Abstract

本发明提供一种极化码的简化译码方法,包括以下步骤:步骤S1,对接收到的极化码建立树图模型;步骤S2,将树图中各个节点分为信息节点、固定值节点和混合节点;步骤S3,对树图进行剪枝,剪去固定值节点和信息节点之间的子树,留下与混合节点相连接的固定值节点和信息节点;步骤S4,对剪枝后的译码树图进行译码运算。本发明能够在确保不对误码性能产生影响的前提下,达到减小译码复杂度和提高译码效率的目的,最后进行了仿真验证,简化后SC译码和BP译码的译码复杂度较原始算法分别降低了36%~65%和41%~67%。

Description

一种极化码的简化译码方法
技术领域
本发明涉及一种极化码的译码方法,尤其涉及一种极化码的简化译码方法。
背景技术
2009年Arlkan提出了信道极化,对于二元离散无记忆信道,通过组合与分离可以使信道容量一部分趋于0,另一部分趋于1,从而形成了信道极化。极化码是一种基于信道极化提出的新型信道编码,根据所需码率在对应信道容量极化于1的部分编入信息,而其他部分编入固定值,能在二元对称信道下逼近香农限。因此,极化码的概念提出以来,在国际上引起了高度关注。目前,极化码的主要研究集中在实际应用中的信道选择和译码器的设计。
极化码的原始译码方法有连续删除(SC)译码和置信传播(BP)译码,SC译码是Arlkan提出的一种针对极化码的译码方法,对信源信号逐位进行估计后删除其冗余的关联信息,并将估计值作为先验信息代入之后的译码运算,码长为N的极化码SC译码复杂度为O(NlogN)。BP译码是一种被广泛应用的迭代译码方法,Arlkan根据对极化码二分图模型的分析,提出了极化码的BP译码算法。BP译码性能较SC译码有一定提高,但码长为N的极化码BP译码复杂度达到了O(tNlogN),t为迭代次数。
极化码原始编译码方法在信道选择和译码器性能上有一定缺陷,使极化码的译码性能较低密度奇偶校验码(LDPC)和Turbo码有一定的差距。为了提高极化码译码性能,基于原始的SC译码方法,Tal等人提出了一种序列SC译码(SCL),提升了误码性能,特别是在CRC方法的辅助下,SCL译码的误码性能可以超过最大似然译码,但同时译码计算复杂度也随之提高至O(lNlogN),l为搜索宽度。Niu等人提出了一种堆栈构造的SC译码(SCS),其误码性能也较原始SC译码有所提高,但其计算复杂度也同样提高为O(lNlogN),l为搜索宽度,实际的算法执行中,高信噪比下搜索宽度会降低,使得SCS译码实际复杂度较SCL译码略低。可见,虽然极化码能够逼近信道香农限,但代价是译码复杂度的提高。
发明内容
本发明所要解决的技术问题是需要提供一种在不降低译码性能的前提下,提高极化码译码效率的极化码的简化译码方法。
对此,本发明提供一种极化码的简化译码方法,包括以下步骤:
步骤S1,对接收到的极化码建立树图模型;
步骤S2,将树图中各个节点分为信息节点、固定值节点和混合节点;
步骤S3,对树图进行剪枝,剪去固定值节点和信息节点之间的子树,留下与混合节点相连接的固定值节点和信息节点;
步骤S4,对剪枝后的译码树图进行译码运算。
本发明的进一步改进在于,所述步骤S1中,对接收到的码长为N=2n的极化码,将其二分图模型表示为一个n+1层的树图形式,该树图中共有2n+1-1个节点,其中处于第n层的各变量节点对应信源比特,而处于第0层的根节点对应编码节点;其中,n为大于1的自然数。
本发明的进一步改进在于,所述步骤S2中,对于第i层的任一节点,如果该层包含的变量节点全部对应固定值比特,则称其为固定值节点;如果该层包含的变量节点全部对应信息比特,则称其为信息节点;如果该层包含的变量节点对应的既有固定值比特又有信息比特,则称其为混合节点。
本发明的进一步改进在于,所述步骤S2中,对于树图第i层左数第m个的节点vim,m∈{1,…,2i},通过进行计算,得到各个节点的分类,其中ceil是向上取整函数,count是重复数计数函数,即对集合中值为m的项进行计数,为构造极化码时挑选的K维信息比特域集;当进而能够对树图每一层的节点进行分类,i为大于1的自然数。
本发明的进一步改进在于,所述步骤S4中,将接收到的信道观测软信息对数似然比L输入根节点,然后在剪枝后的SC译码树图中依次按照访问根节点、遍历左子树和遍历右子树的顺序从根节点开始进行递推译码运算。
本发明的进一步改进在于,所述步骤S4中,在遍历到第i层任一固定值节点vi,f时,i∈1,…,n,根据公式直接译出vi,f包含的各变量节点vn,f[j]所对应的固定值比特un,f[j],其中,j=1,...,2n-i;并根据计算vi,f对应的辅助信息Ui,f并送至vi,f的父节点,其中,vi,f为固定值节点,un,f[j]为第i层任一固定值节点vi,f包含的第n层各变量节点对应的固定值比特,为固定值比特的估计值,是极化码2n-i阶生成矩阵;在遍历到第i层任一信息节点vi,e时,i∈1,…,n,对信息节点vi,e获得的软信息Li,e进行判决计算得到信息节点vi,e对应的辅助信息Ui,e=g(Li,e),使用直接对信息节点vi,e包含的变量节点对应的信息比特un,e[j],j=1,…,2n-i进行估计,并将辅助信息Ui,e送至信息节点vi,e的父节点,其中,un,e[j]为第i层任一信息节点vi,e包含的第n层各变量节点对应的信息比特,为信息比特的估计值;在遍历到第i层任一混合节点vi,m时,节点vi通过接收到i-1层传递的软信息对数似然比Li,v,计算并向处于第i+1层的其左子节点vi+1,l和右子节点vi+1,r所分别传递的软信息Li+1,l和软信息Li+1,r节点vi通过接收到i+1层传递的辅助信息Ui+1,l和Ui+1,r,计算并向处于第i-1层的其父节点vi-1,p传递的辅助信息Ui,v U i , v ( 2 k - 1 ) = U i + 1 , l ( k ) ⊕ U i + 1 , r ( k ) U i , v ( 2 k ) = U i + 1 , r ( k ) , k = 1 , ... , 2 n - i - 1 , 其中,为实数运算符号,为模二和运算符号;处于第n层的任一变量节点vn,j向上提供的辅助信息为其中,是变量节点vn,j对应的信源比特估计值。
本发明的进一步改进在于,步骤S2中,将树图中各个节点分为信息节点、固定值节点和混合节点;对树图进行首次剪枝,在译码运算开始前,将初始值按公式输入根节点;其中,是节点(n+1,j)在第0次迭代过程中向左传递的软信息对数似然比,是原始独立信道道的信道转移概率。
本发明的进一步改进在于,所述步骤S2中,对剪枝后树图上的各混合节点对应的软信息β值计算,按照第n层到第0层层序遍历的顺序进行译码,在遍历到第i层任一混合节点vi,m时,若其左子节点vi+1,l是固定值节点,其右子节点vi+1,r是信息节点,则若其左子节点vi+1,l是固定值节点,其右子节点vi+1,r是混合节点,则 β i , v ( 2 k - 1 ) = α i , v ( 2 k ) + β i + 1 , r ( k ) β i , v ( 2 k ) = α i , v ( 2 k - 1 ) + β i + 1 , r ( k ) , k = 1 , . . . , 2 n - i - 1 ; 若其左子节点vi+1,l是混合节点,其右子节点vi+1,r是信息节点,则若其左子节点vi+1,l和右子节点vi+1,r均是混合节点,则其中,节点vi通过接收到i+1层传递的软信息βi+1,l和软信息βi+1,r,计算并向处于第i-1层的其父节点vi-1,p传递其软信息βi,v;节点vi通过接收到i-1层传递的软信息对数似然比αi,v,计算并向处于第i+1层的其左子节点vi+1,l和其右子节点vi+1,r所分别传递的软信息αi+1,l和软信息αi+1,r
本发明的进一步改进在于,所述步骤S2中,完成各混合节点对应的软信息βi,v计算后,按照第0层到第n层层序遍历的顺序进行各混合节点对应的软信息α值计算,若vi+1,l是固定值节点,vi+1,r是混合节点,则αi+1,r(k)=αi,v(2k-1)+αi,v(2k),k=1,…,2n-i-1;若vi+1,l是混合节点,vi+1,r是信息节点,则若vi+1,l和vi+1,r均是混合节点,则然后将计算结果依次代入下一次迭代过程中,直至前t-1次迭代过程完成,保存第t-1次迭代后的各节点对应的软信息α值结果,然后重新建立极化码的树图模型。
本发明的进一步改进在于,
所述步骤S4中,进行译码运算的第t次迭代过程,计算剪枝后树图上的各混合节点对应的软信息β值,根据已保存的第t-1次迭代后的α值;
首先,在遍历到第i层任一混合节点vi,m时,若其左子节点vi+1,l是固定值节点,其右子节点vi+1,r是信息节点,则 β i , v ( 2 k - 1 ) = α i , v ( 2 k ) β i , v ( 2 k ) = α i , v ( 2 k - 1 ) , k = 1 , ... , 2 n - i - 1 ; 若其左子节点vi+1,l是固定值节点,其右子节点vi+1,r是混合节点,则 β i , v ( 2 k - 1 ) = α i , v ( 2 k ) + β i + 1 , r ( k ) β i , v ( 2 k ) = α i , v ( 2 k - 1 ) + β i + 1 , r ( k ) , k = 1 , ... , 2 n - i - 1 ; 若其左子节点vi+1,l是混合节点,其右子节点vi+1,r是信息节点,则若其左子节点vi+1,l和右子节点vi+1,r均是混合节点,则
然后,按照第n层到第0层层序遍历的顺序进行各节点对应的α值计算,对于遍历到的三类节点,按以下方法进行处理:当遍历到第i层任一固定值节点vi,f时,根据公式直接译出vi,f包含的各变量节点对应的固定值比特;当遍历到第i层任一信息节点vi,e时,根据公式计算vi,e获得的软信息αi,e值,直接对其进行判决计算得到g(αi,e),并对vi,e包含的所有变量节点对应的信息比特进行估计,其中,当遍历到第i层任一混合节点vi,m时,按照第0层到第n层层序遍历的顺序进行各混合节点对应的软信息α值计算;直至译出全部信源比特则完成译码。
与现有技术相比,本发明的有益效果在于:对极化码译码过程进行树图建模分析,将树图上的节点分类,证明了在译码过程中部分节点对应的运算是冗余的;通过对树图的剪枝,删除了这部分冗余的译码运算,从而提出了极化码SC译码和BP译码的简化算法,进而能够在确保不对误码性能产生影响的前提下,达到减小译码复杂度和提高译码效率的目的,最后进行了仿真验证,简化后SC译码和BP译码的译码复杂度较原始算法分别降低了36%~65%和41%~67%。
附图说明
图1是本发明一种实施例的极化码的树图示意图;
图2是本发明一种实施例极化码二分图中处理单元上的软信息传递示意图;
图3是本发明一种实施例码长为8的极化码二分图模型示意图;
图4是本发明一种实施例码长为8的SC简化译码过程示意图;
图5是本发明一种实施例码长为8的BP简化译码的译码前t-1次迭代译码过程示意图;
图6是本发明一种实施例极化码的误码性能与迭代次数关系示意图;
图7是本发明一种实施例的SC简化译码使计算复杂度减少的程度的示意图;
图8是本发明一种实施例的BP简化译码使计算复杂度减少的程度的示意图;
图9是本发明一种实施例SC简化译码中删除概率与误码性能的关系示意图;
图10是本发明一种实施例SC简化译码中误码性能的示意图;
图11是本发明一种实施例BP简化译码中删除概率与误码性能的关系示意图;
图12是本发明一种实施例BP简化译码中误码性能的示意图。
具体实施方式
下面结合附图,对本发明的较优的实施例作进一步的详细说明:
实施例1:
本例提供一种极化码的简化译码方法,包括以下步骤:
步骤S1,对接收到的极化码建立树图模型;
步骤S2,将树图中各个节点分为信息节点、固定值节点和混合节点;
步骤S3,对树图进行剪枝,剪去固定值节点和信息节点之间的子树,留下与混合节点相连接的固定值节点和信息节点;
步骤S4,对剪枝后的译码树图中按顺序从根节点开始进行递推译码运算。本例所述步骤S1中,所述极化码为(N,K,A,uA c)码;所述极化码的简化译码方法也称为基于树图剪枝的SC译码简化方法。
信道极化是将N个二元离散无记忆信道经过信道组合和分离,使部分信道的信道容量趋于1,另一部分信道的信道容量趋于0。极化码的编码原理是对应信道容量趋于1的部分编入信息比特,而对应其他部分编入固定值比特,从而确保信息传输的可靠性。极化码的编码过程可以写为GN,其中是未编码的N个信源比特值(u1,u2,…,uN),GN是极化码N阶极化码生成矩阵,是N个极化码比特值(x1,x2,…,xN)。信源比特由信息比特uA和固定值比特uA c两部分构成,为构造极化码时挑选的K维信息比特域集;信息比特uA是以集合A为指数的信源比特uA={ui,i∈A},固定值比特uA c是以集合A的补集为指数的信源比特。挑选得到K维信息域集合A={a1,…,aK}的原则是使各传输信息比特信道的Bhattacharyya参数 Z ( W N ( i ) ) = Δ Σ y i Σ u j W N ( i ) ( y 1 N , u 1 i - 1 / 0 ) W N ( i ) ( y 1 N , u 1 i - 1 / 1 ) 之和尽可能的小,其中是将N个独立信道组合并分离后,第i个分离信道的信道转移概率,是N个对应的信道观测(y1,y2,…,yN)。在二元删除信道(BEC)情况下,可以利用Bhattacharyya参数的递推公式
Z ( W N ( 2 i ) ) = Z ( W N / 2 ( i ) ) 2 Z ( W N ( 2 i - 1 ) ) = 2 Z ( W N / 2 ( i ) ) - Z ( W N / 2 ( i ) ) 2 计算出每个分离信道的Bhattacharyya参数并从小到大进行排序,选择Bhattacharyya参数小的前K位作为信息比特域A。这种信道选择方法也可以应用于加性高斯白噪声信道(AWGN),具有较好的效果。极化码的生成矩阵为其中BN为N阶翻转矩阵, F = 1 0 1 1 , 为矩阵F的n阶Kronecker积。由此生成的极化码表示为(N,K,A,uA c)码。
极化码的传统译码方法主要有SC译码和BP译码两种。SC译码通过接收得到信道观测的似然函数对信源比特ui,i∈1,…,N逐位进行估计,每次对ui进行估计运算时,将之前已译出的前i-1位信源比特作为已知加入运算。软判决译码过程需要得到信源比特ui对应的软信息对数似然比 L N ( i ) ( y 1 N , u ^ 1 i - 1 ) = l n W N ( i ) ( y 1 N , u ^ 1 i - 1 | 0 ) W N ( i ) ( y 1 N , u ^ 1 i - 1 | 1 ) , 根据转移函数的递推关系,,可以得到软信息对数似然比的递推计算公式
其中,上述公式中,是一种实数运算:对于两个实数a、b,有=2tanh-1(tanh(a/2)·tanh(b/2)),是模二和运算, 定义函数g(·)是一种判决运算如下,对于实数c有 g ( c ) = 0 , c &GreaterEqual; 0 g ( c ) = 1 , c < 0 ; 由于双曲正切函数是奇函数,所以有根据对各独立信道估计得到的软信息对数似然比初始值通过上述软信息对数似然比的递推计算公式依次计算得到信源比特ui对应的软信息对数似然比从而进一步得到信源比特ui的估计值
BP译码是一种被广泛应用的迭代译码算法,极化码BP译码的似然比信息传递方式可由其二分图得到。可以用n阶的二分图来表示码长为N=2n的(N,K,A,uA c)极化码,二分图中共有(n+1)N个节点和nN/2个处理单元。图2表示的是极化码二分图中的处理单元,由四个节点组成,二分图的每一阶中均有N/2个处理单元。如图3所示的是码长为8的极化码二分图模型,每一阶均有4个处理单元。是节点(i,j)在第t次迭代过程中向左传递的软信息对数似然比,是节点(i,j)在第t次迭代过程中向右传递的软信息对数似然比,其中i∈1,…,n+1,j∈1,…,N,初始值为 &beta; 1 , j 0 = 0 , j &Element; A &infin; , j &Element; A c .
每一次迭代,二分图中每个中间节点都向左和向右传递似然信息。每个处理单元中软信息α值和β值传递的过程为在设定的最大迭代次数tmax次迭代后,利用二分图中左侧第一列各节点得到的软信息计算对应的各信源比特uj估计值
本例所述步骤S1中,对接收到的码长为N=2n的极化码,将其二分图模型表示为一个n+1层的树图形式,该树图中共有2n+1-1个节点,其中处于第n层的各变量节点对应信源比特,而处于第0层的根节点对应编码节点;其中,n为大于1的自然数。所述步骤S2中,对于第i层的任一节点,如果该层包含的变量节点全部对应固定值比特,则称其为固定值节点;如果该层包含的变量节点全部对应信息比特,则称其为信息节点;如果该层包含的变量节点对应的既有固定值比特又有信息比特,则称其为混合节点。
如图1所示,在所述SC译码的建模及分析中,码长为N=2n的极化码,可以将其二分图模型表示为一个n+1层的树图形式,该树图中共有2n+1-1个节点,其中处于第n层的各变量节点对应信源比特,而处于第0层的根节点对应编码节点。图1是(8,3,{6,7,8},(0,0,0,0,0))极化码(以下简称(8,3)极化码)的4层树图表示示意,对于第i层的任一节点vi,i∈1,…,n,如果它包含的变量节点全部对应固定值比特,那么称其为固定值节点;如果它包含的变量节点全部对应信息比特,那么称其为信息节点;如果它包含的变量节点对应的既有固定值比特又有信息比特,那么称其为混合节点。于是,固定值节点的子树中均为固定值节点;信息节点的子树中均为信息节点;混合节点的子树包含信息节点和固定值节点,也有可能包含混合节点。
译码过程是根据信道观测得到信源比特估计值的过程,从图1的树图上看,极化码译码过程即软信息从第0层向第n层传递的过程。SC译码逐位对信源比特进行估计,于是软信息对数似然比L需按照“先序遍历”的顺序从根节点递推传递到变量节点。“先序遍历”的顺序是先访问根节点,然后遍历左子树,最后遍历右子树,在遍历左、右子树时,仍然先访问此时的“根节点”,然后遍历其左子树,最后遍历其右子树。而当遍历右子树时,总是需要左子树提供辅助信息
从而可以将信息传输的递推过程改写,处于第i层的节点vi,i∈1,…,n对应的译码运算为:节点vi通过接收到i-1层传递的软信息对数似然比Li,v,计算并向处于第i+1层的其左子节点vi+1,l和右子节点vi+1,r所分别传递的软信息Li+1,l和软信息Li+1,r节点vi通过接收到i+1层传递的辅助信息Ui+1,l和Ui+1,r,计算并向处于第i-1层的其父节点vi-1,p传递的辅助信息Ui,v U i , v ( 2 k - 1 ) = U i + 1 , l ( k ) &CirclePlus; U i + 1 , r ( k ) U i , v ( 2 k ) = U i + 1 , r ( k ) , k = 1 , ... , 2 n - i - 1 .
每当计算右子树的软信息对数似然比L的过程中需要左子树的辅助信息U,U的递推方向同L是相反的。处于第n层的各变量节点vn,j向上提供的辅助信息为其中,是变量节点vn,j对应的信源比特估计值。
值得一提的是,所述SC译码简化算法为,将SC译码树图中的节点进行分类,其中固定值节点和信息节点对应的译码运算是冗余的,不对译码结果产生任何影响,其证明如下。
第i层任一固定值节点vi,f包含的第n层各变量节点对应固定值比特un,f[j],j=1,...,2n-i,其中,i∈1,…,n,固定值比特在译码过程中认为是已知的于是固定值节点vi,f获得的软信息Li,f与固定值比特的估计值无关,固定值节点vi,f对应的译码运算是冗余的。
第i层任一信息节点vi,e包含的第n层各变量节点对应信息比特un,e[j],j=1,…,2n -i,,其中,i∈1,…,n,对信息比特的估计过程即是对各变量节点vn,e[j]获得的软信息Ln,e[j]的判决计算。因为vi,e子树中的译码过程不涉及已知的固定值比特,所以Ln,e[j]只与vi,e获得的软信息Li,e有关,即un,e[j]=g(Ln,e[j])=g(f1(Li,e)),j=1,…,2n-i;其中f1为根据递推公式得到的映射。此时只需用Li,e进行译码,而信息节点vi,e及其子树中各节点对应的译码运算均是冗余的,不对译码结果产生任何影响,证毕。
下面给出一种极化码树图中节点分类的方法:(N,K,A,uA c)极化码的信息域A={a1,…,aK}由初始信道参数和式(1)确定,对于树图第i层左数第m个的节点vim,m∈{1,…,2i},其中ceil是向上取整函数,count是重复数计数函数,即对集合中值为m的项进行计数,为构造极化码时挑选的K维信息比特域集合;若进而可以对树图每一层的节点进行分类。
值得一提的是,根据以上非常严谨的推论和计算,对译码树图上的固定值节点和信息节点进行剪枝处理,从而删去冗余的译码运算,达到简化译码的目的。
因此,本例所述步骤S2中,对于树图第i层左数第m个的节点vim,m∈{1,…,2i},通过进行计算,得到各个节点的分类,其中ceil是向上取整函数,count是重复数计数函数,即对集合中值为m的项进行计数,为构造极化码时挑选的K维信息比特域集合;当进而能够对树图每一层的节点进行分类。
本例所述步骤S4中,将接收到的信道观测软信息对数似然比L输入根节点,然后在剪枝后的SC译码树图中依次按照访问根节点、遍历左子树和遍历右子树的顺序从根节点开始进行递推译码运算。
本例所述步骤S4中,在遍历到第i层任一固定值节点vi,f时,i∈1,…,n,根据公式直接译出vi,f包含的各变量节点vn,f[j]所对应的固定值比特un,f[j],其中,j=1,...,2n-i;并根据计算vi,f对应的辅助信息Ui,f并送至vi,f的父节点,其中,vi,f为固定值节点,un,f[j]为第i层任一固定值节点vi,f包含的第n层各变量节点对应的固定值比特,为固定值比特的估计值,是极化码2n-i阶生成矩阵;在遍历到第i层任一信息节点vi,e时,i∈1,…,n,对信息节点vi,e获得的软信息Li,e进行判决计算得到信息节点vi,e对应的辅助信息Ui,e=g(Li,e),使用直接对信息节点vi,e包含的变量节点vn,e[j]对应的信息比特un,e[j]进行估计,并将辅助信息Ui,e送至信息节点vi,e的父节点,其中,un,e[j]为第i层任一信息节点vi,e包含的第n层各变量节点对应的信息比特,为信息比特的估计值;在遍历到第i层任一混合节点vi,m时,节点vi通过接收到i-1层传递的软信息对数似然比Li,v,计算并向处于第i+1层的其左子节点vi+1,l和右子节点vi+1,r所分别传递的软信息Li+1,l和Li+1,r节点vi通过接收到i+1层传递的辅助信息Ui+1,l和Ui+1,r,计算并向处于第i-1层的其父节点vi-1,p传递的辅助信息Ui,v U i , v ( 2 k - 1 ) = U i + 1 , l ( k ) &CirclePlus; U i + 1 , r ( k ) U i , v ( 2 k ) = U i + 1 , r ( k ) , k = 1 , ... , 2 n - i - 1 , 其中,为实数运算符号,为模二和运算符号;处于第n层的任一变量节点vn,j向上提供的辅助信息为其中,是变量节点vn,j对应的信源比特估计值。
下面举例说明(8,3)码的简化译码算法,如图4所示,首先建立(8,3)的译码树图,如图4所示,并根据信息域A={6,7,8}对节点进行分类。对树图进行剪枝,图4中虚线部分是被剪枝的路径;图4中标示了对剪枝后树图的遍历顺序,树图节点v11接收到的信道观测信息,遍历进行到固定值节点v21,根据计算v21对应的辅助信息U21,由于(8,3,{6,7,8},(0,0,0,0,0))码的固定值比特均为0,所以U21为0,将辅助信息U21送至v11;遍历v11的右子树,进行到混合节点v22,按原始算法进行运算并继续遍历,进行到节点v33;当运算进行到信息节点v46时,根据式译出所含变量节点对应的信息比特并将其作为辅助信息送至v33;当运算进行到节点v34,根据式 ( u ^ n , e [ 1 ] , . . . , u ^ n , e [ 2 n - i ] ) = U i , e G 2 n - i 译出信息比特后,译码运算完成。
复杂度分析如下:码长为N的极化码原始SC译码算法的计算复杂度为O(NlogN)。译码过程包含了N/2·logN次软信息L值操作和N/2·logN次辅助信息U值操作。每一次L值操作中有一次加法运算和一次运算,每一次U值操作中有一次模二和运算。例如(8,3)极化码,原始SC译码算法的运算过程包含了12次加法运算、12次运算和12次模二和运算,而(8,3)极化码简化SC译码运算过程中,只需要进行7次加法运算、2次运算、2次模2和运算及一次“×G2”的运算。“×GN”运算可以等效成极化码编码过程,即N/2·logN次模二和运算。
实施例2:
在实施例1的基础上,与实施例1不同的是,本例主要是通过基于树图剪枝的BP译码简化方法,所述基于树图剪枝的BP译码简化方法也就是所述极化码的简化译码方法其中,BP译码的建模及分析如下。
利用图1中树图模型对极化码BP译码算法进行分析。处于第i层的节点vi,i∈1,…,n对应的译码计算为:节点vi通过接收到i-1层传递的软信息对数似然比αi,v,计算并向处于第i+1层的其左子节点vi+1,l和右子节点vi+1,r传递的软信息αi+1,l和软信息αi+1,r节点vi通过接收到i+1层传递的软信息βi+1,l、βi+1,r,计算并向处于第i-1层的其父节点vi-1,p传递的软信息βi,v
每次迭代中软信息β值的计算需要前一次迭代中得到的软信息α值代入。β值的递推计算顺序是从信源比特逐步向信道观测进行的,在树图上,第n层的变量节点对应信源比特而第0层根节点对应信道观测,于是β值的递推计算过程可以认为是对树图从第n层到第0层层序遍历的过程。从第n层到第0层层序遍历,即从第n层向上逐层从左到右遍历节点,直到第0层。β值的递推计算结束之后,对α值进行更新,α值的递推计算顺序与β值相反,于是传递过程可以认为是对树图从第0层到第n层层序遍历的过程。
运行到迭代次数上限后,对第n层每个变量节点最终得到的软信息α值进行判决,即得到对应信源比特的估计值,完成译码,BP译码简化译码方法为:设定极化码BP译码算法的迭代次数为t,BP译码树图中同样也将节点分为固定值节点、信息节点和混合节点三类,于是根据式 &beta; 1 , j 0 = 0 , j &Element; A &infin; , j &Element; A c 对于第i层不同可能类型节点vi对应的β值计算可写为:当vi是固定值节点,则βi,v(2k-1)=βi,v(2k)=∞,k=1,…,2n-i-1;当vi是信息节点,则βi,v(2k-1)=βi,v(2k)=0,k=1,…,2n-i-1;当vi是混合节点,其左子节点vi+1,l和右子节点vi+1,r可能类型如下:若左子节点vi+1,l是固定值节点,右子节点vi+1,r是信息节点,则 &beta; i , v ( 2 k - 1 ) = &alpha; i , v ( 2 k ) &beta; i , v ( 2 k ) = &alpha; i , v ( 2 k - 1 ) , k = 1 , ... , 2 n - i - 1 ; 若左子节点vi+1,l是固定值节点,右子节点vi+1,r是混合节点,则 &beta; i , v ( 2 k - 1 ) = &alpha; i , v ( 2 k ) + &beta; i + 1 , r ( k ) &beta; i , v ( 2 k ) = &alpha; i , v ( 2 k - 1 ) + &beta; i + 1 , r ( k ) , k = 1 , ... , 2 n - i - 1 ; 若左子节点vi+1,l是混合节点,右子节点vi+1,r是信息节点,则若左子节点vi+1,l和右子节点vi+1,r均是混合节点,则其中,节点vi通过接收到i+1层传递的软信息βi+1,l和βi+1,r,计算并向处于第i-1层的其父节点vi-1,p传递其软信息βi,v
由于在信道选择过程中,每个节点的左子代节点的Bhattacharyya参数大于右子代节点的Bhattacharyya参数,使得左子代节点的信道选择优先级弱于右子代节点,所以混合节点的结构只会存在以上四种情形。此时,只有混合节点对应的β值计算,需要用到上一次迭代得到的α值,而固定值节点和信息节点的β值分别为固定的∞和0。
设定极化码BP译码方法的迭代次数为t,将BP译码树图中的节点进行分类,BP译码树图中固定值节点和信息节点对应的译码计算是冗余的,不对译码结果产生任何影响。其证明如下:对于β值计算,根据式βi,v(2k-1)=βi,v(2k)=∞,k=1,…,2n-i-1和βi,v(2k-1)=βi,v(2k)=0,k=1,…,2n-i-1,固定值节点和信息节点传递的β值是固定的∞和0,所以这两类节点对应的β值计算是冗余的。
由于每次迭代中固定值节点和信息节点传递的β值均是固定的,所以并不需要在前一次迭代中进行固定值节点和信息节点对应的α值计算,于是在前t-1次迭代中,固定值节点和信息节点对应的α值计算不对译码结果产生任何影响。第t次迭代中,固定值比特在译码过程中认为是已知,于是固定值节点对应的α值计算与译码结果无关。信息比特的译码是对第i层信息节点vi,e包含的第n层各变量节点vn,e[j],j=1,...,2n-i得到软信息αn,e[j]的判决计算,由于信息节点的β值为0,所以节点vi,e对应的软信息α值计算可写为:依次递推到第n层,此时软信息αn,e[j]只与αi,e有关, u ^ n , e &lsqb; j &rsqb; = g ( &alpha; n , e &lsqb; j &rsqb; ) = f 2 ( &alpha; i , e ) , j = 1 , ... , 2 n - i ; 其中f2是根据得到的映射。此时只需使用αi,e进行译码,而信息节点vi,e及其子树中各节点对应的α值计算均是冗余的,不对译码结果产生任何影响,证毕。
对于混合节点vi,m,仍需要进行对应的α值计算,对vi,m的左子节点vi+1,l和右子节点vi+1,r分类,若vi+1,l是固定值节点,vi+1,r是混合节点,则αi+1,r(k)=αi,v(2k-1)+αi,v(2k),k=1,…,2n-i-1;若vi+1,l是混合节点,vi+1,r是信息节点,则若vi+1,l和vi+1,r均是混合节点,则然后将计算结果依次代入下一次迭代过程中,直至前t-1次迭代过程完成。
根据上述证明,对极化码BP译码前t-1次迭代和第t次迭代中的树图分别进行剪枝处理,从而删去冗余的译码运算,达到简化译码的目的。由此提出基于树图剪枝的BP简化译码方法,即提出极化码的简化译码方法以下步骤:步骤S1,对接收到的(N,K,A,uA c)极化码建立树图模型;步骤S201,将树图中各节点分为信息节点、固定值节点和混合节点;步骤S202,对树图进行剪枝,剪去全部固定值节点和信息节点,剪枝后译码树图中仅存在混合节点,如图5所示,就是(8,3)码BP译码前t-1次迭代过程译码树图的剪枝示意图;步骤S203,译码运算开始前,将初始值按公式输入根节点;其中,是节点(n+1,j)在第0次迭代过程中向左传递的软信息对数似然比,W1 (1)是原始独立信道的信道转移概率;步骤S204,译码运算开始,对剪枝后树图上的各混合节点对应的β值计算,按照第n层到第0层层序遍历的顺序进行。完成各混合节点对应的软信息β值计算后,按照第0层到第n层层序遍历的顺序进行各混合节点对应的软信息α值计算,将计算结果代入下一次迭代过程中。依次进行,直至前t-1次迭代过程完成;步骤S205,保存第t-1次迭代后的各节点对应的软信息α值结果,重新建立(N,K,A,uA c)极化码树图模型;步骤S3,对树图进行剪枝,剪去固定值节点和信息节点的子树,即只留下与混合节点相连接的固定值节点和信息节点,其剪枝方法与实施例1中的图4中SC译码树图剪枝相同;步骤S4,进行译码运算的第t次迭代过程。计算剪枝后树图上的各混合节点对应的软信息β值,根据已保存的第t-1次迭代后的α值,按照第n层到第0层层序遍历的顺序进行。完成各混合节点的β值计算后,进行剪枝后树图上的各节点对应的α值计算,按照第n层到第0层层序遍历的顺序进行,对于遍历到的三类节点。这里所述步骤S201至步骤S205都是本例所述步骤S2的子步骤,所述步骤S4对剪枝后的译码树图进行译码运算的第t次迭代过程。而实施例1为SC的简化译码方法,实施例2为BP的简化译码方法,即SC和BP剪枝简化译码的区别在于,BP的简化译码方法中,包含了两次剪枝过程,即所述步骤S2中以及所述步骤S3中的剪枝过程。三类节点按以下方法进行处理:若遍历到第i层任一固定值节点vi,f时,i∈1,…,n,直接译出vi,f包含的各变量节点对应的固定值比特;若遍历到第i层任一信息节点vi,e时,根据式计算vi,e获得的软信息αi,e值,直接对其进行判决计算得到g(αi,e),并对vi,e包含的所有变量节点对应的信息比特进行估计 ( u ^ n , e &lsqb; 1 &rsqb; , ... , u ^ n , e &lsqb; 2 n - i &rsqb; ) = g ( &alpha; i , e ) G 2 n - i ; 若遍历到第i层任一混合节点vi,m时,进行其对应的软信息α值计算,直至译出全部信源比特,译码完成。
复杂度分析如下:码长为N的极化码原始BP译码算法的计算复杂度为O(tNlogN),t为迭代次数。译码器中共有N/2·logN个处理单元,每个处理单元在每次迭代过程中都要进行一次α值操作和一次β值操作。每一次α值、β值操作均含有两次运算和两次加法运算。仿真迭代次数t对BP译码性能的影响,如图6所示,AWGN信道下(1024,256)极化码简化BP译码算法的误码性能仿真。可以看出,BP译码迭代次数取20,即可获得较好的译码性能。简化BP译码方法,在前t-1次迭代过程中,软信息只在混合节点之间传递。软信息α值从处于第i层的节点vi传递到第i+1层的混合节点vi+1,m至多需要2n-i-1运算和加法运算,软信息β值从vi向上层传递的运算量亦是如此。例如(8,3)码简化BP译码算法,在前t-1次迭代过程中,每次迭代只需进行4次运算和8次加法运算。在第t次迭代中,α值的运算数量较之前各次迭代增加至多c次运算和加法运算,同时增加了c次“×Gi”运算,c为信息节点树的数量,信息节点树是由信息节点构成的子树。(8,3)码简化BP译码算法,第t次迭代中需进行6次运算、11次加法运算和1次“×G2”运算,而原始BP译码算法每次迭代需进行48次运算和加法运算。当t较大时,可以用前t-1次迭代过程的平均运算量近似得到总体译码计算复杂度,忽略第t次迭代中增加的计算量对总体复杂度的影响。
实施例1和实施例2中,对极化码译码过程进行树图建模分析,将树图上的节点分类,证明了在译码过程中部分节点对应的运算是冗余的;通过对树图的剪枝,删除了这部分冗余的译码运算,从而提出了极化码SC译码和BP译码的简化算法,进而能够在确保不对误码性能产生影响的前提下,达到减小译码复杂度和提高译码效率的目的,最后进行了仿真验证,简化后SC译码和BP译码的译码复杂度较原始算法分别降低了36%~65%和41%~67%。
在实施例1提出的SC译码和BP译码的极化码的简化译码方法的基础上,这里对不同码率下简化前后的SC译码和BP译码,进行计算复杂度仿真比较,得到简化算法使译码复杂度减少的程度。在SC译码和BP译码过程中,由于运算占用了大部分计算量,所以可以用运算次数来近似统计译码复杂度。
取初始信道参数Z(W)=0.6,对码长为210、213、216的极化码,仿真不同码率下对简化译码算法计算较原始译码减少的百分比,如图7和图8所示。仿真结果显示,简化过程使不同码率极化码SC译码和BP译码的计算复杂度均有较大程度的减少。其中,图7SC译码的极化码的简化译码方法的计算复杂度较原始算法减少了36%~65%。图8BP译码的极化码的简化译码方法的计算复杂度较原始算法减少了41%~67%。
误码性能仿真与分析如下:对极化码简化前后的SC译码和BP译码分别在不同信道环境下进行误码性能仿真对比,并与LDPC码[10]的误码性能进行比较,仿真中BP译码算法的迭代次数均为20。
由图9和图12的仿真结果显示,所述极化码的简化译码方法简化前后的SC译码和BP译码在误码性能上没有任何区别,验证了误码性能不改变的前提。图9的仿真结果显示,在BEC下,码率取0.25时,极化码在较大删除概率时的误码性能优于LDPC码,而且此时极化码的编译码复杂度远低于LDPC码。极化码BP译码的误码性能较SC译码在删除概率上提高了约0.02。图10的仿真结果显示,在AWGN下,码率取0.25,极化码的误码性能与LDPC码相差无几,极化码BP译码的误码性能较SC译码提高了约0.2dB。
图11和图12的仿真结果显示,码率取0.5,LDPC码的误码性能优于极化码。在BEC下,极化码BP译码的误码性能较简化SC译码在删除概率上提高了约0.02。在AWGN下,极化码BP译码的误码性能较SC译码提高了约0.2dB。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (10)

1.一种极化码的简化译码方法,其特征在于,包括以下步骤:
步骤S1,对接收到的极化码建立树图模型;
步骤S2,将树图中各个节点分为信息节点、固定值节点和混合节点;
步骤S3,对树图进行剪枝,剪去固定值节点和信息节点之间的子树,留下与混合节点相连接的固定值节点和信息节点;
步骤S4,对剪枝后的译码树图进行译码运算。
2.根据权利要求1所述的极化码的简化译码方法,其特征在于,所述步骤S1中,对接收到的码长为N=2n的极化码,将其二分图模型表示为一个n+1层的树图形式,该树图中共有2n +1-1个节点,其中处于第n层的各变量节点对应信源比特,而处于第0层的根节点对应编码节点;其中,n为大于1的自然数。
3.根据权利要求2所述的极化码的简化译码方法,其特征在于,所述步骤S2中,对于第i层的任一节点,如果该层包含的变量节点全部对应固定值比特,则称其为固定值节点;如果该层包含的变量节点全部对应信息比特,则称其为信息节点;如果该层包含的变量节点对应的既有固定值比特又有信息比特,则称其为混合节点。
4.根据权利要求3所述的极化码的简化译码方法,其特征在于,所述步骤S2中,对于树图第i层左数第m个的节点vim,m∈{1,…,2i},通过进行计算,得到各个节点的分类,其中ceil是向上取整函数,count是重复数计数函数,对集合中值为m的项进行计数,为构造极化码时挑选的K维信息比特域集;进而能够对树图每一层的节点进行分类,i为大于1的自然数。
5.根据权利要求1至4任意一项所述的极化码的简化译码方法,其特征在于,所述步骤S4中,将接收到的信道观测软信息对数似然比L输入根节点,然后在剪枝后的SC译码树图中依次按照访问根节点、遍历左子树和遍历右子树的顺序从根节点开始进行递推译码运算。
6.根据权利要求5所述的极化码的简化译码方法,其特征在于,所述步骤S4中,在遍历到第i层任一固定值节点vi,f时,i∈1,…,n,根据公式直接译出vi,f包含的各变量节点vn,f[j]所对应的固定值比特un,f[j],其中,j=1,…,2n-i;并根据计算vi,f对应的辅助信息Ui,f并送至vi,f的父节点,其中,vi,f为固定值节点,un,f[j]为第i层任一固定值节点vi,f包含的第n层各变量节点对应的固定值比特,为固定值比特的估计值,是极化码2n-i阶生成矩阵;在遍历到第i层任一信息节点vi,e时,i∈1,…,n,对信息节点vi,e获得的软信息Li,e进行判决计算得到信息节点vi,e对应的辅助信息Ui,e=g(Li,e),使用直接对信息节点vi,e包含的变量节点vn,e[j]对应的信息比特un,e[j]进行估计,并将辅助信息Ui,e送至信息节点vi,e的父节点,其中,un,e[j]为第i层任一信息节点vi,e包含的第n层各变量节点对应的信息比特,为信息比特的估计值;在遍历到第i层任一混合节点vi,m时,节点vi通过接收到i-1层传递的软信息对数似然比Li,v,计算并向处于第i+1层的其左子节点vi+1,l和右子节点vi+1,r所分别传递的软信息Li+1,l和软信息Li+1,r节点vi通过接收到i+1层传递的辅助信息Ui+1,l和Ui+1,r,计算并向处于第i-1层的其父节点vi-1,p传递的辅助信息Ui,v U i , v ( 2 k - 1 ) = U i + 1 , l ( k ) &CirclePlus; U i + 1 , r ( k ) U i , v ( 2 k ) = U i + 1 , r ( k ) , k = 1 , ... , 2 n - i - 1 , 其中,为实数运算符号,为模二和运算符号;处于第n层的任一变量节点vn,j向上提供的辅助信息为其中,是变量节点vn,j对应的信源比特估计值。
7.根据权利要求1所述的极化码的简化译码方法,其特征在于,步骤S2中,将树图中各个节点分为信息节点、固定值节点和混合节点;对树图进行首次剪枝,在译码运算开始前,将初始值按公式输入根节点;其中,是节点(n+1,j)在第0次迭代过程中向左传递的软信息对数似然比,是原始独立信道的信道转移概率。
8.根据权利要求7所述的极化码的简化译码方法,其特征在于,所述步骤S2中,对剪枝后树图上的各混合节点对应的软信息β值计算,按照第n层到第0层层序遍历的顺序进行译码,在遍历到第i层任一混合节点vi,m时,若其左子节点vi+1,l是固定值节点,其右子节点vi+1,r是信息节点,则 &beta; i , v ( 2 k - 1 ) = &alpha; i , v ( 2 k ) &beta; i , v ( 2 k ) = &alpha; i , v ( 2 k - 1 ) , k = 1 , ... , 2 n - i - 1 ; 若其左子节点vi+1,l是固定值节点,其右子节点vi+1,r是混合节点,则 &beta; i , v ( 2 k - 1 ) = &alpha; i , v ( 2 k ) + &beta; i + 1 , r ( k ) &beta; i , v ( 2 k ) = &alpha; i , v ( 2 k - 1 ) + &beta; i + 1 , r ( k ) , k = 1 , . . . , 2 n - i - 1 ; 若其左子节点vi+1,l是混合节点,其右子节点vi+1,r是信息节点,则若其左子节点vi+1,l和右子节点vi+1,r均是混合节点,则其中,为实数运算符号,节点vi通过接收到i+1层传递的软信息βi+1,l和软信息βi+1,r,计算并向处于第i-1层的其父节点vi-1,p传递其软信息βi,v;节点vi通过接收到i-1层传递的软信息对数似然比αi,v,计算并向处于第i+1层的其左子节点vi+1,l和其右子节点vi+1,r所分别传递的软信息αi+1,l和软信息αi+1,r
9.根据权利要求8所述的极化码的简化译码方法,其特征在于,所述步骤S2中,完成各混合节点对应的软信息βi,v计算后,按照第0层到第n层层序遍历的顺序进行各混合节点对应的软信息α值计算,若vi+1,l是固定值节点,vi+1,r是混合节点,则αi+1,r(k)=αi,v(2k-1)+αi,v(2k),k=1,…,2n-i-1;若vi+1,l是混合节点,vi+1,r是信息节点,则若vi+1,l和vi+1,r均是混合节点,则然后将计算结果依次代入下一次迭代过程中,直至前t-1次迭代过程完成,保存第t-1次迭代后的各节点对应的软信息α值结果,然后重新建立极化码的树图模型。
10.根据权利要求9所述的极化码的简化译码方法,其特征在于,所述步骤S4中,进行译码运算的第t次迭代过程,计算剪枝后树图上的各混合节点对应的软信息β值,根据已保存的第t-1次迭代后的α值;
首先,在遍历到第i层任一混合节点vi,m时,若其左子节点vi+1,l是固定值节点,其右子节点vi+1,r是信息节点,则 &beta; i , v ( 2 k - 1 ) = &alpha; i , v ( 2 k ) &beta; i , v ( 2 k ) = &alpha; i , v ( 2 k - 1 ) , k = 1 , ... , 2 n - i - 1 ; 若其左子节点vi+1,l是固定值节点,其右子节点vi+1,r是混合节点,则 &beta; i , v ( 2 k - 1 ) = &alpha; i , v ( 2 k ) + &beta; i + 1 , r ( k ) &beta; i , v ( 2 k ) = &alpha; i , v ( 2 k - 1 ) + &beta; i + 1 , r ( k ) , k = 1 , ... , 2 n - i - 1 ; 若其左子节点vi+1,l是混合节点,其右子节点vi+1,r是信息节点,则若其左子节点vi+1,l和右子节点vi+1,r均是混合节点,则
然后,按照第n层到第0层层序遍历的顺序进行各节点对应的α值计算,对于遍历到的三类节点,按以下方法进行处理:当遍历到第i层任一固定值节点vi,f时,根据公式直接译出vi,f包含的各变量节点对应的固定值比特;当遍历到第i层任一信息节点vi,e时,根据公式计算vi,e获得的软信息αi,e值,直接对其进行判决计算得到g(αi,e),并对vi,e包含的所有变量节点对应的信息比特进行估计,其中,当遍历到第i层任一混合节点vi,m时,按照第0层到第n层层序遍历的顺序进行各混合节点对应的软信息α值计算;直至译出全部信源比特则完成译码。
CN201610045755.2A 2016-01-22 2016-01-22 一种极化码的简化译码方法 Active CN105720992B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610045755.2A CN105720992B (zh) 2016-01-22 2016-01-22 一种极化码的简化译码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610045755.2A CN105720992B (zh) 2016-01-22 2016-01-22 一种极化码的简化译码方法

Publications (2)

Publication Number Publication Date
CN105720992A true CN105720992A (zh) 2016-06-29
CN105720992B CN105720992B (zh) 2019-07-23

Family

ID=56153986

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610045755.2A Active CN105720992B (zh) 2016-01-22 2016-01-22 一种极化码的简化译码方法

Country Status (1)

Country Link
CN (1) CN105720992B (zh)

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106487479A (zh) * 2016-09-27 2017-03-08 清华大学深圳研究生院 一种基于多位判决的极化码译码方法
CN106656213A (zh) * 2016-12-22 2017-05-10 东南大学 基于k段分解的低复杂度极化码折叠硬件构架的实现方法
CN106656214A (zh) * 2016-12-22 2017-05-10 东南大学 一种基于串行抵消列表极化码译码的动态分布排序算法
CN107888202A (zh) * 2017-11-29 2018-04-06 哈尔滨工业大学 一种非递归的sc译码似然比确定方法及装置
CN108173624A (zh) * 2018-01-22 2018-06-15 合肥工业大学 一种部分译码的极化码串行抵消译码电路及其方法
CN108304658A (zh) * 2018-02-02 2018-07-20 中国计量大学 一种基于fpga的极化码编码器硬件实现方法
CN108390677A (zh) * 2018-03-20 2018-08-10 山东大学 一种极化码优化的编译码方法
CN108566213A (zh) * 2018-04-25 2018-09-21 东南大学 极化码的串行抵消列表比特翻转译码方法
CN108847848A (zh) * 2018-06-13 2018-11-20 电子科技大学 一种基于信息后处理的极化码的bp译码算法
CN108964671A (zh) * 2017-05-27 2018-12-07 华为技术有限公司 一种译码方法及译码器
CN109194338A (zh) * 2018-11-07 2019-01-11 桂林电子科技大学 一种混合节点多比特处理的极化码译码方法
CN109347600A (zh) * 2018-09-26 2019-02-15 河南科技大学 一种极化码译码方法
WO2019041067A1 (zh) * 2017-08-26 2019-03-07 华为技术有限公司 极化译码的方法与装置
CN109547156A (zh) * 2018-10-18 2019-03-29 北京理工大学 polar短码的非相干检测方法、装置与电子设备
CN109586730A (zh) * 2018-12-06 2019-04-05 电子科技大学 一种基于智能后处理的极化码bp译码算法
CN109725307A (zh) * 2018-12-13 2019-05-07 四川九洲空管科技有限责任公司 一种框架剪切的二次应答数据链译码方法
CN109802690A (zh) * 2017-11-17 2019-05-24 中兴通讯股份有限公司 译码方法、装置和计算机可读存储介质
CN110324111A (zh) * 2018-03-31 2019-10-11 华为技术有限公司 一种译码方法及设备
CN110612738A (zh) * 2017-03-03 2019-12-24 索克普拉科学与工程公司 广义极化码
CN110999095A (zh) * 2017-07-10 2020-04-10 交流通讯有限公司 用于极化码的按块并行冻结位生成
CN111224677A (zh) * 2018-11-27 2020-06-02 华为技术有限公司 编码方法、译码方法及装置
CN111614437A (zh) * 2020-04-30 2020-09-01 北京航空航天大学 一种基于scl译码算法的极化码构造优化方法
CN112564716A (zh) * 2020-08-05 2021-03-26 新疆大学 一种基于剪枝迭代的pc-scma系统联合译码方法
TWI754670B (zh) * 2016-09-20 2022-02-11 南韓商三星電子股份有限公司 極化碼之並列連續消除解碼與連續消除清單解碼的裝置和方法、製造方法及建立方法
CN117375635A (zh) * 2023-11-09 2024-01-09 中国人民解放军军事科学院系统工程研究院 一种卫星通信极化码bp译码的几何表示方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102122966A (zh) * 2011-04-15 2011-07-13 北京邮电大学 基于信道极化的交错结构重复码的编码器及其编译码方法
CN105187073A (zh) * 2015-10-13 2015-12-23 东南大学 一种极化码的bp译码方法及装置
CN105262494A (zh) * 2015-10-13 2016-01-20 东南大学 一种具有迭代早停止机制的极化码bp译码方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102122966A (zh) * 2011-04-15 2011-07-13 北京邮电大学 基于信道极化的交错结构重复码的编码器及其编译码方法
CN105187073A (zh) * 2015-10-13 2015-12-23 东南大学 一种极化码的bp译码方法及装置
CN105262494A (zh) * 2015-10-13 2016-01-20 东南大学 一种具有迭代早停止机制的极化码bp译码方法

Cited By (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI754670B (zh) * 2016-09-20 2022-02-11 南韓商三星電子股份有限公司 極化碼之並列連續消除解碼與連續消除清單解碼的裝置和方法、製造方法及建立方法
CN106487479A (zh) * 2016-09-27 2017-03-08 清华大学深圳研究生院 一种基于多位判决的极化码译码方法
CN106487479B (zh) * 2016-09-27 2019-04-16 清华大学深圳研究生院 一种基于多位判决的极化码译码方法
CN106656213A (zh) * 2016-12-22 2017-05-10 东南大学 基于k段分解的低复杂度极化码折叠硬件构架的实现方法
CN106656214A (zh) * 2016-12-22 2017-05-10 东南大学 一种基于串行抵消列表极化码译码的动态分布排序算法
CN110612738A (zh) * 2017-03-03 2019-12-24 索克普拉科学与工程公司 广义极化码
CN108964671B (zh) * 2017-05-27 2020-12-15 华为技术有限公司 一种译码方法及译码器
CN108964671A (zh) * 2017-05-27 2018-12-07 华为技术有限公司 一种译码方法及译码器
CN110999095A (zh) * 2017-07-10 2020-04-10 交流通讯有限公司 用于极化码的按块并行冻结位生成
CN110999095B (zh) * 2017-07-10 2023-08-04 交流通讯有限公司 用于极化码的按块并行冻结位生成
WO2019041067A1 (zh) * 2017-08-26 2019-03-07 华为技术有限公司 极化译码的方法与装置
US11012092B2 (en) 2017-08-26 2021-05-18 Huawei Technologies Co., Ltd. Polar decoding method and apparatus
CN109802690B (zh) * 2017-11-17 2022-04-22 中兴通讯股份有限公司 译码方法、装置和计算机可读存储介质
CN109802690A (zh) * 2017-11-17 2019-05-24 中兴通讯股份有限公司 译码方法、装置和计算机可读存储介质
CN107888202A (zh) * 2017-11-29 2018-04-06 哈尔滨工业大学 一种非递归的sc译码似然比确定方法及装置
CN108173624B (zh) * 2018-01-22 2020-08-07 合肥工业大学 一种部分译码的极化码串行抵消译码电路及其方法
CN108173624A (zh) * 2018-01-22 2018-06-15 合肥工业大学 一种部分译码的极化码串行抵消译码电路及其方法
CN108304658A (zh) * 2018-02-02 2018-07-20 中国计量大学 一种基于fpga的极化码编码器硬件实现方法
CN108390677A (zh) * 2018-03-20 2018-08-10 山东大学 一种极化码优化的编译码方法
CN110324111A (zh) * 2018-03-31 2019-10-11 华为技术有限公司 一种译码方法及设备
CN110324111B (zh) * 2018-03-31 2020-10-09 华为技术有限公司 一种译码方法及设备
CN108566213B (zh) * 2018-04-25 2022-05-27 东南大学 极化码的串行抵消列表比特翻转译码方法
CN108566213A (zh) * 2018-04-25 2018-09-21 东南大学 极化码的串行抵消列表比特翻转译码方法
CN108847848B (zh) * 2018-06-13 2021-10-01 电子科技大学 一种基于信息后处理的极化码的bp译码算法
CN108847848A (zh) * 2018-06-13 2018-11-20 电子科技大学 一种基于信息后处理的极化码的bp译码算法
CN109347600B (zh) * 2018-09-26 2021-01-05 河南科技大学 一种极化码译码方法
CN109347600A (zh) * 2018-09-26 2019-02-15 河南科技大学 一种极化码译码方法
CN109547156A (zh) * 2018-10-18 2019-03-29 北京理工大学 polar短码的非相干检测方法、装置与电子设备
CN109194338A (zh) * 2018-11-07 2019-01-11 桂林电子科技大学 一种混合节点多比特处理的极化码译码方法
CN111224677A (zh) * 2018-11-27 2020-06-02 华为技术有限公司 编码方法、译码方法及装置
CN111224677B (zh) * 2018-11-27 2021-10-15 华为技术有限公司 编码方法、译码方法及装置
CN109586730A (zh) * 2018-12-06 2019-04-05 电子科技大学 一种基于智能后处理的极化码bp译码算法
CN109586730B (zh) * 2018-12-06 2020-07-07 电子科技大学 一种基于智能后处理的极化码bp译码算法
CN109725307A (zh) * 2018-12-13 2019-05-07 四川九洲空管科技有限责任公司 一种框架剪切的二次应答数据链译码方法
CN111614437A (zh) * 2020-04-30 2020-09-01 北京航空航天大学 一种基于scl译码算法的极化码构造优化方法
CN111614437B (zh) * 2020-04-30 2021-04-30 北京航空航天大学 一种基于scl译码算法的极化码构造优化方法
CN112564716A (zh) * 2020-08-05 2021-03-26 新疆大学 一种基于剪枝迭代的pc-scma系统联合译码方法
CN117375635A (zh) * 2023-11-09 2024-01-09 中国人民解放军军事科学院系统工程研究院 一种卫星通信极化码bp译码的几何表示方法及装置
CN117375635B (zh) * 2023-11-09 2024-05-03 中国人民解放军军事科学院系统工程研究院 一种卫星通信极化码bp译码的几何表示方法及装置

Also Published As

Publication number Publication date
CN105720992B (zh) 2019-07-23

Similar Documents

Publication Publication Date Title
CN105720992A (zh) 一种极化码的简化译码方法
CN109361402B (zh) 编码方法及编码设备
CN105811998A (zh) 一种基于密度演进的极化码构造方法及极化码编译码系统
CN101803206B (zh) 近码字的rom列表解码
CN101656541B (zh) Rs码的译码方法和装置
CN100425000C (zh) 双涡轮结构低密度奇偶校验码解码器及解码方法
CN106936446A (zh) 一种基于多进制ldpc码的高速译码器及其译码方法
US8103945B2 (en) Decoding method and decoding apparatus as well as program
CN104242957B (zh) 译码处理方法及译码器
Esfahanizadeh et al. A matrix completion approach to linear index coding problem
CN101373978B (zh) 一种Turbo码译码方法以及装置
CN110138390A (zh) 一种基于深度学习的极化码sscl算法译码器
CN107453807A (zh) 一种大气光通信信道模型的极化方法、装置及电子设备
Hu et al. A comparative study of polar code decoding algorithms
Yang et al. Nonlinear programming approaches to decoding low-density parity-check codes
CN114157309A (zh) 极化码译码方法、装置及系统
CN110661533A (zh) 优化译码器存储极化码译码性能的方法
Teng et al. Convolutional neural network-aided tree-based bit-flipping framework for polar decoder using imitation learning
CN104052499A (zh) 一种ldpc码的纠删译码方法及系统
CN104158549A (zh) 一种极性码译码方法及译码装置
CN102832950A (zh) 一种低密度奇偶校验码的误帧率估计方法
CN110190857B (zh) 一种crc辅助校验极化码译码方法和智能终端
CN113131950B (zh) 一种极化码的自适应连续消除优先译码方法
CN109525252A (zh) 基于简化三阶关键集合的极化码串行抵消列表译码方法
CN110324111A (zh) 一种译码方法及设备

Legal Events

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