一类的HMM:
prior0_L = normalise(rand(Q_L,1));
transmat0_L = mk_stochastic(rand(Q_L,Q_L));
[mu0_L, Sigma0_L] = mixgauss_init(Q_L*M_L, reshape(E_L(:,:,1:50), [O T*nex]), cov_type);
mu0_L = reshape(mu0_L, [O Q_L M_L]);
Sigma0_L = reshape(Sigma0_L, [O O Q_L M_L]);
mixmat0_L = mk_stochastic(rand(Q_L,M_L));
[LL_L, prior1_L, transmat1_L, mu1_L, Sigma1_L, mixmat1_L] = ...
mhmm_em(E_L(:,:,1:50), prior0_L, transmat0_L, mu0_L, Sigma0_L, mixmat0_L, 'max_iter', 5);