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

传真(Fax)编程

开发平台:

Matlab

  1. function [hard_out,soft_out]=deturbo(in,alphain)
  2. %****************************************************************
  3. % 内容概述:turbo解码器,in是RSC编码器输出
  4. %          利用硬件化的方式实现TURBO码的p-MAX-LOG-MAP译码
  5. %          生成矩阵按照3GPP标准为[1 1 0 1;1 0 1 1]
  6. %          未使用另外一个译码器反馈的外部信息
  7. %          输入为经过高斯信道的RSC软输入,而输出为软输出
  8. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  9. % 单    位:南京邮电大学,通信工程系
  10. % 创建时间:2005年7月17日
  11. % 修改时间:
  12. % 参考文献:《数字通信--基础与应用》
  13. %          《High performace parallelised 3GPP Turbo Decoder》
  14. %          《改进的Turbo码算法及其FPGA实现过程的研究》,天津大学,张宁,赵雅兴
  15. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  16. %****************************************************************
  17. L_seq=length(in);
  18. in1=in(1:2,:);
  19. in2(1,:)=in(1,alphain);
  20. in2(2,:)=in(3,:);
  21. iter_lim=1; % 迭代次数
  22. for it=1:iter_lim
  23.     if it==1
  24.         a_p_first(1:L_seq)=0;
  25.         [so1,ep1] = com_decoder(a_p_first,in1);
  26.     else
  27.         a_p1(alphain)=ep2;
  28.         [so1,ep1] = com_decoder(a_p1,in1);
  29.     end
  30.     a_p2=ep1(alphain);
  31.     [so2,ep2] = com_decoder(a_p2,in2);    
  32. end
  33. % 解码结束,输出--------------------
  34. soft_out(alphain)=so2;
  35. for i=1:L_seq
  36.     if soft_out(i)>=0
  37.         hard_out(i)=1;
  38.     else
  39.         hard_out(i)=-1;
  40.     end
  41. end