论文部分内容阅读
这篇论文根据作者在硕士生期间参加的一项课题撰写而成,旨在研究半分析方法运用于地球卫星轨道长期预报的性能(速度和精度),并给出半分析法针对不同类型轨道可供参考的精度和速度。第一章绪论部分。首先,介绍了本课题的背景,即为了能够充分掌握空间目标的长期演化规律和分布情况,需要对它们的轨道进行长期预报,考虑到速度和精度的实际需求,半分析方法是一种十分合适的预报手段。虽然半分析法不是一种新方法,但是几乎找不到定量评估其性能的文献,这给一些实际应用造成了困难。在这样的前提下,我们定量分析了半分析法在常用地球卫星轨道中的精度和速度,该结果可以帮助我们了解半分析法在长期预报中的适用性。其次,介绍了半分析方法目前国内外的研究情况。最后,简单介绍了论文中用到的时间系统、坐标系统、摄动方程和单位系统。第二章介绍现有卫星预报方法。众所周知,卫星轨道预报实为求解一个常微初值问题的方程组,目前较为流行的有三种做法,即数值法、分析法和半分析法。数值法,即利用数值积分器一步步外推初始根数,直到到达预报时间点。该方法的优点是实现简单,精度高,缺点则是计算耗时较长。分析法,即显式地给出各个根数在受摄二体问题下的分析解,由于已经给出了摄动解,所以在有了初始根数的情况下,可以方便地直接给出预报时间点处的根数。它的优点是速度快,因为不用像数值法一样一步一步外推,缺点则是受到收敛弧段的限制预报时间不宜过长,一般情况下为103弧段(归一化单位下)。半分析法利用数值方法积分平均化后的动力系统。求解则采用数值法外推,因此可以想象其精度优于分析法而速度将优于数值法。半分析法的核心内容就是如何得到平均化右函数,因此着重介绍了分析平均和数值平均两种平均化方法,并简单比较了两者的性能,阐述了最终选择分析平均的理由。第三章各摄动力在半分析法中的处理。本章依次介绍非球形、第三体、光压、大气在半分析法中的具体处理细节。每新增一个摄动源的同时,将半分析法应用于不同高度卫星,以便直观地了解新的摄动源带来的动力学特征以及对于精度和速度的影响。第四章将半分析方法运用于不同类型轨道。首先,通过简单的算例给出半分析法不同轨道高度适合的非球形引力阶次。而后,将完整力模型的半分析方法运用于不同类型地球卫星轨道,研究在不同初始根数a,i的组合下,半分析法在100年尺度长期预报中的精度。地球卫星轨道按高度主要划分为低轨(有大气),低轨(无大气),12小时中轨及24小时高轨地球同步。通过比较大量的算例后发现,相比于数值法,半分析法100年尺度下瞬时根数a,e,i的精度误差在O(10-3)量级,Ω的精度误差在10°之内。虽然右函数精确到二阶的半分析法在大部分情况下都表现出了良好的性能,但是仍然会遇到一些特殊算例,半分析法的表现不尽如人意。这些算例有些伴随着特殊的动力学现象,有些则归咎于处理方法上的原因。前者,我们在本章中就给出了解释;后者,我们将在第五章中提出合理的解决方案。第五章特殊问题讨论。首先解释了第四章中倾角i=90°时半分析精度下降问题,而后给出了大偏心率轨道的处理方法,接着讨论了坐标系附加摄动中的补充项在半分析法中是否有必要添加,最后讨论了低轨根数误差图中的“岛状”结构成因。文章针对不同类型的地球卫星轨道,定量地展示了半分析法100年尺度长期预报的性能。当然这里只是一些基础性的精度和速度的讨论,半分析法的应用范围远不止如此,比如卫星碰撞后碎片分布的快速估计,空间目标编目中的定轨、关联等。希望本文的工作可以给相关领域的工作者起到一定借鉴意义。