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.     

On bridges B(pL, qS) around triangular libration points

When μ is smaller than Routh’s critical value μ ₁ = 0.03852 . . . , two planar Lyapunov families around triangular libration points exist, with the names of long and short period families. There are periodic families which we call bridges connecting these two Lyapunov families. With μ increasing from 0 to 1, how these bridges evolve was studied. The interval (0,1) was divided into six subintervals (0, μ ₅), (μ ₅, μ ₄), (μ ₄, μ ₃), (μ ₃, μ ₂), (μ ₂, μ ₁), (μ ₁, 1), and in each subinterval the families B(pL, qS) were studied, along with the families B(qS, qS′). Especially in the interval (μ ₂, μ ₁), the conclusion that the bridges B(qS, qS′) do not exist was obtained. Connections between the short period family and the bridges B(kS, (k + 1)S) were also studied. With these studies, the structure of the web of periodic families around triangular libration points was enriched.     

两行根数辅助的SLR单站定轨

单站测距资料定轨的困难限制了漫反射SLR(Satellite Laser Ranging)测距资料的应用.为此,提出利用两行根数模拟多站SLR测距资料作为辅助,实现单站SLR测距资料定轨的方法.该方法对卫星Ajisai单站SLR测距资料定轨并生成5 d预报轨道,误差小于40 m,实现利用单站测距资料的轨道改进,验证了方法的可行性.