[发明专利]一种基于自适应EKF的无人侦察机同步定位与构图方法

权利要求书

1.一种基于自适应EKF的无人侦察机同步定位与构图方法，其特征在于，该方法包括如下步骤：

步骤S1：对无人侦察机SLAM系统建模；

1.1地图模型

采用特征地图来构建环境地图模型，将环境中的目标用几何原型来描述；

1.2坐标系统

由于无人侦察机所携带的传感器与环境特征之间的感知行为发生在传感器坐标系中，速度的测量发生在无人侦察机运动坐标系中，而最终环境地图的构建需要表示在全局地图中，所以建立全局坐标系，无人侦察机运动坐标系及传感器坐标系，并表明他们之间的关系，定义GX_GY_GZ_G为全局坐标系，G为地心，X_G指向正北，Y_G指向正东，Z_G垂直于X_GY_G平面并指向G；定义AX_AY_AZ_A为无人侦察机运动坐标系，A为无人侦察机质心，X_A指向机头，Y_A垂直于机头，Z_A垂直于X_AY_A平面并指向A；定义SX_SY_SZ_S为环境感知传感器坐标系，S为传感器质心，X_S指向传感器电磁波发射方向，Y_S垂直于电磁波发射方向，Z_S垂直于X_SY_S平面并指向S；三个坐标系均符合右手定则；

1.3特征模型

根据公式(1)，定义特征的运动学模型为：

Xfi(k)=Xfi(k-1),i=1...n---(1)]]>

其中，k为离散时间变量；

特征在传感器坐标系中一般以极坐标的形式表示：

Xfi(k)=ρiθi(k)]]>

其中，表示第i个特征的状态，即位置；n为特征的个数；

ρ_i为特征i在极坐标系下的极径；θ_i为特征i在极坐标系下的极角；

1.4无人侦察机运动模型

通过建立一个四自由度的常速运动学模型来描述无人侦察机的运动变换趋势，如公式(2)所示：

X_A(k)＝f(X_A(k-1),n(k-1))(2)

无人侦察机X_A＝[xyzψuvwr]^T，[x,y,z,ψ]为全局坐标系中无人侦察机的位置和艏向，[u,v,w,r]表示运动坐标系A中相应的X_A、Y_A、Z_A方向的线速度和转艏角速度，n＝[n_u,n_v,n_w,n_r]表示以加速度的形式作用在速度上的高斯噪声；

根据特征状态和无人侦察机状态X_A可以获得SLAM系统的状态X，包括无人侦察机的状态X_A和特征的状态X_f：

X＝[X_AX_f]^T(3)

xyzψuvwr(k)=x+(uT+nuT22)cosψ-(vT+nvT22)sinψy+(uT+nuT22)sinψ+(vT+nvT22)cosψz+wT+nwT22ψ+rT+nrT22u+nuTv+nvTw+nwTr+nrT(k-1)]]>

其中，T为采样时间；

1.5传感器测量模型

(1)速度测量模型

速度测量传感器可提供X、Y、Z三个方向的速度，测量模型为：

Z_V＝H_VX+s_V(4)

其中，Z_V为速度的测量值；

测量矩阵H_V为：

H_V＝[0_3×4I_3×30_3×10_3×2n]

其中，O为全零矩阵，I为单位矩阵；

速度测量噪声s_V的协方差矩阵为：

RV=σVu2000σVv2000σVw2]]>

其中，为速度u的方差，为速度v的方差，为速度w的方差；

(2)高度测量模型

通过高度测量传感器可以获得无人侦察机的高度，测量模型为：

Z_H＝H_HX+s_H(5)

其中，Z_H为高度的测量值；

测量矩阵H_H为：

H_H＝[001000000_1×2n]

高度测量噪声的s_H协方差矩阵表示为：

RH=σH2]]>

(3)艏向测量模型

通过艏向测量传感器可以获得无人侦察机的艏向，测量模型为：

Z_C＝H_CX+s_C(6)

其中，Z_C为艏向的测量值；

测量矩阵H_C为：

H_C＝[000100000_1×2n]

艏向测量噪声s_C的协方差为：

RC=σC2]]>

(4)环境测量模型

通过环境感知传感器可以获得环境中特征相对于无人侦查机的距离和方位，测量模型为：

ZiA=ρiAθiAT=h(X(k),si)---(7)]]>

其中，s_i为特征测量噪声；

h(X(k),si)=(ρicosθi-x)2+(ρisinθi-y)2arctan(ρisinθi-yρicosθi-x)-ψ+si]]>

测量矩阵为：

