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

传真(Fax)编程

开发平台:

Matlab

  1. function [soft_out,ex_info]=com_decoder_linear_log_map(in,app)
  2. %****************************************************************
  3. % 内容概述:子译码器。
  4. %          利用硬件化的方式实现TURBO码的Linear_LOG-MAP译码
  5. %          线性近似
  6. %          生成矩阵按照3GPP标准为[1 1 0 1;1 0 1 1]
  7. %          输入为经过高斯信道的RSC软输入,而输出为软输出和外部信息
  8. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  9. % 单    位:南京邮电大学,通信工程系
  10. % 创建时间:2005年9月3日
  11. % 修改时间:
  12. % 参考文献:《数字通信--基础与应用》
  13. %          《改进的Turbo码算法及其FPGA实现过程的研究》,天津大学,张宁,赵雅兴
  14. %           后向度量的计算方式和decoder_3GPP_max_th.m有所不同
  15. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  16. %****************************************************************
  17. x=in(1,:);              %输入系统位
  18. y=in(2,:);              %输入校验位
  19. in_length=length(in);
  20. %---初始化&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
  21. Infty = -128;
  22. d=zeros(8,2,in_length);     %分支量度,8种可能状态,输入为0或者1
  23.                             %D(S,i,k)
  24. a=Infty*ones(8,in_length);  %前向分支量度,A(S,k)
  25. a(1,1)=0;                   %寄存器状态由全零开始
  26. b=Infty*ones(8,in_length+1);%后向分支量度,B(S,k)
  27. b(1,in_length)=0;           %寄存器状态由全零结束
  28. %---计算度量和LLR&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
  29. for k=1:in_length
  30.     d(1,2,k)=x(k)+y(k)+app(k);
  31.     d(2,2,k)=d(1,2,k);
  32.     d(7,2,k)=d(1,2,k);
  33.     d(8,2,k)=d(1,2,k);
  34.     
  35.     d(3,2,k)=x(k)+app(k);
  36.     d(4,2,k)=d(3,2,k);
  37.     d(5,2,k)=d(3,2,k);
  38.     d(6,2,k)=d(3,2,k);
  39.     d(3,1,k)=y(k);
  40.     d(4,1,k)=d(3,1,k);
  41.     d(5,1,k)=d(3,1,k);
  42.     d(6,1,k)=d(3,1,k);
  43.     
  44.     if k>1
  45.         a(1,k)=max_lin(a(1,k-1),(a(2,k-1)+d(2,2,k-1)));
  46.         a(2,k)=max_lin((a(4,k-1)+d(4,1,k-1)),(a(3,k-1)+d(3,2,k-1)));
  47.         a(3,k)=max_lin((a(5,k-1)+d(5,1,k-1)),(a(6,k-1)+d(6,2,k-1)));
  48.         a(4,k)=max_lin(a(8,k-1),(a(7,k-1)+d(7,2,k-1)));
  49.         a(5,k)=max_lin(a(2,k-1),(a(1,k-1)+d(1,2,k-1)));
  50.         a(6,k)=max_lin((a(3,k-1)+d(3,1,k-1)),(a(4,k-1)+d(4,2,k-1)));
  51.         a(7,k)=max_lin((a(6,k-1)+d(6,1,k-1)),(a(5,k-1)+d(5,2,k-1)));
  52.         a(8,k)=max_lin(a(7,k-1),(a(8,k-1)+d(8,2,k-1)));       
  53.     end
  54.     
  55.     if k==in_length
  56.         b(1,k)=max_lin(b(1,k+1),(b(5,k+1)+d(1,2,k)));
  57.         b(2,k)=max_lin(b(5,k+1),(b(1,k+1)+d(2,2,k)));
  58.         b(3,k)=max_lin((b(6,k+1)+d(3,1,k)),(b(2,k+1)+d(3,2,k)));
  59.         b(4,k)=max_lin((b(2,k+1)+d(4,1,k)),(b(6,k+1)+d(4,2,k)));
  60.         b(5,k)=max_lin((b(3,k+1)+d(5,1,k)),(b(7,k+1)+d(5,2,k)));
  61.         b(6,k)=max_lin((b(7,k+1)+d(6,1,k)),(b(3,k+1)+d(6,2,k)));
  62.         b(7,k)=max_lin(b(8,k+1),(b(4,k+1)+d(7,2,k)));
  63.         b(8,k)=max_lin(b(4,k+1),(b(8,k+1)+d(8,2,k)));
  64.         %计算LLR--------------------------------------
  65.         l(k)=max_lin_plus([...
  66.             (a(1,k)+d(1,2,k)+b(5,k+1)),(a(2,k)+d(2,2,k)+b(1,k+1)),...
  67.             (a(3,k)+d(3,2,k)+b(2,k+1)),(a(4,k)+d(4,2,k)+b(6,k+1)),...
  68.             (a(5,k)+d(5,2,k)+b(7,k+1)),(a(6,k)+d(6,2,k)+b(3,k+1)),...
  69.             (a(7,k)+d(7,2,k)+b(4,k+1)),(a(8,k)+d(8,2,k)+b(8,k+1))...
  70.             ])-max_lin_plus([...
  71.             (a(1,k)+b(1,k+1)),(a(2,k)+b(5,k+1)),...
  72.             (a(3,k)+d(3,1,k)+b(6,k+1)),(a(4,k)+d(4,1,k)+b(2,k+1)),...
  73.             (a(5,k)+d(5,1,k)+b(3,k+1)),(a(6,k)+d(6,1,k)+b(7,k+1)),...
  74.             (a(7,k)+b(8,k+1)),(a(8,k)+b(4,k+1))...
  75.             ]);
  76.     end
  77. end
  78. for k=in_length-1:-1:1
  79.         b(1,k)=max_lin(b(1,k+1),(b(5,k+1)+d(1,2,k)));
  80.         b(2,k)=max_lin(b(5,k+1),(b(1,k+1)+d(2,2,k)));
  81.         b(3,k)=max_lin((b(6,k+1)+d(3,1,k)),(b(2,k+1)+d(3,2,k)));
  82.         b(4,k)=max_lin((b(2,k+1)+d(4,1,k)),(b(6,k+1)+d(4,2,k)));
  83.         b(5,k)=max_lin((b(3,k+1)+d(5,1,k)),(b(7,k+1)+d(5,2,k)));
  84.         b(6,k)=max_lin((b(7,k+1)+d(6,1,k)),(b(3,k+1)+d(6,2,k)));
  85.         b(7,k)=max_lin(b(8,k+1),(b(4,k+1)+d(7,2,k)));
  86.         b(8,k)=max_lin(b(4,k+1),(b(8,k+1)+d(8,2,k)));
  87.         %计算LLR--------------------------------------
  88.         l(k)=max_lin_plus([...
  89.             (a(1,k)+d(1,2,k)+b(5,k+1)),(a(2,k)+d(2,2,k)+b(1,k+1)),...
  90.             (a(3,k)+d(3,2,k)+b(2,k+1)),(a(4,k)+d(4,2,k)+b(6,k+1)),...
  91.             (a(5,k)+d(5,2,k)+b(7,k+1)),(a(6,k)+d(6,2,k)+b(3,k+1)),...
  92.             (a(7,k)+d(7,2,k)+b(4,k+1)),(a(8,k)+d(8,2,k)+b(8,k+1))...
  93.             ])-max_lin_plus([...
  94.             (a(1,k)+b(1,k+1)),(a(2,k)+b(5,k+1)),...
  95.             (a(3,k)+d(3,1,k)+b(6,k+1)),(a(4,k)+d(4,1,k)+b(2,k+1)),...
  96.             (a(5,k)+d(5,1,k)+b(3,k+1)),(a(6,k)+d(6,1,k)+b(7,k+1)),...
  97.             (a(7,k)+b(8,k+1)),(a(8,k)+b(4,k+1))...
  98.             ]);
  99. end
  100. soft_out=l;
  101. ex_info=soft_out-app-x;