Turbo码的编译码算法仿真讲解.docx

1、Turbo码的编译码算法仿真讲解电子信息类实践课III通信系统仿真题 目Turbo码的编译码算法仿真专 业 学 号 姓 名 日 期 注：本报告仅供参考1、课程设计目的（黑体小三，段前段后个一行）通过完成在在衰落信道下采用不同调制信号进行Turbo码编译码的编程实现，进一步了解了 Turbo码的编码解码过程，以及在不同调制方式不同信道下的性能比见。通过对卷积和交织器的设计，深入了解卷积和交织的作用。以及熟悉了通信 仿真的整体流程。2、课程设计内容具体叙述课程设计的主要内容和原理。表1主要课程设计内容列表课程设计内容负责人Turbo编码译码过程、信道设计、各模块肖雨桐卷积、交织杨佳能不同调制信号的

2、设计熊风在进行本次Turbo仿真时，采用了两种不同编程方式。在程序一中是直接调用matlab Communications System Toolbox中的Turbo编码和解码工具箱，通过配置参数进行仿真。而 在程序二中则根据 Turbo码编译码原理编写。如果程序一更像是一个黑匣子，只能知道通过编解码模块前、后的数据，而具体做了哪些则 不得而 知。a.编码输入信编码器I編码输出息数据图1 Turbo码编码器结构典型的Turbo码编码器结构框图如图所示：由两个反馈的编码器通过一个交织 器并行连接而成。如果必要，由成员编码器输出的序列经过删余阵，从而可以产生一系列不同码率的 码。例如，对于生成矩阵

3、为g-Lgi,g2的（2 , 1, 2）卷积码通过编码后，如果进行删余， 则得到码率为1/2的编码输出序列；如果不进行删余，得到的码率为1/3。一般情况下， Turbo码成员编码器是RSC编码器。原因在于递 归编码器可以改善码的比特误码率性能。通俗理解1/2码率就是信号中有一半都是“无用信号”，这些“无用信号”就是两个分 量编码器的生成的校验码，而删余则是各删除一部分校验码，把剩下的再和信息比特合在一 块，形成编码好的矩阵。b.译码图2 Turbo译码结构Turbo码获得优异性能的根本原因之一是使用了迭代译码，通过与分量编码器 对应的分 量译码器之间软信息的交换来提高译码性能。对于Turbo码

4、这样的并行级联码，如果分量译 码器的输出为硬判决，则不可能实现分量译码器之间软信息的交换，从而限制了系统性能的 进一步提高。从信息论的角度来看，任何硬判决都会损失部分信息，因此，如果分量译码器 能够提供一个反映其输出可靠性的软输出，则其他分量译码器也可以采用软判决译码，从而 系统的性能可以得到进一步提高。为此，人们提出了软输入软输出译码(SISO)的概念和方 法。Turbo码的分 量码SISO译码算法总体上可分为MAP和SOVA两类主要算法。其中MAP算法 是一种最佳后验概率算法。SOVA类算法主要包括软输出的维特比算法(SOVA)和连续列表输 出维特比算法(SLVA ) oSOVA算法实际上

5、就是veterbi算法的扩展，保留的路径数不是两条而是和状态数有关从 所有可行的路径中找最优的是最大似然序列检测。d.衰落信道由于多径和移动台运动等影响因素，使得移动信道对传输信号在时间、频率 和角度上造成了色散，如时间色散、频率色散、角度色散等等。根据不同无线环境，接收信 号包络一般服从几种典型分布，如瑞利分布、莱斯分布和 Nakagami-m分布。在仿真衰落信道时，最重要的参数是多径扩展和多普勒频移。通常在离基站较远、反射物较多的地区，发射机和接收机之间没有直射波路径，存在大量反射波；到达接收天线的方向角随机且在（ 0i）均匀分布；各反射波的幅度和相位都统计独立。3、设计与实现过程信号生成

