基于自适应稳健最小p范数估计的轨道确定方法研究

当测轨数据误差不服从正态分布时,传统的最小二乘(LSE)轨道确定方法将不是最优的.为了获得高精度的定轨结果,一种可行的策略是采用基于最小p范数(Lp)的轨道确定方法.通过分析Lp估计的相关性质,得出普通Lp估计不具有良好的抗差性的结论.为抑制模型误差和异常值的影响,提出了基于数据深度加权的稳健最小p范数估计方法,并证明了相关性质,得出了其崩溃点可以达到1/2的结论.最后,通过残差分析和矩估计法自适应估计相关参数,使得估计达到最大效率.以天基空间目标监视系统为背景进行了仿真试验.结果表明,当观测数据存在系统误差或异常值时,或者当目标动力学模型存在误差或者天基观测平台存在系统误差时,即使观测数据服从正态分布,LSE也不是最优的,在这种意义下自适应稳健Lp估计轨道确定方法比传统轨道确定方法更加稳健,定轨精度也更高.     

抗差估计在星载GPS卫星非差运动学定轨中的应用

针对星载GPS卫星非差运动学定轨的特点,提出将抗差估计应用于星载GPS卫星运动学定轨,一方面尽量保留来自每颗GPS卫星的观测值,保持较强的卫星几何强度,避免轨道求解出现奇异;另一方面,采用等价权思想,有效地降低质量较差的观测值对定轨结果的影响,保证定轨精度．还采用CHAMP卫星的实测GPS数据验证了新方案的可行性和有效性．     

长春站非合作目标激光测距资料的定轨

精密解算了非合作目标的单站激光测距数据.观测数据少、数据弧段分布不好是对非合作目标进行精密定轨的难点.通过定轨过程中对动力学模型的选择及求解参数的选取,使得轨道计算收敛.解算多组圈数的非合作目标数据,将轨道重叠弧段对比作为评判定轨精度的指标;从多组圈数中提取出一圈的观测数据,对其余数据进行定轨处理,将定轨后的轨道结果与提取出的观测数据进行对比,得到在同一时刻的距离偏差,使其作为精密定轨的外符合.结果表明:对非合作目标(4814)进行精密定轨,平均测距残差为1.01 m,在测距方向上,测量数据外符合的平均轨道精度为14.35 m,预报1 d的测距精度为24.60 m.     

GIM在LEO卫星单频GPS定轨中的应用

电离层延迟误差是单频GPS(Global Positioning System)数据最主要的误差源,为提高基于单频GPS数据的LEO(Low Earth Orbiting)卫星定轨精度,必须消除/减弱GPS观测数据中电离层延迟影响.研究了全球电离层模型GIM(Global IonosphericMaps)在基于单频GPS伪距数据的低轨卫星运动学和动力学定轨中的应用,并通过估算电离层尺度因子的方法消除C/A码伪距观测量中电离层延迟影响.由于LEO卫星星载GPS信号受电离层延迟影响与卫星轨道高度相关,选取了轨道高度在300～800 km的CHAMP(CHAllenging Mini-satellite Payload)、GRACE(Gravity Recovery AndClimate Experiment)、TerraSAR-X及SAC-C等LEO卫星C/A码伪距观测量作为试算数据.CHAMP等卫星实测数据计算结果表明:以JPL(Jet Propulsion Laboratory)发布的GIM模型作为背景模型,通过电离层比例因子法能很好地消除C/A码伪距观测量中电离层延迟影响,提高LEO卫星运动学和动力学定轨精度,其中,CHAMP卫星轨道最低,受电离层延迟影响最严重,定轨精度提高最显著,分别为55.6%和47.6%;SAC-C卫星轨道高度最高,受电离层延迟影响最小,相应的定轨精度提高幅度也最低,分别为47.8%和38.2%.     

双星定位系统的近地卫星联合定轨同质异质观测数据加权方法及应用

针对双星定位系统的近地卫星联合定轨中的多源观测数据的融合处理问题,建立了同质观测数据的二步系统误差修正的改进的方差分量估计最优加权方法;分析指出异质观测数据的多源融合测量模型本质为多结构多参数的非线性回归模型,建立了异质观测数据的模型结构特征分析和方差分量估计相结合的最优加权方法.设计了两类观测数据最优加权及联合定轨参数估计的实现算法,并以双星及备份星的距离和同质观测数据以及双星距离和与星敏感器测角的异质观测数据为例,进行了联合定轨仿真实验.理论分析和仿真计算结果表明;对于同质观测数据联合定轨,采用二步系统误差修正的方差分量估计法,可以获得比传统的经验加权算法更优的定轨精度;对于异质观测数据联合定轨,通过引入表征模型结构特征的加权因子,与平均加权方式相比,近地卫星及静地卫星的联合定轨精度均得到一定程度的改善.     

