【技术实现步骤摘要】
本专利技术涉及一种空间目标运动估计方法，属于空间。
技术介绍
1、空间目标运动估计是当前近地空间、地月空间以及深远空间航天活动中的重要任务。空间目标运动估计可以为航天器以及空间碎片、行星、小天体等提供位置信息，而准确快速得估计位置信息对于航天活动中轨迹规划、制导、控制至关重要，在空间碎片清理、地火转移中途修正、地月空间周期轨道保持等航天任务中广泛应用。通常，空间目标运动估计可以通过卡尔曼滤波等方法进行求解。卡尔曼滤波算法包括两部分，一是状态更新，二是测量更新。在状态更新过程中，通常需要进行大量的轨道递推，轨道递推占据大量的计算资源，使得现有的卡尔曼滤波等算法在空间目标运动估计中难以在线应用。针对现有空间目标运动估计方法中，轨道递推耗时巨大的问题，有必要通过训练代理模型，预测轨道递推结果，进而减少算法中轨道递推次数，从而大幅减少计算耗时，实现空间目标实时运动估计。
技术实现思路
1、本专利技术主要目的是提供一种基于代理模型的空间目标运动估计方法，构建用于轨道递推的代理模型，使用代理模型预测的轨道状态代替轨道递推得到的轨道状态，进而预测空间目标状态均值与协方差，降低空间目标运动估计耗时，提升空间目标运动估计精度。本专利技术具有定轨精度高，估计误差小、计算耗时短的优点。
2、本专利技术的目的是通过以下技术方案实现的。
3、本专利技术公开的一种基于代理模型的空间目标运动估计方法，建立空间目标运动估计的动力学方程方程与观测方程，设定空间目标运动估计初始时刻，设定空间目
4、本专利技术公开的一种基于代理模型的空间目标运动估计方法，包括如下步骤：
5、步骤一：建立空间目标运动估计的动力学方程方程与观测方程，设定空间目标运动估计初始时刻t0，设定空间目标初始时刻运动状态估计值设定空间目标初始时刻估计协方差矩阵设定观测间隔dt，设定观测次数n，设定观测噪声协方差矩阵r，将所建立的空间目标运动估计的动力学方程与观测方程转化成离散形式。设定粒子采样数目n，设定迭代参数k←1。
6、建立空间目标运动估计的动力学方程与观测方程如下：
7、
8、其中x为空间目标在时刻t的运动状态，表示空间目标运动的动力学方程，z＝h(x)+v表示空间目标运动估计的观测方程，z为观测矢量，h(x)为测量模型，v为观测噪声。
9、设定空间目标运动估计初始时刻t0，设定空间目标初始时刻运动状态估计值设定空间目标初始时刻估计协方差矩阵设定观测间隔dt，设定观测次数n，设定观测噪声协方差矩阵r。
10、将所建立的空间目标运动估计的动力学方程与观测方程转化成离散形式，有如下表达式：
11、
12、其中下标k与k-1分别表示第k个观测时刻与第k-1个观测时刻对应的状态，xk＝f(xk-1)为离散的空间目标运动估计的动力学方程，表达式为：
13、
14、设定粒子采样数目n，设定迭代参数k←1。
15、步骤二：定义tk表示当前时刻，tk-1表示当前时刻的上一时刻，更新时刻tk＝tk-1+dt，构建用于空间目标运动估计的轨道递推的代理模型fhdmr(xk-1)。
16、定义tk表示当前时刻，tk-1表示当前时刻的上一时刻，更新时刻tk＝tk-1+dt，构建用于空间目标运动估计的轨道递推的代理模型如式(4)所示：
17、
18、其中fhdmr(xk-1)表示空间目标运动估计的动力学方程离散形式的代理模型，与表示输入矢量xk-1的第i个元素和第j个元素，表示代理模型fhdmr(xk-1)的0阶项，表示代理模型fhdmr(xk-1)的1阶项，表示代理模型fhdmr(xk-1)的2阶项。其中代理模型fhdmr(xk-1)的1阶项与代理模型fhdmr(xk-1)的2阶项均为多项式，即有：
19、
20、
21、其中ai表示代理模型fhdmr(xk-1)的1阶项的系数，bij表示代理模型fhdmr(xk-1)的2阶项的系数。
22、步骤三：以运动状态估计值为参考值，根据步骤二构建的轨道递推的代理模型fhdmr(xk-1)，计算用于空间目标运动估计的轨道递推的代理模型fhdmr(xk-1)的系数。
23、以运动状态估计值为参考值，构造代理模型fhdmr(xk-1)的0阶项代理模型fhdmr(xk-1)的1阶项与代理模型fhdmr(xk-1)的2阶项如下：
24、
25、
26、
27、其中表示将的第i个元素替换为递推得到的结果，表示将的第i个元素和第j个元素分别替换为和递推得到的结果。通过式(7)-(9)得用于空间目标运动估计的轨道递推的代理模型的系数如下：
28、
29、
30、式(10)与(11)即为用于空间目标运动估计的轨道递推的代理模型的系数。
31、步骤四：根据运动状态估计值与估计协方差矩阵随机生成用于轨道预测的粒子，利用代理模型预测的轨道状态代替轨道递推得到的轨道状态，根据代理模型预测的轨道状态计算一步预测的均值与协方差。
32、根据运动状态估计值与估计协方差矩阵随机生成n个粒子，记为其中下标p表示第p个粒子。利用构建的用于空间目标运动估计的轨道递推的代理模型fhdmr(xk-1)计算第p个粒子对应的输出收集所有粒子的输出，按照式(12)计算一步预测的均值与协方差：
33、
34、
35、其中为一步预测的均值，为一步预测的协方差矩阵。由于使用大量粒子进行一步预测，一步预测的误差小，空间目标运动估计精度高；另外使用构建的用于空间目标运动估计的轨道递推的代理模型fhdmr(xk-1)计算粒子对应的输出，避免大量轨道递推计算，减少空间目标运动估计计算量。
36、步骤五：根据步骤四中一步预测的均值与协方差、步本文档来自技高网...
【技术保护点】
1.一种基于代理模型的空间目标运动估计方法，其特征在于：包括如下步骤，
2.如权利要求1所述的一种基于代理模型的空间目标运动估计方法，其特征在于：还包括步骤八，基于得到的空间目标运动估计值与估计协方差矩阵，对观测航天器进行规划，观测航天器执行规划得到的制导、控制相关指令，实现对于空间目标的抵近，进而实现对空间目标进行绕飞观测。
3.如权利要求1或2所述的一种基于代理模型的空间目标运动估计方法，其特征在于：步骤一实现方法为，
4.如权利要求3所述的一种基于代理模型的空间目标运动估计方法，其特征在于：步骤二实现方法为，
5.如权利要求4所述的一种基于代理模型的空间目标运动估计方法，其特征在于：步骤三实现方法为，
6.如权利要求5所述的一种基于代理模型的空间目标运动估计方法，其特征在于：步骤四实现方法为，
7.如权利要求6所述的一种基于代理模型的空间目标运动估计方法，其特征在于：步骤五实现方法为，
8.如权利要求7所述的一种基于代理模型的空间目标运动估计方法，其特征在于：步骤六实现方法为，
【技术特征摘要】
1.一种基于代理模型的空间目标运动估计方法，其特征在于：包括如下步骤，
2.如权利要求1所述的一种基于代理模型的空间目标运动估计方法，其特征在于：还包括步骤八，基于得到的空间目标运动估计值与估计协方差矩阵，对观测航天器进行规划，观测航天器执行规划得到的制导、控制相关指令，实现对于空间目标的抵近，进而实现对空间目标进行绕飞观测。
3.如权利要求1或2所述的一种基于代理模型的空间目标运动估计方法，其特征在于：步骤一实现方法为，
4.如权利要求3所述的一种...
【专利技术属性】
技术研发人员：乔栋，周星宇，李翔宇，郑德智，秦同，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。