HMM viterbi

function states = my_dhmm_viterbi(obs, prior, trans, emit)
  % Input
  %  obs   : sequence of observations
  %  prior : initial state probs
  %  trans : transition matrix
  %  emit  : emission matrix
  % Output
  %  states: sequence of states
  %
  % See: http://www.cs.ubc.ca/~murphyk/Software/HMM/hmm_usage.html
  %

  V = multinomial_prob(obs, emit);
  states = viterbi_path(prior, trans, V);
end