CN111338212B - 一种基于极值搜索的欠驱动船舶寻源控制方法 - Google Patents

一种基于极值搜索的欠驱动船舶寻源控制方法 Download PDF

Info

Publication number
CN111338212B
CN111338212B CN202010168726.1A CN202010168726A CN111338212B CN 111338212 B CN111338212 B CN 111338212B CN 202010168726 A CN202010168726 A CN 202010168726A CN 111338212 B CN111338212 B CN 111338212B
Authority
CN
China
Prior art keywords
ship
under
actuated
source
searching
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
CN202010168726.1A
Other languages
English (en)
Other versions
CN111338212A (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN202010168726.1A priority Critical patent/CN111338212B/zh
Publication of CN111338212A publication Critical patent/CN111338212A/zh
Application granted granted Critical
Publication of CN111338212B publication Critical patent/CN111338212B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种基于极值搜索的欠驱动船舶寻源控制方法,包括:建立信号源非线性模型,明确寻找目标—信号源中心的位置表示;建立六自由度欠驱动船舶的控制系统模型,通过船舶线速度控制器控制和滤波器处理,得到设计的损失函数的梯度,通过平均化方法对得到的梯度不断逼近,使得xs和ys逐渐收敛到极值x*和y*,即控制欠驱动船舶(xs,ys)逼近到信号源中心(x*,y*)附近。经过严格数学证明,证实了该方法的有效性,可运用于小型船舶寻源实际控制。

Description

一种基于极值搜索的欠驱动船舶寻源控制方法
技术领域
本发明属于欠驱动船舶的寻源控制技术领域,具体涉及一种基于非模型极值搜索算法的偏航角速度固定欠驱动船舶的寻源控制方法。
背景技术
早在二十世纪初,在车辆寻源控制领域就呈现一种趋势—开发完全自主的车辆,即该车辆不必需要具有观察自身位置的能力。原因有两方面:(1)在水下、冰下或在没有全球定位系统的洞穴等环境的应用需求;(2)高成本的INS系统,难以在较长的时间内保持准确。
现在关于非完整车辆传感器的控制通常有以下几种方法:(1)动态表面控制(DSC)方法,(2)链式系统控制方法,(3)滑模控制方法。然而这些非完整寻源搜索方法计算量大,计算复杂,提高了算法实现成本,而且这些寻源控制在一般陆地需求状况下可以满足要求,然而因为环境约束以及运动体的运动方式不同,类似的非完整寻源搜索便难以满足水环境的寻源要求。
所以一种稳定性好、实时性强、成本低的关于在水环境中欠驱动船舶的寻源控制方法有了实际需求。
发明内容
本发明的目的是克服现有技术的不足,提供一种能应用在类似海洋河流之类的水环境中、计算量小且成本低的基于极值搜索的欠驱动船舶寻源控制方法,例如运用于水上救援,水上清理污染物等作业环境中。
为了解决上述技术问题,本发明采用如下技术方案实现:
一种基于极值搜索的欠驱动船舶寻源控制方法,包括如下步骤:
步骤一、建立信号源非线性模型,明确寻找目标,即未知信号源位置(x*,y*)表示;
步骤二、建立六自由度欠驱动船舶的运动学模型,并简化为浪涌、摇摆、偏航运动控制;
步骤三、设计信号观测器将船舶上传感器采集到的变化的源信号传输到极值搜索算法模块,通过控制器控制和滤波器处理得到损失函数梯度,传感器的位置坐标为(xs,ys);
步骤四、通过平均化方法对得到的梯度不断逼近,使得xs和ys逐渐收敛到极值x*和y*,即欠驱动船舶的传感器坐标(xs,ys)稳定逼近信号源坐标(x*,y*),完成欠驱动船舶寻源追踪的目的。
进一步改进,寻源过程欠驱动船舶沿着信号场最优梯度方向。
进一步改进,确定最优梯度方向的方法为:
利用极值搜索算法对船身传感器位姿与信号源位姿进行处理;
采用平均化方法对得到的梯度不断逼近,同时通过滤波器使输出误差减小到固定较小值。
进一步改进,所述步骤一、建立信号源非线性模型,明确寻找目标,即未知信号源位置(x*,y*)表示:
在寻源控制过程中,由一般规律可知,跟踪信号源分布于一个未知的非线性映射J=f(x,y),该映射具有孤立的局部最大值,即信号源中心f*=f(x*,y*),其表达式为:
J=f(x,y)=f*-qx(x-x*)2-qy(y-y*)2 (1)
其中,(x*,y*)是未知信号源位置坐标,(x,y)为任意时刻欠驱动船舶位置坐标,而qx,qy是未知正数。
进一步改进,所述步骤二、建立六自由度欠驱动船舶的运动学模型,并简化为浪涌、摇摆、偏航运动控制;
建立的欠驱动船舶模型的运动学描述为:
Figure BDA0002408369530000021
Figure BDA0002408369530000022
Figure BDA0002408369530000023
其中,(xc,yc)为船舶中心位置的坐标值,xc是船舶中心位置的横坐标值,yc是船舶中心位置的纵坐标值;
Figure BDA0002408369530000024
是船舶在大地坐标系参考下的偏航角度,且固定初始角度为0,u和v是浪涌线速度和摇摆线速度,r则是偏航角速度;
在距离船舶中心位置(xc,yc)为d的位置配置传感器,则在参考大地坐标系中传感器安装位置的坐标为(xs,ys),可以得到该传感器在控制过程中的位置坐标表示:
Figure BDA0002408369530000025
Figure BDA0002408369530000026
进一步改进,所述步骤三、设计信号观测器将船身传感器采集到的变化的源信号传输到极值搜索算法模块,通过控制器控制和滤波器处理得到损失函数梯度;
损失函数:
Figure BDA0002408369530000031
Figure BDA0002408369530000032
定义输出误差变量:
Figure BDA0002408369530000033
设计浪涌线速度和摇摆线速度控制器:
Figure BDA0002408369530000034
其中,c,α,w,ω0都是输入参数,为正值常数,Δ是滤波后信号,则
Figure BDA0002408369530000035
将传感器收集到的信号源信号送入极值搜索模块,通过滤波器和控制器共同作用,得到损失函数梯度:
Figure BDA0002408369530000036
Figure BDA0002408369530000037
Figure BDA0002408369530000038
其中τ=ωt;
进一步改进,所述步骤四、通过平均化方法对得到的梯度不断逼近,使得xs和ys逐渐收敛到极值x*和y*
设计平均化方法模块如下:
Figure BDA0002408369530000039
Figure BDA00024083695300000310
Figure BDA0002408369530000041
平均化模型(14-16)的平衡点为:
Figure BDA0002408369530000042
由平均化定理得到以下结果:
limτ→∞sup|xs-x*|=O(α+1/ω+r) (17)
limτ→∞sup|ys-y*|=O(α+1/ω+d) (18)
由(17)和(18)可以得到
Figure BDA0002408369530000043
即欠驱动船舶的传感器坐标(xs,ys)稳定逼近信号源坐标(x*,y*),完成欠驱动船舶寻源追踪的目的。
本发明相对于现有控制技术具有以下的优点和效果:
本发明提出的基于极值搜索的控制算法着眼于水环境,克服了船舶运动学与普通非完整机器人运动学的差异,创新性地改进了现有陆地非完整机器人寻源控制方法。选择的参数较为简单、物理含义清晰,误差小,实际运用中成本低;运动过程中船舶运动稳定性好。
附图说明
图1是六自由度的欠驱动船舶坐标系;
图2是寻源控制算法控制框图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述:
如图1、2所示,一种基于极值搜索的欠驱动船舶寻源控制方法,包含以下步骤:
1)、建立欠驱动船舶的控制系统模型:
一般的六自由度欠驱动船舶模型可以简化为浪涌、摇摆和偏航运动控制,则建立的欠驱动船舶模型的运动学描述为:
Figure BDA0002408369530000044
Figure BDA0002408369530000045
Figure BDA0002408369530000046
其中,(xc,yc)为船舶中心位置的坐标值,xc是船舶中心位置的横坐标值,yc是船舶中心位置的纵坐标值;
Figure BDA0002408369530000051
是船舶在大地坐标系参考下的偏航角度,且固定初始角度为0,u和v是浪涌线速度和摇摆线速度,r则是偏航角速度。
在距离船舶中心位置(xc,yc)为d的位置配置传感器,则在参考大地坐标系中传感器安装位置的坐标为(xs,ys),可以得到该传感器在控制过程中的位置坐标表示:
Figure BDA0002408369530000052
Figure BDA0002408369530000053
在寻源控制过程中,由一般规律可知,跟踪信号源分布于一个未知的非线性映射J=f(x,y),该映射具有孤立的局部最大值,即信号源中心f*=f(x*,y*),其表达式为:
J=f(x,y)=f*-qx(x-x*)2-qy(y-y*)2 (6)
其中,(x*,y*)是未知信号源位置坐标,(x,y)为任意时刻欠驱动船舶位置坐标,而qx,qy是未知正数。
2)、设计极值搜索算法模块,包括线速度控制器以及滤波器:
针对上述步骤一中的公式(1-3,6),给出损失函数:
Figure BDA0002408369530000054
Figure BDA0002408369530000055
定义输出误差变量:
Figure BDA0002408369530000056
设计浪涌线速度和摇摆线速度控制器:
Figure BDA0002408369530000057
其中,c,α,w,ω0都是输入参数,为正值常数,Δ是滤波后信号,则
Figure BDA0002408369530000058
将传感器收集到的信号源信号送入极值搜索模块,通过滤波器和控制器共同作用,得到损失函数梯度:
Figure BDA0002408369530000061
Figure BDA0002408369530000062
Figure BDA0002408369530000063
其中τ=ωt;
3)、设计平均化方法模块,对得到的梯度不断逼近,使得xs和ys逐渐收敛到极值x*和y*,即控制欠驱动船舶(xs,ys)逼近到信号源中心(x*,y*)附近,达成寻源目标;设计平均化方法模块如下:
Figure BDA0002408369530000064
Figure BDA0002408369530000065
Figure BDA0002408369530000066
平均化模型(14-16)的平衡点为:
Figure BDA0002408369530000067
由平均化定理得到以下结果:
limτ→∞sup|xs-x*|=O(α+1/ω+r) (17)
limτ→∞sup|ys-y*|=O(α+1/ω+d) (18)
由(17)和(18)可以得到
Figure BDA0002408369530000068
即欠驱动船舶的传感器坐标(xs,ys)稳定逼近信号源坐标(x*,y*),完成欠驱动船舶寻源追踪的目的。该控制算法着眼于水环境,克服了船舶运动学与普通非完整机器人运动学的差异,创新性地改进了现有陆地非完整机器人寻源控制方法。选择的参数较为简单、物理含义清晰,误差小,实际运用中成本低;运动过程中船舶运动稳定性好。
本发明中涉及的未说明部份与现有技术相同或采用现有技术加以实现。

