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

传真(Fax)编程

开发平台:

Matlab

  1. %****************************************************************
  2. % 内容概述:生成随机的输入信号
  3. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  4. % 单    位:南京邮电大学,通信工程系
  5. % 创建时间:2005年7月23日
  6. % 修改时间:
  7. % 参考文献:
  8. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  9. %****************************************************************
  10. clear;
  11. L_seq=1000;
  12. random_in=random('Normal',0,1,1,L_seq-3);
  13. for ii=1:L_seq-3
  14.     if random_in(ii)>=0
  15.         random_in(ii)=1;
  16.     else
  17.         random_in(ii)=0;
  18.     end
  19. end
  20. %random_in=ones(1,L_seq-3);
  21. %random_in=random_in*0.8;
  22. [soft_in,alphain]=turbo(random_in);
  23. soft_in2(1,:)=soft_in(1,alphain);
  24. %soft_in2(2,:)=soft_in(3,:);
  25. %soft_out=decoder_3GPP_MAX_new(soft_in2);
  26. %soft_out(alphain)=soft_out;
  27. soft_out=decoder_3GPP_MAX_new(soft_in);
  28. %[hard_out,soft_out]=deturbo_zhang(soft_in,alphain);
  29. for ii=1:L_seq
  30.     if soft_out(ii)>0
  31.         hard_out(ii)=1;
  32.     else
  33.         hard_out(ii)=-1;
  34.     end
  35. end
  36. counter_ber=0;
  37. for jj=1:L_seq
  38.     if soft_in(1,jj)~=hard_out(jj)
  39.         counter_ber=counter_ber+1;
  40.     end
  41. end
  42. counter_ber
  43. soft_in_out=soft_in(1,:);
  44. soft_in_out(1:10);
  45. hard_out(1:10);
  46. soft_out(1:10)