HMM+algorithm 源代码

function [y,ss] = hmmgendisc(T,HMM) % % Generate T outputs of a HMM with a discrete output distribution % % function y = hmmg…

function f = hmmfupdate(y,alpha,beta,HMM) % % Provide an update to the state output distributions for the HMM model % % funct…

function [A,pi] = hmmApiup(y,alpha,beta,f,HMM) % % Update the A and pi probabilities in the HMM using the forward and % backwa…

function f = hmmgausfup(y,alpha,beta,HMM) % % Update the Gaussian output distribution f of the HMM using the % probabilities a…

function lpy = hmmlpyseq(y,HMM) % % Find the log likelihood of the sequence y[1],y[2],...,y[T], % i.e., compute log P(y|HMM) …

function f = hmmdiscfup(y,alpha,beta,HMM) % % Update the output probability distribution f of the HMM using the forward % and …

function [alpha,beta,f] = hmmab(y,HMM) % % Compute the forward and backward probabilities for the model HMM % and the output p…

function hmmo = hmmupdate(y,HMM) % % Compute updated HMM model from observations (nonnormalized) % % function HMM = hmmupdate…

% Test data for the HMM % Copyright 1999 by Todd K. Moon HMM.A = [0 1/4 .2; 1 1/4 .4; 0 1/2 .4]; HMM.pi = [1;0;0]; HMM.f…

% Test data for the HMM % Copyright 1999 by Todd K. Moon HMM.A = [0 1/4 .2; 1 1/4 .4; 0 1/2 .4]; HMM.pi = [1;0;0]; HMM.f…