CN109347600B - 一种极化码译码方法 - Google Patents

一种极化码译码方法 Download PDF

Info

Publication number
CN109347600B
CN109347600B CN201811125637.8A CN201811125637A CN109347600B CN 109347600 B CN109347600 B CN 109347600B CN 201811125637 A CN201811125637 A CN 201811125637A CN 109347600 B CN109347600 B CN 109347600B
Authority
CN
China
Prior art keywords
function
update
log
matrix
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.)
Active
Application number
CN201811125637.8A
Other languages
English (en)
Other versions
CN109347600A (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.)
Henan University of Science and Technology
Original Assignee
Henan 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 Henan University of Science and Technology filed Critical Henan University of Science and Technology
Priority to CN201811125637.8A priority Critical patent/CN109347600B/zh
Publication of CN109347600A publication Critical patent/CN109347600A/zh
Application granted granted Critical
Publication of CN109347600B publication Critical patent/CN109347600B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045—Arrangements at the receiver end
    • H04L1/0054—Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • 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
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0056—Systems characterized by the type of code used
    • H04L1/0057—Block codes

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Theoretical Computer Science (AREA)
  • Error Detection And Correction (AREA)

Abstract

本发明提供一种极化码译码方法，在极化码SC译码方法的基础上，对f函数节点的计算进行简化处理，用折线逼近函数代替双曲正切函数及其反函数运算，折线逼近函数只需要加法和乘法运算。本发明有益效果：大大降低了运算复杂度，与查表法相比大大节约了内存资源，降低了硬件成本，解决了极化码在对数似然比域译码运算复杂度较高的问题。

Description