Claims (1)

1.一种基于极值搜索的欠驱动船舶寻源控制方法,其特征在于,包括如下步骤:
步骤一、建立信号源非线性模型,明确寻找目标,即未知信号源位置(x*,y*)表示;在寻源控制过程中,由一般规律可知,跟踪信号源分布于一个未知的非线性映射J=f(x,y),该映射具有孤立的局部最大值,即信号源中心f*=f(x*,y*),其表达式为:
J=f(x,y)=f*-qx(x-x*)2-qy(y-y*)2 (1)
其中,(x*,y*)是未知信号源位置坐标,(x,y)为任意时刻欠驱动船舶位置坐标,而qx,qy是未知正数;
步骤二、建立六自由度欠驱动船舶的运动学模型,并简化为浪涌、摇摆、偏航运动控制;建立的欠驱动船舶模型的运动学描述为:
Figure FDA0003595361840000011
Figure FDA0003595361840000012
Figure FDA0003595361840000013
其中,(xc,yc)为船舶中心位置的坐标值,xc是船舶中心位置的横坐标值,yc是船舶中心位置的纵坐标值;
Figure FDA0003595361840000014
是船舶在大地坐标系参考下的偏航角度,且固定初始角度为0,u和v是浪涌线速度和摇摆线速度,r则是偏航角速度;
在距离船舶中心位置(xc,yc)为d的位置配置传感器,则在参考大地坐标系中传感器安装位置的坐标为(xs,ys),可以得到该传感器在控制过程中的位置坐标表示:
Figure FDA0003595361840000015
Figure FDA0003595361840000016
步骤三、设计信号观测器将船舶上传感器采集到的变化的源信号传输到极值搜索算法模块,通过控制器控制和滤波器处理得到损失函数梯度,传感器的位置坐标为(xs,ys);损失函数:
Figure FDA0003595361840000017
Figure FDA0003595361840000018
定义输出误差变量:
Figure FDA0003595361840000019
设计浪涌线速度和摇摆线速度控制器:
Figure FDA00035953618400000110
其中,c,α,w,ω0都是输入参数,为正值常数,Δ是滤波后信号,则
Figure FDA0003595361840000021
将传感器收集到的信号源信号送入极值搜索模块,通过滤波器和控制器共同作用,得到损失函数梯度:
Figure FDA0003595361840000022
Figure FDA0003595361840000023
Figure FDA0003595361840000024
其中τ=ωt;
步骤四、通过平均化方法对得到的梯度不断逼近,使得xs和ys逐渐收敛到极值x*和y*,即欠驱动船舶的传感器坐标(xs,ys)稳定逼近信号源坐标(x*,y*),完成欠驱动船舶寻源追踪的目的;寻源过程欠驱动船舶沿着信号场最优梯度方向,确定最优梯度方向的方法为:
利用极值搜索算法对船身传感器位姿与信号源位姿进行处理;
采用平均化方法对得到的梯度不断逼近,同时通过滤波器使输出误差减小到固定较小值;设计平均化方法模块如下:
Figure FDA0003595361840000025
Figure FDA0003595361840000026
Figure FDA0003595361840000027
平均化模型(14-16)的平衡点为:
Figure FDA0003595361840000028
由平均化定理得到以下结果:
limτ→∞sup|xs-x*|=O(α+1/ω+r) (17)
limτ→∞sup|ys-y*|=O(α+1/ω+d) (18)
由(17)和(18)可以得到
Figure FDA0003595361840000031
即欠驱动船舶的传感器坐标(xs,ys)稳定逼近信号源坐标(x*,y*),完成欠驱动船舶寻源追踪的目的。
CN202010168726.1A 2020-03-12 2020-03-12 一种基于极值搜索的欠驱动船舶寻源控制方法 Active CN111338212B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010168726.1A CN111338212B (zh) 2020-03-12 2020-03-12 一种基于极值搜索的欠驱动船舶寻源控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010168726.1A CN111338212B (zh) 2020-03-12 2020-03-12 一种基于极值搜索的欠驱动船舶寻源控制方法

