2008航天器动力学10-摄动理论解答.ppt

上述方程组原则上可以求解。但由于严重非线性，只可能进行数值计算，没有解析解。 另一方面，由于摄动的影响与万有引力相比及其微小，因此轨道要素的变化是缓慢的，我们关心的是轨道的长期摄动。 的变化缓慢。 的变化快。 当方程中的变量变化快慢明显不同时，理论分析中可考虑采用平均法对方程进行简化。 4、长期摄动影响 在数值计算中，如果变量变化缓慢，计算的步长可以大些；如果变量变化剧烈，计算的步长就要小些。如果两种情况同时存在（特征值相差大） ，方程称为“刚性方程” ，也称为“病态方程” ，需要特殊的算法。 平均法的思路是： 用 统一表示各轨道要素 在卫星运行的一个周期内，近似认为轨道要素没有变化，将方程在一个周期内积分，得到每运行一周轨道要素的微小增量 由于?i不是时间的显函数，而是真点角的显函数，因此… 虽然轨道根数的变化很复杂，但是其对 f 的积分有可能有解析解。经过积分，得 clear all syms mu e p theta f omega J2 Re inc disp [－－－－－符号运算－－积分－－－－－－ ] ; theta f+omega; r p/ 1+e*cos f ; dt_df r^2/sqrt mu*p ; %转换关系 fr -3*mu*J2*Re^2/2/r^4* 1-3* sin inc ^2* sin theta ^2 ; ft -3*mu*