双行根数编目体系轨道误差研究

大量空间目标的真实轨道无法精确知道,目前只能通过跟踪观测的数据进行定轨来得到估计轨道,而估计的轨道就会有误差.双行根数(TLE)是广泛使用的一种特殊编目轨道根数,其配套的轨道模型为Simplified General Perturbations 4(SGP4)/Simplified Deep-space Perturbations 4(SDP4)模型.编目轨道的精度主要依赖于相应的观测模型和动力学模型,这些模型一般都不会非常准确,往往会有误差,有些误差可能直接导致编目定轨结果在局部为有偏估计.通过理论研究和仿真模拟,分析了动力学模型中地球非球形引力位田谐摄动项对编目轨道精度的影响,发现TLE编目轨道中存在随时间周期变化的系统差,该系统误差甚至可以达到千米量级.几何构型较好的测站分布在一定程度上可以削弱编目定轨中产生的系统误差,由于力学模型的限制,无法彻底消除.     

电离层延时对低轨卫星频率传递的影响分析

利用低轨卫星,对受干扰区域的卫星导航接收机进行频率传递,是提升其性能的一种有效手段.而电离层延时误差是影响传递效果的重要因素.对基于低轨卫星的频率传递应用背景进行了介绍,分析了其基本原理,并着重分析了电离层延时的存在对频率传递的影响,提出利用站间差历元差的方法对电离层延时进行修正.最后,通过仿真对理论分析进行了验证.结果表明,通过站间差历元差的手段,对当前电离层变化值进行求解与预测,可以将电离层延时误差的变化控制在一定范围,满足频率传递的要求.     

GPS星历对LEO星载GPS精密定轨精度的影响

目前,越来越多的低轨卫星上都搭载了用于精密定轨的星载GPs接收机,星载GPS已成为低轨卫星精密定轨的主要手段之一.星载GPS精密定轨精度依赖于GPS星历及钟差精度.基于SHORDE-Ⅲ非差动力学定轨功能,以2005年8月1日至8月7日一周的GRACE卫星实测数据为例,采用事后精密轨道(igs)、快速轨道(igr)和超快速轨道(igu)三种GPS星历在同等条件下定轨,估计GPS星历精度对低轨卫星定轨精度的影响,实际计算结果表明igs和igr两类GPS星历定轨精度相当,约为9.5 cm,igu星历定轨精度略低于igs和igr星历,约为10.5cm:高频GPS卫星钟差数据对定轨精度会产生1-6cm影响.     

基于智能优化算法的小天体初轨确定

经典的初轨确定方法包括Laplace方法和Gauss方法以及它们的各种变化形式. 除这些经典方法之外, 基于当今光学观测数据的特点, 学者们也陆续提出了一些其他的初轨确定方法, 包括双r (目标距离观测者的距离)方法和可行域方法. 双r方法的一种实现方式是通过猜测某两个时刻(通常是定轨弧段的首、末时刻)目标离观测者的距离, 结合观测者在空间中的位置矢量, 即可求解相应的Lambert弧段作为目标轨道的初始猜测. 进一步, 以其他观测时刻的RMS (Root Mean Square)为优化变量可以改进初始猜测从而确定初轨. 可行域方法则是针对一组初始观测参数(包括赤经、赤纬及其变率), 根据一些初始假设将目标(离观测者的)距离及其变率约束在可行域内, 并通过三角划分逐步逼近的方式寻找到使观测RMS最小的猜测解. 针对一系列模拟观测数据以及实测数据, 将智能优化算法(粒子群算法)应用于这两种初轨方法, 并将结果与改进的Laplace算法的结果进行比较. 由于双r方法不仅可以用于短弧定轨还可用于长弧关联, 所以进一步给出了针对长弧段数据的关联结果.     

一种初轨计算的稳健方法

光学测角资料的初轨计算在空间目标搜索发现中具有重要作用,当观测资料存在野值时,基于最小二乘的经典初轨计算方法不够稳健.采用最小一乘方法建立了一种初轨计算的稳健方法,方法将初轨计算问题转换为线性规划问题求解,并通过bootstrap方法给出估计精度.数值计算结果表明方法稳健有效,并具有较高的崩溃点.     

A thermodynamic model of high temperature lava vaporization on Io

