CN109960262B - 一种基于几何法的无人艇动态避障方法和系统 - Google Patents

一种基于几何法的无人艇动态避障方法和系统 Download PDF

Info

Publication number
CN109960262B
CN109960262B CN201910225075.2A CN201910225075A CN109960262B CN 109960262 B CN109960262 B CN 109960262B CN 201910225075 A CN201910225075 A CN 201910225075A CN 109960262 B CN109960262 B CN 109960262B
Authority
CN
China
Prior art keywords
obstacle
unmanned
unmanned ship
obstacle avoidance
point
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
CN201910225075.2A
Other languages
English (en)
Other versions
CN109960262A (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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201910225075.2A priority Critical patent/CN109960262B/zh
Publication of CN109960262A publication Critical patent/CN109960262A/zh
Application granted granted Critical
Publication of CN109960262B publication Critical patent/CN109960262B/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, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/0206Control of position or course in two dimensions specially adapted to water vehicles

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种基于几何法的无人艇动态避障方法和系统,属于路径规划领域。本发明在无人艇感知传感器的检测范围内,计算无人艇和每个障碍物的碰撞风险,对于存在碰撞风险的障碍物,提前进行避障动作,从障碍物运动的后方避开障碍物,保障无人艇的安全,在避障过程中,无人艇基本保持直线航向,在避障点处进行航向调整,缩短无人艇的避障航程,该避障方法对静止障碍物和动态障碍物都具有良好可靠的避障效果。本发明基于几何原理计算每个障碍物的避障点,采用LOS跟踪方法对避障点进行跟踪,间接控制无人艇的航向,降低了路径规划对无人艇运动模型的依赖。

Description

一种基于几何法的无人艇动态避障方法和系统
技术领域
本发明属于路径规划领域,更具体地,涉及一种基于几何法的无人艇动态避障方法和系统。
背景技术
水面无人艇(Unmanned Surface Vessel,USV),简称无人艇,是一种水面智能机器人,能够实现自主航行、自主路径规划以及自主作业,并且能够自主感知环境信息。无人艇无论是在军事上还是在民用上都具有广泛的应用,可以完成一定的海上作战任务,对制定区域的探测、搜索、监视任务,以及对特定目标的围捕等和海岛之间物资投送任务等军事任务,还可以完成河湖水资源的检测、堰塞湖探测和水上救援等民用功能。海上环境复杂,不仅存在暗礁、岛屿、灯塔等静止障碍物,同时也存在其他船只和大型漂浮物等动态障碍物,要求无人艇能够在自主完成航行完成任务,则无人艇必须能够独立自主完成路径规划,因此路径规划是无人艇能够安全自主航行的关键。
路径规划分为基于已知地理信息全局路径规划和基于传感器信息的局部路径规划:基于已知地理信息全局路径规划时,全局路径规划在与任务下达给无人艇的时候,无人艇根据出发点和目标点的位置以及已知的地理信息,规划出一条能够避开已知静止障碍物的路径;基于传感器信息的局部路径规划时,在无人艇执行任务的过程中,为了保证无人艇的安全航行,根据艇上装备的实时环境感知传感器信息,在检测新的动态或静态障碍物时,进行局部路径规划,避开新出现的障碍物。
然而,无人艇与其他机器人相比,具有大惯性、长时延、非线性、模型复杂等运动特性,目前无人艇局部路径规划技术极少结合无人艇的运动特性,实际的避障路径不理想;另外,目前无人艇局部路径规划算法基本采用直接控制无人艇航向的方向,直接控制无人艇航向的方法则需要基于无人艇的运动学模型,对不同运动学模型的无人艇不具备通用性。
发明内容
针对现有技术的缺陷,本发明的目的在于解决现有技术无人艇局部路径规划方法避障路径不理想、不具备通用性的技术问题。
为实现上述目的,第一方面,本发明实施例提供了一种基于几何法的无人艇动态避障方法,该方法包括以下步骤:
S0.获取无人艇的目标点位置信息,初始化避障点集合为{目标点},判断检测到的障碍物个数n是否大于0,若是,初始化t为1,进入步骤S1,否则,直接进入步骤S7;
S1.判断t是否大于n,若是,进入步骤S7,否则,进入步骤S2;
S2.获取无人艇与第t个障碍物的位置信息和速度信息;
S3.根据目标点位置信息、无人艇与第t个障碍物的位置信息和速度信息,计算无人艇与第t个障碍物之间位置关系参数、无人艇和目标点之间位置关系参数;
S4.根据无人艇和目标点之间位置关系参数,判断无人艇是否到达目标点,若是,停止动态避障,否则,进入步骤S5;
S5.根据无人艇与第t个障碍物之间位置关系参数、无人艇与第t个障碍物的速度信息、无人艇与目标点之间位置关系参数,判断无人艇与第t个障碍物是否存在碰撞风险,若是,进入步骤S6,否则,t加1,进入步骤S1;
S6.计算该障碍物对应的避障点,加入避障点集合,t加1,进入步骤S1;
S7.对避障点集合的每个避障点,计算其到无人艇的距离,找出距离无人艇最近的避障点;
S8.计算出无人艇坐标到距离无人艇最近的避障点坐标方向与正东方向的夹角,作为无人艇的航向角。
具体地,所述根据目标点位置信息、无人艇与第t个障碍物的位置信息和速度信息,计算无人艇与第t个障碍物之间位置关系参数、无人艇和目标点之间位置关系参数,具体如下:
无人艇坐标到障碍物圆心坐标方向与正东方向的夹角
Figure BDA0002004907330000031
根据无人艇坐标(xa,ya)与障碍物圆心坐标(xo,yo)采用夹角公式计算得到;
无人艇与障碍物圆心和无人艇与目标点连线的夹角
Figure BDA0002004907330000032
Figure BDA0002004907330000033
Figure BDA0002004907330000034
计算得到,计算公式如下:
Figure BDA0002004907330000035
其中,
Figure BDA0002004907330000036
为无人艇坐标到目标点坐标方向与正东方向的夹角;
Figure BDA0002004907330000037
上一个计算周期
Figure BDA0002004907330000038
的值;
无人艇与障碍物圆心和无人艇与避障点连线的夹角
Figure BDA0002004907330000039
计算公式如下:
Figure BDA00020049073300000310
其中,r0为障碍物安全范围圆的半径,d3为无人艇与障碍物圆心的距离;
无人艇与避障点连线和正东方向的夹角
Figure BDA00020049073300000311
Figure BDA00020049073300000312
Figure BDA00020049073300000313
计算得到:
当选择障碍物安全范围圆右边切点为避障点时,计算公式如下:
Figure BDA00020049073300000314
当选择障碍物安全范围圆左边切点为避障点时,计算公式如下:
Figure BDA0002004907330000041
无人艇与目标点之间的距离d0,根据无人艇坐标(xa,ya)与目标点坐标(xt,yt)采用距离公式计算得到;
障碍物圆心到无人艇与目标点连线的距离d1,由d3
Figure BDA0002004907330000042
计算得到,计算公式如下:
Figure BDA0002004907330000043
d′1:上一个计算周期d1的值;
障碍物到避障点的距离d2,也就是障碍物安全范围圆的半径,d2=r0
无人艇与障碍物圆心的距离d3,可根据无人艇与障碍物圆心坐标采用距离公式计算得到;
无人艇到避障点的距离d4,由d3
Figure BDA0002004907330000044
计算得到,计算公式如下:
Figure BDA0002004907330000045
无人艇在保持当前速度航向情况下,与障碍物之间的出现的最短距离dmin,计算公式如下:
Figure BDA0002004907330000046
Figure BDA0002004907330000047
m1=Xo-Xa
Figure BDA0002004907330000048
m2=yo-ya
Figure BDA0002004907330000049
m1、n1、m2、n2为中间常量。
具体地,圆形障碍物半径加上无人艇与障碍物的安全距离,作为该障碍物的安全范围圆的半径ro
具体地,无人艇与障碍物的安全距离为5L,L为无人艇的长度。
具体地,若d0≤L,即无人艇与目标点的距离小于一倍船长,判断无人艇已到达目标点;若d0>L,则无人艇未达到目标点。
具体地,判断以下三个条件任一个是否满足:(1)障碍物的速度大小vo=0;(2)vo≠0且
Figure BDA0002004907330000051
且vo<va
Figure BDA0002004907330000052
若满足,根据障碍物圆心到无人艇与目标点连线的距离d1
Figure BDA0002004907330000053
判断无人艇与障碍物是否存在碰撞风险,判断方式如下:
Figure BDA0002004907330000054
否则,根据dmin和无人艇与障碍物之间的安全距离ro,判断无人艇与障碍物是否存在碰撞风险,判断方式如下:
Figure BDA0002004907330000055
具体地,判断以下三个条件任一个是否满足:(1)障碍物的速度大小vo=0;(2)vo≠0且
Figure BDA0002004907330000056
且vo<va
Figure BDA0002004907330000057
若满足,避障点的选择依据为:
Figure BDA0002004907330000058
时,选取安全范围圆右边切点为避障点,
Figure BDA0002004907330000059
时,选取安全范围圆左边切点为避障点;
否则,避障点的选择依据为:
当(
Figure BDA00020049073300000510
且d1-d′1<0)或者(
Figure BDA00020049073300000511
Figure BDA00020049073300000512
)时,选择安全范围圆左切点为避障点;
当(
Figure BDA00020049073300000513
且d1-d′1<0)或者(
Figure BDA00020049073300000514
Figure BDA00020049073300000515
)时,选择安全范围圆右切点为避障点。
具体地,选择避障点之后,计算避障点坐标(xb,yb),计算公式如下:
Figure BDA0002004907330000061
第二方面,本发明实施例提供了一种基于几何法的无人艇动态避障系统,所述无人艇动态避障系统采用上述第一方面所述的无人艇动态避障方法。
第三方面,本发明实施例提供了一种计算机可读存储介质,该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时实现上述第一方面所述的无人艇动态避障方法。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
1.本发明基于数学几何的原理计算避障点,采用LOS跟踪方法对避障点进行跟踪,通过跟踪避障点的方法间接控制无人艇的航向,而不是直接控制无人艇航,降低了路径规划对无人艇运动模型的依赖。
2.本发明在无人艇感知传感器的检测范围内,计算碰撞风险,提前进行避障动作,从障碍物运动的后方避开障碍物,保障无人艇的安全,在避障过程中,无人艇基本保持直线航向,在避障点处进行航向调整,缩短无人艇的避障航程。
附图说明
图1为本发明实施例提供的一种基于几何法的无人艇动态避障方法流程图;
图2为本发明实施例提供的障碍物、无人艇、目标点与避障点位置示意图;
图3为本发明实施例提供的避障点在障碍物安全范围圆右边切点示意图;
图4为本发明实施例提供的避障点在障碍物安全范围圆左边切点示意图;
图5为本发明实施例提供的MATLAB平台上对静止障碍物的仿真结果图;
图6为本发明实施例提供的MATLAB平台上对动态障碍物的仿真结果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
夹角公式
夹角
Figure BDA0002004907330000071
表示坐标(x1,y1)到坐标(x2,y2)方向与x轴正方向的夹角,s表示两个坐标点的距离公式。
Figure BDA0002004907330000072
如图1所示,一种基于几何法的无人艇动态避障方法,该方法包括以下步骤:
S0.获取无人艇的目标点位置信息,初始化避障点集合为{目标点},判断检测到的障碍物个数n是否大于0,若是,初始化t为1,进入步骤S1,否则,直接进入步骤S7;
S1.判断t是否大于n,若是,进入步骤S7,否则,进入步骤S2;
S2.获取无人艇与第t个障碍物的位置信息和速度信息;
S3.根据目标点位置信息、无人艇与第t个障碍物的位置信息和速度信息,计算无人艇与第t个障碍物之间位置关系参数、无人艇和目标点之间位置关系参数;
S4.根据无人艇和目标点之间位置关系参数,判断无人艇是否到达目标点,若是,停止动态避障,否则,进入步骤S5;
S5.根据无人艇与第t个障碍物之间位置关系参数、无人艇与第t个障碍物的速度信息、无人艇与目标点之间位置关系参数,判断无人艇与第t个障碍物是否存在碰撞风险,若是,进入步骤S6,否则,t加1,进入步骤S1;
S6.计算该障碍物对应的避障点,加入避障点集合,t加1,进入步骤S1;
S7.对避障点集合的每个避障点,计算其到无人艇的距离,找出距离无人艇最近的避障点;
S8.计算出无人艇坐标到距离无人艇最近的避障点坐标方向与正东方向的夹角,作为无人艇的航向角。
步骤S0.获取无人艇的目标点位置信息,初始化避障点集合为{目标点},判断检测到的障碍物个数n是否大于0,若是,初始化t为1,进入步骤S1,否则,直接进入步骤S7。
目标点是指无人艇执行任务的目的地,在整个动态避障过程是固定不变的。获取无人艇的目标点位置信息,目标点位置坐标表示为(xt,yt)。无人艇在避障过程中,遇到的可能不止存在一个障碍物,因此,需要计算每一个障碍物的避障点,最后会输出一个避障点的集合{xb1,yb1;xb2,yb2…xbn,ybn},n为障碍物的数量。如果n=0,即不存在障碍物。
步骤S1.判断t是否大于n,若是,进入步骤S7,否则,进入步骤S2。
步骤S2.获取无人艇与第t个障碍物的位置信息和速度信息。
如图2所示,障碍物用字母O表示,无人艇用字母A表示,目标点用字母T表示,避障点用字母B表示。
根据差分GPS,获得无人艇的经纬度坐标和速度。将无人艇经纬度坐标信息转化为北东地坐标系(x,y)下的位置信息,正东方向为x轴正方向,正北方向为y轴正方向,选取差分GPS基站为坐标原点,障碍物运动方向角以x轴正方向逆时针旋转为正,范围为0~2π。无人艇位置和速度信息用
Figure BDA0002004907330000094
表示,(xa,ya)表示无人艇的位置坐标,va表示无人艇速度大小,
Figure BDA0002004907330000095
表示无人艇坐标到目标点坐标方向与x轴正方向的夹角,可根据无人艇坐标和目标点坐标采用夹角公式计算得到。
障碍物包括静止障碍物和动态障碍物,静态障碍物如小岛、灯塔和暗礁等,动态障碍物如其他运动船只和其他大型漂浮物等。基于视觉或雷达传感器实时检测障碍物信息,获取在传感器检测范围内的所有障碍物的尺寸、位置以及速度信息。障碍物位置和速度信息用
Figure BDA0002004907330000091
表示,(xo,yo)表示障碍物圆心坐标,vo表示障碍物速度大小,
Figure BDA0002004907330000092
表示障碍物运动方向与x轴正方向的夹角,ro表示障碍物安全范围圆的半径。
根据障碍物尺寸的大小,将障碍物建模为大小不同圆形障碍物。根据无人艇的长度L,设定无人艇与障碍物的安全距离,一般设定为5L。圆形障碍物半径加上无人艇与障碍物的安全距离,作为该障碍物的安全范围圆的半径ro。如果两个障碍物的安全范围圆存在相交或者相切的情况,则这两个障碍物合并为一个尺寸更大的障碍物。
步骤S3.根据目标点位置信息、无人艇与第t个障碍物的位置信息和速度信息,计算无人艇与第t个障碍物之间位置关系参数、无人艇和目标点之间位置关系参数。
输入(xt,yt)、
Figure BDA0002004907330000093
计算以下无人艇与障碍物之间位置关系参数、无人艇和目标点之间位置关系参数。
无人艇坐标到障碍物圆心坐标方向与x轴正方向的夹角
Figure BDA0002004907330000101
可根据无人艇坐标(xa,ya)与障碍物圆心坐标(xo,yo)采用夹角公式计算得到。
无人艇与障碍物圆心和无人艇与目标点连线的夹角
Figure BDA0002004907330000102
Figure BDA0002004907330000103
Figure BDA0002004907330000104
计算得到,计算公式如下:
Figure BDA0002004907330000105
其中,
Figure BDA0002004907330000106
为无人艇坐标到目标点坐标方向与x轴正方向的夹角,
Figure BDA0002004907330000107
为无人艇坐标到障碍物圆心坐标方向与x轴正方向的夹角。
Figure BDA0002004907330000108
具有以下作用:
(1)用于计算障碍物圆心到无人艇与目标点连线的距离d1
(2)根据
Figure BDA0002004907330000109
的正负可以判断障碍物在无人艇与目标点连线的左右。
(3)作为判断障碍物与无人艇是否存在碰撞风险参数之一。
Figure BDA00020049073300001010
上一个计算周期
Figure BDA00020049073300001011
的值。
无人艇与障碍物圆心和无人艇与避障点连线的夹角
Figure BDA00020049073300001012
计算公式如下:
Figure BDA00020049073300001013
其中,r0为障碍物安全范围圆的半径,d3为无人艇与障碍物圆心的距离。
无人艇与避障点连线和x轴正方向的夹角
Figure BDA00020049073300001014
用于计算避障点的坐标,由
Figure BDA00020049073300001015
Figure BDA00020049073300001016
计算得到,计算公式如下:
如图3所示,当选择障碍物安全范围圆右边切点为避障点时:
Figure BDA00020049073300001017
如图4所示,当选择障碍物安全范围圆左边切点为避障点时:
Figure BDA0002004907330000111
无人艇与目标点之间的距离d0,根据无人艇(xa,ya)与目标点的坐标(xt,yt)采用距离公式计算得到,用于判断是否已经到达目标点。
障碍物圆心到无人艇与目标点连线的距离d1,由d3
Figure BDA0002004907330000112
计算得到,计算公式如下:
Figure BDA0002004907330000113
d′1:上一个计算周期d1的值。
障碍物圆到避障点的距离d2,避障点为障碍物安全范围圆切点时d2=r0
无人艇与障碍物圆心的距离d3,可根据无人艇与障碍物圆心坐标采用距离公式计算得到。
无人艇到避障点的距离d4,由d3
Figure BDA0002004907330000114
计算得到,用于计算避障点的坐标,计算公式如下:
Figure BDA0002004907330000115
无人艇在保持当前速度航向情况下,与障碍物之间的出现的最短距离dmin,计算公式如下:
Figure BDA0002004907330000116
Figure BDA0002004907330000117
m1=Xo-Xa
Figure BDA0002004907330000118
m2=yo-ya
Figure BDA0002004907330000119
m1、n1、m2、n2为中间常量。
步骤S4.根据无人艇和目标点之间位置关系参数,判断无人艇是否到达目标点,若是,停止动态避障,否则,进入步骤S5。
如果当d0≤L,即无人艇与目标点的距离小于一倍船长,判断无人艇已到达目标点;如果d0>L,则无人艇没有达到目标点。该判断可有效避免检测到的障碍物位于目标点之后的情形。
步骤S5.根据无人艇与第t个障碍物之间位置关系参数、无人艇与第t个障碍物的速度信息、无人艇与目标点之间位置关系参数,判断无人艇与第t个障碍物是否存在碰撞风险,若是,进入步骤S6,否则,t加1,进入步骤S1。
根据障碍物的位置信息和速度信息,将障碍物分为静止、对遇、被超越和其他动态避障模式四种情况。在判断无人艇与障碍物是否存在碰撞风险时,根据障碍物的运动情况,需要分四种情况进行讨论,分别为障碍物在静态障碍物、“对遇型”动态障碍物、“超越型”动态障碍物和其他动态障碍物四种类型。
(1)静态障碍物
根据环境感知传感器输出的信息,如果障碍物的速度大小vo=0,则该障碍物为静止障碍物。对于静止障碍物,根据障碍物圆心到无人艇与目标点连线的距离d1
Figure BDA0002004907330000121
进行判断无人艇与障碍物是否存在碰撞风险,判断方式如下:
Figure BDA0002004907330000122
(2)“对遇型”动态障碍物
根据环境感知传感器输出的信息,vo≠0且
Figure BDA0002004907330000123
时,无人艇与障碍物为“对遇型”情况。在对遇情况下,根据障碍物圆心到无人艇与目标点连线的距离d1
Figure BDA0002004907330000124
进行判断无人艇与障碍物是否存在碰撞风险,判断方式如下:
Figure BDA0002004907330000131
(3)“超越型”动态障碍物
根据环境感知传感器输出的信息,满足以下情况时,为无人艇超越障碍物情况。当vo≠0且vo<va
Figure BDA0002004907330000132
时,为“超越型”动态障碍物,在无人艇超越障碍物情况下,根据障碍物圆心到无人艇与目标点连线的距离d1
Figure BDA0002004907330000133
进行判断无人艇与障碍物是否存在碰撞风险,判断方式如下:
Figure BDA0002004907330000134
(4)其他动态障碍物
在其他动态障碍物情况下,根据无人艇与障碍物当前的位置信息和速度信息计算得到无人艇和障碍保持目前航行状态会出现的最短距离dmin,根据dmin和无人艇与障碍物之间的安全距离ro可以判断无人艇与障碍物是否存在碰撞风险,判断方式如下:
Figure BDA0002004907330000135
检测到存在碰撞风险之后,通过跟踪避障点开始避障动作。
步骤S6.计算该障碍物对应的避障点,加入避障点集合,t加1,进入步骤S1。
下面为存在碰撞风险的4个类型障碍物避障点的选择依据。
1)静止障碍物、“对遇型”障碍物、“超越型”障碍物的避障点的选择依据为:
Figure BDA0002004907330000136
2)其他动态障碍物避障点的选择依据为:
当(
Figure BDA0002004907330000141
且d1-d′1<0)或者(
Figure BDA0002004907330000142
Figure BDA0002004907330000143
)时,选择安全范围圆左切点避障点;
当(
Figure BDA0002004907330000144
且d1-d′1<0)或者(
Figure BDA0002004907330000145
Figure BDA0002004907330000146
)时,选择安全范围圆右切点避障点。
对于其他动态障碍物,选择的避障点位于障碍物运动的后方。
选择避障点之后,计算避障点坐标(xb,yb),计算公式如下:
Figure BDA0002004907330000147
步骤S7.对避障点集合的每个避障点,计算其到无人艇的距离,找出距离无人艇最近的避障点。
输出的避障点的集合{xb1,yb1;xb2,yb2…xbn,ybn},n为障碍物的数量。采用距离公式计算每个避障点(xbi,ybi)到无人艇(xa,ya)之间的距离{s1,s2…sn},n为障碍物的数量。
选出距离无人艇最近的避障点(xbmin,ybmin),即对应的距离最小。
步骤S8.计算出无人艇坐标到距离无人艇最近的避障点坐标方向与正东方向的夹角,作为无人艇的航向角。
采用夹角公式计算出无人艇坐标(xa,ya)到距离无人艇最近的避障点坐标(xbmin,ybmin)方向与x轴正方向的夹角
Figure BDA0002004907330000148
角作为无人艇的航向角。无人艇就根据这个航向角航行就可以顺利避障。
根据LOS(line-of-sight,视线)导航算法,如果使被控无人艇的航向保持对准视线角(LOS角),那么经过适当的控制,就能使被控船舶达到预期的位置。无人艇基本保持直线航向,在避障点处进行航向调整,缩短无人艇的避障航程。
MATLAB平台上对静止障碍物的仿真结果如图5所示,MATLAB平台上对动态障碍物的仿真结果如图6所示。本算法根据无人艇的运动特性进行实时避障点计算和跟踪,通过仿真实验证明,对静止障碍物和动态障碍物都具有良好可靠的避障效果。
以上,仅为本申请较佳的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应该以权利要求的保护范围为准。

