matlab 轨道动力学,轨道动力学数学方法作业.pdf

哈尔滨工业大学

轨道动力学中的数学方法作业

学院 航天学院

班级 1318201

姓名 宋俊霖

学号 1133730115

指导教师：张刚

1. 利用轨道六根数和位置速度的转化关系(或Lagrange 系数法)，编写Matlab

程序求解初值问题，并画出2 个周期内的轨迹曲线。初始时刻卫星轨道参数

E [a , e ,i , ,  , ] [10000km, 0.1,10, 20,30, 60] 。(必选)

2 2 2 2 2 2 2

解：

1、求解思路

1) 求出t 时刻对应平近点角M

2) 求出M 对应偏近点角E

3) 求出E 对应真近点角

4) 利用轨道方程求出r () 在轨道坐标系下坐标

5) 通过坐标旋转变换到地心惯性系下

6) 绘图

2、结果

周期T =9952s

A、轨道三维图

2000

1000

0

-1000

-2000

1

1

0

0.5

x 104 0

-1 -0.5 4

x 10

-1

-2 -1.5

B、x 坐标随时间变化

4

x 10

1

0.5

0

x

-0.5

-1

-1.5

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2

t 4

x 10

C、y 坐标随时间变化

4

x 10

1

0.5

0

y

-0.5

-1

-1.5

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2

t