本发明属于网络可靠性分析技术领域，具体涉及一种基于边扩展图指挥控制网络二端可靠性评估方法。
背景技术：
传统的对网络可靠性的研究通常假设其节点完全可靠，其失效率为零。但是实际网络中这种假设是不成立的，因为往往构成网络节点的物体可能发生故障，如通信网中的电台路由设备、电力网络中的发电输电设备、交通网络中的车辆等。当网络规模很大时，即时节点失效率非常小，由于节点数量的庞大，也会造成对网络连通的巨大影响，因此研究节点不可靠状态下的网络可靠性是非常有必要的。
技术实现要素：
本发明的发明目的是：为了解决现有技术中存在的以上问题，本发明提出了一种基于边扩展图指挥控制网络二端可靠性评估方法。
本发明的技术方案是：一种基于边扩展图指挥控制网络二端可靠性评估方法，包括以下步骤：
a、在节点完全可靠条件下将指挥控制网络抽象为节点和边的拓扑结构图，并进行简化处理；
b、根据步骤a得到的拓扑结构图进行指挥控制网络边扩展处理，生成原始路径函数；
c、在节点不完全可靠条件下对步骤b得到的原始路径函数进行关联边替换处理，构造新路径函数；
d、采用宽度优先搜索方法对步骤a中拓扑结构图的节点和边进行变量排序处理；
e、将步骤c得到的新路径函数按照步骤d中变量排序顺序构造二元决策图；
f、根据步骤e中得到的二元决策图由下往上递归计算指挥控制网络二端网络可靠度。
进一步地，所述步骤b根据步骤a得到的拓扑结构图进行指挥控制网络边扩展处理，生成原始路径函数，具体包括以下分步骤：
b1、从步骤a得到的拓扑结构图中的源点出发，遍历网络中与源点关联的边，分别沿与源点关联的边进行网络遍历，通过点收缩和边删除处理得到多个子网络；
b2、依次对得到的每一个子网络根据其中源点的位置再次遍历其关联边，重复步骤b1，直至拓扑结构图中源点与汇点重合；
b3、构造在节点完全可靠条件下的原始路径函数；表示为：
p(g)＝x1p(g1)+x2p(g2)+…+xkp(gk)
其中，p(g)为原始路径函数，g为指挥控制网络，xi(i＝1,2,…,k)为与源点关联的边，
进一步地，所述步骤c在节点不完全可靠条件下对步骤b得到的原始路径函数进行关联边替换处理，构造新路径函数，具体为：
在节点不完全可靠条件下将拓扑结构图中的边等价为边与其两端节点的交，将步骤b得到的原始路径函数中的边替换为对应边与其两端节点的交，进行布尔运算得到新路径函数。
进一步地，所述步骤d采用宽度优先搜索方法对步骤a中拓扑结构图的节点和边进行变量排序处理，具体为：
从源点开始依次访问与源点直接相连的各个未被访问过的节点，按照被访问的先后顺序对与源点相关联的边进行标号作为该边的边序号，再分别以与源点直接相连的各个未被访问过的节点为起点重复上述操作，直至拓扑结构图中所有节点都被访问过为止，将拓扑结构图的节点和边按照边序号进行排列，得到变量排序顺序。
进一步地，所述步骤e将步骤c得到的新路径函数按照步骤d中变量排序顺序构造二元决策图，具体为：将步骤c得到的新路径函数转化为ite结构形式
并满足
其中，f为转化为ite结构形式的新路径函数，
进一步地，所述步骤f根据步骤e中得到的二元决策图由下往上递归计算指挥控制网络二端网络可靠度的计算公式为：
其中，r(obdd(g))为指挥控制网络二端网络可靠度,
本发明的有益效果是：本发明通过将指挥控制网络抽象为节点和边的拓扑结构图并进行边扩展处理生成路径函数，再采用宽度优先搜索方法构造二元决策图，能够通过识别同构子图来避免冗余计算，这对网络二端可靠性分析研究具有重要的理论和实际的应用价值。
附图说明
图1是本发明的基于边扩展图指挥控制网络二端可靠性评估方法的流程示意图；
图2是本发明实施例中指挥控制网络拓扑结构图；
图3是本发明实施例中简化处理后的指挥控制网络拓扑结构图；
图4是本发明实施例中指挥控制网络边扩展示意图；
图5是本发明实施例中指挥控制网络路径函数构造示意图；
图6是本发明实施例中关联边替换示意图；
图7是本发明实施例中二元决策图；
图8是本发明实施例中指挥控制网络二端连通二元决策结构图；
图9是本发明实施例中指挥控制网络二端可靠性计算示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
如图1所示，为本发明的基于边扩展图指挥控制网络二端可靠性评估方法的流程示意图。一种基于边扩展图指挥控制网络二端可靠性评估方法，包括以下步骤：
a、在节点完全可靠条件下将指挥控制网络抽象为节点和边的拓扑结构图，并进行简化处理；
b、根据步骤a得到的拓扑结构图进行指挥控制网络边扩展处理，生成原始路径函数；
c、在节点不完全可靠条件下对步骤b得到的原始路径函数进行关联边替换处理，构造新路径函数；
d、采用宽度优先搜索方法对步骤a中拓扑结构图的节点和边进行变量排序处理；
e、将步骤c得到的新路径函数按照步骤d中变量排序顺序构造二元决策图；
f、根据步骤e中得到的二元决策图由下往上递归计算指挥控制网络二端网络可靠度。
在本发明的一个可选实施例中，如图2所示，为本发明实施例中指挥控制网络拓扑结构图；其中各节点表示指挥控制网络中的通信车，边表示车与车之间的通信，不同边表示通信方式的不同。本发明将图中两个节点设为需要通信连接的二端，对指挥控制网络的二端可靠性进行研究；由于图中的有些冗余节点和冗余边是可有可无的，并不影响网络通信连接的效果，而网络中有的节点和边的连接关系能够用简单简化规则将其化简，如串联简化和并联简化，因此将指挥控制网络拓扑结构图进行图形处理，采用简化规则和冗余删除操作将其转化为最简形式。如图3所示，为本发明实施例中简化处理后的指挥控制网络拓扑结构图。
最简形式的指挥控制网络中，其节点和边的失效率随着简化处理而改变。将其以图形表示方法表示为g(e,v,j2,j4)，其中e＝(e1,e1,…,e10),v＝(j1,j2,…,j5)分别表示边集和节点集，边数为10，节点数为5。
在本发明的一个可选实施例中，本发明通过生成二端网络的路径函数对网络二端可靠性进行评估。上述步骤b根据步骤a得到的拓扑结构图进行指挥控制网络边扩展处理，生成原始路径函数，具体包括以下分步骤：
b1、从步骤a得到的拓扑结构图中的源点出发，遍历网络中与源点关联的边，分别沿与源点关联的边进行网络遍历，通过点收缩和边删除处理得到多个子网络；
从源点j2出发，遍历网络g中与源点j2关联的边，设为xi(i＝1,2,…,k)。分别沿着xi(i＝1,2,…,k)进行网络遍历，通过点收缩和边删除得到相应的k各子网gi(i＝1,2,…,k),所谓点收缩，是指将源点与关联边xi(i＝1,2,…,k)的另一个端点进行合并，并将合并后的节点作为新的源点，而边删除是指将与原来源点关联的边删除。
b2、依次对得到的每一个子网络根据其中源点的位置再次遍历其关联边，重复步骤b1，直至拓扑结构图中源点与汇点重合；
对得到的每一个子网gi(i＝1,2,…,k)根据其中源点的位置再次遍历其关联边，重复上述操作，依次进行直到最终源点与汇点重合。
b3、构造在节点完全可靠条件下的原始路径函数；表示为：
p(g)＝x1p(g1)+x2p(g2)+…+xkp(gk)
其中，p(g)为原始路径函数，g为指挥控制网络，xi(i＝1,2,…,k)为与源点关联的边，
在指挥控制网络中，假设节点ji，对图2进行边扩展操作，如图4所示，为本发明实施例中指挥控制网络边扩展示意图。根据网络二端通信的路径函数的构造方法，对指挥控制网络二端通信的路径函数进行构造。如图5所示，为本发明实施例中指挥控制网络路径函数构造示意图。从而得到考虑节点完全可靠条件下指挥控制网络二端通信的路径函数，表示为
p(g)＝e2+e1e7+e3e9+e4e8+e1e5e8+e1e6e9+e3e6e7
+e3e8e10+e4e5e7+e4e9e10+e1e5e9e10+e1e6e8e10
+e3e5e6e8+e3e5e7e10+e4e5e6e9+e4e6e7e10
在本发明的一个可选实施例中，上述步骤c在节点不完全可靠条件下对步骤b得到的原始路径函数进行关联边替换处理，构造新路径函数，具体为：
在节点不完全可靠条件下将拓扑结构图中的边等价为边与其两端节点的交，将步骤b得到的原始路径函数中的边替换为对应边与其两端节点的交，进行布尔运算得到新路径函数。
本发明在考虑节点存在失效概率条件下，对关联边替换。如图6所示，为本发明实施例中关联边替换示意图；对于网络拓扑图中的边em，其两端节点分别为va和vb。在考虑边em，节点va和vb均存在失效概率的情况下，可将其等价为三者的交，意味着三者必须均可靠该链路才可靠，其中一个节点或者边失效则该链路失效。因此可将之前节点完全可靠条件下指挥控制网络二端通信的路径函数中的边替换为相应的节点和边的交。指挥控制网络中的关联边替换情况具体为
e1→j2∧e1∧j3,e2→j2∧e2∧j4,e3→j2∧e3∧j5
e4→j2∧e4∧j1,e5→j3∧e5∧j1,e6→j3∧e6∧j5
e7→j3∧e7∧j4,e8→j4∧e8∧j1,e9→j4∧e9∧j5
e10→j1∧e10∧j5
由此替换关系并将其带入原路径函数，进行相关的布尔运算得到如下考虑节点不完全可靠条件下指挥控制网络二端通信的路径函数，表示为
p(g)＝j2e2j4+j2e1j3e7j4+j2e3j5e9j4+j2e4j1e8j4
+j2e1j3e5j1e8j4+j2e1j3e6j5e9j4+j2e3j5e6j3e7j4
+j2e3j5e10j1e8j4+j2e4j1e5j3e7j4+j2e4j1e10j5e9j4
+j2e1j3e5j1e10j5e9j4+j2e1j3e6j5e10j1e8j4+j2e3j5e6j3e5j1e8j4
+j2e3j5e10j1e5j3e7j4+j2e4j1e5j3e6j5e9j4+j2e4j1e10j5e6j3e7j4
在本发明的一个可选实施例中，上述步骤d采用宽度优先搜索方法对步骤a中拓扑结构图的节点和边进行变量排序处理，具体为：
从源点开始依次访问与源点直接相连的各个未被访问过的节点，按照被访问的先后顺序对与源点相关联的边进行标号作为该边的边序号，再分别以与源点直接相连的各个未被访问过的节点为起点重复上述操作，直至拓扑结构图中所有节点都被访问过为止，将拓扑结构图的节点和边按照边序号进行排列，得到变量排序顺序。
本发明从源点j2开始依次访问与j2点直接相连的各个未被访问过的节点v1,v2,…,vn，按照被访问的先后顺序对与源点相关联的边<j2,v1>,…,<j2,vn>进行标号，分别为e1,e2,…,en，此标号即为该边的边序号；接着分别以v1,v2,…,vn为起点访问与v1,v2,…,vn直接相连的节点，重复上述操作，一直到图中所有节点都被访问过为止，图中被标号的变量依据标号排列起来形成了宽度优先搜索排序下的变量顺序。
因此根据宽度优先搜索排序方法对指挥控制网络中节点和边进行排序，得到变量顺序为
j2＜e1＜e2＜e3＜e4＜j3＜j4＜j5＜j1＜e5＜e6＜e7＜e8＜e9＜e10
在本发明的一个可选实施例中，上述步骤e将步骤c得到的新路径函数按照步骤d中变量排序顺序构造二元决策图，具体为：将步骤c得到的新路径函数转化为另外的结构形式
其中，f为转化为ite结构形式的新路径函数，
并满足
其中，g,h分别为两个子函数，<op>为布尔操作“与”、“或”，index(xi)为变量xi的标号，
本发明采用基于shannon分解的转化规则由指挥控制网络二端通信的路径函数构造相应的二元决策图(obbd)，其中
在本发明的一个可选实施例中，上述步骤f根据步骤e中得到的二元决策图由下往上递归计算指挥控制网络二端网络可靠度的计算公式为：
其中，r(obdd(g))为指挥控制网络二端网络可靠度，
本发明根据构造的指挥控制网络二端通信的obdd结构，采用上式递归计算其网络可靠度值。
如图9所示，为本发明实施例中指挥控制网络二端可靠性计算示意图，可以看到，当假设变量可靠度均为0.9时，从obdd结构最底层开始从下向上按照变量顺序的反序依次递归计算可靠度，最后在排序的第一个变量节点出计算所得的可靠度即为整个网络的可靠度。
假设进行图形处理后的指挥控制网络各节点和边的可靠度如表1所示
表1指挥控制网络节点及边可靠度
则由obdd结构图由下至上递归计算指挥控制网络二端可靠性为p＝0.9988。
本发明针对网络可靠性分析，运用有序的bdd方法(obdd)提出了一种更好地评估方法，这种方法利用网络的边扩展(eed)图来构建网络的obdd结构图，能够通过识别同构子图来避免冗余计算。该方法将k端问题分解为若干个两端可靠性问题，在两端可靠性函数基础上利用obdd来进行计算，这种方法相对于因子分解方法来说减少了多余的计算，提高了计算效率。对理论分析以及实际的工程应用方面均有实际的意义。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。