一种极化码译码方法
技术领域
本发明涉及编码技术领域，具体地说是一种极化码译码方法。
背景技术
极化码由土耳其毕尔肯大学
Figure BDA0001812314220000014
教授于2008年首次提出，
Figure BDA0001812314220000015
从理论上严格证明了在二进制对称离散无记忆信道(binary symmetric discrete memorylesschannel，B-DMC)下，极化码的码率可以达到香农极限，并且极化码有着低的编码和译码复杂度。极化码译码的特点是采用SC译码算法，按照自然顺序对发送比特进行逐级判决译码，先判决的比特作为可靠信息参加后译比特的判决。在判决时，对于冻结比特位部分可直接译码，对于信息比特位部分则要根据接收到之前的所有已经被译码的比特和当前计算得到的似然比的值进行判决。
极化码译码过程涉及f函数节点和g函数节点运算，在对数似然比域译码时，f函数节点的运算采用基于双曲正切规则的和积运算。和积算法涉及双曲正切函数及其反函数计算，运算复杂度较高。对于中长码字，和积算法难以保证实时性，一般采用查表法计算双曲正切函数，但需要先将双曲正切函数值存储到只读存储器(ROM)中，然后将输入转化为查表地址得到双曲正切函数的近似值，在一定精度下该方法需要消耗大量的硬件资源。
发明内容
本发明所要解决的技术问题是提供一种极化码译码方法，解决目前极化码在对数似然比域译码运算复杂度较高的问题。
本发明为解决上述技术问题所采用的技术方案是：一种极化码译码方法，其特征在于：极化码参数为
Figure BDA0001812314220000011
其中N表示极化码的码长，K表示信息位的长度，N-K表示冻结位的长度，A表示信息位位置的集合，令Ac为A的补集，则Ac表示冻结位位置的集合，
Figure BDA0001812314220000012
表示长度为N-K的冻结位的二进制向量，一般将冻结位的值设置为0，即令
Figure BDA0001812314220000013
设yi(0≤i＜N)为经过二进制相移键控调制以及AWGN信道后接收端接收到的信号，令矩阵LL表示(n+1)×N的对数似然比矩阵，矩阵B表示(n+1)×N的位矩阵，其中n＝log2 N，并将LL,B设置为非数值(NaN)，矩阵LL的更新函数为Update LLl,i，(0≤l＜log2 N)，其中l对应于译码阶段的列索引，i对应于行索引，LLl,i表示索引为(l,i)的节点的对数似然比的值；矩阵B的更新函数为Update Bl,i，(0≤l＜log2 N)，Bl,i表示索引为(l,i)的节点的位的更新值，部分和
Figure BDA0001812314220000021
则极化码译码的具体步骤为：
步骤一、初始化：接收信号yi(0≤i＜N)的对数似然比LLn,i的值为：
Figure BDA0001812314220000022
EbN0＝10^(SNR/10)，其中n＝log2 N，SNR为AWGN信道的信噪比；设j(1≤j≤N)表示第j次译码循环，并令j＝1；
步骤二、码位倒序计算：i＝bitreversal(j-1)，i表示(j-1)的比特反序，即将十进制(j-1)转化为二进制后，将得到的二进制数进行倒序，然后转化为十进制得到i；
步骤三、更新函数Update LL0,i得到LL0,i和矩阵LL：函数Update LL0,i是函数Update LLl,i(0≤l＜log2 N)在l＝0时的情况，更新函数Update LL0,i即更新函数UpdateLLl,i(0≤l＜log2 N)。
步骤四、判断i是否属于冻结位：若i属于冻结位，则令B0,i＝0；若i不属于冻结位且LL0,i≥0，则B0,i＝0，否则，B0,i＝1。
步骤五、更新函数Update B0,i得到B0,i和矩阵B：函数Update B0,i是函数UpdateBl,i(0≤l＜log2 N)在l＝0时的情况，更新函数Update B0,i即更新函数Update Bl,i(0≤l＜log2 N)。
步骤六、判断j＜N是否成立，若成立，令j＝j+1，重复步骤二～步骤六，否则，执行步骤七。
步骤七：得到译码结果
Figure BDA0001812314220000023
从矩阵B的第一行依次提取出信息位的值作为译码结果
Figure BDA0001812314220000031
本发明所述步骤三中函数Update LLl,i(0≤l＜log2 N)的更新方法为：
步骤3.1、输入索引l,i；
步骤3.2、令s＝2n-l，
Figure BDA0001812314220000032
其中n＝log2 N，
Figure BDA0001812314220000033
表示模2加运算；
步骤3.3、若m＜s/2，且LLl+1,i＝NaN，则更新函数Update LLl+1,i得到LLl+1,i，同时若LLl+1,i+s/2＝NaN，则更新函数Update LLl+1,i+s/2得到LLl+1,i+s/2，然后根据公式(1)计算f函数节点对数似然比的值LLl,i：
LLl,i＝2h2(h1(x1)h1(x2)) (1)；
其中变量
Figure BDA0001812314220000034
变量
Figure BDA0001812314220000035
x1,x2∈(-∞,∞)，函数h1(x)的表达式为公式(2)：
Figure BDA0001812314220000036
函数h2(x′)的表达式为公式(3)：
Figure BDA0001812314220000037
若m＜s/2不成立，则执行步骤(3.4)；
步骤3.4、根据公式(4)计算g函数节点对数似然比的值，此时需要部分和
Figure BDA0001812314220000041
部分和
Figure BDA0001812314220000042
若Bl,i-s/2＝0，则根据公式(5)计算g函数节点对数似然比的值LLl,i，若Bl,i-s/2＝0不成立，则根据公式(6)计算g函数节点对数似然比的值LLl,i，计算f函数或g函数节点对数似然比的值LLl,i后输出更新后的矩阵LL；
Figure BDA0001812314220000043
LLl,i＝LLl+1,i+LLl+1,i-s/2 (5)
LLl,i＝LLl+1,i-LLl+1,i-s/2 (6)。
本发明所述步骤五中函数Update Bl,i(0≤l＜log2 N)的更新方法为：
步骤5.1、输入索引l,i；
步骤5.2、令s＝2n-l，
Figure BDA0001812314220000044
其中n＝log2 N，
Figure BDA0001812314220000045
表示模2加运算；
步骤5.3、若m＜s/2不成立，则
Figure BDA0001812314220000046
Bl+1,i＝Bl,i，更新函数Update Bl+1,i和Update Bl+1,i-s/2，并输出更新后的矩阵B；若m＜s/2成立，则执行步骤六。
本发明的有益效果是：本发明针对f函数节点的计算，用简单的折线逼近算法代替双曲正切函数及其反函数运算，仅需要加法运算与乘法运算，大大降低了运算复杂度，与查表法相比大大节约了内存资源，降低了硬件成本，解决了极化码在对数似然比域译码运算复杂度较高的问题。
附图说明
图1为本发明极化码SC译码方法主要流程示意图；
图2为本发明的f函数的实现结构图；
图3为本发明的折线逼近函数h1(x)的实现结构图；
图4为本发明的折线逼近函数h2(x′)的实现结构图；
图5为本发明g函数的实现结构图；
图6为本发明函数Update Bl,i的更新流程图；
图7为本发明的折线逼近函数h1(x)与双曲正切函数函数的对比图；
图8为本发明的折线逼近函数h2(x′)与反双曲正切函数函数的对比图；
图9为本发明实施例N＝8时极化码的编码图；
图10为本发明实施例N＝8时极化码的译码图。
具体实施方式
下面结合说明书附图对本发明的具体实施方式(实施例)进行描述，使本领域的技术人员能够更好地理解本发明。
如图1所示，一种极化码译码方法，包括以下步骤：
步骤一、初始化：接收信号yi(0≤i＜N)的对数似然比的值LLn,i为：
Figure BDA0001812314220000051
EbN0＝10^(SNR/10)，其中n＝log2 N，SNR为AWGN信道的信噪比；设j(1≤j≤N)表示第j次译码循环，并令j＝1；
步骤二、码位倒序计算：i＝bitreversal(j-1)，i表示(j-1)的比特反序，即将十进制(j-1)转化为二进制后，将得到的二进制数进行倒序，然后转化为十进制得到i；
步骤三、更新函数Update LL0,i得到LL0,i和矩阵LL：函数Update LL0,i是函数Update LLl,i(0≤l＜log2 N)在l＝0时的情况，更新函数Update LL0,i即更新函数UpdateLLl,i(0≤l＜log2 N)，其步骤为：
步骤(a)、输入索引l,i。
步骤(b)、令s＝2n-l，
Figure BDA0001812314220000052
其中n＝log2 N，
Figure BDA0001812314220000053
表示模2加运算。
步骤(c)、若m＜s/2，且LLl+1,i＝NaN，则更新函数Update LLl+1,i得到LLl+1,i，同时若LLl+1,i+s/2＝NaN，则更新函数Update LLl+1,i+s/2得到LLl+1,i+s/2，然后根据公式(1)计算f函数节点对数似然比的值LLl,i。f函数的实现结构用如图2所示的结构实现，包括h1(x1)实现单元、h1(x2)实现单元和h2(x′)实现单元，由于x1,x2∈(-∞,∞)，所以h1(x1)和h1(x2)实现单元都用h1(x)的实现结构。
函数h1(x)的实现结构如图3所示，包括区间检测电路、参数存储器、旁路缓冲器、加法器、乘法器。在区间检测电路的RAM中存入区间转折点的集合X1，X1＝{-∞,-7.0,-3.0,-1.6,-0.8,0.8,1.6,3.0,7.0,∞}，在参数存储器中存入各分段的系数K1(k)和b1(k)，K1＝{0,0.0012,0.0524,0.322,0.83,0.322,0.0524,0.0012,0}，b1＝{-1.0000,-0.9914,-0.8378,-0.4064,0,0.4064,0.8378,0.9914,1.0000}，根据X1中元素值将(-∞,∞)分为9个区间，(X1(k),X1(k+1)],k＝1,…8为第k区间，(X(9),X(10))为第9区间。
函数h1(x)的具体实现步骤为：
(1)将输入值x同时送入旁路缓冲区和区间检测电路，在区间检测电路中与X1中的值逐个比较，若x∈(X1(k),X1(k+1)],k＝1,…8，则x属于第k个区间。若输入值x∈(X(9),X(10))，则x属于第9个区间。
(2)根据k值在参数存储单元中取出K1(k)和b1(k)。
(3)将K1(k)与接收缓冲器中的x相乘并将所得结果与b1(k)相加，得到输出值h1(x)。
函数h1(x)的计算平均需要5次比较操作，1次加法操作和1次乘法操作。
函数h2(x′)的实现结构如图4所示，包括区间检测电路、参数存储器、旁路缓冲器、加法器、乘法器。在区间检测电路的RAM中存入区间转折点的集合X2,X2＝{-∞,-0.999998,-0.9951,-0.9217,-0.6640,0.6640,0.9217,0.9951,0.999998,∞},在参数存储器中存入各分段的系数K2(k′)和b2(k′)，K2＝{0,833.3333,19.0840,3.1056,1.2048,3.1056,19.0840,833.3333,0}，b2＝{-7.0000,826.1667,15.9885,1.2621,0,-1.2621,-15.9885,-826.1667,7.0000}根据X2中元素值将(-∞,∞)分为9个区间，(X2(k′),X2(k′+1)],k′＝1,…8为第k′区间，(X(9),X(10))为第9区间。
函数h2(x′)的具体实现步骤为：
(1)将输入值x′同时送入旁路缓冲区和区间检测电路，在区间检测电路中与X2中的值逐个比较，若x′∈(X2(k′),X2(k′+1)],k′＝1,…8，则x′属于第k′个区间。若输入值x′∈(X(9),X(10))，则x′属于第9个区间。
(2)根据k′值在参数存储单元中取出K2(k′)和b2(k′)。
(3)将K2(k′)与接收缓冲器中的x′相乘并将所得结果与b2(k′)相加，得到输出值h2(x′)，即h2(x′)＝K2(k′)*X2(k′)+b2(k′)。
函数h2(x′)的计算平均需要5次比较操作，1次加法操作和1次乘法操作。
由以上函数h1(x)和函数h2(x′)的计算可以得到f函数节点对数似然比的值：LLl,i＝2h2(h1(x1)h1(x2))。
若m＜s/2不成立，则执行步骤(d)。
步骤(d)、根据公式(4)计算g函数节点对数似然比的值，此时需要部分和
Figure BDA0001812314220000071
部分和
Figure BDA0001812314220000072
若Bl,i-s/2＝0，则根据公式(5)计算g函数节点对数似然比的值LLl,i，若Bl,i-s/2＝0不成立，则根据公式(6)计算g函数节点对数似然比的值LLl,i。g函数的实现结构如图5所示。计算f函数或g函数节点对数似然比的值LLl,i后输出更新后的矩阵LL。
步骤四、判断i是否属于冻结位，若i属于冻结位，则令B0,i＝0；若i不属于冻结位且LL0,i≥0，则B0,i＝0，否则，B0,i＝1。
步骤五、更新函数Update B0,i得到B0,i和矩阵B：函数Update B0,i是函数UpdateBl,i(0≤l＜log2 N)在l＝0时的情况，更新函数Update B0,i即更新函数Update Bl,i(0≤l＜log2 N)，函数Update Bl,i(0≤l＜log2 N)的更新流程图如图6所示，其步骤为：
步骤(a)、输入索引l,i；
步骤(b)、令s＝2n-l，
Figure BDA0001812314220000073
其中n＝log2 N，
Figure BDA0001812314220000074
表示模2加运算；
步骤(c)、若m＜s/2不成立，则
Figure BDA0001812314220000081
Bl+1,i＝Bl,i，更新函数Update Bl+1,i和Update Bl+1,i-s/2，并输出更新后的矩阵B；若m＜s/2成立，则执行步骤六。
步骤六、判断j＜N-1是否成立，若成立，令j＝j+1，重复步骤二～步骤六，否则，执行步骤七。
步骤七：得到译码结果
Figure BDA0001812314220000085
从矩阵B的第一行依次提取出信息位的值作为译码结果
Figure BDA0001812314220000086
本发明简化后的折线逼近函数h1(x)与双曲正切函数tanh(x)的对比图如图7所示，从图7可以看出h1(x)与tanh(x)函数非常逼近，所以可以用折线逼近函数代替双曲正切函数进行运算；简化后的折线逼近函数h2(x′)与反双曲正切函数tanh-1(x′)的对比图如图8所示，从图8可以看出h2(x′)与tanh-1(x′)函数非常逼近，所以可以用折线逼近函数代替反双曲正切函数进行运算。
实施例
如图9所示，若已知极化码的码长为N＝8，信息位的长度为K＝4，冻结位的长度为N-K＝4，冻结比特序号的集合Ac＝{0,1,2,4}，信息比特序号的集合A＝{3,5,6,7}，冻结比特向量
Figure BDA0001812314220000082
信息比特向量uA＝(i0,i1,i2,i3)，则混合后得到
Figure BDA0001812314220000083
编码后的符号为(x0,...,x7)。将编码后的符号经过BPSK调制后送到AWGN信道下，则接收机接收到信息(y0,...,y7)，然后将(y0,...,y7)送到数据缓冲区。N＝8，K＝4时极化码的译码可由图10来表示，简化译码的具体步骤为：
步骤一、初始化：逐位计算接收信号yi(0≤i＜N)的对数似然比的值LLn,i：
Figure BDA0001812314220000084
EbN0＝10^(SNR/10)，其中n＝log2 N，SNR为AWGN信道的信噪比；设j(1≤j≤N)表示第j次译码循环，并令j＝1；
步骤二、码位倒序计算：i＝bitreversal(j-1)。
步骤三、更新函数Update LL0,i得到LL0,i和矩阵LL：更新函数Update LL0,i即更新函数Update LLl,i(0≤l＜log2 N)，其步骤为：
步骤(a)、输入索引l,i；
步骤(b)、令s＝2n-l，
Figure BDA0001812314220000091
其中n＝log2 N，
Figure BDA0001812314220000092
表示模2加运算；
步骤(c)、若m＜s/2，且LLl+1,i＝NaN，则更新函数Update LLl+1,i得到LLl+1,i，同时若LLl+1,i+s/2＝NaN，则更新函数Update LLl+1,i+s/2得到LLl+1,i+s/2，然后根据公式(1)计算f函数节点对数似然比的值LLl,i。若m＜s/2不成立，则执行步骤(d)；
步骤(d)、根据公式(4)计算g函数节点对数似然比的值，此时需要部分和
Figure BDA0001812314220000093
部分和
Figure BDA0001812314220000094
若Bl,i-s/2＝0，则根据公式(5)计算g函数节点对数似然比的值LLl,i，若Bl,i-s/2＝0不成立，则根据公式(6)计算g函数节点对数似然比的值LLl,i。计算f函数或g函数节点对数似然比的值LLl,i后输出更新后的矩阵LL。
步骤四、判断i是否属于冻结位，若i属于冻结位，则令B0,i＝0；若i不属于冻结位且LL0,i≥0，则B0,i＝0，否则，B0,i＝1。
步骤五、更新函数Update B0,i得到B0,i和矩阵B：更新函数Update B0,i即更新函数Update Bl,i(0≤l＜log2 N)，其步骤为：
步骤(a)、输入索引l,i；
步骤(b)、令s＝2n-l，
Figure BDA0001812314220000095
其中n＝log2 N，
Figure BDA0001812314220000096
表示模2加运算；
步骤(c)、若m＜s/2不成立，则
Figure BDA0001812314220000097
Bl+1,i＝Bl,i，更新函数Update Bl+1,i和Update Bl+1,i-s/2，并输出更新后的矩阵B；若m＜s/2成立，则执行步骤六。
步骤六、判断j＜N-1是否成立，若成立，令j＝j+1，重复步骤二～步骤六，否则，执行步骤七。
步骤七：得到译码结果
Figure BDA0001812314220000098
从矩阵B的第一行
Figure BDA0001812314220000099
依次提取出信息位的值作为译码结果
Figure BDA0001812314220000101
Figure BDA0001812314220000102

