论文推荐| 蒋春华：拉格朗日/高斯无奇点卫星运动方程推导与分析|163

《测绘学报》

构建与学术的桥梁拉近与权威的距离

拉格朗日/高斯无奇点卫星运动方程推导与分析

1. 山东大学空间科学研究院, 山东威海 264209;
2. 西安测绘研究所地理信息工程国家重点实验室, 陕西西安 710054;
3. 东南大学仪器科学与工程学院, 江苏南京 210096;
4. 香港理工大学, 香港 999077

收稿日期：2017-02-23；修回日期：2017-07-19

基金项目：国家自然科学基金（41574025）；地理信息工程国家重点实验室开放基金（SKLGIE2016-M-2-1）；国家重点研发计划（2016YFB0501900）

第一作者简介：蒋春华(1988—)，女，博士生，研究方向为轨道理论与GNSS高精度数据处理的理论与方法。E-mail: jiangchunhua@sdu.edu.cn

通信作者：许国昌, E-mail: gcxu@sdu.edu.cn

摘要：针对卫星轨道理论中的奇点问题，对拉格朗日/高斯无奇点卫星运动方程做了深入的分析和探索，从原始拉格朗日和高斯运动方程及其物理意义出发，考虑圆轨道、圆赤道轨道和赤道轨道3种奇点情况，推导了一种新的拉格朗日/高斯无奇点卫星运动方程，并探讨了方程的连续性。该方程消除了零因子，解决了卫星运动方程的奇点问题。

关键词：拉格朗日/高斯无奇点运动方程奇点问题卫星轨道连续性

Derivation and Analysis of Singular-free Lagrangian/Gaussian Equations of Planetary Motion

Abstract: Aiming at the singularity problem in satellite orbit theory, the singularity-free Lagrangian/Gaussian equations of motion is analyzed.Considering the original and physical meaning of the Lagrangian and Gaussian equations of motion, a new Lagrangian/Gaussian singularity-free disturbed equations of motion is proposed and then discussed in three cases:the circular orbit, equatorial orbit, circular and equatorial orbit.Besides, the continuity of these equations is explored.The proposed equations eliminate the zero factor and in this way the singularity problem in the orbital mechanics is solved.

Key words: singularity-free Lagrangian/Gaussian equations of motion singularity problem satellite orbit continuity

当卫星按圆轨道、赤道轨道以及圆赤道轨道运动时，某些卫星开普勒(Kepler)轨道根数发生退化，用于表达卫星轨道的拉格朗日运动方程和源于拉格朗日运动方程的高斯运动方程发生数学奇异，即为奇点问题，它一直是天体力学和N体问题研究中的核心问题之一[1-4]。卫星轨道理论研究主要基于二体问题的摄动理论[5-7]，二阶摄动理论解研究近年来取得了长足进步[8-12]，为奇点问题的进一步研究奠定了基础。众多学者在理论研究方面做了很多努力[13-16]，也取得了不少成果。第1类和第2类奇点问题为解决卫星导航领域的奇点问题提供了可能，但其表达复杂且解决思路和最终表达很难统一[17-18]。文献[19-20]提出参量代换和坐标系旋转等方法，分别从圆轨道，赤道轨道和极轨道的情况，对其方法进行论证和阐述，为奇点问题的解决提供了新的思路。但其忽略了代换参量实际的物理和几何意义，以及转换过程对整个轨道系统求解的可逆性。2012年，文献[21]以不定积分解的形式提出无奇点理论，并于2013年有了以微分方程描述的公式版本[22]。但其公式推导主要是逻辑推理，2014年给出了拉格朗日无奇点运动方程纯数学的严格推导与证明[23]。2015年又给出了高斯无奇点运动方程纯数学的严格推导与证明[24]。上述论文主要是从纯数学角度进行推导证明，参数的物理意义和方程的连续性有待进一步分析。

为了解决上述问题，本文充分考虑参数的实际几何物理意义、总体方程的统一性以及与原始方程的一致性，从拉格朗日无奇点运动方程出发，对奇点问题作了进一步的分析和探索。奇点情况时，某些Kepler轨道根数发生退化，本文基于原始拉格朗日和高斯运动方程及其物理意义，对卫星无奇点方程的判据进行推导与分类，并分别对拉格朗日/高斯无奇点卫星运动方程在圆轨道、赤道轨道和圆赤道轨道情况下存在奇异的轨道参数(升交点赤经、近地点角距和平近点角)赋予新的实际物理和几何意义。在拉格朗日-许与高斯-许卫星运动方程的基础上，根据参数定义对奇异项进行加和变换，从而消除奇异因子，最后给出一种新的拉格朗日/高斯无奇点卫星运动方程的统一表达式。

1 卫星无奇点方程的判据

卫星运动方程的奇点问题源于某些特定情况时Kepler根数的物理意义的模糊，导致拉格朗日运动方程出现小分母或分母为零的问题。拉格朗日卫星运动方程可表示为[7, 14]

