deturbo_new.m
上传用户:hnyfjx
上传日期:2013-06-30
资源大小:2149k
文件大小:2k
源码类别:

传真(Fax)编程

开发平台:

Matlab

  1. function [hard_out,soft_out]=deturbo_new(in,alphain)
  2. %****************************************************************
  3. % 内容概述:turbo解码器,in是RSC编码器输出
  4. %          与一般的解码器相比,这个子解码器多出了对是否第一个子解码器的判断
  5. %          如果是第一个子解码器,输出的外部信息中包含了输入系统信息比特
  6. %          提供给第二个子解码器使用时不需要将系统位删除,只需和后验信息
  7. % 一同交织后交给第二个子解码器,这样就省掉了专门需要对系统信息进行的交织!!
  8. %          利用硬件化的方式实现TURBO码的p-MAX-LOG-MAP译码
  9. %          生成矩阵按照3GPP标准为[1 1 0 1;1 0 1 1]
  10. %          未使用另外一个译码器反馈的外部信息
  11. %          输入为经过高斯信道的RSC软输入,而输出为软输出
  12. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  13. % 单    位:南京邮电大学,通信工程系
  14. % 创建时间:2005年7月18日
  15. % 修改时间:
  16. % 参考文献:《数字通信--基础与应用》
  17. %          《High performace parallelised 3GPP Turbo Decoder》
  18. %          《改进的Turbo码算法及其FPGA实现过程的研究》,天津大学,张宁,赵雅兴
  19. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  20. %****************************************************************
  21. L_seq=length(in);
  22. in1=in(1:2,:);
  23. in2(1,:)=in(1,alphain);
  24. in2(2,:)=in(3,:);
  25. iter_lim=1; % 迭代次数
  26. for it=1:iter_lim
  27.     %---component decoder1
  28.     is_com1=1;
  29.     if it==1
  30.         a_p1(1:L_seq)=0;
  31.         [so1,ep1] = com_decoder_new(a_p1,in1,is_com1);
  32.     else
  33.         a_p1(alphain)=ep2;
  34.         [so1,ep1] = com_decoder_new(a_p1,in1,is_com1);
  35.     end
  36.     %---component decoder2
  37.     is_com1=0;
  38.     a_p2=ep1(alphain);
  39.     [so2,ep2] = com_decoder_new(a_p2,in2,is_com1);    
  40. end
  41. % 解码结束,输出--------------------
  42. soft_out(alphain)=so2;
  43. for i=1:L_seq
  44.     if soft_out(i)>=0
  45.         hard_out(i)=1;
  46.     else
  47.         hard_out(i)=-1;
  48.     end
  49. end