% x4_em_init_gauss - initialize for Gaussian EM algorithm
%
% use initial FIR = [ 0 ... 0 1 0 ... 0]
%
%  FIR = zeros(size(FIR));
%  FIR( (size(FIR,1)+1)/2) = 1;
%
% try initial FIR close to the true value
%
%  FIR = FIR_save .* (1 + 0.1*(2*(rand(size(FIR_save)) > 0.5)-1));
%
% try initial FIR equal to the true value
%
%  FIR = FIR_save;
%
% use known mean of h
%
  FIR = avg_h;

% set FSM state transition cost matrix
%
  trans = FSM_trans( states, FIR);

% initial estimate of G
%
% for Gaussian h with initial 0's in the channel
%
%  G = inv( diag( (N-[0:Memory]) ./ var_noise) + inv_Rh);
%
% for Gaussian h with (-1,1) initial symbols
%
%  G = inv( diag( (N*ones(1,FIR_length)) ./ var_noise) + inv_Rh);
%
% just using Rh
%
  G = inv(inv_Rh);

% run Viterbi and initialize b
%
  vit_C;

  b = get_ex;
