轨道建模是创建数学模型的过程，以模拟一个大质量物体在重力作用下绕另一个大质量物体绕轨道运动时的运动。通常将其他力（例如来自第三体的引力，空气阻力，太阳压力或来自推进系统的推力）建模为次要作用。由于需要将小扰动建模到非常大的轨道，因此直接对轨道进行建模可能会限制机器精度的极限。因此，通常使用摄动方法来对轨道进行建模，以获得更高的精度。轨道模型通常使用特殊的扰动方法在时间和空间中传播。这是通过首先将轨道建模为开普勒轨道来执行的。然后将扰动添加到模型中，以说明影响轨道的各种扰动。特殊的扰动可以应用于天体力学中的任何问题，因为它不仅限于扰动力小的情况。特殊的摄动方法是最精确的机器生成的行星星历的基础，例如喷气推进实验室发展星历。
在这里，我使用以下积分器和力模型来模拟卫星的扰动：
积分器：具有步长控制的可变阶Radau IIA积分器。
力模型：地球重力场（GGM03S模型），太阳系行星的重力（行星的位置由JPLDE430计算），使用NRLMSISE00大气密度模型的阻力效应，使用圆柱模型的太阳辐射压力，固体地球潮汐（IERS） 《 2010年公约》，海潮和广义相对论。
仿真通过运行test_Envisat.m开始。在InitialState.txt中，为您喜欢的卫星设置初始值。在test_Envisat.m中，可以通过将它们设置为1来考虑不同的扰动：
AuxParam.n = 40; % maximum degree of central body's gravitation field
AuxParam.m = 40; % maximum order of central body's gravitation field
AuxParam.sun = 1; % perturbation of the Sun
AuxParam.moon = 1; % perturbation of Moon
AuxParam.planets = 1; % perturbations of planets
AuxParam.sRad = 1; % solar radiation pressure
AuxParam.drag = 1; % atmospheric drag
AuxParam.SolidEarthTides = 1; % solid Earth tides
AuxParam.OceanTides = 1; % ocean tides
AuxParam.Relativity = 1; % general relativity
参考文献：
Montenbruck O.，Gill E .; 卫星轨道：模型，方法和应用；海德堡Springer Verlag；修正了《第三次印刷》（2005年）。
Montenbruck O.，Pfleger T .；个人计算机上的天文学；海德堡Springer Verlag；第四版（2000）。
Seeber G .；卫星大地测量；沃尔特·德·格鲁特（Walter de Gruyter），柏林，纽约；第2次完全修订和扩展版（2003年）。
巴利亚多D.A; 天体动力学与应用基础；纽约，麦格劳-希尔；第三版（2007年）。
尼玛 2000年。国防部，世界大地测量系统，1984年。NIMA-TR 8350.2，第三版，修正案1.华盛顿特区：国家影像和地图局总部。