开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 传真(Fax)编程 > 查看源码
com_decoder_logmap_justforQ.m
资源名称：Turbomatlab.rar　[点击查看]
上传用户：hnyfjx
上传日期：2013-06-30
资源大小：2149k
文件大小：6k
源码类别：

传真(Fax)编程

开发平台：

Matlab

com_decoder_logmap_justforQ.m：源码内容
  1. %function [soft_out,ex_info]=com_decoder_Th(in,app)
  2. %****************************************************************
  3. % 内容概述：子译码器。
  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月3日
  10. % 修改时间：
  11. % 参考文献：《数字通信－－基础与应用》
  12. %          《改进的Turbo码算法及其FPGA实现过程的研究》,天津大学，张宁，赵雅兴
  13. %           后向度量的计算方式和decoder_3GPP_max_th.m有所不同
  14. % 版权声明：任何人均可复制、传播、修改此文件，同时需保留原始版权信息。
  15. %****************************************************************
  16. clear 
  17. o_in=[1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  18.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  19.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  20.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  21.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  22.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  23.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  24.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  25.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  26.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  27.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  28.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  29.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  30.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  31.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  32.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  33.     1 0 1 0 0 0 1 0 1 1 0 1 0 1 0 1 0 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1,...
  34.     1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 1 0 1 0 1 1 0 1 0 0 0 0 0 1 0 1 1];
  35. %o_in=o_in*0.3;
  36. %编码器输入
  38. [encoder_out,alphaout]=turbo(o_in);
  39. %编码器输出
  40. in1=encoder_out(1:2,:);
  41. %选取RSC1编码输入作为译码输入
  42. app=zeros(1,length(o_in)+3);
  43. %定义app为全零
  44. in=in1;    
  45. x=in(1,:);              %输入系统位
  46. y=in(2,:);              %输入校验位
  47. in_length=length(in);
  48. % Kw=0;
  50. %---初始化&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
  51. Infty = -128;
  52. d=zeros(8,2,in_length);     %分支量度，8种可能状态，输入为0或者1
  53.                             %D(S,i,k)
  54. a=Infty*ones(8,in_length);  %前向分支量度，A(S,k)
  55. a(1,1)=0;                   %寄存器状态由全零开始
  56. b=Infty*ones(8,in_length+1);%后向分支量度，B(S,k)
  57. b(1,in_length+1)=0;           %寄存器状态由全零结束
  59. %---计算度量和LLR&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
  60. for k=1:in_length
  61.     d(1,2,k)=x(k)+y(k)+app(k);
  62.     d(2,2,k)=d(1,2,k);
  63.     d(7,2,k)=d(1,2,k);
  64.     d(8,2,k)=d(1,2,k);
  65.     
  66.     d(3,2,k)=x(k)+app(k);
  67.     d(4,2,k)=d(3,2,k);
  68.     d(5,2,k)=d(3,2,k);
  69.     d(6,2,k)=d(3,2,k);
  71.     d(3,1,k)=y(k);
  72.     d(4,1,k)=d(3,1,k);
  73.     d(5,1,k)=d(3,1,k);
  74.     d(6,1,k)=d(3,1,k);
  75.     
  76.     if k>1
  77.         a(1,k)=max_th(a(1,k-1),(a(2,k-1)+d(2,2,k-1)));
  78.         a(2,k)=max_th((a(4,k-1)+d(4,1,k-1)),(a(3,k-1)+d(3,2,k-1)));
  79.         a(3,k)=max_th((a(5,k-1)+d(5,1,k-1)),(a(6,k-1)+d(6,2,k-1)));
  80.         a(4,k)=max_th(a(8,k-1),(a(7,k-1)+d(7,2,k-1)));
  81.         a(5,k)=max_th(a(2,k-1),(a(1,k-1)+d(1,2,k-1)));
  82.         a(6,k)=max_th((a(3,k-1)+d(3,1,k-1)),(a(4,k-1)+d(4,2,k-1)));
  83.         a(7,k)=max_th((a(6,k-1)+d(6,1,k-1)),(a(5,k-1)+d(5,2,k-1)));
  84.         a(8,k)=max_th(a(7,k-1),(a(8,k-1)+d(8,2,k-1)));      
  85.     end
  86.     
  87.     if k==in_length
  88.         b(1,k)=max_th(b(1,k+1),(b(5,k+1)+d(1,2,k)));
  89.         b(2,k)=max_th(b(5,k+1),(b(1,k+1)+d(2,2,k)));
  90.         b(3,k)=max_th((b(6,k+1)+d(3,1,k)),(b(2,k+1)+d(3,2,k)));
  91.         b(4,k)=max_th((b(2,k+1)+d(4,1,k)),(b(6,k+1)+d(4,2,k)));
  92.         b(5,k)=max_th((b(3,k+1)+d(5,1,k)),(b(7,k+1)+d(5,2,k)));
  93.         b(6,k)=max_th((b(7,k+1)+d(6,1,k)),(b(3,k+1)+d(6,2,k)));
  94.         b(7,k)=max_th(b(8,k+1),(b(4,k+1)+d(7,2,k)));
  95.         b(8,k)=max_th(b(4,k+1),(b(8,k+1)+d(8,2,k)));
  97.         %计算LLR--------------------------------------
  98.         l(k)=max_th_plus([...
  99.             (a(1,k)+d(1,2,k)+b(5,k+1)),(a(2,k)+d(2,2,k)+b(1,k+1)),...
  100.             (a(3,k)+d(3,2,k)+b(2,k+1)),(a(4,k)+d(4,2,k)+b(6,k+1)),...
  101.             (a(5,k)+d(5,2,k)+b(7,k+1)),(a(6,k)+d(6,2,k)+b(3,k+1)),...
  102.             (a(7,k)+d(7,2,k)+b(4,k+1)),(a(8,k)+d(8,2,k)+b(8,k+1))...
  103.             ])-max_th_plus([...
  104.             (a(1,k)+b(1,k+1)),(a(2,k)+b(5,k+1)),...
  105.             (a(3,k)+d(3,1,k)+b(6,k+1)),(a(4,k)+d(4,1,k)+b(2,k+1)),...
  106.             (a(5,k)+d(5,1,k)+b(3,k+1)),(a(6,k)+d(6,1,k)+b(7,k+1)),...
  107.             (a(7,k)+b(8,k+1)),(a(8,k)+b(4,k+1))...
  108.             ]);
  109.     end
  110. end
  112. for k=in_length-1:-1:1
  113.         b(1,k)=max_th(b(1,k+1),(b(5,k+1)+d(1,2,k)));
  114.         b(2,k)=max_th(b(5,k+1),(b(1,k+1)+d(2,2,k)));
  115.         b(3,k)=max_th((b(6,k+1)+d(3,1,k)),(b(2,k+1)+d(3,2,k)));
  116.         b(4,k)=max_th((b(2,k+1)+d(4,1,k)),(b(6,k+1)+d(4,2,k)));
  117.         b(5,k)=max_th((b(3,k+1)+d(5,1,k)),(b(7,k+1)+d(5,2,k)));
  118.         b(6,k)=max_th((b(7,k+1)+d(6,1,k)),(b(3,k+1)+d(6,2,k)));
  119.         b(7,k)=max_th(b(8,k+1),(b(4,k+1)+d(7,2,k)));
  120.         b(8,k)=max_th(b(4,k+1),(b(8,k+1)+d(8,2,k)));
  122.         %计算LLR--------------------------------------
  123.         l(k)=max_th_plus([...
  124.             (a(1,k)+d(1,2,k)+b(5,k+1)),(a(2,k)+d(2,2,k)+b(1,k+1)),...
  125.             (a(3,k)+d(3,2,k)+b(2,k+1)),(a(4,k)+d(4,2,k)+b(6,k+1)),...
  126.             (a(5,k)+d(5,2,k)+b(7,k+1)),(a(6,k)+d(6,2,k)+b(3,k+1)),...
  127.             (a(7,k)+d(7,2,k)+b(4,k+1)),(a(8,k)+d(8,2,k)+b(8,k+1))...
  128.             ])-max_th_plus([...
  129.             (a(1,k)+b(1,k+1)),(a(2,k)+b(5,k+1)),...
  130.             (a(3,k)+d(3,1,k)+b(6,k+1)),(a(4,k)+d(4,1,k)+b(2,k+1)),...
  131.             (a(5,k)+d(5,1,k)+b(3,k+1)),(a(6,k)+d(6,1,k)+b(7,k+1)),...
  132.             (a(7,k)+b(8,k+1)),(a(8,k)+b(4,k+1))...
  133.             ]);
  134. end
  135. soft_out=l;
  136. ex_info=soft_out-app-x;