fix hmm demo

bakisit · Mar 10, 2017 · 4db85d6 · 4db85d6
1 parent 9e70e29
commit 4db85d6
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/chapter13/HMM/hmmEm.m b/chapter13/HMM/hmmEm.m
@@ -32,7 +32,7 @@
 %     M-step 
     s = gamma(:,1);                                                                             % 13.18
     A = normalize(A.*(alpha(:,1:n-1)*(beta(:,2:n).*M(:,2:n)./c(2:n))'),2);      % 13.19 13.43 13.65
-    E = bsxfun(@times,gamma*X',1./sum(gamma,2));                                                 % 13.23
+    E = (gamma*X')./sum(gamma,2);                            % 13.23
 end
 model.s = s;
 model.A = A;

diff --git a/demo/ch13/hmm_demo.m b/demo/ch13/hmm_demo.m
@@ -3,12 +3,13 @@
 d = 3;
 k = 2;
 n = 10000;
+[x,model] = hmmRnd(d,k,n);
 %% Viterbi algorithm
 [z, llh] = hmmViterbi(model, x);
 %% HMM filter (forward algorithm)
 [alpha, llh] = hmmFilter(model, x);
 %% HMM smoother (forward backward)
 [gamma,alpha,beta,c] = hmmSmoother(model, x);
 %% Baum-Welch algorithm
-[model, llh] = hmmEm(x,init);
+[model, llh] = hmmEm(x,2);
 plot(llh)