>> 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


trial_count =

   100


trial_count =

   200


trial_count =

   300


trial_count =

   400


trial_count =

   500


trial_count =

   600


trial_count =

   700


trial_count =

   800


trial_count =

   900


trial_count =

        1000

*wrong_where:

ans =

     0   231   223     0
     0   229   218     0
     0   257   244     0
     0   222   212     0
     0   214   216     0
     0   224   216     0
     0   233   216     0
     0   203   193     0
     0   287   271     0
     0   350   354     0


ans =

           10            0            0         0.23      0.22071            0


trial_count =

   100


trial_count =

   200


trial_count =

   300


trial_count =

   400


trial_count =

   500


trial_count =

   600


trial_count =

   700


trial_count =

   800


trial_count =

   900


trial_count =

        1000

*wrong_where:

ans =

     0   172   171     0
     0   179   177     0
     0   170   164     0
     0   166   160     0
     0   164   157     0
     0   162   163     0
     0   169   155     0
     0   171   159     0
     0   233   221     0
     0   306   306     0


ans =

           10            2            0      0.16886      0.16386            0


trial_count =

   100


trial_count =

   200


trial_count =

   300


trial_count =

   400


trial_count =

   500


trial_count =

   600


trial_count =

   700


trial_count =

   800


trial_count =

   900


trial_count =

        1000

*wrong_where:

ans =

     0   107   103     0
     0   122   109     0
     0   136   131     0
     0   108   106     0
     0   108    98     0
     0   104   106     0
     0   111   111     0
     0   134   127     0
     0   182   189     0
     0   295   299     0


ans =

           10            4            0      0.11371      0.10914            0


trial_count =

   100


trial_count =

   200


trial_count =

   300


trial_count =

   400


trial_count =

   500


trial_count =

   600


trial_count =

   700


trial_count =

   800


trial_count =

   900


trial_count =

        1000

*wrong_where:

ans =

     0    75    73     0
     0    85    84     0
     0    69    72     0
     0    67    63     0
     0    77    73     0
     0    76    75     0
     0    82    78     0
     0    88    89     0
     0   142   144     0
     0   227   217     0


ans =

           10            6            0     0.075857        0.074            0


trial_count =

   100


trial_count =

   200


trial_count =

   300


trial_count =

   400


trial_count =

   500


trial_count =

   600


trial_count =

   700


trial_count =

   800


trial_count =

   900


trial_count =

        1000

*wrong_where:

ans =

     0    20    20     0
     0    37    36     0
     0    34    34     0
     0    32    32     0
     0    31    31     0
     0    36    34     0
     0    28    27     0
     0    44    42     0
     0    87    84     0
     0   191   173     0


ans =

           10            8            0     0.031143     0.030571            0


trial_count =

   100


trial_count =

   200


trial_count =

   300


trial_count =

   400


trial_count =

   500


trial_count =

   600


trial_count =

   700


trial_count =

   800


trial_count =

   900


trial_count =

        1000

*wrong_where:

ans =

     0    12    12     0
     0    12    12     0
     0     9     9     0
     0     9     9     0
     0    10    10     0
     0    14    13     0
     0    18    17     0
     0    26    26     0
     0    55    54     0
     0   137   131     0


ans =

           10           10            0        0.012     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.361
        1.283
        1.184
         1.09
        1.039
        1.026


iter_max_g =

     4
     3
     3
     3
     3
     2

EM local min %

ans =

         65.7
         62.6
         58.3
         62.2
         62.6
         70.1

ML local min %

ans =

     0
     0
     0
     0
     0
     0

EM_equal_ML %

ans =

     0
     0
     0
     0
     0
     0


elapsed_time =

       993.56

>> 