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

传真(Fax)编程

开发平台:

Matlab

  1. %****************************************************************
  2. % 内容概述:TURBO译码性能AWGN信道测试(有菜单)
  3. %          达到误帧限即可停止当前SNR点的测试,节省计算量
  4. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  5. % 单    位:南京邮电大学,通信工程系
  6. % 创建时间:2005年9月12日
  7. % 修改时间:2005年10月26日
  8. % 参考文献:
  9. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  10. %****************************************************************
  11. clc;
  12. clear all;
  13. length_interleave = input('交织长度=帧长-尾比特长度【1024】:');
  14. if isempty(length_interleave)
  15.    length_interleave = 1024;
  16. end
  17. iter = input('迭代次数【[1,2,3]】:');
  18. if isempty(iter)
  19.    iter =[1 2 3];
  20. end
  21. ferrlim = input('误帧限(达到此限即可停止当前SNR点的测试)【10】:');
  22. if isempty(ferrlim)
  23.    ferrlim =10;
  24. end
  25. max_EbNo = input('测试最大Eb/No(dB)【2.4】:');
  26. if isempty(max_EbNo)
  27.    max_EbNo =2.4;
  28. end
  29. step_EbNo = input('测试信噪比步长(dB)【0.2】:');
  30. if isempty(step_EbNo)
  31.    step_EbNo=0.2;
  32. end
  33. save_mat = input('是否保存测试结果到MAT文件 【1-保存(缺省),0-不保存】:');
  34. if isempty(save_mat)
  35.    save_mat=1;
  36. end
  37. if save_mat==1
  38.     matFileName = input('MAT文件名 【''临时测试数据.mat''】:');
  39.     if isempty(matFileName)
  40.         matFileName='临时测试数据.mat';
  41.     end
  42. end
  43. fprintf('----------------------------------------------------n'); 
  44. fprintf(' 交织长度=%6dn',length_interleave);
  45. fprintf(' 测试最大Eb/No = %2.1fdB;测试信噪比步长 = %2.1fdBn',max_EbNo,step_EbNo);
  46. if save_mat==1
  47.     fprintf(' 保存测试结果到 = %4sn',matFileName);
  48. end    
  49. fprintf('----------------------------------------------------n'); 
  50. time_begin=datestr(now);
  51. rate=1/3;           %码率
  52. m=3;                    %尾比特数
  53. fading_a=1;             %Fading amplitude
  54. EbNo=0:step_EbNo:max_EbNo;                            %EbNo的采样点
  55. EbNoLinear=10.^(EbNo.*0.1);
  56. num_block_size=length_interleave+m;     %测试的块尺寸,指包含尾比特的软输入系统系列长度
  57. err_counter=zeros(max(iter),length(EbNo));        %初始化错误比特计数器
  58. nferr= zeros(max(iter),length(EbNo));             %初始化错误帧计数器
  59. ber=zeros(max(iter),length(EbNo));                 %初始化错误比特率
  60. random_in=round(rand(1,length_interleave));  %随机数
  61. [turbod_out,alphain]=turbo(random_in);      %编码
  62. for ii=1:length(iter)
  63.     for nEN=1:length(EbNo)
  64.         L_c=4*fading_a*EbNoLinear(nEN)*rate;
  65.         sigma=1/sqrt(2*rate*EbNoLinear(nEN));
  66.         nframe = 0;    % clear counter of transmitted frames
  67.         if nEN==1 | ber(iter(ii),nEN-1)>9.0e-6
  68.             while nferr(iter(ii),nEN)<ferrlim        %nferr:当前迭代次数、EbNo点的错误帧数
  69.                 nframe = nframe + 1; 
  70.                 noice=randn(4,num_block_size);    %噪声
  71.                 soft_in=L_c*(turbod_out+sigma*noice);            %信息噪声叠加
  72.                 [hard_out,soft_out]=decoder(soft_in,alphain,iter(ii)); %译码
  73.                 errs=length(find(hard_out(1:length_interleave)~=random_in));%当前点错误bit数
  74.                 
  75.                 if errs>0 
  76.                     err_counter(iter(ii),nEN)=err_counter(iter(ii),nEN)+errs;
  77.                     nferr(iter(ii),nEN)=nferr(iter(ii),nEN)+1;
  78.                 end
  79.                 fprintf('当前EbNo点:%1.2fdB;已计算:%2.0f帧;其中:%2.0f误帧n',...
  80.                     EbNo(nEN),nframe,nferr(iter(ii),nEN));
  81.             end
  82.             ber(iter(ii),nEN) = err_counter(iter(ii),nEN)/nframe/(length_interleave);%误比特率
  83.             %fer(iter,nEN) = nferr(iter,nEN)/nframe; %误帧率
  84.         else
  85.             ber(iter(ii),nEN)=1.0e-7;
  86.         end
  87.         fprintf('迭代次数:%1.0f;EbNo:%1.2fdB;误码率:%8.4e;n',...
  88.             iter(ii),EbNo(nEN),ber(iter(ii),nEN));
  89.         if save_mat==1
  90.             save (matFileName,'EbNo','ber');
  91.         end
  92.     end
  93. end
  94. %semilogy(EbNo,ber(1,:),EbNo,ber(2,:),EbNo,ber(3,:));
  95. %xlabel('E_b/N_0 (dB)');
  96. %ylabel('Bit Error Rate');
  97. %title('3GPP标准 max-log-map译码算法,1024交织长度,1/3码率');
  98. %legend('1次迭代','2次迭代','3次迭代');
  99. time_end=datestr(now);
  100. disp(time_begin);
  101. disp(time_end);
  102. disp('恭喜你!测试完成!n');