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

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

传真(Fax)编程

开发平台：

Matlab

com_decoder.m：源码内容
  1. function [so,e_p] = com_decoder(a_p,in)
  2. %****************************************************************
  3. % 内容概述：子解码器,输入a_p是先验信息，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. %****************************************************************
  18. x=in(1,:);   %输入系统位 
  19. y=in(2,:);   %输入校验位
  20. %---初始化&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
  21. Infty = -128;
  22. pii=0.5;
  23. % [n,K] = size(g); 
  24. % m = K - 1;
  25. L_seq = length(in);
  26. d(1:2,1:L_seq)=zeros(2,L_seq);  %分支量度，2种可能结果，输入为-1或者1
  27.                                 %D(i,k)
  28. a(1:8,1:L_seq)=Infty*ones(8,L_seq);     %前向分支量度，A(S,k)
  29. a(1,1)=0;                       %寄存器状态由全零开始
  30. b(1:8,1:L_seq+1)=Infty*ones(8,L_seq+1);     %后向分支量度，B(S,k)
  31. b(1,L_seq+1)=0;                       %寄存器状态由全零结束
  32. va1(1:8,1)=zeros(8,1);
  33. va0(1:8,1)=zeros(8,1);
  34. va0_index=[2 1 4 3 6 5 8 7]';
  35. a2a_index=([4 0 1 5 6 2 3 7]+1)';
  36. vb1(1:8,1)=zeros(8,1);
  37. vb0(1:8,1)=zeros(8,1);
  38. vb0_index=[5 6 7 8 1 2 3 4]';
  39. b2b_index=([1 2 5 6 0 3 4 7]+1)';
  41. %初始化结束，开始计算
  42. for k=1:L_seq
  43.     d(1,k)=pii*(a_p(k)+x(k)+y(k));
  44.     d(2,k)=pii*(a_p(k)+x(k)-y(k));
  45. %     d(1,k)=pii*(x(k)+y(k))-log(1+exp(a_p(k)))+a_p(k);
  46. %     d(2,k)=pii*(x(k)-y(k))-log(1+exp(a_p(k)));
  47.     if k>1
  48.         va1=a(:,k-1)+[d(1,k);d(1,k);d(2,k);d(2,k);d(2,k);d(2,k);d(1,k);d(1,k)];
  49.         va0=a(:,k-1)-[d(1,k);d(1,k);d(2,k);d(2,k);d(2,k);d(2,k);d(1,k);d(1,k)];
  50.         va0=va0(va0_index);
  51.         a(:,k)=max(va0,va1);
  52.         a(a2a_index,k)=a(:,k);
  53.     end
  54.     
  55.     if k==L_seq
  56.         %计算后向分支度量
  57.         vb1=b(:,k+1)+[d(1,k);d(2,k);d(2,k);d(1,k);d(1,k);d(2,k);d(2,k);d(1,k)];
  58.         vb0=b(:,k+1)-[d(1,k);d(2,k);d(2,k);d(1,k);d(1,k);d(2,k);d(2,k);d(1,k)];
  59.         vb0=vb0(vb0_index);
  60.         b(:,k)=max(vb0,vb1);
  61.         b(b2b_index,k)=b(:,k);
  62.         %计算LLR；
  63.         llr_a=a(b2b_index,k);
  64.         llr(k)=max(llr_a+vb1)-max(llr_a+vb0);
  65.     end
  66. end
  68. for k=L_seq-1:-1:1
  69.     vb1=b(:,k+1)+[d(1,k);d(2,k);d(2,k);d(1,k);d(1,k);d(2,k);d(2,k);d(1,k)];
  70.     vb0=b(:,k+1)-[d(1,k);d(2,k);d(2,k);d(1,k);d(1,k);d(2,k);d(2,k);d(1,k)];
  71.     vb0=vb0(vb0_index);
  72.     b(:,k)=max(vb0,vb1);
  73.     b(b2b_index,k)=b(:,k);
  74.     %计算LLR；
  75.     llr_a=a(b2b_index,k);
  76.     llr(k)=max(llr_a+vb1)-max(llr_a+vb0);
  77. end
  78. so=llr;
  79. % 软输出
  80. e_p=so-a_p-x;
  81. % 为下一个子解码器提供的外部信息