贝叶斯网络matlab程序

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（1N）;
%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 ] ）;