一种基于视觉雷达的小天体探测器多模型导航方法

1.本发明属于深空探测器自主导航领域，特别适用于深空探测器跟踪附着小天体的过程，具体涉及一种基于视觉雷达的小天体探测器柔性多模型导航方法。


背景技术：

2.近年来，深空探测不断发展，小天体探测也已经成为一项艰巨且具有科学意义的深空探测任务。然而，与行星探测任务相比，小天体探测任务具有探测目标小、所处环境多变、引力小、着陆点不规则等特点。当前小天体探测器在进行小天体附着时，由于单个探测器、单一传感器的限制，以及对目标运动模型建模的不确定，导致探测器无法准确跟踪小天体；由于探测器与小天体接触时发生刚性碰撞，导致小天体探测器附着时发生反弹而失败，探测器无法准确附着小天体，不能完成小天体样品采集的任务。
3.柔性材料发展迅速，柔性材料以其刚柔耦合特性，成为一种极具应用前景的航天探测器材料。但是，柔性材料有复杂的非线性特性，这使得柔性材料结构复杂且极易受到扰动。将雷达与视觉导航方法结合，提高导航精度，增强准确性；与此同时，运用多模型交互导航的方式，引入柔性探测器的不同动力学模型，将柔性材料与多模型交互导航算法结合，增强组合导航的可靠性的同时，大大加强了探测器附着小天体成功的概率。交互式多模型算法作为一种常见的数据融合算法，精度高，可靠性强，用于深空小天体探测任务，与组合导航算法和新型的柔性材料结合，完成深空探测任务，国内对该方面的研究仍在不断进行中。


技术实现要素：

