2021-03-08动力学方程

运动方程
惯性系下，由卫星的摄动理论可知，卫星的运动方程为：

等式右边第一项为中心引力，第二项为摄动力，表示卫星所受的各种摄动力产生的加速度总和。Pd为卫星的动力学模型参数矢量（如光压参数、大气阻力参数等），认为其为常矢量。
卫星所受力可以分为保守力和非保守力。图1给出了卫星上所受各种力产生的加速度大小随轨道高度的变化；图2给出了三类卫星各摄动力的量级（LEO、MEO、GEO）；
可以看出：

N体运动、太阳光压随轨道高度增加而增大；地球中心引力及镶滚摄动力随轨道高度高度增加而减小。
MEO和GEO卫星受到的大气阻力很小，几乎可以忽略。
摄动力
①地球引力
地球引力包括地球中心引力和地球非球型摄动力。地球对外部点的引力位事一个调和函数。EGM2009可以达到2159阶次的位系数，实际定轨时只需要截取一定借此即可。
（重力场模型为地固系，应为至卫星需转换到惯性系）
②N体摄动
太阳、月球及其他行星摄动天体的引力，称为N体摄动。
天体在惯性系下的位置向量可通过查星历表获取，一般选取DE405星历。
见链接：ftp://ssd.jpl.nasa.gov//pub/eph/export/DE405
③大气阻力
产生影响的主要摄动力之一，大气阻力和大气密度的变化相关。
大气阻力系数描述大气与卫星表面的相互作用，受卫星形状、材料等影响，球星卫星可以较准确的确定；非球形难以预测，一般取1.5-3.0。在定轨中一般当作位置参数估计。
大气密度使用数值密度模型计算。（模型很多，如DTM94、Jacchia等）。
④太阳光压
太阳光压摄动力分为两部分：卫星星体部分和太阳帆板部分产生的太阳光压摄动加速度。
计算公式查文件。
⑤低轨卫星经验力模型

开源POD代码：
https://bitbucket.org/geoscienceaustralia/workspace/projects/GACS
计划学习下此开源代码，记录下学习笔记。