CN113325843B - 一种基于多智能体覆盖的污染检测方法 - Google Patents

一种基于多智能体覆盖的污染检测方法 Download PDF

Info

Publication number
CN113325843B
CN113325843B CN202110556830.2A CN202110556830A CN113325843B CN 113325843 B CN113325843 B CN 113325843B CN 202110556830 A CN202110556830 A CN 202110556830A CN 113325843 B CN113325843 B CN 113325843B
Authority
CN
China
Prior art keywords
imu
voronoi
agent
wheel robot
coverage
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.)
Active
Application number
CN202110556830.2A
Other languages
English (en)
Other versions
CN113325843A (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.)
Hangdian Haining Information Technology Research Institute Co ltd
Hangzhou Dianzi University
Original Assignee
Hangdian Haining Information Technology Research Institute Co ltd
Hangzhou Dianzi University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangdian Haining Information Technology Research Institute Co ltd, Hangzhou Dianzi University filed Critical Hangdian Haining Information Technology Research Institute Co ltd
Priority to CN202110556830.2A priority Critical patent/CN113325843B/zh
Publication of CN113325843A publication Critical patent/CN113325843A/zh
Application granted granted Critical
Publication of CN113325843B publication Critical patent/CN113325843B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0212Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
    • G05D1/0217Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory in accordance with energy consumption, time reduction or distance reduction criteria

Abstract

本发明涉及一种多智能体的覆盖与污染检测方法。本发明通过一组携带便携式环境综合检测报警仪的IMU独轮机器人最优地覆盖检测目标区域。通过该智能体定义一个虚拟质心,使用一个控制器来确保虚拟质心渐近地与Voronoi分区质心重合,使智能体绕Voronoi分区质心旋转运动,来增加单个智能体动态覆盖的面积,实现对目标区域的高效覆盖,当出现挥发性有机物时,便携式环境综合检测报警仪实现污染报警功能。

Description

一种基于多智能体覆盖的污染检测方法
技术领域
本发明属于智能体技术领域,涉及一种多智能体的覆盖与污染检测方法。
背景技术
污染通常是指由于人类活动或自然过程导致某些有害物质进入大气中危害环境的现象。在污染检测系统中,可分为前期报警、定性检测和精准检测三级,其中前期报警环节是之后两级检测的基本保证,在减小污染损失和制定治理措施方面,有着非常重要的作用。如何对整个区域进行高效的覆盖检测,实现前期报警,是当前技术面临的难题。
近年来,任务环境的覆盖与检测问题,逐渐受到学术界的重视。传统的传感器检测方法是设计无线传感器网络,计算最优覆盖节点,人工放置污染报警传感器。以这种方式实现传感器网络,需要人为地计算最优传感器位置,配置不灵活,放置位置相对固定。并且单个传感器负责检测的区域必须小于检测半径,覆盖区域相对有限,面对面积较大的任务环境时,需要大量的传感器来覆盖整个区域。一旦单个传感器失效,将无法覆盖负责区域,鲁棒性不强。
引入多智能体系统,通过多智能体携带相应的传感器,自主配置并移动到最优的覆盖节点,能解决上述无线传感器网络无法解决的问题。在传统的多智能体覆盖中根据覆盖形状不同主要分为三种覆盖形式:区域覆盖,点覆盖和栅栏覆盖。区域覆盖问题要求多智能体网络中的每个智能体负责一个子区域,目标是减少覆盖检测范围重叠。实现区域覆盖首先需要对区域环境进行区域划分。早期的区域划分方式有:正方形区域,三角形区域,六边形区域三种。
目前智能体区域划分方式的污染覆盖检测存在的缺点:智能体之间缺乏合作与联系,无法实时地更新智能体位置,每个智能体负责的区域无法根据覆盖检测效率最优地划分,智能体覆盖检测范围有大量的重叠,会造成资源的浪费;单个智能体无法动态地检测区域,负责的区域必须小于传感器的检测半径,所以单个智能体覆盖检测的范围有限,覆盖面积较大的区域时,要求的智能体数量极多,导致成本上升。
发明内容
本发明针对上述覆盖污染问题进行了研究,提出了基于Voronoi分割的多智能体污染覆盖检测方法,可以让携带便携式环境综合检测报警仪的IMU独轮机器人自动地配置到基于覆盖检测效率最优的位置;该智能体之间遵守Voronoi分割规则,保持合理的距离,实现合作和联系,提高了报警仪检测覆盖的效率,节省智能体的使用数量。
本发明的技术方案,通过一组携带便携式环境综合检测报警仪的IMU独轮机器人最优地覆盖检测目标区域。通过该智能体定义一个虚拟质心,使用一个控制器来确保虚拟质心渐近地与Voronoi分区质心重合,使智能体绕Voronoi分区质心旋转运动,来增加单个智能体动态覆盖的面积,实现对目标区域的高效覆盖,当出现挥发性有机物时,便携式环境综合检测报警仪实现污染报警功能。
本发明的具体方法是:
步骤一:定义覆盖区域,分割Voronoi分区。
假设Q为需要覆盖污染检测的凸多边形区域。Q中的任意点表示为q,定义一组智能体{1,…,n},智能体的位置由P=(p1,…,pn)表示,智能体k表示智能体中的任意一个智能体,即k∈{1,…,n},每个智能体在凸多边形区域Q中移动。Q中的任意点q由距离成本最小的智能体负责,从而获得根据覆盖检测效率最优的分割。
智能体k的Voronoi分区的定义如下:
Figure BDA0003077492130000021
智能体{1,…,n}将凸多边形区域Q分割为对应的Voronoi分区,即V1,…,Vn,由每个智能体负责覆盖。
步骤二:求得对应步骤一的Voronoi分区质心。
覆盖范围性能函数的梯度定义如下:
Figure BDA0003077492130000022
求得步骤一中Voronoi分区Vk的梯度H的临界点
Figure BDA0003077492130000023
Figure BDA0003077492130000024
是Voronoi分区质心的位置。接下来将IMU独轮机器人作为本专利考虑的智能体,绕Voronoi分区质心覆盖检测。
步骤三:IMU独轮机器人的动力学方程。
根据IMU独轮机器人的动力学方程,结合步骤二得出该IMU独轮机器人的运动方式。
IMU独轮机器人的动力学方程:
Figure BDA0003077492130000025
Figure BDA0003077492130000026
其中
Figure BDA0003077492130000027
表示IMU独轮机器人k在复平面上的位置,vk为该智能体恒定运动速度,θ为航向角。uk为运动角度控制输入。
Figure BDA0003077492130000028
分别为IMU独轮机器人k在复平面上的位置和航向角的时间导数。
使用复数符号,IMU独轮机器人k运动的虚拟质心zk定义为
Figure BDA0003077492130000029
虚拟质心的动力学方程为
Figure BDA00030774921300000210
其中ω0为非零期望角速度。
Figure BDA00030774921300000211
为虚拟质心位置的时间导数。
步骤四:结合步骤二,步骤三,使用控制器,IMU独轮机器人进入各自对应的区域,实现绕Voronoi分区质心做圆周运动。
使用的控制器为
Figure BDA00030774921300000212
IMU独轮机器人逐渐运动到虚拟质心和Voronoi分区质心重合,一旦虚拟质心到达Voronoi分区质心,智能体将以固定的前进速度和角速度(即稳态圆形轨道)绕Voronoi分区质心做圆周运动。
步骤五:通过步骤四的控制器作用,IMU独轮机器人到达预计轨道,开始覆盖污染检测。
IMU独轮机器人绕着Voronoi分区质心做圆周运动,在一个旋转周期结束后,对各自的Voronoi分区完成一次污染检测,持续地动态覆盖的Voronoi分区。如发生污染泄露,出现挥发性有机物,当智能体运动到附近,进入检测范围,便携式环境综合检测报警仪报警,实现功能。
本发明具有的有益效果为:
1、提出了基于Voronoi分割的多智能体污染覆盖检测方法。针对智能体区域分割方法效率不高的问题,通过Voronoi法分割目标凸边型区域,使区域内的每个点由距离成本最小的智能体负责,根据性能函数梯度,自动得到最优覆盖检测区域,大大减少了智能体之间的重叠覆盖区域,减少了资源的浪费。
2、实现了智能体对负责Voronoi分区的动态污染检测。针对单个智能体无法动态地检测区域的问题,通过智能体绕着Voronoi分区质心做圆周运动,实现动态污染检测,大大增加了单个智能体覆盖检测的面积,从而减少了智能体节点数量,降低了节点布置的成本。
附图说明
图1为本发明流程图。
图2为Voronoi分区与对应质心
Figure BDA0003077492130000032
图3为IMU独轮机器人仿真运动轨迹。
图4为质心重合后某一时刻,其中两个智能体实施污染检测。
具体实施方式
步骤一:定义覆盖区域,划分为Voronoi分区。
假设Q为需要覆盖污染检测的凸多边形区域。Q中的任意点表示为q,定义一组多智能体{1,…,n},智能体的位置由P=(p1,…,pn)表示,智能体k表示多智能体中任意一个智能体,即k∈{1,…,n},每个智能体在凸多边形区域Q中移动。Q中的任意点q由距离成本最小的智能体负责,从而获得根据覆盖检测效率最优的划分凸多边形区域。智能体k的Voronoi分区的定义如下:
Figure BDA0003077492130000031
当两个Voronoi分区Vk和Vj相邻(即它们共享一条边)时,pj称为pk的相邻智能体(反之亦然)。根据Voronoi分区的定义,对于在边界上的任何q,都有mink∈{1,…,n}||q-pk||2=||q-pj||2。智能体{1,…,n}将凸多边形区域Q分割为对应的Voronoi分区,即V1,…,Vn,由每个智能体负责覆盖。
步骤二:求得对应步骤一的Voronoi分区质心。
覆盖范围性能函数定义如下
Figure BDA0003077492130000041
其中Φ(q)是所有智能体已知的分布密度函数,并且假设整个Voronoi分区每个点污染出现的概率是固定相等的,那么Φ(q)取1。可以将||q-pk||2看作是对位于智能体k的对q点覆盖程度的定量。
覆盖范围性能函数的梯度定义如下:
Figure BDA0003077492130000042
求得步骤一中Voronoi分区Vk的梯度H的临界点
Figure BDA0003077492130000043
Figure BDA0003077492130000044
是Voronoi分区质心的位置,见图2,圆点为质心
Figure BDA0003077492130000045
虚线为运动轨迹。接下来将IMU独轮机器人作为本发明考虑的智能体,绕Voronoi分区质心覆盖检测。
步骤三:IMU独轮机器人的动力学方程。
根据IMU独轮机器人的动力学方程,结合步骤二得出该IMU独轮机器人的运动方式。
IMU独轮机器人的动力学方程:
Figure BDA0003077492130000046
Figure BDA0003077492130000047
其中
Figure BDA0003077492130000048
表示IMU独轮机器人k在复平面上的位置,vk为该智能体恒定运动速度,q为航向角。uk为运动角度控制输入。
Figure BDA0003077492130000049
分别为IMU独轮机器人k在复平面上的位置和航向角的时间导数。
使用复数符号,IMU独轮机器人k运动的虚拟质心zk定义为
Figure BDA00030774921300000410
虚拟质心的动力学方程为
Figure BDA00030774921300000411
其中ω0为非零期望角速度。
Figure BDA00030774921300000412
为虚拟质心位置的时间导数。如果uk=ω0,则圆轨道的中心等于虚拟质心,
Figure BDA00030774921300000413
虚拟质心的位置不再移动,IMU独轮机器人做半径为Vk/|ω0|的圆周运动。旋转方向由ω0的符号确定。如果ω0>0,则IMU独轮机器人相对于圆周运动的中心沿逆时针方向旋转,否则为顺时针旋转。
步骤四:结合步骤二、三,使用控制器,使IMU独轮机器人进入对应的区域,实现绕Voronoi分区质心做圆周运动。
结合步骤二、三的动力学方程,设计控制器,控制IMU独轮机器人的航向角,使IMU独轮机器人k的虚拟质心逐渐运动到和Voronoi分区质心重合,控制器为
Figure BDA00030774921300000414
一旦虚拟质心到达Voronoi分区质心,即
Figure BDA00030774921300000415
根据控制器公式,将实现uk=ω0,根据步骤三,智能体将以稳态圆形轨道绕Voronoi分区质心运动(固定的前进速度和角速度),运动轨迹如图3所示,三角为起点,点划线为每个机器人运动轨迹,圆点为质心
Figure BDA00030774921300000416
步骤五:通过步骤四的控制器作用,IMU独轮机器人到达预计轨道,开始覆盖污染检测。
IMU独轮机器人绕着质心做圆周运动,便携式环境综合检测报警仪开启,在一个旋转周期结束后,对各自的Voronoi分区完成一次污染检测,持续地动态覆盖各自Voronoi分区。如果发生污染泄露,出现挥发性有机物,当机器人运动到附近,污染源进入检测范围,便携式环境综合检测报警仪将会报警,实现功能。图4为检测工作后某一时刻,其中两个机器人的污染检测情况。方点为当前两个机器人的位置,黑色细实线为检测范围,扩散点1,2为污染点。
使用本发明的多智能体覆盖污染检测系统,能够高效地分割区域,自动得到基于覆盖目标的最优分区,减少重叠覆盖检测范围;能够动态地检测区域,增加单个智能体覆盖检测的面积,减少智能体节点数量,从而降低节点布置成本。

Claims (1)

1.一种基于多智能体覆盖的污染检测方法,其特征在于该方法包括以下步骤:
步骤一:定义覆盖区域,分割Voronoi分区;
假设Q为需要覆盖污染检测的凸多边形区域;Q中的任意点表示为q,定义一组智能体{1,…,n},智能体的位置由P=(p1,…,pn)表示,智能体k表示智能体中的任意一个智能体,即k∈{1,…,n},每个智能体在凸多边形区域Q中移动;Q中的任意点q由距离成本最小的智能体负责,从而获得根据覆盖检测效率最优的分割;
智能体k的Voronoi分区的定义如下:
Figure FDA0003077492120000011
智能体{1,…,n}将凸多边形区域Q分割为对应的Voronoi分区,即V1,…,Vn,由每个智能体负责覆盖;
步骤二:求得对应步骤一的Voronoi分区质心;
覆盖范围性能函数的梯度定义如下:
Figure FDA0003077492120000012
求得步骤一中Voronoi分区Vk的梯度H的临界点
Figure FDA00030774921200000112
Figure FDA00030774921200000111
作为Voronoi分区质心的位置;接下来将IMU独轮机器人作为智能体,绕Voronoi分区质心覆盖检测;
步骤三:IMU独轮机器人的动力学方程;
根据IMU独轮机器人的动力学方程,结合步骤二得出IMU独轮机器人的运动方式;IMU独轮机器人的动力学方程:
Figure FDA0003077492120000013
Figure FDA0003077492120000014
其中
Figure FDA0003077492120000015
表示IMU独轮机器人k在复平面上的位置,vk为该IMU独轮机器人恒定运动速度,θ为航向角;uk为运动角度控制输入;
Figure FDA0003077492120000016
分别为IMU独轮机器人k在复平面上的位置和航向角的时间导数;
使用复数符号,IMU独轮机器人k运动的虚拟质心zk定义为
Figure FDA0003077492120000017
虚拟质心的动力学方程为
Figure FDA0003077492120000018
其中ω0为非零期望角速度;
Figure FDA0003077492120000019
为虚拟质心位置的时间导数;
步骤四:结合步骤二,步骤三,使用控制器,IMU独轮机器人进入各自对应的区域,实现绕Voronoi分区质心做圆周运动;
使用的控制器为
Figure FDA00030774921200000110
IMU独轮机器人逐渐运动到虚拟质心和Voronoi分区质心重合,一旦虚拟质心到达Voronoi分区质心,IMU独轮机器人将以固定的前进速度和角速度绕Voronoi分区质心做圆周运动;
步骤五:通过步骤四的控制器作用,IMU独轮机器人到达预计轨道,开始覆盖污染检测;
IMU独轮机器人绕着Voronoi分区质心做圆周运动,在一个旋转周期结束后,对各自的Voronoi分区完成一次污染检测,持续地动态覆盖的Voronoi分区;如发生污染泄露,出现挥发性有机物,当IMU独轮机器人运动到附近,进入检测范围,便携式环境综合检测报警仪报警,实现功能。
CN202110556830.2A 2021-05-21 2021-05-21 一种基于多智能体覆盖的污染检测方法 Active CN113325843B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110556830.2A CN113325843B (zh) 2021-05-21 2021-05-21 一种基于多智能体覆盖的污染检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110556830.2A CN113325843B (zh) 2021-05-21 2021-05-21 一种基于多智能体覆盖的污染检测方法

