基于截断边扩展图的工程网络可靠度近似分析方法
【技术领域】
[0001] 本发明属于网络可靠性分析技术领域，具体涉及一种基于截断边扩展图的工程网 络可靠度近似分析方法。
【背景技术】
[0002] 随着科学技术的飞速发展以及计算机的广泛普及，计算机网络已经在通信、银行、 交通、物流、能源、学校、工厂和军事等领域得到了广泛应用。随着计算机网络应用领域的不 断扩大，人们对计算机网络的可靠性要求也越来越高。一旦计算机网络发生故障，将会对人 们的生活、经济、环境、卫生及社会安全造成重大甚至是灾难性的影响，因此，研究网络可靠 性具有极其重要的现实意义。
[0003] 计算机网络可靠性主要指双端可靠性，即源端S到终端T连通的概率。目前，对计算 机网络可靠性的分析方法主要为精确分析方法，虽然其求解值精度较高，但是，由于分析计 算过程非常繁琐，既消耗了大量计算资源，也延长了可靠性分析时间。所以，传统的可靠性 分析方法只适用于中小型网络，无法应用于中大型网络。
【发明内容】
[0004] 针对现有技术存在的缺陷，本发明提供一种基于截断边扩展图的工程网络可靠度 近似分析方法，以牺牲最小可靠性精度为代价，最大程度简化了工程网络可靠度的分析繁 琐性，可适用于对中大型网络进行可靠性分析。
[0005] 本发明采用的技术方案如下：
[0006] 本发明提供一种基于截断边扩展图的工程网络可靠度近似分析方法，包括以下步 骤：
[0007] 步骤1，对于给定网络G = (V，E，S，T);其中，V为节点集合，E为边集合;S为源节点，T 为汇节点；
[0008] 步骤2,当需要计算源节点S到汇节点T之间的网络可靠度值时，首先遍历给定网络 G，对边集合中的各条边进行排序；
[0009] 步骤3,给定截断深度极大值N;其中，N为自然数;给定截断深度变量值P;建立初始 为空的记录表；
[0010] 初始建立仅存在根节点、并且根节点为给定网络G的树模型;此时，给定网络G为第 0层节点；
[0011] 步骤4,设截断深度变量值P = l;
[0012] 步骤5,判断截断深度变量值P是否小于截断深度极大值N，如果小于，则执行步骤 6;如果大于等于，则执行步骤10;
[0013] 步骤6，当前时刻，设树模型第P-1层共存在X个节点，分别记为G1、G2."Gx;
[0014] 步骤7,令j = l;
[0015] 步骤8,树模型扩展出第P层节点的方法为：
[0016] 步骤8.1，对于树模型第p-l层节点中的任意节点Gj，设网络Gj中与源节点S直接连 接的边为xi(i = 1，2，. . .，k)，其中，k为自然数，xi代表边的排序序号值；
[0017] 步骤8.2,首先，令i = l;
[0018] 步骤8.3,对于网络Gj，将源节点S收缩到边xi的另一端点并作为新的源节点，同 时，删除所有与源节点S直接连接的边、度为1的非S节点、与度为1的非S节点直接相连的边， 得到网络Gj的子网G*xi，其中，G*xi为得到的子网标识；同时，判断所述记录表中是否存在 与子网G*xi的网络结构完全相同的子网，如果不存在，则执行步骤8.4;如果存在，则执行步 骤8.5;
[0019] 步骤8.4,将新得到的子网G*xi的子网网络结构、子网标识以及子网在给定网络G 的深度值的信息记录到所述记录表;其中，子网在给定网络G的深度值是指:从给定网络G到 子网G*x i所经过的边的数量值；
[0020]同时，向树模型中增加子网G*xi，并连接网络Gj节点和子网G*xi，子网G*xi形成为 网络G j的一个儿子节点，子网G*xi的深度值即为P;
[0021] 然后，执行步骤9;
[0022] 步骤8.5，假设记录表中已存在子网G*xt，其中，xt关xi ;此时，子网G*xt和子网G* xi的网络结构完全相同，子网G*xt和子网G*xi为同构子图；
[0023] 设子网G*xt的深度值为Nxt，子网G*xi的深度值为Nxi，判断Nxi是否小于N xt，如果小 于，则执行步骤8.6;否则，执行步骤8.7;
[0024] 步骤8.6,将记录表中已记录的子网G*xt的子网网络结构、子网标识以及子网在给 定网络G的深度值的信息删除；同时，将新得到的子网G*xi的子网网络结构、子网标识以及 子网在给定网络G的深度值的信息新增加到记录表；
[0025]同时，向树模型中增加子网G*xi，并连接网络Gj节点和子网G*xi，子网G*xi形成为 网络Gj的一个儿子节点，子网G*xi的深度值即为P;此外，将树模型中的子网G*xt删除，将直 接连接到子网G*xt的边的终点改为连接到子网G*xi ;其中，边的终点是指：当子网G*xt的上 一层相邻节点与子网G*xt存在边的连接关系时，该条边与子网G*xt的交叉点即为边的终 占 .
[0026] 然后执行步骤9;
[0027] 步骤8.7，
[0028]向树模型中再增加一条从节点Gj连接到子网G*xt的边，并且，所增加边的序号为 xi;
[0029] 然后执行步骤9;
[0030] 步骤9，令i = i+Ι，返回步骤8.2，直到i = k+Ι时，跳出循环，令j = j+Ι，返回步骤7， 直到j = x+l时，跳出循环，令P=P+1，返回步骤5;
[0031] 步骤10,将当前树模型最底层的各个节点直接连接到汇节点T;
[0032] 至此，将所生成的树模型称为简化图模型；
[0033] 步骤11，对于所述简化图模型，最底层每个节点自身的子网连通率为已知值，简化 图模型中各个节点之间边的连通率为已知值，按递归公式计算，计算得到网络G的可靠度近 似值：
[0034] 步骤12,判断计算得到的网络G的可靠度近似值的精度是否符合要求，若不符合， 改变给定截断深度极大值N的值，重复执行步骤3到步骤12，如此不断循环，直到精度满足要 求。
[0035]优选的，步骤2中，采用广度优先排列策略，对边集合中的各条边进行排序。
[0036]本发明的有益效果如下：
[0037]本发明提供的基于截断边扩展图的工程网络可靠度近似分析方法具有以下优点： [0038]经实验证明，能够在生成较小边扩展图和等价BDD的基础上得到误差较小的近似 值，而且，以牺牲最小可靠性精度为代价，最大程度简化了工程网络可靠度的分析繁琐性， 可适用于对中大型网络进行可靠性分析。
【附图说明】
[0039] 图1为给定网络G的具体示例图；
[0040] 图2为对图1进彳丁边排序后的网络图；
[0041] 图3为不考虑同构子图识别以及截断深度问题时得到的扩展图；
[0042]图4为N为4时得到的边扩展图；
[0043]图5为N为3时得到的边扩展图；
[0044]图6为河南省电力系统网络图；
[0045] 图7是北京市轨道交通网络图；
[0046] 图8是大学与产业的关系分配网络图。
【具体实施方式】
[0047] 为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白，以下结合 附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以 解释本发明，并不用于限定本发明。
[0048] 本发明提供一种基于截断边扩展图的工程网络可靠度近似分析方法，包括以下步 骤：
[0049]步骤1，对于给定网络G = (V，E，S，T);其中，V为节点集合，E为边集合;S为源节点，T 为汇节点；
[0050]步骤2,当需要计算源节点S到汇节点T之间的网络可靠度值时，首先遍历给定网络 G，对边集合中的各条边进行排序；
[0051 ]本步骤中，可采用现有技术中已知的广度优先边排序策略对边进行排序，当然，也 可以采用其他排序方法，本发明对此并不限制。
[0052]参考图1，为给定网络G=(V，E，S，T)的一个具体示例；参考图2,为对图1进行边排 序后的网络图。
[0053]步骤3,给定截断深度极大值N;其中，N为自然数;给定截断深度变量值P;建立初始 为空的记录表；
[0054] 初始建立仅存在根节点、并且根节点为给定网络G的树模型;此时，给定网络G为第 0层节点；
[0055] 步骤4,设截断深度变量值P = l;
[0056] 步骤5,判断截断深度变量值P是否小于截断深度极大值N，如果小于，则执行步骤 6;如果大于等于，则执行步骤10;
[0057]步骤6,当前时刻，设树模型第P-1层共存在X个节点，分别记为Gl、G2'"Gx;
[0058] 步骤7,令j = l;
[0059] 步骤8,树模型扩展出第P层节点的方法为：
[0060] 步骤8.1，对于树模型第P-ι层节点中的任意节点G j，设网络G j中与源节点S直接连 接的边为xi(i = 1，2，...，k)，其中，k为自然数，xi代表边的排序序号值；
[0061] 步骤8.2，首先，令1 = 1;
[0062] 步骤8.3,对于网络Gj，将源节点S收缩到边xi的另一端点并作为新的源节点，同 时，删除所有与源节点S直接连接的边、度为1的非S节点、与度为1的非S节点直接相连的边， 得到网络Gj的子网G*xi，其中，G*xi为得到的子网标识；同时，判断所述记录表中是否存在 与子网G*xi的网络结构完全相同的子网，如果不存在，则执行步骤8.4;如果存在，则执行步 骤8.5;
[0063] 步骤8.4,将新得到的子网G*xi的子网网络结构、子网标识以及子网在给定网络G 的深度值的信息记录到所述记录表;其中，子网在给定网络G的深度值是指:从给定网络G到 子网G*x i所经过的边的数量值；
[0064]同时，向树模型中增加子网G*xi，并连接网络Gj节点和子网G*xi，子网G*xi形成为 网络G j的一个儿子节点，子网G*xi的深度值即为P;
[0065]然后，执行步骤9;
[0066] 步骤8.5，假设记录表中已存在子网G*xt，其中，xt关xi ;此时，子网G*xt和子网G* xi的网络结构完全相同，子网G*xt和子网G*xi为同构子图；
[0067] 设子网G*xt的深度值为Nxt，子网G*xi的深度值为Nxi，判断Nxi是否小于N xt，如果小 于，则执行步骤8.6;否则，执行步骤8.7;
[0068] 步骤8.6,将记录表中已记录的子网G*xt的子网网络结构、子网标识以及子网在给 定网络G的深度值的信息删除；同时，将新得到的子网G*xi的子网网络结构、子网标识以及 子网在给定网络G的深度值的信息新增加到记录表；
[0069]同时，向树模型中增加子网G*xi，并连接网络Gj节点和子网G*xi，子网G*xi形成为 网络Gj的一个儿子节点，子网G*xi的深度值即为P;此外，将树模型中的子网G*xt删除，将直 接连接到子网G*xt的边的终点改为连接到子网G*xi ;其中，边的终点是指：当子网G*xt的上 一层相邻节点与子网G*xt存在边的连接关系时，该条边与子网G*xt的交叉点即为边的终 占 .
[0070] 然后执行