贝叶斯网络matlab编程例子？

[azyc414]

谢谢谢谢啊

[Tony2251]

那么专业。
有算法自己实现一下就是了。

[HimeTale]

N = 4;
dag = zeros( N, N );
C = 1; 
S = 2;
R = 3;
W = 4;
dag( C, [R S] ) = 1;
dag( R, W ) = 1;
dag( S, W ) = 1;
 
node_sizes = 2*ones(1,N);
%node_sizes = [ 4 2 3 5];
 
%建立有向无环图
onodes = [];
bNet = mk_bnet( dag, node_sizes );
 
%建立条件概率表
bNet.CPD{C} = tabular_CPD( bNet, C, 'CPT', [0.5 0.5] );
bNet.CPD{R} = tabular_CPD( bNet, R,  'CPT',[0.8 0.2 0.2 0.8 ] );
bNet.CPD{S} = tabular_CPD( bNet, S,  'CPT', [0.5 0.9 0.5 0.1] );
bNet.CPD{W} = tabular_CPD( bNet, W, 'CPT', [1 0.1 0.1 0.01 0 0.9 0.9 0.99 ] );
     
 
 
%用于推断
engine = jtree_inf_engine( bNet);
 
evidence = cell( 1,N );
evidence{W} = 2;
 
%添加证据
[engine, loglik] = enter_evidence( engine, evidence );
marg = marginal_nodes( engine, S );
p = marg.T(2);
 
%添加下雨的证据
evidence{R} = 2;
[ engine, loglik ] = enter_evidence( engine, evidence );
marg = marginal_nodes( engine, S );
p = marg.T(2);
 
bar( marg.T );
 
 
% Joint distributions
evidence = cell(1,N);
[engine, ll] = enter_evidence(engine, evidence);
m = marginal_nodes(engine, [S R W]);
m.T

[Lvyiqing]

想学习一下，看看