Issue |
JNWPU
Volume 41, Number 4, August 2023
|
|
---|---|---|
Page(s) | 802 - 811 | |
DOI | https://doi.org/10.1051/jnwpu/20234140802 | |
Published online | 08 December 2023 |
A random error estimation method for satellites being launched into orbit based on ellipsoid theory
基于椭球理论的卫星入轨位置随机误差估计方法
1
School of Computer Science, China University of Geosciences(Wuhan), Wuhan 430074, China
2
Hubei Key Laboratory of Intelligent Geo-Information Processing, China University of Geosciences(Wuhan), Wuhan 430074, China
Received:
16
September
2022
Accurate and fast random error estimation is essential to improve satellite control. In this paper, the random error of the satellite being launched into orbit is explored. According to the uncertainty matrix of six satellite orbit parameters associated with the project, the uncertainty of the orbit position in the three axes of the geocentric inertial frame of reference is calculated via the propagation rate of the covariance matrix. Then, the random error probability density of the satellite position is constructed for normal distribution. The surface of the error with the same probability density is approximated as an ellipsoid, the size is deduced along with the ellipsoid's shape and direction, and the satellite's spatial position probability is obtained. Monte Carlo method and STK software were employed to simulate low-earth orbit micro-nano satellite launching into the same orbital position. It is confirmed that the results obtained via the theoretical model are in line with practice. Therefore, the ellipsoid theory proposed in this paper can be used to describe the satellite random position error.
摘要
准确快捷地估计卫星位置随机误差有利于更好地控制卫星。针对卫星入轨位置的随机误差, 分析了入轨位置随机误差的概率密度函数, 根据卫星入轨六根数不确定性矩阵, 利用协方差阵传播率确定地心惯性坐标系下入轨位置在3个坐标轴方向的不确定性, 把入轨位置等概率密度曲面近似为椭球, 推导出椭球的形状、大小和方向, 进而计算卫星入轨在一定空间范围内的概率。利用蒙特卡罗方法, 借助STK软件模拟低轨微纳卫星在随机因素影响下发射入轨, 仿真实际入轨位置分布, 并对比椭球理论模型下入轨位置的概率分布。仿真结果表明椭球理论误差模型分析结果与模拟实验结果基本一致, 说明该椭球理论模型可以用来估计卫星随机位置误差。
Key words: satellite orbit / initial position / random error / error ellipsoid / probability density surface
关键字 : 卫星轨道 / 初始位置 / 随机误差 / 误差椭球 / 概率密度曲面
© 2023 Journal of Northwestern Polytechnical University. All rights reserved.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
卫星入轨位置精准度受多种因素的影响, 不可避免地存在系统误差和随机误差。系统误差可以在设计或实际操作过程中予以考虑并削减, 而随机误差是由偶然因素引起的, 难以预测和控制。随机误差亦会导致卫星入轨空间位置的不确定, 并在卫星组网和运行过程中传播, 对工作效能产生影响。工程实际中对于高价值的重要航天器, 一般是通过发动机点火等方式对轨道进行校正和调整, 但对于小卫星和微纳卫星, 由于其体积小携带有限, 进行自身轨道修正的代价较高或远超其设计发射成本, 况且大规模星群也没法进行广泛的轨道修正。但不同用途的卫星对位姿精度要求不同, 如果能准确快速地估计其位置偏差, 可以让工作者做出正确的决策。
当前已经有很多学者对卫星定位精度及运行误差进行了相关研究, 并取得了一定成果。如针对卫星的定位精度, Knogl等[1]利用LEO卫星通信信道对GEO卫星进行精确定位。Vighnesam等[2]研究了通过位置差确定IRS运行轨道的系统方法。Zhou等[3]针对北斗卫星定轨提出了一种基于反单点定位的近实时无推力建模方法。Li等[4]实现了利用卫星激光测距分析北斗3号卫星的轨道精度。Larki等[5]研究了基于梯度法确定卫星轨道的方法。白显宗[6]基于历史轨道数据研究了空间目标轨道预报误差与碰撞概率问题。贺若飞等[7]提出一种基于蒙特卡罗卡尔曼滤波的目标定位算法, 并应用于无人机目标定位中, 取得较好效果。阎海峰等[8]针对伪卫星空中基站位置不易精确确定的问题, 建立了组合导航系统和非线性数学模型, 提出了一种抗差自适应高斯混合Sigma点粒子滤波算法。文献[9]基于多项式区间参数分析法, 针对探测器着陆过程动态特性, 提出基于非线性有限元建模的探测器着陆区间分析流程, 计算得到探测器着陆动响应上界和下界。文献[10]对分布式InSAR目标三维定位的空间几何关系, 给出了定位精度指标与系统参数精度指标之间关系的解析表达式。文献[11]针对三星编队系统, 将卫星位置误差源分为体现误差相关性的整体绝对位置误差和体现误差非相关性的星间相对位置误差, 从而对存在相关性误差的卫星编队系统进行定位精度分析。在运行误差分析中, 文献[12]对影响遥感卫星覆盖的不确定性因素进行归一化量化, 对具有不确定性因素影响的指标进行求解, 得到具有统计特性的指标结果。文献[13]提出一种线性回归近似替代不确定性, 将不确定性参数对于效能影响降维为一个线性组合系数进行求解。文献[14]给出了降维表示方法中功率谱密度函数的解析随机误差, 并通过数值算例加以证明。文献[15]提出了一种用神经网络替代不确定性的分析方法, 通过人工神经网络在拟合回归分析上的非线性特性,研究分布式卫星系统概念设计中的不确定性参数。文献[16]引入模拟误差场对不确定进行描述, 从而获得比较精确的估计数值。
综上可知, 目前针对卫星轨道位置精度或误差描述主要有解析法、数值法或通过第三方测量获得; 研究对象主要是中、高轨单个高价值卫星或低轨少量卫星系统。相关研究中的解析法是将某时刻卫星位置用轨道根数的解析函数表示, 为便于建模学者一般采用简化模型或对不确定性因素做近似处理。该方法在卫星数量不多的情况下适用性较好, 但影响计算精度, 也不适用于大规模星群。数值法在抽样数量足够多的前提下, 可实现对样本的无偏估计, 但在卫星规模较大情况下, 又存在计算效率低、耗时长等缺点。通过第三方介入测量, 又存在过程复杂、成本高等问题。显然上述方法不能满足近几年大量卫星发射入轨, 对位置误差多、精、快的评估要求。
本文结合解析法和数值法, 采用概率论的知识和误差椭球理论进行研究。首先针对卫星发射入轨位置随机误差, 依据卫星入轨六根数不确定性矩阵, 利用协方差阵传播率确定地心惯性坐标系下卫星入轨位置在3个坐标轴方向的位置不确定性, 构建了描述卫星入轨位置随机误差的一种椭球理论, 给出卫星入轨在一定空间范围内的概率计算公式。然后利用STK软件和蒙特卡罗方法, 模拟随机因素影响下低轨微纳卫星入轨位置, 验证误差理论模型与实际工况的一致性。实验结果表明本文构建的误差椭球理论可以很好地用于估计入轨位置随机误差, 尤其是估计大规模星群位置误差,本文方法具有高效性, 很好地克服了已有方法的局限性。
1 卫星发射入轨位置表达式
卫星发射入轨位置通常用轨道六根数表示。在航天工程中, 卫星六根数期望值根据效能要求精准确定, 而卫星实际入轨位置的离散程度通常可通过六根数的不确定性矩阵表示。
1.1 入轨位置不确定性矩阵
假设卫星入轨位置六根数用向量X表示, 即:
向量X中Xa, Xe, Xi, XE, Xω, Xf分别表示卫星轨道的半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角和卫星入轨位置的真近点角。
在卫星发射过程中考虑随机因素的影响, 其入轨位置为随机变量。设其轨道参数向量X的数学期望为
卫星轨道位置六根数向量X的方差用DXX表示, 即
DXX的主对角线上的元素为轨道六根数随机变量Xa, Xe, Xi, XE, Xω, Xf的方差, 非主对角线中的元素是描述2个轨道参数之间相关性的协方差。由于在卫星轨道六根数中, 各参数相互独立, 即协方差为零, 故向量X的方差阵可进一步写为
由于卫星入轨位置受许多相互独立的随机因素的影响, 且每个因素所产生的影响都很小, 由中心极限定理, 可假设入轨位置六根数服从正态分布。
1.2 卫星入轨位置随机误差表示
如果卫星在笛卡尔坐标系中希望入轨位置为, 实际入轨位置为(如图 1所示), 则实际入轨位置的偏差r可表示为:
式中, rx, ry, rz为空间位置误差r在x, y, z轴上的分量。显然
所以
式中, E(r2)是卫星希望入轨位置点P误差平方的理论平均值, 并定义为P点的点位方差, 记为σP2, 则
式中: σP2的大小与坐标系的选择无关; σP称为入轨位置误差。
图1 卫星入轨点位置误差示意图 |
2 卫星入轨位置误差模型构建
卫星入轨位置误差σP可以用来评定入轨位置精度, 但不能显示在哪个方向上的位差, 在讨论卫星或星群的空间碰撞安全、对地覆盖或其他效能时, 需要考虑卫星所处轨道位置在哪个方向上的偏差, 为此, 需要确定卫星在三维空间内的位置矢量差。卫星入轨位置等概率密度点的集合可形成一闭合曲面, 该曲面能把各个方向上的误差表示出来(见图 2)。但该曲面不是一种典型的曲面, 作图也不方便, 但其形状与椭球很相似, 为便于分析, 以椭球近似, 称之为等概率密度误差椭球(见图 3), 简称误差椭球。
图2 等概率密度点曲面示意图 |
图3 误差椭球示意图 |
2.1 误差椭球模型函数
在随机误差下, 卫星入轨位置在地心惯性坐标系下呈三维正态分布, 其联合分布密度函数为
式中:Dxx为地心惯性坐标系下卫星入轨位置的不确定性矩阵;|Dxx|为不确定性矩阵的行列式。Dxx的表达式为
式中: σi2(i=x, y, z)为坐标轴向的方差; σij(i, j=x, y, z, i≠j)为各坐标轴之间的协方差。
2.2 不确定性矩阵Dxx的确定
已知轨道六根数的协方差阵DXX, 根据协方差传播率[17]可以计算不确定性矩阵Dxx
式中, 矩阵A为一个3×6的矩阵,可用(11)式表示。
式中: x, y, z是在地心惯性坐标系下的卫星位置表达式; a, e, i, Ω, ω, f分别表示卫星轨道六根数; 是(1)式所示的轨道六根数的期望值。也就是矩阵A中各元素分别是卫星位置x, y, z的表达式对轨道六根数的偏导,并在轨道六根数期望值处的值。
2.3 误差椭球轴长的确定
根据联合分布密度函数(8)式, 该三维正态分布空间内概率密度相同点的表达式为
令x-rx=x′, y-ry=y′, z-rz=z′, 则(12)式可写为
(13) 式是地心惯性坐标系下一个相似椭球族表达式, k为放大因子,给定一个k值, 就得到一个椭球。为方便分析, 进行坐标转换, 求出其在主轴坐标系OUVW(如图 4所示)内的表达式。
图 4中椭球族的主轴为OU, OV, OW。因为Dxx为实对称矩阵, 则存在正交矩阵M, 使得
对(14)式两边同时取逆, 得到
进一步变换形式, 得到
令
得到
展开得到
(19) 式为在主轴坐标系下的误差椭球方程, 其半轴长平方分别为k2λ1, k2λ2, k2λ3, 这里k为前面定义的放大因子, λ1, λ2, λ3为协方差阵Dxx的特征值。已知协方差阵, 给定k值, 即可确定相应的椭球轴长。
图4 地心坐标与主轴坐标示意图 |
2.4 误差椭球轴向的确定
(14) 式中的正交矩阵M=(aij)3×3即是从地心惯性坐标系OXYZ转换到主轴坐标系OUVW的旋转矩阵, 由此可以确定椭球3个轴分别与Z轴的夹角。以轴OU为例, 如图 5所示, OU轴与OZ轴的夹角为α(0°≤α < 180°), OU轴在OXY平面的投影与OX轴的夹角为β(0°≤β < 360°), 由于M为旋转矩阵, 则可以直接根据旋转矩阵的元素计算相应角度, 并转换到α, β规定的范围内, 即
β可由a21, a31通过Matlab中双变量反正切函数确定。
同理, 可以计算OV和OW分别与OZ轴的夹角, 及其在XOY平面的投影与OX轴的夹角。
图5 椭球轴与地心惯性坐标轴的夹角示意图 |
3 卫星入轨位置在一定误差范围内的概率
根据概率密度函数(8)式和误差椭球方程(19), 卫星入轨在误差椭球内的概率可写为[18]
令
代入(21)式得到
(23) 式中Ck为半径等于的圆球, 目标点落入椭球Bk内的概率就相当于落入圆球Ck内的概率。进一步计算可以近似得到卫星入轨到一定误差范围内的概率
(24) 式中k仍为上文提到的椭球放大因子。在协方差阵确定的条件下, 误差椭球3个轴长, , 的比例是确定的, k值不同, 3个轴长度的放大比例不同, 卫星入轨在椭球内的概率就不同, 放大因子k与卫星入轨位置在椭球内的概率P的关系如图 6所示。
从图 6可以看出, 在随机误差下卫星入轨位置基本都在4倍的椭球轴长范围内。当然, 卫星入轨位置在不同的误差范围内的概率亦可得到。
图6 P与k的关系曲线图 |
4 算例分析及仿真验证
计划一颗离地526 km的太阳同步圆形轨道卫星发射升空, 在真近点角60°入轨, 分析其在随机因素影响下卫星入轨位置的随机误差。
4.1 确定相关参数
根据已知要求和给定条件, 确定卫星相关参数。
卫星轨道位置采用六根数表示, 地球半径按6 378.14 km计算, 则其轨道及位置参数如下:
轨道长半轴a=6 904.14 km; 轨道偏心率e=0;轨道倾角i=97.5°; 轨道近地点幅角ω=0°; 轨道升交点赤经Ω=0°; 卫星真近点角f=60°。
根据卫星研究所提供的经验数值, 依据(3)式确定轨道六根数的协方差阵DXX为
卫星入轨位置在地心惯性坐标系下的位置表达式
将轨道参数期望值代入(26)式得到地心惯性坐标系下期望入轨位置
根据(11)式, 对(26)式中x, y, z表达式分别对轨道六根数求偏导, 并计算在近似值(期望值)的结果, 得到偏导数矩阵A
根据(10)式得到地心惯性坐标系下的不确定性矩阵
将Dxx和rx, ry, rz代入(8)式即可得到该卫星入轨位置误差椭球函数。
4.2 误差椭球分析结果
对卫星入轨位置误差椭球函数进行计算, 得到正交矩阵M和对应的对角阵
据此确定卫星入轨位置误差矢量。当放大因子k取1时计算结果见表 1, 当k取2.8时(95%的概率落在椭球范围内)误差椭球如图 7所示。
进一步地, 利用误差椭球理论分析卫星在不同真近点角处的误差椭球、轴长、轴向分别如图 8~10所示。
椭球分析结果
图7 卫星入轨位置椭球仿真结果 |
图8 不同真近点角处误差椭球示意图 |
图9 不同真近点角下误差椭球三根轴长示意图 |
图10 不同真近点角处误差椭球三根轴轴向示意图 |
4.3 蒙特卡罗仿真结果
采用同样的轨道六根数协方差阵(15)式, 通过STK软件模拟该卫星入轨位置仿真10 000次, 则入轨位置分布如图 11所示。
对10 000次入轨位置求平均值, 得到入轨位置x, y, z的平均值(3 452.056 888, -780.432 294 9, 5 927.999 303), 根据(7)式可得到入轨位置点位方差为15.191 1 km2, 入轨位置中误差为3.897 6 km, 卫星入轨在一定误差范围内的概率用频度表示, 其概率密度直方图如图 12所示。
对蒙特卡罗仿真结果通过主成分分析, 得到卫星分布椭球3个轴长及轴向, 如表 2所示。
图11 蒙特卡罗仿真卫星入轨位置 |
图12 在一定误差范围内出现的概率直方图 |
蒙特卡罗结果主成分分析
4.4 2种方法分析结果比较
从表 1和表 2的结果可以看出, 误差椭球分析和蒙特卡罗仿真结果非常接近。
为进一步验证误差椭球理论的合理性, 在此采用提出的椭球模型计算卫星落入的概率, 并与蒙特卡罗实验仿真方法比较, 在特定k值时, 2种分析方法呈现的概率结果如表 3所示。
在放大因子k从0~4整个区间, 误差椭球理论估算的卫星入轨位置概率和模拟仿真的结果比较, 如图 13所示。
从图中可以看出, 2个结果非常接近, 且误差椭球计算概率略低于蒙特卡罗仿真结果, 说明该误差椭球理论模型可以可靠地估计卫星在随机误差下的入轨位置。
卫星落在一定误差范围内的概率
图13 2种方法分析卫星入轨误差概率曲线比较图 |
5 结论
本文构建了描述卫星入轨位置随机误差的一种椭球理论模型并验证了其合理性。根据随机误差下入轨位置六根数的不确定性矩阵, 把入轨位置随机误差等概率密度曲面近似为椭球, 确定椭球三根轴的长度及轴向, 并给出卫星入轨在一定误差范围内概率的计算方法。利用蒙特卡罗方法, 借助STK软件模拟了一颗低轨微纳卫星在随机因素影响下发射入轨的实验, 得到实验卫星入轨位置实际分布, 并与误差椭球模型下入轨位置的概率分布对比, 结果表明椭球模型分析结果与仿真结果一致, 从而说明误差椭球理论可以用于估计入轨位置随机误差, 尤其是高效估计大规模星群位置误差, 对大规模星群设计及发射具有较好参考价值。
精确地估计卫星入轨位置随机误差, 对于计算卫星运行中的碰撞安全、卫星对地覆盖效果等都有很重要的意义, 后期可以结合入轨位置误差椭球理论及误差传递特性进一步研究任意时刻卫星运行轨道的误差及对效能的影响。
References
- KNOGL J S, HENKEL P, GUNTHER C. Precise positioning of a geostationary data relay using LEO satellites[C]//53rd International Symposium Elmar, 2011: 325–328 [Google Scholar]
- VIGHNESAM NV, SONNEY A, SONI P K. CARTOSAT-1 orbit determination system and achieved accuracy during early phase[C]//2007 AAS/AIAA Astrodynamics Specialist Conference, 2007: 511–522 [Google Scholar]
- ZHOU P, DU L, LI X, et al. Near real-time BDS GEO satellite orbit determination and maneuver analysis with reversed point positioning[J]. Advances in Space Research, 2018, 63(5): 1781–1791 [Google Scholar]
- DING Li, XUE Cheng, JIA Xiaolin, et al. Precision analysis of BDS-3 satellite orbit by using slr data[C]//China Satellite Navigation Conference 2019 Proceedings, 2019: 389–399 [Google Scholar]
- LARKI M H A, MABOODI M, BOLANDI H. Satellite orbit determination based on gradient method[C]//Chinese Control and Decision Conference, 2012: 2554–2559 [Google Scholar]
- BAI Xianzong. Research on orbital prediction error and collision probability of space objects[D]. Changsha: National University of Defense Technology, 2013 (in Chinese) [Google Scholar]
- HE Ruofei, TIAN Xuetao, LIU Hongjuan, et al. A UAV target localization approach based on Monte-Carlo Kalman filter[J]. Journal of Northwestern Polytechnical University, 2017, 35(3): 435–441. [Article] (in Chinese) [Google Scholar]
- YAN Haifeng, WEI Wenhui, FENG Zhihua, et al. Research on high-performance algorithm for pseudolite air based positioning[J]. Journal of Northwestern Polytechnical University, 2015, 33(5): 763–769. [Article] (in Chinese) [Google Scholar]
- CHEN Zhaoyue, LIU Li, CHEN Shulin, et al. Interval uncertainty analysis of soft-landing dynamics of lunar lander[J]. Acta Armamentarii, 2019, 40(2): 442–448 [Google Scholar]
- GU Defeng, YI Dongyun, ZHU Jubo, et al. Closed-form solution of distributed InSAR geolocation and its precision analysis[J]. Acta Electronica Sinica, 2007(6): 1026–1031. [Article] (in Chinese) [Google Scholar]
- ZHENG Jinjun, CAO Zhigang. Positioning precision of three satellite constellation using relative locations in the presence of correlated satellite positioning errors[J]. Journal of Tsinghua University, 2004(10): 1398–1401. [Article] [Google Scholar]
- YI Shan. Effectiveness evaluation of remote sensing satellite system based on uncertainty[D]. Harbin: Harbin Institute of Technology, 2021 (in Chinese) [Google Scholar]
- GAO Chen, YANG Zhen, NIU Wenlong. Linear surrogate uncertainty analysis method for distributed satellite system[J]. Chinese Journal of Space Science, 2021, 41(6): 928–935. [Article] (in Chinese) [NASA ADS] [CrossRef] [Google Scholar]
- ZHANG Houle, WU Yongxin, WANG Juncheng. Analysis of stochastic errors and gaussianity of random vector process simulated by dimension-reduction representation method[J]. Applied Mathematical Modelling, 2022, 109: 283–303 [CrossRef] [Google Scholar]
- YANG Chen. Research on performance uncertainty analysis methods for distributed satellite system in conceptual design phase[D]. Beijing: China University of Science, 2020 (in Chinese) [Google Scholar]
- MAGGIONI V, SAPIANO M R P, ADLER R F. Estimating uncertainties in high-resolution satellite precipitation products: systematic or random error[J]. Journal of Hydrometeorology, 2016, 17(4): 1119–1129 [NASA ADS] [CrossRef] [Google Scholar]
- WANG Suihui. Error theory and basis of surveying adjustment[M]. Shanghai: Tongji University Press, 2015 (in Chinese) [Google Scholar]
- TANG Ning, WANF Shaoping, HONG Difeng. A calculation method of minimum distance between wellbore survey error ellipsoid[J]. Acta Petrolei Sinica, 2017, 38(11): 1320–1325. [Article] (in Chinese) [Google Scholar]
All Tables
All Figures
图1 卫星入轨点位置误差示意图 |
|
In the text |
图2 等概率密度点曲面示意图 |
|
In the text |
图3 误差椭球示意图 |
|
In the text |
图4 地心坐标与主轴坐标示意图 |
|
In the text |
图5 椭球轴与地心惯性坐标轴的夹角示意图 |
|
In the text |
图6 P与k的关系曲线图 |
|
In the text |
图7 卫星入轨位置椭球仿真结果 |
|
In the text |
图8 不同真近点角处误差椭球示意图 |
|
In the text |
图9 不同真近点角下误差椭球三根轴长示意图 |
|
In the text |
图10 不同真近点角处误差椭球三根轴轴向示意图 |
|
In the text |
图11 蒙特卡罗仿真卫星入轨位置 |
|
In the text |
图12 在一定误差范围内出现的概率直方图 |
|
In the text |
图13 2种方法分析卫星入轨误差概率曲线比较图 |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.