基于网络嵌入的随机网络拓扑精确重构方法及系统技术方案

本发明专利技术公开了一种基于网络嵌入的随机网络拓扑精确重构方法，包括以下步骤：S1：将高维稀疏的随机网络图G的邻接矩阵M(G)通过特征嵌入表示为低维稠密特征矩阵Φ，其中邻接矩阵M(G)的维度是n

【技术实现步骤摘要】
基于网络嵌入的随机网络拓扑精确重构方法及系统


[0001]本专利技术属于移动通信
，主要涉及一种随机网络拓扑精确重构方法，具体是一种基于网络嵌入及矢量量化理论实现随机网络拓扑精确重构的方法及系统。

技术介绍

[0002]随机网络是一种拓扑结构随机变化的网络，其通信模式的消息转发机制为“存储
‑
携带
‑
转发”，具体来说网络中互相通信的节点在一定范围内随机移动，这些节点都能独立充当终端完成数据的存储转发任务，当两个具备通信能力的节点移动到对方的通信范围内后可以进行数据信息的传输，在无法找到下一个通信节点时，则暂时将数据存储在本地，继续随机移动直到进入其他节点通信范围内再进行消息转发。与传统网络协议MANET不同的是，这种模式无需两个通信的端点之间建立完整通讯路径，允许网络在某时刻处于非全连通的状态，更适合现实网络部署的需求，特别是对由于地震或军事打击导致基础设施被破坏的特殊场景非常具有应用意义。
[0003]随机网络的拓扑结包含网络中很多重要的信息，与网络系统的性能、可靠性及通讯代价等直接相关，因此及时获取精确的随机网络拓扑结构变得十分重要。由于随机网络拓扑结构具有动态变化性，通过路由协议获取网络的拓扑信息，可能前一时刻构建好了全局或者局域网络拓扑，由于网络中节点的移动，在下一时刻拓扑结构已经发生改变了，导致所获取的网络拓扑已经不是真实的网络拓扑，因此路由协议很难维护随机网络拓扑结构。随着对随机网络研究的逐渐深入，已经有了一些成果。
[0004]目前主流的拓扑重构方法为分布式和集中式。其中分布式估算需要实时获取每个节点的速度和方向等信息，需要占用大量通讯资源，从而导致网络拥塞。而在集中式目前针对随机网络拓扑的获取方法，主要是利用信号滤波技术实现随机网络拓扑邻接矩阵特征普参数的实时估算，然后通过公式1来重构网络邻接矩阵。假设随机网络在节点的移动过程中的某段时间内到达稳态，表示该时间切片下的网络快照邻接矩阵，现有的谱分解精确重构方法，是利用特征普参数进行重构：
[0005][0006]其中G代表当前时刻的随机网络，N为网络中节点的数目，特征普参数包括通过相关估算算法得到的网络邻接的特征值和对应的特征向量其中表示的转置。
[0007]在式1的网络精确重构方法中，当网络节点规模较大时，存在两个亟需解决的问题，1)首先其算法的时间复杂度为O(n3)，因此所需要的运算时间消耗难以保证重构网络拓扑的时效性；2)其次对于大规模随机网络全部特征谱参数的预测存在误差干扰，因此对于重构结果的精度有待提高。
[0008]为了解决或者至少部分解决当网络规模较大时，原有的网络精确重构算法计算时
间开销过大和重构精度较低的问题。

技术实现思路

[0009]针对现有随机网络拓扑恢复方法在大规模随机网络重构时，存在运行时间长重构精度低的问题，本专利技术提供一种可实现随机网络拓扑的快速精确重构的方法及系统。
[0010]本专利技术所采用的技术方案是：
[0011]提供一种基于网络嵌入的随机网络拓扑精确重构方法，包括以下步骤：
[0012]S1：将高维稀疏的随机网络图G的邻接矩阵M(G)通过特征嵌入表示为低维稠密特征矩阵Φ，其中邻接矩阵M(G)的维度是n
×
n，低维稠密特征矩阵Φ的维度为n
×
r，且r＜＜n；
[0013]S2：利用低维稠密特征矩阵初步重构邻接矩阵
[0014]S3：将邻接矩阵进行聚类量化，把离散的带有误差的矩阵元素映射到所属簇中心的准确值，消除重构扰动误差，得到随机网络精确重构后的邻接矩阵。
[0015]接上述技术方案，步骤S1具体包括以下步骤：
[0016]S11：假设高维稀疏的随机网络G的邻接矩阵M(G)的维度是n
×
n，在网络某时刻到达稳态后，获取其邻接矩阵的全部特征值λ
i
和对应的特征向量x
i
，按照特征值集合从大到小排列Eig{λ1，λ2…
λ
n
}，其对应的特征向量为{x1，x2…
x
n
}，采用前r个特征分量表示嵌入维度，则特征值矩阵Λ和特征向量矩阵X表示为Λ＝[λ1，λ2…
λ
r
]，X＝[x1，x2…
x
r
]；
[0017]S12：初始化一个低维矩阵并初始化一个酉矩阵Q＝I，用于获取低维嵌入目标矩阵Φ0＝BQ，其中目标矩阵Φ0的维度为n
×
r且r＜＜n，初始化求解矩阵为U＝rand(n，r)，表示一个维度为n
×
r的随机矩阵；
[0018]S13：定义损失函数为迭代求解U＝max(0，BQ)和Q＝VH
T
，其中H，V
T
为对U
T
B进行奇异值分解得到的酉矩阵U
T
B＝H∑V
T
，当满足迭代终止条件时，则求解矩阵U就是低维稠密特征矩阵Φ，ε为精确重构容忍误差上界。
[0019]接上述技术方案，步骤S2具体为：
[0020]利用得到的求解矩阵U初步重构邻接矩阵其中i，j∈1，2，...，n。
[0021]接上述技术方案，步骤S3具体包括以下步骤：
[0022]S31：根据初步重构邻接矩阵最大元素值m，将矩阵元素分为K个簇，设定划分阈值θ＝m/K，表示各簇中心的初始距离，先将邻接矩阵中的n2个元素放入初始集合S0中，并计算获得初始集合中距离最小的两个元素，并将两元素的中心值作为第一个簇中心，然后根据簇中心距离阈值θ，将该距离阈值范围内的所有元素放入集合S1，其余元素放入集合S2，即S2＝S0‑
S1，然后对S2重复该过程，直到初步选出K个簇中心c1，c2，...，c
k
；
[0023]S32：初始化簇中心C＝{c1，c2，...，c
k
}和隶属度矩阵U(0)＝μ
k
(a
ij
)，μ
k
为隶属度，设定聚类目标函数分别对μ
k
和c
k
求偏导，得到簇中心更新表达式和隶属度更新表达式通过迭代求解该表达式，直到簇中心值不发生变化或者元素分配集群不发生改变；
[0024]S33：按照最大隶属度原则，将元素分到不同的子簇，若μ
p
(a
ij
)＞μ
q
(a
ij
)，则将a
ij
分为第p个簇，p，q＝1，2，...，K
’
p≠q，但后用簇中心值对该簇内的元素进行赋值映射，得到最终精确重构的M(G)。
[0025]接上述技术方案，ε具体计算过程如下：
[0026]假设为重构网络矩阵，A为实际网络拓扑的邻接矩阵，则它们之间的误差通过如下公式计算：
[0027][0028]其中，假设重构的网络邻接本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于网络嵌入的随机网络拓扑精确重构方法，其特征在于，包括以下步骤：S1：将高维稀疏的随机网络图G的邻接矩阵M(G)通过特征嵌入表示为低维稠密特征矩阵Φ，其中邻接矩阵M(G)的维度是n
×
n，低维稠密特征矩阵Φ的维度为n
×
r，且r＜＜n；S2：利用低维稠密特征矩阵初步重构邻接矩阵，得到初步重构后的邻接矩阵S3：将初步重构后的邻接矩阵进行聚类量化，把离散的带有误差的矩阵元素映射到所属簇中心的准确值，消除重构扰动误差，得到随机网络精确重构后的邻接矩阵。2.根据权利要求1所述的基于网络嵌入的随机网络拓扑精确重构方法，其特征在于，步骤S1具体包括以下步骤：S11：假设高维稀疏的随机网络G的邻接矩阵M(G)的维度是n
×
n，在网络某时刻到达稳态后，获取其邻接矩阵的全部特征值λ
i
和对应的特征向量x
i
，按照特征值集合从大到小排列Eig{λ1，λ2…
λ
n
}，其对应的特征向量为{x1，x2…
x
n
}，采用前r个特征分量表示嵌入维度，则特征值矩阵Λ和特征向量矩阵X表示为Λ＝[λ1，λ2…
λ
r
]，X＝[x1，x2…
x
r
]；S12：初始化一个低维矩阵并初始化一个酉矩阵Q＝I，用于获取低维嵌入目标矩阵Φ0＝BQ，其中目标矩阵Φ0的维度为n
×
r且r＜＜n，初始化求解矩阵为U＝rand(n，r)，表示一个维度为n
×
r的随机矩阵；S13：定义损失函数为迭代求解U＝max(0，BQ)和Q＝VH
T
，其中H，V
T
为对U
T
B进行奇异值分解得到的酉矩阵U
T
B＝H∑V
T
，当满足迭代终止条件时，则求解矩阵U就是低维稠密特征矩阵Φ，ε为精确重构容忍误差上界。3.根据权利要求2所述的基于网络嵌入的随机网络拓扑精确重构方法，其特征在于，步骤S2具体为：利用得到的求解矩阵U初步重构邻接矩阵其中i，j∈1，2，...，n。4.根据权利要求2所述的基于网络嵌入的随机网络拓扑精确重构方法，其特征在于，步骤S3具体包括以下步骤：S31：根据初步重构后的邻接矩阵最大元素值m，将矩阵元素分为K个簇，设定划分阈值θ＝m/K，表示各簇中心的初始距离，先将邻接矩阵中的n2个元素放入初始集合S0中，并计算获得初始集合中距离最小的两个元素，并将两元素的中心值作为第一个簇中心，然后根据簇中心距离阈值θ，将该距离阈值范围内的所有元素放入集合S1，其余元素放入集合S2，即S2＝S0‑
S1，然后对S2重复该过程，直到初步选出K个簇中心c1，c2，...，c
k
；S32：初始化簇中心C＝{c1，c2，...，c
k
}和隶属度矩阵U(0)＝μ
k
(a
ij
)，μ
k
为隶属度，设定聚类目标函数f分别对μ
k
和c

【专利技术属性】
技术研发人员：颜昕，郭款，雷中杰，
申请(专利权)人：武汉理工大学，
类型：发明
国别省市：