CN108900320A - 一种互联网测试床拓扑结构大比例规模缩减方法及装置 - Google Patents

一种互联网测试床拓扑结构大比例规模缩减方法及装置 Download PDF

Info

Publication number
CN108900320A
CN108900320A CN201810565884.3A CN201810565884A CN108900320A CN 108900320 A CN108900320 A CN 108900320A CN 201810565884 A CN201810565884 A CN 201810565884A CN 108900320 A CN108900320 A CN 108900320A
Authority
CN
China
Prior art keywords
node
bipartite graph
collection
sampling
degree
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.)
Granted
Application number
CN201810565884.3A
Other languages
English (en)
Other versions
CN108900320B (zh
Inventor
焦波
石建迈
张文生
邢立宁
戎海武
何敏藩
于辉
王向东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Foshan University
Original Assignee
Foshan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Foshan University filed Critical Foshan University
Priority to CN201810565884.3A priority Critical patent/CN108900320B/zh
Publication of CN108900320A publication Critical patent/CN108900320A/zh
Priority to PCT/CN2018/124568 priority patent/WO2019233089A1/zh
Application granted granted Critical
Publication of CN108900320B publication Critical patent/CN108900320B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/12—Discovery or management of network topologies
    • 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/14—Network analysis or design
    • H04L41/145—Network analysis or design involving simulating, designing, planning or modelling of a network

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

本发明公开了一种互联网测试床拓扑结构大比例规模缩减方法及装置，通过拓扑二分图结构的分解、采样与合并技术，在保持重要拓扑属性稳定性的条件下，实现自治系统级互联网测试床96％以上的节点删减。本发明可以在保持节点度分布、聚类系数、RichClub连通性和路径长度等重要拓扑属性稳定，即近似保持不变的前提下，大比例地缩减互联网测试床的自治系统级拓扑规模，规模缩减率可以达到96％以上。

Description

一种互联网测试床拓扑结构大比例规模缩减方法及装置
技术领域
本公开涉及数据挖掘应用技术领域，具体涉及一种互联网测试床拓扑结构大比例规模缩减方法及装置。
背景技术
路由协议、资源定位等互联网技术的开发、测试与评估，需要构建相应的仿真测试环境。面对真实互联网网络中每年快速增长的庞大节点规模，相应测试床的构建必须采用规模缩减的方式，以实现构建成本的缩减以及仿真资源利用效率和测试时间运行效率的提高。虽然已有图采样技术已被应用于自治系统级互联网拓扑的规模缩减，但是这些技术的应用范围广泛，没有考虑互联网拓扑独有的结构特征，难以在保持重要拓扑属性稳定性的条件下达到70％以上的规模缩减率(定义为缩减过程删除的节点数占原始拓扑中节点数的比率)。通用性越好的工具，难以适用于专业性较强的领域。因此，本发明关注于自治系统级互联网拓扑独有的二分图结构特征，采用二分图的分解、采样与合并策略，能够在保持重要拓扑属性稳定性的条件下实现真实世界自治系统级互联网拓扑96％以上的规模缩减率。依据压缩规模后的拓扑结构，在小规模互联网测试床上配置拓扑连接关系，可以实现真实大规模拓扑在小规模测试床的等效部署。
发明内容
本公开的目的是针对现有技术的不足，本公开提出一种互联网测试床拓扑结构大比例规模缩减方法及装置，所述方法包括如下步骤：
S1，获取真实世界互联网拓扑结构的探测数据，以简单无向图G＝(V,E)建模自治系统级互联网拓扑结构，其中集合V中节点的物理意义为真实世界互联网的自治系统域，集合E中边的物理意义为不同自治系统域之间的数据通信路径；
S2，采用互联网拓扑结构模型分解互联网拓扑图G＝(V,E)为七个二分图Q1-II、Q2-BI、QI1-RI、Q4-RI1、QI-PI、Q3-PI、Q-P，一个匹配图BI-BI和一个内核图Core；
S3，调用二分图采样方法1对二分图Q1-II、二分图Q2-BI和二分图QI1-RI进行采样，从而抽取部分节点和边构成采样子图；
S4，调用匹配图采样方法对匹配图BI-BI进行采样，抽取部分节点和边构成采样子图；
S5，调用二分图采样方法2对二分图Q4-RI1进行采样，抽取部分节点和边构成采样子图；
S6，调用二分图采样方法3对二分图QI-PI和二分图Q-P进行采样，抽取部分节点和边构成采样子图；
S7，调用二分图采样方法4对二分图Q3-PI进行采样，抽取部分节点和边构成采样子图；
S8，对S3至S7得到的所有采样子图进行合并，得到这些采样子图的合并图；
S9，调用内核图采样方法对内核图Core进行采样，从而抽取部分节点和边构成采样子图；
S10，对S8得到的合并图和S9得到的采样子图进行合并，从而得到真实世界互联网拓扑图G＝(V,E)的采样子图；
S11，依据S10得到的真实世界互联网拓扑图的规模缩减采样子图，配置互联网测试床的拓扑连接关系，从而调度测试床有限的仿真节点资源，生成节点规模大比例缩减的互联网仿真测试环境；
进一步地，对图进行采样，是指从该图中抽取部分节点和边，从而生成由这些被抽取节点和边构成的一个子图。
其中，所述互联网拓扑结构模型分解为八类节点集，其中节点集OI表示系统噪声：
其中，其中V和E分别为互联网拓扑图G＝(V,E)的节点集和边集，v为拓扑图G＝(V,E)中的节点，dv为拓扑图G＝(V,E)中节点v的度；VI,EI：分别为子图GI＝(VI,EI)的节点集和边集；GI＝(VI,EI)为拓扑图G＝(V,E)中，由节点集R(G)＝V/(P∪Q)和两个端点均属于R(G)的所有边构成边集生成的子图；dI(v)为子图GI＝(VI,EI)中节点v的度。
进一步地，互联网拓扑结构模型的八类节点集节点之间的连接关系为：每个P节点仅且必被连接至唯一的Q节点，每个II节点仅且必被连接至至少两个Q节点，每个PI节点仅且必被连接至唯一的QI节点和至少一个Q节点，每个BI节点仅且必被连接至另外一个BI节点和至少一个Q节点，每个RI节点仅且必被连接至至少两个QI节点，此外QI和RI节点可能被连接至Q节点，由连接关系将自治系统级互联网拓扑分解为七个二分图Q1-II,Q-P,Q2-BI,Q3-PI,QI-PI,QI1-RI和Q4-RI1，其中节点集 一个匹配图BI-BI和一个内核图Core。
进一步地，二分图是指图中全部节点可以被划分至两个不同的节点集，且图中每条边的两个端点必定分别属于这两个不同的节点集。二分图Q1-II是指由节点集Q1、节点集II和两个端点分别属于Q1和II的所有边构成边集生成的二分图结构。二分图Q-P,Q2-BI,Q3-PI,QI-PI,QI1-RI和Q4-RI1的物理意义解释方法与二分图Q1-II相同。匹配图BI-BI是由节点集BI和两个端点都属于BI的所有边构成边集生成的图结构。内核图Core是由节点集Q∪QI和连接该节点集内部节点的所有边构成边集生成的图结构。
进一步地，在S2中，互联网拓扑图分解方法包括如下步骤:
S201，读取互联网拓扑图G＝(V,E)的七个节点集P、Q、II、BI、QI、PI和RI；
S202，计算二分图Q1-II：其中且
S203，计算二分图Q2-BI：其中且
S204，计算二分图QI1-RI：其中且
S205，计算二分图Q4-RI1：其中 且
S206，计算二分图QI-PI：GQI-PI＝(QI,PI,EQI-PI)，其中EQI-PI＝{e＝(v,w)|e∈E∧v∈PI∧w∈QI}；
S207，计算二分图Q3-PI：其中且
S208，计算二分图Q-P：GQ-P＝(Q,P,EQ-P)，其中EQ-P＝{e＝(v,w)|e∈E∧v∈P∧w∈Q}；
S209，计算匹配图BI-BI：GBI-BI＝(BI,EBI-BI)，其中EBI-BI＝{e＝(v,w)|e∈E∧v,w∈BI}；
S210，计算内核图Core：GCore＝(VCore,ECore)，其中VCore＝Q∪QI且ECore＝{e＝(v,w)|e∈E∧v,w∈VCore}；
S211，输出七个二分图GQI-PI、GQ-P，一个匹配图GBI-BI和一个内核图GCore。
其中，∈表示“属于”；表示“存在”；∧表示“逻辑与”；∪表示集合的并集；V和E分别表示互联网拓扑图G＝(V,E)的节点集和边集；P、Q、II、BI、QI、PI和RI分别为互联网拓扑结构模型分解得到的七类节点集。
进一步地，在S3中，所述调用二分图采样方法1包括如下步骤：
S301，输入：二分图G＝(V1,V2,E)，其中V1和V2为两个节点集、E为边集，参数T、R，规模缩减率RR；转到步骤S302；
S302，将V2中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集其中这nh个节点中相同度重复数的最大值为T；转到步骤S303；
进一步地，nh的取值由节点集V2中节点的度和参数T唯一确定。
S303，计算并将二分图G分解为和其中且转到步骤S304；
S304，设d1,d2,…,ds为图G中节点集V1的全部节点度，并设dmax为节点集中的最大节点度，计算节点集的度分布{g(d)|d＝1,2,…,dmax}，其中g(d)表示中度为d的节点比率，初始化二分图G′＝(V′1,V′2,E′)←G＝(V1,V2,E)，初始化二分图初始化二分图初始化二分图GLL′中节点集的度分布{g′(d)|d＝1,2,…,dmax}←{g(d)|d＝1,2,…,dmax}，计算需删除边的总数Ed＝||E||·RR，初始化变量t＝0，转到步骤S305；
S305，如果t<Ed，转到步骤S306；否则G′＝(V′1,V′2,E′)，其中且E′＝ELH′∪ELL′，为最终输出结果，转到步骤S311；
S306，如果任意选择集合S中一个节点v，并计算两个可选边集和如果转到步骤S308；如果转到步骤S309；否则且以概率PLL′＝γLL′/(γLH′+γLL′)转到步骤S308，并以概率1-PLL′转到步骤S309，其中，γLL′＝mLL-(||ELL||-||ELL′||)，γLH′＝mLH-(||ELH||-||ELH′||)，mLH＝Ed-mLL；转到步骤S307；
S307，如果S＝Φ，以式定义的离散概率分布pk/Σi＝1,2,…,spi(k＝1,2,…,s)随机采样一个节点度d∈{d1,d2,...,ds}，其中d1>d2>...>ds>0，f(d1),f(d2),...,f(ds)表示节点集V1中度为dk(1≤k≤s)的节点比率；计算两个可选边集 如果转到步骤S308；如果转到步骤S309，否则，以概率PLL′＝γLL′/(γLH′+γLL′)转到步骤S308，并以概率1-PLL′转到步骤S309；
S308，假设其中且计算dh＝argdmaxd∈D{g′(d)-g(d)}，其中D为集合中节点的度组成的集合，并计算的子集从Ea中均等概率分布地采样一条边e，并将该边从图和G′＝(V′1,V′2,E′)中删除，即更新ELL′←ELL′/{e}且E′←E′/{e}，更新GLL′中的节点度分布{g′(d)|d＝1,2,...,dmax}，其中g′(d)表示中度为d的节点比率，转到步骤S310；
S309，从中均等概率分布地采样一条边e，并将该边从二分图和G′＝(V′1,V′2,E′)中删除，即更新ELH′←ELH′/{e}且E′←E′/{e}，转到步骤S310；
S310，更新t←t+1，转到步骤S305；
S311，输出：采样子图G′＝(V′1,V′2,E′)。
其中，||S||表示集合S的势，即集合S包含元素的总数；x←y表示将右边的函数或数值y赋值给左边的变量x；dv表示节点v在图G′＝(V′1,V′2,E′)中的度；A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素；dh＝argdmaxd∈D{g′(d)-g(d)}表示将g′(d)-g(d)取最大值时参数d的取值赋值给dh；∈表示“属于”；表示“不属于”；∧表示“逻辑与”；节点的度是指与该节点相邻边的总数；∪表示集合的并集；输入参数规模缩减率RR被定义为需要从输入的二分图中删除边的总数与输入的二分图包含边数的比率；Φ表示空集，即包含元素个数为零的集合。
进一步地，在S6中，所述调用二分图采样方法3包括如下步骤：
S601，输入：二分图G＝(V1,V2,E)，其中V1和V2为两个节点集、E为边集、节点集V2包含节点的度均为1，节点集参数T，规模缩减率RR；转到步骤S602；
S602，将V1中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集V1 H，其中这nh个节点中相同度重复数的最大值为T；转到步骤S603；
S603，计算V1 L＝V1/V1 H，设d1,d2,…,ds为二分图G中节点集V1的全部节点度，初始化二分图G′＝(V′1,V′2,E′)←G＝(V1,V2,E)，初始化节点集V1 L′←V1 L，初始化节点集V1 H′←V1 H，计算最终需要删除边的总数Ed＝||E||RR，初始化变量t＝0；转到步骤S604；
S604，如果t<Ed，则以概率P＝γL′/(γH′+γL′)转到步骤S605，并以概率1-P转到步骤S608；否则，G′＝(V′1,V′2,E′)，其中V′1＝V1 H′∪V1 L′，为最终输出结果，转到步骤S610；其中，
且mH＝Ed-mL；
S605，如果其中则随机选择集合S中一个节点v，转到步骤S607，否则，转到步骤S606；
S606，以式定义的离散概率分布pk/Σi＝1,2,…,spi(k＝1,2,…,s)，随机采样一个节点度d∈{d1,d2,…,ds}，其中，d1>d2>…>ds>0，f(d1),f(d2),…,f(ds)表示节点集V1在二分图G＝(V1,V2,E)中度为dk(1≤k≤s)的节点比率，计算两个节点集和如果||V2||>0，则随机选择集合V2中一个节点v，否则随机选择集合V1中一个节点v，转到步骤S607；
S607，如果则随机选择一条与节点v相邻的边e，并更新二分图G′＝(V′1,V′2,E′)，E′←E′/{e}，如果dv＝0，则更新V1 L′←V1 L′/{v}，如果则此步不作任何操作；转到步骤S609；
S608，从集合中随机选择一条边e，并更新二分G′＝(V′1,V′2,E′)：E′←E′/{e}；如果dv＝0，则更新V1 H′←V1 H′/{v}；转到步骤S609；
S609，更新t←t+1，转到步骤S604；
S610，输出：采样子图G′＝(V′1,V′2,E′)。
进一步地，输入的节点集被包含于采样子图G′＝(V′1,V′2,E′)，即图G中节点集V2包含节点的度均为1，图G′中节点集V′2包含节点的度均为1；
其中，||S||表示集合S的势，即集合S包含元素的总数；x←y表示将右边的函数或数值y赋值给左边的变量x；当节点v∈V1或v∈V1 H或v∈V1 L，dv表示节点v在图G＝(V1,V2,E)的度，其他情况或没有特定约束下，dv表示节点v在图G′＝(V′1,V′2,E′)的度；A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素；∈表示“属于”；表示“不属于”；∧表示“逻辑与”；∨表示“逻辑或”；节点的度是指与该节点相邻边的总数；表示集合A被包含于集合B；∪表示集合的并集；∩表示集合的交集；输入参数规模缩减率RR被定义为需要从输入的二分图中删除边的总数与输入的二分图包含边数的比率。
进一步地，在S7中，所述调用二分图采样方法4包括如下步骤：
S701，输入：二分图G＝(V1,V2,E)，其中V1和V2为两个节点集、E为边集，多对一映射函数ψ:V1→W，设ψ(v1)为节点集W中与节点v1∈V1对应的唯一节点，对任意节点w∈W预选设定参数L(w)，预选给定节点度分布{g(d)|d＝1,2,…}，预先设定参数T、R，规模缩减率RR；转到步骤S702；
S702，将V2中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集其中这nh个节点中相同度重复数的最大值为T；转到步骤S703；
S703，计算并将二分图G分解为和其中且转到步骤S704；
S704，设d1,d2,…,ds为二分图G中节点集V1的全部节点度，并设dmax为节点集中的最大节点度，初始化二分图G′＝(V′1,V′2,E′)←G＝(V1,V2,E)，初始化二分图初始化二分图计算二分图GLL′中节点集的度分布{g′(d)|d＝1,2,…,dmax}，其中，g′(d)表示中度为d的节点比率，计算最终需要删除边的总数Ed＝||E||RR，其中，初始化变量t＝0，η＝0；转到步骤S705；
S705，如果t<Ed，转到步骤S706；否则二分图G′＝(V′1,V′2,E′)，其中且E′＝ELH′∪ELL′，为最终输出结果，转到步骤S713；
S706，如果则更新η←η-1，计算两个可选边集和
其中且集合N′(ψ(v1))＝{v|v∈V′1∧ψ(v)＝ψ(v1)}，
并转到步骤S709，否则，转到步骤S707；
S707，如果选择集合S中一个节点v，设该节点的度为d＝dv；否则以式定义的离散概率分布pk/Σi＝1,2,…,spi(k＝1,2,…,s)随机采样一个节点度d∈{d1,d2,…,ds}，其中d1>d2>…>ds>0，f(d1),f(d2),…,f(ds)表示节点集V1中度为dk(1≤k≤s)的节点比率，转到步骤S708；
S708，如果则更新η←η+1，转到步骤S712；如果则计算两个可选边集和转到步骤S709；如果d>1，则计算两个可选边集和转到步骤S709；
S709，如果转到步骤S710；如果转到步骤S711；如果 以概率PLL′＝γLL′/(γLH′+γLL′)转到步骤S710，以概率1-PLL′转到步骤S711，
其中，γLL′＝mLL-(||ELL||-||ELL′||)，γLH′＝mLH-(||ELH||-||ELH′||)，
mLH＝Ed-mLL；
S710，假设其中且计算其中D为集合中节点的度组成的集合；如果||Dh||＝0，则更新Dh←{argdmaxd∈D{g′(d)-g(d)}}；
首先，计算集合的一个子集：
假设其中且计算ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}，其中N′(ψ(v))＝{v1|v1∈V′1∧ψ(v1)＝ψ(v)}
然后，计算集合Ea的一个子集，
最后，从集合Eb中随机选择一条边e，并将该边从二分图和G′＝(V′1,V′2,E′)中删除，即更新ELL′←ELL′/{e}且E′←E′/{e}，更新二分GLL′中的节点度分布{g′(d)|d＝1,2,…,dmax}，其中g′(d)表示中度为d的节点比率，转到步骤S712；
S711，假设其中且计算ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}，其中N′(ψ(v))＝{v1|v1∈V′1∧ψ(v1)＝ψ(v)}，
计算集合的一个子集：
并从集合Eb中随机选择一条边e，并将该边从二分图和G′＝(V′1,V′2,E′)中删除，即更新ELH′←ELH′/{e}且E′←E′/{e}，转到步骤S712；
S712，更新t←t+1，转到步骤S705；
S713，输出：采样子图G′＝(V′1,V′2,E′)。
其中，对集合N′(w)＝{v1|v1∈V′1∧ψ(v1)＝w}，其中W为步骤S701的输入，S701至S713期望对满足||N′(w)||接近于L(w)；
其中，||S||表示集合S的势，即集合S包含元素的总数；x←y表示将右边的函数或数值y赋值给左边的变量x；dv表示节点v在图G′＝(V′1,V′2,E′)中的度；A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素；argdmaxd∈D{g′(d)-g(d)}表示g′(d)-g(d)取最大值时参数d的取值；ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}表示将||N′(ψ(v))||/L(ψ(v))取最大值时参数ψ(v)的取值赋值给ψ′；∈表示“属于”；表示“不属于”；∧表示“逻辑与”；∪表示集合的并集；表示“任意”；节点的度是指与该节点相邻边的总数；输入参数规模缩减率RR被定义为需要从输入的二分图中删除边的总数与输入的二分图包含边数的比率。
本发明还提供了一种互联网测试床拓扑结构大比例规模缩减装置，所述装置包括：
建模单元，用于获取真实世界互联网拓扑结构的探测数据，以简单无向图G＝(V,E)建模自治系统级互联网拓扑结构，其中集合V中节点的物理意义为真实世界互联网的自治系统域，集合E中边的物理意义为不同自治系统域之间的数据通信路径；
分解拓扑单元，用于采用互联网拓扑结构模型分解互联网拓扑图G＝(V,E)为七个二分图Q1-II、Q2-BI、QI1-RI、Q4-RI1、QI-PI、Q3-PI、Q-P，一个匹配图BI-BI和一个内核图Core；
二分图采样1单元，用于调用二分图采样方法1对二分图Q1-II、二分图Q2-BI和二分图QI1-RI进行采样，从而抽取部分节点和边构成采样子图；
匹配图采样单元，用于调用匹配图采样方法对匹配图BI-BI进行采样，抽取部分节点和边构成采样子图；
二分图采样2单元，用于调用二分图采样方法2对二分图Q4-RI1进行采样，抽取部分节点和边构成采样子图；
二分图采样3单元，用于调用二分图采样方法3对二分图QI-PI和二分图Q-P进行采样，抽取部分节点和边构成采样子图；
二分图采样4单元，用于调用二分图采样方法4对二分图Q3-PI进行采样，抽取部分节点和边构成采样子图；
子图合并单元，用于对得到的所有采样子图进行合并，得到互联网外围拓扑所有采样子图的合并图；
内核图采样单元，用于调用内核图采样方法对内核图进行采样，从而抽取部分节点和边构成采样子图；
采样子图获取单元，用于对得到的外围拓扑采样子图的合并图和内核拓扑的采样子图进行合并，从而得到真实世界互联网拓扑图的采样子图；
仿真测试单元，依据得到的真实世界互联网拓扑图的规模缩减采样子图，配置互联网测试床的拓扑连接关系，从而调度测试床有限的仿真节点资源，生成节点规模大比例缩减的互联网仿真测试环境。
本公开的有益效果为：本发明公开了一种互联网测试床拓扑结构大比例规模缩减方法，采用面向真实世界互联网拓扑的结构模型，设计测试床的拓扑结构大比例规模缩减技术，在维持拓扑属性稳定不变的前提条件下大比例降低测试床构建所需的仿真资源节点数量，解决有限仿真资源条件下的测试床构建难点问题，从而达到大幅度降低测试床构建成本和测试任务运行时间复杂性等经济效益，本发明可以在保持节点度分布、聚类系数、RichClub连通性和路径长度等重要拓扑属性稳定，即近似保持不变的前提下，大比例地缩减互联网测试床的自治系统级拓扑规模，规模缩减率可以达到96％以上。
附图说明
通过对结合附图所示出的实施方式进行详细说明，本公开的上述以及其他特征将更加明显，本公开附图中相同的参考标号表示相同或相似的元素，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图，在附图中：
图1所示为自治系统级互联网拓扑独有的二分图结构节点分类及它们之间的关系；
图2所示为七个二分图、一个匹配图和一个内核图；
图3所示为互联网测试床拓扑结构大比例规模缩减技术方案的流程图；
图4所示为拓扑规模缩减测试床的组播路由协议测试需求。
具体实施方式
以下将结合实施例和附图对本公开的构思、具体结构及产生的技术效果进行清楚、完整的描述，以充分地理解本公开的目的、方案和效果。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
如图1所示为自治系统级互联网拓扑独有的二分图结构节点分类及它们之间的关系，图2所示为七个二分图、一个匹配图和一个内核图，图3所示为互联网测试床拓扑结构大比例规模缩减技术方案的流程图，图4所示为拓扑规模缩减测试床的组播路由协议测试需求，下面结合图1、图2、图3与图4来阐述根据本公开的实施方式的一种互联网测试床拓扑结构大比例规模缩减方法。
本公开提出一种互联网测试床拓扑结构大比例规模缩减方法，具体包括以下步骤：
S1.获取真实世界互联网拓扑结构的探测数据，并以简单无向图G＝(V,E)建模自治系统级互联网拓扑结构，其中集合V中节点的物理意义为真实世界互联网的自治系统域，集合E中边的物理意义为不同自治系统域之间的数据通信路径。
S2.采用互联网拓扑结构模型，分解互联网拓扑图G＝(V,E)为七个二分图Q1-II、Q2-BI、QI1-RI、Q4-RI1、QI-PI、Q3-PI、Q-P，一个匹配图BI-BI和一个内核图Core。分解方法的具体步骤详见第1.1节。
S3.调用“二分图采样方法1”(具体步骤详见第1.2节)，对二分图Q1-II、二分图Q2-BI和二分图QI1-RI进行采样，从而抽取部分节点和边构成采样子图，具体步骤分别见第1.3节、第1.4节和第1.6节。
S4.调用“匹配图采样方法”，对匹配图BI-BI进行采样，从而抽取部分节点和边构成采样子图，具体步骤详见第1.5节。
S5.调用“二分图采样方法2”，对二分图Q4-RI1进行采样，从而抽取部分节点和边构成采样子图，具体步骤详见第1.7节。
S6.调用“二分图采样方法3”(具体步骤详见第1.8节)，对二分图QI-PI和二分图Q-P进行采样，从而抽取部分节点和边构成采样子图，具体步骤分别见第1.9节和第1.12节。
S7.调用“二分图采样方法4”(具体步骤详见第1.10节)，对二分图Q3-PI进行采样，从而抽取部分节点和边构成采样子图，具体步骤详见第1.11节。
S8.对S3至S7得到的所有采样子图进行合并，得到互联网外围拓扑所有采样子图的合并图，具体步骤详见第1.13节。
S9.调用“内核图采样方法”，对内核图进行采样，从而抽取部分节点和边构成采样子图，具体步骤详见第1.14节。
S10.对S8得到的外围拓扑采样子图的合并图和S9得到的内核拓扑的采样子图进行合并，从而得到真实世界互联网拓扑图的采样子图，具体步骤详见第1.15节。
S11.依据S10得到的真实世界互联网拓扑图的规模缩减采样子图，配置互联网测试床的拓扑连接关系，从而调度测试床有限的仿真节点资源，生成节点规模大比例缩减的互联网仿真测试环境。
自治系统级互联网拓扑可以采用简单无向图G＝(V,E)进行建模，其中V和E分别为节点集和边集。本发明的目标是要得到一个规模大比例缩减(缩减率大于96％)的简单无向图G′，要求G和G′在以下重要拓扑属性上保持相似；其中自治系统级互联网测试床是指该测试床的每个节点为真实互联网中一个自治系统(Autonomous System)的抽象。自治系统是指有权自主确定在本系统中应采用何种路由协议的小型单位。
节点度分布二维属性定义为其中d为节点度，为度大于d的节点数与拓扑图节点总数的比率，f(k)为拓扑图中度为k节点的比率。二维属性在直角坐标系中，是以d为变量且以为函数值的二维函数曲线。
聚类系数二维属性C(d)vs.d：
定义为C(d)＝2Td/d(d-1)，其中Td＝∑i＝1,2,…,tT(vi)/t，d为节点度，C(d)为节点度d的聚类系数，v1,v2,…,vt为拓扑图中所有度为d的节点，t为拓扑图中度为d的节点总数，T(vi)为节点vi任意两个相邻节点之间边的总数。二维属性C(d)vs.d在直角坐标系中，是以d为变量且以C(d)为函数值的二维函数曲线。
RichClub连通性二维属性ρ(r)vs.r：定义ρ(r)为拓扑图中前r个度最大的节点构成子图包含边的总数占据该子图中最大可能边数(即，r(r-1)/2)的比率；r的取值范围为从1到n的全部自然数，其中n为拓扑图包含的节点总数。二维属性ρ(r)vs.r在直角坐标系中，是以r为变量且以ρ(r)为函数值的二维函数曲线。
路径长度分布二维属性μ(l)vs.l：定义μ(l)为拓扑图中相互间最短路径长度为l的节点对总数占据全部节点对总数的比率。二维属性μ(l)vs.l在直角坐标系中，是以l为变量且以μ(l)为函数值的二维函数曲线。
其中vs.表示“相对于”，即表示vs.前的函数是相对于vs.后的变量的函数。例如，vs.d表示是相对于变量d的函数。
为了大比例缩减自治系统级互联网拓扑的节点规模，需要了解该拓扑独有的结构特征。本发明发现该拓扑模型可以被分解为八类节点集：
模型中参数说明如下：
V,E：分别为拓扑图G＝(V,E)的节点集和边集；
dv：为拓扑图G＝(V,E)中节点v的度；
GI＝(VI,EI)：为拓扑图G＝(V,E)中，由节点集R(G)＝V/(P∪Q)和两个端点均属于R(G)的所有边构成边集，生成的子图；
VI,EI：分别为子图GI＝(VI,EI)的节点集和边集；
dI(v)：为子图GI＝(VI,EI)中节点v的度；
∈、∧：分别表示“属于”、“存在”、“任意”和“逻辑与”；
∪：表示集合的并集；
A/B：表示从集合A中删除集合B的所有元素。
上述拓扑模型中OI表示系统噪声节点，其在真实系统占据的比率通常不超过1.9％。为了简化问题，本发明的互联网拓扑规模缩减方法将忽略OI节点。
进一步地，式(1)中七类节点之间的连接关系。如图1所示，每个P节点仅且必被连接至唯一的Q节点，每个II节点仅且必被连接至至少两个Q节点，每个PI节点仅且必被连接至唯一的QI节点和至少一个Q节点，每个BI节点仅且必被连接至另外一个BI节点和至少一个Q节点，每个RI节点仅且必被连接至至少两个QI节点，此外QI和RI节点可能被连接至Q节点。由上述连接关系可以将自治系统级互联网拓扑分解为七个二分图Q1-II,Q-P,Q2-BI,Q3-PI,QI-PI,QI1-RI和Q4-RI1，其中节点集 一个匹配图BI-BI和一个内核图Core，如图2所示。在图2中的箭头表示节点集之间的包含关系，其箭头起点的集合被包含于箭头终点的集合。其中，P节点是指属于集合P的节点，Q节点是指属于集合Q的节点，II节点是指属于集合II的节点，PI节点是指属于集合PI的节点，QI节点是指属于集合QI的节点，BI节点是指属于集合BI的节点，RI节点是指属于集合RI的节点。
其中，二分图是指图中全部节点可以被划分至两个不同的节点集，且图中每条边的两个端点必定分别属于这两个不同的节点集。二分图Q1-II是指由节点集Q1、节点集II和两个端点分别属于Q1和II的所有边构成边集生成的二分图结构。二分图Q-P,Q2-BI,Q3-PI,QI-PI,QI1-RI和Q4-RI1的物理意义解释方法与二分图Q1-II相同。匹配图BI-BI是由节点集BI和两个端点都属于BI的所有边构成边集生成的图结构。内核图Core是由节点集Q∪QI和连接该节点集内部节点的所有边构成边集生成的图结构。
本发明的规模缩减技术，将依据上述自治系统级互联网拓扑独有的二分图结构展开设计。
二分图定义：如果简单无向图G＝(V,E)的节点集V＝V1∪V2，其中V1∩V2＝Φ且对任意一条边e＝(u,w)∈E的两个端点u和w不可能同时属于V1或同时属于V2，则该图可以表示为G＝(V1,V2,E)，且该图被称为二分图。
其中，∪表示集合的并集；∩表示集合的交集；Φ表示不包含任何元素的空集；V和E分别为图G的节点集和边集。
其中，互联网外围拓扑是由互联网外围节点组成的拓扑，互联网内核拓扑是由互联网内核节点组成的拓扑；其中外围节点是指互联网中只进行网络数据的发送与接收，而不进行网络数据转发的节点，内核节点是指互联网中可进行网络数据转发的节点。互联网外围拓扑由其拓扑模型中七个二分图Q1-II、Q2-BI、QI1-RI、Q4-RI1、QI-PI、Q3-PI、Q-P和一个匹配图BI-BI组成，互联网内核图由其拓扑模型中内核图Core构成。
1.1互联网拓扑图分解
输入：原始互联网拓扑图G＝(V,E)。
输出：七个二分图GQI-PI、GQ-P，一个匹配图GBI-BI和一个内核图GCore。
方法的流程与步骤：
S1.依据式(1)计算得到图G的七个节点集P、Q、II、BI、QI、PI和RI。
S2.计算二分图Q1-II：其中且EQ1-II＝{e＝(v,w)|e∈E∧v∈II∧w∈Q1}。
S3.计算二分图Q2-BI：其中且
S4.计算二分图QI1-RI：其中且
S5.计算二分图Q4-RI1：其中 且
S6.计算二分图QI-PI：GQI-PI＝(QI,PI,EQI-PI)，其中EQI-PI＝{e＝(v,w)|e∈E∧v∈PI∧w∈QI}。
S7.计算二分图Q3-PI：其中且
S8.计算二分图Q-P：GQ-P＝(Q,P,EQ-P)，其中EQ-P＝{e＝(v,w)|e∈E∧v∈P∧w∈Q}。
S9.计算匹配图BI-BI：GBI-BI＝(BI,EBI-BI)，其中EBI-BI＝{e＝(v,w)|e∈E∧v,w∈BI}。
S10.计算内核图Core：GCore＝(VCore,ECore)，其中VCore＝Q∪QI且ECore＝{e＝(v,w)|e∈E∧v,w∈VCore}。
其中，∈表示“属于”；表示“存在”；∧表示“逻辑与”；∪表示集合的并集；V和E分别为互联网拓扑图G的节点集和边集。
1.2二分图采样方法1
输入：二分图G＝(V1,V2,E)，其中V1,V2为节点集、E为边集，参数T、R，规模缩减率RR。其中，规模缩减率RR被定义为需要从输入(待采样)的二分图G＝(V1,V2,E)中删除边的总数与该二分图包含边数的比率。
输出：采样子图G′＝(V′1,V′2,E′)。
符号定义：||S||表示集合S的势，即集合S包含元素的总数；
x←y表示将右边的函数或数值y赋值给左边的变量x；
dv表示节点v在图G′＝(V′1,V′2,E′)中的度；
Φ表示不包含任何元素的空集；
A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素；
dh＝argdmaxd∈D{g′(d)-g(d)}表示将g′(d)-g(d)取最大值时参数d的取值赋值给dh；
∈表示“属于”；表示“不属于”；∧表示“逻辑与”；∪表示集合的并集；
节点的度是指与该节点相邻边的总数；
表示集合的势，即集合S包含元素的总数。
方法的流程与步骤：
S1.将节点集V2中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集其中这nh个节点中相同度重复数的最大值为T；也就是说，nh可以由V2包含节点的度和输入参数T唯一确定。转S2。
S2.计算节点集并将二分图G分解为和其中边集且边集转S3。
S3.设d1,d2,…,ds为二分图G中节点集V1的全部节点度，并设dmax为节点集中的最大节点度。计算节点集的度分布{g(d)|d＝1,2,…,dmax}，其中g(d)表示中度为d的节点比率。初始化二分图G′＝(V′1,V′2,E′)←G＝(V1,V2,E)，初始化二分图初始化二分图初始化二分图GLL′中节点集的度分布{g′(d)|d＝1,2,…,dmax}←{g(d)|d＝1,2,…,dmax}。计算需删除边的总数Ed＝||E||RR。初始化变量t＝0。转S4。
S4.如果t<Ed，转S5；否则G′＝(V′1,V′2,E′)(其中且E′＝ELH′∪ELL′)为最终输出结果。
S5.如果任意选择集合S中一个节点v，并计算两个可选边集和如果转S7；如果转S8；否则(即且)，以概率PLL′＝γLL′/(γLH′+γLL′)转S7，并以概率1-PLL′转S8。
其中，γLL′＝mLL-(||ELL||-||ELL′||)，γLH′＝mLH-(||ELH||-||ELH′||)，(2)
S6.如果S＝Φ，以式(4)定义的离散概率分布pk/Σi＝1,2,…,spi(k＝1,2,…,s)随机采样一个节点度d∈{d1,d2,…,ds}，其中d1>d2>…>ds>0，f(d1),f(d2),…,f(ds)表示节点集V1中度为dk(1≤k≤s)的节点比率。
计算两个可选边集 如果转S7；如果转S8；否则，以式(2)和(3)定义的概率PLL′＝γLL′/(γLH′+γLL′)转S7，并以概率1-PLL′转S8。
S7.假设边集其中且计算dh＝argdmaxd∈D{g′(d)-g(d)}，其中D为集合中节点的度组成的集合，并计算边集的一个子集
从Ea中均等概率地采样一条边e，并将该边从二分图和G′＝(V′1,V′2,E′)中删除，即更新边集ELL′←ELL′/{e}和边集E′←E′/{e}。更新二分图GLL′中节点集的度分布{g′(d)|d＝1,2,…,dmax}，其中g′(d)表示中度为d的节点比率。转S9。
S8.从边集中均等概率分布地采样一条边e，并将该边从二分图和G′＝(V′1,V′2,E′)中删除，即更新边集ELH′←ELH′/{e}和边集E′←E′/{e}。转S9。
S9.更新t←t+1，转S4。
方法的应用效果：
本节的二分图采样方法1可以在保证删除Ed＝||E||·RR条边的前提下，使得原始图G＝(V1,V2,E)和采样子图G′＝(V′1,V′2,E′)趋向于拥有相同的节点度分布。其中节点度分布是指图中全部节点的度的概率分布。
1.3二分图采样
输入：第1.1节输出的二分图预先设定参数T(缺省情况下赋值为1)和R(缺省值为0.096)，规模缩减率RR(缺省值为0.96)。
其中缺省情况下参数的取值被称为缺省值，其是系统参数的默认值，即表示一个参数在被修改前的初始值。
输出：采样二分图
方法的流程与步骤：
S1.建立与第1.2节“二分图采样方法1”输入的映射关系：
将赋值于G＝(V1,V2,E)，其中V1＝II，V2＝Q1，
S2.采用“二分图采样方法1”得到的采样子图
S3.如果采样子图的节点集II′中存在度为1的节点，则将这些度为1的节点全部从采样子图中删除。
S4.为最终输出。
1.4二分图采样
输入：第1.1节输出的二分图预先设定参数T(缺省情况下赋值为1)和R(缺省值为0.026)，规模缩减率RR(缺省值为0.96)。
输出：采样二分图
方法的流程与步骤：
S1.建立与第1.2节“二分图采样方法1”输入的映射关系：
将赋值于G＝(V1,V2,E)，其中V1＝BI，V2＝Q2，
S2.调用“二分图采样方法1”得到的采样子图
S3.如果采样子图的节点集BI′中节点总数为奇数，则在BI′中选择一个度最小的节点，并将该节点从采样子图中删除。
S4.为最终输出。
1.5匹配图采样方法(匹配图GBI-BI＝(BI,EBI-BI)采样)
输入：第1.1节输出的二分图和匹配图GBI-BI＝(BI,EBI-BI)，第1.4节输出的采样二分图
输出：采样匹配图G′BI-BI＝(BI′,E′BI-BI)。
符号定义：表示BI节点v在二分图中的度，BI节点v是指属于集合BI的节点v。
表示BI′节点v在采样图中的度，BI′节点v是指属于集合BI′的节点v。
||S||表示集合S的势，即集合S包含元素的总数。
Φ表示空集；表示元素0到x之间最接近于x的整数。
x←y表示将右边的函数或数值y赋值给左边的变量x。
max(x,y)表示取x和y中最大值；min(x,y)表示取x和y中最小值。
∈表示“属于”；∧表示“逻辑与”；表示“任意”；∪表示集合的并集。
方法的流程与步骤：
S1.对每条边e＝(u,v)∈EBI-BI(其中u,v∈BI)，生成一个节点度对(不失一般性，假设)，并将节点度对映射为计算集合并计算节点度对分布函数ρ(x)＝||{y|y∈X∧y＝x}||/||X||，其表示集合X中出现元素x的概率值。
S2.设采样二分图中BI′节点的最大度为dm。初始化边集E′BI-BI＝Φ。将节点度d1以步长1从1增加至dm-1，并将节点度d2以步长1从d1增加至dm：对每个节点度对(d1,d2)，首先计算其期望数量和最大可能数量其中
然后，标准分布(概率均等)地随机生成条边ei＝(ui,vi)(i＝1,2,…,s)，要求对满足最后，更新边集E′BI-BI←E′BI-BI∪{ei}i＝1,2,…,s。
S3.对于节点集标准分布地随机生成t＝||BI′||/2-||E′BI-BI||条边ei＝(ui,vi)(i＝1,2,…,t)，要求对满足ui∈Y∧vi∈Y∧ui≠uj∧ui≠vj∧vi≠vj，并更新边集E′BI-BI←E′BI-BI∪{ei}i＝1,2,…,t。此时，生成的采样匹配图G′BI-BI＝(BI′,E′BI-BI)为本节的输出结果。
1.6二分图采样
输入：第1.1节输出的二分图预先设定参数T(缺省情况下赋值为1)和R(缺省值为0)，规模缩减率RR(缺省值为0.96)。
输出：采样二分图
方法的流程与步骤：
S1.建立与第1.2节“二分图采样方法1”输入的映射关系：
将赋值于G＝(V1,V2,E)，其中V1＝RI，V2＝QI1，
S2.调用“二分图采样方法1”得到的采样子图
S3.如果采样子图的节点集RI′中存在度为1的节点，则将这些度为1的节点全部从采样子图中删除。
S4.为最终输出。
1.7二分图采样方法2(二分图采样)
输入：第1.1节输出的二分图第1.6节输出的采样二分图
输出：采样二分图
符号定义：∈表示“属于”；∧表示“逻辑与”；表示“存在”。
方法的流程与步骤：
S1.计算节点集节点集和边集
S2.则，由节点集RI′1、节点集Q′4和边集生成的二分图即为最终的输出结果。
1.8二分图采样方法3
输入：二分图G＝(V1,V2,E)，其中V1,V2为节点集、E为边集，节点集参数T，规模缩减率RR。其中，规模缩减率RR被定义为需要从输入(待采样)的二分图G＝(V1,V2,E)中删除边的总数与该二分图包含边数的比率。
输出：采样子图G′＝(V′1,V′2,E′)。
约束条件：节点集要求被包含于采样子图，即
二分图G中节点集V2包含节点的度均为1。
采样子图G′中节点集V′2包含节点的度均为1。
符号定义：||S||表示集合S的势，即集合S包含元素的总数。
x←y表示将右边的函数或数值y赋值给左边的变量x。
当节点v∈V1或v∈V1 H或v∈V1 L(其中V1 H和V1 L将在S1、S2中定义)，dv表示v在图G＝(V1,V2,E)的度。其他情况或没有特定约束下，dv表示v在图G′＝(V′1,V′2,E′)的度。∈表示“属于”；表示“不属于”；∧表示“逻辑与”；∨表示“逻辑或”。
A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素。
表示集合的势，即集合S包含元素的总数。
方法的流程与步骤：
S1.将节点集V1中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集V1 H，其中这nh个节点中相同度重复数的最大值为T；也就是说，参数nh可以由节点集V1包含全部节点的度和输入参数T唯一确定。转S2。
S2.计算节点集V1 L＝V1/V1 H，设d1,d2,…,ds为二分图G中节点集V1的全部节点度。初始化二分图G′＝(V′1,V′2,E′)←G＝(V1,V2,E)，初始化节点集V1 L′←V1 L，初始化节点集V1 H′←V1 H。计算最终需要删除边的总数Ed＝||E||·RR。初始化变量t＝0。转S3。
S3.如果t<Ed，则以概率P＝γL′/(γH′+γL′)转S4，并以概率1-P转S7；否则，G′＝(V′1,V′2,E′)(其中V′1＝V1 H′∪V1 L′)为最终输出结果。其中，
且mH＝Ed-mL。(9)
S4.如果其中节点集则随机选择集合S中一个节点v，转S6；否则，转S5。
S5.以式(10)定义的离散概率分布pk/Σi＝1,2,…,spi(k＝1,2,…,s)随机采样一个节点度d∈{d1,d2,…,ds}，其中d1>d2>…>ds>0，f(d1),f(d2),…,f(ds)表示节点集V1在二分图G＝(V1,V2,E)中度为dk(1≤k≤s)的节点比率。
计算两个节点集和如果||V2||>0，则随机选择集合V2中一个节点v，否则随机选择集合V1中一个节点v；转S6.
S6.如果则随机选择一条与节点v相邻的边e，并更新二分图G′＝(V′1,V′2,E′)：E′←E′/{e}；如果dv＝0，则更新节点集V1 L′←V1 L′/{v}。如果则此步不作任何操作。转S8。
S7.从集合中随机选择一条边e，并更新二分图G′＝(V′1,V′2,E′)：E′←E′/{e}；如果dv＝0，则更新节点集V1 H′←V1 H′/{v}。转S8。
S8.更新t←t+1，转S3。
方法的应用效果：
本节的二分图采样方法3可以在保证删除Ed＝||E||·RR条边的前提下，使得原始图G＝(V1,V2,E)和采样子图G′＝(V′1,V′2,E′)趋向于拥有相同的节点度分布，并要求采样子图G′＝(V′1,V′2,E′)保证满足其中节点度分布是指图中全部节点的度的概率分布。
1.9二分图GQI-PI＝(QI,PI,EQI-PI)采样
输入：第1.1节输出的二分图GQI-PI＝(QI,PI,EQI-PI)，第1.6节输出的采样二分图预先设定参数T(缺省情况下赋值为0)，规模缩减率RR(缺省情况下赋值为0.96)。
输出：采样二分图G′QI-PI＝(QI′,PI″,E′QI-PI)。
方法的流程与步骤：
S1.建立与第1.8节“二分图采样方法3”输入G＝(V1,V2,E)和的映射关系：
将GQI-PI＝(QI,PI,EQI-PI)赋值于G＝(V1,V2,E)，其中V1＝QI，V2＝PI，E＝EQI-PI；并设节点集
S2.调用“二分图采样方法3”得到GQI-PI的采样子图G′QI-PI＝(QI′,PI″,E′QI-PI)。
S3.则，G′QI-PI为最终输出。
1.10二分图采样方法4
输入：二分图G＝(V1,V2,E)，其中V1,V2为节点集、E为边集，多对一映射函数ψ:V1→W(设ψ(v1)为节点集W中与节点v1∈V1对应的唯一节点)，对任意节点w∈W预选设定参数L(w)，预选给定节点度分布{g(d)|d＝1,2,…}，预先设定参数T、R，规模缩减率RR。其中，规模缩减率RR被定义为需要从输入(待采样)的二分图G＝(V1,V2,E)中删除边的总数与该二分图包含边数的比率。
输出：采样子图G′＝(V′1,V′2,E′)。
对定义集合N′(w)＝{v1|v1∈V′1∧ψ(v1)＝w}，其中
本节的二分图采样方法4期望对满足||N′(w)||接近于L(w)。
符号定义：||S||表示集合S的势，即集合S包含元素的总数。
x←y表示将右边的函数或数值y赋值给左边的变量x。
dv表示节点v在图G′＝(V′1,V′2,E′)中的度。
∈表示“属于”；表示“不属于”；∧表示“逻辑与”；表示“任意”。
∪表示集合的并集。
A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素。
argdmaxd∈D{g′(d)-g(d)}表示g′(d)-g(d)取最大值时参数d的取值。
ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}表示将||N′(ψ(v))||/L(ψ(v))取最大值时参数ψ(v)的取值赋值给ψ′。
节点的度是指与该节点相邻边的总数。
表示集合的势，即集合S包含元素的总数。
方法的流程与步骤：
S1.将节点集V2中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集其中这nh个节点中相同度重复数的最大值为T；也就是说，参数nh可以由节点集V2包含全部节点的度和输入参数T唯一确定。转S2。
S2.计算节点集并将二分图G分解为和其中边集且边集
S3.设d1,d2,...,ds为二分图G中节点集V1的全部节点度，并设dmax为节点集中的最大节点度。初始化二分图G′＝(V′1,V′2,E′)←G＝(V1,V2,E)，初始化二分图初始化二分图计算节点集的度分布{g′(d)|d＝1,2,...,dmax}，其中g′(d)表示节点集中度为d的节点比率。计算最终需要删除边的总数Ed＝||E||RR。初始化变量t＝0，η＝0。转S4。
S4.如果t<Ed，转S5；否则G′＝(V′1,V′2,E′)(其中且E′＝ELH′∪ELL′)为最终输出结果。
S5.如果则更新η←η-1，计算两个可选边集和其中
且N′(ψ(v1))＝{v|v∈V′1∧ψ(v)＝ψ(v1)}，并转S8；否则，转S6。
S6.如果选择集合S中一个节点v，设该节点的度为d＝dv；否则以式(12)定义的离散概率分布pk/Σi＝1,2,...,spi(k＝1,2,...,s)随机采样一个节点度d∈{d1,d2,…,ds}，其中d1>d2>…>ds>0，f(d1),f(d2),…,f(ds)表示节点集V1中度为dk(1≤k≤s)的节点比率。转S7。
S7.如果则更新η←η+1，转S11；如果则计算两个可选边集和转S8；如果d>1，则计算两个可选边集和转S8。
S8.如果转S9；如果转S10；如果以概率PLL′＝γLL′/(γLH′+γLL′)转S9，以概率1-PLL′转S10。
其中，γLL′＝mLL-(||ELL||-||ELL′||)，γLH′＝mLH-(||ELH||-||ELH′||)，(13)
S9.假设其中且计算其中D为集合中节点的度组成的集合；如果||Dh||＝0，则更新Dh←{argdmaxd∈D{g′(d)-g(d)}}。
首先，计算集合的一个子集：
假设其中且计算ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}，其中N′(ψ(v))＝{v1|v1∈V′1∧ψ(v1)＝ψ(v)}。
然后，计算集合Ea的一个子集：
最后，从集合Eb中随机选择一条边e，并将该边从二分图和G′＝(V′1,V′2,E′)中删除，即更新边集ELL′←ELL′/{e}和边集E′←E′/{e}。更新二分图GLL′中节点集的节点度分布{g′(d)|d＝1,2,…,dmax}，其中g′(d)表示中度为d的节点比率。转S11。
S10.假设其中且计算ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}，其中N′(ψ(v))＝{v1|v1∈V′1∧ψ(v1)＝ψ(v)}。
计算集合的一个子集：
并从集合Eb中随机选择一条边e，并将该边从二分图和G′＝(V′1,V′2,E′)中删除，即更新边集ELH′←ELH′/{e}和边集E′←E′/{e}。转S11。
S11.更新t←t+1，转S4。
方法的应用效果：
本节的二分图采样方法4可以在保证删除Ed＝||E||·RR条边的前提下，使得原始图G＝(V1,V2,E)和采样子图G′＝(V′1,V′2,E′)趋向于拥有相同的节点度分布，并要求采样子图G′＝(V′1,V′2,E′)保证对满足||N′(w)||接近于L(w)。其中节点度分布是指图中全部节点的度的概率分布。
在第1.2节二分图采样方法1基础上设计本节二分图采样方法4的目标，是为了保证图2中二分图QI-PI和二分图Q3-PI拥有相同的PI节点集。
1.11二分图采样
输入：第1.1节输出的二分图GQI-PI＝(QI,PI,EQI-PI)和第1.9节输出的采样二分图G″QI-PI＝(QI′,PI″,E″QI-PI)，预先设定参数T(缺省情况下赋值为1)、R(缺省情况下赋值为1)，规模缩减率RR(缺省情况下赋值为0.96)。
输出：采样二分图和G′QI-PI＝(QI′,PI′,E′QI-PI)。
符号定义：Nv(G)表示节点v在图G的相邻节点集。
dv(G)表示节点v在图G的度。
||S||表示集合S的势，即集合S包含元素的总数。
∈表示“属于”；表示“存在”；∧表示“逻辑与”；∪表示集合的并集。
方法的流程与步骤：
S1.计算节点集节点集和边集定义由这两个节点集和一个边集生成的二分图为
S2.建立与第1.10节“二分图采样方法4”输入的映射关系：
将赋值于G＝(V1,V2,E)，其中V2＝Q″3，设定多对一映射函数为：对ψ(v)∈QI′∧(v,ψ(v))∈EQI-PI；对设定预先参数L(w)＝dw(G″QI-PI)；设定{g(d)|d＝1,2,…}为二分图中节点集Q3包含节点的度分布，其中g(d)表示度为d的节点比率。
S3.调用二分图采样方法4对删除mo·RR-(mo-m″o)条边，其中且并生成采样结果图
S4.为了保证采样结果和二分图G″QI-PI＝(QI′,PI″,E″QI-PI)拥有相同的PI节点集，将二分图G″QI-PI＝(QI′,PI″,E″QI-PI)更新为二分图G′QI-PI＝(QI′,PI′,E′QI-PI)，其中边集E′QI-PI＝{(w,v)∈EQI-PI|w∈QI′∧v∈PI′}。
S5.则，和G′QI-PI为最终的输出结果。
1.12二分图GQ-P＝(Q,P,EQ-P)采样
输入：第1.1节输出的二分图GQ-P＝(Q,P,EQ-P)，第1.3节输出的采样二分图第1.4节输出的采样二分图第1.7节输出的采样二分图第1.11节输出的采样二分图参数T(缺省值为3)，规模缩减率RR(缺省值为0.96)。
输出：采样二分图G′Q-P＝(Q′,P′,E′Q-P)和Q′1、Q′2、Q′3、Q′4的分解节点集Q′11、Q′21、Q′31、Q′41、Q′12、Q′22、Q′32、Q′42。
符号定义：dv(G)表示节点v在图G的度。
∪表示集合的并集。
A/B：如果A和B为集合，则表示从集合A中删除集合B的所有元素。
方法的流程与步骤：
S1.分解节点集Q′1、Q′2、Q′3和Q′4如下：
S2.建立与第1.8节“二分图采样方法3”输入G＝(V1,V2,E)和的映射关系：
将GQ-P＝(Q,P,EQ-P)赋值于G＝(V1,V2,E)，其中V1＝Q，V2＝P，E＝EQ-P。并设节点集
S3.调用二分图采样方法3得到GQ-P的采样子图G′Q-P＝(Q′,P′,E′Q-P)。
S4.则，G′Q-P和式(18)的八个分解节点集Q′11、Q′21、Q′31、Q′41、Q′12、Q′22、Q′32、Q′42为最终输出。
1.13互联网外围拓扑所有采样子图的合并
输入：第1.3节输出的采样二分图第1.4节输出的采样二分图第1.5节输出的采样匹配图G′BI-BI＝(BI′,E′BI-BI)，第1.6节输出的采样二分图第1.7节输出的采样二分图第1.11节输出的采样二分图和G′QI-PI＝(QI′,PI′,E′QI-PI)，第1.12节输出的采样二分图G′Q-P＝(Q′,P′,E′Q-P)和节点集Q′11、Q′21、Q′31、Q′41、Q′12、Q′22、Q′32、Q′42，原始互联网拓扑图G＝(V,E)。
输出：合并子图G′m，内核子图G′Core。
符号定义：||S||表示集合S的势，即集合S包含元素的总数。
dv(G)表示节点v在图G的度。
Nv(G)表示节点v在图G的相邻节点集。
表示元素0到x之间最接近于x的整数。
x←y表示将右边的函数或数值y赋值给左边的变量x。
∈表示“属于”；表示“存在”；∪表示集合的并集；∩表示集合的交集。
A/B：如果A和B为集合，则表示从集合A中删除集合B的所有元素。
表示集合A中的所有元素都被包含于集合B。
Φ表示不包含任何元素的空集。
方法的流程与步骤：
S1.从节点集Q′/Q′11中均等概率地抽取||Q′12||个节点Qs＝{q1,q2,…}，并将采样二分图节点集Q′1包含的节点子集Q′12替换为Qs。类似地，采用从Q′/Q′21、Q′/Q′31和Q′/Q′41抽取的节点集，分别替换的节点子集Q′22、的节点子集Q′32和的节点子集Q′42。通过替换之后，可以满足
S2.计算节点集V″m＝II′∪BI′∪PI′∪RI′∪QI′∪Q′∪P′和边集 则，合并图G″m＝(V″m,E″m)为所有采样二分图和采样匹配图的合并结果。
S3.设V′Q,QI为合并图G″m中所有Q和QI节点组成的节点集(注：第1.1节“互联网拓扑图分解”已经将拓扑图中全部节点标记为P、Q、II、BI、QI、PI和RI七类不同的节点，系统噪声OI节点已被从拓扑图中删除)，并设V′QI为合并图G″m中所有QI节点组成的节点集。定义G′Core＝(V′Q,QI,E′Q,QI)为原始拓扑G中由节点集V′Q,QI及这些节点之间的边生成的子图，并定义G′inner＝(V′QI,E′QI)为原始拓扑G中由节点集V′QI及这些节点之间的边生成的子图。
S4.计算节点集和其中G″inner为图G″m中由所有PI、RI、QI节点及这些节点之间的边生成的子图。将节点集分解为和并将节点集分解为和
S5.初始化边集和对每个节点均等概率地选择一个节点w∈Nv(G′Core)，并更新其中(w,v)∈E′Q,QI。此外，对每个节点均等概率地选择一个节点w∈Nv(G″inner)，并更新其中
S6.均等概率地生成x条连接节点集中2x个节点(定义Vx为这2x个节点的集合)的边其中并更新如果则均等概率地选择一个节点w∈V′Q,QI/{v}，并更新
S7.均等概率地生成y条连接节点集中2y个节点(定义Vy为这2y个节点的集合)的边其中并更新如果则均等概率地选择一个节点w∈V′QI/{v}，并更新
S8.更新G′Core＝(V′Q,QI,E′Q,QI)：并计算则，G′Core和G′m为最终的输出结果。
1.14内核图采样方法
输入：原始互联网拓扑G＝(V,E)，第1.13节输出的合并子图G′m和内核子图G′Core＝(V′Core,E′Core)。
输出：内核边集E″Core。
符号定义：||S||表示集合S的势，即集合S包含元素的总数。
dv(G)表示节点v在图G的度。
∈表示“属于”；max(x,y)表示取x和y中的最大值。
平均节点度是指图中全部节点的度之和除以全部节点的总数。
方法的流程与步骤：
S1.计算原始拓扑G的平均节点度davg，计算合并子图G′m的边数m(G′m)和节点数n(G′m)。将所有边(w,v)∈E′Core按照max(dw(G′Core),dv(G′Core))不降顺序排列。
S2.从边集E′Core中删除排在前面的||E′Core||-davg·n(G′m)/2+m(G′m)条边，并将剩余边组成的集合定义为E″Core(其包含davg·n(G′m)/2-m(G′m)条边)。
S3.E″Core为最终输出的内核边集。
1.15互联网拓扑最终采样结果
输入：第1.13节输出的合并子图G′m，第1.14节输出的内核边集E″Core。
输出：最终采样结果G′。
符号定义：∪表示集合的并集。
方法的流程与步骤：
S1.真实世界互联网拓扑的最终采样结果为G′＝G′m∪E″Core。
实施例以互联网测试床的组播路由协议测试任务为案例。选择该任务需求的原因在于组播路由协议测试指标对测试床拓扑结构的变化比较敏感。实施例以Stanford网络数据池(http://snap.stanford.edu/data/)和ITDK工程(http://www.caida.org/)提供的三个真实世界自治系统级互联网拓扑的探测数据为原始拓扑。这三个原始拓扑分别来源于AS-733(2000年01月探测)、AS-Caida(2007年11月探测)和ITDK(2017年02月探测)。
测试任务首先采用本发明的拓扑采样技术，对真实世界探测得到的原始拓扑进行规模缩减(要求在保持节点度分布、聚类系数、RichClub连通性和路径长度等重要拓扑属性近似稳定不变前提下，实现大于96％的规模缩减率)，得到小规模拓扑结构的互联网测试床；然后在小规模测试床上计算组播路由协议测试指标(包含延时率和费用率)的测量数值；最后通过与真实世界拓扑结构上测试结果的对比，分析本发明(拓扑规模缩减技术)的实用性，即是否能够由小规模测试床上的测试结论等效推演得到组播路由协议在大规模真实网络的运行效果。
1、组播路由协议测试
互联网组播架构可以划分为组播路由协议、组播传输协议和群组管理协议三类。其中组播路由协议构建组播拓扑(即，群组内源节点到多个接收节点的组播数据路径分布)，组播传输协议管控组播拓扑上数据流的传输机制，组播管理协议负责群组成员的管理。因此，组播路由协议紧密依赖于互联网测试床的拓扑结构，组播拓扑通常是测试床拓扑的一个子集。目前，组播路由协议可以划分为基于源的树SBT(Source Based Tree)协议和基于中心的树CBT(Center Based Tree)协议两类。SBT和CBT采用不同的机制生成包含群组全部成员的组播拓扑，其中SBT由许多最短路径树的并集组成(每个树以群组内的一个数据发送源节点为根，形成了从该源节点发送数据的组播路径)；反之，CBT由唯一的最短路径树组成，该树被群组内所有的成员节点共享。
以某个源节点为根的最短路径树，由该源节点至群组中所有接收节点的最短路径组成。一般地，从一个源节点到一个接收节点的延时，采用它们之间的最短路径长度来测量。因此，SBT可以得到从每个源节点至所有接收节点的最小平均延时。然而，费用是评估该协议的另一个重要指标，其被定义为组播分布树(即，组播拓扑)覆盖互联网测试床拓扑的链路(边)总数。易知，CBT可以得到较优的费用指标。虽然如何使平均延时和费用之间达到平衡仍是组播路由领域的一个开放问题，但是CBT选择网络拓扑的“中心”作为其唯一最短路径树的根节点，已成为当前应用较为广泛的组播路由协议。以什么标准定义“网络中心”，以及“网络中心”如何选择，仍是CBT协议簇的一个重要研究方向。
目前，CBT协议的评价通常以SBT协议为基准，因为SBT可以达到最优(小)的平均延时指标。常用的CBT协议评价指标包括：
延时率：被定义为SBT的平均延时与CBT的平均延时的比率。
费用率：被定义为SBT的费用与CBT的费用的比率。
选择“网络中心”(即，唯一最短路径树的根节点)对于CBT协议至关重要。然而，本实施例并不关心“网络中心”的选择方法。精确地计算“网络中心”已被证明是NP完全问题；也就是说，这种精确计算难以被用于大规模网络的协议设计。因此，CBT协议的网络中心通常采用启发式的方法计算。
本实施例重点关心自治系统级互联网测试床拓扑的规模缩减，不关心被试CBT技术的优化。因此，仅考虑一个具体的CBT协议，其“网络中心”的选择策略为：标准分布地随机选择一个源节点作为“网络中心”。标准分布是指分布概率值均等的概率分布。
2、测试任务验证结果
拓扑属性验证结果：
本发明的技术针对AS-733(2000年01月探测，包含6,474节点)、AS-Caida(2007年11月探测，包含26,389节点)和ITDK(2017年02月探测，包含41,626节点)三个真实世界的自治系统级互联网拓扑，都可以在保证节点度分布、聚类系数、RichClub连通性和路径长度分布等重要拓扑属性近似稳定不变的前提下，实现96％以上的节点规模缩减。
组播路由协议测试验证结果：
针对上述AS-733、AS-Caida和ITDK三个真实世界原始拓扑，依据本发明技术采样得到的小规模拓扑结构(节点规模缩减率大于96％)，在互联网测试床上配置拓扑连接关系。通过实验对比验证得到：在小规模测试床上组播路由协议的测试结论与在大规模真实网络上组播路由协议的运行效果基本保持一致。
本公开的实施例提供的一种互联网测试床拓扑结构大比例规模缩减装置，该实施例的一种互联网测试床拓扑结构大比例规模缩减装置包括：处理器、存储器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，例如拓扑程序。所述处理器执行所述计算机程序时实现上述各个方法实施例中的步骤。
所述装置包括：存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序的装置包括：
建模单元，用于获取真实世界互联网拓扑结构的探测数据，以简单无向图G＝(V,E)建模自治系统级互联网拓扑结构，其中集合V中节点的物理意义为真实世界互联网的自治系统域，集合E中边的物理意义为不同自治系统域之间的数据通信路径；
分解拓扑单元，用于采用互联网拓扑结构模型分解互联网拓扑图G＝(V,E)为七个二分图Q1-II、Q2-BI、QI1-RI、Q4-RI1、QI-PI、Q3-PI、Q-P，一个匹配图BI-BI和一个内核图Core；
二分图采样1单元，用于调用二分图采样方法1对二分图Q1-II、二分图Q2-BI和二分图QI1-RI进行采样，从而抽取部分节点和边构成采样子图；
匹配图采样单元，用于调用匹配图采样方法对匹配图BI-BI进行采样，抽取部分节点和边构成采样子图；
二分图采样2单元，用于调用二分图采样方法2对二分图Q4-RI1进行采样，抽取部分节点和边构成采样子图；
二分图采样3单元，用于调用二分图采样方法3对二分图QI-PI和二分图Q-P进行采样，抽取部分节点和边构成采样子图；
二分图采样4单元，用于调用二分图采样方法4对二分图Q3-PI进行采样，抽取部分节点和边构成采样子图；
子图合并单元，用于对得到的所有采样子图进行合并，得到互联网外围拓扑所有采样子图的合并图；
内核图采样单元，用于调用内核图采样方法对内核图进行采样，从而抽取部分节点和边构成采样子图；
采样子图获取单元，用于对得到的外围拓扑采样子图的合并图和内核拓扑的采样子图进行合并，从而得到真实世界互联网拓扑图的采样子图；
仿真测试单元，依据得到的真实世界互联网拓扑图的规模缩减采样子图，配置互联网测试床的拓扑连接关系，从而调度测试床有限的仿真节点资源，生成节点规模大比例缩减的互联网仿真测试环境。
所述一种互联网测试床拓扑结构大比例规模缩减装置可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述一种互联网测试床拓扑接结构大比例规模缩减装置可包括，但不仅限于，处理器、存储器。本领域技术人员可以理解，所述例子仅仅是一种互联网测试床拓扑结构大比例规模缩减装置的示例，并不构成对一种互联网测试床拓扑结构大比例规模缩减装置的限定，可以包括比例子更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述一种互联网测试床拓扑结构大比例规模缩减装置还可以包括输入输出设备、网络接入设备、总线等。
所称处理器可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述处理器是所述一种互联网测试床大比例规模缩减装置的控制中心，利用各种接口和线路连接整个一种互联网测试床拓扑结构大比例规模缩减装置的各个部分。
所述存储器可用于存储所述计算机程序和/或模块，所述处理器通过运行或执行存储在所述存储器内的计算机程序和/或模块，以及调用存储在存储器内的数据，实现所述一种互联网测试床拓扑结构大比例规模缩减装置的各种功能。所述存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等；存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
尽管本公开的描述已经相当详尽且特别对几个所述实施例进行了描述，但其并非旨在局限于任何这些细节或实施例或任何特殊实施例，而是应当将其视作是通过参考所附权利要求考虑到现有技术为这些权利要求提供广义的可能性解释，从而有效地涵盖本公开的预定范围。此外，上文以发明人可预见的实施例对本公开进行描述，其目的是为了提供有用的描述，而那些目前尚未预见的对本公开的非实质性改动仍可代表本公开的等效改动。

Claims (7)

1.一种互联网测试床拓扑结构大比例规模缩减方法，其特征在于，所述方法包括如下步骤：
S1，获取真实世界互联网拓扑结构的探测数据，以简单无向图G＝(V,E)建模自治系统级互联网拓扑结构，其中集合V中节点的物理意义为真实世界互联网的自治系统域，集合E中边的物理意义为不同自治系统域之间的数据通信路径；
S2，采用互联网拓扑结构模型分解互联网拓扑图G＝(V,E)为七个二分图Q1-II、Q2-BI、QI1-RI、Q4-RI1、QI-PI、Q3-PI、Q-P，一个匹配图BI-BI和一个内核图Core；
S3，调用二分图采样方法1对二分图Q1-II、二分图Q2-BI和二分图QI1-RI进行采样，从而抽取部分节点和边构成采样子图；
S4，调用匹配图采样方法对匹配图BI-BI进行采样，抽取部分节点和边构成采样子图；
S5，调用二分图采样方法2对二分图Q4-RI1进行采样，抽取部分节点和边构成采样子图；
S6，调用二分图采样方法3对二分图QI-PI和二分图Q-P进行采样，抽取部分节点和边构成采样子图；
S7，调用二分图采样方法4对二分图Q3-PI进行采样，抽取部分节点和边构成采样子图；
S8，对S3至S7得到的所有采样子图进行合并，得到这些采样子图的合并图；
S9，调用内核图采样方法对内核图Core进行采样，从而抽取部分节点和边构成采样子图；
S10，对S8得到的合并图和S9得到的采样子图进行合并，从而得到真实世界互联网拓扑图G＝(V,E)的采样子图；
S11，依据S10得到的真实世界互联网拓扑图的规模缩减采样子图，配置互联网测试床的拓扑连接关系，从而调度测试床有限的仿真节点资源，生成节点规模大比例缩减的互联网仿真测试环境；
其中，所述互联网拓扑结构模型分解为八类节点集，其中节点集OI表示系统噪声：
其中，V和E分别为互联网拓扑图G＝(V,E)的节点集和边集，v为拓扑图G＝(V,E)中的节点，dv为拓扑图G＝(V,E)中节点v的度；GI＝(VI,EI)为拓扑图G＝(V,E)中由节点集R(G)＝V/(P∪Q)和两个端点均属于R(G)的所有边构成边集生成的子图；VI,EI：分别为子图GI＝(VI,EI)的节点集和边集；dI(v)为子图GI＝(VI,EI)中节点v的度。
2.根据权利要求1所述的一种互联网测试床拓扑结构大比例规模缩减方法，其特征在于，互联网拓扑结构模型的八类节点集节点之间的连接关系为：每个P节点仅且必被连接至唯一的Q节点，每个II节点仅且必被连接至至少两个Q节点，每个PI节点仅且必被连接至唯一的QI节点和至少一个Q节点，每个BI节点仅且必被连接至另外一个BI节点和至少一个Q节点，每个RI节点仅且必被连接至至少两个QI节点，此外QI和RI节点可能被连接至Q节点，由连接关系将自治系统级互联网拓扑分解为七个二分图Q1-II,Q-P,Q2-BI,Q3-PI,QI-PI,QI1-RI和Q4-RI1，其中节点集 一个匹配图BI-BI和一个内核图Core，
其中，二分图是指图中全部节点可以被划分至两个不同的节点集，且图中每条边的两个端点必定分别属于这两个不同的节点集，二分图Q1-II是指由节点集Q1、节点集II和两个端点分别属于Q1和II的所有边构成边集生成的二分图结构，二分图Q-P,Q2-BI,Q3-PI,QI-PI,QI1-RI和Q4-RI1的物理意义解释方法与二分图Q1-II相同；匹配图BI-BI是由节点集BI和两个端点都属于BI的所有边构成边集生成的图结构；内核图Core是由节点集Q∪QI和连接该节点集内部节点的所有边构成边集生成的图结构。
3.根据权利要求1所述的一种互联网测试床拓扑结构大比例规模缩减方法，其特征在于，在S2中，互联网拓扑图分解方法包括如下步骤:
S201，读取互联网拓扑图G＝(V,E)的七个节点集P、Q、II、BI、QI、PI和RI；
S202，计算二分图Q1-II：其中且
S203，计算二分图Q2-BI：其中且
S204，计算二分图QI1-RI：其中且
S205，计算二分图Q4-RI1：其中 且
S206，计算二分图QI-PI：GQI-PI＝(QI,PI,EQI-PI)，其中EQI-PI＝{e＝(v,w)|e∈E∧v∈PI∧w∈QI}；
S207，计算二分图Q3-PI：其中且
S208，计算二分图Q-P：GQ-P＝(Q,P,EQ-P)，其中EQ-P＝{e＝(v,w)|e∈E∧v∈P∧w∈Q}；
S209，计算匹配图BI-BI：GBI-BI＝(BI,EBI-BI)，其中EBI-BI＝{e＝(v,w)|e∈E∧v,w∈BI}；
S210，计算内核图Core：GCore＝(VCore,ECore)，其中VCore＝Q∪QI且ECore＝{e＝(v,w)|e∈E∧v,w∈VCore}；
S211，输出七个二分图GQI-PI、GQ-P，一个匹配图GBI-BI和一个内核图GCore，
其中，∈表示“属于”；表示“存在”；∧表示“逻辑与”；∪表示集合的并集；V和E分别表示互联网拓扑图G＝(V,E)的节点集和边集；P、Q、II、BI、QI、PI和RI分别为互联网拓扑结构模型分解得到的七类节点集。
4.根据权利要求1所述的一种互联网测试床拓扑结构大比例规模缩减方法，其特征在于，在S3中，所述调用二分图采样方法1包括如下步骤：
S301，输入二分图G＝(V1,V2,E)，其中V1和V2为两个节点集、E为边集，参数T、R，规模缩减率RR；转到步骤S302；
S302，将V2中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集其中这nh个节点中相同度重复数的最大值为T；转到步骤S303；
S303，计算并将二分图G分解为和其中且转到步骤S304；
S304，设d1,d2,…,ds为图G中节点集V1的全部节点度，并设dmax为节点集中的最大节点度，计算节点集的度分布{g(d)|d＝1,2,…,dmax}，其中g(d)表示中度为d的节点比率，初始化二分图G′＝(V1′,V2′,E′)←G＝(V1,V2,E)，初始化二分图初始化二分图初始化二分图GLL′中节点集的度分布{g′(d)|d＝1,2,…,dmax}←{g(d)|d＝1,2,…,dmax}，计算需删除边的总数Ed＝||E||·RR，初始化变量t＝0；转到步骤S305；
S305，如果t<Ed，转到步骤S306；否则二分图G′＝(V1′,V2′,E′)，其中且E′＝ELH′∪ELL′，为最终输出结果，转到步骤S311；
S306，如果任意选择集合S中一个节点v，并计算两个可选边集和如果转到步骤S308；如果转到步骤S309；否则且以概率PLL′＝γLL′/(γLH′+γLL′)转到步骤S308，并以概率1-PLL′转到步骤S309，其中，γLL′＝mLL-(||ELL||-||ELL′||)，γLH′＝mLH-(||ELH||-||ELH′||)，mLH＝Ed-mLL；转到步骤S307；
S307，如果S＝Φ，以式定义的离散概率分布pk/Σi＝1,2,…,spi(k＝1,2,…,s)随机采样一个节点度d∈{d1,d2,…,ds}，其中d1>d2>…>ds>0，f(d1),f(d2),…,f(ds)表示节点集V1中度为dk(1≤k≤s)的节点比率；计算两个可选边集 如果转到步骤S308；如果转到步骤S309，否则，以概率PLL′＝γLL′/(γLH′+γLL′)转到步骤S308，并以概率1-PLL′转到步骤S309；
S308，假设其中且计算dh＝argdmaxd∈D{g′(d)-g(d)}，其中D为集合中节点的度组成的集合，并计算的子集从Ea中均等概率分布地采样一条边e，并将该边从二分图和G′＝(V1′,V2′,E′)中删除，即更新ELL′←ELL′/{e}且E′←E′/{e}，更新GLL′中的节点度分布{g′(d)|d＝1,2,…,dmax}，其中g′(d)表示中度为d的节点比率，转到步骤S310；
S309，从中均等概率分布地采样一条边e，并将该边从二分图和G′＝(V1′,V2′,E′)中删除，即更新ELH′←ELH′/{e}且E′←E′/{e}，转到步骤S310；
S310，更新t←t+1，转到步骤S305；
S311，输出：采样子图G′＝(V1′,V2′,E′)；
其中，||S||表示集合S的势，即集合S包含元素的总数；x←y表示将右边的函数或数值y赋值给左边的变量x；dv表示节点v在图G′＝(V1′,V2′,E′)中的度；A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素；dh＝argdmaxd∈D{g′(d)-g(d)}表示将g′(d)-g(d)取最大值时参数d的取值赋值给dh；∈表示“属于”；表示“不属于”；∧表示“逻辑与”；节点的度是指与该节点相邻边的总数；∪表示集合的并集；输入参数规模缩减率RR被定义为需要从输入的二分图中删除边的总数与输入的二分图包含边数的比率，其中，nh的取值由节点集V2中节点的度和参数T唯一确定；Φ表示空集，即元素个数为零的集合。
5.根据权利要求1所述的一种互联网测试床拓扑结构大比例规模缩减方法，其特征在于，在S6中，所述调用二分图采样方法3包括如下步骤：
S601，输入：二分图G＝(V1,V2,E)，其中V1和V2为两个节点集、E为边集、节点集V2包含节点的度均为1，节点集参数T，规模缩减率RR；转到步骤S602；
S602，将V1中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集V1 H，其中这nh个节点中相同度重复数的最大值为T；转到步骤S603；
S603，计算V1 L＝V1/V1 H，设d1,d2,…,ds为二分图G中节点集V1的全部节点度，初始化二分图G′＝(V1′,V2′,E′)←G＝(V1,V2,E)，初始化节点集V1 L′←V1 L，初始化节点集V1 H′←V1 H，计算最终需要删除边的总数Ed＝||E||·RR，初始化变量t＝0；转到步骤S604；
S604，如果t<Ed，则以概率P＝γL′/(γH′+γL′)转到步骤S605，并以概率1-P转到步骤S608；否则，G′＝(V1′,V2′,E′)，其中V1′＝V1 H′∪V1 L′，为最终输出结果，转到步骤S610；其中，
且mH＝Ed-mL；
S605，如果其中则随机选择集合S中一个节点v，转到步骤S607，否则，转到步骤S606；
S606，以式定义的离散概率分布pk/Σi＝1,2,…,spi(k＝1,2,…,s)，随机采样一个节点度d∈{d1,d2,…,ds}，其中，d1>d2>…>ds>0，f(d1),f(d2),…,f(ds)表示节点集V1在二分图G＝(V1,V2,E)中度为dk(1≤k≤s)的节点比率，计算两个节点集和如果||V2||>0，则随机选择集合V2中一个节点v，否则随机选择集合V1中一个节点v，转到步骤S607；
S607，如果则随机选择一条与节点v相邻的边e，并更新二分图G′＝(V1′,V2′,E′)，E′←E′/{e}，如果dv＝0，则更新V1 L′←V1 L′/{v}，如果则此步不作任何操作；转到步骤S609；
S608，从集合中随机选择一条边e，并更新二分G′＝(V1′,V2′,E′)：E′←E′/{e}；如果dv＝0，则更新V1 H′←V1 H′/{v}；转到步骤S609；
S609，更新t←t+1，转到步骤S604；
S610，输出：采样子图G′＝(V1′,V2′,E′)；
其中，输入的节点集被包含于采样子图G′＝(V1′,V2′,E′)，即图G中节点集V2包含节点的度均为1，图G′中节点集V2′包含节点的度均为1；
其中，||S||表示集合S的势，即集合S包含元素的总数；x←y表示将右边的函数或数值y赋值给左边的变量x；当节点v∈V1或v∈V1 H或v∈V1 L，dv表示节点v在图G＝(V1,V2,E)的度，其他情况或没有特定约束下，dv表示节点v在图G′＝(V1′,V2′,E′)的度；A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素；∈表示“属于”；表示“不属于”；∧表示“逻辑与”；∨表示“逻辑或”；节点的度是指与该节点相邻边的总数；表示集合A被包含于集合B；∪表示集合的并集；∩表示集合的交集；输入参数规模缩减率RR被定义为需要从输入的二分图中删除边的总数与输入的二分图包含边数的比率。
6.根据权利要求1所述的一种互联网测试床拓扑结构大比例规模缩减方法，其特征在于，在S7中，所述调用二分图采样方法4包括如下步骤：
S701，输入：二分图G＝(V1,V2,E)，其中V1和V2为两个节点集、E为边集，多对一映射函数ψ:V1→W，设ψ(v1)为节点集W中与节点v1∈V1对应的唯一节点，对任意节点w∈W预选设定参数L(w)，预选给定节点度分布{g(d)|d＝1,2,…}，预先设定参数T、R，规模缩减率RR；转到步骤S702；
S702，将V2中全部节点按度从大到小顺利排列，并抽取前nh个最大度节点构成节点子集其中这nh个节点中相同度重复数的最大值为T；转到步骤S703；
S703，计算并将二分图G分解为和其中且转到步骤S704；
S704，设d1,d2,…,ds为二分图G中节点集V1的全部节点度，并设dmax为节点集中的最大节点度，初始化二分图G′＝(V1′,V2′,E′)←G＝(V1,V2,E)，初始化二分图初始化二分图计算二分图GLL′中节点集的度分布{g′(d)|d＝1,2,…,dmax}，其中，g′(d)表示中度为d的节点比率，计算最终需要删除边的总数Ed＝||E||·RR，其中，初始化变量t＝0，η＝0；转到步骤S705；
S705，如果t<Ed，转到步骤S706；否则二分图G′＝(V1′,V2′,E′)，其中且E′＝ELH′∪ELL′，为最终输出结果，转到步骤S713；
S706，如果则更新η←η-1，计算两个可选边集和
其中且集合N′(ψ(v1))＝{v|v∈V1′∧ψ(v)＝ψ(v1)}，
并转到步骤S709，否则，转到步骤S707；
S707，如果选择集合S中一个节点v，设该节点的度为d＝dv；否则以式定义的离散概率分布pk/Σi＝1,2,…,spi(k＝1,2,…,s)随机采样一个节点度d∈{d1,d2,…,ds}，其中d1>d2>…>ds>0，f(d1),f(d2),…,f(ds)表示节点集V1中度为dk(1≤k≤s)的节点比率，转到步骤S708；
S708，如果则更新η←η+1，转到步骤S712；如果则计算两个可选边集和转到步骤S709；如果d>1，则计算两个可选边集和转到步骤S709；
S709，如果转到步骤S710；如果转到步骤S711；如果 以概率PLL′＝γLL′/(γLH′+γLL′)转到步骤S710，以概率1-PLL′转到步骤S711，
其中，γLL′＝mLL-(||ELL||-||ELL′||)，γLH′＝mLH-(||ELH||-||ELH′||)，
mLH＝Ed-mLL；
S710，假设其中且计算其中D为集合中节点的度组成的集合；如果||Dh||＝0，则更新Dh←{argdmaxd∈D{g′(d)-g(d)}}；
首先，计算集合的一个子集：
假设其中且计算ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}，其中N′(ψ(v))＝{v1|v1∈V1′∧ψ(v1)＝ψ(v)}，
然后，计算集合Ea的一个子集，
最后，从集合Eb中随机选择一条边e，并将该边从二分图和G′＝(V1′,V2′,E′)中删除，即更新ELL′←ELL′/{e}且E′←E′/{e}，更新二分GLL′中的节点度分布{g′(d)|d＝1,2,…,dmax}，其中g′(d)表示中度为d的节点比率，转到步骤S712；
S711，假设其中且计算ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}，其中N′(ψ(v))＝{v1|v1∈V1′∧ψ(v1)＝ψ(v)}，
计算集合的一个子集：
并从集合Eb中随机选择一条边e，并将该边从二分图和G′＝(V1′,V2′,E′)中删除，即更新ELH′←ELH′/{e}且E′←E′/{e}，转到步骤S712；
S712，更新t←t+1，转到步骤S705；
S713，输出：采样子图G′＝(V1′,V2′,E′)；
其中，对集合N′(w)＝{v1|v1∈V1′∧ψ(v1)＝w}，其中W为步骤S701的输入，S701至S713期望对满足||N′(w)||接近于L(w)；
其中，||S||表示集合S的势，即集合S包含元素的总数；x←y表示将右边的函数或数值y赋值给左边的变量x；dv表示节点v在图G′＝(V1′,V2′,E′)中的度；A/B：如果A和B为数值，则表示除法，如果A和B为集合，则表示从集合A中删除集合B的所有元素；argdmaxd∈D{g′(d)-g(d)}表示g′(d)-g(d)取最大值时参数d的取值；ψ′＝argψ(v)maxv∈V{||N′(ψ(v))||/L(ψ(v))}表示将||N′(ψ(v))||/L(ψ(v))取最大值时参数ψ(v)的取值赋值给ψ′；∈表示“属于”；表示“不属于”；∧表示“逻辑与”；∪表示集合的并集；表示“任意”；节点的度是指与该节点相邻边的总数；输入参数规模缩减率RR被定义为需要从输入的二分图中删除边的总数与输入的二分图包含边数的比率。
7.一种互联网测试床拓扑结构大比例规模缩减装置，其特征在于，所述装置包括：
建模单元，用于获取真实世界互联网拓扑结构的探测数据，以简单无向图G＝(V,E)建模自治系统级互联网拓扑结构，其中集合V中节点的物理意义为真实世界互联网的自治系统域，集合E中边的物理意义为不同自治系统域之间的数据通信路径；
分解拓扑单元，用于采用互联网拓扑结构模型分解互联网拓扑图G＝(V,E)为七个二分图Q1-II、Q2-BI、QI1-RI、Q4-RI1、QI-PI、Q3-PI、Q-P，一个匹配图BI-BI和一个内核图Core；
二分图采样1单元，用于调用二分图采样方法1对二分图Q1-II、二分图Q2-BI和二分图QI1-RI进行采样，从而抽取部分节点和边构成采样子图；
匹配图采样单元，用于调用匹配图采样方法对匹配图BI-BI进行采样，抽取部分节点和边构成采样子图；
二分图采样2单元，用于调用二分图采样方法2对二分图Q4-RI1进行采样，抽取部分节点和边构成采样子图；
二分图采样3单元，用于调用二分图采样方法3对二分图QI-PI和二分图Q-P进行采样，抽取部分节点和边构成采样子图；
二分图采样4单元，用于调用二分图采样方法4对二分图Q3-PI进行采样，抽取部分节点和边构成采样子图；
子图合并单元，用于对得到的所有采样子图进行合并，得到互联网外围拓扑所有采样子图的合并图；
内核图采样单元，用于调用内核图采样方法对内核图进行采样，从而抽取部分节点和边构成采样子图；
采样子图获取单元，用于对得到的外围拓扑采样子图的合并图和内核拓扑的采样子图进行合并，从而得到真实世界互联网拓扑图的采样子图；
仿真测试单元，依据得到的真实世界互联网拓扑图的规模缩减采样子图，配置互联网测试床的拓扑连接关系，从而调度测试床有限的仿真节点资源，生成节点规模大比例缩减的互联网仿真测试环境。
CN201810565884.3A 2018-06-04 2018-06-04 一种互联网测试床拓扑结构大比例规模缩减方法及装置 Active CN108900320B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201810565884.3A CN108900320B (zh) 2018-06-04 2018-06-04 一种互联网测试床拓扑结构大比例规模缩减方法及装置
PCT/CN2018/124568 WO2019233089A1 (zh) 2018-06-04 2018-12-28 一种互联网测试床拓扑结构大比例规模缩减方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810565884.3A CN108900320B (zh) 2018-06-04 2018-06-04 一种互联网测试床拓扑结构大比例规模缩减方法及装置

