CN113271643A - 一种多节点协作无人机自组网分簇拓扑重构方法 - Google Patents

Abstract

本发明公开了一种多节点协作无人机自组网分簇拓扑重构方法。该方法首先采用最小节点ID法进行快速分簇，并根据运动的无人机自组网中节点的速度相似度、节点间平均距离、平均链路维持时间、无人机节点剩余能量等因素计算每个无人机节点的权值，进行簇的维护。然后在此基础上提出了一种基于邻域内邻接矩阵的网络中关键节点判断方法，分布式地判断网络中是否有关键节点失效。最后在该分簇算法与关键节点判断方法的基础上针对网络遭受敌方攻击而出现拓扑损坏的问题提出了两阶段拓扑修复方法。无线网络仿真环境EXata中仿真结果表明，与传统分簇算法相比，该方法可以快速、有效地修复损坏的网络拓扑，提升网络的吞吐量。

Description

一种多节点协作无人机自组网分簇拓扑重构方法

技术领域

本发明属于无线网络领域，特别涉及多节点协作无人机自组网分簇拓扑重 构方法。

背景技术

无人机自组织网络具有分布式、自组织、多业务和远距离等优势，被广泛 应用于军事和民用领域。在实际应用过程中，拓扑高动态变化、战场环境复杂 恶劣的特点对无人机自组网的通信质量造成了严重威胁。利用网络中空闲或冗 余节点进行协作传输的方法可以有效提升网络的通信性能，这种方法被称为节 点协作策略。传统的分簇算法在一定程度上解决了无人机自组网的规模扩展性 问题，然而在实际战场环境中，无人机自组织网络具有节点高速移动、拓扑动 态变化的特点，同时也极易遭受敌方的火力攻击，这些都对分簇网络的拓扑重 构问题提出了新的挑战。

附图1给出了一种具体的分簇算法的流程图，其具体步骤如下：

步骤1：节点首先随机退避一段时间，随后发送簇首竞争帧。发送完毕后， 节点将簇首竞争帧中的发送节点ID号设为自身当前簇首，并设置定时器等待 WF_CLUSTER时间，等待该段时间是考虑到网络中有ID号更小的节点存在且 竞选簇首的情况。如果该定时器超时，则说明该节点的一跳邻居范围内没有ID 号更小的节点了，则该节点成功竞选成为簇首，退避一段时间后向周围所有节 点发送簇首广播帧，告知其他节点自身的位置、速度等信息。

步骤2：若该节点在WF_CLUSTER定时器超时之前收到了其他节点广播的 簇首竞选帧，则会比较收到的簇首竞选帧中的节点ID号与自身的大小，如果接 收到的簇首竞选帧中的ID号更小，则将该簇首竞选帧的发送节点设为自身的簇 首节点，重新设置WF_CLUSTER定时器，进入同样的等待过程，当 WF_CLUSTER定时器超时后，节点加入以该节点为簇首的簇。

步骤3：对于在退避阶段准备发送簇首竞选帧的节点，如果在退避完成之 前就收到了其他节点广播的簇首竞选帧，首先会挂起自身的退避定时器，随后 判断帧头部字段中的节点ID号是否比自身要小，如果更小，则取消退避定时器， 并将该簇首竞选帧发送节点设置为簇首，设置定时器等待WF_CLUSTER时间。 如果自身的节点ID号更小，则会忽略该簇首竞选帧，继续进行退避参与竞选簇 首。

随着无人机自组织网络技术的不断发展，人们对于网络拓扑连通性的需求 越来越高，如何设计适合无人机自组织网络的拓扑重构策略具有重要意义。为 了在通信距离远、节点高速移动、拓扑动态变化的无人机自组织网络中维持拓 扑的连通性，本文提出了一种多节点协作无人机自组网分簇拓扑重构方法。

发明内容

本发明的目的是针对无人机自组织网络环境，提出一种多节点协作无人机 自组网分簇拓扑重构方法。为了实现该目的，本发明所采用的步骤是：

步骤1：采用最小节点ID法进行快速分簇，分簇完成后簇首增大发射功率 与邻居簇首进行通信，交换信息。

步骤2：根据运动的无人机自组网中节点的速度相似度、节点间平均距离、 平均链路维持时间、无人机节点剩余能量等因素计算每个无人机节点的权值， 进行簇的维护。

步骤3：提出了一种基于邻域内邻接矩阵的网络中关键节点判断方法，分布 式地判断网络中是否有关键节点失效。

步骤4：在分簇算法与关键节点判断方法的基础上针对网络遭受敌方攻击而 出现拓扑损坏的问题提出了两阶段拓扑修复策略。

本发明提出的基于动态中继选择的无人机自组网协作时分信道接入方法已 在EXata网络仿真环境下得到了验证。仿真场景为边长为20000m的正方形区域， 其中初始位置分布有36个节点，朝着目标战略位置飞行。仿真时长为80秒。 数据包长度固定为2000字节，信道传输速率为20Mbps，传输层采用UDP协议。

附图说明

图1是本发明采用的最小节点ID号方法初始化分簇流程图；

图2是分簇完成后的网络拓扑示意图；

图3是有向图示例图与邻域内邻接矩阵示例图；

图4是二次竞选簇首阶段流程图；

图5是两阶段簇首修复流程关系图；

图6是簇间移动控制拓扑修复方案流程图；

图7是两种节点失效情况下吞吐量随发包速率变化图；

图8是两种节点失效情况下路由协议开销随发包速率变化图。

具体实施方式

下面结合附图和实例对本发明作进一步详细描述。

本发明提出的多节点协作无人机自组网分簇拓扑重构方法已在EXata网络 仿真环境下得到了验证，并通过仿真结果证明了该方法的有效性。下面给出本 发明的具体实施步骤：

步骤1：采用最小节点ID法进行快速分簇，分簇完成后簇首增大发射功率 与邻居簇首进行通信并交换彼此信息。

当进入网络运行阶段，由于节点刚刚初始化配置完成，对于周围节点的信 息掌握不充分，所以先采用最小ID号分簇算法进行分簇。为防止冲突，节点首 先随机退避一段时间，随后发送簇首竞争帧。发送完毕后，节点将簇首竞争帧 中的发送节点ID号设为自身当前簇首，并设置定时器等待WF_CLUSTER时间， 等待该段时间是考虑到网络中有ID号更小的节点存在且竞选簇首的情况。如果 该定时器超时，则说明该节点的一跳邻居范围内没有ID号更小的节点了，则该 节点成功竞选成为簇首，退避一段时间后向周围所有节点发送簇首广播帧，告 知其他节点自身的位置、速度等信息。

若该节点在WF_CLUSTER定时器超时之前收到了其他节点广播的簇首竞 选帧，则会比较收到的簇首竞选帧中的节点ID号与自身的大小，如果接收到的 簇首竞选帧中的ID号更小，则将该簇首竞选帧的发送节点设为自身的簇首节点， 重新设置WF_CLUSTER定时器，进入同样的等待过程，当WF_CLUSTER定时 器超时后，节点加入以该节点为簇首的簇。

同样的，对于在退避阶段准备发送簇首竞选帧的节点，如果在退避完成之 前就收到了其他节点广播的簇首竞选帧，首先会挂起自身的退避定时器，随后 判断帧头部字段中的节点ID号是否比自身要小，如果更小，则取消退避定时器， 并将该簇首竞选帧发送节点设置为簇首，设置定时器等待WF_CLUSTER时间。 如果自身的节点ID号更小，则会忽略该簇首竞选帧，继续进行退避参与竞选簇 首。

在初始化阶段采用最小节点ID号方法完成初步分簇之后，簇首节点周期性 向簇内节点发送簇首广播帧，簇首广播帧内携带发送节点号、帧类型、簇首节 点位置与速度等信息。簇首节点与簇内节点的运动方法采用领航者-跟随者法， 簇首节点以区域中的目的地点为目标，受到来自目的地点的导航力的作用而移 动，簇内节点在建簇完成后，从收到的簇首广播帧中获得簇首节点的速度矢量， 将自身的速度大小、方向保持与簇首节点一致，一起向目标战略位置移动并执 行任务，初始分簇完成后的拓扑示例如附图2所示。

步骤2：根据运动的无人机自组网中节点的速度相似度、节点间平均距离、 平均链路维持时间、无人机节点剩余能量等因素计算每个无人机节点的权值， 进行簇的维护。

在初始化分簇完成后，随着网络的运行与信息的交换，网络中的节点将会 了解自身与邻居节点的距离、邻居节点的速度、方向等信息。在后续的分簇维 护过程中，将利用加权分簇的思想，综合考虑节点的速度相似度与邻居节点的 平均距离、平均通信链路维持时间、节点的剩余能量等因素对每个节点计算自 身的权值W_i：

其中W_V，W_D，W_T，W_P为权值系数，满足W_V+W_D+W_T+W_P＝1，具体大小 根据实际情况进行取值。V_i为节点的速度相似度，D_i为与邻居节点平均距离的 度，T_i为平均通信链路维持时间，P_i为节点的剩余能量。接下来对公式(4.1) 中的节点速度相似度、与邻居节点的平均距离、平均通信链路维持时间进行理 论推导。

对于节点i的速度相似度V_i，设节点i与一个邻居节点j的速度标量分别为 v_i与v_j，v_i与v_j与水平轴的夹角分别为θ_i与θ_j，则节点i与节点j在水平轴(即x轴) 上的速度分量差为：

v_ijx＝v_icosθ_i-v_jcosθ_j (2)

节点i与节点j在纵轴(即y轴)上的速度分量差为：

v_ijy＝v_isinθ_i-v_jsinθ_j (3)

假设节点i有N个邻居节点，则节点i与周围邻居节点在x轴上的平均速度 差为：

节点i与周围邻居节点在y轴上的平均速度差为：

则节点i与周围邻居节点的平均速度差为：

根据上述结果，可以分别求出节点i与周围邻居节点在x轴与y轴上平均速 度差的方差与标准差。节点i与周围邻居节点在x轴上的平均速度差的标准差为：

节点i与周围邻居节点在y轴上的平均速度差的标准差为：

根据所求的节点i与周围邻居节点的平均速度差的标准差，本章将节点i与 周围邻居节点的速度相似度V_i定义为：

其中σ_i表示为节点i与周围邻居节点在x轴与y轴方向速度差的方差的平均 值的开方，即：

对于节点i与邻居节点平均距离的度D_i，可以求出节点i与周围邻居节点的 平均距离为：

根据节点i与周围邻居节点的平均距离可求得距离的标准差，即：

根据节点i与周围邻居节点的平均距离及距离的标准差，将节点i与邻居节 点平均距离的度D_i定义为：

接下来求节点i与周围邻居节点的平均链路通信维持时间。先求单条链路通 信维持时间，以节点i与其一邻居节点j为例进行说明，设节点i所在位置为点 I，节点j所在位置为J，节点I的通信半径为R，节点J与节点I的相对速度为 v_ji，v_ji的大小可由公式(2)与(3)求得，即：

则节点i与节点j的链路生存时间可以表示为：

根据所求的T_ij即可求出节点i与周围邻居节点的平均通信链路生存时间为：

将所求的值代入公式(1)，每个节点即可计算出自身的权值W_i。在后续簇 维护过程中，簇首节点会收集簇内成员节点的权值，每一个簇首选举周期如果 簇首节点没有失效，则会周期性地按照权值大小更新簇首，以保持簇的稳定性。 如果簇首节点失效，则簇内节点会设定退避计时器进行退避，退避完成后首先 发送一次入簇请求帧，随后设置定时器等待一段WF_CLUSTER_RTJ时间，如 果定时器超时，则表明节点一跳范围内没有簇首，节点开始竞选簇首。

为提高网络覆盖率，在分簇完成后，竞选成为簇首的节点在发送数据帧时 会增大发射功率，增大发射功率的值根据网络拓扑需要设定，簇内成员节点依 旧采用默认发射功率。

步骤3：提出了一种基于邻域内邻接矩阵的网络中关键节点判断方法，分布 式地判断网络中是否有关键节点失效。

设节点v的邻居节点集合N_max(v)＝{v₁，v₂，...，v_n|d(v，v_i)＜r(v)，i＝1，2，...，n}。其中d(v，v_i)表示节点v与邻居节点v_i的距离，r(v)表示节点v的通信半径。 设节点v的一个邻域节点为v_u，建立集合X(v_u)＝Φ，遍历N_max(v)中除v 与u之外的节点v_i，将a_ui＝1(即节点u与节点i可达)对应的节点v_i加入集合 X(v_u)。随后对X(v_u)进行判断，如果X(v_u)为空集，即说明节点u与节点v 的其他邻居节点都不连通，即如果节点v失效，会造成网络的分割。随后按照 此流程对集合N_max(v)中的其他节点做如上的操作，判断节点v是否为关键节 点。采用该思想，即可以求出节点v的邻域内邻接矩阵，若邻域内接矩阵中有 一行全为0，则节点v是关键节点。

例如附图3(a)所示的拓扑，对于节点5，可以到达的节点有节点2、节点3、 节点4、节点6、节点7，即N_max(v₅)＝{v₂，v₃，v₄，v₆，v₇}，节点5的邻域拓扑为 R_max(v₅)＝G(N_max(v₅)，E_max(v₅))。

对于节点5的邻居节点集(即邻域)中的节点，按照上述的流程进行遍历， 节点2可以到达邻居节点集中的节点3与节点4，节点3可以到达节点2，节点 4可以到达节点6，节点6可以到达节点4，节点7无法到达其它节点。则节点 5的邻域内邻接矩阵可以用图3(b)表示：

由节点5的邻域内邻接矩阵可以看出，节点7所在行均为0，即表示节点5 为关键节点，节点5的失效会节点7与节点5的其它邻节点无法连通。

根据上述图论的思想，网络中的簇首节点即可以根据两跳邻居节点信息结 合有向图与邻域内邻接矩阵分布式地判断出自身是否为关键节点。具体步骤如 下：

网络中的簇首节点周期性地发送簇首广播帧，簇首广播帧中携带节点位置、 速度、邻居节点信息等信息。

一跳范围内的其它节点接收到该簇首广播帧后，首先判断自身是否为簇内 成员节点或者其它簇首节点，若满足条件则接收该簇首广播帧。簇内成员节点 更新簇首有效时间与位置、速度等信息，随后丢弃该簇首广播帧以释放内存空 间。接收该簇首广播帧的其它簇首节点根据该簇首广播帧中节点携带的邻居节 点表与关键节点判断位信息更新自身的邻接矩阵，并将该节点信息存入自身维 护的邻居节点数据表中，当节点更新邻接矩阵与信息完毕后，同样会丢弃掉该 簇首广播帧。

随后，在该接收簇首节点自身发送簇首广播帧时，会启动关键节点判断流 程，如果邻域内邻接矩阵中存在有一行全为0的情况，即说明该节点为关键节 点，将簇首广播帧头部的关键节点标志位信息置为1，随后向周围邻居簇首节点 进行广播。全网簇首节点都按照这样的流程不断更新周围簇首与自身的信息， 从而可以实现分布式的全网关键簇首节点判断。

以上介绍了簇首节点之间进行失效簇首节点的方法。若失效的簇首节点范 围内还有其他在该簇首失效前维护的簇内成员节点，则簇内成员也会通过簇首 失效检测机制检测到簇首失效。下面将对簇内成员的簇首失效检测机制进行介 绍。

簇内成员节点在每一次接收到簇首广播帧后会根据接收时间设定一个定时 器，定时器超时的时间为当前时间加上簇首失效检测时间间隔，当定时器超时 即可判定簇首节点已经失效，随后根据失效簇首节点信息与所计算的权值进行 新簇首的选举与调整。如果在该定时器超时前收到了新的簇首广播帧，则表明 簇首节点依然有效，节点更新本地维护的簇首节点信息，并取消之前的定时器， 设置新的定时器，继续等待下一个簇首广播帧。

步骤4：在分簇算法与关键节点判断方法的基础上针对网络遭受敌方攻击而 出现拓扑损坏的问题提出了两阶段拓扑修复策略。

如果失效簇首范围内存在原簇内成员，则簇内成员节点会通过簇首失效检 测机制发现该簇首失效，随后通过之前所述加权分簇的方法产生一个新的簇首， 完成拓扑修复，竞选成功的簇首会广播簇首修复告知帧，告知周围的簇首节点 原簇首失效的情况已被修复，并且根据原簇首节点的位置与速度信息进行位置 的调整。周围簇首节点收到该簇首修复告知帧后即知道原失效节点被修复，不 会启动第二阶段拓扑修复过程。

当簇内节点检测到簇首节点失效后会启动簇内修复流程，进入二次竞选簇 首阶段，选举出一个权值最小的节点作为新的簇首，完成拓扑的修复，如附图4 所示。节点首先设置定时器，随机退避一段时间，如果退避定时器超时，即顺 利完成退避，则节点将自身计算的权值添加到帧头部，广播簇首竞选帧，发送 完毕后，节点将自身节点ID号设为当前簇首，随后设置定时器，等待一段 WF_CLUSTER时间，如果在定时器超时前未收到其它簇首竞选帧，则节点判定 自身竞选簇首成功，随后计算原簇首未失效情况下当前时刻的位置，并进行调 整。如果节点在退避完成前收到了其它节点的簇首竞选帧，则会首先查看该簇 首竞选帧中携带的权值数值，并与自身进行比较，如果该簇首竞选帧中权值比 自身要大，则节点将该簇首竞选帧中携带的源节点地址设为自身的簇首节点地 址，重新设置定时器等待WF_CLUSTER时间；否则节点将判断权值是否相等， 如果权值相等，则进一步比较节点ID号，如果收到的簇首竞选帧中的ID号更 小，则将簇首竞选帧中的源节点地址设置为自身的簇首节点地址，重新设置定 时器等待WF_CLUSTER时间，相反，如果自身的ID号更小，则继续进行退避。 如果权值不相等，则显然只有自身权值更大这一种情况，节点会丢弃该簇首竞选帧，继续进行退避。如果经历了这一阶段后节点将簇首更改为收到的簇首竞 选帧中的源节点并进入等待WF_CLUSTER状态，则如果定时器超时则重新分 簇成功，簇首节点为新收到的簇首竞选帧中的源节点；如果在定时器超时前又 收到了新的簇首竞选帧，则重新对簇首竞选帧中的权值进行判断，如果簇首竞 选帧中的权值大于当前节点维护的簇首的权值，则节点将该簇首竞选帧的源节 点更新为簇首，重新设置定时器等待WF_CLUSTER时间；否则按照上述判断 流程，如果权值相等但自身节点号更小或者权值不相等(即自身权值更大)，则丢弃收到的簇首竞选帧，继续等待之前的定时器超时完成分簇。否则就将该簇 首竞选帧的源节点更新为簇首，重新设置定时器等待WF_CLUSTER时间。在 定时器超时前未收到新的簇首竞选帧则判定自身入簇成功。

当新的簇首节点产生后，会广播簇首修复告知帧，告知周围的簇首节点原 失效簇已被修复，同时簇首节点会根据失效簇首节点的位置与速度信息进行位 置的调整。

周围簇首节点在原簇首节点失效后会设置定时器等待一段 WF_FirstStage_Reconstruction时间，如果该定时器超时前收到该失效簇首的簇首 修复告知帧，则更新新的簇首信息后丢弃该帧；如果该定时器超时且并没有收 到簇首修复告知帧，即表明关键簇整体失效，启动第二阶段拓扑修复流程，即 簇间移动控制拓扑修复方案，如附图5所示。

当簇首节点检测到邻居簇首节点失效，并且在WF_FirstStage_Reconstruction时间后仍未收到新簇首修复告知帧，则首先统计失效簇首节点的邻居簇首节点 数量，随后根据失效簇首邻居节点数量与自身簇内成员数量计算协作修复概率， 根据该概率发送拓扑修复请求帧，请求簇内成员节点进行拓扑修复。簇内成员 节点收到该拓扑修复请求帧后，会随机退避一段时间并向簇首节点回复修复应 答帧，在修复应答帧中包括了成员节点当前的位置信息。簇首节点根据簇内成 员节点的位置信息与剩余能量综合选择一个节点进行拓扑修复，并向该节点回 复修复确认帧，该节点收到后即知道自身将承担修复任务，随即开始拓扑修复 进程，如附图6所示。

簇内节点接收到应答确认帧，得知自身被选择参与失效关键簇首节点修复， 即进入簇内成员节点移动修复过程。对于修复节点来说，首先需要确定失效簇 首节点的位置(即如果关键簇首节点未失效，此时应运动到的位置)，修复节点 可以从修复请求帧中获取失效关键簇首节点失效前的位置，用历史位置来推测 当前位置，也可以通过关键簇首节点失效前的速度、方向等信息对当前时刻的 预期位置进行计算，为了获得较为准确的失效关键簇首位置，修复节点采用卡 尔曼滤波法对这两种估计结果进行加权计算随后在虚拟力作用下向确定位置进 行移动。

当修复节点获得较为准确的失效关键簇首位置之后，便会在虚拟力的作用 下向着该位置移动，假设该目的地点位置坐标为D(x_d，y_d，z_d)，修复节点i坐标为 I(x_i，y_i，z_i)，则修复节点受到的吸引力可表示为：

其中M_a为常数，表示吸引力系数，

的方向与节点i到目的位置的方向相 同，且大小与节点i和目的位置之间的距离正相关。在修复节点向目的位置移动 的过程中，还可能因为与周围节点距离过近而受到排斥力，设某一时刻节点i 的邻居节点的集合为N，假设集合N中一个节点为j，坐标为J(x_j，y_j，z_j)，则节点 i受到节点j的排斥力为：

其中M_r为常数，表示排斥力系数，

的方向与节点i到邻居节点的方向相 反，且大小与节点i到邻居节点之间的距离正相关。结合修复无人机受到的来自 障碍物与边界的排斥力，无人机i受到的合力可以表示为：

其中O为无人机节点i附近的障碍物与边界的集合。在按照上述流程计算 出合力后，进一步判断该计算值是否为0，若为0则对节点施加随机扰动，让节 点跳出当前状态。若该计算出的值不为0，则继续比较当前矢量方向与之前的方 向是否相反，若相反，则修复节点会出现往复移动的情况，同样需要施加一次 随机扰动力。若不相反，则节点按照合力的方向进行移动，直至移动到待修复 范围内，完成拓扑修复任务。综上，修复节点受到的合力可表示为：

修复节点的运动方程可表示为：

其中P_i(t)表示t时刻节点坐标，P_i(t+Δt)表示t+Δt时刻节点坐标，F_i/|F_i|表 示当前时刻前进方向。

当修复节点到达修复位置后，为了防止有其它节点已经提前到达修复位置 并且建簇成功，会首先广播一次入簇请求帧，随后设置定时器等待一段 WF_Cluster_RTJ时间，定时器超时后即表明修复节点邻居范围内没有已存在的 簇，随后进入竞选簇首流程，广播一次簇首竞选帧并设置定时器等待一段 WF_Cluster时间，定时器超时后即表明建簇成功，发送簇首广播帧告知周围邻 居簇首节点，拓扑修复完成。如果修复节点在定时器超时前收到了其它簇首回 复的入簇应答帧，则表明已经有其它簇派出的修复节点到达了失效簇首位置， 并且成功建簇，修复节点加入该簇，成为簇内成员节点。

随后，簇内成员节点之间采取基于虚拟力导向的控制方法自主地实现局部 拓扑的优化，节点受到来自簇内其它成员节点之间的吸引力与排斥力、与障碍 物和边界之间的避障力、与新簇首之间的吸引力。修复节点i所受的合力可以表 示为：

其中N表示节点i的邻居节点集，O表示障碍物指示位置与边界的集合，

表示节点受到簇首节点的吸引力。节点可由该式计算出自身受力大小与方向， 从而可以实现动态的分布式的拓扑优化。簇内成员节点在该虚拟力的作用下进 入平衡状态后，根据簇首广播帧中的速度信息调整自身的速度，与簇首节点保 持一致，继续向战略目标位置移动。

在附图7中描述了两种节点失效情况下不同方法的吞吐量与发包速率的关 系，其中方法2为本发明提出的方法。由附图7可以看出，本发明提出的两阶 段拓扑重构方法网络吞吐量均高于不具有重构功能的方法1，其中第一阶段重构 方法可以迅速修复网络拓扑，可获得最高的吞吐量，第二阶段重构方法由于需 要节点移动进行簇间修复，所以吞吐量低于第一阶段重构方法，但仍高于方法1。

在附图8中描述了两种节点失效情况下不同方法的吞吐量与发包速率的关 系，其中方法2为本发明提出的方法。由附图8可以看出，本发明提出的两阶 段拓扑重构方法路由开销均低于不具有重构功能的方法1，其中第一阶段重构方 法可以迅速修复网络拓扑，使网络恢复连通，因此具有最小的路由开销。第二 阶段重构方法由于需要节点移动进行簇间修复，所需时间相比第一阶段更长， 因此需要更多的路由控制帧，路由开销大于第一阶段，但仍低于方法1。

本发明申请书中未作详细描述的内容属于本领域专业技术人员公知的现有 技术。

Claims (5)

1.一种多节点协作无人机自组网分簇拓扑重构方法，所采用的步骤是：

步骤1：采用最小节点ID法进行快速分簇，分簇完成后簇首增大发射功率与邻居簇首进行通信并交换彼此信息。

步骤2：根据运动的无人机自组网中节点的速度相似度、节点间平均距离、平均链路维持时间、无人机节点剩余能量等因素计算每个无人机节点的权值，进行簇的维护。

步骤3：提出了一种基于邻域内邻接矩阵的网络中关键节点判断方法，分布式地判断网络中是否有关键节点失效。

步骤4：在分簇算法与关键节点判断方法的基础上针对网络遭受敌方攻击而出现拓扑损坏的问题提出了两阶段拓扑修复策略。

2.根据权利要求1所述的采用最小节点ID法进行快速分簇，分簇完成后簇首增大发射功率，其特征在于：网络中节点根据一跳范围内所有节点ID号的大小竞选簇首的方法为：

当进入网络运行阶段，由于节点刚刚初始化配置完成，对于周围节点的信息掌握不充分，所以先采用最小ID号分簇算法进行分簇。为防止冲突，节点首先随机退避一段时间，随后发送簇首竞争帧。发送完毕后，节点将簇首竞争帧中的发送节点ID号设为自身当前簇首，并设置定时器等待WF_CLUSTER时间，等待该段时间是考虑到网络中有ID号更小的节点存在且竞选簇首的情况。如果该定时器超时，则说明该节点的一跳邻居范围内没有ID号更小的节点了，则该节点成功竞选成为簇首，退避一段时间后向周围所有节点发送簇首广播帧，告知其他节点自身的位置、速度等信息。

若该节点在WF_CLUSTER定时器超时之前收到了其他节点广播的簇首竞选帧，则会比较收到的簇首竞选帧中的节点ID号与自身的大小，如果接收到的簇首竞选帧中的ID号更小，则将该簇首竞选帧的发送节点设为自身的簇首节点，重新设置WF_CLUSTER定时器，进入同样的等待过程，当WF_CLUSTER定时器超时后，节点加入以该节点为簇首的簇。

同样的，对于在退避阶段准备发送簇首竞选帧的节点，如果在退避完成之前就收到了其他节点广播的簇首竞选帧，首先会挂起自身的退避定时器，随后判断帧头部字段中的节点ID号是否比自身要小，如果更小，则取消退避定时器，并将该簇首竞选帧发送节点设置为簇首，设置定时器等待WF_CLUSTER时间。如果自身的节点ID号更小，则会忽略该簇首竞选帧，继续进行退避参与竞选簇首。

在初始化阶段采用最小节点ID号方法完成初步分簇之后，簇首节点周期性向簇内节点发送簇首广播帧，簇首广播帧内携带发送节点号、帧类型、簇首节点位置与速度等信息。簇首节点与簇内节点的运动方法采用领航者-跟随者法，簇首节点以区域中的目的地点为目标，受到来自目的地点的导航力的作用而移动，簇内节点在建簇完成后，从收到的簇首广播帧中获得簇首节点的速度矢量，将自身的速度大小、方向保持与簇首节点一致，一起向目标战略位置移动并执行任务。

3.根据权利要求1所述的根据运动的无人机自组网中节点的速度相似度、节点间平均距离、平均链路维持时间、无人机节点剩余能量等因素计算每个无人机节点的权值进行簇的维护，其特征在于基于权值计算公式，计算每个节点权值大小的方法为：

在初始化分簇完成后，随着网络的运行与信息的交换，网络中的节点将会了解自身与邻居节点的距离、邻居节点的速度、方向等信息。在后续的分簇维护过程中，将利用加权分簇的思想，综合考虑节点的速度相似度与邻居节点的平均距离、平均通信链路维持时间、节点的剩余能量等因素对每个节点计算自身的权值W_i：

其中W_V，W_D，W_T，W_P为权值系数，满足W_V+W_D+W_T+W_P＝1，具体大小根据实际情况进行取值。V_i为节点的速度相似度，D_i为与邻居节点平均距离的度，T_i为平均通信链路维持时间，P_i为节点的剩余能量。接下来对公式(4.1)中的节点速度相似度、与邻居节点的平均距离、平均通信链路维持时间进行理论推导。

对于节点i的速度相似度V_i，设节点i与一个邻居节点j的速度标量分别为v_i与v_j，v_i与v_j与水平轴的夹角分别为θ_i与θ_j，则节点i与节点j在水平轴(即x轴)上的速度分量差为：

v_ijx＝v_icosθ_i-v_jcosθ_j (2)

节点i与节点j在纵轴(即y轴)上的速度分量差为：

v_ijy＝v_isinθ_i-v_jsinθ_j (3)

假设节点i有N个邻居节点，则节点i与周围邻居节点在x轴上的平均速度差为：

节点i与周围邻居节点在y轴上的平均速度差为：

则节点i与周围邻居节点的平均速度差为：

根据上述结果，可以分别求出节点i与周围邻居节点在x轴与y轴上平均速度差的方差与标准差。节点i与周围邻居节点在x轴上的平均速度差的标准差为：

节点i与周围邻居节点在y轴上的平均速度差的标准差为：

根据所求的节点i与周围邻居节点的平均速度差的标准差，本章将节点i与周围邻居节点的速度相似度V_i定义为：

其中σ_i表示为节点i与周围邻居节点在x轴与y轴方向速度差的方差的平均值的开方，即：

对于节点i与邻居节点平均距离的度D_i，可以求出节点i与周围邻居节点的平均距离为：

根据节点i与周围邻居节点的平均距离可求得距离的标准差，即：

根据节点i与周围邻居节点的平均距离及距离的标准差，将节点i与邻居节点平均距离的度D_i定义为：

接下来求节点i与周围邻居节点的平均链路通信维持时间。先求单条链路通信维持时间，以节点i与其一邻居节点j为例进行说明，设节点i所在位置为点I，节点j所在位置为J，节点I的通信半径为R，节点J与节点I的相对速度为v_ji，v_ji的大小可由公式(2)与(3)求得，即：

则节点i与节点j的链路生存时间可以表示为：

根据所求的T_ij即可求出节点i与周围邻居节点的平均通信链路生存时间为：

将所求的值代入公式(1)，每个节点即可计算出自身的权值W_i。在后续簇维护过程中，簇首节点会收集簇内成员节点的权值，每一个簇首选举周期如果簇首节点没有失效，则会周期性地按照权值大小更新簇首，以保持簇的稳定性。如果簇首节点失效，则簇内节点会设定退避计时器进行退避，退避完成后首先发送一次入簇请求帧，随后设置定时器等待一段WF_CLUSTER_RTJ时间，如果定时器超时，则表明节点一跳范围内没有簇首，节点开始竞选簇首。

为提高网络覆盖率，在分簇完成后，竞选成为簇首的节点在发送数据帧时会增大发射功率，增大发射功率的值根据网络拓扑需要设定，簇内成员节点依旧采用默认发射功率。在网络运行过程中，MAC层发送控制帧与数据帧采用CSMA/CA协议的基本思想，可以在网络中节点数量较多的情况下尽可能地避免冲突，保证一定的网络性能。网络层路由部分，簇首节点之间采用DSR协议，簇内成员节点发送数据包时，首先广播路由请求包(Route Request，RREQ)，并设置等待定时器，如果目的节点在其一跳范围内，则回送路由应答包(RouteReply，RREP)，而如果目的节点不在RREQ发送节点一跳范围内，则节点会判断自身是否为簇首节点，只有簇首节点才会做进一步的路由请求转发。在该簇首节点转发该路由请求后，其他节点接收到该路由请求，同样只有簇首节点会接收该转发路由请求并作进一步转发，簇内成员节点只有接收到自身所在簇的簇首转发的路由请求才会接收和处理，如果自身为路由请求包中的目的节点，则回送路由应答包RREP，由簇首节点转发回源节点，完成路由发现过程。

4.根据权利要求1所述的一种基于邻域内邻接矩阵的网络中关键节点判断方法，其特征在于利用图论中的邻域内邻矩阵判断节点是否为关键节点的方法为：

设节点v的邻居节点集合N_max(v)＝{v₁，v₂，...，v_n|d(v，v_i)＜r(v)，i＝1，2，...，n}。其中d(v，v_i)表示节点v与邻居节点v_i的距离，r(v)表示节点v的通信半径。设节点v的一个邻域节点为v_u，建立集合X(v_u)＝Φ，遍历N_max(v)中除v与u之外的节点v_i，将a_ui＝1(即节点u与节点i可达)对应的节点v_i加入集合X(v_u)。随后对X(v_u)进行判断，如果X(v_u)为空集，即说明节点u与节点v的其他邻居节点都不连通，即如果节点v失效，会造成网络的分割。随后按照此流程对集合N_max(v)中的其他节点做如上的操作，判断节点v是否为关键节点。采用该思想，即可以求出节点v的邻域内邻接矩阵，若邻域内接矩阵中有一行全为0，则节点v是关键节点。

根据上述图论的思想，网络中的簇首节点即可以根据两跳邻居节点信息结合有向图与邻域内邻接矩阵分布式地判断出自身是否为关键节点。具体步骤如下：

网络中的簇首节点周期性地发送簇首广播帧，簇首广播帧中携带节点位置、速度、邻居节点信息等信息。

一跳范围内的其它节点接收到该簇首广播帧后，首先判断自身是否为簇内成员节点或者其它簇首节点，若满足条件则接收该簇首广播帧。簇内成员节点更新簇首有效时间与位置、速度等信息，随后丢弃该簇首广播帧以释放内存空间。接收该簇首广播帧的其它簇首节点根据该簇首广播帧中节点携带的邻居节点表与关键节点判断位信息更新自身的邻接矩阵，并将该节点信息存入自身维护的邻居节点数据表中，当节点更新邻接矩阵与信息完毕后，同样会丢弃掉该簇首广播帧。

随后，在该接收簇首节点自身发送簇首广播帧时，会启动关键节点判断流程，如果邻域内邻接矩阵中存在有一行全为0的情况，即说明该节点为关键节点，将簇首广播帧头部的关键节点标志位信息置为1，随后向周围邻居簇首节点进行广播。全网簇首节点都按照这样的流程不断更新周围簇首与自身的信息，从而可以实现分布式的全网关键簇首节点判断。

以上介绍了簇首节点之间进行失效簇首节点的方法。若失效的簇首节点范围内还有其他在该簇首失效前维护的簇内成员节点，则簇内成员也会通过簇首失效检测机制检测到簇首失效。下面将对簇内成员的簇首失效检测机制进行介绍。

簇内成员节点在每一次接收到簇首广播帧后会根据接收时间设定一个定时器，定时器超时的时间为当前时间加上簇首失效检测时间间隔，当定时器超时即可判定簇首节点已经失效，随后根据失效簇首节点信息与所计算的权值进行新簇首的选举与调整。如果在该定时器超时前收到了新的簇首广播帧，则表明簇首节点依然有效，节点更新本地维护的簇首节点信息，并取消之前的定时器，设置新的定时器，继续等待下一个簇首广播帧。

5.根据权利要求1所述的提出了两阶段拓扑修复策略，其特征在于为每个时帧中的不同阶段设计的内容流程为：

如果失效簇首范围内存在原簇内成员，则簇内成员节点会通过簇首失效检测机制发现该簇首失效，随后通过之前所述加权分簇的方法产生一个新的簇首，完成拓扑修复，竞选成功的簇首会广播簇首修复告知帧，告知周围的簇首节点原簇首失效的情况已被修复，并且根据原簇首节点的位置与速度信息进行位置的调整。周围簇首节点收到该簇首修复告知帧后即知道原失效节点被修复，不会启动第二阶段拓扑修复过程。

当簇内节点检测到簇首节点失效后会启动簇内修复流程，进入二次竞选簇首阶段，选举出一个权值最小的节点作为新的簇首，完成拓扑的修复。节点首先设置定时器，随机退避一段时间，如果退避定时器超时，即顺利完成退避，则节点将自身计算的权值添加到帧头部，广播簇首竞选帧，发送完毕后，节点将自身节点ID号设为当前簇首，随后设置定时器，等待一段WF_CLUSTER时间，如果在定时器超时前未收到其它簇首竞选帧，则节点判定自身竞选簇首成功，随后计算原簇首未失效情况下当前时刻的位置，并进行调整。如果节点在退避完成前收到了其它节点的簇首竞选帧，则会首先查看该簇首竞选帧中携带的权值数值，并与自身进行比较，如果该簇首竞选帧中权值比自身要大，则节点将该簇首竞选帧中携带的源节点地址设为自身的簇首节点地址，重新设置定时器等待WF_CLUSTER时间；否则节点将判断权值是否相等，如果权值相等，则进一步比较节点ID号，如果收到的簇首竞选帧中的ID号更小，则将簇首竞选帧中的源节点地址设置为自身的簇首节点地址，重新设置定时器等待WF_CLUSTER时间，相反，如果自身的ID号更小，则继续进行退避。如果权值不相等，则显然只有自身权值更大这一种情况，节点会丢弃该簇首竞选帧，继续进行退避。如果经历了这一阶段后节点将簇首更改为收到的簇首竞选帧中的源节点并进入等待WF_CLUSTER状态，则如果定时器超时则重新分簇成功，簇首节点为新收到的簇首竞选帧中的源节点；如果在定时器超时前又收到了新的簇首竞选帧，则重新对簇首竞选帧中的权值进行判断，如果簇首竞选帧中的权值大于当前节点维护的簇首的权值，则节点将该簇首竞选帧的源节点更新为簇首，重新设置定时器等待WF_CLUSTER时间；否则按照上述判断流程，如果权值相等但自身节点号更小或者权值不相等(即自身权值更大)，则丢弃收到的簇首竞选帧，继续等待之前的定时器超时完成分簇。否则就将该簇首竞选帧的源节点更新为簇首，重新设置定时器等待WF_CLUSTER时间。在定时器超时前未收到新的簇首竞选帧则判定自身入簇成功。

当新的簇首节点产生后，会广播簇首修复告知帧，告知周围的簇首节点原失效簇已被修复，同时簇首节点会根据失效簇首节点的位置与速度信息进行位置的调整。

周围簇首节点在原簇首节点失效后会设置定时器等待一段WF_FirstStage_Reconstruction时间，如果该定时器超时前收到该失效簇首的簇首修复告知帧，则更新新的簇首信息后丢弃该帧；如果该定时器超时且并没有收到簇首修复告知帧，即表明关键簇整体失效，启动第二阶段拓扑修复流程，即簇间移动控制拓扑修复方案。

当簇首节点检测到邻居簇首节点失效，并且在WF_FirstStage_Reconstruction时间后仍未收到新簇首修复告知帧，则首先统计失效簇首节点的邻居簇首节点数量，随后根据失效簇首邻居节点数量与自身簇内成员数量计算协作修复概率，根据该概率发送拓扑修复请求帧，请求簇内成员节点进行拓扑修复。簇内成员节点收到该拓扑修复请求帧后，会随机退避一段时间并向簇首节点回复修复应答帧，在修复应答帧中包括了成员节点当前的位置信息。簇首节点根据簇内成员节点的位置信息与剩余能量综合选择一个节点进行拓扑修复，并向该节点回复修复确认帧，该节点收到后即知道自身将承担修复任务，随机开始拓扑修复进程。

簇内节点接收到应答确认帧，得知自身被选择参与失效关键簇首节点修复，即进入簇内成员节点移动修复过程。对于修复节点来说，首先需要确定失效簇首节点的位置(即如果关键簇首节点未失效，此时应运动到的位置)，修复节点可以从修复请求帧中获取失效关键簇首节点失效前的位置，用历史位置来推测当前位置，也可以通过关键簇首节点失效前的速度、方向等信息对当前时刻的预期位置进行计算，为了获得较为准确的失效关键簇首位置，修复节点采用卡尔曼滤波法对这两种估计结果进行加权计算随后在虚拟力作用下向确定位置进行移动。

当修复节点获得较为准确的失效关键簇首位置之后，便会在虚拟力的作用下向着该位置移动，假设该目的地点位置坐标为D(x_d，y_d，z_d)，修复节点i坐标为I(x_i，y_i，z_i)，则修复节点受到的吸引力可表示为：

其中M_a为常数，表示吸引力系数，

的方向与节点i到目的位置的方向相同，且大小与节点i和目的位置之间的距离正相关。在修复节点向目的位置移动的过程中，还可能因为与周围节点距离过近而受到排斥力，设某一时刻节点i的邻居节点的集合为N，假设集合N中一个节点为j，坐标为J(x_j，y_j，z_j)，则节点i受到节点j的排斥力为：

其中M_r为常数，表示排斥力系数，

的方向与节点i到邻居节点的方向相反，且大小与节点i到邻居节点之间的距离正相关。结合修复无人机受到的来自障碍物与边界的排斥力，无人机i受到的合力可以表示为：

其中O为无人机节点i附近的障碍物与边界的集合。在按照上述流程计算出合力后，进一步判断该计算值是否为0，若为0则对节点施加随机扰动，让节点跳出当前状态。若该计算出的值不为0，则继续比较当前矢量方向与之前的方向是否相反，若相反，则修复节点会出现往复移动的情况，同样需要施加一次随机扰动力。若不相反，则节点按照合力的方向进行移动，直至移动到待修复范围内，完成拓扑修复任务。综上，修复节点受到的合力可表示为：

修复节点的运动方程可表示为：

其中P_i(t)表示t时刻节点坐标，P_i(t+Δt)表示t+Δt时刻节点坐标，F_i/|F_i|表示当前时刻前进方向。

当修复节点到达修复位置后，为了防止有其它节点已经提前到达修复位置并且建簇成功，会首先广播一次入簇请求帧，随后设置定时器等待一段WF_Cluster_RTJ时间，定时器超时后即表明修复节点邻居范围内没有已存在的簇，随后进入竞选簇首流程，广播一次簇首竞选帧并设置定时器等待一段WF_Cluster时间，定时器超时后即表明建簇成功，发送簇首广播帧告知周围邻居簇首节点，拓扑修复完成。如果修复节点在定时器超时前收到了其它簇首回复的入簇应答帧，则表明已经有其它簇派出的修复节点到达了失效簇首位置，并且成功建簇，修复节点加入该簇，成为簇内成员节点。

随后，簇内成员节点之间采取基于虚拟力导向的控制方法自主地实现局部拓扑的优化，节点受到来自簇内其它成员节点之间的吸引力与排斥力、与障碍物和边界之间的避障力、与新簇首之间的吸引力。修复节点i所受的合力可以表示为：

其中N表示节点i的邻居节点集，O表示障碍物指示位置与边界的集合，

表示节点受到簇首节点的吸引力。节点可由该式计算出自身受力大小与方向，从而可以实现动态的分布式的拓扑优化。簇内成员节点在该虚拟力的作用下进入平衡状态后，根据簇首广播帧中的速度信息调整自身的速度，与簇首节点保持一致，继续向战略目标位置移动。

Images