We modified the MAGMA chemical equilibrium code developed by Fegley and Cameron (1987, Earth Planet. Sci. Lett. 82, 207-222) and used it to model vaporization of high temperature silicate lavas on Io. The MAGMA code computes chemical equilibria in a melt, between melt and its equilibrium vapor, and in the gas phase. The good agreement of MAGMA code results with experimental data and with other computer codes is demonstrated. The temperature-dependent pressure and composition of vapor in equilibrium with lava is calculated from 1700 to 2400 K for 109 different silicate lavas in the ONaKFeSiMgCaAlTi system. Results for five lavas (tholeiitic basalt, alkali basalt, Barberton komatiite, dunite, and a molten type B1 Ca, Al-rich inclusion) are discussed in detail. The effects of continuous fractional vaporization on chemistry of these lavas and their equilibrium vapor are presented. The predicted abundances (relative to Na) of K, Fe, Si, Al, Ca, and Ti in the vapor equilibrated with lavas at 1900 K are lower than published upper limits for Io's atmosphere (which do not include Mg). We predict evaporative loss of alkalis, Fe, and Si during volcanic eruptions. Sodium is more volatile than K, and the Na/K ratio in the gas is decreased by fractional vaporization. This process can match Io's atmospheric Na/K ratio of 10±3 reported by Brown (2001, Icarus 151, 190-195). Silicon monoxide is an abundant species in the vapor above lavas. Spectroscopic searches are recommended for SiO at IR and mm wavelengths. Reactions of metallic vapors with S- and Cl-bearing volcanic gases may form other unusual gases including MgCl₂, MgS, MgCl, FeCl₂, FeS, FeCl, and SiS.     

Instructions for authors

正Research in Astronomy and Astrophysics publishes original research papers and reviews on all branches of astronomy and astrophysics.Reviews are by invitation only.Important new results that require rapid publication can be submitted as a Letter(Letters must be restricted in length to 6 printed pages).Authors who submit a paper are expected to be able to certify that the paper is original work,     

Research in Astronomy and Astrophysics CONTENTS

正~~     

Alkali and halogen chemistry in volcanic gases on Io

We use chemical equilibrium calculations to model the speciation of alkalis and halogens in volcanic gases emitted on Io. The calculations cover wide temperature (500-2000 K) and pressure (10⁻⁶ to 10⁺¹ bars) ranges, which overlap the nominal conditions at Pele (T=1760 K, P=0.01 bars). About 230 compounds of 11 elements (O, S, Li, Na, K, Rb, Cs, F, Cl, Br, I) are considered. The elemental abundances for O, S, Na, K, and Cl are based upon observations. CI chondritic elemental abundances relative to sulfur are used for the other alkalis and halogens (as yet unobserved on Io). We predict the major alkali species in Pele-like volcanic gases and the percentage distribution of each alkali are LiCl (73%), LiF (27%); NaCl (81%), Na (16%), NaF (3%); KCl (91%), K (5%), KF (4%); RbCl (93%), Rb (4%), RbF (3%); CsCl (92%), CsF (6%), Cs (2%). Likewise the major halogen species and the percentage distribution of each halogen are NaF (88%), KF (10%), LiF (2%); NaCl (89%), KCl (11%); NaBr (89%), KBr (10%), Br (1%); NaI (61%), I (30%), KI (9%). We predict the major halogen condensates and their condensation temperatures at P=0.01 bar are NaF (1115 K), LiF (970 K); NaCl (1050 K), KCl (950 K); KBr (750 K), RbBr (730 K), CsBr (645 K); and solid I₂ (200 K). We also model disequilibrium chemistry of the alkalis and halogens in the volcanic plume. Based on this work and our prior modeling for Na, K, and Cl in a volcanic plume, we predict the major loss processes for the alkali halide gases are photolysis and/or condensation onto grains. Their estimated photochemical lifetimes range from a few minutes for alkali iodides to a few hours for alkali fluorides. Condensation is apparently the only loss process for elemental iodine. On the basis of elemental abundances and photochemical lifetimes, we recommend searching for gaseous KCl, NaF, LiF, LiCl, RbF, RbCl, CsF, and CsCl around volcanic vents during eruptions. Based on abundance considerations and observations of brown dwarfs we also recommend a search of Io's extended atmosphere and the Io plasma torus for neutral and ionized Li, Cs, Rb, and F.     

Study on the Stressed Mirror Polishing with a Continuous Polishing Machine for Large Aperture Off-axis Aspheric Mirrors

A special stressed annular polishing technique is proposed to mill the off-axis aspheric sub-mirrors of a large segmented mirror with an annular polishing machine. Based on the basic principle of stressed annular polishing technique, a set of special stressing mechanisms are designed to convert milling the aspheric surfaces of sub-mirrors with different off-axis distances into milling the spherical surfaces with identical radii of curvature, so that they can be pol- ished simultaneously on a continuous polishing machine. It took about contin- uous 40 hours to polish a scaled-down mirror of the planning Chinese Future Giant Telescope (CFGT) using this technique. This mirror has the 330 mm di- ameter, 3.6 m off-axis distance, and the 21.6 m radius of curvature, and its max- imum asphericity is 16 micron. The experiment shows that this method has a high effciency, suits batch manufacturing, especially the batch manufacturing of aspheric sub-mirrors of the segmented primary mirror of an extremely large aperture telescope.     

The climatic and physiographic controls of the eastern Mediterranean over the late Pleistocene climates in the southern Levant and its neighboring deserts

Modern-day synoptic-scale eastern Mediterranean climatology provides a useful context to synthesize the diverse late Pleistocene (60–12 ka) paleohydrologic and paleoenvironmental indicators of past climatic conditions in the Levant and the deserts to its south and east. We first critically evaluate, extract, and summarize paleoenvironmental and paleohydrologic records. Then, we propose a framework of eastern Mediterranean atmospheric circulation features interacting with the morphology and location of the southeast Mediterranean coast. Together they strongly control the spatial distribution of rainfall and wind pattern. This cyclone–physiography interaction enforces the observed rainfall patterns by hampering rainfall generation south and southeast of the latitude of the north Sinai coast, currently at 31°15′.The proposed framework explains the much-increased rains in Lebanon and northern Israel and Jordan as deduced from pollen, rise and maintenance of Lake Lisan, and speleothem formation in areas currently arid and semiarid. The proposed framework also accounts for the southward and eastward transition into semiarid, arid, and hyperarid deserts as expressed in thick loess accumulation at the deserts' margins, dune migration from west to east in the Sinai and the western Negev, and the formation of hyperarid (< 80 mm yr^− 1) gypsic–salic soils in the southern Negev and Sinai. Our climatic synthesis explains the hyperarid condition in the southern Negev, located only 200–250 km south of the much-increased rains in the north, probably reflecting a steeper rainfall gradient than the present-day gradient from the wetter Levant into its bordering southern and eastern deserts.At present, the rainiest winter seasons in Lebanon and northern and central Israel are associated with more frequent (+ 20%), deeper Cyprus Lows traversing the eastern Mediterranean at approximately the latitude of southern Turkey. Even these wettest years in northern Israel do not yield above average annual rainfall amounts in the hyperarid southern Negev. This region is mainly influenced by the Active Red Sea Troughs that produce only localized rains. The eastern Mediterranean Cyprus Lows also produce more dust storms and transport higher amounts of suspended dust to the loess area than any other atmospheric pattern. Concurrent rainfall and dust are essential to the late Pleistocene formation of the elongated thick loess zone along the desert northern margin. Even with existing dust storms, the lack of rain and very sparse vegetation account for the absence of late Pleistocene loess sequences from the southern Negev and the formation of hyperarid soils.When the north Sinai coast shifted 30–70 km northwest due to last glacial global sea level lowering, the newly exposed coastal areas supplied the sand and dust to these active eastern Mediterranean cyclones. This enforced the latitude of the northern boundary of the loess zone to be directly due east of the LGM shoreline. This shift of coast to the northwest inhibited rainfall in the southern Levant deserts and maintained their hyperaridity. Concurrently, frequent deep eastern Mediterranean Cyprus Lows were funneled along the northern Mediterranean increasing (probably doubling) the rains in central and northern Israel, Lebanon, southwestern Syria and northern Jordan. These storms and rains formed lakes, forests, and speleothems only a short distance north of the deserts in the southern Levant.     

Active hot star shells

Interferometry in the visible provides milliarcsecond spatial resolution and thus has been used for studying the circumstellar environment of active hot stars. In this paper I will illustrate how the visibility modulus and phase can be used to better constrain the physics of Be disks through results from the VLA, the MkII and the GI2T interferometers. I will insist on the importance and the potential of coupling high angular resolution with high spectral resolution to the study of Be shells. Finally I will present a possible study of the circumstellar disk of Be stars using the VLTI. This revised version was published online in July 2006 with corrections to the Cover Date.     

Fitting X-ray Afterglow Light curves of Gamma-ray Bursts by Using the Magnetar Energy Injection Model

The central compact object for some gamma-ray bursts (GRBs) may be a strongly magnetized millisecond pulsar. It can inject energy to the outer shock of the GRB by through the magnetic dipole radiation, and therefore causes the shallow decay of the early afterglow. Recently, from a large number of GRB X-ray afterglows observed by Swift/XRT(X-ray telescope), it is revealed that many of them exhibit the shallow decay about 10²∼10⁴ s after the burst prompt emission. We have fitted the X-ray afterglow light curves of 11 GRBs by using the energy injection model of a magnetar with the rotation period in the millisecond order of magnitude. The obtained result shows the validity and universality of the magnetar energy injection model in explaining the shallow decay of afterglows, and simultaneously provides some constraints on the magnetic field strength and rotation period of the central magnetar.