turbo_logmap.txt

常用的编码 差错控制能力很强大的turbo码 用matlab实现的编码。

 Please enter the decoding algorithm. (0:Log-MAP, 1:SOVA)  default 0    
 Please enter the frame size (= info + tail, default: 400)   
 Please enter code generator: ( default: g = [1 1 1; 1 0 1 ] )      
 Please choose punctured / unpunctured (0/1): default 0     
 Please enter number of iterations for each frame: default 5       
 Please enter number of frame errors to terminate: default 15        
 Please enter Eb/N0 in dB : default [2.0]    4


----------------------------------------------------
 === Log-MAP decoder === 
 Frame size =    400
 code generator: 
     1     1     1
     1     0     1
 Punctured, code rate = 1/2 
 iteration number =       5
 terminate frame errors =     15
 Eb / N0 (dB) =       4.00
----------------------------------------------------

+ + + + Please be patient. Wait a while to get the result. + + + +
openvar('r', r);
dbquit
fschange('D:\research\multi_user_dete\turbo_sys_demo.m');
clear turbo_sys_demo
turbo_sys_demo
 Please enter the decoding algorithm. (0:Log-MAP, 1:SOVA)  default 0    
 Please enter the frame size (= info + tail, default: 400)   
 Please enter code generator: ( default: g = [1 1 1; 1 0 1 ] )      
 Please choose punctured / unpunctured (0/1): default 0     
 Please enter number of iterations for each frame: default 5       
 Please enter number of frame errors to terminate: default 15        
 Please enter Eb/N0 in dB : default [2.0]    7.6464


----------------------------------------------------
 === Log-MAP decoder === 
 Frame size =    400
 code generator: 
     1     1     1
     1     0     1
 Punctured, code rate = 1/2 
 iteration number =       5
 terminate frame errors =     15
 Eb / N0 (dB) =       7.65
----------------------------------------------------

+ + + + Please be patient. Wait a while to get the result. + + + +
openvar('sigma', sigma);
openvar('sigma', sigma);
erfc(1/sigma)

ans =

  6.4812e-004

erfc(1/(sigma*1.414))

ans =

    0.0159

erfc(1/(sigma*sqrt(2.0)))

ans =

    0.0159

turbo_sys_demo
 Please enter the decoding algorithm. (0:Log-MAP, 1:SOVA)  default 0    
 Please enter the frame size (= info + tail, default: 400)   
 Please enter code generator: ( default: g = [1 1 1; 1 0 1 ] )      
 Please choose punctured / unpunctured (0/1): default 0     
 Please enter number of iterations for each frame: default 5       
 Please enter number of frame errors to terminate: default 15        
 Please enter Eb/N0 in dB : default [2.0]    5.5383


----------------------------------------------------
 === Log-MAP decoder === 
 Frame size =    400
 code generator: 
     1     1     1
     1     0     1
 Punctured, code rate = 1/2 
 iteration number =       5
 terminate frame errors =     15
 Eb / N0 (dB) =       5.54
----------------------------------------------------

+ + + + Please be patient. Wait a while to get the result. + + + +
Error in ==> D:\research\turbo\logmapo.m
On line 57  ==>          Alpha(k,state2) = log( sum( exp( gamma+Alpha(k-1,:) ) ) );  

Error in ==> D:\research\turbo\turbo_sys_demo.m
On line 145  ==>             L_all = logmapo(rec_s(2,:), g, L_a, 2);  % complete info.  

erfc(1/(sigma*sqrt(2.0)))

ans =

    0.0585

   en = 10^(EbN0db(nEN)/10);      % convert Eb/N0 from unit db to normal numbers
   L_c = 4*a*en*rate; 	% reliability value of the channel
   sigma = 1/sqrt(2*rate*en); 	% standard deviation of AWGN noise
sigma = 1/sqrt(2*10^(7.0343/10));
erfc(1/sigma)

ans =

  6.9518e-006

openvar('sigma', sigma);
sigma = 1/sqrt(2*0.4*10^(7.0343/10));
erfc(1/sigma)