Publications (2)

Publication Number Publication Date
CN108900320A true CN108900320A (zh) 2018-11-27
CN108900320B CN108900320B (zh) 2021-02-09

Family

ID=64344132

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810565884.3A Active CN108900320B (zh) 2018-06-04 2018-06-04 一种互联网测试床拓扑结构大比例规模缩减方法及装置

Country Status (2)

Country Link
CN (1) CN108900320B (zh)
WO (1) WO2019233089A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019233089A1 (zh) * 2018-06-04 2019-12-12 佛山科学技术学院 一种互联网测试床拓扑结构大比例规模缩减方法及装置
CN110708186A (zh) * 2019-09-09 2020-01-17 广东南方海岸科技服务有限公司 互联网测试床拓扑逼真度评价方法、装置、设备及介质
CN111597058A (zh) * 2020-04-17 2020-08-28 微梦创科网络科技（中国）有限公司 一种数据流处理方法及系统
CN112100797A (zh) * 2019-12-16 2020-12-18 佛山科学技术学院 基于结构模型的互联网拓扑仿真生成方法及系统
CN112104474A (zh) * 2019-12-16 2020-12-18 佛山科学技术学院 基于结构模型的互联网拓扑特征提取方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140281026A1 (en) * 2010-05-28 2014-09-18 Microsoft Corporation Virtual Data Center Allocation with Bandwidth Guarantees
CN104521160A (zh) * 2012-08-09 2015-04-15 微软公司 扩展型接入点
US20150381334A1 (en) * 2013-12-04 2015-12-31 Nec Laboratories America, Inc. Interference cancellation
CN106462990A (zh) * 2014-06-05 2017-02-22 微软技术许可有限责任公司 使用图形处理单元的可定制路线规划

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104516995A (zh) * 2013-09-29 2015-04-15 华为技术有限公司 网络流图缩减方法和装置
CN107844415B (zh) * 2017-09-28 2021-02-05 西安电子科技大学 一种基于插值的模型检测路径缩减方法、计算机
CN108900320B (zh) * 2018-06-04 2021-02-09 佛山科学技术学院 一种互联网测试床拓扑结构大比例规模缩减方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140281026A1 (en) * 2010-05-28 2014-09-18 Microsoft Corporation Virtual Data Center Allocation with Bandwidth Guarantees
CN104521160A (zh) * 2012-08-09 2015-04-15 微软公司 扩展型接入点
US20150381334A1 (en) * 2013-12-04 2015-12-31 Nec Laboratories America, Inc. Interference cancellation
CN106462990A (zh) * 2014-06-05 2017-02-22 微软技术许可有限责任公司 使用图形处理单元的可定制路线规划

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JINZHOU ZENG，FANGYUAN LI: "Distributed event-triggered bipartite consensus for multiple agents over signed graph topology", 《2015 34TH CHINESE CONTROL CONFERENCE》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019233089A1 (zh) * 2018-06-04 2019-12-12 佛山科学技术学院 一种互联网测试床拓扑结构大比例规模缩减方法及装置
CN110708186A (zh) * 2019-09-09 2020-01-17 广东南方海岸科技服务有限公司 互联网测试床拓扑逼真度评价方法、装置、设备及介质
CN110708186B (zh) * 2019-09-09 2022-07-12 广东南方海岸科技服务有限公司 互联网测试床拓扑逼真度评价方法、装置、设备及介质
CN112100797A (zh) * 2019-12-16 2020-12-18 佛山科学技术学院 基于结构模型的互联网拓扑仿真生成方法及系统
CN112104474A (zh) * 2019-12-16 2020-12-18 佛山科学技术学院 基于结构模型的互联网拓扑特征提取方法及系统
CN112104474B (zh) * 2019-12-16 2021-06-01 佛山科学技术学院 基于结构模型的互联网拓扑特征提取方法及系统
CN112100797B (zh) * 2019-12-16 2021-06-01 佛山科学技术学院 基于结构模型的互联网拓扑仿真生成方法及系统
CN111597058A (zh) * 2020-04-17 2020-08-28 微梦创科网络科技（中国）有限公司 一种数据流处理方法及系统
CN111597058B (zh) * 2020-04-17 2023-10-17 微梦创科网络科技（中国）有限公司 一种数据流处理方法及系统

