CN105490836B - 一种复杂网络可靠度的蒙特卡罗评估方法 - Google Patents
一种复杂网络可靠度的蒙特卡罗评估方法 Download PDFInfo
- Publication number
- CN105490836B CN105490836B CN201510823631.8A CN201510823631A CN105490836B CN 105490836 B CN105490836 B CN 105490836B CN 201510823631 A CN201510823631 A CN 201510823631A CN 105490836 B CN105490836 B CN 105490836B
- Authority
- CN
- China
- Prior art keywords
- event
- network
- link
- reliability
- failure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/08—Configuration management of networks or network elements
- H04L41/0803—Configuration setting
- H04L41/0823—Configuration setting characterised by the purposes of a change of settings, e.g. optimising configuration for enhancing reliability
- H04L41/0836—Configuration setting characterised by the purposes of a change of settings, e.g. optimising configuration for enhancing reliability to enhance reliability, e.g. reduce downtime
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L43/00—Arrangements for monitoring or testing data switching networks
- H04L43/08—Monitoring or testing based on specific metrics, e.g. QoS, energy consumption or environmental parameters
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Environmental & Geological Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
本发明公开了一种复杂网络可靠度的蒙特卡罗评估方法,实现该方法主要包括以下步骤:首先构建整个网络的无向、无自环图,创建事件驱动模型,对链路状态随机取样K次,为每条链路生成故障事件时刻,并向事件表中插入故障事件,然后对事件表判断,如果为空,求出网络可靠度估计值,否则,置所有链路状态是连通的,获得最优先的时间指针,当随机产生链路故障事件的时间指针值为最优先的时间指针时,则链路是故障状态,并在链路事件表中删除故障事件,获得网络状态,若网络状态是连通状态,则自动累计其次数,并返回对事件表进行判断。本发明提出的复杂网络可靠度评估方法不但比其他算法更具优越性和更高的仿真效率,而且在保持估计精度的前提下,大幅降低了复杂度。
Description
技术领域
本发明属于计算机技术领域,涉及网络可靠度的计算问题,尤其涉及一种网络可靠度的蒙特卡罗评估方法。
背景技术
随着网络的高速发展,网络已成为人类生产、生活不可或缺的一部分,这些复杂网络包括医疗网络、教育网络、电网、输油输气管道网络、无线传感器网络等有线网络和无线网络。随着网络规模的不断扩大,网络的结构越来越复杂。网络的设计可能需要随时根据网络可靠度的计算结果对网络进行调整和优化,如增删一些链路、结点或者改变一些结点的位置,如果缺乏可靠性的准确评估,工程人员设计出的网络将面临很多安全隐患,给网络管理、维护和修复带来很多不确定性,甚至一旦发生故障,将造成重大甚至灾难性影响。网络可靠性是保证网络安全因素的重要方面,尤其是近年来,网络安全问题日益成为国际社会关注焦点,针对我国现状急需加强自身网络建设和可靠性评估,提高基础网络的安全性。网络可靠性评估也是度量网络安全程度的重要手段,通过加强复杂网络可靠性研究能有效提升网络建设和管理能力,有利于提高网络安全系数,有利于提高网络的生存性和抗毁性,有利于增强复杂网络可靠性,减少网络的故障概率以及网络资源浪费,因此针对复杂网络可靠性评估具有极其重要的现实意义。
大型复杂网络系统可靠度精确计算是一个NP-hard问题,如何对其进行快速、精确的计算是可靠度研究的重点与热点。从网络可靠性计算精确性角度分析,现有的网络可靠度计算方法可分为两类:精确计算和近似计算。精确计算方法主要包括因子分解、二元决策图,不交最小路方法、不交最小割方法、排序二分决策图、状态枚举法和容斥原理法等方法,这些精确可靠度求解算法随着网络的增大其复杂程度呈指数增长趋势,这类算法往往大多只能适用于中、小型网络的可靠度分析,对于大型复杂工程网络系统,会因网络复杂性而无法求解或者计算量快速膨胀,算法的计算效率变得低下;另一类是近似估计与抽样仿真的方法,近似计算方法通过计算可靠度边界值来逼近精确值,尤其随着网络规模越来越大的客观现实,采用近似估计的方法就受到重视和青睐,该类方法主要有可靠度上下界、蒙特卡罗(MC,Monte Carlo)法等方法,上下界法在确定较理想的可靠度界值范围时,而且往往以牺牲计算复杂度为代价进行权衡;蒙特卡罗法是以概率统计实验为基础的仿真方法,能够真实地模拟实际物理过程,解决一些系统过于复杂而难以建立精确数学模型的问题,具有直观、复杂度低的优点,广泛应用于大型网络的仿真研究。
发明内容
本发明的目的在于提供一种复杂网络可靠度的评估方法,该方法有利于减少时间复杂度和提高计算效率。
为了实现上述目的,本发明提出了一种复杂网络可靠度的蒙特卡罗评估方法,该方法采用事件驱动的蒙特卡罗评估法EMC(Event-driven MC),为每条链路生成故障事件时刻,通过寻找最优先的故障时刻和对应事件更新网络状态,进而实现复杂网络的可靠性评估。
本发明提出的一种网络可靠度的蒙特卡罗评估方法,具体包括如下步骤:
步骤S1:构建整个网络的无向、无自环图,统计网络中的所有节点和链路(边),为每个节点顺序编号为vi,设节点总数为N,1≤i≤N;为每条边顺序编号为ei,链路总数为M,1≤i≤M。
步骤S2:创建事件驱动模型。
步骤S3:基于事件驱动的概率分布取样。
步骤S4:为每条链路生成故障事件时刻,为所有链路ei生成Yi个TPei,并在链路事件表中插入链路ei事件。
步骤S5:判断事件表是否为空,如果为空,执行步骤S9,否则执行步骤S6。
步骤S6:置所有链路ei状态是连通的,寻找并返回事件表中最优先的时间指针TPmin。
步骤S7:当随机产生链路ei故障事件的时间指针值为TPmin,则链路ei是故障状态,并在链路事件表中删除链路ei故障事件。
步骤S8:获得网络状态X,若X状态是网络连通,则累计其次数,并返回步骤S5。
步骤S9:求出网络可靠度估计值。
本发明针对复杂网络结构可靠度的评价问题,提出了一种基于网络故障事件驱动的蒙特卡罗评估方法。该方法通过产生各链路故障事件的时间指针,来构建网络的事件表。通过更新事件表驱动网络状态的更新,在保持了估计精度的前提下,大幅降低了复杂度,是处理大型网络和可靠度极值问题的有效仿真方法。
附图说明
图1为本发明的一种复杂网络可靠度的蒙特卡罗评估方法流程图。
图2为本发明方法的EMC与CMC取样示意图。
图3为本发明所采用的三种25节点栅格网络。
图4为本发明与其它方法可靠性评估精度对比。
图5为本发明与其它方法可靠性评估消耗时间对比。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图及具体实施例对本发明实施方式作进一步地详细描述。
本发明是一种复杂网络可靠度的蒙特卡罗评估方法,针对复杂网络可靠度的计算问题,采用事件驱动的EMC方法,为每条链路生成故障事件时刻,通过寻找最优先的故障时刻和对应事件更新网络状态,进而驱动仿真,请参阅图1所示,本发明包括以下步骤:
步骤S1:构建无向、无自环图G=(V,E)表示网络,V={v1,v2,…,vi,…vN}为网络中N个节点的集合,E={e1,e2,…,ei,…eM}为M条链路(边)的集合。
其中对无向、无自环图G=(V,E)网络做如下假设:
(1)网络为不可修网络;
(2)节点vi完全可靠;
(3)链路ei存在“连通”和“故障”两种状态,分别表示为x(ei)=1、x(ei)=0。
(4)各链路状态相互独立。
(5)链路ei处于x(ei)=1“连通”状态的概率为p(ei),即p(ei)为ei的连通可靠度;“故障”状态的概率为q(ei)=1-p(ei)。
网络状态向量记为X={x(ei)|ei∈E},可得网络处于状态X的概率为(1)式:
定义1网络可靠度为网络处于连通状态的概率,用R(G)表示。
定义2网络结构函数为Ψ(X),若X状态是网络连通,记为事件Ψ(X)=1,反之为Ψ(X)=0;并设X的状态空间为S,共2M个元素,它是网络所有可能状态的全集。可得网络可靠度为(2)式:
步骤S2:创建事件驱动模型,利于事件驱动的方法对复杂网络的可靠度进行蒙特卡罗仿真,为了实现事件驱动的方法,具体包括如下四个部分:
(1)事件表(ET,Event Table);
(2)用于描述各链路故障事件的时间指针(TP,Time Point);
(3)事件更新,用于更新网络状态和事件表;
(4)时间更新,用于更新最优先事件的时间指针。
基于上述设计和分析,定义如下事件表操作:
1)Insert(ET,ei,TPei):向ET中插入参数为(ei,TPei)的故障事件,事件更新。
2)Remove(ET,ei,TPei):删除ET中参数为(ei,TPei)的故障事件,事件更新。
3)Get-TPmin(ET):寻找并返回ET中最优先的时间指针TPmin,时间更新。
步骤S3:事件驱动的随机分布取样。
若对ei的链路状态x(ei)取样K次,用Yi表示链路故障事件的次数,即x(ei)=0的次数。由于p(ei)恒定,故Yi是服从Binomial分布的随机变量:
取样具体分为两步:
(1)依据K和p(ei)生成Binomial变量Yi;
(2)在[0,K]区间内随机产生Yi个ei故障事件的时间指针TPei。
对于链路ei,传统的蒙特卡罗(CMC,CrudeMC)通过K个随机变量U构造长度为K的状态序列X(ei)。EMC与CMC不同在于,EMC直接通过Yi的分布构造故障事件的时间指针序列TPei,即隐含确定了X(ei),如图2所示。
步骤S4:初始化K,k=0,为所有链路ei生成Yi个TPei,并Insert(ET,ei,TPei),寻找最优先的故障时刻和对应事件,从而更新网络状态。
其中,K是ei的链路状态x(ei)取样K次,k是X网络状态是连通时的累计次数。
步骤S5:判断事件表是否为空,如果为空,执行步骤S9,否则执行步骤S6。
步骤S6:置所有链路ei状态是连通的,即x(ei)=1,TPmin=Get-TPmin(ET),寻找并返回事件表中最优先的时间指针TPmin。
步骤S7:随机产生链路ei故障事件,若TPei=TPmin,则链路ei是故障状态,置x(ei)=0,并Remove(ET,ei,TPei),在链路事件表中删除链路ei故障事件。
步骤S8:获得网络状态X,若Ψ(X)=1,则k++;返回步骤S6。
步骤S9:求出网络可靠度估计值
为验证本发明的复杂网络可靠性评估准确性和效率,以25节点的栅格网络为例进行仿真实验,以比较各种算法的优劣。这种栅格网络生成三种栅格网络,链路总数分别为40、56、72,网络复杂度逐渐增大,如图3所示,网络A、B、C。
在相同硬件环境下,分别采用CMC和EMC方法对网络A进行可靠度仿真。参数:p(ei)=0.99;ei∈E;K=106。分别记录三种方法的估计结果和消耗时间,每种方法共仿真30次,并根据30次统计值计算出各方法的估计标准差,结果如图4和图5所示。估计标准差如公式(4):
其中,Var(R)是方差,公式如(5)所示。
依据因子定理,参考文献(赵彦,张新锋,徐国华.因子定理在计算机集成制造系统网络可靠性分析中的应用[J].计算机集成制造系统,2005,11(12):1621-1625.)可计算网络A可靠度的精确值为R(G)=0.99958。实验结果如图4所示,三种方法的估计均值均收敛于此精确值;同时EMC的标准差要略小。实验结果如图5所示,EMC算法不但消耗时间大幅小于CMC,而且仿真效率显著提高。
为比较各方法效率和网络规模、链路状态之间的关系,采用上述方法分别对图3中的网络A、B、C和p(ei)=0.9、0.99、0.999进行九次仿真,实验结果如表1所示。
表1不同网络状态下的仿真比较
由表1分析可知,当网络复杂度相同时,随着链路可靠度的增加,四种方法的消耗时间均逐步减小,但EMC的减幅要远大于CMC、文献[1](ADBULLAH K.Combining NetworkReductions and Simulation to Estimate Network Reliability[C].Proceedings ofthe 2007Winter Simulation Conference,Washington D.C.USA,Association forComputing Machinery,2007:2301-2305.)和文献[2](MANZI E,LABBE M,et al.Fishman’sSampling Plan for Computing Network Reliability[J]IEEE Trans.Reliability,2001,50(1):41-46.)的方法。这是由于在相同的样本容量下,EMC仿真需要产生更少的随机种子;同时当可靠度较高时,EMC跳过了大量的“连通”事件,获得了更高的效率。同时,由于文献[1]的缩减预处理存在适用局限性,对图3网络拓扑的简化作用有限,导致其时间消耗与CMC基本相当(网络A与B略小);文献[2]在取样之前需要进行网络割集计算,时间消耗最大,且对链路可靠度变化不敏感,并验证了文献[2]取样预处理具有较高的复杂度。
当链路可靠度相同时,随着网络复杂度的增加,五种方法的消耗时间逐步增大,同时估计标准差具有逐步减小的趋势。网络复杂度增大,各方法需要处理的状态向量维数将增大,时间消耗也相应增大。由CMC方法的如下(5)式可以发现,当链路可靠度较大时,p(ei)↑,R(G)↑,StdCMC↓,与表1仿真结果相吻合。
当链路可靠度与网络复杂度均相同时,在这四种方法中,文献[2]估计精度最高,利用割集的预处理,将问题转化为割集的分部仿真,压缩了取样维度,具有较低的估计方差。但是,文献[2]估计精度的提升幅度较小,而时间消耗则比EMC高出近一个数量级。这种通过牺牲效率换取精度提升的策略,需要结合实际应用进行具体权衡。
以上是本发明的较佳实施例,凡在本发明的精神和原则之内,所作的任何修改、等同替换等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种复杂网络可靠度的蒙特卡罗评估方法,其特征在于,该方法包括如下步骤:
步骤S1:构建无向、无自环图G=(V,E)表示网络,统计网络中的所有节点和链路或边,为每个节点顺序编号为vi,设节点总数为N,V={v1,v2,…,vi,…vN}为网络中N个节点的集合,1≤i≤N,为每条边顺序编号为ei,链路总数为M,E={e1,e2,…,ei,…eM}为M条链路或边的集合,1≤i≤M;
步骤S2:创建事件驱动模型;
所述事件驱动模型,包括如下四个部分:
(1)事件表ET,Event Table;
(2)用于描述各链路故障事件的时间指针TP,Time Point;
(3)事件更新,用于更新网络状态和事件表;
(4)时间更新,用于更新最优先事件的时间指针;
所述事件表,事件表的操作定义如下:
1)Insert(ET,ei,TPei):向ET中插入参数为(ei,TPei)的故障事件,事件更新;
2)Remove(ET,ei,TPei):删除ET中参数为(ei,TPei)的故障事件,事件更新;
3)Get-TPmin(ET):寻找并返回ET中最优先的时间指针TPmin,时间更新;
步骤S3:事件驱动的随机分布取样;
对ei的链路状态x(ei)取样K次,用Yi表示链路故障事件的次数,即x(ei)=0的次数;由于p(ei)恒定,故Yi是服从Binomial分布的随机变量:
其中,p(ei)为链路ei处于x(ei)=1“连通”状态的概率,即p(ei)为ei的连通可靠度;q(ei)=1-p(ei)为“故障”状态的概率;
所述的事件驱动的随机分布取样,具体按如下步骤进行:
(3a)依据K和p(ei)生成Binomial变量Yi;
(3b)在[0,K]区间内随机产生Yi个ei故障事件的时间指针TPei;
步骤S4:为所有链路生成故障事件时刻,并向事件表中插入故障事件;
步骤S5:判断事件表是否为空,如果为空,执行步骤S9,否则执行步骤S6;
步骤S6:置所有链路状态是连通的,获得最优先的时间指针;
步骤S7:将故障链路事件从事件表中删除;
步骤S8:统计网络状态X连通状态的次数k,并返回步骤S6;
步骤S9:求出网络可靠度估计值
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510823631.8A CN105490836B (zh) | 2015-11-19 | 2015-11-19 | 一种复杂网络可靠度的蒙特卡罗评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510823631.8A CN105490836B (zh) | 2015-11-19 | 2015-11-19 | 一种复杂网络可靠度的蒙特卡罗评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105490836A CN105490836A (zh) | 2016-04-13 |
CN105490836B true CN105490836B (zh) | 2019-04-19 |
Family
ID=55677580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510823631.8A Active CN105490836B (zh) | 2015-11-19 | 2015-11-19 | 一种复杂网络可靠度的蒙特卡罗评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105490836B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107018011B (zh) * | 2017-03-08 | 2020-10-20 | 北京航空航天大学 | 一种保持网络性能可靠度的网络简化方法 |
CN108037659B (zh) * | 2017-11-07 | 2020-10-16 | 西安电子科技大学 | 基于事件驱动的时变耦合复杂动态网络同步方法 |
CN112464422B (zh) * | 2020-11-30 | 2022-09-13 | 西北工业大学 | 一种基于贪婪算法的提高k-终端网络可靠性的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009232144A (ja) * | 2008-03-24 | 2009-10-08 | Kddi Corp | 故障推定装置 |
CN103501235A (zh) * | 2013-07-15 | 2014-01-08 | 中国航天标准化研究所 | 基于元胞自动机的复杂系统可用性确定方法 |
CN103684827A (zh) * | 2012-09-18 | 2014-03-26 | 中兴通讯股份有限公司 | 基于业务的通信网络评估方法和装置 |
CN104468196A (zh) * | 2014-11-13 | 2015-03-25 | 北京邮电大学 | 基于证据筛选的虚拟网络故障诊断方法及装置 |
CN104506354A (zh) * | 2014-12-24 | 2015-04-08 | 云南电力调度控制中心 | 一种基于半实物仿真的数据网风险评估方法 |
-
2015
- 2015-11-19 CN CN201510823631.8A patent/CN105490836B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009232144A (ja) * | 2008-03-24 | 2009-10-08 | Kddi Corp | 故障推定装置 |
CN103684827A (zh) * | 2012-09-18 | 2014-03-26 | 中兴通讯股份有限公司 | 基于业务的通信网络评估方法和装置 |
CN103501235A (zh) * | 2013-07-15 | 2014-01-08 | 中国航天标准化研究所 | 基于元胞自动机的复杂系统可用性确定方法 |
CN104468196A (zh) * | 2014-11-13 | 2015-03-25 | 北京邮电大学 | 基于证据筛选的虚拟网络故障诊断方法及装置 |
CN104506354A (zh) * | 2014-12-24 | 2015-04-08 | 云南电力调度控制中心 | 一种基于半实物仿真的数据网风险评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105490836A (zh) | 2016-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Automatic partitioning of water distribution networks using multiscale community detection and multiobjective optimization | |
Delelegn et al. | Multi-objective optimisation of cost–benefit of urban flood management using a 1D2D coupled model | |
CN104715282A (zh) | 一种基于改进pso-bp神经网络的数据预测方法 | |
CN105490836B (zh) | 一种复杂网络可靠度的蒙特卡罗评估方法 | |
CN104200045A (zh) | 一种流域大尺度水系分布式水动力模型的并行计算方法 | |
CN112491096B (zh) | 一种用于生成电网仿真分析算例的方法及系统 | |
CN104539601B (zh) | 动态网络攻击过程可靠性分析方法及系统 | |
CN102708327A (zh) | 一种基于谱优化的网络社区发现方法 | |
CN113807520A (zh) | 基于图神经网络的知识图谱对齐模型的训练方法 | |
Zhang | Global behavior of a computer virus propagation model on multilayer networks | |
CN111709350B (zh) | 基于fcm聚类的低频振荡模态参数识别方法及系统 | |
CN104579815A (zh) | 网络防护策略的概率模型检测方法 | |
CN109858616A (zh) | 基于神经网络的功率放大器行为级建模系统及方法 | |
CN104794332B (zh) | 一种高层建筑风致响应分析模型的不确定性分析方法 | |
CN106355091B (zh) | 基于生物智能的传播源定位方法 | |
CN104950261B (zh) | 电池的硬件在环仿真测试方法及系统 | |
Lu et al. | Flow simulation system based on high order space-time extension of flux reconstruction method | |
CN105302979B (zh) | 两相流体网络模型中阀门组的建模方法和系统 | |
CN110674581B (zh) | 一种数字孪生模型一致性精准判定方法及系统 | |
Tong et al. | Finite iterative algorithm for solving coupled Lyapunov equations appearing in discrete-time Markov jump linear systems | |
Buccafusca et al. | Maximizing power in wind turbine arrays with variable wind dynamics | |
Lu et al. | High order large eddy simulations for a transonic turbine blade using hybrid unstructured meshes | |
CN104037756B (zh) | 一种含复杂电力设备模型的电力系统稳定评估方法 | |
CN106934729A (zh) | 建筑检测鉴定方法及装置 | |
Amrit et al. | Efficient multi-objective aerodynamic optimization by design space dimension reduction and co-kriging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |