近地卫星的运动受到各种力的影响。其中一个力是地球的中央着引力，其他力量被称为扰动。这些扰动被分为重力和非引力力。在这种情况下，运动方程可以写为：
r̈= - （gm / r ^ 3）* r +Γ_p
γ_P是扰动力诱导的附加加速度的向量。
Γ_p=r̈__e+r̈_s+r̈_m+ r̈_p+r̈_e+r̈_o+r̈_d+r̈_sp+r̈__a+r̈_emp
r̈_e=由于地球内的非球形和不均匀质量分布而导致的加速度（中央体）
R̈_S，r̈_m，r̈_p=由于其他天体（太阳，月亮和行星）引起的加速
R̈_E，r̈_o=由于地球和海洋潮汐引起的加速
r̈_d=由于大气阻力导致的加速度
R̈_SP，R̈_A=由于直接和地球反射的太阳辐射压力而导致的加速度
r̈_emp=由于未确定力引起的加速度
在这里，我使用以下集成器和力模型来模拟卫星的扰动运动：
Integrator：具有梯级控制的可变阶层IIA Integrator
力量模型：
- 地球的重力场（GGM03S模型）
- 太阳系行星的重力（通过JPLDE436计算行星的位置）
- 使用Jacchia-Bowman 2008，NRLMSISE-00，MSIS-86，Jacchia 70或修改后的Harris-Priester大气密度模型（在Accel.m.M中取消注释您最喜欢的型号）
- 使用几何或圆柱暗影模型的太阳辐射压力
- 固体地球潮（IERS约定2010）
- 海潮
- 一般相对性
- 使用IAU 2006分辨率的ECEF2ECI和ECI2ECEF转换
仿真通过运行HPOP.M开始。在InitialState.txt中为您最喜欢的卫星设置初始值。在HPOP.M中，您可以通过将其设置为如下方式来考虑不同的扰动：
auxparam.n = 70;中央体重的最大程度的最大程度
Auxparam.m = 70;中央体重引力领域的最低阶数
auxparam.sun = 1;％扰动太阳
auxparam.moon = 1;％扰动月亮
auxparam.planets = 1;％的行星扰动
auxparam.srad = 1;％太阳辐射压力
auxparam.drag = 1;％大气拖累
auxparam.solidearthtiates = 1;％固体地球潮汐
auxparam.oceantides = 1;％海洋潮汐
auxparam.relativity = 1;％相对论
参考：
Montenbruck O.，Gill E .;卫星轨道：型号，方法和应用;斯普林克·韦尔格拉格，海德堡;纠正第3次印刷（2005）。
Montenbruck O.，Pfleger T。个人计算机上的天文学;斯普林克·韦尔格拉格，海德堡;第四版（2000）。
Seeber G .;卫星测地;沃尔特德格莱特，柏林，纽约;第二次完全修订和扩展版（2003）。
Vallado D.A;actrocumics和申请的基础;麦格劳山，纽约;第3版（2007）。
尼玛。2000年。国防部世界大地测量系统1984。Nima-TR 8350.2，3 ed ed，修正1.华盛顿特区：总部，国家图像和映射机构。
http://sol.spacenvironment.net/jb2008/
http://ssd.jpl.nasa.gov/?ephemerides.
http://celestrak.com/spacedata/
引用
Meysam Mahooti（2021）。高精度轨道传播者（//www.tianjin-qmedu.com/matlabcentral/fileexchange/55167-high-precision-Orbit-propagator），Matlab中央文件兑换。检索到。
