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

传真(Fax)编程

开发平台:

Matlab

  1. function [hard_out,soft_out]=decoder_SemiTh(in,alphain,num_iterate)
  2. %****************************************************************
  3. % 内容概述:turbo解码器,in是RSC编码器输出
  4. %          利用硬件化的方式实现TURBO码的MAX-LOG-MAP译码
  5. %          生成矩阵按照3GPP标准为[1 1 0 1;1 0 1 1]
  6. %          输入为经过高斯信道的RSC软输入,而输出为软、硬输出
  7. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  8. % 单    位:南京邮电大学,通信工程系
  9. % 创建时间:2005年9月12日
  10. % 修改时间:2005年10月29日
  11. % 参考文献:《数字通信--基础与应用》
  12. %          《High performace parallelised 3GPP Turbo Decoder》
  13. %          《改进的Turbo码算法及其FPGA实现过程的研究》,天津大学,张宁,赵雅兴
  14. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  15. %****************************************************************
  16. m=3;    %寄存器个数
  17. L_seq=length(in);
  18. in1=in(1:2,:);
  19. in2=in(3:4,:);
  20. e_p=zeros(1,L_seq);
  21. for it=1:num_iterate
  22.     a_p(alphain)=e_p(1:L_seq-m);  %解交织
  23.     a_p(L_seq-m+1:L_seq)=e_p(L_seq-m+1:L_seq);  %尾比特不变化
  24.     [so,e_p] = constituent_decoder_SemiTh(in1,a_p);
  25.     
  26.     a_p(1:L_seq-m)=e_p(alphain);  %交织
  27.     a_p(L_seq-m+1:L_seq)=e_p(L_seq-m+1:L_seq);  %尾比特不变化
  28.     [so,e_p] = constituent_decoder_SemiTh(in2,a_p);    
  29. end
  30. % 解码结束,输出--------------------
  31. soft_out(alphain)=so(1:L_seq-m);
  32. hard_out=(sign(soft_out)+1)/2;
  33. end