Publications (2)

Publication Number Publication Date
CN113325843A CN113325843A (zh) 2021-08-31
CN113325843B true CN113325843B (zh) 2022-03-18

Family

ID=77416212

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110556830.2A Active CN113325843B (zh) 2021-05-21 2021-05-21 一种基于多智能体覆盖的污染检测方法

Country Status (1)

Country Link
CN (1) CN113325843B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114755373B (zh) * 2022-06-16 2023-03-31 西安工业大学 一种基于多机器人编队的空气污染源预警定位方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102427584A (zh) * 2011-12-30 2012-04-25 江南大学 实现多智能体系统针对性覆盖及追踪的方法
CN103529847A (zh) * 2013-10-22 2014-01-22 南京邮电大学 一种基于Voronoi图的多机器人污染控制方法
CN106484955A (zh) * 2016-09-14 2017-03-08 中国科学院计算技术研究所 一种基于Voronoi图的室内布局评价方法及系统
CN108681323A (zh) * 2018-05-14 2018-10-19 浙江工业大学 基于Voronoi图的群组机器人自组装路径规划方法
CN108770079A (zh) * 2018-05-31 2018-11-06 深圳汇通智能化科技有限公司 一种基于水下机器人的水域环境监测系统
CN108981714A (zh) * 2018-08-21 2018-12-11 北航(四川)西部国际创新港科技有限公司 一种无人机路径规划和可行性分析方法
CN110320930A (zh) * 2019-06-17 2019-10-11 中国工程物理研究院电子工程研究所 基于Voronoi图的多无人机编队队形可靠变换方法
CN110703747A (zh) * 2019-10-09 2020-01-17 武汉大学 一种基于简化广义Voronoi图的机器人自主探索方法
CN111750869A (zh) * 2020-08-14 2020-10-09 中国电子科技集团公司第五十四研究所 一种实时重构Voronoi图的无人机路径规划方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102427584A (zh) * 2011-12-30 2012-04-25 江南大学 实现多智能体系统针对性覆盖及追踪的方法
CN103529847A (zh) * 2013-10-22 2014-01-22 南京邮电大学 一种基于Voronoi图的多机器人污染控制方法
CN106484955A (zh) * 2016-09-14 2017-03-08 中国科学院计算技术研究所 一种基于Voronoi图的室内布局评价方法及系统
CN108681323A (zh) * 2018-05-14 2018-10-19 浙江工业大学 基于Voronoi图的群组机器人自组装路径规划方法
CN108770079A (zh) * 2018-05-31 2018-11-06 深圳汇通智能化科技有限公司 一种基于水下机器人的水域环境监测系统
CN108981714A (zh) * 2018-08-21 2018-12-11 北航(四川)西部国际创新港科技有限公司 一种无人机路径规划和可行性分析方法
CN110320930A (zh) * 2019-06-17 2019-10-11 中国工程物理研究院电子工程研究所 基于Voronoi图的多无人机编队队形可靠变换方法
CN110703747A (zh) * 2019-10-09 2020-01-17 武汉大学 一种基于简化广义Voronoi图的机器人自主探索方法
CN111750869A (zh) * 2020-08-14 2020-10-09 中国电子科技集团公司第五十四研究所 一种实时重构Voronoi图的无人机路径规划方法

