地球卫星三维运行轨道MATLAB仿真

1、问题的描述

3 轨道上运行的地球卫星，根据牛顿第二定律F=ma以及万有引力定律F=-GmME*r/r，

3可得a=-GME*r/r，即

x''= -GME*x/r3 3 y''= -GME*y/r; z''= -GM*z/r3

E (1)

式中，(x，y，z)表示卫星的三维坐标，为G=6.672*10(N m/kg)引力常数，-1122

ME=5.97*1024(kg)是地球的质量。假定卫星的三个方向的初始位置和速度如下

[x(0),y(0),z(0),x'(0),y'(0),z'(0)]=[2043922.17 8186504.63 4343461.71 -5379.54 -407.10 3516.05]。

该卫星轨道求解过程实际上是求解一个二阶常微分方程，可首先将该方程转换为一阶常微分方程，令X=[x,y,z,x',y',z']，故公式(1)可转化为 T

X(4) X(5) X(6) X'(t) = A X(1) A X(2) A X(3) A=-GME/r3 (2)

初始条件即为X(0)=[2043922.17 8186504.63 4343461.71 -5379.54 -407.10 3516.05]。

2、MATLAB仿真代码

分两段程序：(1)子程序将二阶微分方程转换为一阶微分方程，代码如下

function fy=vdp(t,x)

r=x(1)^2+x(2)^2+x(3)^2;

G=3.986005e14;

A=-G/r^(3/2);

fy=[x(4)