computevard.m

来自「Nonlinear dynamical factor analysis Matl」· M 代码 · 共 23 行

function [dcp_dsvn, newac] = computevard(dcp_dsv, dcp_dsvn, ac, mv)
% COMPUTEVARD is used by by ndfa_iter.  This function is the feedback version
% of COMPUTEVAR.

% Copyright (C) 2002 Harri Valpola and Antti Honkela.
%
% This package comes with ABSOLUTELY NO WARRANTY; for details
% see License.txt in the program package.  This is free software,
% and you are welcome to redistribute it under certain conditions;
% see License.txt for details.

newac=zeros(size(ac));
d=size(newac,2);
for i=1:size(mv,1)
  newac(i,2:end)=reshape(mv(i,i,:), [1 d-1]);
end

for i=d:-1:2
  newac(:,i) = newac(:,i) .* dcp_dsvn(:,i) ./ (dcp_dsvn(:,i) + dcp_dsv(:,i));
  dcp_dsvn(:,i) = dcp_dsvn(:,i) + dcp_dsv(:,i);
  dcp_dsv(:,i-1) = dcp_dsv(:,i-1) + newac(:,i).^2 .* dcp_dsv(:,i);
end
dcp_dsvn(:,1) = dcp_dsvn(:,1) + dcp_dsv(:,1);