ans =

    0.0045

sigma = 1/sqrt(2*10^(7.6464/10));
erfc(1/sigma)

ans =

  1.4116e-006

sigma = 1/sqrt(2*0.4*10^(7.6464/10));
erfc(1/sigma)

ans =

    0.0023

sigma = 1/sqrt(2*10^(7.6464/10));
erfc(1/sigma)

ans =

  1.4116e-006

en= 10^(7.6464/10)

en =

    5.8162

sigma = 1/sqrt(2*en);
erfc(1/sigma)

ans =

  1.4116e-006

erfc(1/0.9114)

ans =

    0.1207

erfc(1/0.9194)

ans =

    0.1240

erfc(1/0.9194/0.9194)

ans =

    0.0943

erfc(1/0.9194/1.414)

ans =

    0.2767

0.2767/2

ans =

    0.1384

 sigma = 1/sqrt(2*10^(7.6464/10));
0.5*erfc(1/sigma/1.414)

ans =

  3.2345e-004

sigma = 1/sqrt(2*0.4*10^(7.6464/10));
0.5*erfc(1/sigma/1.414)

ans =

    0.0155

sigma = 1/sqrt(2*0.4*10^(5.5383/10));
0.5*erfc(1/sigma/1.414)

ans =

    0.0453

exit
 Please enter the decoding algorithm. (0:Log-MAP, 1:SOVA)  default 0    0
 Please enter the frame size (= info + tail, default: 400)   1024
 Please enter code generator: ( default: g = [1 1 1; 1 0 1 ] )      [1 1 1 1 1;1 0 0 0 1]
 Please choose punctured / unpunctured (0/1): default 0     1
 Please enter number of iterations for each frame: default 5       5
 Please enter number of frame errors to terminate: default 15        10
 Please enter Eb/N0 in dB : default [2.0]    0


----------------------------------------------------
 === Log-MAP decoder === 
 Frame size =   1024
 code generator: 
     1     1     1     1     1
     1     0     0     0     1
 Unpunctured, code rate = 1/3 
 iteration number =       5
 terminate frame errors =     10
 Eb / N0 (dB) =       0.00
----------------------------------------------------

+ + + + Please be patient. Wait a while to get the result. + + + +
************** Eb/N0 =  0.00 db **************
Frame size = 1024, rate 1/3. 
3 frames transmitted, 3 frames in error.
Bit Error Rate (from iteration 1 to iteration 5):
1.3105e-001    1.0556e-001    8.7908e-002    7.8758e-002    6.6667e-002    
Frame Error Rate (from iteration 1 to iteration 5):
1.0000e+000    1.0000e+000    1.0000e+000    1.0000e+000    1.0000e+000    
***********************************************

************** Eb/N0 =  0.00 db **************
Frame size = 1024, rate 1/3. 
6 frames transmitted, 6 frames in error.
Bit Error Rate (from iteration 1 to iteration 5):
1.3709e-001    1.0752e-001    8.4314e-002    7.5327e-002    6.6340e-002    
Frame Error Rate (from iteration 1 to iteration 5):
1.0000e+000    1.0000e+000    1.0000e+000    1.0000e+000    1.0000e+000    
***********************************************

************** Eb/N0 =  0.00 db **************
Frame size = 1024, rate 1/3. 
9 frames transmitted, 9 frames in error.
Bit Error Rate (from iteration 1 to iteration 5):
1.3780e-001    1.1078e-001    8.8780e-002    7.6362e-002    6.4052e-002    
Frame Error Rate (from iteration 1 to iteration 5):
1.0000e+000    1.0000e+000    1.0000e+000    1.0000e+000    1.0000e+000    
***********************************************

