test_logmap.m
上传用户: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年9月5日
  10. % 修改时间:
  11. % 参考文献:
  12. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  13. %****************************************************************
  14. clear;
  15. clc;
  16. time_begin=datestr(now);
  17. rate=1/3;           %码率
  18. m=3;                    %尾比特数
  19. fading_a=1;             %Fading amplitude
  20. snr=0:0.2:2.4;                            %snr的采样点
  21. EbNoLinear=10.^(snr.*0.1);
  22. iter=[1 2 3];                                %迭代次数
  23. ferrlim=10;                              %误帧限,达到此限即可停止当前SNR点的测试
  24. %num_frame=round(10.^(snr+1));           %测试的帧数;
  25. num_block_size=1024;                    %测试的块尺寸,指包含尾比特的软输入系统系列长度
  26. err_counter=zeros(max(iter),length(snr));        %初始化错误比特计数器
  27. nferr= zeros(max(iter),length(snr));             %初始化错误帧计数器
  28. ber=zeros(max(iter),length(snr));                 %初始化错误比特率
  29. random_in=round(rand(1,num_block_size-m));  %随机数
  30. [turbod_out,alphain]=turbo(random_in);      %编码
  31. for ii=1:length(iter)
  32.     for nEN=1:length(snr)
  33.         L_c=4*fading_a*EbNoLinear(nEN)*rate;
  34.         sigma=1/sqrt(2*rate*EbNoLinear(nEN));
  35.         nframe = 0;    % clear counter of transmitted frames
  36.         if nEN==1 | ber(iter(ii),nEN-1)>9.0e-6
  37.             while nferr(iter(ii),nEN)<ferrlim        %nferr:当前迭代次数、snr点的错误帧数
  38.                 nframe = nframe + 1; 
  39.                 noice=randn(3,num_block_size);    %噪声
  40.                 soft_in=L_c*(turbod_out+sigma*noice);            %信息噪声叠加
  41.                 [hard_out,soft_out]=deturbo_logmap(soft_in,alphain,iter(ii)); %译码
  42.                 errs=length(find(hard_out(1:num_block_size-m)~=random_in));%当前点错误bit数
  43.                 if errs>0 
  44.                     err_counter(iter(ii),nEN)=err_counter(iter(ii),nEN)+errs;
  45.                     nferr(iter(ii),nEN)=nferr(iter(ii),nEN)+1;
  46.                 end
  47.             end
  48.             ber(iter(ii),nEN) = err_counter(iter(ii),nEN)/nframe/(num_block_size-m);%误比特率
  49.             %fer(iter,nEN) = nferr(iter,nEN)/nframe; %误帧率
  50.         else
  51.             ber(iter(ii),nEN)=1.0e-7;
  52.         end
  53.         fprintf('迭代次数:%1.0f;snr:%1.2f;误码率:%8.4e;n',...
  54.             iter(ii),snr(nEN),ber(iter(ii),nEN));
  55.         %save cap算法06_准Th_log_map.mat snr ber;
  56.     end
  57. end
  58. semilogy(snr,ber(1,:),snr,ber(2,:),snr,ber(3,:));
  59. xlabel('SNR(dB)');
  60. ylabel('Bit Error Rate');
  61. title('3GPP标准 准Threshold-log-map译码算法 译码性能图,1024交织长度,WYF噪声加法');
  62. legend('1次迭代','2次迭代','3次迭代');
  63. time_end=datestr(now);
  64. disp(time_begin);
  65. disp(time_end);