近地卫星的运动受到各种力量的影响。其中一种力是地球的中心引力，其他的力称为摄动。这些扰动可分为引力和非引力。在这种情况下，运动方程可以写成:
r̈= - (GM / r ^ 3) * r +γ_p
γ_p是扰动力引起的附加加速度的矢量。
γ_p = r̈_E + r̈_ + 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 =由未建模力引起的加速度
这里我使用了以下积分器和力模型来模拟卫星的摄动:
积分器:带步长控制的变阶Radau IIA积分器
力模型:
-地球重力场(GGM03S模型)
-太阳系行星的重力(行星的位置由JPLDE436计算)
-阻力效应使用Jacchia- bowman 2008, NRLMSISE-00, MSIS-86, Jacchia 70或改进的哈里斯-普里斯特大气密度模型(在加速。你可以取消评论你最喜欢的模型)
-太阳辐射压力使用几何或圆柱形阴影模型
-固体地球潮汐(IERS公约2010)
——海洋潮汐
——广义相对论
-使用IAU 2006决议的ECEF2ECI和ECI2ECEF转换
通过运行HPOP.m开始模拟。在InitialState.txt中设置你最喜欢的卫星的初始值。在高档。M你可以通过将它们设置为1来考虑不同的摄动:
AuxParam。n = 70;中心物体重力场的最大程度%
AuxParam。m = 70;中心物体重力场最小阶数%
AuxParam。太阳= 1;太阳的扰动
AuxParam。月亮= 1;%月球摄动
AuxParam。行星= 1;%行星的扰动
AuxParam。sRad = 1;%太阳辐射压力
AuxParam。阻力= 1;%大气阻力
AuxParam。SolidEarthTides = 1;固体潮%
AuxParam。OceanTides = 1;%海洋潮汐
AuxParam。相对论= 1;%广义相对论
引用:
蒙特布鲁克O.，吉尔E.;卫星轨道:模型、方法和应用;施普林格-海德堡;更正第三版(2005)。
蒙滕布鲁克，弗莱格;个人计算机上的天文学;施普林格-海德堡;第四版(2000年)。
Seeber g;卫星大地测量学;Walter de Gruyter，柏林，纽约;第二完全修订和扩展版(2003年)。
Vallado d;天体动力学基础与应用;麦格劳-希尔,纽约;第3版(2007年)。
尼玛。2000.1984年美国国防部世界大地测量系统。NIMA-TR 8350.2，第三版，修正案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中央文件交换。检索．
