ASO-S卫星工程科学应用系统的数据库设计

先进天基太阳天文台(Advanced Space-based Solar Observatory, ASO-S)卫星是我国首颗太阳观测卫星, 主要观测太阳耀斑和日冕物质抛射以及产生它们的磁场结构. ASO-S卫星的科学应用系统是科学卫星工程的6大系统之一, 它连接科学用户和卫星数据, 为将卫星的科学数据转化为科学成果提供保障. 科学应用系统的数据库是连接软件与海量数据的枢纽, 为科学数据生产和用户服务及运行提供数据层的支撑. 介绍了科学应用系统的数据库架构设计、数据库的选择以及数据库性能优化和表样例. 这里的数据库包括观测计划、工程参数、运维日志、科学数据、定标数据和特征事件识别等数据库. 这些数据库的建设将为ASO-S卫星工程科学应用系统的顺利运行提供数据支撑, 也可以为未来其他科学卫星类似数据库的搭建提供参考和借鉴.     

遗传算法优化的BP神经网络卫星钟差预报

针对BP (Back Propagation)神经网络模型预测卫星钟差中权值和阈值的最优化问题, 提出了基于遗传算法优化的BP神经网络卫星钟差短期预报模型, 给出了遗传算法优化BP神经网络的基本思想、具体方法和实施步骤. 为验证该优化模型的有效性和可行性, 利用北斗卫星导航系统(BeiDou navigation satellite system, BDS)卫星钟差数据进行钟差预报精度分析, 并将其与灰色模型(GM(1,1))和BP神经网络模型预报的结果比较分析. 结果表明: 该模型在短期钟差预报中具有较好的精度, 优于GM(1,1)模型和BP神经网络模型.     

基于斩波轮技术的K波段接收机噪声校准研究

接收机是射电天文中用于探测微弱射电信号的重要接收设备.接收机的强度校准就是将接收机对射电源的响应转换为天文意义上的流量密度.常规方法就是使用经典的冷热负载法,将接收机自身的强度响应转换为一个等效的温度值,之后再据此对射电源做进一步标定.通过搭建基于斩波轮技术的K波段接收机强度校准平台,使用斩波轮法测试K波段常温接收机的噪声温度,并与传统冷热负载法的测试结果进行比对.结果显示,在晴好天气条件下,斩波轮法在30°、90°仰角下噪声温度的最大测试误差为7.5%和8.4%,可以很好地应用于实际噪声温度测试中;但在5°仰角测试中,由于过低仰角引入了地面噪声,使得斩波轮法的测试误差上升至20%–30%之间而无法使用.希望在此基础上进一步开展K波段天空亮温度的理论计算与实测,从而完善斩波轮技术的应用,使之可以满足在不同气象条件下的噪声校准测试需求.     

基于Stacking集成学习的恒星/星系分类研究

机器学习在当今诸多领域已经取得了巨大的成功,但是机器学习的预测效果往往依赖于具体问题.集成学习通过综合多个基分类器来预测结果,因此,其适应各种场景的能力较强,分类准确率较高.基于斯隆数字巡天(Sloan Digital Sky Survey,SDSS)计划恒星/星系中最暗源星等集分类正确率低的问题,提出一种基于Stacking集成学习的恒星/星系分类算法.从SDSS-DR7(SDSS Data Release 7)中获取完整的测光数据集,并根据星等值划分为亮源星等集、暗源星等集和最暗源星等集.仅针对分类较为复杂且困难的最暗源星等集展开分类研究.首先,对最暗源星等集使用10折嵌套交叉验证,然后使用支持向量机(Support Vector Machine,SVM)、随机森林(Random Forest,RF)、XGBoost(eXtreme Gradient Boosting)等算法建立基分类器模型;使用梯度提升树(Gradient Boosting Decision Tree,GBDT)作为元分类器模型.最后,使用基于星系的分类正确率等指标,与功能树(Function Tree,FT)、SVM、RF、GBDT、XGBoost、堆叠降噪自编码(Stacked Denoising AutoEncoders,SDAE)、深度置信网络(Deep Belief Network,DBN)、深度感知决策树(Deep Perception Decision Tree,DPDT)等模型进行分类结果对比分析.实验结果表明,Stacking集成学习模型在最暗源星等集分类中要比FT算法的星系分类正确率提高了将近10%.同其他传统的机器学习算法、较强的提升算法、深度学习算法相比,Stacking集成学习模型也有较大的提升.     

An Automatic Method for Detecting Transients and Variable Sources in AST3 Survey Based on Image Subtraction and Random Forest

