fwd1.m

来自「麻省理工学院的人工智能工具箱,很珍贵,希望对大家有用!」· M 代码 · 共 28 行

function [f, logscale] = fwd1(engine, ev, t)
% Forwards pass for slice 1.

bnet = bnet_from_engine(engine);
ss = bnet.nnodes_per_slice;

CPDpot = cell(1,ss);
for n=1:ss
  fam = family(bnet.dag, n, 1);
  e = bnet.equiv_class(n, 1);
  CPDpot{n} = convert_to_table(bnet.CPD{e}, fam(:), ev);
end       
f.evidence = ev;
f.t = t;

pots = CPDpot;
slice1 = 1:ss;
%CPDclqs = engine.clq_ass_to_node1(slice1);

[f.clpot, f.seppot] =  init_pot(engine.jtree_engine1, slice1, CPDpot);
[f.clpot, f.seppot] = collect_evidence(engine.jtree_engine1, f.clpot, f.seppot);
for c=1:length(f.clpot)
   [f.clpot{c}, lik] = normalise(f.clpot{c});
   ll(c) = log(lik + (lik==0)*eps);
end
logscale = ll(engine.root1);