贝叶斯网络模型代码.doc

贝叶斯网络模型代码贝叶斯网络模型代码贝叶斯网络模型代码addpath(genpathKPM(pwd））N=4;dag=zeros(N，N)；C=1;S=2；R=3；W=4；dag（C，［RS］)=1;dag（R，W)=1;dag（S,W）=1;discrete_nodes=1：N;node_sizes=2*ones（1,N）;bnet=mk_bnet(dag,node_sizes，'discrete'，discrete_nodes）；onodes=[];bnet=mk_bnet(dag,node_sizes，'discrete’，discrete_nodes，’observed’，onodes）；bnet=mk_bnet（dag,node_sizes,’names’,{’cloudy’,'S’，'R’,’W'｝，’discrete',1：4)；C=bnet.names('cloudy’)；％bnet.names是一个关联数组；bnet.CPD{C｝=tabular_CPD（bnet,C，[0.50.5]);CPT=zeros（2,2,2);CPT(1，1,1）=1.0;CPT(2,1,1）=0.1；CPT=reshape（[10.10.10。0100。90。90.99],[222］）；bnet.CPD｛W｝=tabular_CPD(bnet，W，'CPT’，[10。10。10。0100。90.90。99］）；bnet。CPD{C｝=tabular_CPD（bnet,C，［0。50.5］)；bnet。CPD{R｝=tabular_CPD(bnet，R，［0.80.20。20。8］）；bnet。CPD｛S}=tabular_CPD（bnet，S，［0。50.90.50.1］）;bnet.CPD{W}=tabular_CPD(bnet，W，［10.10。10。0100.90.90.99］）;figuredraw_graph（dag)