卫星轨道预报的一种分析方法

人造地球卫星的轨道预报是空间环境监测和实时跟踪测量中一个重要环节，由于监测对象众多，要求精度也不太高，通常采用分析法预报．在已有分析法得到t时刻平均根数的基础上给出一种轨道预报方法，由t时刻的平均根数给出该时刻卫星的位置和速度，在此基础上将地球非球形引力摄动的周期项直接用卫星直角坐标的位置和速度分量表示，这样可以避免在计算轨道根数变化的周期项时出现的奇点问题，从而对根数的选择无特殊要求，可适用于各种轨道，简化预报程序和相应的软件，提高预报效率。     

Behaviour of the SMF method for the numerical integration of satellite orbits

The SMF algorithms were recently developed by the authors as a multistep generalization of the ScheifeleG-functions one-step method. Like the last, the proposed codes integrate harmonic oscillations without truncation error and the perturbing parameter appears as a factor of that error when integrating perturbed oscillations. Therefore they seemed to be convenient for the accurate integration of orbital problems after the application of linearizing transformations, such as KS or BF. In this paper we present several numerical experiments concerning the propagation of Earth satellite orbits, that illustrate the performance of the the SMF method. In general, it provides greater accuracy than the usual standard algorithms for similar computational cost.     

精密定轨用地球大气模型误差的补偿方法

本文提出了一种用于精密轨道确定的地球大气模型误差的补偿方法，这种方法采用改进 大气周日峰角的技术，有效地补偿了大气密度模型的误差，并具有较好的物理意义．与传统 的每圈改进一次加速度的方法比较，轨道确定达到了相当的径向精度，且避免了轨道改进参 数间的相关性，同时解算参数的数量也少了一半．     

Dynamics of Two Planets in the 2/1 Mean-Motion Resonance

The dynamics of two planets near a first-order mean-motion resonance is modeled in the domain of the general three-body planar problem. The system studied is the pair Uranus-Neptune (2/1 resonance). The phase space of the resonance and near-resonance regions is studied by means of surfaces of section and spectral analysis techniques. After a thorough investigation of the topology of the phase space, we find that several regimes of motion are possible for the Uranus-Neptune system, and the regions of transition between the regimes of motion are the seats of chaotic motion. This revised version was published online in July 2006 with corrections to the Cover Date.     

A survey of orbits of co-orbitals of Mars

Many asteroids with a semimajor axis close to that of Mars have been discovered in the last several years. Potentially some of these could be in 1:1 resonance with Mars, much as are the classic Trojan asteroids with Jupiter, and its lesser-known horseshoe companions with Earth. In the 1990s, two Trojan companions of Mars, 5261 Eureka and 1998 VF₃₁, were discovered, librating about the L₅ Lagrange point, 60° behind Mars in its orbit. Although several other potential Mars Trojans have been identified, our orbital calculations show only one other known asteroid, 1999 UJ₇, to be a Trojan, associated with the L₄ Lagrange point, 60° ahead of Mars in its orbit. We further find that asteroid 36017 (1999 ND₄₃) is a horseshoe librator, alternating with periods of Trojan motion. This asteroid makes repeated close approaches to Earth and has a chaotic orbit whose behavior can be confidently predicted for less than 3000 years. We identify two objects, 2001 HW₁₅ and 2000 TG₂, within the resonant region capable of undergoing what we designate “circulation transition”, in which objects can pass between circulation outside the orbit of Mars and circulation inside it, or vice versa. The eccentricity of the orbit of Mars appears to play an important role in circulation transition and in horseshoe motion. Based on the orbits and on spectroscopic data, the Trojan asteroids of Mars may be primordial bodies, while some co-orbital bodies may be in a temporary state of motion.     

Effective Hamiltonian for the D'Alembert Planetary Model Near a Spin/Orbit Resonance

The D'Alembert model for the spin/orbit problem in celestial mechanics is considered. Using a Hamiltonian formalism, it is shown that in a small neighborhood of a p:q spin/orbit resonance with (p,q) different from (1,1) and (2,1) the 'effective' D'Alembert Hamiltonian is a completely integrable system with phase space foliated by maximal invariant curves; instead, in a small neighborhood of a p:q spin/orbit resonance with (p,q) equal to (1,1) or (2,1) the 'effective' D'Alembert Hamiltonian has a phase portrait similar to that of the standard pendulum (elliptic and hyperbolic equilibria, separatrices, invariant curves of different homotopy). A fast averaging with respect to the 'mean anomaly' is also performed (by means of Nekhoroshev techniques) showing that, up to exponentially small terms, the resonant D'Alembert Hamiltonian is described by a two-degrees-of-freedom, properly degenerate Hamiltonian having the lowest order terms corresponding to the 'effective' Hamiltonian mentioned above.     

Kepler Equation solver

Kepler's Equation is solved over the entire range of elliptic motion by a fifth-order refinement of the solution of a cubic equation. This method is not iterative, and requires only four transcendental function evaluations: a square root, a cube root, and two trigonometric functions. The maximum relative error of the algorithm is less than one part in 10¹⁸, exceeding the capability of double-precision computer arithmetic. Roundoff errors in double-precision implementation of the algorithm are addressed, and procedures to avoid them are developed.     

Numerical Exploration of Chermnykh’s Problem

We study the equilibrium points and the zero-velocity curves of Chermnykh’s problem when the angular velocity ω varies continuously and the value of the mass parameter is fixed. The planar symmetric simple-periodic orbits are determined numerically and they are presented for three values of the parameter ω. The stability of the periodic orbits of all the families is computed. Particularly, we explore the network of the families when the angular velocity has the critical value ω = 2√2 at which the triangular equilibria disappear by coalescing with the collinear equilibrium point L₁. The analytic determination of the initial conditions of the family which emanate from the Lagrangian libration point L₁ in this case, is given. Non-periodic orbits, as points on a surface of section, providing an outlook of the stability regions, chaotic and escape motions as well as multiple-periodic orbits, are also computed. Non-linear stability zones of the triangular Lagrangian points are computed numerically for the Earth–Moon and Sun–Jupiter mass distribution when the angular velocity varies.     

Synthetically adaptive robust filtering for satellite orbit determination

The precise orbit information is the premise of the application of satellite, and the precise orbit deter- mination includes dynamic and geometric methods. The dynamic method estimates the position and veloc-ity of a satellite at one epoch by extending the obser-vational arc in order to include more observations of different epochs in it. It is realized by integrating the dynamic equation of the satellite so that the observa-tions at different epochs are related to the particular state. The dy…     

GPS卫星轨道数值积分与广播星历及IGS精密星历的比较

本文采用作者自编的SPPORB IT程序,对GPS卫星轨道的运动方程进行Adam s数值积分求解,同时利用广播星历计算卫星轨道坐标,然后将两者结果同IGS精密星历提供的卫星坐标进行比较,并探讨其轨道误差,计算结果显示广播星历与精密星历差值在2m左右,而数值积分与精密星历的差值在2 cm左右,进一步的分析表明前者误差较大是没有考虑卫星所受的太阳光压、日月引力等影响,而后者考虑了这些影响。鉴于IGS提供的是地固系坐标,而本文数值积分是在惯性系坐标系下进行的,因此本文还举例对惯性坐标系和地固系之间的坐标转换进行了描述。最后,通过实例说明SPPORB IT程序的稳定性以及Adam s数值积分方法的有效性。