Also Published As

Publication number Publication date
CN113325843A (zh) 2021-08-31

Similar Documents

Publication Publication Date Title
CN108958289B (zh) 基于相对速度障碍的集群无人机避碰方法
WO2017202141A1 (zh) 基于邻域约束的空地异构机器人系统路径规划方法
Tan et al. Three dimensional collision avoidance for multi unmanned aerial vehicles using velocity obstacle
CN102158988B (zh) 一种有向传感器网络强栅栏覆盖判定方法
CN113325843B (zh) 一种基于多智能体覆盖的污染检测方法
US20210114622A1 (en) Movement control
Kim et al. UAV path planning for maximum visibility of ground targets in an urban area
CN106339004A (zh) 一种植保无人机的作业方法和装置
CN111766784B (zh) 障碍物环境下多机器人图案构成的迭代优化方法
US11662726B2 (en) Controlling movement of a device
CN108212719A (zh) 一种飞机整机自动化喷涂片区划分和机器人站位规划方法
CN115202349A (zh) 基于通信干扰的多移动机器人协同编队控制方法、装置、设备及存储介质
CN107703970B (zh) 无人机集群环绕追踪方法
Lu et al. Fixed-wing UAV path planning in a dynamic environment via dynamic RRT algorithm
CN112596548B (zh) 一种环形任务模式下的多无人机一致性编队控制方法
CN113359831A (zh) 基于任务逻辑调度的集群四旋翼无人机路径生成方法
CN113885565A (zh) 一种多旋翼无人机圆弧转弯的控制方法
Yaguchi et al. A waypoint navigation method with collision avoidance using an artificial potential method on random priority
Han et al. Research on UAV indoor path planning algorithm based on global subdivision grids
CN110216688A (zh) 一种办公区送件服务机器人及其控制方法
Tripathy et al. Guidance of an autonomous agent for coverage applications using range only measurement
Ruangwiset Path generation for ground target tracking of airplane-typed UAV
KR20200062874A (ko) 날개의 각도가 꺾인 6개의 회전날개를 가지는 드론
CN115857543A (zh) 一种基于先验信息的多无人机协同搜索方法
Motonaka et al. 3-dimensional kinodynamic motion planning for an X4-Flyer using 2-dimensional harmonic potential fields

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant