高精度轨道传播器：精确模拟卫星的扰动（特殊扰动方法）-matlab开发

近地卫星的运动受到各种力的影响。 这些力之一是地球的中心引力，其他力称为扰动。 这些扰动分为重力和非重力。 在这种情况下，运动方程可以写为： 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 = 由于未建模的力引起的加速度在这里，我使用了以下积分器和力模型来模拟卫星的扰动运动： 积分