综合评价线性多步积分公式轨道积分性能的两项新指标

概述了制约线性多步积分公式轨道积分状态的多种因素。提出了综合评价线性多步积分公式积分性能的两项新指标。建议在对数值计算有较高精度要求的科研项目中，应将构造并选择适合研究项目的线性多步积分公式以及高效的积分方式列为课题前期工作的重要部分。     

综合评价线性多步积分公式轨道积分性能的两项新指标

概述了制约线性多步积分公式轨道积分状态的多种因素 ;提出了综合评价线性多步积分公式积分性能的两项新指标。建议在对数值计算有较高精度要求的科研项目中 ,应将构造并选择适合研究项目的线性多步积分公式以及高效的积分方式列为课题前期工作的重要部分     

推荐数值求解y‘’=f（x,y）的几组织分系数

为求解特殊二阶常微分方程ｙ＇＇＝ｆ（ｘ，ｙ）的初值问题，本文采用最大阶算子方法构造了一类线性多步积分公式，并与Ｃｏｗｅｌｌ方法的同阶公式作了大量的平等计算，通过对不同轨道类型、不同步长、不同积分间隔时的计算结果的全面仔细地分析比较，我们从八阶、十阶、十二阶、十四阶中各自选定一组积分系数，推荐给同行计算使用，结果表明，采用本文推荐的积分方法计算天体轨道是有益的，因为它的积分精度以及积分过程中误差累积     

推荐一组对长、短积分间隔皆优的线性多步积分公式

推荐了一组对短、中、长积分间隔以及带耗散力或较大偏心率等多种类型卫星轨道数值积分皆可获取较高计算精度的线性多步积分公式（MTM）,供读者选择使用。     

可用于带耗散力卫星轨道长间隔积分的数值积分方法

本文分析了对称方法不适合带耗散力的卫星轨道长间隔积分的缺陷和本质原因,并针对这个问题,采用文所介绍的方法,构造并推荐了一组积分公式。     

推荐数值求解y~″＝f（x，y）的几组积分系数

为求解特殊二阶常微分方程ｙ″＝（ｘ，ｙ）的初值问题，本文采用最大阶算子方法构造了一类线性多步积分公式，并与Ｃｏｗｅｌｌ方法的同阶公式作了大量的平行计算，通过对不同轨道类型、不同步长、不同积分间隔时的计算结果的全面仔细地分析比较，我们从八阶、十阶、十二阶、十四阶中各自选定一组积分系数，推荐给同行计算使用，结果表明，采用本文推荐的积分方法计算天体轨道是有益的，因为它的积分精度以及积分过程中误差累积的方式都十分明显地好于同阶的Ｃｏｗｅｌｌ方法。     

利用IGS星历预报GPS卫星轨道

在动力学轨道拟合以及轨道积分的基础上,提出了基于IGS精密星历的GPS卫星轨道预报方法。该方法首先利用已知的IGS精密星历作为虚拟观测值,采用动力学方法拟合出GPS卫星的初始轨道和动力学参数,然后再通过积分来预报GPS卫星的轨道。主要讨论了基于不同弧段的IGS星历时,该方法对GPS卫星轨道的拟合和预报情况。研究结果显示:对于6 d弧段以内的IGS精密星历,其拟合轨道与IGS精密星历差值的三维RMS值均优于4 cm,随着拟合弧段的增加,拟合残差变大;当利用2～6 d弧段的IGS星历来预报GPS轨道时,大部分卫星第1天、第7天和第30天的三维预报精度可优于0.1 m、3 m和100 m。其中,2d弧段的IGS星历对GPS卫星第1天和第7天的预报结果最好,5 d弧段的IGS星历对GPS卫星第30天的预报结果最好。     

卫星跟踪技术的发展

本文介绍了卫星跟踪技术发展的历史、现状、卫星轨道分类方法、几种观测方法和它们的精度，卫星轨道测定和运动微分方程积分方法、几个著名的应用卫星系统、卫星技术对世界科技发展产生的巨大贡献、以及中国卫星跟踪技术发展的历史和现状。     

LAGEOS卫星轨道改进的一种简易方法

分析表明，卫星轨道偏差对预报的影响主要表现在M和Ω这两个轨道根数上。本给出了M、Ω的修正公式，并结合LAGEOS测距数据，对卫星的短期预报进行修正，取得了满意的结果。     

月球卫星轨道力学综述

月球探测器的运动通常可分为3个阶段，这3个阶段分别对应3种不同类型的轨道：近地停泊轨道、向月飞行的过渡轨道与环月飞行的月球卫星轨道。近地停泊轨道实为一种地球卫星轨道；过渡轨道则涉及不同的过渡方式(大推力或小推力等)；环月飞行的月球卫星轨道则与地球卫星轨道有很多不同之处，它决不是地球卫星轨道的简单克隆。针对这一点，全面阐述月球卫星的轨道力学问题，特别是环月飞行中的一些热点问题，如轨道摄动解的构造、近月点高度的下降及其涉及的卫星轨道寿命、各种特殊卫星(如太阳同步卫星和冻结轨道卫星等)的轨道特征、月球卫星定轨等。     

Two multiderivate multistep (MDMS) integrators

A preliminary survey of multiderivative multistep integrators is carried out. It is found that all of them are much more accurate than the classical linear multistep methods, but most of them have poor stability. After parameter adjustment, two of them (called MDMS I and MDMS II by us) are competitive with or superior to the classical methods in some aspects, such as accuracy and stability. MDMS I behaves especially well in all the cases which have been studied.     

Hamilton系统数值计算的新方法

系统地介绍了近年来对Ｈａｍｉｌｔｏｎ系统数值计算新建立的辛算法和线性对称多步法，并对它们在动力天文中的应用作了一简要回顾。     

Orbitally stable multistep methods

Lambert and Watson (1976) examine the family of symmetric linear multistep methods for the special second-order initial value problem, and connect the property of symmetry with a property of periodicity. The problems of celestial mechanics may be formulated as second-order initial value problems, but these frequently incorporate the first derivative explicity. It is common for such equations to be reduced to a system of first-order equations. Thus motivated, we utilize ideas from the aforementioned paper to determine the family of linear multistep methods for first-order initial value problems that possess an analogous property of periodicity. This family of orbitally stable methods is illustrated by examining the regularized equations of motion of an artificial earth satellite in an oblate atmosphere.     

Adams—Cowell方法与KSG积分器的比较

在人造地球卫星精密定轨中,有摄星历等量的计算常采用Adams-Cowell方法,美国Texas大学空间研究中心(CSR)的定轨软件中则采用了一种有别于Adams-Cowell方法的KSG积分器。本文对这两种线性多步法作了全面比较,并用典型算例作了数值验证,列出了两种方法中卫星轨道沿迹误差的状况,以此表明为什么人们常采用Adams-Cowell方法。     

A new family of multistep numerical integration methods based on Hermite interpolation

In this paper, a new family of explicit and implicit multistep methods is presented both for the error-controlled and uncontrolled modes. The main concept is to replace the Newton interpolation with the Hermite interpolation, where the Hermite polynomial is fitted to the function values and its derivatives. This idea is very useful in the numerical solution of problems (e.g., orbit propagation problem) where higher-order derivatives can easily be computed. In addition to the theoretical concept, the stability regions of the proposed methods are determined. The new methods are more stable than the well-known multistep numerical integrators (i.e., Adams–Bashforth and Adams–Bashforth–Moulton) in the explicit, implicit, and predictor–corrector forms. Using the second-order derivatives gives smaller error constants in the proposed method. The new integrators are numerically tested for a few examples, and the solutions are compared with those of the well-known multistep methods. Moreover, the CPU time and absolute integration error are compared in the satellite orbit propagation problem using various integration methods. The CHAMP mission, i.e., a German small-satellite mission for geoscientific and atmospheric research and applications, is considered as a case study for comparing the achievable accuracy of the proposed method with the existing method for solving the two-body problem.     

动力天文中常用多步方法的比较（Ⅱ）

本文从几个方面给出了动力天文中三种常用的多步积分方法的比较结果,为积分器的选择提供了依据。     

Symplectic methods and their application to the motion of small bodies in the Solar System

Symplectic methods have been widely used in Solar System dynamics. This paper discusses both single step and multistep symplectic methods. For single step methods we point out that the modified algorithm (Wisdom et al., 1991, Kinoshita et. al., 1991) can be executed in the mass center coordinate system and in the Jacobian coordinate system. For multistep methods we describe the connections between symmetric and symplectic methods.     

Symplectic methods and their application to the motion of small bodies in the Solar System

Symplectic methods have been widely used in Solar System dynamics. This paper discusses both single step and multistep symplectic methods. For single step methods we point out that the modified algorithm (Wisdom et al., 1991, Kinoshita et. al., 1991) can be executed in the mass center coordinate system and in the Jacobian coordinate system. For multistep methods we describe the connections between symmetric and symplectic methods.     

Multistep methods of numerical integration using back-corrections

A new class of linear multistep methods is proposed for the solution of the equations of motion of certain dynamical systems encountered in celestial mechanics and astrodynamics. These methods are distinguished from the classical predictor-corrector methods in that they permit back-corrections of the solution to be made. As the integration advances in time, the numerical solution is corrected or improved at certain points in the past. The enhanced numerical stability of these methods allows the meaningful application of high-order algorithms. Consequently, stepsizes larger than those attainable with the classical methods may be adopted and thus greater over-all efficiency may be realized. The application of these methods to the problem of determining the orbit of an artificial satellite is accomplished and the results are compared with those obtained using classical methods.