% x4_em_init_det - initialize for deterministic 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, same as Gaussian initialization
%
  FIR = avg_h;

  G = [];

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

% run Viterbi and initialize b
%
  vit_C;

  b = get_ex;
