cos E e cos f 1 e cos E
见章仁为“卫星轨道姿 态动力学与控制”，p5 －7
根据上式可由平近点角 M 迭代求出偏近点角 E 、 再求出真近点角 f。 从而确定航天器的运动。
a(1 e ) r 1 e cos f
2
2018年10月8日星期一
因此，利用轨道根数可以很直观地 表示航天器的运动，并且只需求解 代数方程。
p h2

πab 1 A h T 2 2π ab T p
p a(1 e2 ) b 1 e2
T 2π
a3

2π 因此轨道平均角速度 n 为： n f 3 T a
2018年10月8日星期一 Page 27
定义：
平近点角M ：航天器从 近地点开始按平均角速 度 n 转过的角度。
航天器的开普勒三大定律
谐和定律：航天器轨道半长 轴的三次方同轨道周期的平 方成正比。
a k 2 T
3
a
T
是轨道半长轴 是航天器的运行周期
k
是与轨道无关的常数
a
a
Page 4
2018年10月8日星期一
轨道的几何描述
O为地球的质心， 也是椭圆的一个焦点. S为航天器的质心.
S
b A
p
r
O
P
P 是近地点 (perigee） A 是远地点 (apogee） a 是半长轴 (semi-major axis) b 是半短轴 (semi-minor axis) p 是半通径 (semi-parameter) e 是偏心率 (eccentricity) c 是半焦距 (semi-focus)
2018年10月8日星期一 Page 20
1、问题的提出
如果用航天器的初始位置和速度 ( x0 , y0 , z0 , x 0 , y 0 , z 0 ) 来描述航天器的运动，则在任一时刻，需要求解 微分方程才能确定航天器的位置，不方便。 另一方面，我们已知航天器在某一个平面内的运动 轨迹为圆锥曲线，如果已知： （1）轨道平面在空间惯性坐标系中的方位； （2）圆锥曲线的方向（长半轴方向）； （3）在某一时刻航天器在轨道的某一个点上， 则可以通过求解代数方程确定任一时刻航天器的位 置。
d2r r 3 2 dt r
x x r3 0 y y 3 0 r z z r3 0
r x y z
2 2
2
如果给定初始条件：
x0 , y0 , z0 , x0 , y0 , z0
就可以计算出以后任意时刻航天器的位置和速度。
2018年10月8日星期一 Page 7
再利用
p

p3

可得到
t


0
dt (1 e cos ) 2
r
E
S
（过程略）
e
物理意义： 为积分常数，表示矢径 r 与 e 重合的 时刻，称为过近地点时间。
2018年10月8日星期一 Page 19
§1.4 航天器的轨道要素
前面介绍了航天器轨道的特点及积分情况， 导出了一些积分常数（ E , h, e, ），根据轨道 运动方程，只有六个参数是独立的。 原则上，要唯一确定航天器的轨道，六个独 立的参数可以有多种选取方法，比如取航天器的 初始位置和速度：( x0 , y0 , z0 , x 0 , y 0 , z 0 ) ，也可以取 E, h, e, 。 但在航天领域，一般习惯用另外的六个独立 参数来描述轨道的状况。
2018年10月8日星期一 Page 21
解决方案：
已知航天器的运动轨迹为圆锥曲线，而圆锥 曲线的统一方程为：
p r 1 e cos
除了p、e 外再引入四个量，可以用这六个独立变 量来描述航天器的轨道运动。 在航天领域，一般习惯用下面的六个独立参数 来描述轨道的状况: i、Ω 、ω 、p、e、τ 。 这些 量称为轨道要素，或轨道根数。

r3
[rr ] r r 2 r
Page 15
v h

r
3
[rr ] r r r
2
积分后为
e的方向
1

(v h
r
r 1 v h [ 2 r r ] 0 r r
？) e r
？
1 r e h (v h) h (r v ) 0 r

E
（3）如果E＝0，临界情况，满足 v p
抛物线
由于已经知道航天器的轨道是圆锥曲线，根据 第（2）点，E<0时r有界，因此是椭圆轨道。 根据第（1）点，E>0时r可以无界，因此是 双曲线轨道。
2018年10月8日星期一 Page 13
2、动量矩积分
d 2r r 3 2 dt r
0 方程两边叉乘 r： r v
数学概念：微分方程的定解由初始条件确定；而方程的积分常数 是初始条件的某种组合。因此方程独立的积分常数数目不超过初 始条件的数目。在轨道问题中，积分常数不超过6个。
因此可能还存在另1个积分常数…
2018年10月8日星期一 Page 18
4、时间积分
利用 h r
2
及 r
h2
p 1 e cos
M n(t )

a
3
(t )
E O f
S
偏近点角 E ：椭圆轨道 存在内、外接圆，航天 器在内、外接圆上的投 影点与椭圆中心对应的 夹角。如图。
2018年10月8日星期一
Page 28
各角度的关系
M

a
3
(t )
求微分方程 与求代数方 程的比较？
M E e sin E
e的大小
e2 e e 1
2 Eh2
2
所以 e 在轨道平面内，且只有一个独立的量。 物理意义此处还不太明确。
关于e的大小，你有何直觉？ 椭圆轨道： E 0 e [0, 1)
2018年10月8日星期一 Page 16
e
1

(v h
r
r
)
v r
E
S e
e的物理意义 两边叉乘r
算例
为解决这 些问题， 需要对轨 道进行深 入研究
问题： （1）如果参数不适当，航天器可能会撞上地球! （2）如何得到希望的轨道？
2018年10月8日星期一 Page 8
一些尝试
假设引力公式为
G ms m r F r r
其中η 不一定为2；Gη为相应的引力常数。 你估计会出现什么现象？
er 1

(v h) r
er 0 可以看出，在一般情况下，
但如果r与v垂直，则 e r 0 所以，e平行于椭圆长轴方向，再根据其大小，e 指向近地点。
2018年10月8日星期一 Page 17
思考
我们已找到了5个积分常数E, h, e。 问题是：当我们求出常数E，h，并为其中所 使用的技巧而得意时，拉普拉斯利用更复杂的技巧 又找到了一个积分常数e…… 那么我们是否求出了微分方程全部的积分常数？ 难到这些微分方程的积分常数会没完没了吗？
2018年10月8日星期一 Page 1
航天器的开普勒三大定律
椭圆定律：航天器绕地球运 动的轨道为一椭圆，地球位 于椭圆的一个焦点上。
2018年10月8日星期一
Page 2
航天器的开普勒三大定律
面积定律：航天器与地球中 心的连线在相同的时间内扫 过的面积相等。
2018年10月8日星期一
Page 3
2018年10月8日星期一
c
a
p a(1 e2 ) b 1 e2 c ae
Page 5
轨道的微分描述
设 Oxyz 为参考坐标系，O为 地球中心，xyz 指向三颗恒星。 设 me 为地球质量，m为航天器 质量，r为航天器的矢径。
Gme m r d2r ma m 2 F 2 dt r r
Page 29
算例
a(1 e2 ) r 1 e cos f
2018年10月8日星期一 Page 22
2、地球赤道惯性坐标系
为了定义轨道根数，有必要先介绍地球赤道惯性坐标系。
定义地球赤道惯性坐标系OXYZ：O在地球中心，X 轴沿地球赤道面与黄道面的交线，Z轴指向北极星。
Z (指向北极星)
黄道面
太 阳
X
春分点方向：春天的头一天地心 与日心的连线
2018年10月8日星期一
x z
F rySE NhomakorabeaO
d2r r 3 2 dt r
G 6.67 1011 m3 / kg s 2
万有引力常数
Gme 3.99 105 km3 / s 2 地心引力常数
这就是航天器绕地球运动的运动微分方程。
2018年10月8日星期一 Page 6
如果在直角坐标系中进行计算：
2018年10月8日星期一
Z
h i
O
r
ω
N
S
e t=τ
Y
Ω
X
Page 25
4、轨道要素描述的公式及计算方法
p r 1 e cos( )
p h
2
S

b
p
O
r
根据几何关系有
a(1 e ) r 1 e cos f
2

ae a
f

N
其中 f 是真近点角：航天器相对于椭圆长轴的极角。 真近点角 f 的变化就是航天器的轨道角速度。
Z
h i
O
r ω
N
S
e t=τ
Y
Ω