Commit dcec34b0 authored by Sylvain Calinon's avatar Sylvain Calinon

Commented alpha param in semi-tied GMM

parent 162b66a9
...@@ -77,6 +77,6 @@ view(-40,6); axis equal; ...@@ -77,6 +77,6 @@ view(-40,6); axis equal;
%print('-dpng','graphs/demo_semitiedGMM01.png'); %print('-dpng','graphs/demo_semitiedGMM01.png');
pause; %pause;
close all; %close all;
...@@ -56,6 +56,10 @@ end ...@@ -56,6 +56,10 @@ end
if ~isfield(model,'params_nbVariationSteps') if ~isfield(model,'params_nbVariationSteps')
model.params_nbVariationSteps = 50; model.params_nbVariationSteps = 50;
end end
% if ~isfield(model,'params_alpha')
% model.params_alpha = 1.0;
% end
if ~isfield(model,'B') if ~isfield(model,'B')
model.B = eye(model.nbVar) * model.params_Bsf; model.B = eye(model.nbVar) * model.params_Bsf;
model.InitH = pinv(model.B) + eye(model.nbVar) * model.params_diagRegFact; model.InitH = pinv(model.B) + eye(model.nbVar) * model.params_diagRegFact;
...@@ -100,6 +104,7 @@ for nbIter=1:model.params_nbMaxSteps ...@@ -100,6 +104,7 @@ for nbIter=1:model.params_nbMaxSteps
model.H = pinv(model.B) + eye(model.nbVar) * model.params_diagRegFact; model.H = pinv(model.B) + eye(model.nbVar) * model.params_diagRegFact;
for i=1:model.nbStates for i=1:model.nbStates
model.Sigma(:,:,i) = model.H * model.SigmaDiag(:,:,i) * model.H'; %Eq.(3) model.Sigma(:,:,i) = model.H * model.SigmaDiag(:,:,i) * model.H'; %Eq.(3)
%model.Sigma(:,:,i) = model.params_alpha*(model.H*model.SigmaDiag(:,:,i)*model.H') + (1 - model.params_alpha)*model.S(:,:,i); % Eq. (10)
end end
%Compute average log-likelihood %Compute average log-likelihood
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment