wolffd@0
|
1 function bnet = mk_alarm_bnet()
|
wolffd@0
|
2
|
wolffd@0
|
3 % Written by Qian Diao <qian.diao@intel.com> on 11 Dec 01
|
wolffd@0
|
4
|
wolffd@0
|
5 N = 37;
|
wolffd@0
|
6 dag = zeros(N,N);
|
wolffd@0
|
7 dag(21,23) = 1 ;
|
wolffd@0
|
8 dag(21,24) = 1 ;
|
wolffd@0
|
9 dag(1,24) = 1 ;
|
wolffd@0
|
10 dag(1,23) = 1 ;
|
wolffd@0
|
11 dag(2,26) = 1 ;
|
wolffd@0
|
12 dag(2,25) = 1 ;
|
wolffd@0
|
13 dag(2,24) = 1 ;
|
wolffd@0
|
14 dag(2,13) = 1 ;
|
wolffd@0
|
15 dag(2,23) = 1 ;
|
wolffd@0
|
16 dag(13,30) = 1 ;
|
wolffd@0
|
17 dag(30,31) = 1 ;
|
wolffd@0
|
18 dag(3,14) = 1 ;
|
wolffd@0
|
19 dag(3,19) = 1 ;
|
wolffd@0
|
20 dag(4,36) = 1 ;
|
wolffd@0
|
21 dag(14,35) = 1 ;
|
wolffd@0
|
22 dag(32,33) = 1 ;
|
wolffd@0
|
23 dag(32,35) = 1 ;
|
wolffd@0
|
24 dag(32,34) = 1 ;
|
wolffd@0
|
25 dag(32,36) = 1 ;
|
wolffd@0
|
26 dag(15,21) = 1 ;
|
wolffd@0
|
27 dag(5,31) = 1 ;
|
wolffd@0
|
28 dag(27,30) = 1 ;
|
wolffd@0
|
29 dag(28,31) = 1 ;
|
wolffd@0
|
30 dag(28,29) = 1 ;
|
wolffd@0
|
31 dag(26,28) = 1 ;
|
wolffd@0
|
32 dag(26,27) = 1 ;
|
wolffd@0
|
33 dag(16,31) = 1 ;
|
wolffd@0
|
34 dag(16,37) = 1 ;
|
wolffd@0
|
35 dag(23,26) = 1 ;
|
wolffd@0
|
36 dag(23,29) = 1 ;
|
wolffd@0
|
37 dag(23,25) = 1 ;
|
wolffd@0
|
38 dag(6,15) = 1 ;
|
wolffd@0
|
39 dag(7,27) = 1 ;
|
wolffd@0
|
40 dag(8,21) = 1 ;
|
wolffd@0
|
41 dag(19,20) = 1 ;
|
wolffd@0
|
42 dag(19,22) = 1 ;
|
wolffd@0
|
43 dag(31,32) = 1 ;
|
wolffd@0
|
44 dag(9,14) = 1 ;
|
wolffd@0
|
45 dag(9,17) = 1 ;
|
wolffd@0
|
46 dag(9,19) = 1 ;
|
wolffd@0
|
47 dag(10,33) = 1 ;
|
wolffd@0
|
48 dag(10,34) = 1 ;
|
wolffd@0
|
49 dag(11,16) = 1 ;
|
wolffd@0
|
50 dag(12,13) = 1 ;
|
wolffd@0
|
51 dag(12,18) = 1 ;
|
wolffd@0
|
52 dag(35,37) = 1 ;
|
wolffd@0
|
53
|
wolffd@0
|
54 node_sizes = 2*ones(1,N);
|
wolffd@0
|
55 node_sizes(2) = 3;
|
wolffd@0
|
56 node_sizes(6) = 3;
|
wolffd@0
|
57 node_sizes(14) = 3;
|
wolffd@0
|
58 node_sizes(15) = 4;
|
wolffd@0
|
59 node_sizes(16) = 3;
|
wolffd@0
|
60 node_sizes(18) = 3;
|
wolffd@0
|
61 node_sizes(19) = 3;
|
wolffd@0
|
62 node_sizes(20) = 3;
|
wolffd@0
|
63 node_sizes(21) = 4;
|
wolffd@0
|
64 node_sizes(22) = 3;
|
wolffd@0
|
65 node_sizes(23) = 4;
|
wolffd@0
|
66 node_sizes(24) = 4;
|
wolffd@0
|
67 node_sizes(25) = 4;
|
wolffd@0
|
68 node_sizes(26) = 4;
|
wolffd@0
|
69 node_sizes(27) = 3;
|
wolffd@0
|
70 node_sizes(28) = 3;
|
wolffd@0
|
71 node_sizes(29) = 4;
|
wolffd@0
|
72 node_sizes(30) = 3;
|
wolffd@0
|
73 node_sizes(32) = 3;
|
wolffd@0
|
74 node_sizes(33) = 3;
|
wolffd@0
|
75 node_sizes(34) = 3;
|
wolffd@0
|
76 node_sizes(35) = 3;
|
wolffd@0
|
77 node_sizes(36) = 3;
|
wolffd@0
|
78 node_sizes(37) = 3;
|
wolffd@0
|
79
|
wolffd@0
|
80 bnet = mk_bnet(dag, node_sizes);
|
wolffd@0
|
81
|
wolffd@0
|
82 bnet.CPD{1} = tabular_CPD(bnet, 1,[0.96 0.04 ]);
|
wolffd@0
|
83 bnet.CPD{2} = tabular_CPD(bnet, 2,[0.92 0.03 0.05 ]);
|
wolffd@0
|
84 bnet.CPD{3} = tabular_CPD(bnet, 3,[0.8 0.2 ]);
|
wolffd@0
|
85 bnet.CPD{4} = tabular_CPD(bnet, 4,[0.95 0.05 ]);
|
wolffd@0
|
86 bnet.CPD{5} = tabular_CPD(bnet, 5,[0.8 0.2 ]);
|
wolffd@0
|
87 bnet.CPD{6} = tabular_CPD(bnet, 6,[0.01 0.98 0.01 ]);
|
wolffd@0
|
88 bnet.CPD{7} = tabular_CPD(bnet, 7,[0.01 0.99 ]);
|
wolffd@0
|
89 bnet.CPD{8} = tabular_CPD(bnet, 8,[0.95 0.05 ]);
|
wolffd@0
|
90 bnet.CPD{9} = tabular_CPD(bnet, 9,[0.95 0.05 ]);
|
wolffd@0
|
91 bnet.CPD{10} = tabular_CPD(bnet, 10,[0.9 0.1 ]);
|
wolffd@0
|
92 bnet.CPD{11} = tabular_CPD(bnet, 11,[0.99 0.01 ]);
|
wolffd@0
|
93 bnet.CPD{12} = tabular_CPD(bnet, 12,[0.99 0.01 ]);
|
wolffd@0
|
94 bnet.CPD{13} = tabular_CPD(bnet, 13,[0.95 0.95 0.05 0.1 0.1 0.01 0.05 0.05 0.95 0.9 0.9 0.99 ]);
|
wolffd@0
|
95 bnet.CPD{14} = tabular_CPD(bnet, 14,[0.05 0.95 0.5 0.98 0.9 0.04 0.49 0.01 0.05 0.01 0.01 0.01 ]);
|
wolffd@0
|
96 bnet.CPD{15} = tabular_CPD(bnet, 15,[0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.97 ]);
|
wolffd@0
|
97 bnet.CPD{16} = tabular_CPD(bnet, 16,[0.3 0.98 0.4 0.01 0.3 0.01 ]);
|
wolffd@0
|
98 bnet.CPD{17} = tabular_CPD(bnet, 17,[0.99 0.1 0.01 0.9 ]);
|
wolffd@0
|
99 bnet.CPD{18} = tabular_CPD(bnet, 18,[0.05 0.01 0.9 0.19 0.05 0.8 ]);
|
wolffd@0
|
100 bnet.CPD{19} = tabular_CPD(bnet, 19,[0.05 0.98 0.01 0.95 0.9 0.01 0.09 0.04 0.05 0.01 0.9 0.01 ]);
|
wolffd@0
|
101 bnet.CPD{20} = tabular_CPD(bnet, 20,[0.95 0.04 0.01 0.04 0.95 0.29 0.01 0.01 0.7 ]);
|
wolffd@0
|
102 bnet.CPD{21} = tabular_CPD(bnet, 21,[0.97 0.97 0.01 0.97 0.01 0.97 0.01 0.97 0.01 0.01 0.97 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.01 ]);
|
wolffd@0
|
103 bnet.CPD{22} = tabular_CPD(bnet, 22,[0.95 0.04 0.01 0.04 0.95 0.04 0.01 0.01 0.95 ]);
|
wolffd@0
|
104 bnet.CPD{23} = tabular_CPD(bnet, 23,[0.97 0.97 0.97 0.97 0.97 0.97 0.01 0.95 0.97 0.97 0.01 0.95 0.01 0.4 0.97 0.97 0.01 0.5 0.01 0.3 0.97 0.97 0.01 0.3 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.03 0.01 0.01 0.97 0.03 0.01 0.58 0.01 0.01 0.01 0.48 0.01 0.68 0.01 0.01 0.01 0.68 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.97 0.01 ]);
|
wolffd@0
|
105 bnet.CPD{24} = tabular_CPD(bnet, 24,[0.97 0.97 0.97 0.97 0.97 0.97 0.01 0.01 0.4 0.1 0.01 0.01 0.01 0.01 0.2 0.05 0.01 0.01 0.01 0.01 0.2 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.49 0.58 0.84 0.9 0.29 0.01 0.01 0.75 0.25 0.01 0.01 0.01 0.01 0.7 0.15 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.3 0.01 0.05 0.08 0.3 0.97 0.08 0.04 0.25 0.38 0.08 0.01 0.01 0.09 0.25 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.2 0.01 0.01 0.01 0.4 0.01 0.9 0.01 0.45 0.6 0.9 0.97 0.97 0.01 0.59 0.97 0.97 ]);
|
wolffd@0
|
106 bnet.CPD{25} = tabular_CPD(bnet, 25,[0.97 0.97 0.97 0.01 0.6 0.01 0.01 0.5 0.01 0.01 0.5 0.01 0.01 0.01 0.01 0.97 0.38 0.97 0.01 0.48 0.01 0.01 0.48 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.01 0.97 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.01 0.97 ]);
|
wolffd@0
|
107 bnet.CPD{26} = tabular_CPD(bnet, 26,[0.97 0.97 0.97 0.01 0.01 0.03 0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.01 0.01 0.97 0.97 0.95 0.01 0.01 0.94 0.01 0.01 0.88 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.97 0.04 0.01 0.01 0.1 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.97 0.01 ]);
|
wolffd@0
|
108 bnet.CPD{27} = tabular_CPD(bnet, 27,[0.98 0.98 0.98 0.98 0.95 0.01 0.95 0.01 0.01 0.01 0.01 0.01 0.04 0.95 0.04 0.01 0.01 0.01 0.01 0.01 0.01 0.04 0.01 0.98 ]);
|
wolffd@0
|
109 bnet.CPD{28} = tabular_CPD(bnet, 28,[0.01 0.01 0.04 0.9 0.01 0.01 0.92 0.09 0.98 0.98 0.04 0.01 ]);
|
wolffd@0
|
110 bnet.CPD{29} = tabular_CPD(bnet, 29,[0.97 0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.97 0.01 0.01 0.01 0.01 0.97 0.97 0.97 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.97 0.97 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.97 0.97 0.43 ]);
|
wolffd@0
|
111 bnet.CPD{30} = tabular_CPD(bnet, 30,[0.98 0.98 0.01 0.98 0.01 0.69 0.01 0.01 0.98 0.01 0.01 0.3 0.01 0.01 0.01 0.01 0.98 0.01 ]);
|
wolffd@0
|
112 bnet.CPD{31} = tabular_CPD(bnet, 31,[0.05 0.01 0.05 0.01 0.05 0.01 0.05 0.01 0.05 0.01 0.05 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.1 0.01 0.95 0.01 0.95 0.05 0.1 0.01 0.95 0.01 0.95 0.05 0.1 0.01 0.3 0.01 0.3 0.01 0.95 0.01 0.99 0.05 0.95 0.05 0.95 0.01 0.99 0.05 0.99 0.05 0.3 0.01 0.99 0.01 0.3 0.01 0.95 0.99 0.95 0.99 0.95 0.99 0.95 0.99 0.95 0.99 0.95 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.9 0.99 0.05 0.99 0.05 0.95 0.9 0.99 0.05 0.99 0.05 0.95 0.9 0.99 0.7 0.99 0.7 0.99 0.05 0.99 0.00999999 0.95 0.05 0.95 0.05 0.99 0.01 0.95 0.01 0.95 0.7 0.99 0.01 0.99 0.7 0.99 ]);
|
wolffd@0
|
113 bnet.CPD{32} = tabular_CPD(bnet, 32,[0.1 0.01 0.89 0.09 0.01 0.9 ]);
|
wolffd@0
|
114 bnet.CPD{33} = tabular_CPD(bnet, 33,[0.98 0.33333334 0.01 0.33333334 0.01 0.33333334 0.01 0.33333334 0.98 0.33333334 0.01 0.33333334 0.01 0.33333334 0.01 0.33333334 0.98 0.33333334 ]);
|
wolffd@0
|
115 bnet.CPD{34} = tabular_CPD(bnet, 34,[0.98 0.33333334 0.01 0.33333334 0.01 0.33333334 0.01 0.33333334 0.98 0.33333334 0.01 0.33333334 0.01 0.33333334 0.01 0.33333334 0.98 0.33333334 ]);
|
wolffd@0
|
116 bnet.CPD{35} = tabular_CPD(bnet, 35,[0.98 0.95 0.3 0.95 0.04 0.01 0.8 0.01 0.01 0.01 0.04 0.69 0.04 0.95 0.3 0.19 0.04 0.01 0.01 0.01 0.01 0.01 0.01 0.69 0.01 0.95 0.98 ]);
|
wolffd@0
|
117 bnet.CPD{36} = tabular_CPD(bnet, 36,[0.98 0.98 0.01 0.4 0.01 0.3 0.01 0.01 0.98 0.59 0.01 0.4 0.01 0.01 0.01 0.01 0.98 0.3 ]);
|
wolffd@0
|
118 bnet.CPD{37} = tabular_CPD(bnet, 37,[0.98 0.98 0.3 0.98 0.1 0.05 0.9 0.05 0.01 0.01 0.01 0.6 0.01 0.85 0.4 0.09 0.2 0.09 0.01 0.01 0.1 0.01 0.05 0.55 0.01 0.75 0.9 ]);
|