Claims (9)

1.一种基于几何法的无人艇动态避障方法,其特征在于,该方法包括以下步骤:
S0.获取无人艇的目标点位置信息,初始化避障点集合为{目标点},判断检测到的障碍物个数n是否大于0,若是,初始化t为1,进入步骤S1,否则,直接进入步骤S7;
S1.判断t是否大于n,若是,进入步骤S7,否则,进入步骤S2;
S2.获取无人艇与第t个障碍物的位置信息和速度信息;
S3.根据目标点位置信息、无人艇与第t个障碍物的位置信息和速度信息,计算无人艇与第t个障碍物之间位置关系参数、无人艇和目标点之间位置关系参数;
S4.根据无人艇和目标点之间位置关系参数,判断无人艇是否到达目标点,若是,停止动态避障,否则,进入步骤S5;
S5.根据无人艇与第t个障碍物之间位置关系参数、无人艇与第t个障碍物的速度信息、无人艇与目标点之间位置关系参数,判断无人艇与第t个障碍物是否存在碰撞风险,若是,进入步骤S6,否则,t加1,进入步骤S1;
S6.计算该障碍物对应的避障点,加入避障点集合,t加1,进入步骤S1;
S7.对避障点集合的每个避障点,计算其到无人艇的距离,找出距离无人艇最近的避障点;
S8.计算出无人艇坐标到距离无人艇最近的避障点坐标方向与正东方向的夹角,作为无人艇的航向角;
所述根据目标点位置信息、无人艇与第t个障碍物的位置信息和速度信息,计算无人艇与第t个障碍物之间位置关系参数、无人艇和目标点之间位置关系参数,具体如下:
无人艇坐标到障碍物圆心坐标方向与正东方向的夹角
Figure FDA0002345658850000021
根据无人艇坐标(xa,ya)与障碍物圆心坐标(xo,yo)采用夹角公式计算得到;
无人艇与障碍物圆心和无人艇与目标点连线的夹角
Figure FDA0002345658850000022
Figure FDA0002345658850000023
Figure FDA0002345658850000024
计算得到,计算公式如下:
Figure FDA0002345658850000025
其中,
Figure FDA0002345658850000026
为无人艇坐标到目标点坐标方向与正东方向的夹角;
无人艇与障碍物圆心和无人艇与避障点连线的夹角
Figure FDA0002345658850000027
计算公式如下:
Figure FDA0002345658850000028
其中,r0为障碍物安全范围圆的半径,d3为无人艇与障碍物圆心的距离;
无人艇与避障点连线和正东方向的夹角
Figure FDA0002345658850000029
Figure FDA00023456588500000210
Figure FDA00023456588500000211
计算得到:
当选择障碍物安全范围圆右边切点为避障点时,计算公式如下:
Figure FDA00023456588500000212
当选择障碍物安全范围圆左边切点为避障点时,计算公式如下:
Figure FDA0002345658850000031
无人艇与目标点之间的距离d0,根据无人艇坐标(xa,ya)与目标点坐标(xt,yt)采用距离公式计算得到;
障碍物圆心到无人艇与目标点连线的距离d1,由d3
Figure FDA0002345658850000032
计算得到,计算公式如下:
Figure FDA0002345658850000033
障碍物到避障点的距离d2,也就是障碍物安全范围圆的半径,d2=r0
无人艇与障碍物圆心的距离d3,可根据无人艇与障碍物圆心坐标采用距离公式计算得到;
无人艇到避障点的距离d4,由d3
Figure FDA0002345658850000034
计算得到,计算公式如下:
Figure FDA0002345658850000035
无人艇在保持当前速度航向情况下,与障碍物之间的出现的最短距离dmin,计算公式如下:
Figure FDA0002345658850000036
Figure FDA0002345658850000037
m1=xo-xa
Figure FDA0002345658850000038
m2=yo-ya
Figure FDA0002345658850000041
m1、n1、m2、n2为中间常量,vo为障碍物速度大小,va为无人艇速度大小,
Figure FDA0002345658850000042
为障碍物运动方向与X轴正方向的夹角。
2.如权利要求1所述的无人艇动态避障方法,其特征在于,圆形障碍物半径加上无人艇与障碍物的安全距离,作为该障碍物的安全范围圆的半径ro
3.如权利要求2所述的无人艇动态避障方法,其特征在于,无人艇与障碍物的安全距离为5L,L为无人艇的长度。
4.如权利要求1所述的无人艇动态避障方法,其特征在于,若d0≤L,即无人艇与目标点的距离小于一倍船长,判断无人艇已到达目标点;若d0>L,则无人艇未达到目标点。
5.如权利要求1所述的无人艇动态避障方法,其特征在于,判断以下三个条件任一个是否满足:(1)障碍物的速度大小vo=0;(2)vo≠0且
Figure FDA0002345658850000043
(3)vo≠0且vo<va
Figure FDA0002345658850000044
若满足,根据障碍物圆心到无人艇与目标点连线的距离d1
Figure FDA0002345658850000045
判断无人艇与障碍物是否存在碰撞风险,判断方式如下:
Figure FDA0002345658850000046
否则,根据dmin和无人艇与障碍物之间的安全距离ro,判断无人艇与障碍物是否存在碰撞风险,判断方式如下:
Figure FDA0002345658850000047
6.如权利要求1所述的无人艇动态避障方法,其特征在于,判断以下三个条件任一个是否满足:(1)障碍物的速度大小vo=0;(2)vo≠0且
Figure FDA0002345658850000051
(3)vo≠0且vo<va
Figure FDA0002345658850000052
若满足,避障点的选择依据为:
Figure FDA0002345658850000053
时,选取安全范围圆右边切点为避障点,
Figure FDA0002345658850000054
时,选取安全范围圆左边切点为避障点;
否则,避障点的选择依据为:
当(
Figure FDA0002345658850000055
且d1-d′1<0)或者(
Figure FDA0002345658850000056
Figure FDA0002345658850000057
Figure FDA0002345658850000058
)时,选择安全范围圆左切点为避障点,
Figure FDA0002345658850000059
为上一个计算周期
Figure FDA00023456588500000510
的值,d′1为上一个计算周期d1的值;
当(
Figure FDA00023456588500000511
且d1-d′1<0)或者(
Figure FDA00023456588500000512
Figure FDA00023456588500000513
Figure FDA00023456588500000514
)时,选择安全范围圆右切点为避障点。
7.如权利要求6所述的无人艇动态避障方法,其特征在于,选择避障点之后,计算避障点坐标(xb,yb),计算公式如下:
Figure FDA00023456588500000515
8.一种基于几何法的无人艇动态避障系统,其特征在于,所述无人艇动态避障系统采用如权利要求1至7任一项所述的无人艇动态避障方法。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述的无人艇动态避障方法。
CN201910225075.2A 2019-03-25 2019-03-25 一种基于几何法的无人艇动态避障方法和系统 Active CN109960262B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910225075.2A CN109960262B (zh) 2019-03-25 2019-03-25 一种基于几何法的无人艇动态避障方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910225075.2A CN109960262B (zh) 2019-03-25 2019-03-25 一种基于几何法的无人艇动态避障方法和系统

