一种时变网络中时效稳定社团的检测方法和装置制造方法及图纸

本发明专利技术属于复杂网络分析技术领域，具体为一种时变网络中时效稳定社团的检测方法和装置。本发明专利技术方法包括：获取节点之间的连接随时间变化的情况，构建时变网络；使用波动率定量刻画网络中连边的动态变化程度；初始化社团结构，结合波动率计算时变网络的动态模块度；优化动态模块度，其最大值对应的社团结构即为时变网络的时效稳定社团。本发明专利技术使用波动率定量刻画网络动态变化，并将网络动态变化特点与社团检测方法结合起来，给出了识别时变网络中稳定社团的方法，提高了社团检测的准确性和可靠性，在社交网络、生物网络、交通网络等不同领域具有广泛的应用前景，同时为理解实际生活中复杂系统中社团结构的功能、动态过程提供新的视角。

【技术实现步骤摘要】
一种时变网络中时效稳定社团的检测方法和装置
本专利技术属于复杂网络分析
，具体涉及一种时变网络中的时效稳定社团检测方法和装置。
技术介绍
网络是刻画复杂系统中海量元素之间复杂关系的有效工具，实际生活中的很多系统都可以使用网络建模分析，如社交网络，蛋白质相互作用网络，交通网络等。基于网络科学的研究已经发现了实际系统许多重要的性质，包括小世界属性，富人俱乐部以及社团结构等。网络中的社团是指节点的聚类，相同社团的节点连接紧密，不同社团的节点连接稀疏。基于网络结构将节点划分为社团的过程被称为社团检测。社团结构广泛地存在于实际系统当中，如社会学、生物学、计算机科学、经济学、交通系统和电力系统等。网络的拓扑结构与节点属性、网络功能的联系通过社团得以体现，社团结构对网络中的传播过程产生影响。检测和研究社团结构有助于我们理解网络的结构特征和行为变化。同时，社团检测在实际生活中有丰富的应用场景，如线上购物网站中的商品推荐，社交软件中的广告投放，生物学研究中蛋白质和基因功能的预测等。由于社团结构的普遍存在和重要意义，近年来社团检测算法已成为网络科学的重要研究问题之一。图割和聚类等传统算法可以直接用于网络的社团检测，如层次聚类和谱聚类，其中层次聚类算法中最为著名的是Newman和Girvan提出的基于连边介数中心性的分裂算法。Newman和Girvan提出模块度(modularity)来衡量社团划分的质量，模块度值越大表明社团结构越明显。给定网络结构，最大化模块度是检测网络中社团结构的有效手段。由于这是一个NP难问题，一系列启发式算法被提出寻找近似解，如Blondel等人提出的Louvain算法适用于加权网络。基于模块度最优化的方法是应用最广泛的社团检测算法，但是该方法存在精度限制，无法找到相对网络而言尺寸较小的社团。除了直接分析网络的拓扑结构，还可以基于网络中的动态过程检测社团，包括自旋模型，随机游走和同步。已有的社团检测算法被广泛地应用到静态网络中，但是多数实际系统是随时间不断变化的。例如，当大脑受到外界信号的刺激时，脑区的活跃程度和不同脑区之间的联系会相应地改变；节假日期间，城市之间的交通流量会与平日不同。从网络拓扑的角度而言，系统的动态变化体现在节点的增加和减少，连边的出现和消失，以及连边权重的变化三个方面。为了保留网络结构在时间维度上的信息且模型的复杂度又不至于过高，多种关于时变网络的建模方式被提出。结构的动态变化和建模方式复杂度的增加为时变网络中的社团检测带来一系列挑战。关于时变网络中社团结构的研究能够为我们理解、控制和预测网络结构的演化、网络结构特征的出现和形成，网络中的传播过程提供更加深入和全面的视角。近年来，时变网络中的社团检测引起了学者的广泛关注。使用网络序列表示结构随时间的变化，并用静态网络的方法检测序列中每个网络的社团是一个简单直观的思路，这类方法被称为序列检测，根据是否考虑社团结构变化的平滑性可以分为独立检测和演化聚类两类方法。这类方法的关键问题在于排除噪声等因素对社团结构变化的干扰，以及确定相邻时间社团的对应关系。而基于多层网络的检测算法被称为整体检测方法，但目前层间连边局限在相邻网络的同一节点之间，非相邻网络不同节点之间的联系还有待研究。无论跟踪社团结构的变化还是寻找时变网络中不变的社团，对网络结构时效变化的研究还不够充分，已有的研究多数局限在对网络整体结构在短时间内连续变化这一特点的考虑，缺乏对网络在更长时间范围内的时效变化特征，网络中连边的时效变化，以及网络结构非连续变化的讨论。
技术实现思路
本专利技术目的在于提供一种时变网络中时效稳定社团的检测方法和装置，以解决现有社团检测方法局限于静态网络，忽视了实际系统动态变化特点的问题，从而提高社团检测的准确性，深化对实际系统时效变化的理解。本专利技术提供的时变网络中时效稳定社团的检测方法，具体步骤为：(1)获取一段时间内节点之间的连接随时间变化的情况，构建时变网络，使用邻接矩阵的集合表示各个时刻的网络结构；(2)基于步骤(1)得到的时变网络计算连边的波动率矩阵，定量刻画连边的动态变化程度；(3)初始化社团结构，结合波动率计算时变网络的动态模块度，用于定量刻画时变网络中的时效稳定社团；(4)优化所述步骤(3)的动态模块度；动态模块度最大值对应的社团划分即为时变网络中的时效稳定社团。本专利技术步骤(1)中，所述时变网络Gd表示N个节点之间的连边随时间变化的情况，具体地，Gd＝(G(t),t＝1…T)，其中G(t)＝(V,E(t))表示时刻t对应的网络结构，共T个时间点，V是N个节点构成的集合，各个时刻均相同，E(t)是t时刻连边构成的集合，随时间发生变化；将每个时刻的网络结构用邻接矩阵表示，即A(t)＝[aij(t)]N×N，如果在时刻t节点i和节点j之间不存在连边，则aij(t)＝0，否则aij(t)等于连边的权重，网络G(t)为无权无向网络，所以A(t)为对称矩阵。本专利技术步骤(2)中，所述定量刻画连边的动态变化程度，是使用波动率表示在一段时间内连边权重的波动程度，用于定量刻画网络的动态变化；节点i和j之间连边eij的波动率Vij按照如下方法计算：其中，std(·)表示标准差运算符，Δij(t)表示连边eij权重在时刻t的增长率；时变网络Gd中所有连边的波动率组成波动率矩阵V＝[Vij]N×N。本专利技术步骤(3)中，所述结合波动率计算时变网络的动态模块度，将通过动态模块度定量刻画时变网络中的时效稳定社团；时效稳定社团具体表现为：节点关于社团的隶属关系不随时间改变，所以不同时刻网络中对应的时效稳定社团结构是相同的，用向量C＝[C1,C2,…,CN]表示各个节点所属的社团；在不同时刻，同一个社团内节点之间的连接比不同社团节点之间的连接密切；不同社团之间连边的动态变化程度比同一社团内连边的变化更剧烈；动态模块度的计算方法如下：其中，si(t)是在时刻t网络中与节点i相连的连边权重之和，ω(t)为所有连边的权重之和，Ci表示节点i所属的社团，如果节点i和j属于同一个社团，则δ(Ci,Cj)＝1，否则为0，mt为t时刻网络中的连边数目；动态模块度是结合波动率，将Girvan-Newman模块度从静态网络向时变网络的拓展，用社团内连边的平均权重表示连接的密切程度，用连边的平均波动率表示连边的动态变化程度。本专利技术步骤(4)中，使用优化动态模块度的方法寻找时效稳定社团；当改变网络中的社团划分结构时，动态模块度发生变化；最大动态模块度对应的社团结构就是时变网络的时效稳定社团。本专利技术提供的时变网络中时效稳定社团的检测装置包括：时变网络构建模块、波动率计算模块、动态模块度优化模块、时效稳定社团输出模块，其中：所述时变网络构建模块，用于获取不同时刻节点之间的连接情况，构建包含N个节点，T个时间点的时变网络，并以邻接矩阵集合的形式表示网络在各个时刻拓扑结构；所述波动率计算模块，用于计算连边权重在一段时间内的波动率，并用波动率矩阵存储所有节点对之间连边的波动率；所述动态模块度优化模块，用于优化网络的动态模块度，寻找动态模块度的最大值；所述时效稳定社团输出模块，用于给出最大动态模块度对应的社团划分，即时变网络的时效稳定社团结构。本专利技术提供的技术方案具有如下优点：本专利技术首次将社团本文档来自技高网...