Claims (1)

1.一种极化码译码方法，其特征在于：极化码参数为(N,K,A,
Figure FDA0002762348520000015
)，其中N表示极化码的码长，K表示信息位的长度，N-K表示冻结位的长度，A表示信息位位置的集合，令Ac为A的补集，则Ac表示冻结位位置的集合，
Figure FDA0002762348520000016
表示长度为N-K的冻结位的二进制向量，一般将冻结位的值设置为0，即令
Figure FDA0002762348520000017
设yi(0≤i＜N)为经过二进制相移键控调制以及AWGN信道后接收端接收到的信号，令矩阵LL表示(n+1)×N的对数似然比矩阵，矩阵B表示(n+1)×N的位矩阵，其中n＝log2 N，并将LL,B设置为非数值(NaN)，矩阵LL的更新函数为Update LLl,i，(0≤l＜log2 N)，其中l对应于译码阶段的列索引，i对应于行索引，LLl,i表示索引为(l,i)的节点的对数似然比的值；矩阵B的更新函数为Update Bl,i，(0≤l＜log2 N)，Bl,i表示索引为(l,i)的节点的位的更新值，部分和
Figure FDA0002762348520000011
则极化码译码的具体步骤为：
步骤一、初始化：接收信号yi(0≤i＜N)的对数似然比LLn,i的值为：
Figure FDA0002762348520000012
EbN0＝10^(SNR/10)，其中n＝log2 N，SNR为AWGN信道的信噪比；设j(1≤j≤N)表示第j次译码循环，并令j＝1；
步骤二、码位倒序计算：i＝bitreversal(j-1)，i表示(j-1)的比特反序，即将十进制(j-1)转化为二进制后，将得到的二进制数进行倒序，然后转化为十进制得到i；
步骤三、更新函数Update LL0,i得到LL0,i和矩阵LL：函数Update LL0,i是函数UpdateLLl,i(0≤l＜log2 N)在l＝0时的情况，更新函数Update LL0,i即更新函数Update LLl,i(0≤l＜log2 N)；
所述步骤三中函数Update LLl,i(0≤l＜log2 N)的更新方法为：
步骤3.1、输入索引l,i；
步骤3.2、令s＝2n-l，
Figure FDA0002762348520000013
其中n＝log2 N，
Figure FDA0002762348520000014
表示模2加运算；
步骤3.3、若m＜s/2，且LLl+1,i＝NaN，则更新函数Update LLl+1,i得到LLl+1,i，同时若LLl+1,i+s/2＝NaN，则更新函数Update LLl+1,i+s/2得到LLl+1,i+s/2，然后根据公式(1)计算f函数节点对数似然比的值LLl,i：
LLl,i＝2h2(h1(x1)h1(x2)) (1)；
其中变量
Figure FDA0002762348520000021
变量
Figure FDA0002762348520000022
x1,x2∈(-∞,∞)，函数h1(x)的表达式为公式(2)：
Figure FDA0002762348520000023
函数h2(x′)的表达式为公式(3)：
Figure FDA0002762348520000024
若m＜s/2不成立，则执行步骤(3.4)；
步骤3.4、根据公式(4)计算g函数节点对数似然比的值，此时需要部分和
Figure FDA0002762348520000025
部分和
Figure FDA0002762348520000026
若Bl,i-s/2＝0，则根据公式(5)计算g函数节点对数似然比的值LLl,i，若Bl,i-s/2＝0不成立，则根据公式(6)计算g函数节点对数似然比的值LLl,i，计算f函数或g函数节点对数似然比的值LLl,i后输出更新后的矩阵LL；
Figure FDA0002762348520000031
LLl,i＝LLl+1,i+LLl+1,i-s/2 (5)
LLl,i＝LLl+1,i-LLl+1,i-s/2 (6)；
步骤四、判断i是否属于冻结位：若i属于冻结位，则令B0,i＝0；若i不属于冻结位且LL0,i≥0，则B0,i＝0，否则，B0,i＝1；
步骤五、更新函数Update B0,i得到B0,i和矩阵B：函数Update B0,i是函数Update Bl,i(0≤l＜log2 N)在l＝0时的情况，更新函数Update B0,i即更新函数Update Bl,i(0≤l＜log2N)；
所述步骤五中函数Update Bl,i(0≤l＜log2 N)的更新方法为：
步骤5.1、输入索引l,i；
步骤5.2、令s＝2n-l，
Figure FDA0002762348520000032
其中n＝log2 N，
Figure FDA0002762348520000033
表示模2加运算；
步骤5.3、若m＜s/2不成立，则
Figure FDA0002762348520000034
Bl+1,i＝Bl,i，更新函数UpdateBl+1,i和Update Bl+1,i-s/2，并输出更新后的矩阵B；若m＜s/2成立，则执行步骤六；
步骤六、判断j＜N是否成立，若成立，令j＝j+1，重复步骤二～步骤六，否则，执行步骤七；
步骤七：得到译码结果
Figure FDA0002762348520000035
从矩阵B的第一行依次提取出信息位的值作为译码结果
Figure FDA0002762348520000036
CN201811125637.8A 2018-09-26 2018-09-26 一种极化码译码方法 Active CN109347600B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811125637.8A CN109347600B (zh) 2018-09-26 2018-09-26 一种极化码译码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811125637.8A CN109347600B (zh) 2018-09-26 2018-09-26 一种极化码译码方法

