CN115865740B - 一种基于网络结构的关键链路识别方法及装置 - Google Patents
一种基于网络结构的关键链路识别方法及装置 Download PDFInfo
- Publication number
- CN115865740B CN115865740B CN202310186777.0A CN202310186777A CN115865740B CN 115865740 B CN115865740 B CN 115865740B CN 202310186777 A CN202310186777 A CN 202310186777A CN 115865740 B CN115865740 B CN 115865740B
- Authority
- CN
- China
- Prior art keywords
- network
- edge
- matrix
- edges
- representing
- 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
- 238000000034 method Methods 0.000 title claims abstract description 126
- 238000013145 classification model Methods 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000010606 normalization Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 110
- 238000004364 calculation method Methods 0.000 claims description 41
- FESBVLZDDCQLFY-UHFFFAOYSA-N sete Chemical compound [Te]=[Se] FESBVLZDDCQLFY-UHFFFAOYSA-N 0.000 claims description 39
- 240000005523 Peganum harmala Species 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 5
- 238000012512 characterization method Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000010187 selection method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
本发明公开了一种基于网络结构的关键链路识别方法及装置。该方法包括:遍历网络的所有边,计算所述网络中任一边的归一化介数值,得到边归一化介数值集;根据预设的边分类模型,将所述网络中的所有边进行分类,得到第一关键链路集和待选边集;根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集中的边进行处理,得到第二关键链路集;将所述第一关键链路集和所述第二关键链路集进行合并,得到关键链路集。可见,本发明通过对网络中的边进行分类,解决了仅靠介数值确定关键链路的片面性问题,通过设置合理的阈值参数,结合边的介数值确定其余关键链路,解决了介数法中科学性的问题。
Description
技术领域
本发明涉及网络技术领域,尤其涉及一种基于网络结构的关键链路识别方法及装置。
背景技术
关键链路识别是网络技术领域的基础性问题,其目标是发现网络中的关键链路,期望以最小的关键链路集最大程度上实现网络功能的正常运行。以路由网络为例,网络中的路由器为网络节点,路由器之间的链路为网络中的边。给定路由网络,承担路由器之间信息交互的链路的重要性不尽相同,由于在网络中的位置不同,有的链路承载了网络中的大量信息流量,有的链路仅承载了很少一部分的信息流量。为了维持整个路由网路的正常运行,系统维护人员需要确定哪些链路在网络中起着关键核心作用。基于网络中链路的重要性差异,用户可以在系统维护过程中对关键链路进行针对性的安全加固,从而最大程度地保障路由网络的健壮运行。考虑到路由网络中链路数目较多,对多数链路同时加固往往又面临资源不足的问题,因此,系统维护人员总希望找到尽量少的链路进行维护同时保障路由网络最大程度上正常运行。这也就是关键链路识别需要解决的问题,发现使得整个网络系统最大程度正常运行的最少的链路集,即关键链路集。
传统的关键链路识别方法往往基于单一指标方法来识别,比如基于H-index的链路预测方法、基于网络演化模型的关键链路识别方法,这些方法中对于关键链路集规模的估计常常采用单指标排序,取其中排序靠前的一定数量的链路。在实际应用中,仅依据单指标排序进行关键链路的识别,存在片面性与精准性问题。
发明内容
鉴于上述网络中关键链路识别方法存在的问题,本发明提出一种基于网络结构的关键链路识别方法及装置,所述方法中提出了综合利用两种指标识别关键链路,分别对应于第一关键链路集和第二关键链路集。第一种指标重点关注于网络全局的连通性,第二种指标关注链路在网络中的“中介作用”,即介数法,基于介数中心性恒等式给出了关键链路集构建的阈值选取方法,为关键链路识别问题的解决提供了科学性与精准性俱佳的解决方案。
为了解决上述技术问题,本发明实施例第一方面公开了一种基于网络结构的关键链路识别方法,所述方法包括:
遍历网络G的所有边,计算所述网络G中任一边的归一化介数值,得到边归一化介数值集;所述边表征所述网络G中节点之间的链路;
根据预设的边分类模型,将所述网络G的中的所有边进行分类,得到第一关键链路集E 1 和待选边集E 2 ;
根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集E 2 中的边进行处理,得到第二关键链路集E 3 ;
将所述第一关键链路集E 1 和所述第二关键链路集E 3 进行合并,得到关键链路集E。
作为一种可选的实施方式,在本发明实施例第一方面中,所述计算所述网络G中任一边的归一化介数值,包括:
对于网络G中的任一边,设待计算归一化介数值的边为e;
计算所述边e的介数中心值;介数中心值计算公式为:
式中,Bet(e)表示边e的介数中心值,V表示网络G中节点集合,v i ,v j 表示V中任意两个节点,σ(v i, v j│ e)表示节点v i 到节点v j 经过边e的最短路径的数目,σ(v i, v j )表示节点v i 到节点v j 的最短路径的数目;
对所述边e的介数中心值进行归一化处理,得到边e的归一化介数值;所述归一化介数值计算公式为:
式中,bet(e)表示边e的归一化介数值,Bet(e)表示边e的介数中心值,n表示网络中含有的节点的总数。
作为一种可选的实施方式,在本发明实施例第一方面中,所述预设的边分类模型,包括:
设待判断的边为e;
构建网络G的邻接矩阵A=(a i,j );所述a i,j 表示所述邻接矩阵A的第i行第j列元素,若所述网络G中节点i和j之间有边相连,则a i,j =1,否则,
a i,j =0;
将所述网络G中所述边e删除,得到网络G 2 ,构建所述网络G 2 的邻接矩阵A e ;
计算所述邻接矩阵A的n-1次幂,得到矩阵B;所述n表征网络G中节点V的个数;所述B=A n-1 ;
计算所述矩阵B中非零分块矩阵的个数,得到第一矩阵数b;
计算所述邻接矩阵A e 的n-1次幂,得到矩阵B e ;所述B e =A e n-1 ;
计算所述矩阵B e 中非零分块矩阵的个数,得到第二矩阵数b e ;
若所述第一矩阵数b与所述第二矩阵数b e 不相等,则所述边e为关键链路候选边,将所述边e加入第一关键链路集E 1 ,否则,所述边e为待选边,将所述边e加入待选边集E 2 。
作为一种可选的实施方式,在本发明实施例第一方面中,所述根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集E 2 中的边进行处理,得到第二关键链路集E 3 ,包括:
利用预设的阈值参数计算模型,得到阈值参数t;
计算所述网络G的平均最短路径长度l;
计算所述第一关键链路集E 1 中所有边的归一化介数值的和,得到S;
基于所述边归一化介数值集,按边归一化介数值从大到小顺序,依次将所述待选边集E 2 内的待选边输入阈值判断模型,若模型不等式成立,则将所述待选边加入第二关键链路集E 3 ,若模型不等式不成立,则停止输入,得到最终第二关键链路集E 3 ;
所述阈值判断模型为:
式中,E 3 表示所述第二关键链路集E 3 ,e i 表示所述第二关键链路集E 3 中第i条边,bet(e i )表示边e i 的介数值,t表示阈值参数,l表示所述网络G的平均最短路径长度;S表示E 1 中所有边的归一化介数值和。
作为一种可选的实施方式,在本发明实施例第一方面中,所述利用预设的阈值参数计算模型,得到阈值参数t,包括:
计算所述网络G的平均度数c,计算公式为:
式中,c表示所述网络G的平均度数,|E|表示所述网络G中所有边的个数,|V|表示所述网络G中所有节点的个数;
构建平均度数为所述平均度数c的随机Erdos-Renyi网络,所述随机Erdos-Renyi网络的度分布为:
式中,e表示自然常数e,c表示网络G的平均度数;
基于所述网络G和所述随机Erdos-Renyi网络的度分布向量的相关度,利用预设的阈值参数计算模型,得到阈值参数t;所述阈值参数t取值范围满足0≤t≤1;所述预设的阈值参数计算模型为:
式中,向量表示网络G中节点的度数分布,所述向量,其中代表网络G中度数为x的节点出现的概率,所述;向量表示所述随机Erdos-Renyi网络中节点的度数分布,所述向量,其中,代表网络G对应的等效随机Erdos-Renyi网络中度数为y的节点出现的概率,所述;表示向量与向量的内积;,分别代表向量和向量的欧几里得范数;k max 表示所述网络G=(V,E)与所述随机Erdos-Renyi网络中节点的最大度数;向量与向量的维数为所述k max 。
作为一种可选的实施方式,在本发明实施例第一方面中,所述计算网络G的平均最短路径长度l,包括:
基于所述最短距离矩阵L,计算所述网络G的平均最短路径长度l;所述平均最短路径长度l的计算公式为:
式中,n表示网络G=(V,E)中节点的个数;l i,j 表示所述最短距离矩阵L中第i行第j列的元素。
本发明实施例第二方面公开了一种基于网络结构的关键链路识别装置,所述装置包括:
边介数值计算模块;用于计算网络中所有边的归一化介数值,得到边归一化介数值集;
边判断模块;用于根据预设的边分类模型,将所述网络中的所有边进行分类,得到第一关键链路集E 1 和待选边集E 2 ;
阈值判断模块;根据预设的阈值判断模型,对所述待选边集中的边进行处理,得到第二关键链路集E 3 ;
关键链路集生成模块;用于将所述第一关键链路集E 1 和所述第二关键链路集E 3 进行合并,得到关键链路集E。
作为一种可选的实施方式,在本发明实施例第二方面中,所述计算网络中所有边的归一化介数值,包括:
对于网络中的任一边,设待计算归一化介数值的边为e;
计算所述边e的介数中心值;所述介数中心值计算公式为:
式中,Bet(e)表示边e介数中心值,V表示网络G中节点集合,v i, v j 表示V中任意两个节点,σ(v i, v j│ e)表示节点v i 到节点v j 经过边e的最短路径的数目,σ(v i, v j )为节点v i 到节点v j 的最短路径的数目;
对所述边e的介数中心值进行归一化处理,得到边e的归一化介数值;所述归一化计算公式为:
式中,bet(e)表示边e的归一化介数值,Bet(e)表示边e的介数中心值,n表示网络中含有的节点的总数。
作为一种可选的实施方式,在本发明实施例第二方面中,所述预设的边分类模型,包括:
设待判断的边为e;
构建网络G的邻接矩阵A=(a i,j );所述a i,j 表示所述邻接矩阵A的第i行第j列元素,若所述网络G中节点i和j之间有边相连,则a i,j =1,否则,a i,j =0;
将所述网络G中所述边e删除,得到网络G 2 ,构建所述网络G 2 的邻接矩阵A e ;
计算所述邻接矩阵A的n-1次幂,得到矩阵B;所述n表征网络G中节点V的个数;所述B=A n-1 ;
计算所述矩阵B中非零分块矩阵的个数,得到第一矩阵数b;
计算所述邻接矩阵A e 的n-1次幂,得到矩阵B e ;所述B e =A e n-1 ;
计算所述矩阵B e 中非零分块矩阵的个数,得到第二矩阵数b e ;
若所述第一矩阵数b与所述第二矩阵数b e 不相等,则所述边e为关键链路候选边,将所述边e加入第一关键链路集E1,否则,所述边e为待选边,将所述边e加入待选边集E 2 。
作为一种可选的实施方式,在本发明实施例第二方面中,所述根据预设的阈值判断模型,对所述待选边集中的边进行处理,得到第二关键链路集E 3 ,包括:
利用预设的阈值参数计算模型,得到阈值参数t;
计算所述网络G的平均最短路径长度l;
计算所述第一关键链路集E 1 中所有边的归一化介数值的和,得到S;
基于所述边归一化介数值集,按边归一化介数值从大到小顺序,依次将所述待选边集E 2 内的待选边输入阈值判断模型,若模型不等式成立,则将所述待选边加入第二关键链路集E 3 ,若模型不等式不成立,则停止输入,得到最终第二关键链路集E 3 ;
所述阈值判断模型为:
式中,E 3 表示所述第二关键链路集E 3 ,e i 表示所述第二关键链路集E 3 中第i条边,bet(e i )表示边e i 的介数值,t表示阈值参数,l表示所述网络G的平均最短路径长度;S表示E 1 中所有边的归一化介数值和。
可选的,所述利用预设的阈值参数计算模型,得到阈值参数t,包括:
计算所述网络G的平均度数c,计算公式为:
式中,c表示所述网络G的平均度数,|E|表示所述网络G中所有边的个数,|V|表示所述网络G中所有节点的个数;
构建平均度数为所述平均度数c的随机Erdos-Renyi网络,所述随机Erdos-Renyi网络的度分布为:
式中,e表示自然常数e,c表示网络G的平均度数;
基于所述网络G和所述随机Erdos-Renyi网络的度分布向量的相关度,利用预设的阈值参数计算模型,得到阈值参数t;所述阈值参数t取值范围满足0≤t≤1;所述预设的阈值参数计算模型为:
式中,向量表示网络G中节点的度数分布,所述向量,其中代表网络G中度数为x的节点出现的概率,所述;向量表示所述随机Erdos-Renyi网络中节点的度数分布,所述向量,其中,代表网络G对应的等效随机Erdos-Renyi网络中度数为y的节点出现的概率,所述;表示向量与向量的内积;,分别代表向量和向量的欧几里得范数;k max 表示所述网络G=(V,E)与所述随机Erdos-Renyi网络中节点的最大度数;向量与向量的维数为所述k max 。
可选的,所述计算网络G的平均最短路径长度l,包括:
基于所述最短距离矩阵L,计算所述网络G的平均最短路径长度l;所述平均最短路径长度l的计算公式为:
式中,n表示网络G=(V,E)中节点的个数;l i,j 表示所述最短距离矩阵L中第i行第j列的元素。
本发明实施例第三方面公开了另一种基于网络结构的关键链路识别装置,所述装置包括:
存储有可执行程序代码的存储器;
与所述存储器耦合的处理器;
所述处理器调用所述存储器中存储的所述可执行程序代码,执行本发明实施例第一方面公开的基于网络结构的关键链路识别方法中的部分或全部步骤。
本发明实施例第四方面公开了又一种计算机存储介质,所述计算机存储介质存储有计算机指令,所述计算机指令被调用时,用于执行本发明实施例第一方面公开的基于网络结构的关键链路识别方法中的部分或全部步骤。
与现有技术相比,本发明实施例具有以下有益效果:
本发明实施例中,遍历网络的所有边,计算所述网络中任一边的归一化介数值,得到边归一化介数值集;所述边表征所述网络中节点之间的链路;根据预设的边分类模型,将所述网络中的所有边进行分类,得到第一关键链路集和待选边集;根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集中的边进行处理,得到第二关键链路集;将所述第一关键链路集和所述第二关键链路集进行合并,得到关键链路集。可见,本发明通过先对网络中的边进行分类,解决了仅靠介数值确定关键链路的片面性问题,通过设置合理的阈值参数,结合边的介数值确定关键链路,解决了介数法中科学性的问题。
附图说明
图1是本发明实施例公开的一种基于网络结构的关键链路识别方法流程示意图;
图2是本发明实施例公开的一种网络结构示意图;
图3是本发明实施例公开的一种基于网络结构的关键链路识别装置结构示意图;
图4是本发明实施例公开的另一种基于网络结构的关键链路识别装置结构示意图;
图5是本发明实施例公开的又一种基于网络结构的关键链路识别装置结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别不同对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、装置、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
在本文中提及“实施例”意味着,结合实施例描述的特定特征、结构或特性可以包含在本发明的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例,也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是,本文所描述的实施例可以与其它实施例相结合。
本发明公开了基于网络结构的关键链路识别方法及装置,首先找出能将连通的网络分隔为不连通的子网络的边,作为第一关键链路集,再利用阈值判断模型,基于边的介数值确定剩余边中的第二关键链路集,将第一关键链路集和第二关键链路集确定为最终关键链路集,提高关键链路提取的全面性和科学性。
实施例一
请参阅图1,图1是本发明实施例公开的一种基于网络结构的关键链路识别方法流程示意图。其中,图1所描述基于网络结构的关键链路识别方法可应用于信息网络、社交网络、生物网络和交通网络,本发明实施例不做限定。如图1所示,该基于网络结构的关键链路识别方法可以包括以下操作:
101、遍历网络的所有边,计算所述网络中任一边的归一化介数值,得到边归一化介数值集。
本发明实施例中,上述边表征所述网络中节点之间的链路。
102、根据预设的边分类模型,将所述网络中的所有边进行分类,得到第一关键链路集和待选边集。
103、根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集中的边进行处理,得到第二关键链路集。
104、将所述第一关键链路集和所述第二关键链路集进行合并,得到关键链路集。
可见,实施本发明实施例所描述的基于网络结构的关键链路识别方法,能够首先找出能将连通的网络分隔为不连通的子网络的边,作为第一关键链路集,再利用阈值判断模型,基于边的介数值确定剩余边中的第二关键链路集,将第一关键链路集和第二关键链路集合并为最终关键链路集,提高关键链路提取的全面性和准确性。
在一个可选的实施例中,上述计算网络中边的归一化介数值,包括:
设待计算归一化介数值的边为e;
计算所述边e的介数中心值;所述介数中心值计算公式为:
式中,Bet(e)表示边e的介数中心值,V表示网络G中节点集合,v i ,v j 表示V中任意两个节点,σ(v i, v j│ e)表示节点v i 到节点v j 经过边e的最短路径的数目,σ(v i, v j )表示节点v i 到节点v j 的最短路径的数目。
对所述边e的介数中心值进行归一化处理,得到边e的归一化介数值;所述归一化计算公式为:
式中,bet(e)表示边e的归一化介数值,Bet(e)表示边e的介数中心值,n表示网络中含有的节点的总数。
在另一个可选的实施例中,上述预设的边分类模型,具体包括:
设待判断的边为e;
构建网络G的邻接矩阵A=(a i,j );所述a i,j 表示所述邻接矩阵A的第i行第j列元素,若所述网络G中节点i和j之间有边相连,则a i,j =1,否则,
a i,j =0;
将所述网络G中所述边e删除,得到网络G 2 ,构建所述网络G 2 的邻接矩阵A e ;
计算所述邻接矩阵A的n-1次幂,得到矩阵B;所述n表征网络G中节点V的个数;所述B=A n-1 ;
计算所述矩阵B中非零分块矩阵的个数,得到第一矩阵数b;
计算所述邻接矩阵A e 的n-1次幂,得到矩阵B e ;所述B e =A e n-1 ;
计算所述矩阵B e 中非零分块矩阵的个数,得到第二矩阵数b e ;
若所述第一矩阵数b与所述第二矩阵数b e 不相等,则所述边e为关键链路候选边,将所述边e加入第一关键链路集E1,否则,所述边e为待选边,将所述边e加入待选边集E 2 。
可选的,如图2中所描述的一种网络结构示意图,所述网络结构中包括节点1、节点2、节点3、节点4、节点5,节点3和节点4之间的边为e,判断边e是否为关键链路的方法如下:
构建所述网络结构的邻接矩阵A,计算所述A的2~(n-1)次幂,所述n=6;
构建删除边e的子网络对应的邻接矩阵A e ,计算所述A e 的2~(n-1)次幂,所述n=6。
从矩阵A 5和的结果上可以看出,矩阵A 5具有1个非零分块矩阵,矩阵具有2个非零分块矩阵。也就是说,网络G删除边e之后,从一个全连通网络变成了两个连通分支构成的网络,连通分支数目增加,所以边e为关键链路候选边,应将所述边e加入第一关键链路集。类似地,可以判断所述网络结构中除边e之外的其余边均不是关键链路候选边,属于待选边,应加入待选边集。
可见,预设的边分类模型重点关注了关系到网络全局的连通性的边,找出了能将连通的网络分隔为不连通的子网络的边,如果删除该边,将导致网络被拆分两个隔离的网络,说明该边对网络连通性至关重要,因此确定该类边作为第一关键链路集。
在又一个可选的实施例中,上述所述根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集中的边进行处理,得到第二关键链路集,包括,具体包括:
利用预设的阈值参数计算模型,得到阈值参数t;
计算所述网络G的平均最短路径长度l;
基于所述边归一化介数值集,按介数值从大到小顺序,依次将所述待选边集E 2 内的待选边输入阈值判断模型,若模型不等式成立,则将所述待选边加入第二关键链路集E 3 ,若模型不等式不成立,则停止输入,得到最终第二关键链路集E 3 ;
所述阈值判断模型为:
式中,E 3 表示所述第二关键链路集E 3 ,e i 表示所述第二关键链路集E 3 中第i条边,bet(e i )表示边e i 的介数值,t表示阈值参数,l表示所述网络G的平均最短路径长度;S表示E 1 中所有边的归一化介数值和;
可见,预设的阈值判断模型对介数法进行了改进,边的介数代表了该链路在所有节点对的交换过程中的重要程度,通过利用阈值判断模型,基于介数中心性恒等式,全面地、科学地选取出第二关键链路集。
可选的,上述利用预设的阈值参数计算模型,得到阈值参数t,具体包括:
计算所述网络G的平均度数c,计算公式为:
式中,c表示所述网络G的平均度数,|E|表示所述网络G中所有边的个数,|V|表示所述网络G中所有节点的个数。
构建平均度数为所述平均度数c的随机Erdos-Renyi网络,所述随机Erdos-Renyi网络的度分布为:
式中,e表示自然常数e,c表示网络G的平均度数;
基于所述网络G和所述随机Erdos-Renyi网络的度分布向量的相关度,利用预设的阈值参数计算模型,得到阈值参数t;所述阈值参数t取值范围满足0≤t≤1;所述预设的阈值参数计算模型为:
式中,向量表示网络G中节点的度数分布,所述向量,其中代表网络G中度数为x的节点出现的概率,所述;向量表示所述随机Erdos-Renyi网络中节点的度数分布,所述向量,其中,代表网络G对应的等效随机Erdos-Renyi网络中度数为y的节点出现的概率,所述;表示向量与向量的内积;,分别代表向量和向量的欧几里得范数;k max 表示所述网络G=(V,E)与所述随机Erdos-Renyi网络中节点的最大度数;向量与向量的维数为所述k max 。
可选的,上述计算网络G的平均最短路径长度l,包括:
基于所述最短距离矩阵L,计算所述网络G的平均最短路径长度l;所述平均最短路径长度l的计算公式为:
式中,n表示网络G=(V,E)中节点的个数;l i,j 表示所述最短距离矩阵L中第i行第j列的元素。
可见,通过基于介数中心性恒等式给出了关键链路集构建的阈值选取方法,为关键链路识别问题的解决提供了科学性与精准性的解决方案。
实施例二
请参阅图3,图3是本发明实施例公开的一种基于网络结构的关键链路识别装置结构示意图。其中,图3所描述基于网络结构的关键链路识别装置可应用于信息网络、社交网络、生物网络和交通网络,本发明实施例不做限定。如图3所示,该装置可以包括:
边介数值计算模块301;用于计算网络中所有边的归一化介数值,得到边归一化介数值集;
边判断模块302;用于根据预设的边分类模型,将所述网络中的所有边进行分类,得到第一关键链路集E 1 和待选边集E 2 ;
阈值判断模块303;用于根据预设的阈值判断模型,对所述待选边集中的边进行处理,得到第二关键链路集E 3 ;
关键链路集生成模块304;用于将所述第一关键链路集E 1 和所述第二关键链路集E 3 进行合并,得到关键链路集E。
可见,实施图3所描述的基于网络结构的关键链路识别装置,首先计算网络中所有边的归一化介数值,得到边介数值集;通过预设的边分类模型,找出能将连通的网络分隔为不连通的子网络的边,作为第一关键链路集,再利用阈值判断模型,基于边的介数值集,确定第二关键链路集,将第一关键链路集和第二关键链路集合并为最终关键链路集,不仅提高关键链路提取的全面性,而且更加科学。
在一个可选的实施例中,上述计算网络中所有边的归一化介数值,包括:
对于网络中的任一边,设待计算归一化介数值的边为e;
计算所述边e的介数中心值;所述介数中心值计算公式为:
式中,Bet(e)表示边e的介数中心值,V表示网络G中节点集合,v i ,v j 表示V中任意两个节点,σ(v i, v j│ e)表示节点v i 到节点v j 经过边e的最短路径的数目,σ(v i, v j )表示节点v i 到节点v j 的最短路径的数目;
对所述边e的介数中心值进行归一化处理,得到边e的归一化介数值;所述归一化计算公式为:
式中,bet(e)表示边e的归一化介数值,Bet(e)表示边e的介数中心值,n表示网络中含有的节点的总数。
在另一个可选的实施例中,上述预设的边分类模型,包括:
设待判断的边为e;
构建网络G的邻接矩阵A=(a i,j );所述a i,j 表示所述邻接矩阵A的第i行第j列元素,若所述网络G中节点i和j之间有边相连,则a i,j =1,否则,
a i,j =0;
将所述网络G中所述边e删除,得到网络G 2 ,构建所述网络G 2 的邻接矩阵A e ;
计算所述邻接矩阵A的n-1次幂,得到矩阵B;所述n表征网络G中节点V的个数;所述B=A n-1 ;
计算所述矩阵B中非零分块矩阵的个数,得到第一矩阵数b;
计算所述邻接矩阵A e 的n-1次幂,得到矩阵B e ;所述B e =A e n-1 ;
计算所述矩阵B e 中非零分块矩阵的个数,得到第二矩阵数b e ;
若所述第一矩阵数b与所述第二矩阵数b e 不相等,则所述边e为关键链路候选边,将所述边e加入第一关键链路集E1,否则,所述边e为待选边,将所述边e加入待选边集E 2 。
在又一个可选的实施例中,上述根据预设的阈值判断模型,对所述待选边集中的边进行处理,得到第二关键链路集E 3 ,包括:
基于所述边归一化介数值集,利用阈值判断模型,从待选边集E 2 内选出关键链路,得到第二关键链路集E 3 。
在又一个可选的实施例中,如图4所示,阈值判断模块303包括阈值参数计算模块3031、平均最短路径计算模块3032和边检验模块3033,其中:
平均最短路径计算模块3032,用于计算网络的平均最短路径l;
边检验模块3033,用于从待选边集E 2 内选出关键链路,得到第二关键链路集E 3 。
在又一个可选的实施例中,上述计算阈值参数t,包括:
计算所述网络G的平均度数c,计算公式为:
式中,c表示所述网络G的平均度数,|E|表示所述网络G中所有边的个数,|V|表示所述网络G中所有节点的个数;
构建平均度数为所述平均度数c的随机Erdos-Renyi网络,所述随机Erdos-Renyi网络的度分布为:
式中,e表示自然常数e,c表示网络G的平均度数;
基于所述网络G和所述随机Erdos-Renyi网络的度分布向量的相关度,利用预设的阈值参数计算模型,得到阈值参数t;所述阈值参数t取值范围满足0≤t≤1;所述预设的阈值参数计算模型为:
式中,向量表示网络G中节点的度数分布,所述向量,其中代表网络G中度数为x的节点出现的概率,所述;向量表示所述随机Erdos-Renyi网络中节点的度数分布,所述向量,其中,代表网络G对应的等效随机Erdos-Renyi网络中度数为y的节点出现的概率,所述;表示向量与向量的内积;,分别代表向量和向量的欧几里得范数;k max 表示所述网络G=(V,E)与所述随机Erdos-Renyi网络中节点的最大度数;向量与向量的维数为所述k max 。
在又一个可选的实施例中,上述计算网络的平均最短路径长度l,包括:
基于所述最短距离矩阵L,计算所述网络G的平均最短路径长度l;所述平均最短路径长度l的计算公式为:
式中,n表示网络G=(V,E)中节点的个数;l i,j 表示所述最短距离矩阵L中第i行第j列的元素。
在又一个可选的实施例中,上述从待选边集E 2 内选出关键链路,得到第二关键链路集E 3 ,包括:
计算所述第一关键链路集E 1 中所有边的介数值的和,得到S;
基于所述边归一化介数值集,按归一化介数值从大到小顺序,依次将所述待选边集E 2 内的待选边输入阈值判断模型,若模型不等式成立,则将所述待选边加入第二关键链路集E 3 ,若模型不等式不成立,则停止输入,得到最终第二关键链路集E 3 。
所述阈值判断模型为:
式中,E 3 表示所述第二关键链路集E 3 ,e i 表示所述第二关键链路集E 3 中第i条边,bet(e i )表示边e i 的介数值,t表示阈值参数,l表示所述网络G的平均最短路径长度;S表示E 1 中所有边的归一化介数值和。
可见,实施本实施例所描述的基于网络结构的关键链路识别装置,基于介数中心性恒等式给出了关键链路集构建的阈值选取方法,确定阈值参数,计算得出网络的平均最短路径长度l,利用阈值判断模型,基于所述边归一化介数值集,从待选边集E 2 内选出第二关键链路集,提高关键链路提取的科学性。
实施例三
请参阅图5,图5是本发明实施例公开的又一种基于网络结构的关键链路识别装置结构示意图。其中,图5所描述的装置可应用于信息网络、社交网络、生物网络和交通网络,本发明实施例不做限定。如图5所示,该装置可以包括:
存储有可执行程序代码的存储器501;
与存储器501耦合的处理器502;
处理器502调用存储器501中存储的可执行程序代码,用于执行实施例一所描述的基于网络结构的关键链路识别方法中的步骤。
实施例四
本发明实施例公开了一种计算机读存储介质,其存储用于电子数据交换的计算机程序,其中,该计算机程序使得计算机执行实施例一所描述的基于网络结构的关键链路识别方法中的步骤。
实施例五
本发明实施例公开了一种计算机程序产品,该计算机程序产品包括存储了计算机程序的非瞬时性计算机可读存储介质,且该计算机程序可操作来使计算机执行实施例一所描述的基于网络结构的关键链路识别方法中的步骤。
最后应说明的是:本发明实施例公开的基于网络结构的关键链路识别方法及装置所揭露的仅为本发明较佳实施例而已,仅用于说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解;其依然可以对前述各项实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或替换,并不使相应的技术方案的本质脱离本发明各项实施例技术方案的精神和范围。
Claims (8)
1.一种基于网络结构的关键链路识别方法,其特征在于,包括:
遍历网络G的所有边,计算所述网络G中任一边的归一化介数值,得到边归一化介数值集;所述边表征所述网络G中节点之间的链路;
根据预设的边分类模型,将所述网络G的中的所有边进行分类,得到第一关键链路集E 1 和待选边集E 2 ;所述待选边集E 2 由网络G中除所述第一关键链路集E 1 之外的边组成;具体的:
设待判断的边为e;
构建网络G的邻接矩阵A=(a i,j );
将所述网络G中所述边e删除,得到网络G 2 ,构建所述网络G 2 的邻接矩阵A e ;
计算所述邻接矩阵A的n-1次幂,得到矩阵B;所述n表征网络G中节点V的个数;所述B=A n -1 ;
计算所述矩阵B中非零分块矩阵的个数,得到第一矩阵数b;
计算所述邻接矩阵A e 的n-1次幂,得到矩阵B e ;所述B e =A e n-1 ;
计算所述矩阵B e 中非零分块矩阵的个数,得到第二矩阵数b e ;
若所述第一矩阵数b与所述第二矩阵数b e 不相等,则所述边e为关键链路候选边,将所述边e加入第一关键链路集E 1 ,否则,所述边e为待选边,将所述边e加入待选边集E 2 ;
根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集E 2 中的边进行处理,得到第二关键链路集E 3 ;
将所述第一关键链路集E 1 和所述第二关键链路集E 3 进行合并,得到关键链路集E。
2.根据权利要求1所述的基于网络结构的关键链路识别方法,其特征在于,所述计算所述网络G中任一边的归一化介数值,包括:
对于网络G中的任一边,设待计算归一化介数值的边为e;
计算所述边e的介数中心值;
对所述边e的介数中心值进行归一化处理,得到边e的归一化介数值。
3.根据权利要求1所述的基于网络结构的关键链路识别方法,其特征在于,所述根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集E 2 中的边进行处理,得到第二关键链路集E 3 ,包括:
利用预设的阈值参数计算模型,得到阈值参数t;
计算所述网络G的平均最短路径长度l;
计算所述第一关键链路集E 1 中所有边的归一化介数值的和,得到S;
基于所述边归一化介数值集,按边归一化介数值从大到小顺序,依次将所述待选边集E 2 内的待选边输入阈值判断模型,若模型不等式成立,则将所述待选边加入第二关键链路集E 3 ,若模型不等式不成立,则停止输入,得到最终第二关键链路集E 3 ;
所述阈值判断模型为:
计算所述网络G的平均度数c;
构建平均度数为所述平均度数c的随机Erdos-Renyi网络,所述随机Erdos-Renyi网络的度分布为:
式中,e表示自然常数e,c表示网络G的平均度数;
基于所述网络G和所述随机Erdos-Renyi网络的度分布向量的相关度,利用预设的阈值参数计算模型,得到阈值参数t;所述阈值参数t取值范围满足0≤t≤1;
所述预设的阈值参数计算模型为:
6.一种基于网络结构的关键链路识别装置,其特征在于,所述装置包括:
边介数值计算模块,用于计算网络中所有边的归一化介数值,得到边归一化介数值集;
边判断模块,用于根据预设的边分类模型,将所述网络中的所有边进行分类,得到第一关键链路集E 1 和待选边集E 2 ;所述待选边集E 2 由网络G中除所述第一关键链路集E 1 之外的边组成;
所述预设的边分类模型为:
设待判断的边为e;
构建网络G的邻接矩阵A=(a i,j );
将所述网络G中所述边e删除,得到网络G 2 ,构建所述网络G 2 的邻接矩阵A e ;
计算所述邻接矩阵A的n-1次幂,得到矩阵B;所述n表征网络G中节点V的个数;所述B=A n -1 ;
计算所述矩阵B中非零分块矩阵的个数,得到第一矩阵数b;
计算所述邻接矩阵A e 的n-1次幂,得到矩阵B e ;所述B e =A e n-1 ;
计算所述矩阵B e 中非零分块矩阵的个数,得到第二矩阵数b e ;
若所述第一矩阵数b与所述第二矩阵数b e 不相等,则所述边e为关键链路候选边,将所述边e加入第一关键链路集E 1 ,否则,所述边e为待选边,将所述边e加入待选边集E 2 ;
阈值判断模块,用于根据边归一化介数值集,利用预设的阈值判断模型,对所述待选边集E 2 中的边进行处理,得到第二关键链路集E 3 ;
关键链路集生成模块,用于将所述第一关键链路集E 1 和所述第二关键链路集E 3 进行合并,得到关键链路集E。
7.一种基于网络结构的关键链路识别装置,其特征在于,所述装置包括:
存储有可执行程序代码的存储器;
与所述存储器耦合的处理器;
所述处理器调用所述存储器中存储的所述可执行程序代码,执行如权利要求1-5任一项所述的基于网络结构的关键链路识别方法。
8.一种计算机存储介质,其特征在于,所述计算机存储介质存储有计算机指令,所述计算机指令被调用时,执行如权利要求1-5任一项所述的基于网络结构的关键链路识别方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310186777.0A CN115865740B (zh) | 2023-03-02 | 2023-03-02 | 一种基于网络结构的关键链路识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310186777.0A CN115865740B (zh) | 2023-03-02 | 2023-03-02 | 一种基于网络结构的关键链路识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115865740A CN115865740A (zh) | 2023-03-28 |
CN115865740B true CN115865740B (zh) | 2023-05-12 |
Family
ID=85659605
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310186777.0A Active CN115865740B (zh) | 2023-03-02 | 2023-03-02 | 一种基于网络结构的关键链路识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115865740B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1530324A1 (en) * | 2003-11-06 | 2005-05-11 | Siemens Aktiengesellschaft | A method and a network node for routing a call through a communication network |
CN110135092A (zh) * | 2019-05-21 | 2019-08-16 | 江苏开放大学(江苏城市职业学院) | 基于半局部中心性的复杂加权交通网关键节点识别方法 |
CN113225215A (zh) * | 2021-05-19 | 2021-08-06 | 华北电力大学 | 一种sdn架构下区分服务网络关键链路识别方法及系统 |
CN113271226A (zh) * | 2021-05-19 | 2021-08-17 | 湖南经研电力设计有限公司 | 电力通信网关键链路识别方法及其保护方法 |
CN115440392A (zh) * | 2022-08-19 | 2022-12-06 | 青海师范大学 | 一种基于删后拉普拉斯矩阵的重要超边识别方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8948053B2 (en) * | 2011-09-12 | 2015-02-03 | Honeywell International Inc. | Apparatus and method for detecting critical nodes and critical links in a multi-hop network |
US8842520B2 (en) * | 2011-09-12 | 2014-09-23 | Honeywell International Inc. | Apparatus and method for identifying optimal node placement to form redundant paths around critical nodes and critical links in a multi-hop network |
CN106027399B (zh) * | 2016-07-26 | 2018-11-13 | 华北电力大学(保定) | 一种通信网关键链路识别方法 |
CN107623319B (zh) * | 2017-08-17 | 2021-03-02 | 广东电网有限责任公司惠州供电局 | 一种基于多评价指标的电网关键线路辨识方法 |
US10716048B2 (en) * | 2018-09-10 | 2020-07-14 | Telefonaktiebolaget Lm Ericsson (Publ) | Detecting critical links in bluetooth mesh networks |
US10536366B1 (en) * | 2018-09-13 | 2020-01-14 | Charter Communication Operating, LLC | Methods and apparatus for controlling and making link bundle advertisements to support routing decisions |
CN109921939B (zh) * | 2019-03-18 | 2022-04-15 | 中电科大数据研究院有限公司 | 一种通信网络中关键节点的选取方法及系统 |
CN110138437B (zh) * | 2019-04-03 | 2021-04-20 | 西安建筑科技大学 | 基于时延性能的卫星通信网络关键链路序列检测方法 |
CN111428323B (zh) * | 2020-04-16 | 2023-06-23 | 太原理工大学 | 复杂网络中利用广义折扣度与k-shell识别一组关键节点的方法 |
CN115189949B (zh) * | 2022-07-11 | 2023-06-16 | 中国电子科技集团公司第十五研究所 | 一种基于全局信息的网络安全控制关键节点识别方法 |
CN115622902B (zh) * | 2022-12-19 | 2023-04-07 | 中国人民解放军国防科技大学 | 基于网络结构和节点价值的电信网节点重要度计算方法 |
-
2023
- 2023-03-02 CN CN202310186777.0A patent/CN115865740B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1530324A1 (en) * | 2003-11-06 | 2005-05-11 | Siemens Aktiengesellschaft | A method and a network node for routing a call through a communication network |
CN110135092A (zh) * | 2019-05-21 | 2019-08-16 | 江苏开放大学(江苏城市职业学院) | 基于半局部中心性的复杂加权交通网关键节点识别方法 |
CN113225215A (zh) * | 2021-05-19 | 2021-08-06 | 华北电力大学 | 一种sdn架构下区分服务网络关键链路识别方法及系统 |
CN113271226A (zh) * | 2021-05-19 | 2021-08-17 | 湖南经研电力设计有限公司 | 电力通信网关键链路识别方法及其保护方法 |
CN115440392A (zh) * | 2022-08-19 | 2022-12-06 | 青海师范大学 | 一种基于删后拉普拉斯矩阵的重要超边识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115865740A (zh) | 2023-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111565205B (zh) | 网络攻击识别方法、装置、计算机设备和存储介质 | |
CN111090807B (zh) | 一种基于知识图谱的用户识别方法及装置 | |
CN115577858B (zh) | 基于区块链的碳排放量预测方法、装置及电子设备 | |
CN111260220A (zh) | 群控设备识别方法、装置、电子设备和存储介质 | |
CN111639230B (zh) | 一种相似视频的筛选方法、装置、设备和存储介质 | |
CN116596095A (zh) | 基于机器学习的碳排放量预测模型的训练方法及装置 | |
CN115830649A (zh) | 一种网络资产指纹特征的识别方法、装置及电子设备 | |
CN111598711A (zh) | 目标用户账号识别方法、计算机设备及存储介质 | |
CN114494753A (zh) | 聚类方法、装置、电子设备及计算机可读存储介质 | |
CN110309154B (zh) | 基于图谱的实体特征选择方法、装置、设备和存储介质 | |
CN114492601A (zh) | 资源分类模型的训练方法、装置、电子设备及存储介质 | |
CN114142923A (zh) | 光缆故障定位方法、装置、设备及可读介质 | |
WO2019184325A1 (zh) | 基于平均互信息的社区划分质量评价方法及系统 | |
CN117807245A (zh) | 网络资产图谱中节点特征提取方法及相似节点搜索方法 | |
CN109492844B (zh) | 业务策略的生成方法和装置 | |
CN117390480A (zh) | 一种信息提取方法、装置、设备及存储介质 | |
CN115865740B (zh) | 一种基于网络结构的关键链路识别方法及装置 | |
CN111784402A (zh) | 基于多通路的下单率预测方法、设备及可读存储介质 | |
CN114978878B (zh) | 定位方法、装置、电子设备及计算机可读存储介质 | |
CN109993338B (zh) | 一种链路预测方法及装置 | |
CN111917589B (zh) | 一种电力通信网络资源备份方法及相关装置 | |
KR102409160B1 (ko) | 복수개의 커뮤니티를 포함하는 네트워크에서 커뮤니티 재구성 방법 및 이를 위한 전자 장치 | |
CN110032596B (zh) | 流量异常用户识别方法及系统 | |
CN107766373A (zh) | 图片所属类目的确定方法及其系统 | |
CN113946717A (zh) | 一种子图指标特征获得方法、装置、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |