[go: up one dir, main page]

CN103336527A - Device node formation method, device and system - Google Patents

Device node formation method, device and system Download PDF

Info

Publication number
CN103336527A
CN103336527A CN2013102257773A CN201310225777A CN103336527A CN 103336527 A CN103336527 A CN 103336527A CN 2013102257773 A CN2013102257773 A CN 2013102257773A CN 201310225777 A CN201310225777 A CN 201310225777A CN 103336527 A CN103336527 A CN 103336527A
Authority
CN
China
Prior art keywords
node
distance
azimuth
formation
individual
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.)
Pending
Application number
CN2013102257773A
Other languages
Chinese (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.)
Hanshan Normal University
Original Assignee
Hunan University of Science and Engineering
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 Hunan University of Science and Engineering filed Critical Hunan University of Science and Engineering
Priority to CN2013102257773A priority Critical patent/CN103336527A/en
Publication of CN103336527A publication Critical patent/CN103336527A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Manipulator (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明提供了一种节点编队方法、装置及系统,该节点编队方法包括:步骤S102,检测第一节点与第二节点之间的第一距离以及第一方位角,以及检测第一节点与第三节点之间的第二距离以及第二方位角,并计算出第二节点与第三节点之间的第三距离;步骤S104,判断第一距离、第二距离以及第三距离是否均小于或等于预设的第一阈值长度,如果是,则执行步骤S108,如果否,则执行步骤S106;步骤S106,计算出平均位置点,并驱动第一节点向平均位置点移动,并以预定时间间隔返回执行步骤S102;步骤S108,计算出目标位置点,并驱动第一节点向目标位置点移动,并以预定时间间隔执行步骤S110;步骤S110,判断第一节点是否到达目标位置点,如果否,则返回执行步骤S102,如果是,则第一节点停止编队过程。

The present invention provides a node formation method, device and system, the node formation method includes: step S102, detecting the first distance and the first azimuth between the first node and the second node, and detecting the first node and the second node The second distance and the second azimuth between the three nodes, and calculate the third distance between the second node and the third node; step S104, judge whether the first distance, the second distance and the third distance are all less than or Equal to the preset first threshold length, if yes, execute step S108, if not, execute step S106; step S106, calculate the average position point, and drive the first node to move to the average position point, and at a predetermined time interval Return to step S102; step S108, calculate the target location point, and drive the first node to move to the target location point, and execute step S110 at predetermined time intervals; step S110, judge whether the first node reaches the target location point, if not, Then return to step S102, if yes, the first node stops the formation process.

Description

节点编队方法、装置及系统Node formation method, device and system

技术领域 technical field

本发明涉及自动控制领域,具体而言,涉及一种节点编队方法、装置及系统。  The present invention relates to the field of automatic control, in particular to a node formation method, device and system. the

背景技术 Background technique

目前用于多智能体编队行为的分布式控制方法主要有势函数法、人工物理法、基于行为法,网络图论法等。这些方法从某种程度上都能够实现三近邻智能体的聚集行为。  At present, the distributed control methods for multi-agent formation behavior mainly include potential function method, artificial physics method, behavior-based method, network graph theory method, etc. These methods can achieve the aggregation behavior of three-nearest neighbor agents to some extent. the

但是,当前的这些方法出发点并不是为了形成多成员的规则队形结构,势函数法,人工物理法存在着多控制参数,编队行为对参数敏感;而基于行为法存在行为控制精度的偏差和较大的行为响应延迟;网络图论法虽然能形成指定的队形,然而各个体动作存在着先后次序而不能并发执行,需花费时间较多。  However, the starting point of these current methods is not to form a multi-member regular formation structure. The potential function method and the artificial physics method have multiple control parameters, and the formation behavior is sensitive to parameters; while the behavior-based method has deviations in behavior control accuracy and relatively low Large behavior response delay; although the network graph theory method can form a specified formation, however, each individual action has a sequence and cannot be executed concurrently, which takes a lot of time. the

针对相关技术中多智能体编队行为的分布式控制所存在的上述问题,目前尚未提出有效的解决方案。  For the above-mentioned problems in the distributed control of multi-agent formation behavior in related technologies, no effective solution has been proposed yet. the

发明内容 Contents of the invention

本发明提供了一种节点移动编队方法、装置及系统,以至少解决上述现有技术中多智能体编队行为的分布式控制所存在的问题。  The present invention provides a node mobile formation method, device and system to at least solve the problems existing in the distributed control of multi-agent formation behavior in the prior art. the

根据本发明的一个方面,提供了一种节点编队方法,包括:步骤S102,检测第一节点与第二节点之间的第一距离以及第一方位角,以及检测第一节点与第三节点之间的第二距离以及第二方位角,并根据所述第一距离、所述第一方位角、所述第二距离和所述第二方位角计算出所述第二节点与所述第三节点之间的第三距离;步骤S104,判断第一距离、第二距离及第三距离是否均小于或等于预设的第一阈值长度,如果是,则执行步骤S108,如果否,则执行步骤S106;步骤S106,根据第一距离、第一方位角、第二距离和第二方位角计算出第一节点、第二节点和第三节点之间的平均位置点,并驱动第一节点向平均位置点移动,并以预定时间间隔返回执行步骤S102;步骤S108,根据当前第一距离、当前第一方位角、当前第二距离和当前第二方位角计算出目标位置点,其中,目标位置点为到第二节点以及到第三节点的距离均为第一阈值长度的一半并且到第一节点的距离最小的位置点,并驱动第一节点向目标位置点移动,并以预定时间间隔执行步骤S110;步骤S110,判断所述第一节点是否到达所述目标位置点,如果否,则返回执行所述步骤S102,如果是,则第一节点停止编队过程。  According to one aspect of the present invention, a node formation method is provided, including: step S102, detecting the first distance and the first azimuth between the first node and the second node, and detecting the distance between the first node and the third node The second distance between and the second azimuth, and calculate the second node and the third azimuth according to the first distance, the first azimuth, the second distance and the second azimuth The third distance between nodes; step S104, judge whether the first distance, the second distance and the third distance are all less than or equal to the preset first threshold length, if yes, then execute step S108, if not, then execute step S106; step S106, according to the first distance, the first azimuth, the second distance and the second azimuth, calculate the average position point between the first node, the second node and the third node, and drive the first node to the average The position point moves, and returns to step S102 at a predetermined time interval; step S108, calculates the target position point according to the current first distance, the current first azimuth angle, the current second distance and the current second azimuth angle, wherein the target position point It is the position point whose distance to the second node and the third node are both half of the first threshold length and the distance to the first node is the smallest, and drives the first node to move to the target position point, and executes steps at predetermined time intervals S110; Step S110, judging whether the first node has reached the target point, if not, return to step S102, if yes, the first node stops the formation process. the

优选地,在步骤S102之后还包括:当根据当前第一方位角和第二方位角判断出第一节点、第二节点和第三节点处于同一直线上,并且第一节点位于第二节点和第三节点的同一侧时,停止驱动第一节点。  Preferably, after step S102, it also includes: when it is judged according to the current first azimuth angle and the second azimuth angle that the first node, the second node and the third node are on the same straight line, and the first node is located between the second node and the second node When the same side of the three nodes, stop driving the first node. the

优选地,在步骤S102之后还包括:当根据当前第一方位角和第二方位角判断出第一节点、第二节点和第三节点处于同一直线上,并且第一节点位于第二节点和第三节点的中间时,驱动第一节点沿该直线的垂直线的任意一侧移动。  Preferably, after step S102, it also includes: when it is judged according to the current first azimuth angle and the second azimuth angle that the first node, the second node and the third node are on the same straight line, and the first node is located between the second node and the second node In the middle of the three nodes, drive the first node to move along any side of the vertical line of the straight line. the

优选地,平均位置点为第一节点、第二节点和第三节点所组成的三角形的重心或内心。  Preferably, the average position point is the center of gravity or incenter of the triangle formed by the first node, the second node and the third node. the

优选地,在上述过程中,同时,与同样的方式分别独立地驱动第二节点和第三节点的移动。  Preferably, during the above process, at the same time, the movement of the second node and the third node are independently driven in the same manner. the

优选地,第一节点、第二节点和第三节点均为智能机器人。  Preferably, the first node, the second node and the third node are all intelligent robots. the

根据本发明的另一方面,提供了一种节点编队装置,位于第一节点中,该节点编队装置包括:检测单元,用于检测第一节点与第二节点之间的第一距离以及第一方位角,以及检测第一节点与第三节点之间的第二距离以及第二方位角;计算单元,用于判断第一距离、第二距离及第三距离是否均小于或等于预设的第一阈值长度,如果否,则根据第一距离、第一方位角、第二距离和第二方位角计算出第一节点、第二节点和第三节点之间的平均位置点;如果是,则根据第一距离、第一方位角、第二距离和第二方位角计算出目标位置点,其中,目标位置点为到第二节点以及到第三节点的距离均为第一阈值长度的一半并且到第一节点的距离最小的位置点;驱动单元,用于根据计算单元的计算结果驱动第一节点向平均位置点或目标位置点移动。  According to another aspect of the present invention, a node formation device is provided, located in a first node, the node formation device includes: a detection unit, configured to detect a first distance between a first node and a second node and a first Azimuth, and detecting the second distance and the second azimuth between the first node and the third node; the calculation unit is used to judge whether the first distance, the second distance and the third distance are all less than or equal to the preset first A threshold length, if not, then calculate the average position point between the first node, the second node and the third node according to the first distance, the first azimuth, the second distance and the second azimuth; if yes, then The target position point is calculated according to the first distance, the first azimuth, the second distance and the second azimuth, wherein the target position point is half of the first threshold length to the second node and to the third node, and The position point with the smallest distance to the first node; the driving unit, used to drive the first node to move to the average position point or the target position point according to the calculation result of the calculation unit. the

优选地,计算单元还用于根据当前第一方位角和第二方位角判断第一节点、第二节点和第三节点是否处于同一直线上以及它们之间的相邻关系。  Preferably, the calculation unit is also used to judge whether the first node, the second node and the third node are on the same straight line and the adjacent relationship among them according to the current first azimuth angle and the second azimuth angle. the

优选地,驱动单元还用于在计算单元判断出第一节点、第二节点和第三节点处于同一直线上并且第一节点位于第二节点和第三节点的同一侧的情况下,控制第一节点停止移动。  Preferably, the drive unit is also used to control the first node, the second node and the third node to be on the same straight line and the first node to be on the same side of the second node and the third node when the calculation unit judges that The node stops moving. the

优选地,驱动单元还用于在计算单元判断出第一节点、第二节点和第三节点处于同一直线上,并且第一节点位于第二节点和第三节点之间的情况下,驱动第一节点沿该直线的垂直线的任意一侧移动。  Preferably, the driving unit is also used to drive the first Nodes move along either side of a vertical line to that line. the

优选地,平均位置点为第一节点、第二节点和第三节点所组成的三角形的重心或内心。  Preferably, the average position point is the center of gravity or incenter of the triangle formed by the first node, the second node and the third node. the

优选地,检测单元包括:红外传感器,用于检测第一距离和第二距离;数字指南针,用于检测第一方位角和第二方位角;  Preferably, the detection unit includes: an infrared sensor for detecting the first distance and the second distance; a digital compass for detecting the first azimuth and the second azimuth;

根据本发明的又一方面,提供了一种节点编队系统,包括:前文的第一节点、第二节点和第三节点,其中,第二节点和第三节点分别具有与第一节点相同的检测单元、计算单元和驱动单元。  According to yet another aspect of the present invention, a node formation system is provided, including: the first node, the second node, and the third node mentioned above, wherein the second node and the third node respectively have the same detection unit, calculation unit and drive unit. the

优选地,第一节点、第二节点和第三节点均为智能机器人。  Preferably, the first node, the second node and the third node are all intelligent robots. the

在本发明中,通过在节点之间的初始距离过大时,以驱动节点向平均位置点靠近,直到节点满足距离要求时,才直接转入位置调整,解决了相关技术中多智能体编队行为的分布式控制所存在的控制精度不高以及较大的响应延迟的问题,进而达到了提高多智能体编队的控制精度和加快响应的效果。  In the present invention, when the initial distance between the nodes is too large, the nodes are driven closer to the average position point, until the nodes meet the distance requirements, the position adjustment is directly transferred, and the multi-agent formation behavior in the related art is solved. The problems of low control precision and large response delay in the distributed control of the distributed control system have achieved the effect of improving the control precision of the multi-agent formation and speeding up the response. the

附图说明 Description of drawings

此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:  The accompanying drawings described here are used to provide a further understanding of the present invention and constitute a part of the application. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations to the present invention. In the attached picture:

图1是根据本发明实施例的节点编队方法流程图;  Fig. 1 is a node formation method flowchart according to an embodiment of the present invention;

图2是根据本发明实施例的节点编队装置的模块结构示意图;  Fig. 2 is a schematic diagram of a module structure of a node formation device according to an embodiment of the present invention;

图3是根据本发明实施例的节点编队系统分布示意图;  Fig. 3 is a schematic diagram of the distribution of a node formation system according to an embodiment of the present invention;

图4是根据本发明实施例一的个体状态转移模型图;  Fig. 4 is an individual state transition model diagram according to Embodiment 1 of the present invention;

图5是根据本发明实施例的个体Ai获取邻近个体局部位置信息的过程示意图;  Fig. 5 is a schematic diagram of the process of obtaining local position information of adjacent individuals by individual A i according to an embodiment of the present invention;

图6a至6d是根据本发明实施例的计算pig的示意图;  6a to 6d are schematic diagrams of calculating p ig according to an embodiment of the present invention;

图7是根据本发明实施例一的相对距离及个体Ai理想目的地位置确定示意图;  Fig. 7 is a schematic diagram of determining the relative distance and the ideal destination position of an individual A i according to Embodiment 1 of the present invention;

图8是根据本发明实施例一的个体Ai分别执行单纯CA及AA算法时的状态流程图;  Fig. 8 is a state flow diagram of individual A i respectively executing simple CA and AA algorithms according to Embodiment 1 of the present invention;

图9是根据本发明实施例一的个体Ai执行TFS交互控制方法时的状态流程图;  Fig. 9 is a state flow chart when an individual A i executes a TFS interactive control method according to Embodiment 1 of the present invention;

图10是根据本发明实施例一的三近邻个体位置共线场景示意图;  Fig. 10 is a schematic diagram of a collinear scene of three nearest neighbor individual positions according to Embodiment 1 of the present invention;

图11是根据本发明实施例一的修正子算法流程图;  Fig. 11 is a flow chart of the correction sub-algorithm according to Embodiment 1 of the present invention;

图12是根据本发明实施例一的转换将发生时的临界位置示意图;  Fig. 12 is a schematic diagram of a critical position when conversion will occur according to Embodiment 1 of the present invention;

图13是根据本发明实施例一的三近邻个体CF-2分布时的E编队过程示意图;以及  Fig. 13 is a schematic diagram of the E formation process when three neighbor individuals CF-2 are distributed according to Embodiment 1 of the present invention; and

图14是根据本发明实施例一的典型初始分布下三近邻个体的平均边长误差与时间的变化关系示意图。  Fig. 14 is a schematic diagram of the relationship between the average side length error and time of three neighboring individuals under a typical initial distribution according to Embodiment 1 of the present invention. the

具体实施方式 Detailed ways

下文中将参考附图并结合实施例来详细说明本发明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。  Hereinafter, the present invention will be described in detail with reference to the drawings and examples. It should be noted that, in the case of no conflict, the embodiments in the present application and the features in the embodiments can be combined with each other. the

图1是根据本发明实施例的节点编队方法流程图。如图1所示,该节点编队方法包括以下步骤:  Fig. 1 is a flowchart of a node formation method according to an embodiment of the present invention. As shown in Figure 1, the node formation method includes the following steps:

步骤S102,检测第一节点与第二节点之间的第一距离以及第一方位角,以及检测第一节点与第三节点之间的第二距离以及第二方位角,并根据所述第一距离、所述第一方位角、所述第二距离和所述第二方位角计算出所述第二节点与所述第三节点之间的第三距离;  Step S102, detecting the first distance and the first azimuth between the first node and the second node, and detecting the second distance and the second azimuth between the first node and the third node, and according to the first The distance, the first azimuth, the second distance and the second azimuth calculate a third distance between the second node and the third node;

步骤S104,判断第一距离、第二距离及第三距离是否均小于或等于预设的第一阈值长度,如果是,则执行步骤S108,如果否,则执行S106;  Step S104, judging whether the first distance, the second distance and the third distance are all less than or equal to the preset first threshold length, if yes, execute step S108, if not, execute S106;

步骤S106,根据第一距离、第一方位角、第二距离和第二方位角计算出第一节点、第二节点和第三节点之间的平均位置点,并驱动第一节点向平均位置点移动,并以预定时间间隔返回执行步骤S102;  Step S106, calculate the average position point between the first node, the second node and the third node according to the first distance, the first azimuth angle, the second distance and the second azimuth angle, and drive the first node to the average position point Move, and return to execute step S102 at predetermined time intervals;

步骤S108,根据当前第一距离、当前第一方位角、当前第二距离和当前第二方位角计算出目标位置点,其中,目标位置点为到第二节点以及到第三节点的距离均为第一阈值长度的一半并且到第一节点的距离最小的位置点,并驱动第一节点向目标位置点移动,并以预定时间间隔执行步骤S110;  Step S108, calculate the target position point according to the current first distance, the current first azimuth angle, the current second distance and the current second azimuth angle, wherein the target position point is the distance to the second node and the distance to the third node is Half of the first threshold length and the minimum distance to the first node, and drive the first node to move to the target position, and perform step S110 at predetermined time intervals;

步骤S110,判断第一节点是否到达目标位置点,如果否,则返回执行步骤S102,如果是,则第一节点停止编队过程。  Step S110, judging whether the first node has reached the target point, if not, return to step S102, if yes, the first node stops the formation process. the

在本实施例中,通过在节点之间的初始距离过大时,以驱动节点向平均位置点靠近,直到节点满足距离要求时,才驱动节点向目标位置点移动,解决了相关技术中多智能体编队行为的分布式控制所存在的控制精度不高以及较大的响应延迟的问题,进而达到了提高多智能体编队的控制精度和加快响应的效果。  In this embodiment, when the initial distance between the nodes is too large, the nodes are driven to move closer to the average position point, and the nodes are not driven to move to the target position point until the nodes meet the distance requirements, which solves the problem of multi-intelligence in related technologies. The distributed control of agent formation behavior has the problems of low control precision and large response delay, and then achieves the effect of improving the control accuracy and speeding up the response of multi-agent formation. the

在上述方法中,除第一节点外,第二节点和第三节点也同时异步执行与第一节点同样的节点编队策略,即,第一、第二和第三节点的位置调整是动态进行的。  In the above method, in addition to the first node, the second node and the third node also execute the same node formation strategy as the first node asynchronously at the same time, that is, the position adjustment of the first, second and third nodes is carried out dynamically . the

其中,在步骤S102之后还包括:当根据当前第一方位角和第二方位角判断出第一节点、第二节点和第三节点处于同一直线上,并且第一节点位于第二节点和第三节点的同一侧时,第一节点暂时停留在原位置;如果第一节点、第二节点和第三节点处于同一直线上,并且第一节点位于第二节点和第三节点的中间时,则驱动第一节点沿该直线的垂直线的任意一侧移动,这样的话,各个节点则可继续执行上述移动策略。  Wherein, after step S102, it also includes: when judging that the first node, the second node and the third node are on the same straight line according to the current first azimuth angle and the second azimuth angle, and the first node is located between the second node and the third node When the nodes are on the same side, the first node temporarily stays at the original position; if the first node, the second node and the third node are on the same straight line, and the first node is in the middle of the second node and the third node, the drive A node moves along either side of the vertical line of the line, so that each node can continue to implement the above moving strategy. the

其中,平均位置点可以选择为第一节点、第二节点和第三节点所组成的三角形的重心或内心。  Wherein, the average position point can be selected as the center of gravity or incenter of the triangle formed by the first node, the second node and the third node. the

其中,第一节点、第二节点和第三节点为需要进行编队移动的个体,例如,智能机器人。  Wherein, the first node, the second node and the third node are individuals that need to move in formation, for example, intelligent robots. the

图2是根据本发明实施例的节点编队装置模块结构示意图。如图2所示,该节点编队装置包括以下功能模块:  Fig. 2 is a schematic structural diagram of a node formation device module according to an embodiment of the present invention. As shown in Figure 2, the node formation device includes the following functional modules:

检测单元120,用于检测第一节点100与第二节点200之间的第一距离以及第一方位角,以及检测第一节点100与第三节点300之间的第二距离以及第二方位角;  A detection unit 120, configured to detect a first distance and a first azimuth between the first node 100 and the second node 200, and detect a second distance and a second azimuth between the first node 100 and the third node 300 ;

计算单元140,用于判断第一距离、第二距离及第三距离是否均小于或等于预设的第一阈值长度,如果否,则根据第一距离、第一方位角、第二距离和第二方位角计算出第一节点100、第二节点200和第三节点300之间的平均位置点;如果是,则根据第一距离、第一方位角、第二距离和第二方位角计算出目标位置点,其中,目标位置点为到第二节点200以及到第三节点300的距离均为第一阈值长度的一半并且到第一节点100的距离最小的位置点;  Calculation unit 140, for judging whether the first distance, the second distance and the third distance are all less than or equal to the preset first threshold length, if not, according to the first distance, the first azimuth, the second distance and the first distance The second azimuth calculates the average position point between the first node 100, the second node 200 and the third node 300; if yes, then calculates according to the first distance, the first azimuth, the second distance and the second azimuth Target position point, wherein, the target position point is the position point that the distance to the second node 200 and the distance to the third node 300 is half of the first threshold length and the distance to the first node 100 is the smallest;

驱动单元160,用于根据计算单元140的计算结果驱动第一节点100向平均位置点或目标位置点移动。  The driving unit 160 is configured to drive the first node 100 to move to the average position point or the target position point according to the calculation result of the calculation unit 140 . the

其中,检测单元120可包括:红外传感器,用于检测第一距离和第二距离;数字指南针,用于检测第一方位角和第二方位角;  Wherein, the detection unit 120 may include: an infrared sensor, for detecting the first distance and the second distance; a digital compass, for detecting the first azimuth and the second azimuth;

其中,计算单元140还用于根据当前第一方位角和第二方位角判断第一节点100、第二节点200和第三节点300是否处于同一直线上以及它们之间的相邻关系。  Wherein, the computing unit 140 is also used to judge whether the first node 100 , the second node 200 and the third node 300 are on the same straight line and the adjacent relationship among them according to the current first azimuth angle and the second azimuth angle. the

其中,驱动单元160还用于在计算单元140判断出第一节点100、第二节点200和第三节点300处于同一直线上并且第一节点100位于第二节点200和第三节点300的同一侧的情况下,控制第一节点100停止移动,以及第一节点100、第二节点200和第三节点300处于同一直线上并且第一节点100位于第二节点200和第三节点300之间的情况下,驱动第一节点100沿该直线的垂直线的任意一侧移动。  Wherein, the driving unit 160 is also used to judge in the calculation unit 140 that the first node 100, the second node 200 and the third node 300 are on the same straight line and the first node 100 is located on the same side of the second node 200 and the third node 300 In the case of controlling the first node 100 to stop moving, and the first node 100, the second node 200 and the third node 300 are on the same straight line and the first node 100 is located between the second node 200 and the third node 300 Next, drive the first node 100 to move along any side of the vertical line of the straight line. the

图3是根据本发明实施例的节点编队系统分布示意图。如图3所示,该节点编队系统包括前文所述的第一节点100、第二节点200和第三节点300,其中,第二节点200和第三节点300分别具有与第一节点100相同的检测单元120、计算单元140和驱动单元160。即,该系统中,第一、第二和第三节点的位置调整是动态进行的。在实际的应用中,第一、第二和第三节点可以是智能机器人,当然根据实际应用的不同,这些节点也可以是其它具有移动和控制功能的个体,例如,无人机、传感器等。通过本实施的所提供的节点编队系统可以实现三近邻节点的编队。  Fig. 3 is a schematic diagram of distribution of a node formation system according to an embodiment of the present invention. As shown in FIG. 3 , the node formation system includes the aforementioned first node 100, second node 200 and third node 300, wherein the second node 200 and the third node 300 have the same A detection unit 120 , a calculation unit 140 and a drive unit 160 . That is, in this system, the position adjustment of the first, second and third nodes is performed dynamically. In practical applications, the first, second and third nodes can be intelligent robots, of course, according to different practical applications, these nodes can also be other individuals with movement and control functions, such as drones, sensors and so on. The formation of three adjacent nodes can be realized through the node formation system provided by this implementation. the

实施例一  Embodiment one

本实施例详细描述了一个根据本发明实现一个三近邻个体机器人E结构编队的实际应用。  This embodiment describes in detail a practical application of realizing an E-structure formation of three adjacent individual robots according to the present invention. the

1.本实例的个体状态模型:  1. The individual state model of this instance:

在本实施例中,实际个体机器人的结构由四个硬件模块组成:(1)红外(近距)传感器,用于观测邻近个体相对自己的距离;(2)数字指南针,用于观测邻近个体相对自己的方位角;(3)计算单元,利用前两个单元获取到的局部数据,按预先设计好的控制算法计算出下一时刻的目的地位置;(4)驱动器,将控制算法得到的结果转换成相应的可执行运动状态,调整个体运动方向,驱使自己向目的地移动。前两个模块结合在一起负责完成邻近个体的局部位置观测。  In this embodiment, the structure of the actual individual robot is composed of four hardware modules: (1) infrared (proximity) sensor, used to observe the distance of adjacent individuals relative to itself; (2) digital compass, used to observe the relative distance of adjacent individuals own azimuth; (3) the calculation unit, using the local data obtained by the first two units, calculates the destination position at the next moment according to the pre-designed control algorithm; (4) the driver, the result obtained by the control algorithm Convert to the corresponding executable motion state, adjust the individual motion direction, and drive yourself to move to the destination. The first two modules are combined to complete the local position observation of neighboring individuals. the

为更好的将分析集中在控制协作机制本身,结合个体硬件及实际任务执行时的动作特点给出如下个体状态抽象。  In order to better focus the analysis on the control cooperation mechanism itself, the following individual state abstraction is given in combination with the individual hardware and the action characteristics of the actual task execution. the

假设:个体机器人仅具有局部范围(观测半径rs)的环境信息获取能力,拥有观测、计算及移动三种可执行状态,运行期间依次在这三种状态之间循环(如图4所示),状态执行及各状态之间的转移没有时间延迟和误差干扰。  Assumption: the individual robot only has the ability to obtain environmental information in a local range (observation radius r s ), has three executable states of observation, calculation and movement, and cycles between these three states during operation (as shown in Figure 4) , state execution and transition between states without time delay and error interference.

实际上,个体机器人的每个可执行状态都需要时间,状态之间的转换也存在时间上的延迟。例如,在观测阶段,个体需要做旋转扫描利用近距传感器来获取各邻居相对自己的距离,同时,通过数字指南针来获取对应邻居的局部方位角;在计算阶段,个体根据观测阶段所得到的距离和局部方位角信息来确定邻居的局部位置,进而通过执行预先设计的交互控制方法求得下一时刻的目的地位置,这些都需要花费时间。任意个体Ai利用直接观测到的邻居径向 距离和方位角来确定其相对自己的局部位置的过程如下所述。图5表示了个体Ai利用红外传感器从局部坐标系的x正半轴开始逆时针方向旋转扫描一周从而获得两个邻近个体Ai1、Ai2的局部位置的情形。个体Ai1相对个体Ai的距离为di1,局部方位角度为θi1,个体Ai2相对个体Ai的距离为di2,局部方位角度为θi2,而个体Ai3分布在个体Ai的可观测范围之外。个体Ai会按公式(1)计算出近邻个体Ai1相对自己的局部位置pi1=(pi1(x),pi1(y))。类似的,可计算出Ai2的局部位置。而对于Ai3,由于Ai1无法感知到它的存在,故无法获得它的位置信息。  In fact, each executable state of an individual robot takes time, and transitions between states also have time delays. For example, in the observation stage, the individual needs to do a rotation scan to obtain the distance of each neighbor relative to itself by using the proximity sensor, and at the same time, obtain the local azimuth of the corresponding neighbor through the digital compass; in the calculation stage, the individual uses the distance obtained in the observation stage It will take time to determine the neighbor’s local position based on the local azimuth information, and then execute the pre-designed interactive control method to obtain the destination position at the next moment. The process of any individual A i using the directly observed radial distance and azimuth of its neighbors to determine its local position relative to itself is as follows. Fig. 5 shows the situation where the individual A i uses the infrared sensor to scan counterclockwise from the positive semi-axis of the local coordinate system for one circle to obtain the local positions of two adjacent individuals A i1 and A i2 . The distance between individual A i1 and individual Ai is d i1 , and the local azimuth angle is θ i1 . The distance between individual A i2 and individual A i is d i2 , and the local azimuth angle is θ i2 . outside the observation range. The individual A i will calculate the local position p i1 of the neighbor individual A i1 relative to itself according to formula (1) =(p i1 (x), p i1 (y)). Similarly, the local position of A i2 can be calculated. As for A i3 , since A i1 cannot perceive its existence, its location information cannot be obtained.

pp ii 11 (( xx )) == dd ii 11 coscos θθ ii 11 pp ii 11 (( ythe y )) == dd ii 11 sinsin θθ ii 11 -- -- -- (( 11 ))

显然,由于个体的观测范围是有限的,三个体实现编队E结构的前提是初始分布时三个体互相可观测。因此,三个体分布时应该充分靠近,以满足彼此可观测的条件,此处分析过程假设了初始分布三邻近个体互相可观测。  Obviously, due to the limited observation range of individuals, the prerequisite for the three entities to realize the formation E structure is that the three entities are mutually observable in the initial distribution. Therefore, the distribution of the three individuals should be sufficiently close to meet the condition of mutual observability. The analysis process here assumes that the initial distribution of the three adjacent individuals is mutually observable. the

个体通过观测和简单运算得到其余个体相对自己的局部位置信息,然后执行预先按照任务需求所设计好的交互控制方法,计算出下一时刻的目的地位置,并将此目的地位置带入运动控制方程从而得到下一时刻所需运动状态(即速度,包括大小和方向)。  The individual obtains the local position information of other individuals relative to itself through observation and simple calculation, and then executes the interactive control method designed in advance according to the task requirements, calculates the destination position at the next moment, and brings this destination position into the motion control The equation thus obtains the required motion state (ie speed, including size and direction) at the next moment. the

下面说明后续用到的符号所表示的意义以及相关定义和问题陈述。  The following explains the meanings of the symbols used later, as well as related definitions and problem statements. the

对于三个邻近个体,任意个体及其位置分别记为Ai和pi,其余两个个体及其位置分别记为Ai1、Ai2和pi1、pi2,并记NPi={pi1,pi2}为邻居位置集合。三邻近个体的索引号集合记为ID={i,i1,i2}。  For three adjacent individuals, any individual and its position are recorded as A i and p i respectively, and the other two individuals and their positions are respectively recorded as A i1 , A i2 and p i1 , p i2 , and NP i ={p i1 ,p i2 } is the set of neighbor positions. The set of index numbers of three adjacent individuals is recorded as ID={i,i1,i2}.

定义1(三点结构,T):由三近邻个体分布位置决定的几何结构(通常为三角形),记为Ti={pi,pi1,pi2},其中,T中三元素可处于同一直线。  Definition 1 (three-point structure, T): a geometric structure (usually triangular) determined by the distribution positions of three neighbors, recorded as T i ={p i ,p i1 ,p i2 }, where the three elements in T can be in the same line.

定义2(等边三角形,E):当任意两元素之间的距离都为D时的三点结构称为等边三角形,记为E(此时D称为E的边长)。  Definition 2 (equilateral triangle, E): When the distance between any two elements is D, the three-point structure is called an equilateral triangle, denoted as E (at this time D is called the side length of E). the

在T和E定义的基础上,定义三邻近个体的E结构编队行为如下。  On the basis of the definition of T and E, define the formation behavior of the E structure of three adjacent individuals as follows. the

三近邻个体E结构编队:三近邻个体从任意初始分布(满足互相可观测)自组织编队成E结构队形的过程,即以三近邻个体为顶点的T结构随着个体间的运动协调转换为E结构的过程。  Three-nearest-neighbor individual E-structure formation: the process in which three-nearest-neighbor individuals self-organize into an E-structure formation from any initial distribution (satisfying mutual observability), that is, the T-structure with the three-nearest-neighbor individuals as vertices transforms into The process of E structure. the

2.本实施例的运动控制方程  2. The motion control equation of this embodiment

在以上状态转移模型假设的基础上,我们提出一种以连续离散时间表示的个体运动控制方程。个体以离散时间表示的形式在每一采样时刻计算出下一时刻的目的地位置,两个采样时刻的间隔内则保持目的地位置不变。下面给出个体以连续时间表示的运动控制方程:  Based on the assumptions of the above state transition model, we propose a control equation for individual motion expressed in continuous discrete time. The individual calculates the destination position of the next time at each sampling time in the form of discrete time representation, and keeps the destination position unchanged during the interval between two sampling times. The motion control equation of the individual expressed in continuous time is given below:

pp ·· ii == ff (( pp ii -- pp igig )) -- -- -- (( 22 ))

f(pi-pig)=-C·(pi-pig)    (3)  f(p i -p ig )=-C·(p i -p ig ) (3)

其中,pi和pig分别表示个体i的当前位置和下一时刻的目的地位置;负号表示个体i朝着目的地位置的方向运动;常数C由实际个体的最大运动速率|vmax|以及观测半径rs共同决定,具体可由如下推导过程得到。  Among them, p i and p ig respectively represent the current position of individual i and the destination position at the next moment; the negative sign indicates that individual i is moving towards the destination position; the constant C is determined by the maximum movement rate of the actual individual |v max | and the observation radius r s are jointly determined, which can be obtained by the following derivation process.

由式(2)及式(3)可知,个体在任何时候运动,都满足以下关系式  From formula (2) and formula (3), it can be seen that the individual moves at any time, and it satisfies the following relationship

CC ·· || pp ii -- pp igig || == || pp ·&Center Dot; ii ||

对某个固定常数C,上式左侧C·|pi-pig|的最大值是C·rs,表示所期望的最大个体即时速率,而上式右侧为机器人的实际即时速率

Figure BDA00003321227500072
显然,前者期望即时速率不能大于实际即时速率的最大值,记为
Figure BDA00003321227500073
因此,合理的常数C值应该满足下式:  For a fixed constant C, the maximum value of C·|p i -p ig | on the left side of the above formula is C·r s , which represents the expected maximum individual instant speed, while the right side of the above formula is the actual instant speed of the robot
Figure BDA00003321227500072
Obviously, the former expected instant rate cannot be greater than the maximum value of the actual instant rate, denoted as
Figure BDA00003321227500073
Therefore, a reasonable constant C value should satisfy the following formula:

CC ≤≤ || vv maxmax || rr sthe s -- -- -- (( 44 ))

确定好控制常数C之后,由运动控制方程可知,下一时刻运动状态取决于目的地位置。而此目的地位置是由个体在每一时刻通过执行交互控制方法来获得的。因此,交互控制方法是个体整个行为决策的关键所在,它应依据三近邻个体的具体全局任务来进行设计。下面将E结构编队作为三近邻个体的全局任务,给出相应的交互控制方法的详细实现过程。  After the control constant C is determined, it can be seen from the motion control equation that the motion state at the next moment depends on the destination position. And this destination position is obtained by the individual at each moment by executing the interactive control method. Therefore, the interactive control method is the key to the individual's overall behavior decision-making, and it should be designed according to the specific global tasks of the three nearest neighbors. In the following, the E-structure formation is regarded as the global task of the three nearest neighbors, and the detailed implementation process of the corresponding interactive control method is given. the

3.本实施的三近邻个体E结构编队  3. The three-nearest neighbor individual E structure formation in this implementation

尽管每个个体中执行的是同一个交互控制方法且它们的行为是异步的,然而三近邻个体要想从任意初始分布通过彼此间的运动协调编队成E结构,每一时刻每个个体的动作应该趋向于有助于完成全局任务目标(E结构编队)。这符合群机器人学中群体中个体的动作应赋予达成全局任务完成的趋向性特征。这里用D表示所期望的E结构边长,显然,D的值应小于个体的观测半径rs(即D<rs),否则会由于彼此观测不到对方而造成无法编队成所希望的E结构。按照以上思想,任意个体Ai应采用的一个很自然的行为决策是:在每一时刻将到个体Ai1和Ai2距离相等且长度为D的位置pig作为自己下一时刻的目的地位置。就E结构编队来说,位置pig是个体Ai下一时刻最理想的去处,这是因为,如果个体Ai处于位置pig就能够与个体Ai1、个体Ai2构成腰长为D的等腰三角形结构。显然满足要求的pig如果存在的话会有两个,这里我们统一的将距离个体Ai较近的pig作为个体Ai下一时刻的目的地位置。每一个时刻保证所有个体都有明确的目的地位置pig的条件是三邻近个体构成的T结构的边长都不大于2D,即需要满足下式:  Although each individual executes the same interactive control method and their behavior is asynchronous, if the three neighbors want to form an E structure from any initial distribution through mutual motion coordination, each individual’s action at each moment Should tend to contribute to the achievement of global mission objectives (E-structure formations). This is consistent with the swarm robotics in which the actions of individuals in the swarm should be endowed with the tendency to achieve global task completion. Here, D is used to represent the desired side length of the E structure. Obviously, the value of D should be smaller than the observation radius r s of the individual (ie D< rs ), otherwise it will not be possible to form the desired E structure because they cannot observe each other. structure. According to the above ideas, a very natural behavior decision for any individual A i is: at each moment, take the position p ig which is equal in distance to individuals A i1 and A i2 and has a length of D as its destination position at the next moment . As far as the E-structure formation is concerned, the position p ig is the most ideal place for the individual A i to go at the next moment, because if the individual A i is in the position p ig , it can form a waist-length D formation with the individual A i1 and the individual A i2 Isosceles triangle structure. Obviously, if there are two p igs that meet the requirements, here we uniformly take the p ig that is closer to the individual A i as the destination position of the individual A i at the next moment. The condition to ensure that all individuals have a clear destination position p ig at each moment is that the side length of the T structure composed of three adjacent individuals is not greater than 2D, that is, the following formula needs to be satisfied:

dj,k≤2D j≠k∈ID    (5)  d j,k ≤2D j≠k∈ID (5)

这样一来,如果后续时刻的T结构仍然满足边长都不大于2D的条件,三个体会在观测、计算、移动三个执行状态之间周而复始的循环,通过彼此间的运动协调最终完成E结构编队,我们将此过程称为调整过程(APr)。在调整过程中,每个个体中执行的是相同的单纯调整算法(AA)。  In this way, if the T structure at the subsequent time still satisfies the condition that the side length is not greater than 2D, the three experiences will cycle through the three execution states of observation, calculation, and movement, and finally complete the E structure through mutual motion coordination Formation, we refer to this process as the adjustment process (APr). During the adjustment process, the same simple adjustment algorithm (AA) is executed in each individual. the

3.1.调整过程(APr)  3.1. Adjustment Process (APr)

图6给出了满足式(5)的所有典型T结构类型和此时任意个体Ai相对应的pig的计算方法(算法1),虚线箭头表示个体Ai下一时刻的期望运动方向。  Figure 6 shows all typical T structure types that satisfy formula (5) and the calculation method of p ig corresponding to any individual A i at this time (algorithm 1). The dotted arrow indicates the expected movement direction of the individual A i at the next moment.

当某采样时刻,三近邻个体之间的距离都一致的满足式(5)时,所有个体都能为下一时刻计算出一个明确的目的地位置pig。以下两式用数学语言描述了个体Ai的目的地位置pig的含义。式中符号||·||为2-范数,表示平面中任意两个位置之间的距离。表1给出了与调整过程相对应的算法伪代码。  When the distances among the three nearest neighbors all satisfy the formula (5) at a certain sampling moment, all individuals can calculate a definite destination position p ig for the next moment. The following two formulas describe the meaning of the destination position p ig of individual A i in mathematical language. In the formula, the symbol ||·|| is a 2-norm, which means the distance between any two positions in the plane. Table 1 gives the algorithm pseudo code corresponding to the adjustment process.

P={p|||p-pi1||=D且||p-pi2||=D,p∈R2}    (6)  P={p|||pp i1 ||=D and ||pp i2 ||=D,p∈R 2 } (6)

pp igig == argarg pp &Element;&Element; PP [[ minmin (( || || pp -- pp ii || || )) ]] -- -- -- (( 77 ))

表1  Table 1

Figure BDA00003321227500082
Figure BDA00003321227500082

3.2内聚过程(CPr)  3.2 Cohesive process (CPr)

在某采样时刻,当个体Ai通过观测和计算后获知其它两近邻个体之间的距离满足式(5)时,个体Ai能够为下一时刻计算出一个明确的目的地位置pig,此位置距离个体Ai最近且到其余两个个体的距离同为D。由于此目的地位置与另外两个个体能构成腰长为D的等腰三角形结构,个体Ai朝此目的地位置方向移动显然有利于形成Ei结构。然而,同一采样时刻并不能保证其它两个个体也都能够获得各自下一时刻明确的目的地位置。如图7所示,描述了某时刻三近邻个体Ai、Ai1及Ai2的位置分布情况,di12<2D,个体Ai的目的地位置应为pig(注意:不是p′ig),因为只有pig满足pig,i1=pig,i2=D且距离个体Ai最近。由于个体Ai和Ai1之间相距较远,个体Ai2无法找到满足条件pig,i1=pig,i2=D的明确目的地位置pig。个体Ai利用近距传感器直接获得di,i1、di,i2,利用数字指南针获得方位角θi,i1和θi,i2,就可通过式(9)计算出个体Ai1和Ai2之间的距离di1,i2 At a certain sampling moment, when the individual A i knows that the distance between other two neighbor individuals satisfies the formula (5) through observation and calculation, the individual A i can calculate a definite destination position p ig for the next moment, where The location is the closest to individual A i and the distance to the other two individuals is the same as D. Since this destination location and two other individuals can form an isosceles triangle structure with waist length D, it is obviously beneficial for individual A i to move towards this destination location to form a structure E i . However, the same sampling moment does not guarantee that the other two individuals will also be able to obtain their respective definite destination positions at the next moment. As shown in Figure 7, it describes the position distribution of three neighbor individuals A i , A i1 and A i2 at a certain moment, d i12 <2D, the destination position of individual A i should be p ig (note: not p′ ig ) , because only p ig satisfies p ig,i1 =p ig,i2 =D and is closest to individual A i . Due to the large distance between individuals A i and A i1 , individual A i2 cannot find a definite destination p ig satisfying the condition p ig,i1 =p ig,i2 =D. The individual A i directly obtains d i,i1 and d i,i2 by using the proximity sensor, and obtains the azimuth angle θ i,i1 and θ i,i2 by using the digital compass, then the individual A i1 and A i2 can be calculated by formula (9) The distance between d i1,i2

dd ii 11 ,, ii 22 == dd ii ,, ii 11 22 ++ dd ii ,, ii 22 22 -- 22 &CenterDot;&CenterDot; dd ii ,, ii 11 &CenterDot;&CenterDot; dd ii ,, ii 22 &CenterDot;&Center Dot; coscos (( &theta;&theta; ii ,, ii 22 -- &theta;&theta; ii ,, ii 11 )) -- -- -- (( 88 ))

对于三近邻个体之间的间距不全都满足式(5)时,如何让所有个体在下一时刻都有恰当的目的地位置去向,我们提出如下内聚策略:在某采样时刻,个体如果发现三个体彼此之间的距离不全都满足式(5)时,便停止对明确的目的地位置的计算,转而计算三者的平均位置并将其做为下一时刻的近似目的地位置。表2给出此算法的伪代码。内聚策略本身其实对三个体的位置分布情况并没有特别的要求,因为三个体的平均位置总能够利用三个体的位置计算来获得。这里是将内聚策略视为整个交互控制策略的一部分,给内聚策略加上了一个位置分布执行条件。这主要是为了在三近邻个体彼此间的距离不是全都满足式(5)的情况下实现三者的充分靠近,保证彼此之间的间距最终全都满足式(5)。个体从此会转入调整过程,三近邻个体通过持续的行为协调最终实现E结构编队。  When the distances between the three neighbors do not all satisfy the formula (5), how to make all the individuals have the appropriate destination at the next moment, we propose the following cohesion strategy: at a certain sampling time, if the individual finds three individuals When the distances between them do not all satisfy the formula (5), the calculation of the definite destination position is stopped, and the average position of the three is calculated instead, and it is used as the approximate destination position at the next moment. Table 2 gives the pseudo code of this algorithm. The cohesion strategy itself has no special requirements on the position distribution of the three bodies, because the average position of the three bodies can always be obtained by calculating the positions of the three bodies. Here, the cohesion strategy is regarded as a part of the entire interaction control strategy, and a location distribution execution condition is added to the cohesion strategy. This is mainly to realize the closeness of the three neighbors when the distances between the three neighbors do not all satisfy the formula (5), so as to ensure that the distances between each other finally all satisfy the formula (5). Individuals will then enter the adjustment process, and the three-nearest neighbors will finally realize the E-structure formation through continuous behavioral coordination. the

表2  Table 2

Figure BDA00003321227500091
Figure BDA00003321227500091

由以上可知,内聚过程作为整个交互控制过程的一部分,主要是在当三近邻个体初始位置彼此拉开很大时用以实现三者向其平均位置充分靠近。当三个体分布已足够靠拢以至于彼此之间的间距全都满足式(5)时,个体便直接转入调整过程。调整过程直接面向目标结构,以使三近邻机器人编队成期望的E结构。将相应的两个算法独立开来看,其实内聚算法(CA)并没有对机器人位置分布有特定要求,执行单纯CA算法的三近邻个体会聚集到同一位置。调整算法(AA)有对位置分布的特定要求即三个体任意两者之间的距离都一致的满足式(5)时才会被执行,并驱使三近邻个体直接编队成E结构,否则会发生目的地位置无法计算所导致的严重错误或者让算法退出执行。图8给出了个体Ai分别执行两个单纯算法时的状态流程示意图。  It can be seen from the above that the cohesion process, as part of the entire interactive control process, is mainly used to achieve the three neighbors sufficiently close to their average positions when the initial positions of the three neighbors are far apart from each other. When the distribution of the three bodies is close enough so that the distances between them all satisfy the formula (5), the individual will directly enter the adjustment process. The adjustment process is directly oriented to the target structure, so that the three-nearest neighbor robots form into the desired E-structure. Looking at the corresponding two algorithms independently, in fact, the cohesion algorithm (CA) does not have specific requirements for the robot position distribution, and the three nearest neighbors who execute the simple CA algorithm will gather at the same position. The adjustment algorithm (AA) has a specific requirement for the position distribution, that is, the distance between any two of the three bodies is the same, and it will be executed only when the formula (5) is satisfied, and it will drive the three nearest neighbors to form an E structure directly, otherwise it will happen A critical error that causes the destination location to fail to be calculated or cause the algorithm to abort execution. Fig. 8 shows a schematic diagram of the state flow when individual A i respectively executes two simple algorithms.

3.3交互控制方法(TFS)  3.3 Interactive control method (TFS)

个体根据当前自己与其它近邻个体的相对位置结构来选择首先进入内聚过程还是调整过程。个体如果首先进入的是面向目标结构的调整过程,则个体将不会再进入内聚过程,反之如果首先进入的是内聚过程,则在三近邻个体取得充分内聚之后还是会进入调整过程,以最终完成E结构编队。我们将个体在整个编队过程中的个体所采用的行为决策称为三角形编队策略(Triangle Formation Strategy),所对应的交互控制方法记为TFS。所有个体执行的是同一 个交互控制方法,它们的动作及在其上执行TFS算法的过程是异步的。为便于对TFS算法的进一步描述,下面给出两种相对位置结构的定义用以概括所有可能的三点结构类型。  Individuals choose whether to first enter the cohesion process or the adjustment process according to the current relative position structure between themselves and other neighbors. If the individual first enters the adjustment process facing the target structure, the individual will not enter the cohesion process again; on the contrary, if the individual enters the cohesion process first, it will still enter the adjustment process after the three neighbors have achieved sufficient cohesion. To finally complete the E structure formation. We refer to the behavior decisions adopted by individuals in the entire formation process as Triangle Formation Strategy, and the corresponding interactive control method is denoted as TFS. All individuals execute the same interactive control method, and their actions and the process of executing the TFS algorithm on them are asynchronous. In order to facilitate the further description of the TFS algorithm, the definitions of two relative position structures are given below to summarize all possible three-point structure types. the

定义(内聚结构,CF):至少存在一条大于2D边长的T结构。  Definition (cohesive structure, CF): There exists at least one T structure larger than the 2D side length. the

定义(调整结构,AF):边长都不大于2D的T结构。  Definition (Adjustment Structure, AF): T structure whose side length is not greater than 2D. the

由以上定义可知,任意时刻三近邻个体构成的T结构只可能是属于以上定义的两种类型结构之一,即内聚结构(CF)或调整结构(AF)。个体位于CF结构中时会执行CA算法,而处于AF结构中时会执行AA算法,但两种算法的共同目的都是为了实现三近邻个体的E结构编队。个体通过执行CA算法向三个体的平均位置移动以实现彼此间的充分内聚,从而使得执行AA算法的位置分布条件得以满足。换句话说,执行内聚算法的目的是保证三个体能够转入AA算法的执行。要想达到目标E结构的成功编队,AA算法是每个个体最终必需执行的TFS算法程序段。而对于TFS算法的CA算法部分,个体未必会执行,这取决于当时三近邻个体位置分布结构是否为CF结构。  From the above definition, it can be known that the T structure formed by three neighbors at any time can only be one of the two types of structures defined above, that is, cohesive structure (CF) or adjusted structure (AF). The CA algorithm will be executed when the individual is in the CF structure, and the AA algorithm will be executed when it is in the AF structure, but the common purpose of the two algorithms is to realize the formation of the E structure of the three nearest neighbors. Individuals move to the average position of the three bodies by executing the CA algorithm to achieve sufficient cohesion among each other, so that the location distribution conditions for executing the AA algorithm are satisfied. In other words, the purpose of executing the cohesion algorithm is to ensure that the three bodies can be transferred to the execution of the AA algorithm. In order to achieve the successful formation of the target E structure, the AA algorithm is the program segment of the TFS algorithm that each individual must eventually execute. For the CA algorithm part of the TFS algorithm, individuals may not be able to execute it, which depends on whether the distribution structure of the three nearest neighbors is a CF structure. the

因此,用于三近邻个体自组织E结构编队的TFS交互控制方法综合了单纯CA算法和单纯AA算法,个体选择哪个算法程序开始执行,要视当时三者的位置分布结构。表3和表4分别给出了TFS交互控制方法的伪代码及相应的程序段(用于目的地位置pig的求取)。个体执行TFS算法时的状态流程如图9所示,流程图中虚线框所表示的修正子算法部分在TFS算法用于实体系统时是需要加以考虑的。  Therefore, the TFS interactive control method for the self-organized E-structure formation of three-neighboring individuals combines the simple CA algorithm and the simple AA algorithm. Which algorithm program the individual chooses to start execution depends on the position distribution structure of the three at that time. Table 3 and Table 4 respectively give the pseudo-code of the TFS interactive control method and the corresponding program segment (for the calculation of the destination position p ig ). Figure 9 shows the state flow of an individual executing the TFS algorithm. The modified sub-algorithm part indicated by the dotted box in the flow chart needs to be considered when the TFS algorithm is used in the entity system.

表3  table 3

Figure BDA00003321227500101
Figure BDA00003321227500101

表4  Table 4

Figure BDA00003321227500102
Figure BDA00003321227500102

Figure BDA00003321227500111
Figure BDA00003321227500111

Figure BDA00003321227500121
Figure BDA00003321227500121

Figure BDA00003321227500131
Figure BDA00003321227500131

对于实体机器人,当三个体位置共线分布时,处于两端的个体由于受中间个体遮挡的影响,彼此之间无法观测到对方,都只能观测到中间个体。而位于中间的个体则能同时观测到位于两端的个体,如图10所示。值得提醒的是,这里只讨论当三个体位置共线时,它们仍处于彼此之间的观测能力区域以内的情况。图中阴影条状区域就是由于中间个体遮挡所造成的个体Ai1的观测盲区,Ai2正好处于此盲区内,Ai1观测不到Ai2。类似的原因,Ai2也观测不到Ai1。  For a physical robot, when the positions of the three bodies are distributed collinearly, the individuals at both ends cannot observe each other due to the occlusion of the middle individual, and can only observe the middle individual. While the individual in the middle can observe the individuals at both ends at the same time, as shown in Figure 10. It is worth reminding that here we only discuss the situation that when the positions of the three bodies are collinear, they are still within the observation capability area of each other. The shaded strip area in the figure is the observation blind area of individual A i1 caused by the occlusion of the middle individual. A i2 is just in this blind area, and A i1 cannot observe A i2 . For similar reasons, A i2 cannot observe A i1 either.

因此,可以对在实体机器人中执行的TFS交互控制方法引入如图11所示的修正子算法。个体当观测到单个个体时,就暂时停留在原位置不动(注意:系统仍在运行,执行观测和计算过程)。当观测到两个邻居且发现与它们的位置共线时,可以判断出自己处于这两邻居的中间位置,就沿垂直共线的任意一侧方向移动。这样一来位于两端的个体就能彼此感知到,从而引起个体对TFS交互控制方法(的CA算法或AA算法部分)的继续执行。  Therefore, the modified sub-algorithm shown in Figure 11 can be introduced into the TFS interactive control method implemented in the physical robot. When an individual observes a single individual, it temporarily stays in place (note: the system is still running, performing observation and calculation processes). When two neighbors are observed and found to be collinear with their positions, you can judge that you are in the middle of the two neighbors, and move along any side of the vertical collinear. In this way, the individuals at both ends can perceive each other, thereby causing the individuals to continue to execute the TFS interactive control method (the CA algorithm or the AA algorithm part). the

4.E结构编队稳定性分析  4. E-structure formation stability analysis

以上详细描述了实现三近邻个体协调行为E结构编队的TFS交互控制方法的设计过程及有关算法的伪代码和程序段。下面就三近邻个体利用TFS交互控制方法进行E结构编队的有关性质给出相关的结论:  The above describes in detail the design process of the TFS interactive control method to realize the E-structure formation of the three-nearest neighbor individual coordinated behavior, as well as the pseudo-code and program segments of the related algorithm. The following conclusions are given on the relevant properties of the E-structure formation of the three nearest neighbors using the TFS interactive control method:

引理4.1:三近邻个体在内聚过程中的平均位置

Figure BDA00003321227500132
不变,即有:  Lemma 4.1: The average position of three neighbors in the cohesion process
Figure BDA00003321227500132
unchanged, that is:

pp &OverBar;&OverBar; &CenterDot;&CenterDot; == 00

证明:  prove:

三个体参与内聚过程时执行的是CA算法,如图8所示,由算法定义有

Figure 174709DEST_PATH_GDA00003465689600141
并由运动控制方程(式(2)和式(3))可知  When the three bodies participate in the cohesion process, the CA algorithm is executed, as shown in Figure 8, defined by the algorithm
Figure 174709DEST_PATH_GDA00003465689600141
And from the motion control equation (type (2) and type (3)), it can be known that

pp &CenterDot;&CenterDot; ii == -- CC &CenterDot;&CenterDot; (( pp ii -- pp igig )) == -- CC &CenterDot;&CenterDot; (( pp ii -- pp &OverBar;&OverBar; ))

pp &CenterDot;&CenterDot; ii 11 == -- CC &CenterDot;&CenterDot; (( pp ii 11 -- pp ii 11 gg )) == -- CC &CenterDot;&CenterDot; (( pp ii 11 -- pp &OverBar;&OverBar; ))

pp &CenterDot;&CenterDot; ii 22 == -- CC &CenterDot;&CenterDot; (( pp ii 22 -- pp ii 22 gg )) == -- CC &CenterDot;&CenterDot; (( pp ii 22 -- pp &OverBar;&OverBar; ))

由于对此式两边求导后,将以上式子代入可得  because After deriving both sides of this formula, we can substitute the above formula into

pp &OverBar;&OverBar; &CenterDot;&Center Dot; == (( pp &CenterDot;&Center Dot; ii ++ pp &CenterDot;&Center Dot; ii 11 ++ pp &CenterDot;&Center Dot; ii 22 33 )) == (( -- CC &CenterDot;&Center Dot; (( pp ii -- pp &OverBar;&OverBar; )) -- CC &CenterDot;&Center Dot; (( pp ii 11 -- pp &OverBar;&OverBar; )) -- CC &CenterDot;&Center Dot; (( pp ii 22 -- pp &OverBar;&OverBar; )) 33 ))

== -- CC &CenterDot;&CenterDot; (( pp ii ++ pp ii 11 ++ pp ii 22 33 -- pp &OverBar;&OverBar; )) == -- CC &CenterDot;&CenterDot; (( pp &OverBar;&OverBar; -- pp &OverBar;&OverBar; )) == 00

命题得证。  The proposition is proved. the

定理4.1:通过单纯CA算法三近邻个体经过时间T后将稳定内聚在以平均位置

Figure BDA00003321227500147
为中心的η邻域内:  Theorem 4.1: Through the simple CA algorithm, the three nearest neighbors will be stable and cohesive at the average position after time T
Figure BDA00003321227500147
In the η neighborhood of the center:

BB &eta;&eta; (( pp &OverBar;&OverBar; )) == {{ pp || || || pp -- pp &OverBar;&OverBar; || || << &eta;&eta; }}

TT == 11 CC lnln (( 22 EE. &eta;&eta; ))

其中η为任意正数, E = max ( | | p i ( 0 ) - p &OverBar; | | : i &Element; ID ) . Wherein n is any positive number, E. = max ( | | p i ( 0 ) - p &OverBar; | | : i &Element; ID ) .

证明:  prove:

由引理4.1可知,执行CA算法的三个体的平均位置满足

Figure BDA000033212275001411
From Lemma 4.1, we can see that the average position of the three bodies performing the CA algorithm satisfies
Figure BDA000033212275001411

对机器人i,i∈ID,取误差变量为ei(t)=pi-pig,取能量函数为

Figure BDA000033212275001412
并对能量函数求导,可得:  For robot i, i∈ID, the error variable is e i (t)=p i -p ig , and the energy function is
Figure BDA000033212275001412
And deriving the energy function, we can get:

VV &CenterDot;&CenterDot; ii (( tt )) == (( ee &CenterDot;&Center Dot; ii (( tt )) )) TT ee ii (( tt ))

== (( pp &CenterDot;&Center Dot; ii -- pp &CenterDot;&Center Dot; igig )) TT (( pp ii -- pp igig ))

= ( p &CenterDot; i - p &OverBar; &CenterDot; ) T ( p i - p &OverBar; ) (由内聚过程有 p ig = p &OverBar; = ( p &Center Dot; i - p &OverBar; &Center Dot; ) T ( p i - p &OverBar; ) (by the cohesive process has p ig = p &OverBar; )

== (( pp &CenterDot;&Center Dot; ii )) TT (( pp ii -- pp &OverBar;&OverBar; ))

== -- CC (( pp ii -- pp &OverBar;&OverBar; )) TT (( pp ii -- pp &OverBar;&OverBar; ))

== -- CC || || pp ii -- pp &OverBar;&OverBar; || || 22

&le;&le; 00

由于对任意正数η,当

Figure BDA00003321227500157
时,恒有
Figure BDA00003321227500158
由Lyapunov稳定性理论可知,个体经过足够长的时间T必会稳定内聚到以
Figure BDA00003321227500159
为中心的η邻域内。  Since for any positive number η, when
Figure BDA00003321227500157
always have
Figure BDA00003321227500158
From the Lyapunov stability theory, it can be known that after a long enough time T, the individual will be stable and cohesive to the following
Figure BDA00003321227500159
In the η neighborhood of the center.

又由 V &CenterDot; i ( t ) = - C | | p i - p &OverBar; | | 2 = - C ( e i ( t ) ) T e i ( t ) = - 2 C V i ( t ) , 求解此微分方程可得Vi(t)=Vi(0)e-2Ct。  And by V &Center Dot; i ( t ) = - C | | p i - p &OverBar; | | 2 = - C ( e i ( t ) ) T e i ( t ) = - 2 C V i ( t ) , Solving this differential equation gives V i (t)=V i (0)e -2Ct .

对任意正数η,当

Figure BDA000033212275001511
时,恒有当||ei(t)||=η时,表示个体i已经进入了
Figure BDA000033212275001513
的η邻域,用ti表示个体i内聚到此η邻域内所需要的时间,则有Vi(ti)=η2/2,所以
Figure BDA000033212275001514
解得
Figure BDA000033212275001515
因此,三近邻个体都内聚到以
Figure BDA000033212275001516
为中心的η邻域内所需要的时间为:  For any positive number η, when
Figure BDA000033212275001511
always have When ||e i (t)||=η, it means that individual i has entered
Figure BDA000033212275001513
η neighborhood, use t i to represent the time required for individual i to cohere into this η neighborhood, then V i (t i )=η 2 /2, so
Figure BDA000033212275001514
Solutions have to
Figure BDA000033212275001515
Therefore, the three nearest neighbors are all cohesive to the following
Figure BDA000033212275001516
The time required in the η neighborhood of the center is:

TT == maxmax (( tt ii :: ii &Element;&Element; IDID ))

== maxmax (( 11 CC lnln (( 22 VV ii (( 00 )) &eta;&eta; )) :: ii &Element;&Element; IDID ))

== maxmax (( 11 CC lnln (( 22 || || pp ii (( 00 )) -- pp &OverBar;&OverBar; || || &eta;&eta; )) :: ii &Element;&Element; IDID ))

= 1 C ln ( 2 E &eta; ) , 其中 E = max ( | | p i ( 0 ) - p &OverBar; | | : i &Element; ID ) = 1 C ln ( 2 E. &eta; ) , in E. = max ( | | p i ( 0 ) - p &OverBar; | | : i &Element; ID )

命题得证。  The proposition is proved. the

引理4.2:三近邻个体通过CA算法能从任意CF结构转变到AF结构且对于TFS交互控制方法不可逆。  Lemma 4.2: The three nearest neighbors can transform from any CF structure to AF structure through the CA algorithm, and it is irreversible for the TFS interactive control method. the

证明:  prove:

由定理4.1可知,存在一个正数η以让三个体充分内聚使得彼此之间的距离足够小以至于max(||pi-pj||:i,j∈ID且i≠j)≤2D成立,因此,命题前半部分得证。  According to Theorem 4.1, there is a positive number η to make the three bodies sufficiently cohesive so that the distance between them is small enough that max(||p i -p j ||:i,j∈ID and i≠j)≤ 2D holds, therefore, the first half of the proposition is proved.

对于TFS算法,AF结构不会转换为CF结构。三邻近个体位于AF结构中时,参与的是调整过程,执行的是AA算法。AF结构要想转换为CF结构,必须经过一个临界结构,如图12所示,此时有max(||pi-pj||:i,j∈ID且i≠j)=2D。考虑到对称性,只分析个体Ai位于由两圆及连接个体Ai1和Ai2的直线所围成区域中的情况(注:图中以Ai1、Ai2为中心半径为2D的圆边界是用于辅助分析的而不是表示个体的观测边界)。  For the TFS algorithm, AF structures are not converted to CF structures. When the three adjacent individuals are in the AF structure, they participate in the adjustment process and execute the AA algorithm. If the AF structure is to be transformed into a CF structure, it must go through a critical structure, as shown in Figure 12. At this time, max(||p i -p j ||:i,j∈ID and i≠j)=2D. Considering the symmetry, only analyze the case that the individual A i is located in the area surrounded by two circles and the straight line connecting the individual A i1 and A i2 (Note: In the figure, A i1 and A i2 are the circle boundaries with a radius of 2D is used to aid in the analysis rather than represent the individual observation boundaries).

显然,如果AF结构要想经历这种临界结构转换到CF结构,个体Ai1和Ai2必须要发生背离运动。同样考虑到对称性特点,这里只对个体Ai1作讨论,个体Ai2的情况可类似讨论。如果背离运动发生,个体Ai1的运行方向必须是朝向水平线l1上方(注:l1与Ai1Ai2垂直),也就是说Ai1的目的地位置必须位于此水平线以上区域。这里假设在水平线l1上方区域的任意位置pi1g就是个体Ai1的目的地位置。从上图可知,目的地位置pi1g到个体Ai2的距离必然大于2D,由于此时个体还是处于调整结构,个体进行的是调整过程,执行的是AA算法。又由调整过程定义知,目的地位置pi1g到个体Ai2的距离不会大于等于2D,引起矛盾。从而得出,个体Ai1的目的地位置pi1g位于水平线l1上方区域的假设不成立,即个体Ai1的目的地位置pi1g必处于水平线l1下方区域。类似分析可得,个体Ai2的目的地位置pi2g必处于水平线l2上方区域。因此,个体Ai1和Ai2不会发生背离运动,即在TFS算法下,AF结构到CF结构的转换不会发生,命题后半部分得证。  Apparently, if the AF structure wants to undergo this transition from the critical structure to the CF structure, the individual A i1 and A i2 must move away from each other. Considering the characteristic of symmetry, here only individual A i1 is discussed, and the situation of individual A i2 can be discussed similarly. If the departure movement occurs, the running direction of the individual A i1 must be above the horizontal line l1 (note: l1 is perpendicular to A i1 and A i2 ), that is to say, the destination position of A i1 must be located above the horizontal line. Here it is assumed that any position p i1g in the area above the horizontal line l1 is the destination position of individual A i1 . It can be seen from the above figure that the distance between the destination position p i1g and the individual A i2 must be greater than 2D. Since the individual is still in the adjustment structure at this time, the individual is performing the adjustment process and the AA algorithm is executed. It is also known from the definition of the adjustment process that the distance from the destination position p i1g to the individual A i2 will not be greater than or equal to 2D, causing a contradiction. Therefore, it can be concluded that the assumption that the destination position p i1g of individual A i1 is located in the area above the horizontal line l1 is not valid, that is, the destination position p i1g of individual A i1 must be in the area below the horizontal line l1 . Similar analysis shows that the destination position p i2g of individual A i2 must be in the area above the horizontal line l 2 . Therefore, individual A i1 and A i2 will not deviate, that is, under the TFS algorithm, the conversion from AF structure to CF structure will not occur, and the second half of the proposition is proved.

命题得证。  The proposition is proved. the

定理4.2:三近邻个体通过AA算法能从任意AF结构实现E结构编队。  Theorem 4.2: The three nearest neighbors can realize E-structure formation from any AF structure through AA algorithm. the

证明:  prove:

对任意个体i,i∈ID,由于在采样间隔内个体的目的地位置不变,即有

Figure BDA00003321227500161
取误差函数为ei(t)=pi-pig,取能量函数为并对此能量函数求导,可得:  For any individual i, i∈ID, since the destination position of the individual does not change within the sampling interval, that is,
Figure BDA00003321227500161
Take the error function as e i (t)=p i -p ig , and take the energy function as And taking the derivative of this energy function, we can get:

VV &CenterDot;&Center Dot; ii (( tt )) == (( ee &CenterDot;&Center Dot; ii (( tt )) )) TT ee ii (( tt ))

== (( pp &CenterDot;&Center Dot; ii -- pp &CenterDot;&Center Dot; igig )) TT (( pp ii -- pp igig ))

== (( pp &CenterDot;&Center Dot; ii )) TT (( pp ii -- pp igig ))

== -- CC (( pp ii -- pp igig )) TT (( pp ii -- pp igig ))

== -- CC || || pp ii -- pp igig || || 22

&le;&le; 00

可知,当||pi-pig||2>0,恒有

Figure BDA00003321227500172
按照Lyapunov稳定性理论,pi必渐进收敛到pig,即pi=pig。  It can be seen that when ||p i -p ig || 2 >0, there is always
Figure BDA00003321227500172
According to the Lyapunov stability theory, p i must asymptotically converge to p ig , that is, p i = p ig .

又由调整过程的定义,可知  From the definition of the adjustment process, it can be seen that

||pig-pj||=D,j∈ID\{i}  ||p ig -p j ||=D,j∈ID\{i}

进一步,将pi=pig代入上式,可得  Further, substituting p i =p ig into the above formula, we can get

||pi-pj||=D,j∈ID\{i}  ||p i -p j ||=D,j∈ID\{i}

即任意个体i到其两个邻居的距离均为D。因此,三近邻个体组成的AF结构最终转换为E结构。  That is, the distance from any individual i to its two neighbors is D. Therefore, the AF structure composed of three neighbors is finally transformed into the E structure. the

命题得证。  The proposition is proved. the

定理4.3:三近邻个体通过TFS算法能从任意初始分布实现E结构编队。  Theorem 4.3: The three-nearest neighbor individuals can realize the E-structure formation from any initial distribution through the TFS algorithm. the

证明:  prove:

如果三近邻个体初始位置分布为CF结构,按照TFS交互控制方法的设计,个体会执行CA算法进入内聚过程,由引理4.2可知三近邻个体通过一个内聚过程会从CF结构转换为AF结构,且这种转换是不可逆的。AF结构一旦形成,个体会立刻进入调整过程,由定理4.2可知三近邻个体必会编队成E结构。如果三近邻个体初始位置分布为AF结构,由TFS交互控制方法流程(见图9)可知,个体会直接执行AA算法进入调整过程。又由定理4.2可知,进入调整过程的三近邻个体会直接编队成E结构。  If the initial position distribution of the three neighbors is a CF structure, according to the design of the TFS interactive control method, the individuals will execute the CA algorithm to enter the cohesion process. From Lemma 4.2, it can be seen that the three neighbors will convert from the CF structure to the AF structure through a cohesion process , and this conversion is irreversible. Once the AF structure is formed, the individuals will immediately enter the adjustment process. From Theorem 4.2, it can be seen that the three nearest neighbors will form an E structure. If the initial position distribution of the three nearest neighbors is an AF structure, it can be seen from the flow of the TFS interactive control method (see Figure 9) that the individual will directly execute the AA algorithm to enter the adjustment process. It can also be seen from Theorem 4.2 that the three nearest neighbors entering the adjustment process will directly form an E structure. the

命题得证。  The proposition is proved. the

5.仿真实验  5. Simulation

下面对三近邻个体通过TFS交互控制方法编队E结构的行为进行仿真分析,以验证TFS算法对于三近邻个体E结构编队的有效性。此处,所有实验程序的编制及其运行都是在Breve仿真平台中进行。Breve是基于3D环境的用于分布式系统和人工生命的多智能体仿真开发平台,使用面向对象语言Steve做为平台的仿真程序设计语言。Breve平台中包含一个常用类的类库,这为实验人员实现快速模型仿真及算法验证提供了方便。  Next, the simulation analysis of the behavior of the three-nearest-neighboring individuals forming the E-structure through the TFS interactive control method is carried out to verify the effectiveness of the TFS algorithm for the formation of the E-structure of the three-nearest-neighboring individuals. Here, the compilation and operation of all experimental programs are carried out on the Breve simulation platform. Breve is a multi-agent simulation development platform for distributed systems and artificial life based on 3D environment. It uses the object-oriented language Steve as the simulation programming language of the platform. The Breve platform contains a class library of commonly used classes, which provides convenience for experimenters to realize fast model simulation and algorithm verification. the

仿真参数设置:E结构边长D=10;个体观测半径rs=50;最大运动速率|vmax|=10,由式(7)我们取运动控制方程中的常数C=vmax/rs=0.2;仿真开始时,三个体随机分布在以环境坐标原点为中心,半径为rd=rs/2的圆形区域内,在此范围内的个体都能彼此观测到对方。另外,对仿真运行的时长进行了限制,当机器人的速率|v|满足|v|<0.01时停止运动,系统结 束仿真。从运动控制方程可知,这种对速率的限制其实是指定了个体与其目的地位置的接近程度,即目的地位置距离误差。此处相当于将特定的目的地距离误差值作为了个体移动及E结构编队仿真的停止条件。  Simulation parameter setting: E structure side length D=10; individual observation radius r s =50; maximum motion velocity |v max |=10, from formula (7) we take the constant C=v max /r s in the motion control equation =0.2; at the beginning of the simulation, the three bodies are randomly distributed in a circular area centered on the origin of the environment coordinates and the radius is r d =r s /2, and the individuals within this range can observe each other. In addition, the running time of the simulation is limited. When the velocity |v| of the robot satisfies |v|<0.01, the robot stops moving and the system ends the simulation. It can be seen from the motion control equation that this limitation on the speed actually specifies the closeness of the individual to its destination position, that is, the distance error of the destination position. Here, it is equivalent to using a specific destination distance error value as the stopping condition of individual movement and E-structure formation simulation.

对三近邻个体初始分布为CF和AF的两种情形,分别进行E编队行为仿真和平均边长误差的收敛性测试。平均边长误差由下式定义,用于表示三近邻个体趋近E结构的程度与时间的关系。  For the two cases where the initial distribution of the three nearest neighbors is CF and AF, the E formation behavior simulation and the convergence test of the average side length error are carried out respectively. The average side length error is defined by the following formula, which is used to represent the relationship between the degree of three neighbors approaching the E structure and time. the

Figure BDA00003321227500181
Figure BDA00003321227500181

其中,size(ID)表示ID中的元素个数。  Among them, size(ID) indicates the number of elements in ID. the

对如下典型初始分布进行了E结构编队行为仿真实验:(1)CF分布,CF-1表示初始内聚结构中只有一条边大于2D,具体初始边长为(17.4716.0630.42);CF-2表示初始内聚结构中有两条边大于2D,具体初始边长为(36.6615.7530.56);CF-3表示初始内聚结构中三条边长都大于2D,具体初始边长为(21.9424.4927.00)。(2)AF分布,表示初始调整结构即三条边都满足式(8),具体初始边长为(11.6015.6911.61)。仿真表明:三近邻个体在任意位置分布条件下都能编队成以其为顶点的E结构。图13给出了三近邻个体从CF-2初始分布编队成目标E结构所经历的轨迹情况,图中左上角数字表示以秒为单位所用时间。  The simulation experiment of E-structure formation behavior was carried out for the following typical initial distributions: (1) CF distribution, CF-1 means that only one side in the initial cohesive structure is larger than 2D, and the specific initial side length is (17.4716.0630.42); CF-2 means Two sides in the initial cohesive structure are larger than 2D, and the specific initial side length is (36.6615.7530.56); CF-3 indicates that the three side lengths in the initial cohesive structure are all larger than 2D, and the specific initial side length is (21.9424.4927.00). (2) AF distribution, which means the initial adjustment structure, that is, all three sides satisfy the formula (8), and the specific initial side length is (11.6015.6911.61). The simulation shows that the three nearest neighbors can be formed into an E structure with it as the vertex under any position distribution condition. Figure 13 shows the trajectories experienced by the three neighbors from the initial distribution of CF-2 to the target E structure, and the number in the upper left corner of the figure indicates the time taken in seconds. the

图14给出了在以上几种典型初始分布下三近邻个体动态逼近E结构的过程中平均边长误差随时间的变化情况。从平均误差收敛的曲线可知,无论哪种初始位置分布,T结构所有的边长最后都能收敛到D。仿真实验表明,三近邻个体通过执行TFS交互控制方法能从任意初始分布位置实现E结构编队。  Figure 14 shows the variation of the average side length error over time in the process of the three-neighbor individuals dynamically approaching the E structure under the above typical initial distributions. From the average error convergence curve, it can be seen that no matter what kind of initial position distribution, all the side lengths of the T structure can converge to D in the end. The simulation experiment shows that the three nearest neighbors can realize E-structure formation from any initial distribution position by implementing the TFS interactive control method. the

在上述的实施例中,详细的描述了如何实现三近邻智能机器人编队的控制,通过上述实施例所提供的技术手段,解决了相关技术中多智能体编队行为的分布式控制所存在的控制精度不高以及较大的响应延迟的问题,进而达到了提高多智能体编队的控制精度和加快响应的效果。  In the above-mentioned embodiment, it is described in detail how to realize the control of the three-nearest neighbor intelligent robot formation. Through the technical means provided by the above-mentioned embodiment, the control precision existing in the distributed control of multi-agent formation behavior in the related art is solved. The problem of low and large response delay, and thus achieve the effect of improving the control accuracy of the multi-agent formation and speeding up the response. the

在另外一个实施例中,还提供了一种交互控制的软件,该软件用于执行上述实施例中描述的节点E结构编队的方法。  In another embodiment, interactive control software is also provided, and the software is used to execute the method for node E structure formation described in the above embodiment. the

在另外一个实施例中,还提供了一种存储介质,该存储介质中存储有上述软件,该存储介质包括但不限于光盘、软盘、硬盘、可擦写存储器等。  In another embodiment, there is also provided a storage medium in which the above software is stored, and the storage medium includes but not limited to an optical disk, a floppy disk, a hard disk, a rewritable memory, and the like. the

显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。  Obviously, those skilled in the art should understand that each module or each step of the above-mentioned present invention can be realized by a general-purpose computing device, and they can be concentrated on a single computing device, or distributed in a network formed by multiple computing devices Alternatively, they may be implemented in program code executable by a computing device so that they may be stored in a storage device to be executed by a computing device, and in some cases in an order different from that shown here The steps shown or described are carried out, or they are separately fabricated into individual integrated circuit modules, or multiple modules or steps among them are fabricated into a single integrated circuit module for implementation. As such, the present invention is not limited to any specific combination of hardware and software. the

以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。  The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and changes. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included within the protection scope of the present invention. the

Claims (14)

1.一种节点编队方法,其特征在于,包括:1. A node formation method, characterized in that, comprising: 步骤S102,检测第一节点与第二节点之间的第一距离以及第一方位角,以及检测所述第一节点与第三节点之间的第二距离以及第二方位角,并根据所述第一距离、所述第一方位角、所述第二距离和所述第二方位角计算出所述第二节点与所述第三节点之间的第三距离;Step S102, detecting the first distance and the first azimuth between the first node and the second node, and detecting the second distance and the second azimuth between the first node and the third node, and according to the calculating a third distance between the second node and the third node by the first distance, the first azimuth, the second distance, and the second azimuth; 步骤S104,判断所述第一距离、所述第二距离以及所述第三距离是否均小于或等于预设的第一阈值长度,如果是,则执行步骤S108,如果否,则执行步骤S106;Step S104, judging whether the first distance, the second distance and the third distance are all less than or equal to the preset first threshold length, if yes, execute step S108, if not, execute step S106; 步骤S106,根据所述第一距离、所述第一方位角、所述第二距离和所述第二方位角计算出所述第一节点、所述第二节点和所述第三节点之间的平均位置点,并驱动所述第一节点向所述平均位置点移动,并以预定时间间隔返回执行所述步骤S102;Step S106, calculating the distance between the first node, the second node and the third node according to the first distance, the first azimuth, the second distance and the second azimuth and drive the first node to move to the average position point, and return to execute the step S102 at predetermined time intervals; 步骤S108,根据当前所述第一距离、当前所述第一方位角、当前所述第二距离和当前所述第二方位角计算出目标位置点,其中,所述目标位置点为到所述第二节点以及到所述第三节点的距离均为所述第一阈值长度的一半并且到所述第一节点的距离最小的位置点,并驱动所述第一节点向所述目标位置点移动;Step S108, calculating a target position point according to the current first distance, the current first azimuth angle, the current second distance and the current second azimuth angle, wherein the target position point is to the The distance between the second node and the third node is half of the first threshold length and the minimum distance to the first node, and drives the first node to move to the target position point ; 步骤S110,以预定时间间隔判断所述第一节点是否到达所述目标位置点,如果否,则返回执行所述步骤S102,如果是,则所述第一节点结束编队过程。Step S110, judging whether the first node has reached the target point at predetermined time intervals, if not, return to step S102, and if yes, the first node ends the formation process. 2.根据权利要求1所述的节点编队方法,其特征在于,在所述步骤S102之后还包括:2. The node formation method according to claim 1, further comprising: after the step S102: 当根据当前所述第一方位角和第二方位角判断出所述第一节点、所述第二节点和所述第三节点处于同一直线上,并且所述第一节点位于所述第二节点和所述第三节点的同一侧时,停止驱动所述第一节点。When it is judged according to the current first azimuth angle and the second azimuth angle that the first node, the second node and the third node are on the same straight line, and the first node is located at the second node When it is on the same side as the third node, stop driving the first node. 3.根据权利要求1所述的节点编队方法,其特征在于,在所述步骤S102之后还包括:3. The node formation method according to claim 1, further comprising: after the step S102: 当根据当前所述第一方位角和所述第二方位角判断出所述第一节点、所述第二节点和所述第三节点处于同一直线上,并且所述第一节点位于所述第二节点和所述第三节点的中间时,驱动所述第一节点沿该直线的垂直线的任意一侧移动。When it is determined according to the current first azimuth and the second azimuth that the first node, the second node and the third node are on the same straight line, and the first node is located at the first When the second node is in the middle of the third node, the first node is driven to move along any side of the vertical line of the straight line. 4.根据权利要求1所述的节点编队方法,其特征在于,所述平均位置点为所述第一节点、所述第二节点和所述第三节点所组成的三角形的重心或内心。4. The node formation method according to claim 1, wherein the average position point is the center of gravity or the center of a triangle formed by the first node, the second node and the third node. 5.根据权利要求1所述的节点编队方法,其特征在于,在控制所述第一节点移动的过程中,以同样的方式分别独立地控制所述第二节点和第三节点的移动,以对所述第一、第二和第三节点形成编队。5. The node formation method according to claim 1, characterized in that, in the process of controlling the movement of the first node, the movement of the second node and the third node are independently controlled in the same manner, so as to A formation is formed for the first, second and third nodes. 6.根据权利要求1-5任意一项所述的节点编队方法,其特征在于,所述第一节点、所述第二节点和所述第三节点均为智能机器人。6. The node formation method according to any one of claims 1-5, wherein the first node, the second node and the third node are all intelligent robots. 7.一种节点编队装置,位于第一节点中,其特征在于,包括:7. A node formation device located in the first node, characterized in that it comprises: 检测单元,用于检测所述第一节点与第二节点之间的第一距离以及第一方位角,以及检测所述第一节点与第三节点之间的第二距离以及第二方位角;a detection unit, configured to detect a first distance and a first azimuth between the first node and a second node, and detect a second distance and a second azimuth between the first node and a third node; 计算单元,用于根据所述第一距离、所述第一方位角、所述第二距离和所述第二方位角计算出所述第二节点与所述第三节点之间的第三距离,并判断所述第一距离、所述第二距离以及所述第三距离是否均小于或等于预设的第一阈值长度,如果否,则根据所述第一距离、所述第一方位角、所述第二距离和所述第二方位角计算出所述第一节点、所述第二节点和所述第三节点之间的平均位置点;如果是,则根据所述第一距离、所述第一方位角、所述第二距离和所述第二方位角计算出目标位置点,其中,所述目标位置点为到所述第二节点以及到所述第三节点的距离均为所述第一阈值长度的一半并且到所述第一节点的距离最小的位置点;a calculation unit, configured to calculate a third distance between the second node and the third node according to the first distance, the first azimuth, the second distance, and the second azimuth , and judge whether the first distance, the second distance and the third distance are all less than or equal to the preset first threshold length, if not, according to the first distance, the first azimuth angle , the second distance and the second azimuth to calculate the average position point between the first node, the second node and the third node; if so, according to the first distance, The first azimuth, the second distance, and the second azimuth calculate a target location point, wherein the target location point is a distance from the second node and a distance from the third node. a location point that is half the length of the first threshold and has the smallest distance to the first node; 驱动单元,用于根据所述计算单元的计算结果驱动所述第一节点向所述平均位置点或所述目标位置点移动。A driving unit, configured to drive the first node to move to the average position point or the target position point according to the calculation result of the calculation unit. 8.根据权利要求7所述的节点编队装置,其特征在于,8. The node formation device according to claim 7, characterized in that, 所述计算单元还用于根据当前所述第一方位角和当前所述第二方位角判断所述第一节点、所述第二节点和所述第三节点是否处于同一直线上以及它们之间的相邻关系。The calculation unit is also used to judge whether the first node, the second node, and the third node are on the same straight line and whether they are on the same straight line according to the current first azimuth angle and the current second azimuth angle. neighbor relationship. 9.根据权利要求8所述的节点编队装置,其特征在于,9. The node formation device according to claim 8, characterized in that, 所述驱动单元还用于在所述计算单元判断出所述第一节点、所述第二节点和所述第三节点处于同一直线上并且所述第一节点位于所述第二节点和所述第三节点的同一侧的情况下,停止驱动所述第一节点。The driving unit is further configured to determine that the first node, the second node and the third node are on the same straight line and the first node is located between the second node and the In the case of the same side of the third node, stop driving the first node. 10.根据权利要求8所述的节点编队装置,其特征在于,10. The node formation device according to claim 8, characterized in that, 所述驱动单元还用于在所述计算单元判断出所述第一节点、所述第二节点和所述第三节点处于同一直线上,并且所述第一节点位于所述第二节点和所述第三节点之间的情况下,驱动所述第一节点沿该直线的垂直线的任意一侧移动。The driving unit is further configured to determine that the first node, the second node, and the third node are on the same straight line when the computing unit determines that the first node is located between the second node and the third node. In the case of being between the third node, drive the first node to move along any side of the vertical line of the straight line. 11.根据权利要求7所述的节点编队装置,其特征在于,所述平均位置点为所述第一节点、所述第二节点和所述第三节点所组成的三角形的重心或内心。11 . The node formation device according to claim 7 , wherein the average position point is the center of gravity or incenter of a triangle formed by the first node, the second node and the third node. 12.根据权利要求7-11任一项所述的节点编队装置,其特征在于,所述检测单元包括:12. The node formation device according to any one of claims 7-11, wherein the detection unit comprises: 红外传感器,用于检测所述第一距离和所述第二距离;an infrared sensor for detecting the first distance and the second distance; 数字指南针,用于检测所述第一方位角和所述第二方位角。A digital compass, used to detect the first azimuth and the second azimuth. 13.一种节点编队系统,其特征在于,包括权利要求7-12任一项所述的第一节点、第二节点和第三节点,其中,所述第二节点和第三节点分别具有与所述第一节点相同的所述检测单元、所述计算单元和所述驱动单元。13. A node formation system, characterized in that it comprises the first node, the second node and the third node according to any one of claims 7-12, wherein the second node and the third node respectively have The detecting unit, the calculating unit and the driving unit of the first node are the same. 14.根据权利要求13所述的节点编队系统,其特征在于,所述第一节点、所述第二节点和所述第三节点均为智能机器人。14. The node formation system according to claim 13, wherein the first node, the second node and the third node are all intelligent robots.
CN2013102257773A 2013-06-07 2013-06-07 Device node formation method, device and system Pending CN103336527A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2013102257773A CN103336527A (en) 2013-06-07 2013-06-07 Device node formation method, device and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2013102257773A CN103336527A (en) 2013-06-07 2013-06-07 Device node formation method, device and system

Publications (1)

Publication Number Publication Date
CN103336527A true CN103336527A (en) 2013-10-02

Family

ID=49244718

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2013102257773A Pending CN103336527A (en) 2013-06-07 2013-06-07 Device node formation method, device and system

Country Status (1)

Country Link
CN (1) CN103336527A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106547272A (en) * 2016-10-26 2017-03-29 北京京东尚科信息技术有限公司 The method and apparatus for determining equipment mobile route
CN109917781A (en) * 2017-12-13 2019-06-21 北京京东尚科信息技术有限公司 For dispatching the method, apparatus and system of automated guided vehicle
CN112098945A (en) * 2020-05-13 2020-12-18 苏州触达信息技术有限公司 A method and smart device for finding vehicles
CN115061367A (en) * 2022-07-20 2022-09-16 电子科技大学 A swarm control method based on the Henneberg constraint under azimuth only

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5329450A (en) * 1991-05-10 1994-07-12 Shinko Electric Co., Ltd. Control method for mobile robot system
WO2000041000A1 (en) * 1998-12-30 2000-07-13 L-3 Communications Corporation Close/intra-formation positioning collision avoidance system and method
CN1301001A (en) * 1999-12-23 2001-06-27 李善伯 Method for realizing vehicle queue operation in road traffic
CN101127657A (en) * 2007-07-16 2008-02-20 江南大学 Dynamic Modeling and Control Techniques for Autonomous Mobile Sensor Networks
CN101515179A (en) * 2009-02-17 2009-08-26 浙江大学 Multi- robot order switching method
CN201348824Y (en) * 2008-12-22 2009-11-18 昆明理工大学 Pre-formation none-waiting traffic flow control device
CN102331711A (en) * 2011-08-12 2012-01-25 江苏合成物联网科技有限公司 Formation control method for mobile autonomous robots
CN102662377A (en) * 2012-05-17 2012-09-12 哈尔滨工业大学 Formation system and formation method of multi-mobile robot based on wireless sensor network

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5329450A (en) * 1991-05-10 1994-07-12 Shinko Electric Co., Ltd. Control method for mobile robot system
WO2000041000A1 (en) * 1998-12-30 2000-07-13 L-3 Communications Corporation Close/intra-formation positioning collision avoidance system and method
CN1301001A (en) * 1999-12-23 2001-06-27 李善伯 Method for realizing vehicle queue operation in road traffic
CN101127657A (en) * 2007-07-16 2008-02-20 江南大学 Dynamic Modeling and Control Techniques for Autonomous Mobile Sensor Networks
CN201348824Y (en) * 2008-12-22 2009-11-18 昆明理工大学 Pre-formation none-waiting traffic flow control device
CN101515179A (en) * 2009-02-17 2009-08-26 浙江大学 Multi- robot order switching method
CN102331711A (en) * 2011-08-12 2012-01-25 江苏合成物联网科技有限公司 Formation control method for mobile autonomous robots
CN102662377A (en) * 2012-05-17 2012-09-12 哈尔滨工业大学 Formation system and formation method of multi-mobile robot based on wireless sensor network

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李翔: "群机器人分布式编队控制及建模研究", 《中国博士学位论文全文数据库 信息科技辑》, no. 12, 15 December 2012 (2012-12-15) *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106547272A (en) * 2016-10-26 2017-03-29 北京京东尚科信息技术有限公司 The method and apparatus for determining equipment mobile route
CN106547272B (en) * 2016-10-26 2019-12-03 北京京东尚科信息技术有限公司 The method and apparatus for determining equipment movement routine
CN109917781A (en) * 2017-12-13 2019-06-21 北京京东尚科信息技术有限公司 For dispatching the method, apparatus and system of automated guided vehicle
CN112098945A (en) * 2020-05-13 2020-12-18 苏州触达信息技术有限公司 A method and smart device for finding vehicles
CN112098945B (en) * 2020-05-13 2024-03-22 苏州触达信息技术有限公司 A method and smart device for finding vehicles
CN115061367A (en) * 2022-07-20 2022-09-16 电子科技大学 A swarm control method based on the Henneberg constraint under azimuth only
CN115061367B (en) * 2022-07-20 2025-01-24 电子科技大学 A cluster control method based on Henneberg constraints only in azimuth

Similar Documents

Publication Publication Date Title
Ge et al. A scalable adaptive approach to multi-vehicle formation control with obstacle avoidance
Li et al. Variable time headway policy based platoon control for heterogeneous connected vehicles with external disturbances
CN104865960B (en) A kind of multiple agent approach to formation control based on plane
Robuffo Giordano et al. A passivity-based decentralized strategy for generalized connectivity maintenance
Pan et al. Multi-robot obstacle avoidance based on the improved artificial potential field and PID adaptive tracking control algorithm
CN105138006B (en) A cooperative tracking control method for time-delay nonlinear multi-agent systems
CN104154917B (en) Planning method and device of robot collision prevention path
CN103336527A (en) Device node formation method, device and system
CN102207736A (en) Robot path planning method and apparatus thereof based on Bezier curve
Soltero et al. Collision avoidance for persistent monitoring in multi-robot systems with intersecting trajectories
CN103412564A (en) Unmanned system distributed consistency formation control method and system thereof
CN106483958A (en) A kind of man-machine coordination based on obstacle figure and potential field method is formed into columns and is followed and barrier-avoiding method
CN105353768A (en) Unmanned plane locus planning method based on random sampling in narrow space
CN104181813B (en) There is the Lagrange system self-adaptation control method of connective holding
CN112327839A (en) Formation control method, device, equipment and medium for multi-robot system
Tnunay et al. Distributed collision-free coverage control of mobile robots with consensus-based approach
CN114326810A (en) An obstacle avoidance method for unmanned aerial vehicles in complex dynamic environments
CN114296473A (en) Multi-agent self-adaptive formation control method for avoiding collision and communication interruption
CN114995466B (en) A method and system for generating three-dimensional space-time motion corridors of multiple unmanned vehicles
Gong et al. A virtual spring strategy for cooperative control of connected and automated vehicles at signal-free intersections
CN111007848A (en) Multi-agent cooperative operation control method based on bounded space
Qin et al. Formation control of robotic swarm using bounded artificial forces
Wang et al. Zero Touch Coordinated UAV Network Formation for $360^{\circ} $ Views of a Moving Ground Target in Remote VR Applications
Ma et al. Finite-time flocking control of a swarm of cucker-smale agents with collision avoidance
Tang et al. On non-escape search for a moving target by multiple mobile sensor agents

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: HANSHAN NORMAL COLLEGE

Free format text: FORMER OWNER: HU'NAN UNIVERSITY OF SCIENCE AND ENGINEERING

Effective date: 20140709

C41 Transfer of patent application or patent right or utility model
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Li Xiang

Inventor after: Chen Hongcai

Inventor before: Li Xiang

COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 425101 YONGZHOU, HUNAN PROVINCE TO: 521000 CHAOZHOU, GUANGDONG PROVINCE

Free format text: CORRECT: INVENTOR; FROM: LI XIANG TO: LI XIANG CHEN HONGCAI

TA01 Transfer of patent application right

Effective date of registration: 20140709

Address after: 521000 Qiaodong, Xiangqiao District, Guangdong, Chaozhou

Applicant after: Hanshan Normal University

Address before: Yang Zitang road 425101 Yongzhou city of Hunan province Lingling District No. 130

Applicant before: Hunan University of Science and Engineering

RJ01 Rejection of invention patent application after publication

Application publication date: 20131002

RJ01 Rejection of invention patent application after publication