>> tic; x4_mc; toc

trials =

        1000

using N-FIR_length bits to compute BER
uncorrelated h's
GM h
Gauss EM and ML

alpha =

          0.8

Using seed from gm.04

seed =

      210675


FIR_length =

     3


N =

    10


avg_h =

          0.2
          0.2
          0.2


var_h =

          0.1


T =

     1     0     0
     1     1     0
     1     1     1


inv_T =

     1     0     0
    -1     1     0
     0    -1     1


inv_Rh =

    30    20    10
    20    20    10
    10    10    10


Rh =

          0.1         -0.1            0
         -0.1          0.2         -0.1
            0         -0.1          0.2


SNR_array =

     0     2     4     6     8    10

*wrong_where:

ans =

     0   258   223     0
     0   265   218     0
     0   298   244     0
     0   273   212     0
     0   245   216     0
     0   269   216     0
     0   272   216     0
     0   257   193     0
     0   308   271     0
     0   347   354     0


ans =

           10            0            0      0.26857      0.22071            0

*wrong_where:

ans =

     0   200   171     0
     0   253   177     0
     0   247   164     0
     0   239   160     0
     0   236   157     0
     0   232   163     0
     0   253   155     0
     0   271   159     0
     0   282   221     0
     0   320   306     0


ans =

           10            2            0      0.23714      0.16386            0

*wrong_where:

ans =

     0   170   103     0
     0   196   109     0
     0   234   131     0
     0   242   106     0
     0   201    98     0
     0   222   106     0
     0   201   111     0
     0   239   127     0
     0   259   189     0
     0   320   299     0


ans =

           10            4            0      0.20943      0.10914            0

*wrong_where:

ans =

     0   137    73     0
     0   165    84     0
     0   194    72     0
     0   186    63     0
     0   198    73     0
     0   196    75     0
     0   219    78     0
     0   218    89     0
     0   259   144     0
     0   265   217     0


ans =

           10            6            0        0.185        0.074            0

*wrong_where:

ans =

     0    78    20     0
     0   136    36     0
     0   164    34     0
     0   175    32     0
     0   194    31     0
     0   189    34     0
     0   190    27     0
     0   210    42     0
     0   227    84     0
     0   298   173     0


ans =

           10            8            0      0.16086     0.030571            0

*wrong_where:

ans =

     0    74    12     0
     0   121    12     0
     0   151     9     0
     0   161     9     0
     0   177    10     0
     0   175    13     0
     0   187    17     0
     0   200    26     0
     0   233    54     0
     0   271   131     0


ans =

           10           10            0      0.14943     0.011714            0

N, SNR, det, g, known, ML

iter_avg_det =

     0
     0
     0
     0
     0
     0


iter_max_det =

     0
     0
     0
     0
     0
     0


iter_avg_g =

        1.055
        1.037
        1.018
        1.008
        1.004
        1.002


iter_max_g =

     3
     3
     2
     2
     2
     2

EM local min %

ans =

         75.5
         71.3
         70.4
           72
         66.8
         69.7

ML local min %

ans =

     0
     0
     0
     0
     0
     0

EM_equal_ML %

ans =

     0
     0
     0
     0
     0
     0


elapsed_time =

       1257.7

>>