Publications (2)

Publication Number Publication Date
CN109347600A CN109347600A (zh) 2019-02-15
CN109347600B true CN109347600B (zh) 2021-01-05

Family

ID=65306719

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811125637.8A Active CN109347600B (zh) 2018-09-26 2018-09-26 一种极化码译码方法

Country Status (1)

Country Link
CN (1) CN109347600B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112152639A (zh) 2019-06-27 2020-12-29 深圳市中兴微电子技术有限公司 一种极化码的译码方法、装置、存储介质和终端

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105099622A (zh) * 2015-07-01 2015-11-25 北京邮电大学 极化编码调制中信道可靠度的确定方法及装置
CN105187073A (zh) * 2015-10-13 2015-12-23 东南大学 一种极化码的bp译码方法及装置
CN105720992A (zh) * 2016-01-22 2016-06-29 哈尔滨工业大学深圳研究生院 一种极化码的简化译码方法
CN105939160A (zh) * 2016-01-04 2016-09-14 西南大学 LTE-Advanced标准中低存储容量Turbo码译码器及设计方法
CN106656212A (zh) * 2016-12-05 2017-05-10 东南大学 基于极化码的自适应连续消除译码方法及架构
CN107769894A (zh) * 2017-11-16 2018-03-06 东南大学 极化码编码的mimo系统的联合检测译码方法
CN107888202A (zh) * 2017-11-29 2018-04-06 哈尔滨工业大学 一种非递归的sc译码似然比确定方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180198560A1 (en) * 2017-01-09 2018-07-12 Qualcomm Incorporated Generalized polar code based on polarization of linear block codes and convolutional codes

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105099622A (zh) * 2015-07-01 2015-11-25 北京邮电大学 极化编码调制中信道可靠度的确定方法及装置
CN105187073A (zh) * 2015-10-13 2015-12-23 东南大学 一种极化码的bp译码方法及装置
CN105939160A (zh) * 2016-01-04 2016-09-14 西南大学 LTE-Advanced标准中低存储容量Turbo码译码器及设计方法
CN105720992A (zh) * 2016-01-22 2016-06-29 哈尔滨工业大学深圳研究生院 一种极化码的简化译码方法
CN106656212A (zh) * 2016-12-05 2017-05-10 东南大学 基于极化码的自适应连续消除译码方法及架构
CN107769894A (zh) * 2017-11-16 2018-03-06 东南大学 极化码编码的mimo系统的联合检测译码方法
CN107888202A (zh) * 2017-11-29 2018-04-06 哈尔滨工业大学 一种非递归的sc译码似然比确定方法及装置

