基于声学测量网络的水下多AUV协同定位编队拓扑结构优化方法与流程

本发明属于声学测量领域，具体涉及一种基于声学测量网络的水下多auv协同定位编队拓扑结构优化方法。

背景技术：

协同定位是目前中间层区域多自主式水下航行器最有效的导航方法之一。航行器编队拓扑结构会对整个系统的可观测性产生直接影响，从而影响多自主水下航行器协同定位精度。因此，考虑水声测距误差与距离的相关性，对协同定位系统的编队拓扑结构进行优化，使系统的可观测性达到最优，是协同定位方向需要深入研究的课题。目前针对多领航、多随从auv系统，多采用fisher信息矩阵推导评价函数，采取基于偏导数的梯度算法进行位置的迭代。然而，该方法的求解过程依赖于初始位置的选取，若初始位置选取不当，往往存在陷入局部最优解的危险，而无法得到最优拓扑结构。

本发明考虑和距离相关的水声噪声，采用引入退火思想的步进递推策略对多auv协同定位系统的最优编队拓扑构型的优化方法进行了研究。现在常用的令对随从auv取fisher信息矩阵的行列式作为评价函数的方法，在auv编队规模日益扩大的情况很难通过计算机解算得到最优编队拓扑结构；基于偏导数的迭代策略虽然可以在一定程度上避免多解或无解情况的产生，但在初始位置选取不当时往往受限于局部最优解而无法得到真正最优的拓扑结构。

针对目前在设计多auv协同定位系统编队构型的求解困难的问题，本发明在基于fisher信息矩阵理论的基础上，综合考虑多随从auv、随从auv位置不确定以及与距离相关的水声噪声来推导评价函数，提出一种基于退火思想的步进递推策略，最终实现快速递推得到最优编队拓扑结构。本发明可用于多auv协同定位系统的编队构型优化领域。

技术实现要素：

本发明的目的在于提供一种能够在随从auv位置不确定情形下，领航auv从任意初始位置出发，完成多auv协同定位系统编队构型快速优化设计的方法。

本发明的目的是这样实现的：

基于声学测量网络的水下多auv协同定位编队拓扑结构优化方法，该方法包括以下步骤：

步骤1：建立多auv协同定位模型；

步骤2：推算各随从auv的fisher信息矩阵；

步骤3：根据各随从auv的fisher信息矩阵，构建系统评价函数；

步骤4：根据系统评价函数，计算最优编队拓扑结构约束条件；

步骤5：引入退火思想完成步进递推过程；

步骤5.1：初始化；

步骤5.2：在邻域内产生新解；

步骤5.3：通过得到的新的领航auv的位置参数xinew和yinew计算最优评价函数fnew；

步骤6：根据约束条件验证最终编队拓扑结构，找到最优位置布局。

本发明还包括：

步骤1所述的多auv协同定位模型为：

(i＝1,2,3...n，j＝1,2,3,...m)

式中，rij为第i个领航auv与第j个随从auv之间的真实距离；lij为第i个领航auv与第j个随从auv之间的量测距离；pi＝[xi,yi,zi]^t为第i个领航auv在导航坐标系中的位置，xi为第i个领航auv在x轴方向的坐标，yi为第i个领航auv在y轴方向的坐标，zi为第i个领航auv在z轴方向的坐标；qj＝[xj,yj,zj]^t为第j个随从auv在导航坐标系中的位置，xj为第j个随从auv在x轴方向的坐标，yj为第j个随从auv在y轴方向的坐标，zj为第j个随从auv在z轴方向的坐标；n为领航auv的数目，m为ωij为随从auv的数目；噪声信号，且ωij＝σ(1+ηrij^γ)，σ为量测噪声常数；η为量测噪声相关的常数，γ为距离相关的常数。

步骤2所述的推算各随从auv的fisher信息矩阵为：

式中：向量

步骤3所述的系统评价函数表达式为：

式中：f(qj)为概率密度函数，其具体形式会根据随从auv的工作任务决定；dj为随从auv可能分布的范围，根据概率密度函数决定，用蒙特卡洛方法化简为如下形式：

式中，第一参数aijv、第二参数bijv、第三参数mijv表示为：

jυ为一定区域内，第υ个可能的随从auv位置，在该域内根据构建的可观测性评价函数，解得从艇j的最优位置。

步骤4所述的计算最优编队拓扑结构约束条件为：

用蒙特卡洛方法化简成如下形式：

步骤5.1所述的初始化包括以下步骤：

(5.1.1)初始化温度t；

(5.1.2)初始化解状态；

(5.1.3)初始化每个t值时的迭代次数k；

对每一个随从auv来说，通过当前各个领航auv的横坐标xi和纵坐标yi计算其fisher信息矩阵的行列式|fim(qj)|，之后将fisher信息矩阵的行列式取对数并求和，即得到在当前步骤的评价函数，记为f[t]，作为初始化的解状态；其中温度下降率趋于1。

步骤5.3所述计算最优评价函数fnew，之后会产生两种情况：

情况1：fnew＞fold，则说明新的领航auv位置能够提升系统的可观测性，将新的领航auv编队拓扑结构作为基础，回到步骤2进行下一步计算；

