本发明公开了一种基于简化双二元Turbo码度量值的译码方法，主要解决现有译码方法复杂度大效率低的问题。本发明的步骤为：接收数据；初始化；设置译码迭代次数；计算前半轮迭代的后向度量值；计算前半轮迭代的后验信息和前半轮迭代的外信息；获得后半轮迭代的先验信息；计算后半轮迭代的后向度量值；计算后半轮迭代的后验信息和后半轮迭代的外信息；获得前半轮迭代的先验信息；译码终止判决；判决对数似然比信息。本发明对不同分支度量表达式之间的相同部分只进行了一次计算，对前向度量中间值和后向度量中间值先比较再求解，从而减少译码计算量，降低了译码复杂度，提高了译码效率。
1.一种基于简化双二元Turbo码度量值的译码方法，包括如下步骤：(1)接收数据：
(1a)依据双二元Turbo码交织器的交织深度，按照时序，接收两路系统位软信息，每路系统位软信息的长度与交织深度相等；按照时序，接收两路第一校验位软信息和第二校验位软信息，每路验位软信息的长度与交织深度相等；
(1b)利用双二元Turbo码交织器，对两路的系统位软信息进行交织，得到两路交织后的系统位软信息；
(2b)将前向度量和后向度量初始化为 其中，log(·)表示以自然常数e为底的对数操作；L表示双二元Turbo码分量编码器中寄存器的总数；
设置译码最大迭代次数为M，M∈{1,2，...,10}，∈表示属于操作；将当前迭代次数初始化为1；
(4a)合并接收的两路系统位软信息和两路第一校验位软信息，得到两个合并后的系统位软信息和两个合并后的校验位软信息；
(4b)由合并后的系统位软信息和合并后的校验位软信息，得到8个信道观测值；
(4c)由信道观测值和前半轮迭代的先验信息，计算16个分支度量值；
(4d)将起始状态相同且结束状态不同的所有后向度量中间值，以四个中间值为一组，划分成多个后向度量中间值小组；
(4e)利用下式，在每个后向度量中间值小组中，分别计算后向度量差值和分支度量差值：X＝B1-C1
其中，X表示后向度量差值；B1表示带有对应于码字(0,1)的先验信息的后向度量中间值中的后向度量；C1表示带有对应于码字(1,0)的先验信息的后向度量中间值中的后向度量；P表示分支度量差值；C2表示带有对应于码字(1,0)的先验信息的后向度量中间值中的分支度量；B2表示带有对应于码字(0,1)的先验信息的后向度量中间值中的分支度量；
(4f)判断后向度量差值X是否大于分支度量差值P，若是，则将带有对应于码字(0,1)的先验信息的后向度量中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,
0)的先验信息的后向度量中间值作为第一个待选后向度量中间值；
(4g)使用与步骤(4e)和步骤(4f)中相同的方法，由带有对应于码字(0,0)的先验信息的后向度量中间值和带有对应于码字(1,1)的先验信息的后向度量中间值，获得第二个待选后向度量中间值；
(4h)取每个后向度量中间值小组中两个待选后向度量中间值中最大的一个，作为该后向度量中间值小组选取的后向度量；
(5)计算前半轮迭代的后验信息和前半轮迭代的外信息：
(5a)将结束状态相同且起始状态不同的所有前向度量中间值，以四个中间值为一组，划分成多个前向度量中间值小组；
(5b)利用下式，在每个前向度量中间值小组中，分别计算前向度量差值和分支度量差值：Y＝A1-D1
其中，Y表示前向度量差值；A1表示带有对应于码字(0,0)的先验信息的前向度量中间值中的前向度量；D1表示带有对应于码字(1,1)的先验信息的前向度量中间值中的前向度量；Q表示分支度量差值；D2表示带有对应于码字(1,1)的先验信息的前向度量中间值中的分支度量；A2表示带有对应于码字(0,0)的先验信息的前向度量中间值中的分支度量；
(5c)判断前向度量差值Y是否大于分支度量差值Q，若是，则将带有对应于码字(0,0)的先验信息的前向度量的中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,1)的先验信息的前向度量中间值作为第一个待选前向度量中间值；
(5d)使用与步骤(5b)和步骤(5c)中相同的方法，由带有对应于码字(0,1)的先验信息的前向度量中间值和带有对应于码字(1,0)的先验信息的前向度量中间值，获得第二个待选前向度量中间值；
(5e)取每个前向度量中间值小组中两个待选前向度量中间值中最大的一个，作为该前向度量中间值小组选取的前向度量；
(5f)利用后验信息公式，计算前半轮迭代的后验信息；
利用双二元Turbo码交织器对前半轮迭代的外信息进行交织，获得后半轮迭代的先验信息；
(7a)合并两路交织后的系统位软信息和两路第二校验位软信息，得到两个合并后的系统位软信息和两个合并后的校验位软信息；
(7b)由合并后的系统位软信息和合并后的校验位软信息，合并得到8个信道观测值；
(7c)由信道观测值和先验信息，计算得到16个分支度量值；
(7d)将起始状态相同且结束状态不同的所有后向度量中间值，以四个中间值为一组，划分成多个后向度量中间值小组；
(7e)利用下式，在每个后向度量中间值小组中，分别计算后向度量差值和分支度量差值：X＝B1-C1
其中，X表示后向度量差值；B1表示带有对应于码字(0,1)的先验信息的后向度量中间值中的后向度量；C1表示带有对应于码字(1,0)的先验信息的后向度量中间值中的后向度量；P表示分支度量差值；C2表示带有对应于码字(1,0)的先验信息的后向度量中间值中的分支度量；B2表示带有对应于码字(0,1)的先验信息的后向度量中间值中的分支度量；
(7f)判断后向度量差值X是否大于分支度量差值P，若是，则将带有对应于码字(0,1)的先验信息的后向度量中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,
0)的先验信息的后向度量中间值作为第一个待选后向度量中间值；
(7g)使用与步骤(7e)和步骤(7f)中相同的方法，由带有对应于码字(0,0)的先验信息的后向度量中间值和带有对应于码字(1,1)的先验信息的后向度量中间值，获得第二个待选后向度量中间值；
(7h)取每个后向度量中间值小组中两个待选后向度量中间值中最大的一个，作为该后向度量中间值小组选取的后向度量；
(8)计算后半轮迭代的后验信息和后半轮迭代的外信息：
(8a)将结束状态相同且起始状态不同的后半轮迭代的前向度量中间值，以四个中间值为一组，划分成多个前向度量中间值小组；
(8b)利用下式，在每个前向度量中间值小组中，分别计算前向度量差值和分支度量差值：Y＝A1-D1
其中，Y表示前向度量差值；A1表示带有对应于码字(0,0)的先验信息的前向度量中间值中的前向度量；D1表示带有对应于码字(1,1)的先验信息的前向度量中间值中的前向度量；Q表示分支度量差值；D2表示带有对应于码字(1,1)的先验信息的前向度量中间值中的分支度量；A2表示带有对应于码字(0,0)的先验信息的前向度量中间值中的分支度量；
(8c)判断前向度量差值Y是否大于分支度量差值Q，若是，则将带有对应于码字(0,0)的先验信息的前向度量的中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,1)的先验信息的前向度量中间值作为第一个待选前向度量中间值；
(8d)使用与步骤(8b)和步骤(8c)中相同的方法，由带有对应于码字(0,1)的先验信息的前向度量中间值和带有对应于码字(1,0)的先验信息的前向度量中间值，获得第二个待选前向度量中间值；
(8e)取每个前向度量中间值小组中两个待选前向度量中间值中最大的一个，作为该前向度量中间值小组选取的前向度量；
(8f)利用后验信息公式，计算后半轮迭代的后验信息；
利用双二元Turbo码解交织器，对后半轮迭代的外信息进行解交织，获得前半轮迭代的先验信息；
(10)判断当前迭代次数是否达到最大迭代次数，若是，则执行步骤(11)，否则，将当前迭代次数加1后执行步骤(4)；
(11a)使用双二元Turbo码解交织器，对后半轮迭代的后验信息进行解交织，得到对数似然比；
(11b)使用双二元Turbo码的符号判决规则，对对数似然比进行硬判决，译码结束。
2.根据权利要求1所述的基于简化双二元Turbo码度量值的译码方法，其特征在于，步骤(4a)中所述的合并两路交织后的系统位软信息和两路第一校验位软信息是按照下式得到的：S1＝Ark+Brk
其中，S1和S2分别表示两个合并后的系统位软信息；P1和P2分别表示两个合并后的校验位软信息；Ark表示k时刻所接收的两路系统位软信息中的一路系统位软信息；Brk表示k时刻所接收的两路系统位软信息中的另一路系统位软信息；Yrk表示k时刻所接收的两路第一校验位软信息中的一路第一校验位软信息；Wrk表示k时刻所接收的两路第一校验位软信息中的另一路第一校验位软信息，k依次取值1,2,…,N，N表示双二元Turbo码交织器的交织深度。
3.根据权利要求1所述的基于简化双二元Turbo码度量值的译码方法，其特征在于，步骤(4b)、步骤(7b)中所述的8个信道观测值是按照下式得到的：SP11＝S1+P1
其中，SP11、SP12、SP13、SP14、SP21、SP22、SP23、SP24分别表示由合并后的系统位软信息和合并后的校验位软信息，计算得到的8个信道观测值；S1和S2分别表示两个合并后的系统位软信息；P1和P2分别表示两个合并后的校验位软信息。
4.根据权利要求1所述的基于简化双二元Turbo码度量值的译码方法，其特征在于，步骤(4c)、步骤(7c)中所述的16个分支度量值是按照下式得到的：其中， 表示k时刻第j个分支度量值；k表示译码时刻，k∈{1,2,...,NL}，NL表示双二元Turbo码编码码长；j表示k时刻分支度量值序号，j∈{0,1,...,15}；SP11、SP12、SP13、SP14、SP21、SP22、SP23、SP24分别表示由合并后的系统位软信息和合并后的校验位软信息计算得到的8个信道观测值； 分别表示k时刻先验信息的3个数据。
5.根据权利要求1所述的基于简化双二元Turbo码度量值的译码方法，其特征在于，步骤(5f)、步骤(8f)中所述的后验信息公式如下：其中，Lk表示k时刻的后验信息，k的取值范围为[1,2,…,N]，N表示双二元Turbo码交织器的交织深度；max(·)表示取最大值操作；s′和s分别表示双二元Turbo码编码器状态状态转移图上的两个状态节点；Auk和Buk分别表示双二元Turbo码编码器k时刻两个输入的比特；
表示k-1时刻状态s′的前向度量； 表示k时刻第j个分支度量值，j表示分支度量序号，j∈{0,1,...,15}； 表示k时刻状态s的后向度量。
6.根据权利要求1所述的基于简化双二元Turbo码度量值的译码方法，其特征在于，步骤(5g)中所述的外信息公式如下：其中， 表示k时刻对应码字为(0,1)的前半轮迭代的外信息，k的取值范围为[1,
2,…,N]，N表示双二元Turbo码交织器的交织深度； 表示k时刻对应码字为(0,1)的前半轮迭代的后验信息；Brk表示k时刻所接收的两路系统位软信息中的一路系统位软信息；
表示k时刻对应码字为(0,1)的前半轮迭代的先验信息； 表示k时刻对应码字为(1,0)的前半轮迭代的外信息； 表示k时刻对应码字为(1,0)的前半轮迭代的后验信息；Ark表示k时刻所接收的两路系统位软信息中的另一路系统位软信息； 表示k时刻对应码字为(1,
0)的前半轮迭代的先验信息； 表示k时刻对应码字为(1,1)的前半轮迭代的外信息；
表示k时刻对应码字为(1,1)的前半轮迭代的后验信息； 表示k时刻对应码字为(1,1)的前半轮迭代的先验信息。
7.根据权利要求1所述的基于简化双二元Turbo码度量值的译码方法，其特征在于，步骤(7a)中所述的合并两路交织后的系统位软信息和两路第二校验位软信息是按照下式得到的：S1＝Ark+Brk
其中，S1和S2分别表示两个合并后的系统位软信息；P1和P2分别表示两个合并后的校验位软信息；Ark表示k时刻两路交织后的系统位软信息中的一路交织后的系统位软信息；Brk表示k时刻两路交织后的系统位软信息中的另一路交织后的系统位软信息；Yrk表示k时刻所接收的两路第二校验位软信息中的一路第二校验位软信息；Wrk表示k时刻所接收的两路第二校验位软信息中的另一路第二校验位软信息，k依次取值1,2,…,N，N表示双二元Turbo码交织器的交织深度。
8.根据权利要求1所述的基于简化双二元Turbo码度量值的译码方法，其特征在于，步骤(8g)中所述的外信息计算如下：其中， 表示k时刻对应码字为(0,1)的后半轮迭代的外信息，k的取值范围为[1,
2,…,N]，N表示双二元Turbo码交织器的交织深度； 表示k时刻对应码字为(0,1)的前半轮迭代的后验信息；Brk表示k时刻两路交织后的系统位软信息中的一路交织后的系统位软信息； 表示k时刻对应码字为(0,1)的后半轮迭代的先验信息； 表示k时刻对应码字为(1,0)的后半轮迭代的外信息； 表示k时刻对应码字为(1,0)的后半轮迭代的后验信息；Ark表示k时刻两路交织后的系统位软信息中的一路交织后的系统位软信息； 表示k时刻对应码字为(1,0)的后半轮迭代的先验信息； 表示k时刻对应码字为(1,1)的后半轮迭代的外信息； 表示k时刻对应码字为(1,1)的后半轮迭代的后验信息； 表示k时刻对应码字为(1,1)的后半轮迭代的先验信息。
基于简化双二元Turbo码度量值的译码方法
技术领域
[0001] 本发明属于无线通信技术领域，更进一步涉及信道编码技术领域中的基于简化双二元Turbo码度量值的译码方法。本发明可用于双二元Turbo码译码。
背景技术
[0002] Turbo码作为一种优秀的信道编码方案，性能接近香农限，自其提出以来就备受关注，但是由于Turbo码译码过程复杂，硬件实现难度较大，所以经过二十余年的研究，才被广泛的应用到各种通信标准之中。双二元Turbo码在经典Turbo码基础上发展而来，它具有更高的编码效率，更好的纠错性能，而且误码平层不明显，国际数字视频广播回传信道DVB-RCS(Digital Video Broadcasting-Return Channel via Satellite)和全球微波互联接入WIMAX(Worldwide Interoperability for Microwave Access)等标准采用了二元Turbo码作为信道编码方案。
[0003] 广州杰赛科技股份有限公司在其申请的专利文献“双二元卷积Turbo码译码方法和装置”(申请日：2009年11月02日，申请号：200910193580.X，公开号：CN101697492A)中公开了一种双二元卷积Turbo码译码方法。该方法在增强的Max-Log-MAP算法特有步骤处，并列地包括了常数Log-MAP算法相应的特有步骤，使得用户可以根据译码算法模式编号选择对应的译码算法。该方法存在的不足之处是，前向度量中间值和后向度量中间值只进行了计算，却没有对计算结果加以利用，造成了不必要的运算，从而增加了译码复杂度。
[0004] 李江鹏等人在其发表的论文“WIMAX系统双二进制Turbo码的译码算法及VLSI实现”([J].微电子学与计算机，2011,(03):43-47)中公开了一种双二元Turbo码译码方法。该译码方法在增加少量复杂度的基础上，对经典Max-Log-MAP译码算法进行了改进，提高了译码性能。但是，该方法仍然存在的不足之处是，在分支度量的计算过程中，存在重复的运算，由此增加了运算量，造成了资源的浪费。
发明内容
[0005] 本发明的目的在于克服上述现有技术的不足，提出一种基于简化双二元Turbo码度量值的译码方法。本发明可实现对双二元Turbo码译码，有效的降低了增强的Max-Log-MAP算法译码过程中度量值的计算复杂度。
[0006] 实现本发明目的的具体步骤包括如下：
[0007] (1)接收数据：
[0008] (1a)依据双二元Turbo码交织器的交织深度，按照时序，接收两路系统位软信息，每路系统位软信息的长度与交织深度相等；按照时序，接收两路第一校验位软信息和第二校验位软信息，每路验位软信息的长度与交织深度相等；
[0009] (1b)利用双二元Turbo码交织器，对两路的系统位软信息进行交织，得到两路交织后的系统位软信息；
[0010] (2)初始化：
[0011] (2a)将前半轮迭代的先验信息初始化为0；
[0012] (2b)将前向度量和后向度量初始化为 其中，log(·)表示以自然常数e为底的对数操作；L表示双二元Turbo码分量编码器中寄存器的总数；
[0013] (3)设置译码迭代次数：
[0014] 设置译码最大迭代次数为M，M∈{1,2，...,10}，∈表示属于操作；将当前迭代次数初始化为1；
[0015] (4)计算前半轮迭代的后向度量值：
[0016] (4a)合并接收的两路系统位软信息和两路第一校验位软信息，得到两个合并后的系统位软信息和两个合并后的校验位软信息；
[0017] (4b)由合并后的系统位软信息和合并后的校验位软信息，得到8个信道观测值；
[0018] (4c)由信道观测值和前半轮迭代的先验信息，计算16个分支度量值；
[0019] (4d)将起始状态相同且结束状态不同的所有后向度量中间值，以四个中间值为一组，划分成多个后向度量中间值小组；
[0020] (4e)利用下式，在每个后向度量中间值小组中，分别计算后向度量差值和分支度量差值：
[0021] X＝B1-C1
[0022] P＝C2-B2
[0023] 其中，X表示后向度量差值；B1表示带有对应于码字(0,1)的先验信息的后向度量中间值中的后向度量；C1表示带有对应于码字(1,0)的先验信息的后向度量中间值中的后向度量；P表示分支度量差值；C2表示带有对应于码字(1,0)的先验信息的后向度量中间值中的分支度量；B2表示带有对应于码字(0,1)的先验信息的后向度量中间值中的分支度量；
[0024] (4f)判断后向度量差值X是否大于分支度量差值P，若是，则将带有对应于码字(0,1)的先验信息的后向度量中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,0)的先验信息的后向度量中间值作为第一个待选后向度量中间值；
[0025] (4g)使用与步骤(4e)和步骤(4f)中相同的方法，由带有对应于码字(0,0)的先验信息的后向度量中间值和带有对应于码字(1,1)的先验信息的后向度量中间值，获得第二个待选后向度量中间值；
[0026] (4h)取每个后向度量中间值小组中两个待选后向度量中间值中最大的一个，作为该后向度量中间值小组选取的后向度量；
[0027] (5)计算前半轮迭代的后验信息和前半轮迭代的外信息：
[0028] (5a)将结束状态相同且起始状态不同的所有前向度量中间值，以四个中间值为一组，划分成多个前向度量中间值小组；
[0029] (5b)利用下式，在每个前向度量中间值小组中，分别计算前向度量差值和分支度量差值：
[0030] Y＝A1-D1
[0031] Q＝D2-A2
[0032] 其中，Y表示前向度量差值；A1表示带有对应于码字(0,0)的先验信息的前向度量中间值中的前向度量；D1表示带有对应于码字(1,1)的先验信息的前向度量中间值中的前向度量；Q表示分支度量差值；D2表示带有对应于码字(1,1)的先验信息的前向度量中间值中的分支度量；A2表示带有对应于码字(0,0)的先验信息的前向度量中间值中的分支度量；
[0033] (5c)判断前向度量差值Y是否大于分支度量差值Q，若是，则将带有对应于码字(0,0)的先验信息的前向度量的中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,1)的先验信息的前向度量中间值作为第一个待选前向度量中间值；
[0034] (5d)使用与步骤(5b)和步骤(5c)中相同的方法，由带有对应于码字(0,1)的先验信息的前向度量中间值和带有对应于码字(1,0)的先验信息的前向度量中间值，获得第二个待选前向度量中间值；
[0035] (5e)取每个前向度量中间值小组中两个待选前向度量中间值中最大的一个，作为该前向度量中间值小组选取的前向度量；
[0036] (5f)利用后验信息公式，计算前半轮迭代的后验信息；
[0037] (5g)利用外信息公式，计算前半轮迭代的外信息；
[0038] (6)获得后半轮迭代的先验信息：
[0039] 利用双二元Turbo码交织器对前半轮迭代的外信息进行交织，获得后半轮迭代的先验信息；
[0040] (7)计算后半轮迭代的后向度量值：
[0041] (7a)合并两路交织后的系统位软信息和两路第二校验位软信息，得到两个合并后的系统位软信息和两个合并后的校验位软信息；
[0042] (7b)由合并后的系统位软信息和合并后的校验位软信息，合并得到8个信道观测值；
[0043] (7c)由信道观测值和先验信息，计算得到16个分支度量值；
[0044] (7d)将起始状态相同且结束状态不同的所有后向度量中间值，以四个中间值为一组，划分成多个后向度量中间值小组；
[0045] (7e)利用下式，在每个后向度量中间值小组中，分别计算后向度量差值和分支度量差值：
[0046] X＝B1-C1
[0047] P＝C2-B2
[0048] 其中，X表示后向度量差值；B1表示带有对应于码字(0,1)的先验信息的后向度量中间值中的后向度量；C1表示带有对应于码字(1,0)的先验信息的后向度量中间值中的后向度量；P表示分支度量差值；C2表示带有对应于码字(1,0)的先验信息的后向度量中间值中的分支度量；B2表示带有对应于码字(0,1)的先验信息的后向度量中间值中的分支度量；
[0049] (7f)判断后向度量差值X是否大于分支度量差值P，若是，则将带有对应于码字(0,1)的先验信息的后向度量中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,0)的先验信息的后向度量中间值作为第一个待选后向度量中间值；
[0050] (7g)使用与步骤(7e)和步骤(7f)中相同的方法，由带有对应于码字(0,0)的先验信息的后向度量中间值和带有对应于码字(1,1)的先验信息的后向度量中间值，获得第二个待选后向度量中间值；
[0051] (7h)取每个后向度量中间值小组中两个待选后向度量中间值中最大的一个，作为该后向度量中间值小组选取的后向度量；
[0052] (8)计算后半轮迭代的后验信息和后半轮迭代的外信息：
[0053] (8a)将结束状态相同且起始状态不同的后半轮迭代的前向度量中间值，以四个中间值为一组，划分成多个前向度量中间值小组；
[0054] (8b)利用下式，在每个前向度量中间值小组中，分别计算前向度量差值和分支度量差值：
[0055] Y＝A1-D1
[0056] Q＝D2-A2
[0057] 其中，Y表示前向度量差值；A1表示带有对应于码字(0,0)的先验信息的前向度量中间值中的前向度量；D1表示带有对应于码字(1,1)的先验信息的前向度量中间值中的前向度量；Q表示分支度量差值；D2表示带有对应于码字(1,1)的先验信息的前向度量中间值中的分支度量；A2表示带有对应于码字(0,0)的先验信息的前向度量中间值中的分支度量；
[0058] (8c)判断前向度量差值Y是否大于分支度量差值Q，若是，则将带有对应于码字(0,0)的先验信息的前向度量的中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,1)的先验信息的前向度量中间值作为第一个待选前向度量中间值；
[0059] (8d)使用与步骤(8b)和步骤(8c)中相同的方法，由带有对应于码字(0,1)的先验信息的前向度量中间值和带有对应于码字(1,0)的先验信息的前向度量中间值，获得第二个待选前向度量中间值；
[0060] (8e)取每个前向度量中间值小组中两个待选前向度量中间值中最大的一个，作为该前向度量中间值小组选取的前向度量；
[0061] (8f)利用后验信息公式，计算后半轮迭代的后验信息；
[0062] (8g)利用外信息公式，计算后半轮迭代的外信息；
[0063] (9)获得前半轮迭代的先验信息：
[0064] 利用双二元Turbo码解交织器，对后半轮迭代的外信息进行解交织，获得前半轮迭代的先验信息；
[0065] (10)判断当前迭代次数是否达到最大迭代次数，若是，则执行步骤(11)，否则，将当前迭代次数加1后执行步骤(4)；
[0066] (11)判决对数似然比信息：
[0067] (11a)使用双二元Turbo码解交织器，对后半轮迭代的后验信息进行解交织，得到对数似然比；
[0068] (11b)使用双二元Turbo码的符号判决规则，对对数似然比进行硬判决，译码结束。
[0069] 本发明与现有技术相比具有以下优点：
[0070] 第一，由于本发明在双二元Turbo译码的分支度量计算过程中，对不同分支度量表达式之间的相同部分只进行了一次计算，从而克服了现有技术在双二元Turbo译码的分支度量的计算过程中存在重复运算的缺点，使得本发明具有减少双二元Turbo译码计算量，提高译码效率的优点。
[0071] 第二，由于本发明使用的后向度量计算方法和前向度量计算方法，利用了译码过程中不同前向度量中间值和后向度量中间值的共性，对在双二元Turbo码的前向度量中间值和后向度量中间值先比较再求解，从而克服了现有技术在译码过程中获得后向度量和前向度量时计算量过大的缺点，使得本发明具有降低译码复杂度的优点。
附图说明
[0072] 图1是本发明的流程图；
[0073] 图2是本发明的仿真图。
具体实施方式
[0074] 下面结合附图对本发明做进一步的描述。
[0075] 参照图1，本发明实现的具体步骤如下：
[0076] 步骤1，接收数据。
[0077] 依据双二元Turbo码交织器的交织深度，按照时序，接收两路系统位软信息，每路系统位软信息的长度与交织深度相等；按照时序，接收两路第一校验位软信息和第二校验位软信息，每路验位软信息的长度与交织深度相等。
[0078] 利用双二元Turbo码交织器，对两路的系统位软信息进行交织，得到两路交织后的系统位软信息。
[0079] 步骤2，初始化。
[0080] 将前半轮迭代的先验信息初始化为0。
[0081] 将前向度量和后向度量初始化为 其中，log(·)表示以自然常数e为底的对数操作；L表示双二元Turbo码分量编码器中寄存器的总数。
[0082] 步骤3，设置译码迭代次数。
[0083] 设置译码最大迭代次数为M，M∈{1,2，...,10}，∈表示属于操作；将当前迭代次数初始化为1。
[0084] 步骤4，计算前半轮迭代的后向度量值。
[0085] 第1步，合并接收的两路系统位软信息和两路第一校验位软信息，得到两个合并后的系统位软信息和两个合并后的校验位软信息。
[0086] 所述的合并两路交织后的系统位软信息和两路第一校验位软信息是按照下式得到的：
[0087] S1＝Ark+Brk
[0088] S2＝Ark-Brk
[0089] P1＝Yrk+Wrk
[0090] P2＝Yrk-Wrk
[0091] 其中，S1和S2分别表示两个合并后的系统位软信息；P1和P2分别表示两个合并后的校验位软信息；Ark表示k时刻所接收的两路系统位软信息中的一路系统位软信息；Brk表示k时刻所接收的两路系统位软信息中的另一路系统位软信息；Yrk表示k时刻所接收的两路第一校验位软信息中的一路第一校验位软信息；Wrk表示k时刻所接收的两路第一校验位软信息中的另一路第一校验位软信息，k依次取值1,2,…,N，N表示双二元Turbo码交织器的交织深度。
[0092] 第2步，由合并后的系统位软信息和合并后的校验位软信息，得到8个信道观测值。
[0093] 所述的8个信道观测值是按照下式得到的：
[0094] SP11＝S1+P1
[0095] SP12＝S1-P1
[0096] SP13＝S1+P2
[0097] SP14＝S1-P2
[0098] SP21＝S2+P2
[0099] SP22＝S2-P2
[0100] SP23＝S2+P1
[0101] SP24＝S2-P1
[0102] 其中，SP11、SP12、SP13、SP14、SP21、SP22、SP23、SP24分别表示由合并后的系统位软信息和合并后的校验位软信息，计算得到的8个信道观测值；S1和S2分别表示两个合并后的系统位软信息；P1和P2分别表示两个合并后的校验位软信息。
[0103] 第3步，由信道观测值和前半轮迭代的先验信息，计算16个分支度量值。
[0104] 所述的16个分支度量值是按照下式得到的：
[0105]
[0106]
[0107] 其中， 表示k时刻第j个分支度量值；k表示译码时刻，k∈{1,2,...,NL}；NL表示双二元Turbo码编码码长；j表示k时刻分支度量值序号，j∈{0,1,...,15}；SP11、SP12、SP13、SP14、SP21、SP22、SP23、SP24分别表示由合并后的系统位软信息和合并后的校验位软信息计算得到的8个信道观测值； 分别表示k时刻先验信息的3个数据。
[0108] 第4步，将起始状态相同且结束状态不同的所有后向度量中间值，以四个中间值为一组，划分成多个后向度量中间值小组。
[0109] 第5步，利用下式，在每个后向度量中间值小组中，分别计算后向度量差值和分支度量差值：
[0110] X＝B1-C1
[0111] P＝C2-B2
[0112] 其中，X表示后向度量差值；B1表示带有对应于码字(0,1)的先验信息的后向度量中间值中的后向度量；C1表示带有对应于码字(1,0)的先验信息的后向度量中间值中的后向度量；P表示分支度量差值；C2表示带有对应于码字(1,0)的先验信息的后向度量中间值中的分支度量；B2表示带有对应于码字(0,1)的先验信息的后向度量中间值中的分支度量。
[0113] 第6步，判断后向度量差值X是否大于分支度量差值P，若是，则将带有对应于码字(0,1)的先验信息的后向度量中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,0)的先验信息的后向度量中间值作为第一个待选后向度量中间值。
[0114] 第7步，使用与本步骤第5步、第6步相同的方法，由带有对应于码字(0,0)的先验信息的后向度量中间值和带有对应于码字(1,1)的先验信息的后向度量中间值，获得第二个待选后向度量中间值。
[0115] 第8步，取每个后向度量中间值小组中两个待选后向度量中间值中最大的一个，作为该后向度量中间值小组选取的后向度量。
[0116] 步骤5，计算前半轮迭代的后验信息和前半轮迭代的外信息。
[0117] 第1步，将结束状态相同且起始状态不同的所有前向度量中间值，以四个中间值为一组，划分成多个前向度量中间值小组。
[0118] 第2步，利用下式，在每个前向度量中间值小组中，分别计算前向度量差值和分支度量差值：
[0119] Y＝A1-D1
[0120] Q＝D2-A2
[0121] 其中，Y表示前向度量差值；A1表示带有对应于码字(0,0)的先验信息的前向度量中间值中的前向度量；D1表示带有对应于码字(1,1)的先验信息的前向度量中间值中的前向度量；Q表示分支度量差值；D2表示带有对应于码字(1,1)的先验信息的前向度量中间值中的分支度量；A2表示带有对应于码字(0,0)的先验信息的前向度量中间值中的分支度量。
[0122] 第3步，判断前向度量差值Y是否大于分支度量差值Q，若是，则将带有对应于码字(0,0)的先验信息的前向度量的中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,1)的先验信息的前向度量中间值作为第一个待选前向度量中间值。
[0123] 第4步，使用与本步骤第2步、第3步相同的方法，由带有对应于码字(0,1)的先验信息的前向度量中间值和带有对应于码字(1,0)的先验信息的前向度量中间值，获得第二个待选前向度量中间值；
[0124] 第5步，取每个前向度量中间值小组中两个待选前向度量中间值中最大的一个，作为该前向度量中间值小组选取的前向度量。
[0125] 第6步，利用后验信息公式，计算前半轮迭代的后验信息。
[0126] 所述的后验信息公式如下：
[0127]
[0128] 其中，Lk表示k时刻的后验信息，k的取值范围为[1,2,…,N]，N表示双二元Turbo码交织器的交织深度；max(·)表示取最大值操作；s′和s分别表示双二元Turbo码编码器状态状态转移图上的两个状态节点；Auk和Buk分别表示双二元Turbo码编码器k时刻两个输入的比特； 表示k-1时刻状态s′的前向度量； 表示k时刻第j个分支度量值，j表示分支度量序号，j∈{0,1,...,15}； 表示k时刻状态s的后向度量。
[0129] 第7步，利用外信息公式，计算前半轮迭代的外信息。
[0130] 所述的外信息公式如下：
[0131]
[0132] 其中， 表示k时刻对应码字为(0,1)的前半轮迭代的外信息，k的取值范围为[1,2,…,N]，N表示双二元Turbo码交织器的交织深度； 表示k时刻对应码字为(0,1)的前半轮迭代的后验信息；Brk表示k时刻所接收的两路系统位软信息中的一路系统位软信息；
表示k时刻对应码字为(0,1)的前半轮迭代的先验信息； 表示k时刻对应码字为(1,0)的前半轮迭代的外信息； 表示k时刻对应码字为(1,0)的前半轮迭代的后验信息；Ark表示k时刻所接收的两路系统位软信息中的另一路系统位软信息； 表示k时刻对应码字为(1,
0)的前半轮迭代的先验信息； 表示k时刻对应码字为(1,1)的前半轮迭代的外信息；
表示k时刻对应码字为(1,1)的前半轮迭代的后验信息； 表示k时刻对应码字为(1,1)的前半轮迭代的先验信息。
[0133] 步骤6，获得后半轮迭代的先验信息。
[0134] 利用双二元Turbo码交织器对前半轮迭代的外信息进行交织，获得后半轮迭代的先验信息。
[0135] 步骤7，计算后半轮迭代的后向度量值。
[0136] 第1步，合并两路交织后的系统位软信息和两路第二校验位软信息，得到两个合并后的系统位软信息和两个合并后的校验位软信息。
[0137] 所述的合并两路交织后的系统位软信息和两路第二校验位软信息是按照下式得到的：
[0138] S1＝Ark+Brk
[0139] S2＝Ark-Brk
[0140] P1＝Yrk+Wrk
[0141] P2＝Yrk-Wrk
[0142] 其中，S1和S2分别表示两个合并后的系统位软信息；P1和P2分别表示两个合并后的校验位软信息；Ark表示k时刻两路交织后的系统位软信息中的一路交织后的系统位软信息；Brk表示k时刻两路交织后的系统位软信息中的另一路交织后的系统位软信息；Yrk表示k时刻所接收的两路第二校验位软信息中的一路第二校验位软信息；Wrk表示k时刻所接收的两路第二校验位软信息中的另一路第二校验位软信息，k依次取值1,2,…,N，N表示双二元Turbo码交织器的交织深度。
[0143] 第2步，由合并后的系统位软信息和合并后的校验位软信息，合并得到8个信道观测值。
[0144] 所述的8个信道观测值是按照下式得到的：
[0145] SP11＝S1+P1
[0146] SP12＝S1-P1
[0147] SP13＝S1+P2
[0148] SP14＝S1-P2
[0149] SP21＝S2+P2
[0150] SP22＝S2-P2
[0151] SP23＝S2+P1
[0152] SP24＝S2-P1
[0153] 其中，SP11、SP12、SP13、SP14、SP21、SP22、SP23、SP24分别表示由合并后的系统位软信息和合并后的校验位软信息，计算得到的8个信道观测值；S1和S2分别表示两个合并后的系统位软信息；P1和P2分别表示两个合并后的校验位软信息。
[0154] 第3步，由信道观测值和先验信息，计算得到16个分支度量值。
[0155] 所述的16个分支度量值是按照下式得到的：
[0156]
[0157]
[0158] 其中， 表示k时刻第j个分支度量值；k表示译码时刻，k∈{1,2,...,N}；N表示双二元Turbo码编码码长；j表示k时刻分支度量值序号，j∈{0,1,...,15}；SP11、SP12、SP13、SP14、SP21、SP22、SP23、SP24分别表示由合并后的系统位软信息和合并后的校验位软信息计算得到的8个信道观测值； 分别表示k时刻先验信息的3个数据。
[0159] 第4步，将起始状态相同且结束状态不同的所有后向度量中间值，以四个中间值为一组，划分成多个后向度量中间值小组。
[0160] 第5步，利用下式，在每个后向度量中间值小组中，分别计算后向度量差值和分支度量差值：
[0161] X＝B1-C1
[0162] P＝C2-B2
[0163] 其中，X表示后向度量差值；B1表示带有对应于码字(0,1)的先验信息的后向度量中间值中的后向度量；C1表示带有对应于码字(1,0)的先验信息的后向度量中间值中的后向度量；P表示分支度量差值；C2表示带有对应于码字(1,0)的先验信息的后向度量中间值中的分支度量；B2表示带有对应于码字(0,1)的先验信息的后向度量中间值中的分支度量。
[0164] 第6步，判断后向度量差值X是否大于分支度量差值P，若是，则将带有对应于码字(0,1)的先验信息的后向度量中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,0)的先验信息的后向度量中间值作为第一个待选后向度量中间值。
[0165] 第7步，使用与本步骤第5步、第6步相同的方法，由带有对应于码字(0,0)的先验信息的后向度量中间值和带有对应于码字(1,1)的先验信息的后向度量中间值，获得第二个待选后向度量中间值；
[0166] 第8步，取每个后向度量中间值小组中两个待选后向度量中间值中最大的一个，作为该后向度量中间值小组选取的后向度量。
[0167] 步骤8，计算后半轮迭代的后验信息和后半轮迭代的外信息。
[0168] 第1步，将结束状态相同且起始状态不同的所有前向度量中间值，以四个中间值为一组，划分成多个前向度量中间值小组。
[0169] 第2步，利用下式，在每个前向度量中间值小组中，分别计算前向度量差值和分支度量差值：
[0170] Y＝A1-D1
[0171] Q＝D2-A2
[0172] 其中，Y表示前向度量差值；A1表示带有对应于码字(0,0)的先验信息的前向度量中间值中的前向度量；D1表示带有对应于码字(1,1)的先验信息的前向度量中间值中的前向度量；Q表示分支度量差值；D2表示带有对应于码字(1,1)的先验信息的前向度量中间值中的分支度量；A2表示带有对应于码字(0,0)的先验信息的前向度量中间值中的分支度量。
[0173] 第3步，判断前向度量差值Y是否大于分支度量差值Q，若是，则将带有对应于码字(0,0)的先验信息的前向度量的中间值作为第一个待选前向度量中间值，否则，将带有对应于码字(1,1)的先验信息的前向度量中间值作为第一个待选前向度量中间值。
[0174] 第4步，使用与本步骤第2步、第3步相同的方法，由带有对应于码字(0,1)的先验信息的前向度量中间值和带有对应于码字(1,0)的先验信息的前向度量中间值，获得第二个待选前向度量中间值。
[0175] 第5步，取每个前向度量中间值小组中两个待选前向度量中间值中最大的一个，作为该前向度量中间值小组选取的前向度量。
[0176] 第6步，利用后验信息公式，计算后半轮迭代的后验信息。
[0177] 所述的后验信息公式如下：
[0178]
[0179] 其中，Lk表示k时刻的后验信息，k的取值范围为[1,2,…,N]，N表示双二元Turbo码交织器的交织深度；max(·)表示取最大值操作；s′和s分别表示双二元Turbo码编码器状态状态转移图上的两个状态节点；Auk和Buk分别表示双二元Turbo码编码器k时刻两个输入的比特； 表示k-1时刻状态s′的前向度量； 表示k时刻第j个分支度量值，j表示分支度量序号，j∈{0,1,...,15}； 表示k时刻状态s的后向度量。
[0180] 第7步，利用外信息公式，计算后半轮迭代的外信息。
[0181] 所述的外信息计算如下：
[0182]
[0183] 其中， 表示k时刻对应码字为(0,1)的后半轮迭代的外信息，k的取值范围为[1,2,…,N]，N表示双二元Turbo码交织器的交织深度； 表示k时刻对应码字为(0,1)的前半轮迭代的后验信息；Brk表示k时刻两路交织后的系统位软信息中的一路交织后的系统位软信息； 表示k时刻对应码字为(0,1)的后半轮迭代的先验信息； 表示k时刻对应码字为(1,0)的后半轮迭代的外信息； 表示k时刻对应码字为(1,0)的后半轮迭代的后验信息；Ark表示k时刻两路交织后的系统位软信息中的一路交织后的系统位软信息； 表示k时刻对应码字为(1,0)的后半轮迭代的先验信息； 表示k时刻对应码字为(1,1)的后半轮迭代的外信息； 表示k时刻对应码字为(1,1)的后半轮迭代的后验信息； 表示k时刻对应码字为(1,1)的后半轮迭代的先验信息。
[0184] 步骤9，获得前半轮迭代的先验信息。
[0185] 利用双二元Turbo码解交织器，对后半轮迭代的外信息进行解交织，获得前半轮迭代的先验信息。
[0186] 步骤10，判断当前迭代次数是否达到最大迭代次数，若是，则执行步骤11，否则，将当前迭代次数加1后执行步骤4。
[0187] 步骤11，判决对数似然比信息。
[0188] 使用双二元Turbo码解交织器，对后半轮迭代的后验信息进行解交织，得到对数似然比。
[0189] 使用双二元Turbo码的符号判决规则，对对数似然比进行硬判决，译码结束。
[0190] 下面结合仿真图对本发明的效果做进一步的描述。
[0191] 1、仿真实验条件：
[0192] 本发明的仿真实验软件环境：Windows 7，MATLAB R2014a。
[0193] 本发明的仿真实验采用第二代数字视频广播回传信道系统DVB-RCS2(Digital Video Broadcasting-Return Channel via Satellite 2)标准的Turbo码，设置译码最大迭代次数为8，并与一种现有技术进行了对比(增强Max-Log-MAP译码算法)。
[0194] 2、仿真内容与仿真结果分析：
[0195] 本发明的仿真实验是采用本发明的方法和现有技术的增强Max-Log-MAP译码算法，分别对1/2码率、408比特码长的双二元Turbo码和5/6码率、1504比特码长的双二元Turbo码进行译码，根据误比特率和比特信噪比的对应关系绘制图2。图2中纵坐标表示误比特率，横坐标表示比特信噪比，单位是dB。图2中以圆形标识的虚线，代表采用本发明方法对1/2码率、408比特码长的双二元Turbo码的译码仿真结果曲线。图2中以正方形标识的实线，代表采用现有技术的增强Max-Log-MAP译码算法对1/2码率、408比特码长的双二元Turbo码的译码仿真结果曲线。图2中以五角星标识的虚线，代表采用本发明方法对5/6码率、1504比特码长的双二元Turbo码的译码仿真结果曲线。图2中以菱形标识的实线，代表采用现有技术的增强Max-Log-MAP译码算法对5/6码率、1504比特码长的双二元Turbo码的译码仿真结果曲线。图2中以上三角标识的实线，代表未编码情况的误码仿真曲线。
[0196] 从图2的仿真结果可以看出，采用本发明方法与现有技术的增强Max-Log-MAP译码算法，对1/2码率、408比特码长的双二元Turbo码译码的仿真结果曲线比较接近。采用本发明方法与现有技术的增强Max-Log-MAP译码算法，对5/6码率、1504比特码长的双二元Turbo码译码的仿真结果曲线比较接近。与未编码的方法相比，本发明方法的所有仿真曲线均发生了左移，说明了本发明方法具有优秀的译码性能。
[0197] 统计本发明方法和现有技术的增强Max-Log-MAP译码算法，获得一次双二元Turbo的前向度量计算和后向度量过程中，所使用的加法次数和比较次数，根据统计的加法次数和比较次数绘制表1。
[0198] 表1译码复杂度表
[0199]
[0200] 由表1所示的译码复杂度可以看到，采用现有技术的增强Max-Log-MAP译码算法，获得一次双二元Turbo的前向度量计算和后向度量过程中，使用加法的总次数为248，使用比较的总次数为96。采用本发明方法，获得一次双二元Turbo的前向度量计算和后向度量过程中，使用加法的总次数为134，使用比较的总次数为96。与现有技术的增强Max-Log-MAP译码算法相比，在获得一次双二元Turbo的前向度量计算和后向度量过程中，本发明方法使用加法次数占现有技术的增强Max-Log-MAP译码算法使用加法次数的54％，由此说明本发明具有很低的译码复杂度。