式中，函数R称为摄动位函数，包括了除中心引力位以外的全部项，6个开普勒根数a、e、ω、i、Ω、M，即半长轴、偏心率、近地点角距、倾角、升交点赤经、平近点角。文献[21]将其改写成

式中，δ及其下标表示式(1)的右边各项。当e或者sin i趋近于0时，即可能发生奇点问题。奇点问题解决的前提是需要对奇点域做具体的判断。首先，定义径向和赤道方向的误差如下[22]

此处径向误差ae是轨道几何中心O′与地心O的距离，如图 1所示，也是e引起的最大径向误差。赤道偏差a sin i≈ai表示卫星由i引起的垂直于赤道的最大距离(倾角i是卫星的最大纬度)。在积分时刻t0=0，e或者/和i接近于0，则在积分时刻t，偏差e和i就为Δe和Δi。式(3)可写成

图 1 卫星的扁心率和近点角Fig. 1 The eccentricity and anomalies of a satellite

图选项

圆轨道、赤道轨道以及圆赤道轨道的轨道偏差定义为

式中，ε0表示在卫星运动方向上由ΔΩ、Δω、ΔM引起的最大偏差。径向偏差、赤道方向偏差和轨道方向偏差是轨道最大误差的3个组成部分。3个偏差表示(e, ω, i, Ω, M)扰动对轨道产生的最大误差域。实际应用中，假定轨道3个方向偏差εe、εi和ε0分别为1 m，即这里定义的εe、εi和ε0可以作为判据标准。对于不同轨道，判据条件如表 1所示，用于对奇点问题的判断与进一步分析。

表 1 奇点情况的判断条件Tab. 1 Criteria of singularity

奇点类型	判断条件
非奇点	ae≥εe=1, asin i≥εi=1
圆轨道奇点	ae < εe=1, asin i≥εi=1
赤道轨道奇点	ae≥εe=1, asin i < εi=1
圆赤道轨道奇点	ae < εe=1, asin i < εi=1

表选项

2 拉格朗日无奇点卫星运动方程推导分析

拉格朗日无奇点运动方程的数学推导已在文献[21]中详细给出。对于通常非奇点情况拉格朗日方程[15, 22]可写为

对于奇点情况还需进一步分析。由轨道参数的原始定义有：Ω表示在地球赤道平面上，升交点与春分点之间的地心夹角称为升交点赤经；ω表示在轨道平面上近地点与升交点之间的地心角距称为近地点角距；f表示卫星到近地点之间的地心角距，如图 2所示。为计算真近点角，二体问题中引入偏近点角和平近点角的定义，3种近点角(真近点角f、平近点角M、偏近点角E)均是时间的函数，用来表述卫星在ECI坐标系中的位置随着时间的变化。其中，平近点角M采用平均角速度n来描述卫星绕地球的轨道运动，广泛用于解析轨道的理论推导和应用，本文的推导也是基于此。

图 2 轨道几何图Fig. 2 Orbital geometry

图选项

奇点情况分为圆轨道、赤道轨道和圆赤道轨道3种情况。当为圆轨道时，e=0近地点失去意义，无法确定。此时开普勒轨道根数发生退化，平近点角M和近地点角距ω意义发生改变。因此，从卫星实际运动情况分析，可定义新的平近点M′表示卫星与升交点的夹角，即

当为赤道轨道时，sin i=0，升交点失去意义，无法确定。此时开普勒轨道根数发生退化，升交点赤经Ω和近地点角距ω意义发生改变。因此，可定义新的近地点角距ω′表示近地点与春分点的夹角，即

当为圆赤道轨道时，e=0且sin i=0，近地点和升交点均失去意义，无法确定。此时开普勒轨道根数发生退化，平近点角M和近地点角距ω以及升交点赤经Ω意义均发生变化。因此，定义一个新的平近点M″表示卫星与春分点的夹角，即

2.1 圆轨道情况

圆轨道情况ae < εe=1、a sin i≥εi=1时，从物理意义分析，近地点无法确定，相关轨道参数对时间的偏导数也无法确定。从数学角度分析e除到式(2)右边，分母存在零因子，即发生奇点。由式(7)，可得

因此，圆轨道奇点情况的运动方程由式(6)中一、二、三、四式，以及下式组成

2.2 赤道轨道情况

赤道轨道情况ae≥εe=1、asin i < εi=1时，从物理意义分析，升交点无法确定，相关轨道参数对时间的偏导数也无法确定。从数学角度分析，sin i除到式(2)右边，分母存在零因子，即发生奇点。由式(8)可得

因此，赤道轨道奇点情况的运动方程由式(6)中一、二、三、六式，以及下式组成

2.3 圆赤道轨道情况

圆赤道轨道情况ae < εe=1、asin i < εi=1时，从物理意义分析，升交点和近地点无法确定，相关轨道参数对时间的偏导数也无法确定。从数学角度分析，sin i、e除到式(2)右边，分母存在零因子，即发生奇点。由式(9)，可得