4.本发明的目的是解决小天体跟踪精度低、小天体与探测器刚性碰撞导致附着失败等问题，提出了一种基于视觉雷达小天体多模型柔性探测器导航方法。将视觉/雷达探测器与柔性材料相结合，建立小天体柔性探测器物理模型；同时根据柔性材料的刚柔耦合特性，建立不同位置状态下的多种简化动力学模型；研究视觉/雷达数据融合目标导航算法，设计了将视觉和雷达信息作为观测量的卡尔曼滤波器，基于该滤波器和不同状态下的动力学模型，运用了一种视觉/雷达的多模型导航算方法，提高导航精度，进而解决柔性探测器在小天体表面附着等方面问题。
5.本发明是通过下述技术方案实现的。
6.本发明公开的柔性探测器多模型导航方法，根据刚柔耦合特性对柔性探测器上的多个探测器点建立多种动力学观测模型，将视觉雷达融合信息作为输入量，引入多种动力学模型，实现整体柔性探测器对小天体表面特征点的导航：将视觉/雷达探测器与柔性材料相结合，建立柔性探测器物理模型；基于柔性体的刚柔耦合特性建立不同位置状态下的多种动力学模型；运用视觉/雷达数据融合目标跟踪算法，提出了将视觉和雷达信息作为观测量的卡尔曼滤波器，基于该滤波器和不同状态下柔性探测器的动力学模型，提出了一种视觉/雷达的小天体柔性探测器多模型导航方法，提高导航精度，进而解决柔性探测器在小天体表面附着等方面问题。
7.本发明公开的一种基于视觉雷达的小天体探测器柔性多模型导航方法，包括以下步骤：
8.步骤一、搭建柔性探测器物理模型：
9.将柔性材料与探测装置结合，设每个质量聚集体为一个节点，可以设置n个节点。其中柔性材料基本由硅橡胶与有机硅组成，根据材料构成不同，柔性材料有不同的柔变系数；每个探测装置由视觉相机、激光雷达组成，用于定位与导航。
10.步骤二、构建柔性探测器动力学模型：
11.根据柔性探测器的刚柔耦合特性，将柔性探测器节点之间的连接等效为弹性连接，得到探测器节点的非线性动力学模型。
12.本发明采用混合坐标系法，其中包括描述单个节点与着陆点处的距离、姿态倾斜角关系的聚集体平动系o-x
i-y
i-zi(i＝1,2,...,n)和描述整个柔性探测器与着陆点处地面的距离、姿态倾斜角关系的柔性探测器平动系o
t-x
t-y
t-z
t
，其中，i表示为柔性探测器上第i个探测器节点。
13.柔性探测器节点绕x轴转过的角度为θ，则其转动角速度与角加速度可分别表示为和记聚集体平动系原点在柔性探测器平动系下的位置矢量为r0,记聚集体平动系到柔性探测器平动系的空间转换矩阵为a。可得柔性探测器中的第i个探测器节点在柔性探测器平动系下的位置矢量ri为：
14.ri＝r0 a(ri nq)
15.其中，ri为节点i在聚集体平动系下未发生运动前的位置矢量，nq为节点i在聚集体平动系下的弹性变形。从而，节点i在柔性探测器平动系下的速度与加速度分别为：
[0016][0017][0018]
单元动能可表示为：
[0019][0020]
对单元动能进行变分，并在时间上进行积分，通过哈密顿原理得到柔性探测器节点间刚柔耦合非线性动力学模型为：
[0021][0022]
其中，c为阻尼等效矩阵，g为刚柔耦合引起的附加质量矩阵，m为节点质量矩阵，q
t
为刚性运动加速度引起的外载荷，q
π
为预应变产生的载荷列阵。
[0023]
柔性探测器节点在不同位置下，建立不同刚柔耦合非线性模型，得到各个柔性探测器节点的独立状态方程和独立观测方程。
[0024]
步骤三、通过视觉/雷达目标导航算法进行单模型探测器导航：
[0025]
通过柔性探测器搭载相机拍摄的图像进行解算得到目标的深度信息，雷达信息经过杂质过滤，将视觉信息与雷达数据进行数据融合，设计一种将融合后的像素坐标作为观测量的卡尔曼滤波，将融合结果输入到滤波器中，得到目标特征点的位置估计，原理如图4所示。
[0026]
其中，针对传统的视觉/雷达融合的目标导航方法，即将视觉与雷达解算的位置信
息直接作为卡尔曼滤波的观测量的方法，视觉解算结果误差为非高斯白噪声，导致滤波精度不高，所以提出了一种将像素坐标作为观测量的视觉/雷达信息融合的卡尔曼滤波器，以提高目标跟踪精度。
[0027]
将柔性探测器节点的状态变量选取为：x＝[x,y,z,v
x
,vy,vz,a
x
,ay,az]
t
[0028]
其中，x,y,z为跟踪目标的位置信息，v
x
,vy,vz为跟踪目标三个方向的速度分量，a
x
,ay,az为跟踪目标三个方向的加速度信息。
[0029]
则其目标状态方程为：x(k 1)＝φx(k) γw(k)
[0030]
其中φ为状态转移矩阵，g为过程噪声矩阵，w为过程噪声。观测方程为：
[0031]
其中，(u1,v1)与(u2,v2)分别为柔性探测器节点1与柔性探测器节点2搭载相机拍摄到目标天体部分特征点的像素坐标，(xr,yr,zr)为雷达测量到的小天体目标特征点的三维坐标。
[0032]
步骤四、运用最近临界法进行特征点关联：
[0033]
雷达和视觉方法采集的小天体表面的数据复杂，其中包括不同的特征点和干扰数据，本发明采用最近临界法进行特征点匹配和关联，提高探测器对小天体表面数据采集的准确度。
[0034]
其中，假设检测到小天体表面的特征点有c个类别ω1，ω2，
…
，ωc，每类有ni(i＝1,2,
…
c)个标明类别的样本。则规定ωi类的判别函数为：
[0035][0036]
其中，i表示ωi类，k表示ωi类ni个样本中的第k个。
[0037]
最近邻法应用到小天体表面多目标导航系统中的过程如图5所示。假设存在三个小天体目标特征点，相应传感器在各时刻分别可以得到三组数据(a
x(t)
,a
y(t)
,a
z(t)
)、(b
x(t)
,b
y(t)
,b
z(t)
)与(c
x(t)
,c
y(t)
,c
z(t)
)，将属于同一类的数据串联起来。
[0038]
步骤五、进行柔性探测器多模型轨迹导航：
[0039]
步骤三运用的是单一模型的导航算法，现在运用多模型算法，根据各个柔性探测器节点之间的位置参数不同，得到的不同动力学模型，将各种状态下的模型引入系统中，使各个运动模型尽可能准确的描述柔性探测器不同的运动状态，经过输入交互、状态滤波、模型概率更新和状态融合输出四个步骤，完成精确导航。多模型算法原理如图6所示。
[0040]
本发明采用多模型算法将柔性探测器的各种运动模型引入系统中，使各个运动模型尽可能准确的描述柔性探测器的运动状态。
[0041]
多模型算法主要包括输入交互、状态滤波、模型概率更新和状态融合四个步骤。通过输入交互，可以得到系统内每一个运动模型滤波器的初始状态以及初始状态协方差的估计值，将通过输入交互得到的估计值作为各自模型滤波器的输入值将通过输入交互得到的初始状态估计值与初始状态协方差估计值作为滤波器当前输入，估计出下一时刻的状态与状态协方差；经过滤波器更新后，可以得到第j个滤波器的新息和新息协方差，结合预测模型转移概率μ
ij
(k 1∣k)，可以得到当前时刻的模型概率估计；将求得的模型概率作为加权值，将状态估计与状态协方差估计的加权和作为整个算法的最终输出值：
[0042][0043][0044]
其中与分别为k时刻第j个模型的状态估计值与状态协方差估计值,u
ij
(k∣k 1)为预测模型i转移到模型j的转移概率,可由下式求得：
[0045][0046]
通过上述过程可以分别得到各个柔性探测器节点对目标小天体表面特征点的导航信息，将这些信息经过卡尔曼滤波器可实现对各轨迹的融合，得到最终的导航轨迹估计信息。
[0047]
与现有技术相比，本发明的优点和积极效果在于：
[0048]
1、本发明公开的一种视觉雷达的小天体探测器柔性多模型导航方法，将探测器柔性化，根据柔性材料的柔变系数不同，将柔性体内力、内力矩等效为弹簧、扭簧、阻尼器的作用力及力矩，建立各个探测器之间的不同动力学模型，将视觉与雷达信息进行融合，将融合后的数据作为卡尔曼滤波的输入之一，用于导航轨迹预测，提高导航精度。
[0049]
2、本发明公开的一种视觉雷达的小天体探测器柔性多模型导航方法，利用雷达和视觉得到的小天体特征点和各个探测器节点的位置信息，引入不同位置下不同的柔性探测器模型，各个探测器之间数据交互，得到不同的观测方程，通过滤波，提高对小天体表面特征点导航的准确性，完成探测器对小天体表面的着陆导航。
附图说明
[0050]
图1为本发明的总体流程图；
[0051]
图2为本发明柔性探测器物理模型俯视图；
[0052]
图3为本发明柔性探测器物理模型仰视图；
[0053]
图4为本本发明的视觉雷达单一模型跟踪原理示意图；
[0054]
图5为本发明航迹关联算法示意图；
[0055]
图6为本发明多模型算法示意图；
[0056]
图7为本发明仿真轨迹图；
[0057]
图8为本发明仿真位置误差图；
[0058]
图9为本发明多模型概率估计曲线。
具体实施方式
[0059]
为了能够更加清楚地理解本发明的目的、特征和优点，下面结合附图及实例对本发明做进一步说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用不同于在此描述的其他方式来实施，因此，本发明并不限于下面公开的具体实施例。
[0060]
一种基于视觉雷达的小天体探测器柔性多模型导航方法，包括以下步骤：
[0061]
步骤一、构建柔性探测器物理模型；
[0062]
步骤二、构建柔性探测器动力学模型；
[0063]
步骤三、通过视觉/雷达目标导航算法进行单模型柔性探测器导航；
[0064]
步骤四、运用最近临界法进行特征点关联；
[0065]
步骤五、进行多模型轨迹导航。
[0066]
具体的：
[0067]
步骤一、构建柔性探测器物理模型：
[0068]
本发明将柔性材料与探测装置结合，每个探测装置为一个节点，以4个探测节点为例。其中柔性材料基本由硅橡胶与有机硅组成，根据材料不同，柔性材料有不同的柔变系数；每个探测装置由视觉相机、激光雷达成，用于定位与导航，柔性探测器物理模型结构俯视和仰视图如图2和图3所示。
[0069]
步骤二、构建柔性探测器动力学模型：
[0070]
将探测器节点之间的连接等效为弹性连接，根据柔性探测器的刚柔耦合特性，得到探测器节点的非线性动力学模型。
[0071]
本发明采用混合坐标系法，其中包括描述单个节点与着陆点处的距离、姿态倾斜角关系的聚集体平动系o-x
i-y
i-zi(i＝1,2,...,n)和描述整个柔性探测器与着陆点处地面的距离、姿态倾斜角关系的探测器平动系o
t-x
t-y
t-z
t
，其中i为柔性探测器上不同的探测器节点。
[0072]
柔性探测器节点绕x轴转过的角度为θ，转动角速度与角加速度分别为和记聚集体平动系原点在柔性探测器平动系下的位置矢量为r0,记聚集体平动系到柔性探测器平动系的空间转换矩阵为a。可得柔性探测器中的第i个柔性探测器节点在柔性探测器平动系下的位置矢量ri为：
[0073]ri
＝r0 a(ri nq)
[0074]
其中，ri为节点i在聚集体平动系下未发生运动前的位置矢量，nq为节点i在聚集体平动系下的弹性变形。从而，节点i在柔性探测器平动系下的速度与加速度分别为：
[0075][0076][0077]
单元动能可表示为：
[0078]
[0079]
对单元动能进行变分，并在时间上进行积分，通过哈密顿原理得到柔性探测器节点间刚柔耦合非线性动力学模型为：
[0080][0081]
其中，c为阻尼等效矩阵，g为刚柔耦合引起的附加质量矩阵，m为节点质量矩阵，q
t
为刚性运动加速度引起的外载荷，q
π
为预应变产生的载荷列阵。根据刚柔耦合非线性模型，得到各个探测器节点的独立状态方程和独立观测模型。
[0082]
步骤三、通过视觉/雷达目标导航算法进行单模型柔性探测器导航：
[0083]
通过探测器搭载相机拍摄的图像进行解算得到目标的数据信息，雷达信息经过杂质过滤，将视觉信息与雷达数据进行数据融合，设计一种将融合后的像素坐标作为观测量的卡尔曼滤波，将融合结果输入到滤波器中，得到目标特征点的位置估计，原理如图4所示。
[0084]
其中，针对传统的视觉/雷达融合的目标导航方法，即将视觉与雷达解算的位置信息直接作为卡尔曼滤波的观测量的方法，视觉解算结果误差为非高斯白噪声，导致滤波精度不高，所以提出了一种将融合后的像素坐标作为观测量的视觉/雷达信息融合的卡尔曼滤波器，以提高目标跟踪精度。
[0085]
将柔性探测器节点的状态变量选取为：x＝[x,y,z,v
x
,vy,vz,a
x
,ay,az]
t
[0086]
其中，x,y,z为跟踪目标的位置信息，v
x
,vy,vz为跟踪目标三个方向的速度分量，a
x
,ay,az为跟踪目标三个方向的加速度信息。
[0087]
则其目标状态方程为：x(k 1)＝φx(k) γw(k)
[0088]
其中φ为状态转移矩阵，g为过程噪声矩阵，w为过程噪声，观测方程为：
[0089]
其中，(u1,v1)与(u2,v2)分别为探测器节点1与探测器节点2搭载相机拍摄到目标天体部分特征点的像素坐标，(xr,yr,zr)为雷达测量到的小天体目标特征点的三维坐标。
[0090]
步骤四、运用最近临界法进行特征点关联；
[0091]
雷达和视觉采集的小天体表面的数据信息复杂，其中包括不同的特征点和干扰数据，本发明采用最近临界法进行特征点匹配和关联，提高探测器对小天体表面数据采集的准确度。
[0092]
其中，假设检测到小天体表面有c个类别ω1，ω2，
…
，ωc，每类有ni(i＝1,2,
…
c)个标明类别的样本。则规定ωi类的判别函数为：
[0093][0094]
其中，i表示ωi类，k表示ωi类ni个样本中的第k个。
[0095]
最近邻法应用到小天体表面多目标导航系统中的过程如图5所示。假设存在三个小天体目标特征点，相应传感器在各时刻分别可以得到三组数据(a
x(t)
,a
y(t)
,a
z(t)
)、(b
x(t)
,b
y(t)
,b
z(t)
)与(c
x(t)
,c
y(t)
,c
z(t)
)，将属于同一类的数据串联起来。
[0096]
步骤五、进行多模型轨迹导航：
[0097]
步骤三运用的是单一模型的跟踪算法，现在运用多模型算法将柔性探测器节点不同受力情况下的运动模型引入系统中，使各个运动模型可以尽可能准确的描述各柔性探测器节点的运动状态，经过输入交互、状态滤波、模型概率更新和状态融合输出四个步骤，得到各探测器点的最终运动轨迹，算法详细原理如图6所示。
[0098]
采用多模型算法将柔性探测器节点的各种运动模型引入系统中，尽可能准确的描述柔性探测器节点的运动状态。
[0099]
多模型算法主要包括输入交互、状态滤波、模型概率更新和状态融合四个步骤。具体过程如下：
[0100]
输入交互：通过输入交互，可以得到系统内每一个运动模型滤波器的初始状态以及初始状态协方差的估计值，将通过输入交互得到的估计值作为各自模型滤波器的输入值。输入交互主要包括状态交互与状态协方差交互两个方面。
[0101]
状态交互：
[0102]
状态协方差交互：
[0103]
其中，与分别为k时刻第j个模型的状态估计值与状态协方差估计值，u
ij
(k∣k 1)为预测模型i转移到模型j的转移概率：
[0104][0105]
状态滤波：将通过输入交互得到的初始状态估计值与初始状态协方差估计值作为滤波器当前输入，估计出下一时刻的状态与状态协方差，其过程与单模型滤波过程相同。
[0106]
预测：
[0107]
滤波增益计算：
[0108]
状态更新：
[0109]
模型概率更新：经过滤波器更新后，可以得到第j个滤波器的新息和新息协方差，结合预测模型转移概率μ
ij
(k 1∣k)，可以得到当前时刻的模型概率估计。
[0110]
新息和新息协方差：
[0111]
模型概率更新：
[0112]
其中，π
ij
为模型i转移到模型j的转移概率。
[0113]
状态融合输出：将求得的模型概率作为加权值，将状态估计与状态协方差估计的加权和作为整个算法的最终输出值。
[0114]
通过上述过程可以分别得到各个探测器节点对小天体目标表面特征点的导航信息，控制柔性探测器整体导航附着小天体。
[0115]
为验证本方法的有效性，以三个柔性探测器节点为例，对小天体表面一个特征点进行导航估计，通过模型受力分析，设定模型一(m1)为匀速运动，模型二(m2)为顺时针运动，模型三(m3)为逆时针运动。
[0116]
其中模型一：
[0117]
模型二与模型三：
[0118]
假设各柔性探测器节点的运动仿真条件如表1所示，小天体目标在1～150s、251～400s做匀速直线运动，151～250s以的速度做顺时针运动，以的速度401～500s做逆时针运动。在上述仿真条件下进行100次蒙特卡洛仿真，得到柔性探测器对小天体目标特征点的导航曲线、位置误差、模型概率估计曲线等，仿真结果如图7至9所示。
[0119]
探测器节点运动仿真条件：表1
[0120]
从图7可以看出，该方法可以实现对小天体表面特征点的有效导航，从图8的位置误差图可以看出，该算法在三个方向的误差均小于1.5米，根据图9的模型概率更新来看，本发明提出的方法可以对模型改变做出及时响应，精度更高，可靠性更强。
[0121]
以上所述，仅为本发明的较佳实例，并非是对本发明做其他限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例应用于其它领域，但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所做的任何修改、等同变化与改型，仍属于本发明技术方案的保护范围。