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

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

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数值积分方法的有效性。     

利用C波段转发数据实现GEO卫星轨道确定

为了提高我国卫星导航系统GEO卫星定轨能力,针对我国卫星导航系统所特有的新数据和新方法,利用CAPS系统的C波段转发式测距数据进行了单GEO卫星定轨实验.在分析该方法技术特点的基础上,从定轨残差和轨道重叠弧段比较等方面,对定轨结果进行了分析.由于设备转发时延是影响该方法定轨精度的重要因素,因此同时还分析了系统偏差对定轨精度的影响,得出了一些有益的结论.     

Fog Detection over China’s Adjacent Sea Area by using the MTSAT Geostationary Satellite Data

A fog threshold method for the detection of sea fog from Multi-function Transport Satellite (MTSAT1R) infrared (IR) channel data is presented.This method uses principle component analysis (PCA),texture analysis,and threshold detection to extract sea fog information.A heavy sea fog episode that occurred over China’s adjacent sea area during 7 8 April 2008 was detected,indicating that the fog threshold method can effectively detect sea fog areas nearly 24 hours a day.MTSAT-1R data from March 2006,June 2007,and April 2008 were processed using the fog threshold method,and sea fog coverage information was compared with the meteorological observation report data from ships.The hit rate,miss rate,and false alarm rate of sea fog detection were 66.1%,27.3%,and 33.9%,respectively.The results show that the fog threshold method can detect the formation,evolution,and dissipation of sea fog events over period of time and that the method has superior temporal and spatial resolution relative to conventional ship observations.In addition,through MTSAT-1R data processing and a statistical analysis of sea fog coverage information for the period from 2006 to 2009,the monthly mean sea fog day frequency,spatial distribution and seasonal variation characteristics of sea fog over China’s adjacent sea area were obtained.     

The Data Depth-weight-kernel Estimation of the Model Error of Satellite Orbit Determination

It is an objective fact that there exists error in the satellite dynamic model and it will be transferred to satellite orbit determination algorithm, forming a part of the connotative model error. Mixed with the systematic error and random error of the measurements, they form the unitive model error and badly restrict the precision of the orbit determination. We deduce in detail the equations of orbit improvement for a system with dynamic model error, construct the parametric model for the explicit part of the model and nonparametric model for the error that can not be explicitly described. We also construct the partially linear orbit determination model, estimate and fit the model error using a two-stage estimation and a kernel function estimation, and finally make the corresponding compensation in the orbit determination. Beginning from the data depth theory, a data depth weight kernel estimator for model error is proposed for the sake of promoting the steadiness of model error estimation. Simulation experiments of SBSS are performed. The results show clearly that the model error is one of the most important effects that will influence the precision of the orbit determination. The kernel function method can effectively estimate the model error, with the window width as a major restrict parameter. A data depth-weight-kernel estimation, however, can improve largely the robustness of the kernel function and therefore improve the precision of orbit determination.