一种统计意义下最优接近轨道快速搜索算法

设计空间监视卫星时,需要分析设计卫星与大量被监视空间飞行器的接近能力。通过具体仿真实例提出了一种基于摄动模型的快速最优接近轨道搜索方法,可以大大提高仿真速度,仿真结果证实了该方法的有效性,统计参数通过了精确模型的确认,也验证了算法的可信性。     

基于双星系统定轨的数学方法

通过对双星系统定轨原理的分析,结合卫星的动力学模型,采用联合多时刻观测系统的数据融合方法,得到了一种解算卫星轨道状态的解析方法,同时还给出了误差传递方式。该方法不需要标称轨道,计算方便。理论分析和仿真计算表明,方法确实可行,且定轨精度有所提高。     

Boundary layer aerosol retrieval from thermal infrared nadir sounding – Preliminary results

A non-empirical algorithm is presented to retrieve the optical depth in the 750–1250 cm⁻¹ spectral range, of aerosol located in the boundary layer over the ocean, from nadir high-resolution radiance spectra in the thermal infrared. The algorithm is based on a line-by-line radiative transfer forward model and used the Optimal Estimation Method for the retrieval. Its performance strongly depends on the quality of the a priori temperature and H₂O atmospheric profiles. To demonstrate the relevance of the algorithm, distributions of maritime aerosol parameters have been retrieved from IMG/ADEOS data for December 1996, using the algorithm with the LBLRTM radiative transfer code, and ERA40 (ECMWF) a priori atmospheric profiles and surface conditions.     

Hard X-Ray Modulation Telescope(HXMT) Mission

The HXMT mission concept consists of a slat-collimated hard X-ray detector assembly sensitive in 20～250 keV with a collection area of about 5000 cm2. Based on the reconstruction technique by direct demodulation developed in recent years, HXMT is mainly devoted to performing a hard X-ray all-sky imaging survey with both high sensitivity and high spatial resolution. It can also be used to make pointed observations of X-ray sources to study their spectroscopic and temporal properties in details. The main detector of HXMT consists of 18 individual cylindrical NaI(T1)/CsI(Na) phoswich modules, each with anarea of 283.5 cm2 and a field of view of 5.7°× 1.1° (FWHM). Its spatial resolution and position accuracy are 5′ and 1′ by using the direct demodulation in 1994, and in 2000 its feasibility and technical demonstration study was selected as a project under the Major State Basic Research Program of China. In October 2005, this project entered the full design phase and was listed as a candidate for the first dedicated astronomy satellite around 2010. We are now also considering secondary low energy instruments for this satellite.     

卫星测控天线日照温度场极端工况预测方法

针对卫星测控天线日照温度场产生的热变形影响天线电性能问题,提出了一种预测天线极端温度工况的方法.该方法将太阳照射参数与天线面板方程有效结合起来,求得天线面板的太阳直射点运动轨迹,从而确定天线姿态参数.在Ⅰ-DEAS中建立该姿态的模型,同时引入环境因素的影响,进行瞬态仿真分析,预测天线的极端温度工况.通过试验天线,验证理论计算太阳直射点运动轨迹与极端温度预测方法的正确性.     

Optimal H∞ robust output feedback control for satellite formation in arbitrary elliptical reference orbits

A two degree-of-freedom signal-based optimal H_∞ robust output feedback controller is designed for satellite formation in an arbitrary elliptical reference orbit. Based on high-fidelity linearized dynamics of relative motion, uncertainties introduced by non-zero eccentricity and gravitational J2 perturbation are separated to construct a robust control model. Furthermore, a distributed robust control model is derived by modifying the perturbed robust control model of each satellite with the eigenvalues of the Laplacian matrix of the communication graph, which represent uncertainty in the communication topology. A signal-based optimal H_∞ robust controller is then designed primarily. Considering that the uncertainties involved in the distributed robust control model have a completely diagonal structure, the corresponding analyses are made through structured singular value theory to reduce the conservativeness. Based on simulation results, further designs including increasing the degrees of freedom of the controller, modifying the performance and control weighted functions, adding a post high-pass filter according to the dynamic characteristics, and reducing the control model are made to improve the control performance. Nonlinear simulations demonstrate that the resultant optimal H_∞ robust output feedback controller satisfies the robust performance requirements under uncertainties caused by non-zero eccentricity, J2 perturbation, and varying communication topology, and that 5 m accuracy in terms of stable desired formation configuration can be achieved by the presented optimal H_∞ robust controller. In addition to considering the widely discussed uncertainties caused by the orbit of each satellite in a formation, the optimal H_∞ robust output feedback control model presented in the current work considers the uncertainties caused by varying communication topology in the satellite formation that works in a cooperative way. Other new improvements include adopting a new method to more accurately describe and analyze the effects of the higher-order J2 perturbation, combining all the uncertainties into a diagonal structure, and utilizing a structured singular value to synthesize and analyze the controller.     

