obslike.m

来自「CHMMBOX, version 1.2, Iead Rezek, Oxford」· M 代码 · 共 56 行

function [B] = obslike (X,T,n,hmm)

% function [B] = obslike (X,T,n,hmm)
%
% Evaluate likelihood of data given observation model
% for hmm.obsmodel = 'GaussCom','Gauss','AR' or 'LIKE'
% 
% X          N by p data matrix
% T          length of series to learn
% n          block index (time series data can be split into many blocks)
% hmm        hmm data structure
%
% B          Likelihood of N data points

p=length(X(1,:));
N=length(X(:,1));
K=hmm.K;

B=zeros(T,K);
k1=(2*pi)^(-p/2);
switch hmm.obsmodel
  case 'GaussCom',  
    iCov=inv(hmm.state(1).Cov);  % All Covs are the same
    k2=k1/sqrt(det(hmm.state(1).Cov));
    for i=1:T
      for l=1:K
	d=hmm.state(l).Mu-X((n-1)*T+i,:);
	B(i,l)=k2*exp(-0.5*d*iCov*d');
      end
    end
  case 'Gauss',
    for l=1:K  
      state(l).iCov=inv(hmm.state(l).Cov);
      state(l).k2=k1/sqrt(det(hmm.state(l).Cov));
    end
    for i=1:T
      for l=1:K
	d=hmm.state(l).Mu-X((n-1)*T+i,:);
	B(i,l)=state(l).k2*exp(-0.5*d*state(l).iCov*d');
	%if isnan(B(i,l))
	%  keyboard
	%end
      end
    end
  case 'LIKE',
    % The observations are themselves likelihoods
    for l=1:K
        B(:,l)=X(:,l);
    end    
  otherwise
    disp('Unknown observation model');
end