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

传真(Fax)编程

开发平台:

Matlab

  1. %****************************************************************
  2. % 内容概述:生成随机的输入信号
  3. % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
  4. % 单    位:南京邮电大学,通信工程系
  5. % 创建时间:2005年7月21日
  6. % 修改时间:
  7. % 参考文献:
  8. % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
  9. %****************************************************************
  10. clear;
  11. L_seq=100;
  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. %[hard_out,soft_out]=deturbo(soft_in,alphain);
  24. [hard_out,soft_out]=deturbo_new(soft_in,alphain);
  25. a_p(1:L_seq)=0;
  26. soft_in1=soft_in(1:2,:);
  27. [so,e_p] = com_decoder_new(a_p,soft_in,1)
  28. for i=1:L_seq
  29.     if so(i)>=0
  30.         hard_out(i)=1;
  31.     else
  32.         hard_out(i)=-1;
  33.     end
  34. end 
  35. counter_ber=0;
  36. for jj=1:L_seq
  37.     if soft_in(1,jj)~=hard_out(jj)
  38.         counter_ber=counter_ber+1;
  39.     end
  40. end
  41. counter_ber
  42. soft_in_out=soft_in(1,:);
  43. soft_in_out(1:10)
  44. hard_out(1:10)
  45. soft_out(1:10)