Also Published As

Publication number Publication date
CN109347600A (zh) 2019-02-15

Similar Documents

Publication Publication Date Title
CN107241106B (zh) 基于深度学习的极化码译码算法
Xu et al. Improved polar decoder based on deep learning
CN104079382B (zh) 一种基于概率计算的极化码译码器和极化码译码方法
US5471500A (en) Soft symbol decoding
CN108847848B (zh) 一种基于信息后处理的极化码的bp译码算法
EP0660534B1 (en) Error correction systems with modified viterbi decoding
CN107612560B (zh) 基于部分信息比特似然比的极化码早期迭代停止方法
CN110071779B (zh) 一种低复杂度的极化码多级编码调制方法
KR20070080246A (ko) 하이브리드 서바이버 메모리 구조체
CN110278001B (zh) 基于深度学习的极化码分区译码方法
JP2009232476A (ja) 情報処理方法および情報処理装置
KR101583139B1 (ko) 높은 처리량과 낮은 복잡성을 갖는 연속 제거 극 부호 복호 장치 및 그 방법
CN109347600B (zh) 一种极化码译码方法
Huang et al. Latency reduced method for modified successive cancellation decoding of polar codes
CN107438985B (zh) 处理快于奈奎斯特的信令接收信号
CN111835364B (zh) 一种极化码的低复杂度神经bp译码方法
JP3728171B2 (ja) 信頼性情報計算方法
CN109639290B (zh) 一种半随机分组叠加编码及译码方法
JP5859147B2 (ja) ビット尤度演算装置およびビット尤度演算方法
CN110855298B (zh) 基于子信道冻结条件的低迭代次数极化码bp译码方法
Zhang et al. Adaptive iterative decoding for expediting the convergence of unary error correction codes
CN114257342A (zh) 基于非二进制极化码的两用户多址接入系统编码及译码方法
Xiong et al. Deep learning method of polar codes under colored noise
CN115378443A (zh) 一种低精度的深度神经网络极化码sc译码算法
CN109167651B (zh) 一种传输高斯信源的方法

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