Publications (2)

Publication Number Publication Date
CN109960262A CN109960262A (zh) 2019-07-02
CN109960262B true CN109960262B (zh) 2020-05-19

Family

ID=67024765

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910225075.2A Active CN109960262B (zh) 2019-03-25 2019-03-25 一种基于几何法的无人艇动态避障方法和系统

Country Status (1)

Country Link
CN (1) CN109960262B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110703765B (zh) * 2019-11-07 2022-05-06 大连海事大学 一种无人船的碰撞自规避方法及系统
CN110837255B (zh) * 2019-11-08 2022-07-15 哈尔滨工程大学 一种适用于高速水面无人艇的自主危险规避方法
CN110849370A (zh) * 2019-11-14 2020-02-28 中国船舶重工集团公司第七0七研究所 一种基于水面无人艇的动态航线规划方法
CN110906934B (zh) * 2019-11-29 2022-04-15 华中科技大学 一种基于碰撞危险系数的无人艇避障方法和系统
TWI756647B (zh) * 2020-03-18 2022-03-01 財團法人船舶暨海洋產業研發中心 基於人工勢場法的船舶避碰方法及系統
CN113110460B (zh) * 2021-04-20 2022-09-02 哈尔滨工程大学 一种动态环境下获取水面无人艇艏向可行区间的方法
CN113625709A (zh) * 2021-07-22 2021-11-09 中国舰船研究设计中心 一种水面无人艇避障方法
CN115182408B (zh) * 2022-08-10 2024-03-15 江苏科技大学 一种内河水域清理器避障方法及自主移动方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8849483B2 (en) * 2011-04-13 2014-09-30 California Institute Of Technology Target trailing with safe navigation with colregs for maritime autonomous surface vehicles
CN106845716B (zh) * 2017-01-25 2020-09-08 东南大学 一种基于导航误差约束的水面无人艇局部分层路径规划方法
CN107168335B (zh) * 2017-06-28 2020-06-23 大连海事大学 一种考虑混合多目标避障的水面无人艇路径跟踪制导方法
CN108319291B (zh) * 2017-12-13 2020-12-11 中国人民解放军空军工程大学 一种基于安全边界分析的无人机认知防碰撞控制方法
CN108196451B (zh) * 2017-12-29 2021-01-22 中国电子科技集团公司信息科学研究院 一种仿生鱼群避障行为控制方法
CN108445879B (zh) * 2018-03-12 2021-02-23 上海大学 一种基于预测碰撞危险区域的无人艇避障方法
CN109298708B (zh) * 2018-08-31 2021-08-17 中船重工鹏力(南京)大气海洋信息系统有限公司 一种融合雷达与光电信息的无人艇自主避障方法
CN109358615A (zh) * 2018-09-07 2019-02-19 上海大学 基于传感器检测障碍物速度误差的无人水面艇vo避障方法