AST3-2 (the second Antarctic Survey Telescope) is located in Antarctic Dome A, the loftiest ice dome on the Antarctic Plateau. It produces a huge amount of observational data which require a more efficient data reduction program to be developed. Also the data transmission in Antarctica is much difficult, thus it is necessary to perform data reduction and detect variable and transient sources remotely and automatically in Antarctica, but this attempt is restricted by the unsatisfactory performance of the low power consumption computer in Antarctica. For realizing this purpose, to develop a new method based on the existing image subtraction method and random forest algorithm, taking the AST3-2 2016 dataset as the test sample, becomes an alternative choice. This method performs image subtraction on the dataset, then applies the principle component analysis to extract the features of residual images. Random forest is used as a machine learning classifier, and in the test a recall rate of 97% is resulted for the positive sample. Our work has verified the feasibility and accuracy of this method, and finally found out a batch of candidates for variable stars in the AST3-2 2016 dataset.     

Formation of Carbon-Enhanced Metal-Poor RR Lyrae Stars

Carbon-enhanced metal-poor (CEMP) stars are considered to be related to the first generation of stars, and responsible for the chemical evolution of the early Galaxy. More than half of them are in binaries, and could be explained by the binary evolution, but the formation channel of them is still not fully understood. Among the hundreds of CEMP stars, there are nine CEMP RR Lyrae stars identified, and at least seven of which are very likely not binaries. The usual binary star evolution channel is difficult to produce such a single star, particularly that of carbon enrichment. One way in which such a single star might be produced is the merger of a helium white dwarf with a Hertzsprung gap (HG) star. We use a stellar evolution program to calculate the models of the merger remnants, and find that the models can reproduce the observed distribution of these CEMP single RR Lyrae stars in terms of surface temperature, gravity, and carbon abundance. Hence, it is extremely possible that the helium white dwarf and HG star merger model is one of the formation channels of the metal-poor carbon-rich RR Lyrae stars.     

On the “rotational angular momentum” of the oceans and the corresponding polar motion

Periodic polar motions caused by ocean tides are predicted. In the Liouville equations for rotational motion the complete excitation functions for the ocean tides have to be used. This does not depend on the fact that hydrodynamical ocean tide models do not consider the centrifugal acceleration. The observable polar motion of the Celestial Ephemeris Pole CEP (more exactly: the terrestrial location of the CEP) is tabulated for the ten ocean tides M₂, S₂, N₂, K₁, O₁, P₁, M f, M f′, M m, Ssa. Typical amplitudes for the largest ocean tides are 0.4 milliarcseconds. This is within the reach of geodetic VLBI and SLR observations.     

A multiple shooting approach for the numerical treatment of stellar structure and evolution

We present a new numerical method for solving the system of partial differential equations describing the structure and evolution of a spherically symmetric star. As usual, we employ the transversal method of lines in order to split the equations into a coupled spatial and temporal part. The novel features of the algorithm are the following: (a) Instead of using the Lagrangian picture we formulate the system of partial differential equations in the Eulerian picture. (b) We reformulate the equations of stellar structure as a multipoint boundary-value problem. By means of this reformulation the rather clumsy iterative matching procedure of stellar atmosphere and interior is avoided. (c) The multipoint boundary-value problem is solved by the multiple shooting method. This approach not only ensures a high accuracy of the stellar models calculated at each time step but also allows the free boundaries inside the star due to different energy transport mechanisms to be located exactly. (d) The time derivatives involved in the stellar-structure equations are discretized implicitly to second order accuracy. Moreover, at each time step, the chemical abundances are determined by using a sophisticated update procedure. In this way, a high accuracy is achieved with respect to the integration in time. The algorithm has turned out to be exceedingly reliable and numerically accurate. This is shown by the evolution of a 1 M_⊙ star up to the hydrogen-shell burning phase. In this example, the virial theorem, the law of mass conservation, and the law of energy conservation is fulfilled to a hitherto unattainable degree of accuracy. Since the multiple shooting method, which is at the heart of our approach, is a perfect example of a parallel algorithm, the computational speed of the algorithm might be substantially improved provided easy-to-program, high-performance parallel computers with sufficiently many processors become available in the near future.     

The outer core rotation about an inclined axis from geomagnetic secular variations (Mitteilungen des Zentralinstituts für Physik der Erde,Nr. 1832)

Non-axisymmetric motions of the outer core of the Earth are important for the dynamo problem and the excitation of the decade variations of the polar motion. The components of the vector of a rigid rotation of the outer core about an inclined axis were estimated by a first-order approximation of the frozen- field theory of the geomagnetic secular variation from 1903.5 to 1975.5. The trends and quasi-periodic constituents of these quantities were computed. It was shown that the position and time behaviour of the rotational pole of the outer core differ considerably from the well-know co-ordinates of the dipole axis. Some periods of the equatorial components of the rotational vector are comparable with those of the axial component previously derived for a pure axial rotation. Additionally, the time behaviour of the pole path shows events like the well-known Markowitz wobble but naturally with other extent. These and other results suggested that the investigations are worth to be continued in future by some physical interpretations.