【技术保护点】
1.一种时变网络中时效稳定社团的检测方法，其特征在于，具体步骤为：(1)获取一段时间内节点之间的连接随时间变化的情况，构建时变网络，使用邻接矩阵的集合表示各个时刻的网络结构；(2)基于步骤(1)得到的时变网络计算连边的波动率矩阵，定量刻画连边的动态变化程度；(3)初始化社团结构，结合波动率计算时变网络的动态模块度；用于定量刻画时变网络中的时效稳定社团；(4)优化所述步骤(3)的动态模块度；动态模块度最大值对应的社团划分即为时变网络中的时效稳定社团。

【技术特征摘要】
1.一种时变网络中时效稳定社团的检测方法，其特征在于，具体步骤为：(1)获取一段时间内节点之间的连接随时间变化的情况，构建时变网络，使用邻接矩阵的集合表示各个时刻的网络结构；(2)基于步骤(1)得到的时变网络计算连边的波动率矩阵，定量刻画连边的动态变化程度；(3)初始化社团结构，结合波动率计算时变网络的动态模块度；用于定量刻画时变网络中的时效稳定社团；(4)优化所述步骤(3)的动态模块度；动态模块度最大值对应的社团划分即为时变网络中的时效稳定社团。2.根据权利要求1所述的方法，其特征在于，所述步骤(1)中的时变网络Gd表示N个节点之间的连边随时间变化的情况，具体地，Gd＝(G(t),t＝1…T)，其中G(t)＝(V,E(t))表示时刻t对应的网络结构，共T个时间点，V是N个节点构成的集合，各个时刻均相同，E(t)是t时刻连边构成的集合，随时间发生变化；将每个时刻的网络结构用邻接矩阵表示，即A(t)＝[aij(t)]N×N，如果在时刻t节点i和节点j之间不存在连边，则aij(t)＝0，否则aij(t)等于连边的权重，网络G(t)为无权无向网络，所以A(t)为对称矩阵。3.根据权利要求1所述的方法，其特征在于，所述步骤(2)使用波动率表示在一段时间内连边权重的波动程度，用于定量刻画网络的动态变化；节点i和j之间连边eij的波动率Vij按照如下方法计算：其中，std(·)表示标准差运算符，Δij(t)表示连边eij权重在时刻t的增长率；时变网络Gd中所有连边的波动率组成波动率矩阵V＝[Vij]N×N。4.根据权利要求1所述的方法，其特征在于，所述步骤(3)中通过动态模块度定量刻画时变网络中的时效稳定社团；时效稳定社团具体表现为...

【专利技术属性】
技术研发人员：李翔，王文婧，
申请(专利权)人：复旦大学，
类型：发明
国别省市：上海,31