其中，和是点特征在运动坐标系下的参数表示，ρ_i和θ_i是点特征在全局坐标系下的参数表示；

步骤S2：基于步骤S1所建立的无人侦察机SLAM系统模型，执行自适应EKF的无人侦察机同步定位与构图操作；

其中，建立无人侦察机SLAM系统的状态方程和观测方程如下所示，被估计状态以及观测变量与状态向量的关系呈现非线性：

X(k)＝f(X(k-1),m(k))

(8)

Z(k)＝h(X(k),n(k))

式中，X(k)是无人侦察机SLAM系统的状态向量；Z(k)是系统的观测向量；m(k)是系统的过程噪声序列；n(k)是系统的观测噪声序列；f()和h()表示系统的非线性状态方程和非线性观测方程；

针对上式所示的非线性SLAM系统，将系统过程噪声以及观测噪声建模为带有时变均值及协方差矩阵的相互独立的正态白噪声序列：

E[m(k)]=q(k),E[(m(k)-q(k))(m(j)-q(j))T]=Q(k)δkjE[n(k)]=r(k),E[(n(k)-r(k))(n(j)-r(j))T]=R(k)δkjE[(m(k)-q(k))(n(j)-r(j))T]=0]]>

EKF利用泰勒级数展开的思想通过求解非线性状态方程和观测方程的偏导数来将SLAM系统线性化；若k-1时刻X(k-1)的最优状态估计则X(k)的估计可按Kalman滤波基本方程求解；进而按照下述带时变噪声估值器的自适应EKF-SLAM基本步骤流程执行：

步骤A：执行航位推算的无人侦察机随着不断航行，利用所携带的环境感知传感器不断感知到环境中的特征，并提取所获取图像中的特征；

步骤B：每次提取到特征需要判断当前特征是新特征或是无人侦察机已经观测过并融入地图中的旧特征，即完成数据关联过程；

步骤C：如果数据关联结果表明所观测特征为地图中已存在特征，则根据对该特征的实际测量值与预测测量值之间的差异来实现无人侦察机状态及整个地图中所有特征位置的校正；具体如下：

SLAM过程中包含了系统状态向量的前二阶矩：

E[X(k)]=X^(k)]]>

E[(X(k)-X^(k))(X(k)-X^(k))T]=P(k)]]>

其中，表示无人侦察机及地图中特征状态的估计值，P(k)表示无人侦察机及特征状态估计误差的协方差矩阵；

步骤C-1：状态预测过程，利用状态方程对系统的状态估计及协方差矩阵进行预测：

X^(k|k-1)=F(k)X(k-1)+q^(k-1)---(9)]]>

P(k|k-1)=F(k)P(k-1)F(k)T+W(k)Q^(k-1)W(k)T---(10)]]>

其中，F(k)=FV(k)08×2n02n×8I2n×2n,W(k)=WN(k)02n×4,]]>F_V(k)是非线性运动方程f对状态X_A在当前最优估计处的雅克比矩阵；W(k)是非线性状态方程f对噪声m(k)的雅克比矩阵；由时变噪声统计特性估值器递推获得：

q^(k)=1kΣi=0k-1[X^(i+1|i+1)-F(k)X^(i|i)]---(11)]]>

Q^(k)=1kΣi=0k-1[K(i+1)υ(i+1)υ(i+1)TK(i+1)T+P(i+1|i+1)-F(i)P(i|i)F(i)T]---(12)]]>

步骤C-2：SLAM的状态校正通过速度、高度、艏向及环境感知传感器的测量，来实现无人侦查机状态及地图中特征状态的更新，测量更新方程为：

X^(k)=X^(k|k-1)+Kυ(k)---(13)]]>

P(k)＝[I-K(k)H(k)]P(k|k-1)(14)

K(k)＝P(k|k-1)H(k)^T(S(k))^-1S(k)(15)

υ(k)=Z(k)-H(k)X^(k|k-1)-r^(k)---(16)]]>

S(k)=H(k)P(k|k-1)H(k)T+R^(k-1)---(17)]]>

其中，H是非线性传感器观测方程h对X的偏导数的雅克比矩阵，H根据测量传感器的不同而具有相应的不同形式；和是由时变噪声统计估值器递推获得：

r^(k)=1kΣi=0k-1[Z(i+1)-H(i+1)X^(i+1|i)]---(18)]]>

R^(k)=1kΣi=0k-1[υ(i+1)υ(i+1)T-H(i+1)P(i+1|i)H(i+1)T]---(19)]]>

将式(9)～(19)迭代计算即可实现噪声统计特性的递推估计，进而实现带时变噪声统计特性估值器的自适应EKF-SLAM过程。