删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

低复杂度极化码SCL译码算法

本站小编 哈尔滨工业大学/2019-10-24

低复杂度极化码SCL译码算法

刘士平,马林华,胡星,黄天宇

(空军工程大学 航空航天工程学院,西安710038)



摘要:

极化码的串行抵消列表(SCL)译码的逐次逐比特进行判决过程与路径度量值的计算筛选过程是整个译码系统复杂度与延迟的主要来源.在分析现有SCL及多比特判决译码的优缺点基础上,针对SCL译码造成的译码系统复杂度高和延时大的问题,将每组多比特码字(多个独立信道)视作一个整体,并在译码时通过信道合成构造为一个虚拟多比特信道,从而可以对多比特码字进行同步传输并根据信道递归公式进行同步判决译码.由此基于SCL译码的码树构造提出一种构造多比特虚拟信道的SCL译码算法,并结合设置译码码树节点阈值减少码树节点分裂次数的方法进一步增强了上述算法.在AWGN信道下的分别对虚拟2、4和8比特信道SCL译码的误码率及误帧率性能进行仿真.仿真结果表明在虚拟8比特信道情况下,预设阈值S=30时的译码性能接近传统SCL算法,且总历经节点数降低了63.7%,总加法次数是8比特同时判决译码算法的17%.此算法降低了译码算法的计算复杂度及硬件存储复杂度,更适合于硬件实现, 具有一定的实用价值.

关键词:  极化码  串行抵消列表译码  译码复杂度  阈值

DOI:10.11918/j.issn.0367-6234.201706116

分类号:TN911.22

文献标识码:A

基金项目:国家自然科学基金资助项目(61472442); 综合业务网国家重点实验室(西安电子科技大学)开放研究课题资助项目(INS-15-13); 航空科学基金资助项目(20155896025)



Alow complexity SCL decoding algorithm for polar codes

LIU Shiping,MA Linhua,HU Xing,HUANG Tianyu

(Institute of Aeronautics and Astronautics Engineering, Air Force Engineering University, Xi’an 710038, China)

Abstract:

The successive cancellation list (SCL) decoding algorithm of polar code is of bitwise, which is the main source of the decoding complexity and latency. In view of the analysis of the advantages and disadvantages of the existing SCL and multi-bit decision decoding, each set of multi-bit code words (multiple independent channels) is treated as a whole to construct a multi-bit virtual channel, which can transmit multi-bit code words and decode synchronously. Therefore a construction of multi-bit virtual channel based on the SCL algorithm under the AWGN channel is proposed to against the problem of high complexity and large delay for decoding system caused by SCL decoding. It is also enhanced by combining with the establishment threshold of the decoder tree nodes to reduce the number of tree nodes' split. Under the AWGN channel, the bit error rate and the frame error rate performance of the virtual 2, 4 and 8 bit channels are simulated. The simulation results show that the decoding performance is close to that of the traditional SCL algorithm when the default threshold S=30 in the case of a virtual 8-bit channel. The total number of nodes is reduced by 63.7% and the total number of additions is 17% of the 8-bit decision decoding algorithm. The proposed algorithm can effectively reduce the calculation complexity of decoding and the hardware storage. It is more suitable for hardware implementation and has a certain practical value.

Key words:  polar code  successive cancellation list algorithm  decoding complexity  threshold


刘士平, 马林华, 胡星, 黄天宇. 低复杂度极化码SCL译码算法[J]. 哈尔滨工业大学学报, 2018, 50(5): 160-165. DOI: 10.11918/j.issn.0367-6234.201706116.
LIU Shiping, MA Linhua, HU Xing, HUANG Tianyu. Alow complexity SCL decoding algorithm for polar codes[J]. Journal of Harbin Institute of Technology, 2018, 50(5): 160-165. DOI: 10.11918/j.issn.0367-6234.201706116.
基金项目 国家自然科学基金资助项目(61472442);综合业务网国家重点实验室(西安电子科技大学)开放研究课题资助项目(INS-15-13);航空科学基金资助项目(20155896025) 作者简介 刘士平(1994—),男,博士研究生;
马林华(1965—),男,教授,博士生导师 通信作者 马林华,Land_max@126.com 文章历史 收稿日期: 2017-06-19



Contents            -->Abstract            Full text            Figures/Tables            PDF


低复杂度极化码SCL译码算法
刘士平, 马林华, 胡星, 黄天宇    
空军工程大学 航空航天工程学院,西安 710038

收稿日期: 2017-06-19
基金项目: 国家自然科学基金资助项目(61472442);综合业务网国家重点实验室(西安电子科技大学)开放研究课题资助项目(INS-15-13);航空科学基金资助项目(20155896025)
作者简介: 刘士平(1994—),男,博士研究生;
马林华(1965—),男,教授,博士生导师
通信作者: 马林华,Land_max@126.com


摘要: 极化码的串行抵消列表(SCL)译码的逐次逐比特进行判决过程与路径度量值的计算筛选过程是整个译码系统复杂度与延迟的主要来源.在分析现有SCL及多比特判决译码的优缺点基础上,针对SCL译码造成的译码系统复杂度高和延时大的问题,将每组多比特码字(多个独立信道)视作一个整体,并在译码时通过信道合成构造为一个虚拟多比特信道,从而可以对多比特码字进行同步传输并根据信道递归公式进行同步判决译码.由此基于SCL译码的码树构造提出一种构造多比特虚拟信道的SCL译码算法,并结合设置译码码树节点阈值减少码树节点分裂次数的方法进一步增强了上述算法.在AWGN信道下的分别对虚拟2、4和8比特信道SCL译码的误码率及误帧率性能进行仿真.仿真结果表明在虚拟8比特信道情况下,预设阈值S=30时的译码性能接近传统SCL算法,且总历经节点数降低了63.7%,总加法次数是8比特同时判决译码算法的17%.此算法降低了译码算法的计算复杂度及硬件存储复杂度,更适合于硬件实现, 具有一定的实用价值.
关键词: 极化码    串行抵消列表译码    译码复杂度    阈值    
Alow complexity SCL decoding algorithm for polar codes
LIU Shiping, MA Linhua, HU Xing, HUANG Tianyu    
Institute of Aeronautics and Astronautics Engineering, Air Force Engineering University, Xi'an 710038, China