情况2：fnew＜fold，则说明新的领航auv位置不一定能够提升系统的可观测性，此时不能排除局部最优解的情况，因此以metropolis准则接受新解作为当前解，同样回到步骤2进行下一步计算；

计算得到新的最优评价函数fnew导致温度变化，使得领航auv的位置参数发生变化如下式所示：

xi_new＝xi_old+μ·t

yi_new＝yi_old+μ·t

其中，μ可选择为(-1,1)上任意一个常数，此时产生新解的邻域的大小与温度高低直接相关，随着温度的降低，产生新解的邻域范围逐渐缩小。

本发明的有益效果在于：

1.本发明在多随从auv情况下，考虑了声学测距误差与距离的相关性，具有更高的实用价值；

2.本发明针对位置信息的不确定性，根据相应的概率密度采用蒙特卡洛方法对可能分布区域内的编队构型进行优化设计，针对实际中可能出现的问题有较好的应对方案；

3.本发明采用基于退火思想的步进递推的策略，递推过程的进行与初始值的选取无关，也就是说，不管主艇初始位置在何处，均可经过迭代步骤找到其最优位置布局；

4.本发明引入metropolis准则作为判断是否接受新解作为当前解的准则之一，可有效改善当局部最优解出现时迭代不再继续进行的情况，是一种以概率l收敛于全局最优解的全局优化算法；

5.本发明的迭代过程产生新解邻域的大小与温度高低直接相关，随着温度的降低，产生新解的邻域逐渐缩小至一个极小的范围，这样可以在递推的过程中逐步减小领航auv位置递推移动的距离，增加最终结果的精确性。

附图说明

图1为本发明的流程图。

图2为随从auv位置不确定时交叉计算示意图。

图3为角度投影示意图。

具体实施方式

下面结合附图对本发明做进一步描述：

步骤一：建立多auv协同定位模型；

建立多auv协同定位模型时，采用的模型为：

式中：

ωij＝σ(1+ηrij^γ),

(i＝1,2,3...n，j＝1,2,3,...m)

lij为第i个领航auv与第j个随从auv之间的量测距离；pi＝[xi,yi,zi]^t为第i个领航auv在导航坐标系中的位置；qj＝[xj,yj,zj]^t为第j个随从auv在导航坐标系中的位置；ωij为噪声信号；rij为第i个领航auv与第j个随从auv之间的真实距离；σ为量测噪声常数；η和γ均为量测噪声和距离相关的常数。

步骤二：推算各随从auv的fisher信息矩阵；

推算各随从auv的fisher信息矩阵时，采用的算法模型为：

式中：

步骤三：推导系统评价函数表达式；

考虑随从auv位置的不确定性，推导系统评价函数表达式时采用的算法模型为：

式中：f(qj)为概率密度函数，其具体形式会根据随从auv的工作任务决定；dj为随从auv可能分布的范围，根据概率密度函数决定。用蒙特卡洛方法可以化简为如下形式：

式中：

jυ为一定区域内，第υ个可能的随从auv位置，如令υ＝10。在该域内根据构建的可观测性评价函数，解得从艇j的最优位置。交叉计算示意图如图2所示。

aij＝cos(αij)cos(βij)

bij＝sin(αij)cos(βij)

αij为第j个随从auv到第i个领航auv的向量在水平面的投影与x轴在逆时针方向上的夹角，其范围为0°到360°，βij为第j个随从auv到第i个领航auv的向量在垂直于水平面的二维平面上投影与水平面的夹角，其范围为0°到90°。角度投影示意图如图3所示。

步骤四：计算最优编队拓扑结构约束条件；

计算最优编队拓扑结构约束条件时，采用的算法模型为：

用蒙特卡洛方法可以化简成如下形式：

步骤五：引入退火思想完成步进递推过程；

具体实施步骤如下：

(1)初始化：

a.初始化温度t；

b.初始化解状态；

对每一个随从auv来说，通过当前各个领航auv的横坐标xi和纵坐标yi计算其fisher信息矩阵的行列式|fim(qj)|，之后将fisher信息矩阵的行列式取对数并求和，即可以得到在当前步骤的评价函数，记为f[t]。

c.初始化每个t值时的迭代次数k。

其中，初始温度应足够大，温度下降率趋于1，保证温度下降足够慢，同时每个t值时的迭代次数k足够大，从而保证解状态跑遍整个区域。

(2)在邻域内产生新解；

(3)通过得到的新的领航auv的位置参数xinew和yinew计算最优评价函数fnew。此时可能会产生两种情况：

·如果fnew＞fold，则说明新的领航auv位置能够提升系统的可观测性。那么将新的领航auv编队拓扑结构作为基础，回到步骤(2)进行下一步计算。

·如果fnew＜fold，则说明新的领航auv位置不一定能够提升系统的可观测性。由于此时不能排除局部最优解的情况，因此以metropolis准则接受新解作为当前解，同样回到步骤(2)进行下一步计算。

温度变化使得领航auv的位置参数发生变化如下式所示：

xi_new＝xi_old+μ·t

yi_new＝yi_old+μ·t

其中，μ可选择为(-1,1)上任意一个常数，此时产生新解的邻域的大小与温度高低直接相关。随着温度的降低，产生新解的邻域范围逐渐缩小。

步骤六：根据约束条件验证最终编队拓扑结构。