example_3.m.svn-base

来自「Probabilistic graphical models in matlab」· SVN-BASE 代码 · 共 15 行

my_network = new_model('factor', 6 );
my_network = set_all_variables_sizes(my_network, 2);
M = [0.2 0.1; 0.7 0.9];
M (:,:,2) = [0.4 0.3; 0.6 0.4];
my_network = set_potential(my_network, [1 3 5], M);
my_network = set_potential(my_network, [1 2], [2 0.6; 1.1 1.3] );
M = [0.8 0.6; 0.4 0.2];
M (:,:,2) = [0.4 0.5; 0.5 0.05];
my_network = set_potential(my_network, [3 4 6], M );
my_network = set_observation(my_network, 2 , 2 );
my_network = use_timer(my_network);
my_network = compute_inference(my_network, 'gibbs', 5);
my_network.marginals{:}
my_network = compute_inference(my_network, 'tree_mcmc', 2);
my_network.marginals{:}