Optimal two-impulse rendezvous using constrained multiple-revolution Lambert solutions

A solution to the fixed-time minimum-fuel two-impulse rendezvous problem for the general non-coplanar elliptical orbits is provided. The optimal transfer orbit is obtained using the constrained multiple-revolution Lambert solution. Constraints consist of lower bound for perigee altitude and upper bound for apogee altitude. The optimal time-free two-impulse transfer problem between two fixed endpoints implies finding the roots of an eighth order polynomial, which is done using a numerical iterative technique. The set of feasible solutions is determined by using the constraints conditions to solve for the short-path and long-path orbits semimajor axis ranges. Then, by comparing the optimal time-free solution with the feasible solutions, the optimal semimajor axis for the two fixed-endpoints transfer is identified. Based on the proposed solution procedure for the optimal two fixed-endpoints transfer, a contour of the minimum cost for different initial and final coasting parameters is obtained. Finally, a numerical optimization algorithm (e.g., evolutionary algorithm) can be used to solve this global minimization problem. A numerical example is provided to show how to apply the proposed technique.     

Research on the Method of Orbit Determination Based on the Self-adaptive Stable Least p-norm Estimate

The traditional least square estimation (LSE) method for orbit determination will not be optimal if the error of observational data does not obey the Gaussian distribution. In order to solve this problem, the least p-norm (Lp) estimation method is presented in this paper to deal with the non-Gaussian distribution cases. We show that a suitable selection of parameter p may guarantee a reasonable orbit determination result. The character of Lp estimation is analyzed. It is shown that the traditional Lp estimation method is not a robust method. And a stable Lp estimating based on data depth weighting is put forward to deal with the model error and outlier. In the orbit determination process, the outlier of observational data and coarse model error can be quantitatively described by their weights. The farther is the data from the data center, the smaller is the value of data depth and the smaller is the weighted value accordingly. The result of the new Lp method is stabler than that of the traditional Lp estimation and the breakdown point could be up to 1/2. In addition, the orbit parameter is adaptively estimated by residual analysis and matrix estimation method, and the estimation efficiency is enhanced. Finally, by taking the Space-based Space Surveillance System as an example and performing simulation experiments, we show that if there are system error or abnormal value in the observational data or system error in satellite dynamical model and space-based observation platform, LSE will not be optimal even though the observational data obeys the Gaussian distribution, and the orbit determination precision by the self-adaptive robust Lp estimation method is much better than that by the traditional LSE method.     

Classification of stellar spectra with SVM based on within-class scatter and between-class scatter

Support Vector Machine (SVM) is a popular data mining technique, and it has been widely applied in astronomical tasks, especially in stellar spectra classification. Since SVM doesn’t take the data distribution into consideration, and therefore, its classification efficiencies can’t be greatly improved. Meanwhile, SVM ignores the internal information of the training dataset, such as the within-class structure and between-class structure. In view of this, we propose a new classification algorithm-SVM based on Within-Class Scatter and Between-Class Scatter (WBS-SVM) in this paper. WBS-SVM tries to find an optimal hyperplane to separate two classes. The difference is that it incorporates minimum within-class scatter and maximum between-class scatter in Linear Discriminant Analysis (LDA) into SVM. These two scatters represent the distributions of the training dataset, and the optimization of WBS-SVM ensures the samples in the same class are as close as possible and the samples in different classes are as far as possible. Experiments on the K-, F-, G-type stellar spectra from Sloan Digital Sky Survey (SDSS), Data Release 8 show that our proposed WBS-SVM can greatly improve the classification accuracies.     

The Unusual Absorption Line Spectrum of Quasar SDSS J2220+0109

We report the Balmer broad absorption lines (BALs) in the quasar SDSS J2220 + 0109 discovered from the SDSS data, and present a detailed analysis of the peculiar absorption line spectrum, including the He I* multiplet at λλ3189, 3889 arising from the metastable 2³s-state helium and the Balmer Hα and Hβ lines from the excited hydrogen H I of n = 2 level, which are rarely seen in quasar spectra, as well as many absorption lines arising from the excited Fe II* of the levels 7 955 cm⁻¹, 13 474 cm⁻¹ and 13 673 cm⁻¹ in the wavelength range 3100∼3300 Å. Ca II H, K absorption line doublets also clearly appear in the SDSS spectrum. All absorption lines show a similar blueshifted velocity structure of Δv ≈ − 1500 ∼ 0 km·s⁻¹ relative to the quasar's systematic redshift determined from the emission lines. Detailed analysis suggests that the Balmer absorption lines should arise from the partially ionized region with a column density of N_HI ≈ 10²¹ cm⁻² for an electron density of ne ∼ 10⁶ cm⁻³; and that the hydrogen n = 2 level may be populated via collisional excitation with Lyα pumping.     

