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

传真(Fax)编程

开发平台:

Matlab

  1. %****************************************************************
  2. % 内容概述:AWGN信道测试
  3. %          每个SNR采样点的测试帧数是随snr指数上升的,
  4. %          这样不仅减少了计算量而且可以保证精度。
  5. %          每帧的输入是相同的,(不是每帧都不一样,这样可以大幅降低计算量)
  6. %          在帧比特足够多的情况下,应该可以保证随机性(未理论证实)
  7. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  8. % 单    位:南京邮电大学,通信工程系
  9. % 创建时间:2005年7月24日
  10. % 修改时间:
  11. % 参考文献:
  12. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  13. %****************************************************************
  14. clear;
  15. clc;
  16. time_begin=datestr(now);
  17. snr=0::3;                          %snr的采样点
  18. rate=1/3;                               %码率
  19. num_frame=round(3.2.^(snr+1));           %测试的帧数;
  20. num_block_size=1024;                    %测试的块尺寸,指包含尾比特的软输入系统系列长度
  21. err_counter(3,1:length(snr))=0;         %初始化错误比特计数器
  22. guage=0;                                %初始化计算量
  23. counter_guage=0;                        %初始化计算量计数器
  24. for num_it=1:1:3
  25.     for ii=1:length(snr)
  26.         guage=guage+num_frame(ii);      %统计计算量
  27.     end
  28. end
  29. random_in=round(rand(1,num_block_size-3));  %随机数
  30. [turbod_out,alphain]=turbo(random_in);      %编码
  31. for num_it=1:1:3
  32.     for ii=1:length(snr)
  33.         for frame=1:num_frame(ii)
  34.             soft_in=awgn(turbod_out,snr(ii)+10*log10(rate),'measured');    %通过AWGN信道,加入噪声
  35.             counter_guage=counter_guage+1;
  36.             fprintf('总帧数:%6.0f;完成帧数:%6.0f;当前snr点总帧数:%4.0f;当前snr点完成帧数:%4.0f;计算进度:%2.1f%% n',...
  37.                 guage,counter_guage,num_frame(ii),frame,counter_guage*100/guage);            
  38.             [hard_out,soft_out]=deturbo_cap(soft_in,alphain,num_it);
  39.             err_counter(num_it,ii)=err_counter(num_it,ii)+...
  40.             length(find(hard_out(1:num_block_size-3)~=random_in));
  41.         end
  42.         ber(num_it,ii)=err_counter(num_it,ii)/((num_block_size-3)*num_frame(ii));
  43.     end
  44. end
  45. semilogy(snr,ber(1,:),snr,ber(2,:),snr,ber(3,:));
  46. xlabel('SNR(dB)');
  47. ylabel('BER');
  48. title('3GPP标准 max-log-map译码算法 译码性能图');
  49. legend('1次迭代','2次迭代','3次迭代');
  50. %save cap算法06.mat snr ber soft_out;
  51. time_end=datestr(now);
  52. disp(time_begin);
  53. disp(time_end);
  54. %plot(snr,ber);
  55. %counter_ber
  56. %soft_in(1,1:10);
  57. %hard_out(1:10);
  58. %soft_out(1:10);
  59. %--------------------------
  60. %soft_in2(1,:)=soft_in(1,alphain);
  61. %soft_in2(2,:)=soft_in(3,:);
  62. %soft_out=decoder_3GPP_MAX_new(soft_in2);
  63. %soft_out(alphain)=soft_out;
  64. %----------------------------
  65. %soft_out=decoder_3GPP_MAX_new(soft_in);
  66. %for ii=1:num_block_size
  67. %    if soft_out(ii)>0
  68. %        hard_out(ii)=1;
  69. %    else
  70. %        hard_out(ii)=-1;
  71. %    end
  72. %end
  73. %----------------------------