Publications (2)

Publication Number Publication Date
CN111338212A CN111338212A (zh) 2020-06-26
CN111338212B true CN111338212B (zh) 2022-06-14

Family

ID=71184297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010168726.1A Active CN111338212B (zh) 2020-03-12 2020-03-12 一种基于极值搜索的欠驱动船舶寻源控制方法

Country Status (1)

Country Link
CN (1) CN111338212B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102681548A (zh) * 2012-05-16 2012-09-19 东南大学 一种用于轮式机器人搜索核辐射源的运动控制方式
CN109116856A (zh) * 2018-09-28 2019-01-01 上海海事大学 一种基于扰动观测器的欠驱动船舶路径跟踪控制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102681548A (zh) * 2012-05-16 2012-09-19 东南大学 一种用于轮式机器人搜索核辐射源的运动控制方式
CN109116856A (zh) * 2018-09-28 2019-01-01 上海海事大学 一种基于扰动观测器的欠驱动船舶路径跟踪控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
田玉平 等.一种改进的基于分数幂的极值搜索控制方法.《第三十二届中国控制会议论文集(D卷)》 》.2013,第5433-5437页. *
罗冰.时变气流环境下飞行机器人搜寻气味源关键问题研究.《中国博士学位论文全文数据库 工程科技Ⅱ辑》.2019,(第12期),C031-1. *