Also Published As

Publication number Publication date
CN108900320B (zh) 2021-02-09
WO2019233089A1 (zh) 2019-12-12

Similar Documents

Publication Publication Date Title
CN108900320A (zh) 一种互联网测试床拓扑结构大比例规模缩减方法及装置
De Floriani et al. Morse complexes for shape segmentation and homological analysis: discrete models and algorithms
Comellas et al. Recursive graphs with small-world scale-free properties
CN103810321B (zh) 滤波器设计工具
CN110166312A (zh) 使用随机测试床自动创建网络装置模型
CN109948641A (zh) 异常群体识别方法及装置
CN102810113B (zh) 一种针对复杂网络的混合型聚类方法
Ahmed et al. Nested chinese restaurant franchise process: Applications to user tracking and document modeling
CN103218692B (zh) 基于活动间依赖关系分析的工作流挖掘方法
CN108667895A (zh) 基于二维码的文件共享方法、服务器及终端设备
CN103593393B (zh) 基于微博互动关系挖掘社交圈的方法及装置
CN107611962B (zh) 电网系统支路搜索方法、系统及电子设备
CN112087379A (zh) 业务链的编排方法及装置、存储介质和电子装置
CN108259637A (zh) 一种基于决策树的nat设备识别方法及装置
CN106452858A (zh) 一种识别网络用户的方法、装置和计算设备
CN108880871A (zh) 一种无线传感器网拓扑资源分配方法及装置
CN110188151A (zh) 一种基于cad和gis生成室内地图的方法及电子装置
CN108880872B (zh) 一种互联网测试床拓扑结构分解方法及装置
Budel et al. Random hyperbolic graphs in $ d+ 1$ dimensions
CN105005586A (zh) 一种基于度特征替换策略的流式图抽样方法
CN111709102B (zh) 一种基于层次聚类的供水管网分区方法
CN105242958A (zh) 一种虚拟试验系统与hla仿真系统数据交换方法
CN109063223A (zh) Bim模型的轻量化方法和装置以及bim模型的处理方法和系统
CN108683547A (zh) 一种软件定义的无线传感网络配置方法
CN115543428A (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