Abstract: The successive cancellation list (SCL) decoding algorithm of polar code is of bitwise, which is the main source of the decoding complexity and latency. In view of the analysis of the advantages and disadvantages of the existing SCL and multi-bit decision decoding, each set of multi-bit code words (multiple independent channels) is treated as a whole to construct a multi-bit virtual channel, which can transmit multi-bit code words and decode synchronously. Therefore a construction of multi-bit virtual channel based on the SCL algorithm under the AWGN channel is proposed to against the problem of high complexity and large delay for decoding system caused by SCL decoding. It is also enhanced by combining with the establishment threshold of the decoder tree nodes to reduce the number of tree nodes' split. Under the AWGN channel, the bit error rate and the frame error rate performance of the virtual 2, 4 and 8 bit channels are simulated. The simulation results show that the decoding performance is close to that of the traditional SCL algorithm when the default threshold S=30 in the case of a virtual 8-bit channel. The total number of nodes is reduced by 63.7% and the total number of additions is 17% of the 8-bit decision decoding algorithm. The proposed algorithm can effectively reduce the calculation complexity of decoding and the hardware storage. It is more suitable for hardware implementation and has a certain practical value.
Key words: polar code    successive cancellation list algorithm    decoding complexity    threshold    
目前,香农在有噪信道编码理论中指出,存在可以达到香农限的码字[1].极化码(Polar Code)是2008年由E. Arikan提出的一种基于信道极化定理进行设计的新型信道编码,是第一种能够通过严格的数学方法证明达到信道容量的构造性编码方案[2-5].因此极化编码激起了各国学者的广泛研究并开始应用于各个领域,拥有良好的发展前景.

串行抵消(SC)译码算法是极化码的基本译码算法, 当码长趋近于无穷时,可获得很好的渐进性能[2].但是在中短码长时,SC译码的性能并不理想,相同条件下的编码增益不如Turbo码及LDPC码[6].为提高SC译码性能,文献提出了串行抵消列表(SCL)算法,增加了每层译码搜索宽度, 降低了误码率.在此基础上,文献[9-10]提出基于对数似然比(LLR)的SCL译码算法,对译码码树中每一个节点的判决计算次数缩减了一半,降低了计算及硬件存储复杂度.但上述两种方法都是逐比特串行译码,运算量及运算复杂度随着码长的增加呈指数增长.针对此类问题,文献[11-13]中多比特译码算法,将传统SCL的单比特译码增强为多比特同时译码.此算法尽管减少了译码器判决次数降低了译码延时,但由于其译码算法结构,随着每次译码的比特数增加,运算次数及硬件存储复杂度都会随之增加.基于上述问题,本文提出了基于构造虚拟多比特信道的SCL译码算法.将传统的多个单比特信道构造成多比特虚拟信道,降低了译码运算及存储复杂度,并结合文献[14]中提出的基于信道可靠度的阈值判决的译码方案,进一步增强了此算法,提高译码效率.

1 背景算法描述 1.1 极化码SCL译码算法(N, K)极化编码输入的N比特信息u1N=(u1, u2, …, un)由K比特的信息比特和(N-K)位冻结比特(一般为全0)两部分组成.根据信道极化理论,分别对N=2n个独立的二进制离散无记忆信道(B-DMC)W进行信道合成及信道分解的操作后,可得到N个独立等价子信道WN(i)i =1, 2, …, N.其中一部分子信道的对称容量I(WN(i))接近于1,视作无噪信道,用于传输信息比特; 而另一部分子信道的对称容量I(WN(i))接近于0,视作全噪信道,用于传输冻结比特.

根据文献[2]中提到的有关信道可靠性的计算方法,u1N中的K比特消息通过可靠性高的子信道进行传输,而(N-K)位冻结比特通过可靠度低的子信道进行传输.其次,依据编码式(1),完成极化编码,其中GN为生成矩阵.

$x_1^N = u_1^N{\mathit{\boldsymbol{G}}_N} = {\mathit{\boldsymbol{u}}_A}{\mathit{\boldsymbol{G}}_N}\left( \mathit{\boldsymbol{A}} \right) \oplus {\mathit{\boldsymbol{u}}_{{A^c}}}{\mathit{\boldsymbol{G}}_N}\left( {{\mathit{\boldsymbol{A}}^c}} \right).$ (1)

式中:A为信息信道下标的集合; Ac为其补集即冻结信道的下标集合; 符号$ \oplus $为模二加运算.生成矩阵GN

${\mathit{\boldsymbol{G}}_N} = {\mathit{\boldsymbol{B}}_N}{\mathit{\boldsymbol{F}}^{ \otimes n}}.$ (2)

式中:符号$ \otimes $为张量积运算; ${\mathit{\boldsymbol{B}}_N} = {\mathit{\boldsymbol{R}}_N}({\mathit{\boldsymbol{I}}_2} \otimes {\mathit{\boldsymbol{B}}_{N/2}}) $,实现索引位逆转操作; RN是比特翻转矩阵即实现元素的奇偶分离重排.

由此可见,极化码也是线性分组码的一种,(N, K)极化码的码率R=N/K.极化编码可视作根据信道可靠度的计算选取合适的信道,确定编码参数(N, K, A, uAc)的GN陪集编码过程.

在仿真实验中,极化码的码字序列经二进制相移键控(BPSK)调制为符号x1N∈{-1, 1}N的序列,经过AWGN信道后得到的接收序列y1N=(y1, y2, …, yN).采用串行抵消(SC)译码可看作逐个依次地根据接收端计算出的对数似然比(LLR)进行硬判决. ui的似然函数为$W_n^{(i)}(y_1^N, \hat u_1^{i - 1}|{u_i}) $,LLR为

$L\left( {{u_i}} \right) = \ln \left( {\frac{{W_n^{\left( i \right)}\left( {y_1^N,\hat u_1^{i - 1}\left| 0 \right.} \right)}}{{W_n^{\left( i \right)}\left( {y_1^N,\hat u_1^{i - 1}\left| 1 \right.} \right)}}} \right).$ (3)

式中$ \hat u_1^{i - 1}$是接收端对前i-1个比特的估计序列.根据上述得到的LLR值结合硬判决即可完成对ui的估计,具体判决依据如下:

${{\hat u}_1} = \left\{ \begin{array}{l}0,i \in \mathit{\boldsymbol{A}},L\left( {{u_i}} \right) > 0;\\1,i \in \mathit{\boldsymbol{A}},L\left( {{u_i}} \right) < 0;\\{u_i},i \in {\mathit{\boldsymbol{A}}^c}.\end{array} \right.$ (4)

在上述基于LLR的SC译码算法的基础上,根据译码码树,若在每次按位搜索时将搜索宽度扩大为L,即保留L条候选路径.然后基于L条路径依次扩展,直到叶子节点即最后一层.从最后一层保留的L条路径中选择一个路径度量值最大的路径作为译码输出,此种方法即为基于LLR的SCL译码算法.因此,SCL算法可看作是多路径同时进行的SC译码算法.

1.2 多比特判决译码算法在极化码的译码算法中,路径度量的判决筛选过程是整个系统延迟的主要来源.无论是SC译码还是SCL译码算法都依据硬判决式(4)逐次逐位进行判决.当码长增加时,该算法的系统延迟增加,降低了译码效率.针对此问题,文献[13]提出了多比特判决的译码算法(2Kb-rSCL算法),译码时可每次同时计算出M比特消息的路径度量值进行判决筛选.如图 1所示,当M=1时等价于传统SCL译码,当M时,该算法即为最大似然(ML)译码.

Figure 1
图 1 传统单比特及多比特译码示意 Figure 1 Traditional singal-bit and multi-bit decoding schematic


假设W: XY为一个B-DMC信道,则N个并列单比特信道为WN(i): XYN× Xi-1, 1≤iN.在多比特同步译码算法中,M(M=2K)比特码字的路径可靠性度量过程视作由M个独立的单比特信道进行信道合成得到,即$W_{N, M}^{(j)}:{\rm{ }}\mathit{\boldsymbol{X}}{^M} \to {\rm{ }}\mathit{\boldsymbol{Y}}{^N} \times {\rm{ }}\mathit{\boldsymbol{X}}{^{jM - M}}, 1 \le j \le {\rm{ }}\frac{N}{M} $,其转移概率为$W_{N, M}^{(j)}(y_1^N, x_1^{jM - M}|x_{jM - M + 1}^{jM}) $.基于信道与信道的输出码字u1N的独立性,上述算法可描述为:当M比特码字同时译码判决时,信道转移概率$W_{N, M}^{(i)}(y_1^N, \hat u_1^{iM - M}|u_{iM - M + 1}^{iM}) $转化成M个独立的单比特信道之积,即

$\begin{array}{l}W_{N,M}^{\left( i \right)}\left( {y_1^N,\hat u_1^{iM - M}\left| {u_{iM - M + 1}^{iM}} \right.} \right) = \prod\limits_{j = 0}^{M - 1} {W_{\frac{N}{M}}^{\left( i \right)}\left( {y_{j\frac{N}{M} + 1}^{\left( {j + 1} \right)\frac{N}{M}},} \right.} \\\;\;\;\;\left. {\hat \alpha _{1 + j\frac{N}{M}}^{\left( {i + 1} \right) + j\frac{N}{M}}\left| {{\alpha _{i + j\frac{N}{M}}}} \right.} \right).\end{array}$ (5)

式中$ u_{iM - M + 1}^{iM} = \left( {{\alpha _i}, {\alpha _{i + \frac{N}{M}}}, \ldots, {\alpha _{i + N - \frac{N}{M}{\rm{ }}}}} \right){\mathit{\boldsymbol{G}}_M}$; 生成矩阵${\mathit{\boldsymbol{G}}_M} = {\rm{ }}{\mathit{\boldsymbol{B}}_M}\mathit{\boldsymbol{F}}{^{ \otimes K}} $.根据得到的似然值$W_{N, M}^{(j)}(y_1^N, \hat u_1^{iM - M}|u_{iM - M + 1}^{iM}) $,结合SCL译码算法,通过路径度量模块(MCU)进行筛选保留L条候选路径以及通迫零算法模块(ZFU)将冻结位至“0”.通过增加了每次译码比特数从而提高译码效率.

2 基于虚拟多比特信道的SCL译码多比特译码算法的核心是对似然值$W_{N, M}^{(i)}(y_1^N, \hat u_1^{iM - M}|u_{iM - M + 1}^{iM}) $的计算.在对图 1(b)中每组M比特码字进行译码时,上述算法将其转化为传统单比特信道似然值之积.然而该算法只是通过路径度量上的多位判决降低了译码延迟,随着码长以及M的增加,计算复杂度以及进行路径度量判决(MCU)的查表复杂度也随之增加.并没有从信道转移概率这个根本问题上解决多比特运算的运算及存储复杂度.

将每组M比特码字即M个独立信道视作一个整体.在进行译码判决时,可将M个独立单比特信道通过信道合成构造为一个独立的M比特虚拟信道.此虚拟信道可以同时对M比特码字进行传输并根据信道递归公式进行同步判决,而不是M个单比特信道之积.这个虚拟信道为$W_N^{(i)}:\mathit{\boldsymbol{ \boldsymbol{\varLambda} }} \to {\mathit{\boldsymbol{Y}}^N} \times {\mathit{\boldsymbol{ \boldsymbol{\varLambda} }}^{i - 1}}, 1 \le i \le \frac{N}{M} $,信道转移概率为$W_n^{(i)}(y_1^N, \lambda _1^{i - 1}|{\lambda _i}), {\lambda _i} = ({u_{jM - M + 1}}, \ldots, {u_{jM}}) $.根据虚拟信道的性质,对上述转移概率$W_{N, M}^{(i)}(y_1^N, \hat u_1^{iM - M}|u_{iM - M + 1}^{iM}) $的计算可通过信道递归计算有进一步地简化,并有下述定理1.

定理1??假设u1N中每一个比特ui分别独立的等概率的判决为0或1,那么给定0 < mnN=2nM=2m,对任意的1 < km,0≤ω < n,令K=2kΩ=2ω,并且$0 \le i < \frac{{2\mathit{\Omega} }}{K} $,则K比特虚拟信道W2Ω, K(i+1)可以由两个独立的虚拟$ \frac{K}{2}$比特信道WΩ, K/2(i+1)通过单步信道变换得到,即

$\left( {W_{\mathit{\Omega },K/2}^{\left( {i + 1} \right)},W_{\mathit{\Omega },K/2}^{\left( {i + 1} \right)}} \right) \mapsto W_{2\mathit{\Omega },K}^{\left( {i + 1} \right)}$

且信道转移概率满足

$\begin{array}{*{20}{c}}{W_{2\mathit{\Omega },K}^{\left( {i + 1} \right)}\left( {y_1^{2\mathit{\Omega }},u_1^{iK}\left| {u_{iK + 1}^{iK + K}} \right.} \right) = }\\{W_{\mathit{\Omega },K/2}^{\left( {i + 1} \right)}\left( {y_1^\mathit{\Omega },u_{1,o}^{iK} \oplus u_{1,e}^{iK}\left| {u_{iK + 1,o}^{iK + K} \oplus u_{iK + 1,e}^{iK + K}} \right.} \right) \cdot }\\{W_{\mathit{\Omega },K/2}^{\left( {i + 1} \right)}\left( {y_{\mathit{\Omega } + 1}^{2\mathit{\Omega }},u_{1,e}^{iK}\left| {u_{iK + 1,e}^{iK + K}} \right.} \right).}\end{array}$ (6)

证明??由条件概率公式Pr(B |A)=Pr(A, B)/Pr(A)可得

$W_{2\mathit{\Omega}, K}^{(i + 1)}\left( {y_1^{2\mathit{\Omega} }, u_1^{iK}|u_{iK + 1}^{iK + K}} \right) = \frac{{W_{2\mathit{\Omega}, 1}^{(iK + K)}\left( {y_1^{2\mathit{\Omega} }, u_1^{iK + K - 1}|{u_{iK + K}}} \right)}}{{{\rm{Pr}}\left( {u_{iK + 1}^{iK + K - 1}|{u_{iK + K}}} \right)}} $由于每一个比特ui分别独立的等概率的判决为0或1,则

$\Pr \left( {u_1^{iK + K - 1}\left| {{u_{iK + K}}} \right.} \right) = \Pr \left( {u_1^{iK + K - 1}} \right) = {2^{ - \left( {K - 1} \right)}}$

因此

$\begin{array}{*{20}{c}}{W_{2\mathit{\Omega },K}^{\left( {i + 1} \right)}\left( {y_1^{2\mathit{\Omega }},u_1^{iK}\left| {u_{iK + 1}^{iK + K}} \right.} \right) = }\\{{2^{\left( {K - 1} \right)}}W_{2\mathit{\Omega },K}^{\left( {i + 1} \right)}\left( {y_1^{2\mathit{\Omega }},u_1^{iK + K - 1}\left| {{u_{iK + K}}} \right.} \right).}\end{array}$ (7)

根据文献[2]中的式(23)可得

$\begin{array}{l}W_{2\mathit{\Omega },1}^{\left( {iK + K} \right)}\left( {y_1^{2\mathit{\Omega }},u_1^{iK + K - 1}\left| {{u_{iK + K}}} \right.} \right) = \frac{1}{2}W_{\mathit{\Omega },1}^{\left( {\frac{{iK + K}}{2}} \right)}\left( {y_1^\mathit{\Omega },u_{1,o}^{iK + K - 2} \oplus } \right.\\\;\;\;\left. {u_{1,e}^{iK + K - 2}\left| {{u_{iK + K - 1}} \oplus {u_{iK + K}}} \right.} \right) \cdot \\\;\;\;W_{\mathit{\Omega },1}^{\left( {\frac{{iK + K}}{2}} \right)}\left( {y_{\mathit{\Omega } + 1}^{2\mathit{\Omega }},u_{1,e}^{iK + K - 2}\left| {{u_{iK + K}}} \right.} \right)\end{array}$ (8)

可以推出

$\begin{array}{l}W_{\mathit{\Omega },1}^{\left( {\frac{{iK + K}}{2}} \right)}\left( {y_1^\mathit{\Omega },u_{1,o}^{iK + K - 2} \oplus u_{1,e}^{iK + K - 2}\left| {{u_{iK + K - 1}}} \right. \oplus {u_{iK + K}}} \right) = \\{2^{ - \left( {\frac{K}{2} - 1} \right)}}W_{\mathit{\Omega },K/2}^{\left( {i + 1} \right)}\left( {y_1^\mathit{\Omega },u_{1,o}^{iK} \oplus u_{1,e}^{iK}\left| {u_{iK + 1,o}^{iK + K}} \right. \oplus u_{iK + 1,e}^{iK + K}} \right)\end{array}$ (9)

同理

$\begin{array}{*{20}{c}}{W_{\mathit{\Omega },1}^{\left( {\frac{{iK + K}}{2}} \right)}\left( {y_{\mathit{\Omega } + 1}^{\mathit{2\Omega }},u_{1,e}^{iK + K - 2}\left| {{u_{iK + K}}} \right.} \right) = }\\{{2^{ - \left( {\frac{K}{2} - 1} \right)}}W_{\mathit{\Omega },K/2}^{\left( {i + 1} \right)}\left( {y_{\mathit{\Omega } + 1}^{\mathit{2\Omega }},u_{1,e}^{iK}\left| {u_{1K + 1,e}^{iK + K}} \right.} \right)}\end{array}$ (10)

综上,由式(7)~(10),即可推出(6),证毕.

当利用上述多比特虚拟信道算法进行译码时,译码码树每层对K位信息uiK+1iK+K的判决有2K条候选路径,因此基于LLR的二元判决模式不再实用.针对此问题,结合上述基于阈值度量的算法采取了对数似然(LL)结构,由上述定理1可得到每条候选路径的LL值即$ {\rm{log}}W_{2\mathit{\Omega}, K}^{(i + 1)}\left( {y_1^{2\mathit{\Omega}}, u_1^{iK}|u_{iK + 1}^{iK + K}} \right)$.

随着虚拟信道数K的增加,译码判决候选路径数2K也呈指数增长.为解决此问题,引入AWGN信道下基于高斯近似的低分裂次数优化算法,进一步增强上述译码算法,降低路径度量运算以及存储复杂度.文献[14]提出的ESR-SCL译码算法中,引入Pe(ui)为在译码过程中节点ui的平均错误概率,即在子信道$W_n^{(i)}(y_1^N, \hat u_1^{i - 1}|{u_i}) $中,当u1i-1判决全部正确时ui被错误译码判决的概率.则1-Pe(ui)为ui节点的译码可靠性度量值即置信度,并作为译码节点的阈值. ${P_e}\left( {{u_i}} \right) = Q(\sqrt {E\left[{L\left( {{u_i}} \right)} \right]/2} {\rm{ }}) $其中$ Q\left( x \right) = \frac{1}{{\sqrt {2{\rm{ \mathsf{ π} }}} }}\int_x^{ + \infty } {{e^{ - \frac{{{t^2}}}{2}}}{\rm{d}}t} $.根据阈值的预先计算提出了上述算法的核心,即只有当$W_N^{(i)}(y_1^N, \hat u_1^{i - 1}|{u_i} = 0) > 1 - {P_e}({u_i}) $$W_N^{(i)}(y_1^N, \hat u_1^{i - 1}|{u_i} = 1) > 1 - {P_e}({u_i}) $时,译码码树在此层不分裂直接判为0或1,从而减少了历经节点数,降低了复杂度.

由高斯近似可知,当任意i∈[1, N]时,ui取1或0的L(ui)值具有相同的概率密度函数[15].因此K比特信息uiK+1iK+K的置信度为(1-Pe(uj))K,其中iK+1≤jiK+K.然而由上述公式,Pe(ui)中变量L(ui)的计算需要得到已经译码判决出的u1i-1.根据多比特虚拟信道的译码码树构造可知,在对${u_{iK + 1}^{iK + K}} $进行判决时,$\hat u_1^{iK} $的判决结果是已知的,故只有uiK+1的LLR值L(uiK+1)可以通过计算得到.因此上述计算中取j=iK+1,则${u_{iK + 1}^{iK + K}} $置信度可定义为(1-Pe(uiK+1))K.其中:

$\begin{array}{l}{P_e}\left( {{u_{iK + 1}}} \right) = {\rm{Q}}\left( {\sqrt {E\left[ {L\left( {{u_{iK + 1}}} \right)} \right]/2} } \right),\\\;\;\;L\left( {{u_{iK + 1}}} \right) = \log \frac{{W_N^{\left( i \right)}\left( {y_1^N,\hat u_1^{iK}\left| 0 \right.} \right)}}{{W_N^{\left( i \right)}\left( {y_1^N,\hat u_1^{iK}\left| 1 \right.} \right)}}.\end{array}$ (11)

结合虚拟信道的信道转移概率LL值以及上述给出的可靠性度量的阈值算法,定义通过上述虚拟信道传输算法译码时只有当${\rm{log}}\left( {W_{2\mathit{\Omega}, K}^{(i + 1)}\left( {y_1^{2\mathit{\Omega} }, u_1^{iK}|u_{iK + 1}^{iK + K} = \psi _{iK + 1}^{iK + K}} \right)} \right) > K{\rm{log}}(1 - {P_e}({u_{iK + 1}})) $时,译码码树在此层不分裂直接判为$\psi _{iK + 1}^{iK + K} $,其中$\psi _{iK + 1}^{iK + K} $$u_{iK + 1}^{iK + K} $2K种可能码字中任意一个码字.其余情况下结合多比特虚拟信道SCL译码算法正常分裂.即

$\hat u_{iK + 1}^{iK + K} = \left\{ \begin{array}{l}\psi _{iK + 1}^{iK + K},\log \left( {W_{2\mathit{\Omega },K}^{\left( {i + 1} \right)}\left( {y_1^{2\mathit{\Omega }},u_1^{iK}\left| {u_{iK + 1}^{iK + K} = \psi _{iK + 1}^{iK + K}} \right.} \right)} \right) \ge \\\;\;\;\;\;\;\;\;K\log \left( {1 - {P_e}\left( {{u_{iK + 1}}} \right)} \right)\\{\rm{split}},{\rm{else}}.\end{array} \right.$ (12)

基于可靠性度量的计算可以推断,在路径正常分裂后,若出现${\rm{log}}\left( {W_{2\mathit{\Omega}, K}^{(i + 1)}\left( {y_1^{2\mathit{\Omega} }, u_1^{iK}|u_{iK + 1}^{iK + K} = \psi _{iK + 1}^{iK + K}} \right)} \right) \le K{\rm{log}}(1 - {P_e}({u_{iK + 1}})) $则该路径默认为错误路径不再保留或扩展,直接删除.

基于上述算法的译码路径可分为两类,第一类是基本上不分裂,此类路径以很高的概率包含正确的码字; 第二类是进行频繁地分裂,此类路径译码正确的概率相对较低.因此,定义变量Sl(i)为第l条路径在到达译码树第i层时中间经历的分裂次数.即假如第l条路径在第i层扩展到第i+1层时没有分裂,那么Sl(i+1)=Sl(i)+1.正确的译码路径很少分裂,而错误的译码路径会频繁地分裂,定义阈值S.Sl(i)≥S,则此路径很大可能为正确译码,予以保留,其余Sl(i) < S的路径予以抛弃,从而减少保留的路径数量,以较慢的速度达到SCL译码中定义的L值.

综上所述,提出基于阈值判决的多比特虚拟信道SCL译码算法,描述步骤如下:

第1步:根据式(6)的计算,得到译码输出每组K比特信息的LL值,共计2k个值.

第2步:将上一步得到的2k个LL值依据式(12)和(13)对其进行判决,筛选出不分裂的、被删除的以及正常分裂的子节点,Sl(i)值也随时保持更新.

第3步:在第2步筛选基础上,若此时保留路径数量P>L,则删除Sl(i) < S的路径; 若所有路径Sl(i) < S,或删除Sl(i) < S的路径后P仍大于L,则保留剩余路径中路径度量值最大的L条.

第4步:若$i < \frac{N}{K} $,则令i=i+1并重复第一步至第3步; 若$i = \frac{N}{K} $,则在最后一层保留路径中取出路径度量值最大的一条路径作为译码输出.

3 仿真结果及性能分析在AWGN信道下选取码长为1024,码率为1/2的(1024, 512)极化码,经过二进制频移键控(BPSK)调制,分别对上述多比特虚拟信道的SCL译码及基于不同阈值判决的虚拟信道SCL译码算法进行1 000次蒙特卡洛仿真.

图 2为(1024, 512)极化码的传统的SCL译码算法与K分别取2、4、8的多比特虚拟信道SCL译码算法的误码率及误帧率的对比.从FER/BER曲线可以看出,当虚拟信道的比特数增大时,该算法的译码性能由略微的增加但是并不明显,性能与传统的SCL译码算法基本保持一致.

Figure 2
图 2 (1 024, 512)极化码不同参数下SCL译码算法译码性能的对比 Figure 2 The comparison of SCL decoding performance of (1024, 512) polar code under different parameters


以4比特同时译码为例,由式(6)对$ {\rm{log}}(W_{N, 4}^{(j)}(y_1^N, \hat u_1^{4i - 4}|u_{4i - 4 + 1}^{4i}))$的计算共需要经历(24+2×22)共计24次加法运算,而由式(5)对${\rm{log}}(W_{N, 4}^{(j)}(y_1^N, \hat u_1^{4i - 4}|u_{4i - 4 + 1}^{4i})) $的计算需进行(24×(4-1))共计48次加法运算,以此类推,得到如下复杂度对比表.

由图 2和表 1的分析可得出,上述多比特虚拟信道的译码算法与文献[13]中多比特判决译码相比,在K=4与K=8时可分别减少50%及83%的加法运算次数.在降低了译码计算复杂度的同时,译码性能与传统SCL算法基本保持一致.

表 1
表 1 计算${\rm{log}}({\rm{ }}W_{2\mathit{\Omega}, K}^{(i + 1)}(y_1^{2\mathit{\Omega} }, u_1^{iK}|u_{iK + 1}^{iK + K})) $所需加法次数 Table 1 Numbers of additions to calculate ${\rm{log}}({\rm{ }}W_{2\mathit{\Omega}, K}^{(i + 1)}(y_1^{2\mathit{\Omega} }, u_1^{iK}|u_{iK + 1}^{iK + K})) $ K值 虚拟信道算法 2Kb-SCL算法 本文算法优化率%

K=2 4 4 0

K=4 24 48 50

K=8 304 1 792 83

K=16 66 144 1 179 648 94.4



表 1 计算${\rm{log}}({\rm{ }}W_{2\mathit{\Omega}, K}^{(i + 1)}(y_1^{2\mathit{\Omega} }, u_1^{iK}|u_{iK + 1}^{iK + K})) $所需加法次数 Table 1 Numbers of additions to calculate ${\rm{log}}({\rm{ }}W_{2\mathit{\Omega}, K}^{(i + 1)}(y_1^{2\mathit{\Omega} }, u_1^{iK}|u_{iK + 1}^{iK + K})) $


图 3、4为当(1 024, 512)极化码多比特虚拟信道译码时的K值分别取4、8时,在不同译码码树节点阈值S下的译码性能曲线.由不同阈值下的对比可知,根据阈值判决的虚拟比特信道译码由于减少了码树分裂次数即历经节点数导致译码性能略差.当所选取的阈值越大越接近传统SCL译码或虚拟比特信道SCL译码,随着阈值S的减小,译码性能也随之下降.在阈值S达到30时,性能略优于文献[14]中的ESR-SCL译码算法,S=20和25时译码性能较ESR-SCL算法略差.

Figure 3
图 3 (1 024, 512)极化码虚拟4比特信道译码时的不同阈值S下的译码性能 Figure 3 The decoding performance of (1 024, 512) polar code at different thresholds S under 4-bit virtual channel


Figure 4
图 4 (1 024, 512)极化码虚拟8比特信道译码时的不同阈值S下的译码性能 Figure 4 The decoding performance of (1 024, 512) polar code at different thresholds S under 8-bit virtual channel


由上述实验得到,在K=8、信噪比为2.5 dB时误码率可达到10-5.图 5为当信噪比Eb/N0=2.5 dB时,无预设阈值与不同阈值下虚拟8比特信道SCL译码码树每层历经的节点总数.从图 3、4观察对比可以发现,设定阈值后,虽然译码性能略有下降,但每层历经的节点数量明显减少,在图 5中可观察不同阈值下历经节点数减少的具体情况.表 3记录了信噪比Eb/N0=2.5 dB时不同参数下1 000次实验译码树平均历经的节点总数.

Figure 5
图 5 (1 024,512)极化码虚拟8比特信道译码时的不同阈值S下的平均历经节点总数 Figure 5 The decoding performance of (1 024, 512) polar code at different thresholds S under 8-bit virtual channel


表 2
表 2 不同阈值参数平均历经节点总数 Table 2 The average experience nodes under different thresholds S值 无 30 25 20

历经节点数 1007 365.13 319.00 289.07



表 2 不同阈值参数平均历经节点总数 Table 2 The average experience nodes under different thresholds


仿真实验表明,当阈值S=30时,虚拟8比特信道SCL译码性能接近于传统SCL译码,但是总历经节点数降低了63.7%,所需加法计算次数与8比特同时判决译码算法相比减少了83%.降低了译码算法的复杂度及历经节点次数,降低了硬件存储复杂度.

4 结论本文回顾了极化码编码及SCL译码的原理,分析论证了现有多比特判决译码的优缺点.在此基础上提出了基于构造多比特虚拟信道的SCL译码算法并给出了多比特虚拟信道的递归公式,并结合上述算法的译码码树的构造结合了基于节点阈值判决的减少译码树分裂次数的算法,最后对上述算法进行了仿真实验与复杂度分析.结果表明,采用虚拟8比特信道SCL译码,取预设阈值S=30时译码性能接近传统SCL算法,在信噪比Eb/N0=2.5 dB时误码率可达到10-5.总历经节点数降低了63.7%,总加法计算次数与8比特同时判决译码算法相比减少了83%,降低了译码算法的计算复杂度及硬件存储复杂度,因此,该算法更适合于硬件实现, 更有应用价值.


参考文献
[1]SHANNON C E. A mathematical theory of communication[J].Bell System Technical Journal, 1948, 19(4): 271-285.DOI: 10.1002/j.1538-7305.1948.tb01338.x


[2] ARICAN E. Channel polarization: a method for constructing capacity achieving codes for symmetric binary-input memoryless channels[J]. IEEE Trans. Inform Theory, 2009, 55(7): 3051-3073. DOI: 10.1109/TIT.2009.2021379.


[3]ARICAN E. A performance comparison of polar codes and Reed-Muller codes[J].IEEE Communications Letters, 2008, 12(6): 447-449.DOI: 10.1109/LCOMM.2008.080017


[4] ARICAN E. Channel combining and splitting for cutoff rate improvement[J]. IEEE Trans. Inf. Theory, 2006, 52(2): 628-639. DOI: 10.1109/ISIT.2005.1523420.


[5] ARICAN E, KIM H, MARKARIAN G. Performance of short polar codes under ML decoding[C] // Paul Cunningham and Miriam Cunningham(Eds), ⅡMC International Information Management Corporation, 2009: 1-6.


[6]CHEN Kai, NIU Kai and LIN Jiaru. Improved successive cancellation decoding of polar codes[J].IEEE Transaction on Communications, 2013, 61(8): 3100-3107.DOI: 10.1109/TCOMM.2013.070213.120789


[7]I. TAL, A. VARDY. List decoding of polar codes[J].IEEE Transactions on Information Theory, 2011, 61(5): 2213-2226.DOI: 10.1109/TIT.2015.2410251


[8]LI Bin, SHEN Hui, TSE D. An adaptive successive cancellation list decoder for polar codes with cyclic redundancy check[J].IEEE Communications Letters, 2012, 16(12): 2044-2047.DOI: 10.1109/LCOMM.2012.111612.121898


[9] A BALATSOUKAS-STIMMING, A RAYMOND W GROSS, A BURG. Hardware architecture for list successive cancellation decoding of polar codes[J]. IEEE Trans. Circuits Syst. Ⅱ Express Briefs, 2014, 61(8): 609-613. DOI: 10.1109/TCSⅡ.2014.2327336.


[10] A. BALATSOUKAS-STIMMING, M. B. PARIZI, A. BURG. LLR-based successive cancellation list decoding of polar codes[J]. IEEE Trans. Signal Process, 2015, 63(19): 5165-5179. DOI: 10.1109/ICASSP.2014.6854333.


[11] LI Bin, SHEN Hui, TSE D. Parallel decoders of polar codes. sep. 2013[EB/OL]. 4 sep. 2013, Available: http://arxiv.org/abs/1309.1026.


[12] LI Bin, SHEN Hui, TSE D, TONG Wen. Low-latency polar codes via hybrid decoding[J]. Proc. 8th Int. Symp. Turbo Codes and Iterative Inf. Process, Aug. 2014: 2223-2227. DOI: 10.1109/ISTC.2014.6955118.


[13] YUAN B, PHARHI K K. Low-latency successive-cancellation list decoders for polar codes with multi-bit decision[J]. IEEE Trans. Very Large Scale Integr, 2015, 23(10): 2268-2280. DOI: 10.1109/TVLSI.2014.2359793.


[14]ZHANG Zhaoyang, ZHANG Liang, WANG Xianbin, ZHONG Caijun. A split-reduced successive cancellation list decoder for polar codes[J].IEEE Journal on Selected Areas in Communications, 2016, 34(2): 292-302.DOI: 10.1109/JSAC.2015.2504321


[15] CHUNG S Y, RICHARSON T J, URBANKE R L. Analysis of sum-product decoding of low-density parity-check codes using a gaussian approximation[J] IEEE Trans. Inf. Theory, 2001, 47(2): 657-670. DOI: 10.1109/18.910580.



相关话题/系统 过程 计算 译码 航空航天

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 工业纯钛金属织构标准极图的计算及分析
    工业纯钛金属织构标准极图的计算及分析陈亮维,刘状,虞澜,胡劲,易健宏(昆明理工大学材料科学与工程学院,昆明650093)摘要:工业纯钛中的金属织构会引起各向异性,获得织构信息及分析其演变规律对钛材加工与应用非常重要.本文利用单晶钛的晶体结构数据、乌氏网、极图与织构的定义,建立了纯钛的织构与特定晶面极 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 热轧过程中高纯钴微观组织及织构演变
    热轧过程中高纯钴微观组织及织构演变李震1,宋克睿1,韩彦鹏1,肖柱1,贺昕2,陈志永1(1.中南大学材料科学与工程学院,长沙410083;2.北京有色金属研究总院有研亿金新材料股份有限公司,北京102200)摘要:金属钴具有同素异构转变特性。为探究热轧工艺对高纯钴的微观组织及织构演变规律的影响,对纯 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • AZ31B镁合金带材热轧过程组织均匀性及性能研究
    AZ31B镁合金带材热轧过程组织均匀性及性能研究曹东东1,2,梅瑞斌1,2,包立1,侯铮1,黄芸1(1.东北大学秦皇岛分校资源与材料学院,河北秦皇岛066004;2.东北大学材料科学与工程学院,沈阳110819)摘要:本文开展了变形温度为300、350、400℃和总压下率分别为15%、30%、45% ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 抗癌药物作用预测计算方法的研究现状与展望
    抗癌药物作用预测计算方法的研究现状与展望顾兆伟1,张立忠2,刘晓峰3,谭先4(1.长春中医药大学附属第三临床医院脑病康复科,长春130000;2.长春市朝阳区清和社区卫生服务中心,长春130000;3.空军杭州特勤疗养中心康复理疗科,杭州310000;4.东北师范大学信息科学与技术学院,长春1300 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 液压轮毂马达辅助驱动系统控制策略实车验证
    液压轮毂马达辅助驱动系统控制策略实车验证曾小华,崔臣,张轩铭,宋大凤,李立鑫(汽车仿真与控制国家重点实验室(吉林大学),长春130025)摘要:为充分提升重型牵引车辆通过不良路面的能力,对国内某款重型牵引车在传统结构的基础上加装了前轴液压轮毂马达辅助驱动系统,并针对该混合动力系统,开发了工程化的控制 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 考虑高温水泄漏影响的凝水泵启动过程仿真分析
    考虑高温水泄漏影响的凝水泵启动过程仿真分析覃海波,金家善,倪何(海军工程大学动力工程学院,武汉430033)摘要:针对某船用给水机组凝水泵在备用状态和备用转换启动过程中,由于高温除氧水进入凝水泵及其吸入管路而在机组启动后发生汽蚀,导致凝水泵出口压力长时间达不到规定要求的问题,在基于差异演化算法辨识得 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 铁锰复合阴极MFC-EF耦合系统产电及降解RhB效能
    铁锰复合阴极MFC-EF耦合系统产电及降解RhB效能史珂1,2,赵庆良1,2,王维业2,王琨1,2(1.城市水资源与水环境国家重点实验室(哈尔滨工业大学),哈尔滨150090;2.哈尔滨工业大学环境学院,哈尔滨150090)摘要:为实现微生物燃料电池(MFC)微电的原位利用,结合电芬顿(EF)技术的 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 污水生物处理过程中溶解性有机氮分布和转化特征
    污水生物处理过程中溶解性有机氮分布和转化特征王小东1,陈明飞1,王子文1,王燕1,王硕1,2,3,李激1,2,3(1.江南大学环境与土木工程学院,江苏无锡214122;2.江苏省厌氧生物技术重点实验室(江南大学),江苏无锡214122;3.江苏省高校水处理技术与材料协同创新中心,江苏苏州215009 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 桥梁颤振临界风速的概率密度演化计算
    桥梁颤振临界风速的概率密度演化计算姜保宋,周志勇,唐峰(土木工程防灾国家重点实验室(同济大学),上海200092)摘要:针对桥梁结构自身特性以及外部环境的随机性(如刚度、质量、阻尼比、气动导数等因素)所造成的桥梁的颤振临界风速不确定,难以衡量桥梁颤振稳定性问题.将概率密度演化方法与桥梁颤振多模态耦合 ...
    本站小编 哈尔滨工业大学 2020-12-05
  • 桥梁伸缩缝跳车冲击荷载计算方法与模型实验
    桥梁伸缩缝跳车冲击荷载计算方法与模型实验丁勇1,2,王佩1,游玖昂1,诸葛萍1(1.宁波大学土木工程系,浙江宁波315211;2.桥梁工程结构动力学国家重点实验室(重庆交通科研设计院),重庆400067)摘要:为实测移动车辆对桥梁伸缩缝的冲击荷载,防止桥梁伸缩缝在这种冲击荷载作用下发生早期损坏,制作 ...
    本站小编 哈尔滨工业大学 2020-12-05