>> tic; x4_mc; toc

trials =

        1000

using N-FIR_length bits to compute BER
GM update using alpha and avg_h and h_1
using bad EM init
uncorrelated h's
GM h

alpha =

        0.995

Using seed from gm.04

seed =

      210675


FIR_length =

     3


N =

    10


avg_h =

     0
     0
     0


var_h =

         0.01


T =

     1     0     0
     0     1     0
     0     0     1


inv_T =

     1     0     0
     0     1     0
     0     0     1


inv_Rh =

   100     0     0
     0   100     0
     0     0   100


Rh =

         0.01            0            0
            0         0.01            0
            0            0         0.01


SNR_array =

     0     2     4     6     8    10

*wrong_where:

ans =

     0   302   211     0
     0   727   235     0
     0   702   255     0
     0   746   238     0
     0   745   230     0
     0   745   241     0
     0   725   246     0
     0   758   232     0
     0   727   289     0
     0   659   359     0


ans =

           10            0            0      0.67029      0.23657            0

*wrong_where:

ans =

     0   231   134     0
     0   752   174     0
     0   751   192     0
     0   769   197     0
     0   760   200     0
     0   784   182     0
     0   780   198     0
     0   763   212     0
     0   737   257     0
     0   695   305     0


ans =

           10            2            0      0.68957      0.18243            0

*wrong_where:

ans =

     0   177   106     0
     0   756   138     0
     0   770   153     0
     0   796   147     0
     0   809   126     0
     0   792   133     0
     0   807   121     0
     0   793   161     0
     0   737   196     0
     0   708   261     0


ans =

           10            4            0        0.701        0.132            0

*wrong_where:

ans =

     0   146    60     0
     0   766    74     0
     0   808    84     0
     0   781    86     0
     0   825    94     0
     0   796    95     0
     0   785   109     0
     0   786   137     0
     0   750   151     0
     0   743   213     0


ans =

           10            6            0        0.701        0.086            0

*wrong_where:

ans =

     0    78    18     0
     0   776    39     0
     0   809    42     0
     0   818    31     0
     0   792    43     0
     0   822    38     0
     0   809    41     0
     0   797    65     0
     0   758    92     0
     0   760   164     0


ans =

           10            8            0      0.70057        0.036            0

*wrong_where:

ans =

     0    39     7     0
     0   807     7     0
     0   804     3     0
     0   830     7     0
     0   804    10     0
     0   826    17     0
     0   810    20     0
     0   806    47     0
     0   778    51     0
     0   753    94     0


ans =

           10           10            0      0.70286     0.010143            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 =

        2.257
        2.263
        2.276
        2.309
        2.323
        2.336


iter_max_g =

     5
     5
     5
     5
     4
     4

EM local min %

ans =

            0
            0
            0
            0
          0.2
          0.2

ML local min %

ans =

     0
     0
     0
     0
     0
     0

EM_equal_ML %

ans =

     0
     0
     0
     0
     0
     0


elapsed_time =

       1214.1

>> 