Temporal and spatial distribution of GPS-TEC anomalies prior to the strong earthquakes

Earthquakes are one of the most destructive and harmful natural disasters, especially in recent years, the 2008/5/12 Wenchuan M7.9 earthquake, the 2011/3/11 Tohoku M9.0 earthquake and the 2012/4/11 Sumatra M8.6 earthquake have caused a significant impact to the human life. In this paper, we make a study of the temporal and spatial distribution of the Global Positioning System Total Electron Content (GPS TEC) anomalies prior to the three strong earthquakes by the method of statistical analysis. Our results show that the pre-earthquake ionospheric anomalies are mainly positive anomalies and take the shape of a double-crest structure with a trough near the epicenter. The ionospheric anomalies do not coincide with the vertical projection of the epicenter of the subsequent earthquake, but mainly localize in the near-epicenter region and corresponding ionospheric anomalies are also simultaneously observed in the magnetic conjugate region prior to the three earthquakes. In addition, the amplitude and scale-size of the ionospheric ΔTEC are different with the magnitude of the earthquake, and the horizontal scale-size of the greatest anomalies before the Tohoku M9.0 earthquake is ～30^° in longitude and ～10^° in latitude, with the maximum amplitude of TEC disturbances reaching ～20 TECu relative to the background. The peak of anomaly enhancement usually occurs in the afternoon to sunset (i.e. between 14:00 and 18:00 local time) which lasts for approximate 2 hours. Possible causes of these anomalies are discussed, and after eliminating the effect of solar activities and magnetic storms it can be concluded that the detected obvious and regular anomalous behavior in TEC within just a few days before the earthquakes is related with the forthcoming earthquakes with high probability.     

Simulation of evolutionary mechanism of planetary rings

We devised a numerical model of planetary ring in which the inelastic collision and the gravitation between ring particles are considered. We adopt Hill's equations and the differential algorithm of circular mesh for computation of the particle orbits. The evolutional processes are presented for different coefficient of restitution and dynamical optical depth . The results show that the semi-major axis and eccentricity of the ring particles are changed with and . We compute the average energies transferred and loss in inelastic collisions for various values of the parameters. The dynamical equilibrium properties are discussed in the different cases.     

NGC 7538—IRS1致密HⅡ区研究

本文在分析研究NGC7538-IRS1致密HⅡ区H_2CO和OH脉泽辐射VLBI观测结果的基础上,指出该HⅡ区合理的模型是:HⅡ区表面为厚的尘埃层包围,尘埃层两极已被突破,并形成双极流;HⅡ区外面有一个环形转动气体-尘埃云,存在由环向HⅡ区表面的物质下落;包括环和HⅡ区在内的整个系统视向速度为-61km/s,该系统居于视向速度为-57km/s的更大分子云中。H_2CO和OH脉泽发生在HⅡ区两极附近离HⅡ区表面小于0.2R_(HⅡ)的区域内。利用上述模型,还讨论了H_2O脉泽及其他分子吸收线和发射线的发生区域。     

W3(OH)致密HⅡ区磁场研究

本文以W3(OH)为实例,建立了与OH脉泽成协的HII区气体-尘埃壳层的磁场模型,由此可见,致密HII区所在分子云核心中的磁场强度,与分子云核心的分子数密度之间存在指数α=1/2的幂律关系;而在HII区气体-尘埃壳层中,磁场强度与分子数密度之间存在α=1的正比关系,根据讨论可知,与OH脉泽成协的HII区,其气体-尘埃壳层的分子数密度为10~6cm~(-3)量级,磁场强度为几个mG;一旦分子数密度达10~7cm~(-3)量级,OH脉泽便将熄灭。