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

传真(Fax)编程

开发平台:

Matlab

  1. %****************************************************************
  2. % 内容概述:TURBO译码计算速对测试(有菜单)
  3. %          设定计算时长,检测单位时间内的译码bit数
  4. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  5. % 单    位:南京邮电大学,通信工程系
  6. % 创建时间:2005年10月29日
  7. % 修改时间:
  8. % 参考文献:
  9. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  10. %****************************************************************
  11. clc;
  12. clear all;
  13. timerLimit = input('测试时长【60秒】:');
  14. if isempty(timerLimit)
  15.    timerLimit = 60;
  16. end
  17. length_interleave = input('交织长度=帧长-尾比特长度【1024】:');
  18. if isempty(length_interleave)
  19.    length_interleave = 1024;
  20. end
  21. iter = input('迭代次数【10】:');
  22. if isempty(iter)
  23.    iter =10;
  24. end
  25. algorithm = input('译码算法【1:LOG-MAP(缺省),2:MAX-LOG-MAP,3:SEMITH-LOG-MAP)】:');
  26. if isempty(algorithm)
  27.    algorithm =1;
  28. end
  29. save_mat = input('是否保存测试结果到MAT文件 【1-保存,0-不保存(缺省)】:');
  30. if isempty(save_mat)
  31.    save_mat=0;
  32. end
  33. if save_mat==1
  34.     matFileName = input('MAT文件名 【''运算时长临时测试数据.mat''】:');
  35.     if isempty(matFileName)
  36.         matFileName='运算时长临时测试数据.mat';
  37.     end
  38. end
  39. fprintf('----------------------测试参数----------------------n'); 
  40. fprintf(' 交织长度=%4dbit;迭代次数=%2dn',length_interleave,iter);
  41. fprintf(' 运算时长=%4d秒n',timerLimit);
  42. switch algorithm
  43.     case 1
  44.         fprintf(' 译码算法:LOG-MAPn');
  45.     case 2
  46.         fprintf(' 译码算法:MAX-LOG-MAPn');
  47.     case 3
  48.         fprintf(' 译码算法:门限MAX-LOG-MAPn');
  49. end
  50. if save_mat==1
  51.     fprintf(' 保存测试结果到 : %4sn',matFileName);
  52. else
  53.     fprintf(' 不保存测试结果到文件n');
  54. end    
  55. fprintf('----------------------------------------------------n'); 
  56. rate=1/3;           %码率
  57. m=3;                    %尾比特数
  58. fading_a=1;             %Fading amplitude
  59. EbNo=1.0;                            %EbNo的采样点
  60. EbNoLinear=10.^(EbNo.*0.1);
  61. num_block_size=length_interleave+m;     %测试的块尺寸,指包含尾比特的软输入系统系列长度
  62. random_in=round(rand(1,length_interleave));  %随机数
  63. [turbod_out,alphain]=turbo(random_in);      %编码
  64. L_c=4*fading_a*EbNoLinear*rate;
  65. sigma=1/sqrt(2*rate*EbNoLinear);
  66. nframe = 0;    % clear counter of transmitted frames
  67. time_begin=clock;
  68. while etime(clock,time_begin)<timerLimit        %nferr:当前迭代次数、EbNo点的错误帧数
  69.     nframe = nframe + 1; 
  70.     noice=randn(4,num_block_size);    %噪声
  71.     soft_in=L_c*(turbod_out+sigma*noice);            %信息噪声叠加
  72.     switch algorithm
  73.         case 1
  74.             [hard_out,soft_out]=decoder_logmap(soft_in,alphain,iter); %log-map译码
  75.         case 2
  76.             [hard_out,soft_out]=decoder(soft_in,alphain,iter); %max-log-map译码
  77.         case 3
  78.             [hard_out,soft_out]=decoder_SemiTh(soft_in,alphain,iter); %门限max-log-map译码
  79.     end
  80. end
  81. if save_mat==1
  82.     save (matFileName,'length_interleave','iter','nframe','timerLimit');
  83. end
  84. fprintf('----------------------测试结果----------------------n');  
  85. fprintf(' 交织长度=%4dbit; 迭代次数=%2dn',length_interleave,iter);
  86. fprintf(' 运算帧数=%2d;运算bit数=%6dn',nframe,nframe*length_interleave);
  87. fprintf(' 运算时长=%4d秒;运算bit速率=%6.2e bit/sn',timerLimit,nframe*length_interleave/timerLimit);
  88. switch algorithm
  89.     case 1
  90.         fprintf(' 译码算法:LOG-MAPn');
  91.     case 2
  92.         fprintf(' 译码算法:MAX-LOG-MAPn');
  93.     case 3
  94.         fprintf(' 译码算法:门限MAX-LOG-MAPn');
  95. end
  96. if save_mat==1
  97.     fprintf(' 保存测试结果到 : %4sn',matFileName);
  98. else
  99.     fprintf(' 未保存测试结果到文件n');
  100. end    
  101. fprintf('----------------------------------------------------n');