CN110888104B - 一种接近信标轨迹下的水下机器人定位方法 - Google Patents
一种接近信标轨迹下的水下机器人定位方法 Download PDFInfo
- Publication number
- CN110888104B CN110888104B CN201911064151.2A CN201911064151A CN110888104B CN 110888104 B CN110888104 B CN 110888104B CN 201911064151 A CN201911064151 A CN 201911064151A CN 110888104 B CN110888104 B CN 110888104B
- Authority
- CN
- China
- Prior art keywords
- distance
- underwater robot
- displacement
- observation
- observation distance
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000006073 displacement reaction Methods 0.000 claims abstract description 66
- 238000012549 training Methods 0.000 claims abstract description 41
- 230000008569 process Effects 0.000 claims abstract description 26
- 230000002787 reinforcement Effects 0.000 claims abstract description 12
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 230000009471 action Effects 0.000 claims description 44
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000000875 corresponding effect Effects 0.000 claims description 12
- 230000033001 locomotion Effects 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 3
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 claims description 2
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 claims description 2
- 230000008901 benefit Effects 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 4
- 230000033228 biological regulation Effects 0.000 abstract description 2
- 238000009472 formulation Methods 0.000 abstract description 2
- 239000000203 mixture Substances 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
- G01S1/72—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using ultrasonic, sonic or infrasonic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/14—Systems for determining distance or velocity not using reflection or reradiation using ultrasonic, sonic, or infrasonic waves
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Manipulator (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及水下机器人定位技术,旨在提供一种接近信标轨迹下的水下机器人定位方法。包括:在保持水下机器人匀速率运动且速率大小已知的情况下,按照设定测量周期计算得到当前位移方向与当前观测距离方向之间的夹角;然后通过强化学习和训练,获得艏向角的调整策略并用于调整艏向角,使水下机器人朝向接近信标的方向运动;在接近信标的过程中,利用扩展卡尔曼滤波位置估计方程计算水下机器人的位置以实现定位。本发明在初始估计误差较大的情况下,相比于圆形轨迹,收敛速度更快。用强化学习的方法来调整艏向角,使水下机器人做接近信标的运动。避免了复杂的调整角度规则的制定。不需要准确的初始位置估计值,运算量小,计算简单,定位稳定可靠。
Description
技术领域
本发明涉及水下机器人定位技术,特别涉及一种接近信标轨迹下的水下机器人定位方法。
背景技术
水下机器人在水下执行一段时间任务后,累积误差会随着时间增大。这时候,一般需要浮出水面,接收GPS信号重新定位。对于深水作业的水下机器人,浮出水面需要消耗大量能量,且隐蔽性不好。
单信标测距定位技术由于安装简单、成本低,成为近几年水下定位技术新的发展方向。但是单个测距信息不能完全确定水下机器人的位置,因此需要设计一定的轨迹和相应的算法来实现定位。圆形轨迹和扩展卡尔曼滤波算法是常用的机动轨迹和滤波算法。其中,扩展卡尔曼滤波是一个非线性的滤波器,对初始估计值要求较高,当滤波初始估计值误差较大时,往往会导致收敛速度慢,甚至发散。在初始估计误差较大时,采用圆形轨迹,收敛速度慢。
发明内容
本发明要解决的技术问题是,克服现有技术中的不足,提供一种接近信标轨迹下的水下机器人定位方法。该方法是在单信标测距***下,通过调整艏向角,使水下机器人向信标运动,同时在该轨迹中用扩展卡尔曼滤波进行位置估计,实现定位。
为解决技术问题,本发明的解决方案是:
提供一种接近信标轨迹下的水下机器人定位方法,包括以下步骤:
(1)在保持水下机器人匀速率运动且速率大小已知的情况下,按照设定测量周期,以单个声信标测量水下机器人与信标之间的距离值,以水下机器人自带的罗盘测量艏向角,由前一个时刻的位移距离和前一时刻的观测距离,和当前时刻的观测距离,通过已知三边长度的余弦定理计算得到当前位移方向与当前观测距离方向之间的夹角;然后通过强化学习和训练,获得艏向角的调整策略;
(物理定义中匀速运动是指速度大小和方向不变,本发明所述匀速率是指速度大小不变,但通过调整艏向角的角度改变方向。)
(观测距离方向并不可测量。但通过三边长度已知的余弦定理,可以求出当前位移方向与观测距离方向的夹角。观测距离即由信标测量得到的水下机器人与信标之间的距离值,观测距离也就是量测距离,由于水声测距带有噪声,会与实际距离存在误差。)
(2)利用强化学习训练得到的策略来调整艏向角,使水下机器人朝向接近信标的方向运动;
(3)在接近信标的过程中,利用扩展卡尔曼滤波位置估计方程计算水下机器人的位置以实现定位。
本发明中,所述步骤(1)具体包括:
首先,建立8行×16列的R表格,表格内容是对应于各情况下的平均奖励值,表示选择相应动作时对接近信标的有利程度;其中,
行名包括下述八种情况:情况1:观测距离变大,变化量大于一半的位移距离且位移方向与观测距离方向的夹角变大;情况2:观测距离变大,变化量小于等于一半的位移距离,且位移方向与观测距离方向的夹角变大;情况3:观测距离变小,变化量小于一半的位移距离,且位移距离与观测距离方向的夹角变大;情况4:观测距离变小,变化量大于等于一半的位移距离,且位移距离与观测距离方向的夹角变大;情况5:观测距离变大,变化量大于一半的位移距离,且位移方向与观测距离方向的夹角变小;情况6:观测距离变大,变化量小于等于一半的位移距离,且位移方向与观测距离方向的夹角变小;情况7:观测距离变小,变化量小于一半的位移距离,且位移距离与观测距离方向的夹角变小;情况8:观测距离变小,变化量大于等于一半的位移距离,且位移距离与观测距离方向的夹角变小;
列名包括下述十六个调整动作:顺时针旋转15°、30°、45°、60°、75°、90°、105°、120°和逆时针旋转15°、30°、45°、60°、75°、90°、105°、120°;
其次,以经过强化学习训练的R表格作为艏向角的调整策略,具体训练的过程包括:
(1.1)初始化R表格,确定训练次数;
(1.2)设置水下机器人的起始位置后,使其保持匀速率运动;在第一个时刻里,水下机器人朝任意方向移动;
(1.3)按设定时间间隔,测量当前时刻的观测距离r(t)、上一时刻的观测距离r(t-1)、上一时刻的位移距离d(t-1)和用已知三边长度下的余弦定理,计算位移方向与当前观测距离方向之间的夹角α;如果r(t)<r(t-1)且夹角α小于45°,保持直行;否则将水下机器人的行进方向逆时针旋转5°后直行,并测量获得下一个时刻的观测距离、位移和夹角数据;
(1.4)由连续三个时刻的观测距离和两次位移距离的信息,计算出后两次观测距离的差值r(t+1)-r(t)和两次夹角的差值α2-α1,并在细分后进行归类;
(1.5)以ε-greedy的探索策略选取设计好的十六个动作;执行动作后,根据观测距离变化情况获得对应奖励值Re;当观测距离减小时Re=+1,当观测距离不变或变大时Re=-1;并更新R表格;更新公式为:
R(s,a)=(R(s,a)×(N(s,a)-1)+Re)/N(s,a)
s为情况s1~s8中的一个,a为动作a1~a16,具体为得到奖励值Re前,所归类的情况s,和所执行的动作a;
N(s,a)为训练过程中N(S,A)表格中的一项,表示在情况s下执行了动作a的次数;N(S,A)是一个8×16的表格;更新规则为N(s,a)=N(s,a)+1;
R(s,a)为训练过程中R(S,A)表格中的一项,表示在情况s状态下执行动作a获得奖励值的平均值,R(S,A)是一个8×16的表格;
(1.6)当观测距离小于设定的距离值D;本轮训练结束,返回步骤(1.3);当训练达到设定次数后结束训练,保留训练结果R表格。
(ε-greedy的探索策略是现有技术,常用于强化学习的一种探索策略,以ε概率随机选取一个动作(包括使得R值最大的动作),以1-ε的概率选择R值最大的动作,这样在训练过程中,相比于随机策略,R值最大的动作被更多的选取,被更充分的训练。)
本发明中,所述步骤(2)包括:
(2.1)按设定时间间隔,测量当前时刻的观测距离r(t+1)、上一时刻的观测距离r(t)、上一时刻的位移距离d(t)和用已知三边长度下的余弦定理,计算位移方向与当前观测距离方向之间的夹角α2;如果r(t+1)<r(t)且夹角α2小于45°,保持直行;否则将水下机器人的行进方向逆时针旋转5°后直行,并测量获得下一个时刻的观测距离、位移和夹角数据;
(2.2)调整艏向角后,继续直行,并重新进行测量和计算;重复以上过程,直到水下机器人到信标的距离小于设定的距离值D;然后围绕信标方向随机地运动一段设定时间T,定位结束。
本发明中,所述步骤(3)包括:
(1)设定扩展卡尔曼滤波位置估计方程为:
P(k+1|k)=P(k|k)+Q
K(k+1)=P(k+1|k)·HT(k+1)·[H(k+1)·P(k+1|k)·HT(k+1)+R]-1
P(k+1|k+1)=[I-K(k+1)·H(k+1)]·P(k+1|k)
上式中,k表示k时刻,X(k)=[xk yk]T上标“T”表示向量或者矩阵的转置;xk,yk分别是水下机器人在k时刻以信标为原点的东向和北向坐标;字母上带“^”符号的表示状态量的预测值或估计值;
U(k)=[v·sinθk v·cosθk]T,θk为k时刻的水下机器人的艏向角,v为水下机器人的航速;
I为单位矩阵;P(k|k)是k时刻的误差的协方差矩阵;P(k+1|k)是k+1时刻预测的误差协方差矩阵;P(k+1|k+1)是量测更新后的协方差矩阵;
V(k+1)为k+1时刻过程噪声,W(k+1)为k+1时刻观测噪声,满足零均值高斯分布,Q是过程噪声的方差,R是观测噪声的方差;
Z(k+1)=h(X(k+1))+W(k+1),是即k+1时刻观测到的带有噪声的水下机器与信标之间的距离值;
与现有技术相比,本发明的有益效果是:
1.本发明在初始估计误差较大的情况下,相比于圆形轨迹,收敛速度更快。
2.本发明用强化学习的方法来调整艏向角,使水下机器人做接近信标的运动。避免了复杂的调整角度规则的制定。
3.本发明不需要准确的初始位置估计值,运算量小,计算简单,定位稳定可靠。
附图说明
图1是本发明中信标测距和计算夹角的示意图。
图2是用于八种情况分类的测量示意图。
图3是R表格训练的流程图。
具体实施方式
本发明利用单信标测距导航,在水下机器人初始位置未知的情况下,通过使水下机器人接近信标运动,并用扩展卡尔曼滤波进行位置估计实现定位的。接近信标轨迹下的扩展卡尔曼滤波具有较快的收敛速度,且对初始估计误差大小不敏感。在初始位置未知的情况下,做接近信标的运动需要一定的艏向角调整策略。本发明通过强化学习的方法训练得到艏向角调整策略的结果。
本发明的控制策略以R表格形式存储,通过查询表格选择某一情况下较为合适的角度进行调整。
R表格的结构示例如下:
R(S,A) | a1 | a2 | a3 | … | a16 |
s1 | R(s1,a1) | R(s1,a2) | R(s1,a3) | … | R(s1,a16) |
s2 | R(s2,a1) | R(s2,a2) | R(s2,a3) | … | R(s2,a16) |
s3 | R(s3,a1) | R(s3,a2) | R(s3,a3) | … | R(s3,a16) |
s4 | R(s4,a1) | R(s4,a2) | R(s4,a3) | … | R(s4,a16) |
s5 | R(s5,a1) | R(s5,a2) | R(s5,a3) | … | R(s5,a16) |
s6 | R(s6,a1) | R(s6,a2) | R(s6,a3) | … | R(s6,a16) |
s7 | R(s7,a1) | R(s7,a2) | R(s7,a3) | … | R(s7,a16) |
s8 | R(s8,a1) | R(s8,a2) | R(s8,a3) | … | R(s8,a16) |
表格中,S为情况集合S=<s1,s2,s3,s4....s8>,A为动作集合A=<a1,a2,a3...a16>。
R表格的训练过程示例如下:
例如假设某次得到的R值训练表格如上所示(为清楚展示内容,只显示了表格的一部分),有标记的是该行值最大的一格。
那么当训练完成后,根据表格挑选动作,选择值最大的所对应的动作。
当情况归为s1时候,选择的动作就是a8,情况归为s2的时候动作就是a1。
如果是训练的时候,用ε-greedy随机探索策略时候,有标记的格子所对应的动作被选取的概率大于其他动作,更新R(S,A)表格还需要用的N(S,A)表格(为清楚展示内容,只显示了表格的一部分)。
N(S,A) | a1 | a2 | a3 | a4 | a5 | a6 | a7 | a8 | a9 |
s1 | 7.00 | 57.00 | 73.00 | 89.00 | 58.00 | 67.00 | 211.00 | 8615.00 | 51.00 |
s2 | 2955.00 | 58.00 | 51.00 | 41.00 | 34.00 | 32.00 | 39.00 | 42.00 | 232.00 |
s3 | 168.00 | 34.00 | 48.00 | 59.00 | 74.00 | 28.00 | 20.00 | 32.00 | 2954.00 |
s4 | 342.00 | 7.00 | 17.00 | 2.00 | 9.00 | 9.00 | 3.00 | 1.00 | 3.00 |
s5 | 13.00 | 10.00 | 10.00 | 11.00 | 21.00 | 7.00 | 1397.00 | 61.00 | 12.00 |
s6 | 34.00 | 26.00 | 16.00 | 88.00 | 30.00 | 45.00 | 3448.00 | 27.00 | 21.00 |
s7 | 50.00 | 123.00 | 55.00 | 7468.00 | 78.00 | 72.00 | 54.00 | 73.00 | 62.00 |
s8 | 365.00 | 536.00 | 307.00 | 290.00 | 298.00 | 281.00 | 290.00 | 281.00 | 41265.00 |
比如说情况s1下,执行动作a1后,观测距离减小了,则Re=+1。
更新N(s1,a1)=N(s1,a1)+1=8
更新R(s1,a1)=(R(s1,a1)*(N(s1,a1)-1)+Re)/N(s1,a1)=(-0.22*(8-1)+1)/8=-0.0675。
如图2所示,图中只有7个测量的量(即3+2+2):三个时刻的观测距离值(无方向)、这三个时刻间的两个位移距离,以及这三个时刻间的两个位移方向。按照距离变化(变大变小,变大程度与位移/2之间的大小,变小程度与位移/2大小可分为4种),以及位移与距离夹角的角度(变大与变小)变化来划分出的4x2=8种情况,夹角α是通过三边测距余弦定理测得的。
图2中,后两次观测距离的差值tenzhen1与两次夹角的差值tenzhen2是两个用于情况划分的标志量,α1是前一时刻位移距离d(t-1)与观测距离r(t-1)的夹角,通过d(t-1),r(t-1),r(t)三个距离长度,通过余弦定理计算得到,α2是旋转5°后,位移距离d(t)与观测距离r(t)之间的夹角,d是位移距离值。水下机器人做匀速率运动,则d=d(t)=d(t-1)。tezhen1=r(t+1)-r(t),tezhen2=α2-α1;情况分类的依据是,根据tezhen1大于等于d/2、小于d/2且大于等于0、小于等于0且大于-d/2、小于等于-d/2将所有情况划分为四种情况。再根据tezhen2是否大于0,对这四种情况继续细分,这样所有情况一共被分为八种。为了更好的描述八种情况的分类依据,用以下情况分类表示意。
情况分类表
情况集合S | 分类依据 |
s1 | r(t+1)>r(t)+d(t)/2,α2≥α1 |
s2 | r(t)<r(t+1)≤r(t)+d(t)/2,α2≥α1 |
s3 | r(t)-d(t)/2<r(t+1)≤r(t),α2≥α1 |
s4 | r(t+1)≤r(t)-d(t)/2,α2≥α1 |
s5 | r(t+1)>r(t)+d(t)/2,α2<α1 |
s6 | r(t)<r(t+1)≤r(t)+d(t)/2,α2<α1 |
S7 | r(t)-d(t)/2<r(t+1)≤r(t),α2<α1 |
s8 | r(t+1)≤r(t)-d(t)/2,α2<α1 |
动作集合A里面的a1,a2...a8分别为顺时针旋转15°,30°...120°。a9,a10...a16分别为逆时针旋转15°,30°,...120°。
艏向角调整在t+1时刻调整,调整角度为a1,a2....a16中的一个。训练过程中动作选取策略为ε-greedy策略。ε是一个0到1之间的数,越大说明采取随机动作的概率越大,这里设置为0.9,使得在训练过程中可以充分探索各个动作。调整前,先生成一个0到1之间的随机数,若这个随机数大于ε,则随机选取A集合中的一个动作。小于ε则根据表格,通过情况的分类找到在表格对应的行名,查找到这一行中值最大的那一格,选取那一格列名所对应的动作。
选取动作后,执行该调整角度的动作,并直行。测得直行后,水下机器人到信标的距离r(t+2),计算奖励值,用来作为一个判断该情况下采取该动作好坏的一个信号。目的是使得水下机器人做接近信标的运动。因此当距离变小时,奖励值为+1,距离变大时,奖励值为-1。即当r(t+2)≤r(t+1)时,奖励值Re=+1;当r(t+2)≥r(t+1)时,奖励值Re=-1。后时刻比前时刻距离小的话奖励值为正。目的是为了接近信标,距离减小是好的,给正的奖励值。奖励值用于更新R表格,表格的更新使艏向角调整策略也会改变。更新的公式为R(s,a)=(R(s,a)×(N(s,a)-1)+Re)/N(s,a).s为情况s1~s8中的一个,a为动作a1~a16,具体为得到奖励值Re前,所归类的情况s,和所执行的动作a。
N(s,a)为N(S,A)表格中的一项,表示训练过程中情况s下执行了动作a的次数。N(S,A)是一个8x16的表格。更新规则为N(s,a)=N(s,a)+1。
R(s,a)为R(S,A)表格中的一项,表示情况s下执行动作a获得奖励值的平均值。R(S,A)是一个8x16的表格。
本发明在动作选取策略上采取了ε-greedy随机策略,在训练过程中,每个动作都会被充分选择,这样R(s,a)的值会趋于该值的数学期望值。
训练过程如图3流程图所示。
步骤1:初始化R表格为8x16的一张空表格。设置训练次数
步骤2:设置水下机器人的起始位置。
在每个时刻t测得观测距离r(t),结合上一个时刻的位移d(t-1),和上一个时刻的观测距离r(t-1),如图1所示,进行判断是否要调整角度。若r(t)<r(t-1),且夹角α小于45°,保持直行。
步骤3:不然先逆时针旋转5度,直行,获得下一个时刻的观测距离和位移r(t+1)。
步骤4:如图2所示。由三个时刻的观测距离和两个位移信息,计算特征1和2,由特征将所有情况归为八类。对每一个情况按ε-greedy策略选取动作。选取动作后,执行该动作,直行,并获得下一个时刻的距离r(t+2)。
步骤5:距离变小,奖励值Re为+1,距离变大,奖励值Re为-1。获得奖励值Re,并更新表格。
步骤6:当观测距离小于D。本轮训练结束,返回步骤2。训练达到设定次数,训练结束。保留训练结果R表格。
扩展卡尔曼滤波位置估计方程:
P(k+1|k)=P(k|k)+Q
K(k+1)=P(k+1|k)·HT(k+1)·[H(k+1)·P(k+1|k)·HT(k+1)+R]-1
P(k+1|k+1)=[I-K(k+1)·H(k+1)]·P(k+1|k)
上式中,k表示k时刻,X(k)=[xk yk]T上标“T”表示向量或者矩阵的转置xk,yk分别是水下机器人在k时刻以信标为原点的东向和北向坐标;宁母上带“^”符号的表示该状态量的预测值或估计值。方程中是k时刻的状态估计值,是k+1时刻的状态预测值,是根据观测值校正后的k+1时刻的状态估计值,
U(k)=[v·sinθkv·cosθk]T θk为k时刻的水下机器人的艏向角,v为水下机器人的航速;
I为单位矩阵P(k|k)是k时刻的误差的协方差矩阵;P(k+1|k)是k+1时刻预测的误差协方差矩阵;P(k+1|k+1)是量测更新后的协方差矩阵。
V(k+1)为k+1时刻过程噪声,W(k+1)为k+1时刻观测噪声,满足零均值高斯分布,Q是过程噪声的方差,R是观测噪声的方差;
Z(k+1)=h(X(k+1))+W(k+1),是即k+1时刻观测到的带有噪声的水下机器与信标之间的距离值;
根据R表格,对于情况si(i=1,2,3...,8),比较R(si,a1),R(si,a2),R(si,a3)...R(si,a16)的值,选则其中最大R值所对应的动作,作为需要调整的角度,使水下机器人向信标运动。在水下机器人向信标运动过程中,用扩展卡尔曼滤波进行位置估计。
水下机器人匀速率运动,在每个时刻测得到单信标的观测距离。根据速度信息,可获得两个时刻之间的位移信息。在开始时刻,向任意方向运动。由两个时刻测得的到信标的观测距离r(t-1)、r(t)和这两个时刻之间的位移距离d(t-1),构成一个三角形。进行一次判定,若r(t)<r(t-1),且三角形中边d(t-1)与边r(t-1)之间夹角小于45度,则保持直行。否则,在当前时刻逆时针旋转5°并直行;进入下一个时刻,测得下一个时刻到单信标的观测距离r(t+1)和下一个时刻与当前时刻之间的位移距离d(t)。
由三个时刻的到信标的距离和三个时刻间的两个位移距离信息可以组成两个三角形。对这两个三角形根据距离变化和夹角变化分为八种情况。对每一种情况,根据R表格,选择一个调整的角度。R表格为强化学习训练的结果。调整角度后,直行,获得下一个时刻到信标的距离r(t+2)。
重复以上步骤,直到水下机器人到信标的距离小于D。之后在信标附近运动一段时间T,定位结束。
实施例:
对本发明提供的一种接近信标轨迹下的水下机器人定位方法进行仿真分析。
水下机器人实际起始位置为(400m,-300m),扩展卡尔曼滤波估计的起始位置为(200m,-100m)。水下机器人做2m/s的匀速率运动。测距周期为1s。
量测噪声方差V=[1]m
每个时刻测得观测距离,结合位移信息,判断是否进行角度调整。需要调整角度时,对情况进行分类,并根据R表格做角度调整。直到观测距离小于2m时,再在信标附近做一段时间运动。
Claims (3)
1.一种接近信标轨迹下的水下机器人定位方法,其特征在于,包括以下步骤:
(1)在保持水下机器人匀速率运动且速率大小已知的情况下,按照设定测量周期,以单个声信标测量水下机器人与信标之间的距离值,以水下机器人自带的罗盘测量艏向角,由前一个时刻的位移距离和前一时刻的观测距离,和当前时刻的观测距离,通过已知三边长度的余弦定理计算得到当前位移方向与当前观测距离方向之间的夹角;然后通过强化学习和训练,获得艏向角的调整策略;
(2)利用强化学习训练得到的策略来调整艏向角,使水下机器人朝向接近信标的方向运动;具体包括:
(2.1)按设定时间间隔,测量当前时刻的观测距离r(t+1)、上一时刻的观测距离r(t)、上一时刻的位移距离d(t)和用已知三边长度下的余弦定理,计算位移方向与当前观测距离方向之间的夹角α2;如果r(t+1)<r(t)且夹角α2小于45°,保持直行;否则将水下机器人的行进方向逆时针旋转5°后直行,并测量获得下一个时刻的观测距离、位移和夹角数据;
(2.2)调整艏向角后,继续直行,并重新进行测量和计算;重复以上过程,直到水下机器人到信标的距离小于设定的距离值D;然后围绕信标方向随机地运动一段设定时间T,定位结束;
(3)在接近信标的过程中,利用扩展卡尔曼滤波位置估计方程计算水下机器人的位置以实现定位。
2.根据权利要求1所述的方法,其特征在于,所述步骤(1)具体包括:
首先,建立8行×16列的R表格,表格内容是对应于各情况下的平均奖励值,表示选择相应动作时对接近信标的有利程度;其中,
行名包括下述八种情况:情况1:观测距离变大,变化量大于一半的位移距离且位移方向与观测距离方向的夹角变大;情况2:观测距离变大,变化量小于等于一半的位移距离,且位移方向与观测距离方向的夹角变大;情况3:观测距离变小,变化量小于一半的位移距离,且位移距离与观测距离方向的夹角变大;情况4:观测距离变小,变化量大于等于一半的位移距离,且位移距离与观测距离方向的夹角变大;情况5:观测距离变大,变化量大于一半的位移距离,且位移方向与观测距离方向的夹角变小;情况6:观测距离变大,变化量小于等于一半的位移距离,且位移方向与观测距离方向的夹角变小;情况7:观测距离变小,变化量小于一半的位移距离,且位移距离与观测距离方向的夹角变小;情况8:观测距离变小,变化量大于等于一半的位移距离,且位移距离与观测距离方向的夹角变小;
列名包括下述十六个调整动作:顺时针旋转15°、30°、45°、60°、75°、90°、105°、120°和逆时针旋转15°、30°、45°、60°、75°、90°、105°、120°;
其次,以经过强化学习训练的R表格作为艏向角的调整策略,具体训练的过程包括:
(1.1)初始化R表格,确定训练次数;
(1.2)设置水下机器人的起始位置后,使其保持匀速率运动;在第一个时刻里,水下机器人朝任意方向移动;
(1.3)按设定时间间隔,测量当前时刻的观测距离r(t)、上一时刻的观测距离r(t-1)、上一时刻的位移距离d(t-1)和用已知三边长度下的余弦定理,计算位移方向与当前观测距离方向之间的夹角α;如果r(t)<r(t-1)且夹角α小于45°,保持直行;否则将水下机器人的行进方向逆时针旋转5°后直行,并测量获得下一个时刻的观测距离、位移和夹角数据;
(1.4)由连续三个时刻的观测距离和两次位移距离的信息,计算出后两次观测距离的差值r(t+1)-r(t)和两次夹角的差值α2-α1,并在细分后进行归类;
(1.5)以ε-greedy的探索策略选取设计好的十六个动作;执行动作后,根据观测距离变化情况获得对应奖励值Re;当观测距离减小时Re=+1,当观测距离不变或变大时Re=-1;并更新R表格;更新公式为:
R(s,a)=(R(s,a)×(N(s,a)-1)+Re)/N(s,a)
s为情况s1~s8中的一个,a为动作a1~a16,具体为得到奖励值Re前,所归类的情况s,和所执行的动作a;
N(s,a)为训练过程中N(S,A)表格中的一项,表示在情况s下执行了动作a的次数;N(S,A)是一个8×16的表格;更新规则为N(s,a)=N(s,a)+1;
R(s,a)为训练过程中R(S,A)表格中的一项,表示在情况s状态下执行动作a获得奖励值的平均值,R(S,A)是一个8×16的表格;
(1.6)当观测距离小于设定的距离值D;本轮训练结束,返回步骤(1.3);当训练达到设定次数后结束训练,保留训练结果R表格。
3.根据权利要求1所述位方法,其特征在于,所述步骤(3)包括:
(1)设定扩展卡尔曼滤波位置估计方程为:
K(k+1)=P(k+1|k)·HT(k+1)·[H(k+1)·P(k+1|k)·HT(k+1)+R]-1
P(k+1|k+1)=[I-K(k+1)·H(k+1)]·P(k+1|k)
上式中,k表示k时刻,X(k)=[xk yk]T上标“T”表示向量或者矩阵的转置;xk,yk分别是水下机器人在k时刻以信标为原点的东向和北向坐标;字母上带“^”符号的表示状态量的预测值或估计值;
U(k)=[v·sinθk v·cosθk]T,θk为k时刻的水下机器人的艏向角,v为水下机器人的航速;
I为单位矩阵;P(k|k)是k时刻的误差的协方差矩阵;P(k+1|k)是k+1时刻预测的误差协方差矩阵;P(k+1|k+1)是量测更新后的协方差矩阵;
V(k+1)为k+1时刻过程噪声,W(k+1)为k+1时刻观测噪声,满足零均值高斯分布,Q是过程噪声的方差,R是观测噪声的方差;
Z(k+1)=h(X(k+1))+W(k+1),是即k+1时刻观测到的带有噪声的水下机器与信标之间的距离值;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911064151.2A CN110888104B (zh) | 2019-11-04 | 2019-11-04 | 一种接近信标轨迹下的水下机器人定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911064151.2A CN110888104B (zh) | 2019-11-04 | 2019-11-04 | 一种接近信标轨迹下的水下机器人定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110888104A CN110888104A (zh) | 2020-03-17 |
CN110888104B true CN110888104B (zh) | 2022-03-22 |
Family
ID=69746784
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911064151.2A Active CN110888104B (zh) | 2019-11-04 | 2019-11-04 | 一种接近信标轨迹下的水下机器人定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110888104B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112197765B (zh) * | 2020-09-14 | 2021-12-10 | 浙江大学 | 一种实现水下机器人精细导航的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103376452A (zh) * | 2012-04-18 | 2013-10-30 | 中国科学院沈阳自动化研究所 | 一种用单台声信标修正水下机器人位置误差的方法 |
CN104280024A (zh) * | 2013-07-05 | 2015-01-14 | 中国科学院沈阳自动化研究所 | 一种深水机器人组合导航装置和方法 |
CN105823480A (zh) * | 2016-03-18 | 2016-08-03 | 中国海洋大学 | 基于单信标的水下移动目标定位算法 |
CN107102644A (zh) * | 2017-06-22 | 2017-08-29 | 华南师范大学 | 基于深度强化学习的水下机器人轨迹控制方法及控制*** |
CN107741745A (zh) * | 2017-09-19 | 2018-02-27 | 浙江大学 | 一种实现移动机器人自主定位与地图构建的方法 |
CN109782289A (zh) * | 2018-12-26 | 2019-05-21 | 中国电子科技集团公司第二十研究所 | 一种基于基线几何结构约束的水下航行器定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9213077B2 (en) * | 2012-07-31 | 2015-12-15 | Adaptive Methods, Inc. | GPS assisted torpedo recovery system |
-
2019
- 2019-11-04 CN CN201911064151.2A patent/CN110888104B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103376452A (zh) * | 2012-04-18 | 2013-10-30 | 中国科学院沈阳自动化研究所 | 一种用单台声信标修正水下机器人位置误差的方法 |
CN104280024A (zh) * | 2013-07-05 | 2015-01-14 | 中国科学院沈阳自动化研究所 | 一种深水机器人组合导航装置和方法 |
CN105823480A (zh) * | 2016-03-18 | 2016-08-03 | 中国海洋大学 | 基于单信标的水下移动目标定位算法 |
CN107102644A (zh) * | 2017-06-22 | 2017-08-29 | 华南师范大学 | 基于深度强化学习的水下机器人轨迹控制方法及控制*** |
CN107741745A (zh) * | 2017-09-19 | 2018-02-27 | 浙江大学 | 一种实现移动机器人自主定位与地图构建的方法 |
CN109782289A (zh) * | 2018-12-26 | 2019-05-21 | 中国电子科技集团公司第二十研究所 | 一种基于基线几何结构约束的水下航行器定位方法 |
Non-Patent Citations (1)
Title |
---|
用固定单信标修正水下机器人导航误差;刘健等;《控制与决策》;20100930;第25卷(第9期);第1354-1357页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110888104A (zh) | 2020-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109211276B (zh) | 基于gpr与改进的srckf的sins初始对准方法 | |
CN101464152B (zh) | 一种sins/gps组合导航***自适应滤波方法 | |
CN109724599B (zh) | 一种抗野值的鲁棒卡尔曼滤波sins/dvl组合导航方法 | |
CN109633590B (zh) | 基于gp-vsmm-jpda的扩展目标跟踪方法 | |
CN108120438B (zh) | 一种基于imu和rfid信息融合的室内目标快速跟踪方法 | |
CN103644903B (zh) | 基于分布式边缘无味粒子滤波的同步定位与地图构建方法 | |
CN109521454A (zh) | 一种基于自学习容积卡尔曼滤波的gps/ins组合导航方法 | |
CN103439731A (zh) | 基于无迹卡尔曼滤波的gps/ins组合导航方法 | |
CN110726415B (zh) | 一种自适应的水下多波束同步定位与建图方法 | |
CN101819041A (zh) | 自进化anfis与ukf结合的gps/mems-ins组合定位误差动态预测方法 | |
CN110673148A (zh) | 一种主动声纳目标实时航迹解算方法 | |
CN111508282B (zh) | 低空无人机农田作业飞行障碍物冲突检测方法 | |
CN110954132A (zh) | Grnn辅助自适应卡尔曼滤波进行导航故障识别的方法 | |
CN107179693A (zh) | 基于Huber估计的鲁棒自适应滤波和状态估计方法 | |
CN112083457B (zh) | 一种神经网络优化的imm卫星定位导航方法 | |
CN110763872A (zh) | 一种多普勒测速仪多参数在线标定方法 | |
CN110888104B (zh) | 一种接近信标轨迹下的水下机器人定位方法 | |
CN114777812B (zh) | 一种水下组合导航***行进间对准与姿态估计方法 | |
CN108226887A (zh) | 一种观测量短暂丢失情况下的水面目标救援状态估计方法 | |
CN102508217B (zh) | 建立雷达测量误差标定模型的方法 | |
CN108507587B (zh) | 一种基于坐标变换的三维车载定位导航方法 | |
CN116859356B (zh) | 基于旋转矩阵优化求解的车载4d毫米波雷达自标定方法 | |
CN113503891B (zh) | 一种sinsdvl对准校正方法、***、介质及设备 | |
CN114061592B (zh) | 基于多模型的自适应鲁棒auv导航方法 | |
CN105549003A (zh) | 一种汽车雷达目标跟踪方法 |
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 |