Generation of satellite tracking profile: Problems and validation algorithms

Recently, as a satellite mission becomes complicated, it has been required to generate the schedule of satellite antenna movements automatically without relying upon operator’s ad hoc knowledge. To generate the satellite antenna schedule autonomously, this paper first addresses geometrical problems associated with the antenna scheduling and mission planning problems that can be formulated from satellite navigation and antenna orientation information. Then, based on the solutions of the geometrical problems, a set of antenna azimuth and elevation angles that enables the antenna to point towards the desired ground station is obtained systematically. Using the computed azimuth and elevation angles, the satellite tracking profile (TP) is generated, and to validate it, TP validation algorithms are developed.     

Precise orbit determination using the batch filter based on particle filtering with genetic resampling approach

In this study, genetic resampling (GRS) approach is utilized for precise orbit determination (POD) using the batch filter based on particle filtering (PF). Two genetic operations, which are arithmetic crossover and residual mutation, are used for GRS of the batch filter based on PF (PF batch filter). For POD, Laser-ranging Precise Orbit Determination System (LPODS) and satellite laser ranging (SLR) observations of the CHAMP satellite are used. Monte Carlo trials for POD are performed by one hundred times. The characteristics of the POD results by PF batch filter with GRS are compared with those of a PF batch filter with minimum residual resampling (MRRS). The post-fit residual, 3D error by external orbit comparison, and POD repeatability are analyzed for orbit quality assessments. The POD results are externally checked by NASA JPL’s orbits using totally different software, measurements, and techniques. For post-fit residuals and 3D errors, both MRRS and GRS give accurate estimation results whose mean root mean square (RMS) values are at a level of 5 cm and 10–13 cm, respectively. The mean radial orbit errors of both methods are at a level of 5 cm. For POD repeatability represented as the standard deviations of post-fit residuals and 3D errors by repetitive PODs, however, GRS yields 25% and 13% more robust estimation results than MRRS for post-fit residual and 3D error, respectively. This study shows that PF batch filter with GRS approach using genetic operations is superior to PF batch filter with MRRS in terms of robustness in POD with SLR observations.     

On the effects of each term of the geopotential perturbation along the time I: Quasi-circular orbits

This paper provides a useful new method to determine minimum and maximum range of values for the degree and order of the geopotential coefficients required for simulations of orbits of satellites around the Earth. The method consists in a time integration of the perturbing acceleration coming from each harmonic of the geopotential during a time interval T. More precisely, this integral represents the total velocity contribution of a specific harmonic during the period T  . Therefore, for a pre-fixed minimum contribution, for instance 1×10^-8$1\times {10}^{-8}$ m/s during the period of time T, any harmonic whose contribution is below this value can, safely, be neglected. This fact includes some constraints in the degree and order of the terms which are present in the geopotential formula, saving computational efforts compared to the integration of the full model. The advantage of this method is the consideration of other perturbations in the dynamics (we consider the perturbations of the Sun, the Moon, and the direct solar radiation pressure with eclipses), since these forces affect the value of the perturbation of the geopotential, because these perturbations depend on the trajectory of the spacecraft, that is dependent on the dynamical model used. In this paper, we work with quasi-circular orbits and we present several simulations showing the bounds for the maximum degree and order (M) that should be used in the geopotential for different situations, e. g., for a satellite near 500 km of altitude (like the GRACE satellites at the beginning of their mission) we found 35?M?198$35?M?198$ for T=1$T=1$ day. We analyzed the individual contribution of the second order harmonic (_J2$J_2$) and we use its behavior as a parameter to determine the lower limit of the number of terms of the geopotential model. In order to test the accuracy of our truncated model, we calculate the mean squared error between this truncated model and the “full” model, using the CBERS (China-Brazil Earth Resources Satellite) satellite in this test.     

A new approach to the telescope operation method for satellite tracking using a time synchronization technique

For the development of a telescope that is capable of precisely tracking satellites and high-speed operation such as satellite laser ranging, a special method of telescope operation is required. This study aims to propose a new telescope operation method and system configuration for the independent development of a mount and an operation system which includes the host computer. Considering that the tracking of a satellite is performed in real time, communication and synchronization between the two independent subsystems are important. Therefore, this study applied the concept of time synchronization, which is used in various fields of industry, to the communication between the command computer and the mount. In this case, communication delays do not need to be considered in general, and it is possible to cope with data loss. Above all, when the mount is replaced in the future, only the general communication interface needs to be modified, and thus, it is not limited by replacement in terms of the overall system management. The performance of the telescope operation method developed in this study was verified by applying the method to the first mobile SLR system in Korea. This study is significant in that it proposed a new operation method and system configuration, to which the concept of time synchronization was applied, for the observation system that requires an optical telescope.