初始轨道确定吉布斯误差分析范文

《中国惯性技术学报》2016年第二期

摘要：

现代战争中的空间武器平台捕获敌方的空间目标，确定其运行轨道后实施军事打击，所以对轨道确定参数的误差进行分析对于精确打击具有非常重要的意义。为了准确分析吉布斯三位置矢量定轨法的轨道根数的确定精度，根据吉布斯方法的轨道确定模型，采用向量求导的方法，详细推导了轨道根数对于观测位置误差的灵敏度矩阵，并给出了轨道根数误差与观测位置误差之间的关系，明确初始轨道的确定精度及误差传播规律。最后对吉布斯方法及其误差分析进行仿真分析验证，验证了方法的正确性。

关键词：

初始轨道确定；吉布斯方法；轨道根数；灵敏度矩阵

在轨道力学中，对空间目标进行轨道确定包括两个过程：利用短弧段观测数据的初始轨道确定以及长弧段下的精密定轨[1]。初始轨道确定一般都是采用二体模型，经典的初始轨道确定方法都是利用测角数据来进行计算的，主要有Laplace法和Gauss法两种。长期以来，国内外的众多学者都针对以上两种方法进行了更加深入的研究，并提出了多种初始轨道确定的方法，解决了一些使用传统方法会出现不收敛或不稳定等问题。文献[2]基于单站单圈测量数据研究了低轨卫星初轨确定的方法，对遗传算法和最小二乘法在不同观测条件下确定的初轨信息以及轨道预测信息进行了分析和比较，确定了两种方法的优缺点和各自的适用条件。文献[3]采用最小一乘方法建立了一种初轨计算的稳健方法，将初轨计算问题转换为线性规划问题求解，并通过bootstrap方法给出估计精度，数值计算结果表明该方法稳健有效，并具有较高的崩溃点。文献[4]将采用角度和角速率数据进行初始轨道确定与只采用角度数据进行了对比，采用角度和角速率数据有着更好的平稳性和可靠性。文献[5]提出了一种采用单时刻GPS测量的方法来对地球静止卫星进行初始轨道的确定，并采用EKF对该算法进行了验证。文献[6]利用遗传算法来对极短弧的光学观测进行初始轨道确定，并利用蒙特卡洛仿真来验证算法的可靠性。文献[7]将天基角度测量与遗传算法相结合来确定初始轨道，得到了满意的结果，并且多种群遗传算法（MPGA）可以有效解决角度测量方法的迭代过程中出现的问题。文献[8]中根据开普勒运动定律推导出确定斜距的方法，5个数据一组进行初始轨道的确定。除此之外，还可以利用其他方法来确定初始轨道，如利用高斯混合模型求解[9]。

以上是对初始轨道确定方法研究现状的概述，对轨道确定方法进行误差分析的相关文献相对较少。文献[10]中分析了星载GPS接收机的定位误差对卫星轨道根数的影响，只推导了长半轴，轨道倾角和升交点赤经的误差与位置、速度误差的关系式。本文将主要对吉布斯三位置矢量定轨法进行误差分析。首先给出轨道根数的计算方法和吉布斯三位置矢量定轨法。利用3个位置矢量就能确定轨道根数，但目前没有文献给出位置矢量本身的测量误差、夹角的大小等对于6个轨道根数确定误差的影响程度。在整个轨道确定的过程中，初始轨道确定的精度大多数情况下无法满足精度要求，需要进行精密定轨，而初始轨道确定的精度往往会影响精密定轨初值的选取，高精度的初始轨道的确定可以减少整个轨道确定的时间。所以本文将对吉布斯方法进行误差分析，推导出6个轨道根数的灵敏度矩阵，明确误差传播特性，并对其进行仿真分析与验证。

1轨道根数与初始轨道确定

1.1轨道根数轨道根数是用来描述天体在其轨道运行状态的一组参数。6个轨道根数如图1所示，其中，h为比角动量的模。

1.2初始轨道确定——吉布斯三位置矢量定轨法如果测量出空间目标轨道上P1、P2、P3位置的位置矢量1r、2r和3r后，便可用吉布斯三位置矢量定轨法来计算6个轨道根数。如图2所示。

2初始轨道确定方法的误差分析

对吉布斯三位置矢量定轨法进行误差分析，要根据6个轨道根数的计算公式写成关于3个位置矢量1r、2r和3r的微分表达式，求取关于1r、2r和3r的偏导，转换成矩阵的形式，得到6个轨道根数a、i、、e、ω、θ的灵敏度矩阵，并且h作为重要的中间变量同样需要求取灵敏度矩阵。在求解过程中，主要利用向量对向量求偏导的相关计算公式以及如下几个常用的求取偏导公式。标量对位置矢量的灵敏度矩阵均为19维矩阵，假设各个位置矢量测量的不确定度相等，均为，并假设各次测量均为独立测量，再结合不确定度的合成公式，即当所求取轨道根数的灵敏度矩阵为A时，那么轨道根数的不确定度。

3初始轨道确定方法的仿真分析

利用STK软件绘制轨道，设定轨道参数如下对吉布斯三位置矢量定轨法进行Matlab编程，经过计算可以得到6个轨道根数，长半轴a8002.04km，轨道倾角i59.9992，偏心率e0.100132，升交点赤39.9932，近地点幅角30.1203，真近点角49.8836（相对于位置矢量2r）。把以上建立的运行轨道模型中的位置信息当作基准数据，在基准数据中加入0.1km的位置随机误差，多次取值，求取6个轨道根数的不确定度，并在此基础之上逐步增大选取的3个位置矢量之间的夹角，最大到120°左右，观察6个轨道根数不确定度的变化趋势，绘制出变化曲线。令位置不确定度0.1km，根据前面得到的6个轨道根数的灵敏度矩阵，计算6个轨道根数的不确定度，与由仿真数据计算得到的不确定度进行比较，仿真结果如图3~8所示。从上述仿真结果可以看到，利用吉布斯三位置矢量定轨法计算得到的轨道根数与预设的轨道根数接近，由实际数据计算得到的轨道根数的不确定度与由灵敏度矩阵计算得到的轨道根数的不确定度，两者的结果十分接近，验证了前面所求取的轨道根数的灵敏度矩阵的准确性。从轨道根数的不确定度的变化趋势可以看到，增大选取位置矢量之间的夹角，6个轨道根数的不确定度在减小，所以尽量增大选取的3个位置矢量之间的夹角有利于减小轨道根数的误差。但实际确定初始轨道时，增大3个位置的夹角可能受到所需确定时间的限制，可根据具体所需的轨道确定精度和效率进行取舍。

4结论

针对吉布斯方法确定的6个轨道根数的误差，推导了6个轨道根数误差相对于位置矢量测量误差的灵敏度矩阵，用Matlab软件编程计算了6个轨道根数的不确定度；然后，利用STK软件绘制轨道，验证了空间目标初始轨道确定方法的有效性；最后，对选取位置矢量之间的夹角做出改变进行了仿真分析，绘制了轨道根数不确定度的变化曲线，将真实数据计算得到的不确定度与由灵敏度矩阵计算得到的不确定度进行比较，验证了前面所求取的灵敏度矩阵的准确性。由仿真结果可知，为减小轨道根数的误差，可以通过增大选取的3个位置矢量之间的夹角来实现。在实际的初始轨道确定中，可根据本文给出的6个轨道根数灵敏度矩阵可以计算出所求解的轨道的轨道根数误差。

作者：刘庆博 任顺清 单位：哈尔滨工业大学 空间控制与惯性技术研究中心