6、图4、程序一框图具体实现过程：按照流程图中的各方面(模块)内容进行代码级的详细说明，例如：衰落信道的设计在进行仿真的过程中尝试使用了两种不同的编写方式：a.Create Rayleigh fading channel objectchan_ray = rayleighcha n(1/10000, 100);fadedSig 二 filter(cha n_ ray, modSig nal) ; % Apply the cha nnel effectshChan 二 comm. AWGNCha nn el(NoiseMethod, Sig nal to n oise ratio (SNR); hC

7、ha n. SNR 二 EbNo_db( n);fadedSig 二 filter (cha n_ ray, modSig nal); % Apply the cha nnel effects receivedSig nal 二 step(cha n, fadedSig) ; % Apply the cha nnel effectsb.调用改进JAKES模型产生单径平坦型瑞利衰落信道子程序n samp 二 8; %脉冲抽样点数 ts = 1/(num*n samp) ;%抽样时间间隔t 二(0: num* nsamp-l)*ts;% 抽样时间序列h = rayleigh(10,t);%调用瑞利

8、衰落子程序，输入为(最大多普勒频移，抽样时间序列)modSig nal 二 rectpulse(modSig nal, nsamp) ;%巨形脉冲形成modSig nal = h. *modSig nal;% 通过瑞利信道 receivedSig nal = in tdump (receivedSig nal, nsamp) : %匹配滤波相干解调使用mat lab函数rayleighchan程序一：a调用turbo编码器解码器encoder =comm. TurboE ncoder( 1 nterleaverl ndicesSource，, 1 nputport);decoder =comm

9、. TurboDecoder C 1 nterleaverl ndicesSource，, 1 nputport,Numlteratio ns, 4);b.调用AWGN信道cha n = comm. AWGNCha nn el( EbNo, EbNo_db, BitsPerSymbolJ, log2(M);c.进行编码解码Turbo 编码： encodedDa ta 二 st ep (encoder, msg, index);调制编码过宿道 receivedSignal 二 st ep (chan, modSignal):解调Turbo 译码： receivedB its 二 st ep (d

10、ecoder, -demodSignal, index):程序二code_le ngth = 1024;%码长rate = 1/2;% 码率niter = 4; %迭代次数a信道：AWGNEbN0_db = 2:6;en 二 10A(EbN0_db( nEN)/10);L_c二4*en *rate;%信道置信度sigma = 1/sqrt (2*rate*en) ;%AWGN 信道标准并 过信道：r = en coder_out + sigmara ndn (1, code_le ngth*(2);b.Turbo编码器%第一个分量RSC编码output1 二 rsc_en code(Gmsg,

11、 1);%1*2048y(l, :) = output 1 (1:2:2*code_le ngth) ;% 系统比特y(2, :) = output 1 (2:2:2*code_le ngth) ;% 校验比特%第二个分量RSC编码a = reshape(y(1, :),32,32);y_inv 二 reshape(a5 , 1, 1024);output2 二 rsc_en code (Gv_i nv, -1) ;% 输入 1024 输出 1*2048y(3, :) = output2(2:2:2*code_le ngth) ;% 校验比特输出删余生成码率为1/2的码en coder_out

12、(l:2:e nd)二 y(1,:);en coder_out(2:4:e nd)二 y(2, l:2:e nd);en coder_out (4:4:e nd) = y (3, 2:2:e nd) ;%删余，奇为系统比特c.Turbo译码器通过解复用生成每个分量译码器的译码输入数据；初始化外部信息L_e(l:code_length) 后。根据迭代次数，循环依次进行第一分量译码、第二分量译码。%第一分量译码器L_a(i ndex) = L_e; %先验信息L_all 二 sova(rec_s (1, :), G, L_a, 1);L_e = L_all2*rec_s(1, 1:2:2*code

13、_le ngth)L_a;% 计算外部信息%第二分量译码器L_a 二 L_e(i ndex);L_all 二 sova (rec_s (2, :), G, L_a, 2);L_e 二 L_all-2*rec_s(2, 1:2:2*code_le ngth)_L_a;最后再估计信息比特值。d.软输入Viterbi算法初始化时，建立网格图，网格中的零状态初始化度量为 0,其他为负无穷。输入先验信息L、RSC码生成矩阵以及数据；前向跟踪计算所有路径度量值；第一分量译码 器，从零状态开始向后跟踪；第二分量译码器从最可能的状态开始后向跟踪；后向跟踪获得 估计比特及最大似然路径；找到相应于不同信息比特估计

14、的竞争路径的最小delta值，得到 软输出。(具体可参看刘向东、向良军编著的信道编码与matlab仿真除SOVA算法外里面 还给出了 Log_MAP算法)4、结论在进行仿真时，发现上述两种衰落信道误码率基本都是 50%左右，也就意味着基本上全部都是错的。在进行了分析之后发现，在不加其他模块进行仿真的时候，光对比 经过衰落信道前后的比特流基本就能错一半左右。O210 话0山一ffi10JTurbo Uncoded10”0 Eb/No (dB)图5衰落信道下程序1仿真结果ber =0M956 0 + 4S03 Q. 47B2 0. 4694图6仿真得出ber值因无法仿真出较好的衰落信道，因此两个程序都放弃衰落信道改用加性高斯白噪声信道。图7程序1结果在16QAM调制信号时在1dB便能达到接近10_6o程序二是1/2码率，采用BPSK调制、sova译码算法，需要到9dB时，误码率大概能达到两者程序出现如此大的结果差异，最主要的可能与Turbo译码方式相关。成绩评定表成绩评疋项目比例得分报告成绩（百分制）20%答辩成绩（百分制）80%总评成绩（百分制）100%评语：评审教师： 时 间:注：课题设计成绩评定表单独一页，装订在报告的最后