************** Eb/N0 =  0.00 db **************
Frame size = 1024, rate 1/3. 
10 frames transmitted, 10 frames in error.
Bit Error Rate (from iteration 1 to iteration 5):
1.4353e-001    1.1775e-001    9.8529e-002    8.7549e-002    7.5980e-002    
Frame Error Rate (from iteration 1 to iteration 5):
1.0000e+000    1.0000e+000    1.0000e+000    1.0000e+000    1.0000e+000    
***********************************************

 Please enter the decoding algorithm. (0:Log-MAP, 1:SOVA)  default 0    0
 Please enter the frame size (= info + tail, default: 400)   1024
 Please enter code generator: ( default: g = [1 1 1; 1 0 1 ] )      [1 1 1 1 1;1 0 0 0 1]
 Please choose punctured / unpunctured (0/1): default 0     1
 Please enter number of iterations for each frame: default 5       5
 Please enter number of frame errors to terminate: default 15        10
 Please enter Eb/N0 in dB : default [2.0]    0.5


----------------------------------------------------
 === Log-MAP decoder === 
 Frame size =   1024
 code generator: 
     1     1     1     1     1
     1     0     0     0     1
 Unpunctured, code rate = 1/3 
 iteration number =       5
 terminate frame errors =     10
 Eb / N0 (dB) =       0.50
----------------------------------------------------

+ + + + Please be patient. Wait a while to get the result. + + + +
************** Eb/N0 =  0.50 db **************
Frame size = 1024, rate 1/3. 
3 frames transmitted, 1 frames in error.
Bit Error Rate (from iteration 1 to iteration 5):
1.0098e-001    2.9739e-002    5.2288e-003    6.5359e-004    6.5359e-004    
Frame Error Rate (from iteration 1 to iteration 5):
1.0000e+000    1.0000e+000    1.0000e+000    3.3333e-001    3.3333e-001    
***********************************************

************** Eb/N0 =  0.50 db **************
Frame size = 1024, rate 1/3. 
6 frames transmitted, 1 frames in error.
Bit Error Rate (from iteration 1 to iteration 5):
1.0000e-001    2.9739e-002    5.0654e-003    3.2680e-004    3.2680e-004    
Frame Error Rate (from iteration 1 to iteration 5):
1.0000e+000    1.0000e+000    1.0000e+000    1.6667e-001    1.6667e-001    
***********************************************

Error in ==> D:\turbo\turbo\logmapo.m
On line 72  ==>      if(sum(exp(gamma+Beta(k+1,:)))<1e-300)

Error in ==> D:\turbo\turbo\turbo_sys_demo.m
On line 136  ==>             L_all = logmapo(rec_s(1,:), g, L_a, 1);  % complete info.

turbo_sys_demo
 Please enter the decoding algorithm. (0:Log-MAP, 1:SOVA)  default 0    0
 Please enter the frame size (= info + tail, default: 400)   1024
 Please enter code generator: ( default: g = [1 1 1; 1 0 1 ] )      [1 1 1 1 1;1 0 0 0 1]
 Please choose punctured / unpunctured (0/1): default 0     1
 Please enter number of iterations for each frame: default 5       2
 Please enter number of frame errors to terminate: default 15        1
 Please enter Eb/N0 in dB : default [2.0]    1


----------------------------------------------------
 === Log-MAP decoder === 
 Frame size =   1024
 code generator: 
     1     1     1     1     1
     1     0     0     0     1
 Unpunctured, code rate = 1/3 
 iteration number =       2
 terminate frame errors =      1
 Eb / N0 (dB) =       1.00
----------------------------------------------------

+ + + + Please be patient. Wait a while to get the result. + + + +
************** Eb/N0 =  1.00 db **************
Frame size = 1024, rate 1/3. 
3 frames transmitted, 0 frames in error.
Bit Error Rate (from iteration 1 to iteration 2):
4.2157e-002    0.0000e+000    
Frame Error Rate (from iteration 1 to iteration 2):
1.0000e+000    0.0000e+000    
***********************************************

************** Eb/N0 =  1.00 db **************
Frame size = 1024, rate 1/3. 
4 frames transmitted, 1 frames in error.
Bit Error Rate (from iteration 1 to iteration 2):
4.4608e-002    1.7157e-003    
Frame Error Rate (from iteration 1 to iteration 2):
1.0000e+000    2.5000e-001    
***********************************************