因此，圆赤道轨道奇点情况的运动方程由式(6)中一、二、三式，以及下式组成

3 高斯无奇点卫星运动方程推导分析

非奇点情况下高斯运动方程为[24]

式中，fa、fh、fr分别表示轨道坐标系的3个方向的摄动力。对于奇点情况下的高斯运动方程做进一步推导。

3.1 圆轨道情况

ae < εe=1、asin i≥εi=1，根据拉格朗日圆轨道运动方程，得

因此，圆轨道奇点情况的运动方程由式(16)中一、二、三、五式，以及下式组成

3.2 赤道轨道情况

ae≥εe=1、asin i < εi=1，根据拉格朗日赤道轨道运动方程，得

因此，赤道轨道奇点情况的运动方程由式(16)中一、二、三、六式，以及下式组成

3.3 圆赤道轨道情况

ae < εe=1, asin i < εi=1，根据拉格朗日圆赤道轨道运动方程，得

因此，圆赤道轨道奇点情况的运动方程由式(16)中一、二、三式，以及式(22)组成

4 无奇点卫星运动方程的连续性分析

总结上述非奇点与奇点的各种情况，最终统一的拉格朗日无奇点卫星运动方程为

最终统一的高斯无奇点卫星运动方程为

卫星的运动和力模型是连续的，进一步分析卫星运动方程的连续性可知，式(23)、(24)对于任何一种奇点情况(圆轨道奇点、赤道轨道奇点、圆赤道轨道奇点)方程是非间断的，奇点情况下运动方程中重新定义的组合量及非奇点情况下运动方程的轨道参数在其对应的判据域内均是连续的。因此，新的拉格朗日和高斯无奇点卫星运动方程具有连续性。

5 仿真与分析

(1) 圆轨道情况。轨道坐标系下状态矢量进行一次旋转，分析角度物理意义为

(2) 赤道轨道情况。轨道坐标系下旋转矩阵，分析角度物理意义可知

式中，R1、R3为x和z中的旋转矩阵[16]，圆赤道轨道两点兼有。

下面利用本文推导的拉格朗日无奇点卫星运动方程描述卫星运动，轨道初值如表 2所示(其中后3个可表示新定义的根数)。首先基于式(23)求出各时刻的卫星的轨道根数，然后利用式(25)、(26)将其转换成卫星状态矢量，最后利用自主编写的卫星轨道三维可视化仿真软件分别绘制二体运动圆轨道、赤道轨道以及圆赤道轨道轨迹[25]，如图 3所示。将计算轨道与积分轨道对比，轨道差在厘米级。高斯方程同理。

图 3 圆轨道、赤道轨道以及圆赤道轨道卫星轨迹Fig. 3 The circular/equatorial/circular and equatorial satellite orbit

图选项

表 2 圆轨道、赤道轨道及圆赤道轨道的轨道初值Tab. 2 Initial orbit of circular/equatorial/circular and equatorial satellite orbit

轨道类型	a/km	e/(°)	i/(°)	Ω/(°)	ω/(°)	M/(°)
圆轨道	8000	0	30	0	—	80
赤道轨道	12 000	0.3	0	—	60	20
圆赤道轨道	11 000	0	0	—	—	10

表选项

从上述过程与结果可知，本文推导的无奇点卫星运动方程，可描述二体运动中的圆轨道、赤道轨道以及圆赤道轨道卫星运动，且计算精度与非奇点二体运动相当，可用于解决二体运动中的奇点问题，从而在一定程度上证明本文推导公式的正确性和可用性。

6 结论

本文推导了卫星拉格朗日/高斯无奇点运动方程，无任何近似，与原始拉格朗日/高斯运动方程具有逻辑一致性；奇点情况下，对参数含义重新定义，所得方程在奇点情况下仍具有几何和物理意义；圆轨道、赤道轨道和圆赤道轨道下的卫星运动方程的表达更加简洁，便于理解和计算；轨道方程右端完全消除了零因子，从根本上解决卫星运动方程的奇点问题，并且该卫星运动方程具有连续性。

【引文格式】蒋春华, 徐天河, 乔晶, 等. 拉格朗日/高斯无奇点卫星运动方程推导与分析[J]. 测绘学报，2018，47(4)：455-464. DOI: 10.11947/j.AGCS.2018.20170082

“子午仪”卫星导航系统——世界上第一个卫星导航系统

“2017年度中国遥感领域十大事件”揭晓

同行评议的“丑陋面”！

涨知识！全球有多少CORS站？

2018年泰晤士世界大学200强榜单，中国入榜12所，快来看看都有谁？

杨元喜院士: 北斗三号试验系统性能分析【附PPT】

数学能力是中兴与华为的唯一区别

关于在学术论文署名中常见问题或错误的诚信提醒

新时代之光 | 武大六院士当选“荆楚楷模”

足球大亨许家印1000亿杀入信息产业，重点布局商业航天、智能交通、AI

宁津生院士：测绘界的十几位院士正在研究这个问题