Also Published As

Publication number Publication date
CN109960262A (zh) 2019-07-02

Similar Documents

Publication Publication Date Title
CN109960262B (zh) 一种基于几何法的无人艇动态避障方法和系统
Polvara et al. Obstacle avoidance approaches for autonomous navigation of unmanned surface vehicles
CN108445879B (zh) 一种基于预测碰撞危险区域的无人艇避障方法
WO2020253028A1 (zh) 一种基于航迹重规划的水面无人艇动态避碰方法
US10429845B2 (en) System and method for controlling a position of a marine vessel near an object
Peng et al. Development of the USV ‘JingHai-I’and sea trials in the Southern Yellow Sea
Larson et al. Advances in autonomous obstacle avoidance for unmanned surface vehicles
US11600184B2 (en) Marine propulsion control system and method
US11804137B1 (en) Marine propulsion control system and method with collision avoidance override
CN109460045B (zh) 动态障碍在线感知下usv基于改进蚁群优化的避碰规划方法
WO2021230356A1 (ja) 船舶の自動誘導方法、船舶の自動誘導プログラム、船舶の自動誘導システム、及び船舶
CN102999050A (zh) 一种智能水下机器人的自主避障方法
EP3657214B1 (en) Proximity sensing system and method for a marine vessel
JP2019094046A (ja) 物体付近の船舶のポジションを制御するためのシステムおよび方法
US11702178B2 (en) Marine propulsion control system, method, and user interface for marine vessel docking and launch
CN114387824B (zh) 一种符合国际海上避碰规则的避碰转向判定方法
JP6686249B2 (ja) 船舶の自動操縦システム、船舶、及び船舶の自動操縦方法
US11403955B2 (en) Marine propulsion control system and method with proximity-based velocity limiting
Noguchi et al. Tracking omnidirectional surfaces using a low-cost autonomous underwater vehicle
CN114061565A (zh) 一种无人船舶slam及其应用方法
CN114384917B (zh) 一种基于场论的船舶实时避碰方法
Mahipala et al. Model Predictive Control for Path Following and Collision-Avoidance of Autonomous Ships in Inland Waterways
US11794865B1 (en) Proximity sensing system and method for a marine vessel
CN114408122B (zh) 一种船舶防碰撞控制系统的设计方法
Praczyk et al. High–Level Control System for Biomimetic Autonomous Under-water Vehicle

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
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20190702

Assignee: JIUJIANG DIVISION OF NO.707 INSTITUTE, CHINA SHIPBUILDING INDUSTRY Corp.

Assignor: HUAZHONG University OF SCIENCE AND TECHNOLOGY

Contract record no.: X2022420000075

Denomination of invention: A Geometric-based Dynamic Obstacle Avoidance Method and System for Unmanned Vehicles

Granted publication date: 20200519

License type: Common License

Record date: 20220824