Turbo码

第6章系统分析Turbo码•13.1Turbo码的提出•13.2Turbo码编码器的组成•13.3Turbo码的译码•13.4Turbo码的分量码、交织器与性能限•13.5Turbo码在实际通信系统(3GPP)中的应用•习题第6章系统分析§13.1Turbo码的提出•Turbo码，又称并行级联卷积码(PCCC)，是由C.Berrou等在ICC'93会议上提出的。它巧妙地将卷积码和随机交织器结合在一起，实现了随机编码的思想；同时，采用软输出迭代译码来逼近最大似然译码。模拟结果表明，如果采用大小为65535的随机交织器，并且进行18次迭代，则在Eb/N0≥0.7dB时，码率为1/2的Turbo码在AWGN信道上的误比特率（BER）≤10-5，达到了近Shannon限的性能（1/2码率的Shannon限是0dB）。第6章系统分析•因此，这一超乎寻常的优异性能，立即引起信息与编码理论界的轰动。图13-1中给出了Turbo码及其它编码方案的性能比较，从中可以看出Turbo编码方案的优越性。第6章系统分析图13-1AWGN信道中的码率与Shannon限Shannon限QPSK限1.210.80.60.40.20－20246810容量(比特/符号)Eb/N0(dB)1.4BCH(255,23)cc(2,7)并用ML解码RS(255,223)＋cc(2,1,7)Turbo码第6章系统分析•由于Turbo码的上述优异性能并不是从理论研究的角度给出的，而仅是计算机仿真的结果。因此，Turbo码的理论基础还不完善。后来经过不少人的重复性研究与理论分析，发现Turbo码的性能确实是非常优异的。因此，turbo码的发现，标志着信道编码理论与技术的研究进入了一个崭新的阶段，它结束了长期将信道截止速率作为实际容量限的历史。第6章系统分析•需要说明的是，由于原Turbo编译码方案申请了专利，因此在有关Turbo码的第一篇文章中，作者没有给出如何进行迭代译码的实现细节，只是从原理上加以说明。此后，P.Robertson对此进行了探讨，对译码器的工作原理进行了详细说明。人们依此进行了大量的模拟研究。第6章系统分析•Turbo码的提出，更新了编码理论研究中的一些概念和方法。现在人们更喜欢基于概率的软判决译码方法，而不是早期基于代数的构造与译码方法，而且人们对编码方案的比较方法也发生了变化，从以前的相互比较过渡到现在的均与Shannon限进行比较。同时，也使编码理论家变成了实验科学家。第6章系统分析§13.2Turbo码编码器的组成•Turbo码编码器是由两个反馈的系统卷积码编码器通过一个随机交织器并行连接而成的，编码后的校验位经过删余阵，从而产生不同码率的码字，见图13-2。第6章系统分析图13-2Turbo码编码器结构框图分量码编码器(RSC2)交织器u1分量码编码器(RSC1)删余复用信息序列uXp2Xp1XpXsX第6章系统分析•图13-2所示的是典型的Turbo码编码器结构框图，信息序列u={u1,u2,…,uN}经过一个N位交织器，形成一个新序列u1={u′1,u′2,…,u′N}（长度与内容没变，但比特位置经过重新排列）。u与u1分别传送到两个分量码编码器（RSC1与RSC2）。一般情况下，这两个分量码编码器结构相同，生成序列Xp1与Xp2。为了提高码率，序列Xp1与Xp2需要经过删余器，采用删余（puncturing）技术从这两个校验序列中周期地删除一些校验位，形成校验位序列Xp。第6章系统分析•Xp与未编码序列Xs经过复用调制后，生成了Turbo码序列X。例如，假定图13-2中两个分量编码器的码率均是1/2，为了得到1/2码率的Turbo码，可以采用这样的删余矩阵：P=［10,01］，即删去来自RSC1的校验序列Xp1的偶数位置比RSC2的校验序列Xp2的奇数位置比特。第6章系统分析•例13.1一个码率为1/3的Turbo码编码器的组成框图如图13-3所示。图13-3一个码率为1/3的Turbo码编码器v2Interleaverv1v0cc~第6章系统分析•图13-3所示的是基于(2,1,4)RSC(递归卷积系统码)的Turbo码编码器。分量码是码率为1/2的寄存器级数为4的(2,1,4)RSC码。其生成矩阵为432411,1)(DDDDDDG(13.2.1)第6章系统分析•我们假设输入序列为•c=(1011001)(13.2.2)•则第一个分量码的输出序列为•v0=(1011001)•v1=(1110001)(13.2.3)•假设经过交织器后信息序列变为•c=(1101010)(13.2.4)第6章系统分析•第二个分量码编码器所输出的校验位序列为•v2=(1000000)(13.2.5)•则Turbo码序列为•v=(111,010,110,100,000,000,110)(13.2.6)第6章系统分析§13.3Turbo码的译码•一、Turbo码的迭代译码原理•由于Turbo码是由两个或多个分量码经过不同交织后对同一信息序列进行编码，对任何单个传统编码，通常在译码器的最后得到硬判决译码比特，然而Turbo码译码算法不应局限于在译码器中通过的是硬判决信息。为了更好的利用译码器之间的信息，译码算法所用的应当是软判决信息而不是硬判决。第6章系统分析•一个由两个分量码构成Turbo码的译码器是由两个与分量码对应的译码单元和交织器与解交织器组成的，将一个译码单元的软输出信息作为下一个译码单元的输入；为了获得更好的译码性能，将此过程迭代数次。这就是Turbo码译码器的基本的工作原理。•第6章系统分析•二、Turbo码译码器的组成•Turbo码译码器的基本结构如图13-4所示。它由两个软输入软输出（SISO）译码器DEC1和DEC2串行级联组成，交织器与编码器中所使用的交织器相同。译码器DEC1对分量码RSC1进行最佳译码，产生关于信息序列u中每一比特的似然信息，并将其中的“新信息”经过交织送给DEC2，译码器DEC2将此信息作为先验信息，对分量码RSC2进行最佳译码，产生关于交织后的信息序列中每一比特的似然比信息，然后将其中的“外信息”经过解交织送给DEC1，进行下一次译码。第6章系统分析•这样，经过多次迭代，DEC1或DEC2的外信息趋于稳定，似然比渐近值逼近于对整个码的最大似然译码，然后对此似然比进行硬判决，即可得到信息序列u的每一比特的最佳估值序列。uˆ第6章系统分析图13-4Turbo码译码器的结构＋解交织判决解交织交织交织延时L(un)Le21w(ys)Le12y2py1pypysLe21w()~Le21~uk^软输入软输出译码器DEC1软输入软输出译码器DEC2多路信号分离器第6章系统分析•假定Turbo码译码器的接收序列为y=(ys,yp)，冗余信息yp经解复用后，分别送给DEC1和DEC2。于是，两个软输出译码器的输入序列分别为：•dec1:y1=(ys,y1p)•dec2:y2=(ys,y2p)第6章系统分析•三、分量码的最大后验概率译码（MAP算法）•考虑图13-5所示的软输入软输出（SISO）译码器，它能为每一译码比特提供对数似然比输出。第6章系统分析图13-5软输入软输出译码器框图L(uk)Le(uk)MAP译码器yksykp第6章系统分析•图13-5中，MAP译码器的输入序y=yN1=(y1,y2,…,yk,…,yN)。其中，yk=(ysk,ypk)。Le(uk)是关于uk的先验信息，L(uk)是关于uk的对数似然比。它们的定义如下：)|0()|1(ln)()0()1(ln)(11NkNkkkkkeyuPyuPuLuPuPuL(13.3.1)(13.3.2)第6章系统分析•假定发送端RSC编码器的存储级数为v，约束长度为K，编码器在k时刻的Sk=(ak,ak-1,…,•ak-v+1)，编码输出序列为x=(xs,xp)。传输信道模型如图13-6所示。•从图13-6可知：pkspkpkpkpkpkskssksksksksknExacaynExacay)12()12((13.3.3)(13.3.4)第6章系统分析•式中，ask和apk为信道衰落因子，对于AWGN信道，ask=apk=1;nsk和npk是两个独立同分布的高斯噪声样值，它们的均值为0，方差σ2=N0/2。图13-6信道模型QPSK调制器xksxkpcksckpaksnksakpnkpyksykp第6章系统分析•MAP译码器的任务就是求解式（13.3.3），然后按照下列规则进行判决：0)(0)(01ˆkkkuLuLu(13.3.5)下面利用BCJR算法对式（13.3.2）的计算方法进行推导。第6章系统分析•根据Bayes规则，式（13.3.2）可以写为)(/),,()(/),,(ln)(/),0()(/),1(ln)(1110),(1110),(1111NNkkussNNkkussNNkNNkkypysSsSpypysSsSpypyuPypyuPuLkk(13.3.6)第6章系统分析•式中，求和是对所有由uk=1（或uk=0）引起的Sk-1→Sk的状态转移进行的。根据BCJR算法［12］，•p(Sk-1=s′,Sk=s,yN1)可以按下式计算：•p(s′,s,yN1)=p(s′,yk-11)·p(s,yk|s′)·p(yNk+1|s)•=αk-1(s′)·γk(s′,s)·βk(s)(13.3.7)第6章系统分析•式中：•αk(s)≡p(Sk=s,yk1)为前向递推；•βk(s)≡p(yNk+1|Sk=s)为后向递推;•γk(s′,s)≡p(Sk=s,yk|Sk-1=s′)为s′和s之间的分支转移概率。第6章系统分析•下面求αk(s),βk(s)和γk(s′,s)。由定义得),|,(),(),,()(11111111kkkkkkskkkskysSysSpysSpysSsSps(13.3.8)第6章系统分析•考虑到RSC编码器等效于一个马尔可夫源，在状态Sk-1已知时，在k-1时刻以后发生的事件与以前输入无关。因此，从式（13.3.8）可得前向递推公式),()()|,()()(111ssssSysSpsskkskkkksk(13.3.9)第6章系统分析),()()|,()|()|,()(11111ssssSysSpsSypsSysSpskkskkkkNkskNkksk同样，βk(s)可按下式反向递推得到(13.3.10)第6章系统分析•至于分支转移概率γk(s′,s)，从其定义可得•γk(s′,s)=p(Sk=s|Sk-1=s′)·p(yk|Sk=s,Sk-1=s′)•=P(uk)·p(yk|uk)(13.3.11)•式中，P(uk)是uk的先验概率，p(yk|uk)由信道转移概率决定。第6章系统分析•考虑到式（13.3.11）是从连续随机变量的概率密度计算得到的，γk(s′,s)的值可能大于1，这会使得式（13.3.9）和式（13.3.10）产生溢出，导致整个算法不稳定。因此，有必要对αk(s)、βk(s)进行归一化处理。•令)()()(~1kkkypsas(13.3.12)第6章系统分析•因为，所以),()(11kkskysSpyp)()()(~sasaskskk(13.3.13)将式（13.3.9）代入上式，并且分子分母同除以p(yk1)，得到),()(~),()(~)(/),()()(/),()()(~11111111ssssssypsssypssssskksskkkskkskskkk(13.3.14)第6章系统分析•对于，考虑到•，于是有)(~sk)(/)()|()|(1111111kkkNkkNkypypyypyyp),()(