Also Published As

Publication number Publication date
CN111338212A (zh) 2020-06-26

Similar Documents

Publication Publication Date Title
CN109032178B (zh) 全驱动auv回收控制系统及自主回收方法
CN107168312B (zh) 一种补偿uuv运动学和动力学干扰的空间轨迹跟踪控制方法
CN107957727B (zh) 水下机器人控制系统及动力定位方法
CN110032197B (zh) 一种基于有限时间制导和控制的无人船路径跟踪控制方法
CN108444478B (zh) 一种用于水下航行器的移动目标视觉位姿估计方法
CN111487966B (zh) 一种基于航路点的水面无人艇自适应路径跟踪控制方法
CN108267953B (zh) 一种基于领航者-跟踪者水下机器人位置跟踪方法
CN112036099B (zh) 一种基于帆船风帆攻角确定最小回转半径的方法
CN109579850B (zh) 基于对水速度辅助惯导的深水智能导航方法
CN112710304B (zh) 一种基于自适应滤波的水下自主航行器导航方法
CN112327835A (zh) 一种无人艇滑模编队控制系统及其方法
CN109946976B (zh) 一种宽航速auv运动控制方法
CN113093742B (zh) 一种自动规避多障碍物的无人船路径跟踪系统
CN110926466A (zh) 一种面向无人船组合导航信息融合的多尺度数据分块算法
CN113467231A (zh) 基于侧滑补偿ilos制导律的无人艇路径跟踪方法
CN112015086A (zh) 一种欠驱动水面船有限时间路径跟踪输出反馈控制方法
CN114061565B (zh) 一种无人船舶slam及其应用方法
CN107315348A (zh) 一种基于惩罚式小波网络的无人潜航器位姿控制装置及方法
CN111338212B (zh) 一种基于极值搜索的欠驱动船舶寻源控制方法
CN110456786B (zh) 基于地形特征提取的uuv海底地形跟踪系统及跟踪方法
Kim et al. Development of imaging sonar based autonomous trajectory backtracking using AUVs
CN115031726A (zh) 一种数据融合导航定位方法
CN112241170A (zh) 基于并联六自由度平台的无人艇自稳定系统
CN113110458A (zh) 一种无人船虚拟目标跟踪控制系统
CN110057383A (zh) 一种auv推位导航系统杆臂误差标校方法

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