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

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

LP估计在星载GPS运动学定轨中的应用及精度分析

为得到高精度的星载GPS运动学定轨,必须利用观测精度高的相位观测值,但是相位观测值预处理后,仍然存在残余小周跳．在残差服从正态分布情况下LS法是最佳参数解算方法,但该方法不能解决资料的系统误差消除问题,LP估计是处理资料残差分布含有系统误差的有效方法之一．基于LS、LP方法的有效条件和GPS数据预处理的特性,将LP估计方法引入星载GPS运动学定轨数据处理中,以CHAMP卫星资料为例,研究了LP估计在星载GPS运动学定轨中的应用及其精度分析．实践表明:在处理含有残余小周跳的相位观测值时,LP估计比LS更有效,提高了星载GPS运动学定轨精度,但随着残余周跳的进一步修复,LP估计相对于LS估计的优越性越来越弱,在资料完全没有系统误差,残差服从正态分布的情况下,LP估计不能很好地体现其优越性,精度反而低于LS估计．     

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

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

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

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

利用VLBI数据确定"探测一号"卫星的轨道

双星计划的“探测一号”是中国首颗真正严格意义上的科学实验卫星,其运行轨道为中国迄今所发射的卫星中距地球最远,远地点地心距达7．8万公里．采用射电天文的VLBI技术可以对“探测一号”以及更远的深空目标,如探月飞行器实现跟踪．为了验证VLBI技术在我国探月计划中的作用,上海天文台组织了国内目前仅有的上海、乌鲁木齐和昆明3个台站对“探测一号”进行试跟踪,利用对“探测一号”约两天的VLBI观测数据,确定“探测一号”卫星的轨道,对VLBI的定轨能力做初步的探讨．按照测控部门提供的初轨 (其精度仅保证跟踪)推算的轨道与VLBI时延的拟合误差平均约2 km,时延率的拟合误差平均约15 cm／s．而利用VLBI数据定轨后的拟合程度相对于初轨有了很大的改善,结果表明,单独利用VLBI时延定轨,时延的拟合精度约5．5 m,作为外部检核的VLBI时延率的拟合精度在2 cm／s左右．单独利用VLBI时延率定轨,时延率的拟合精度约为1．3 cm／s,作为外部检核的VLBI时延的拟合精度约为29 m．而若将时延和时延率数据联合定轨,采用其内符精度加权,VLBI时延和时延率的残差分别为5．5 m和 2 cm／s．为了合理地评估VLBI定轨的真实精度,利用模拟数据进行误差协方差分析,结果表明VLBI定轨精度受动力学模型误差的影响较大,由于"探测一号”卫星的动力学模型难以精确确定,所以利用两天弧段的VLBI数据确定“探测一号”卫星轨道的位置误差为km量级,而速度误差可达cm／s量级．模拟计算还表明, VLBI和USB数据联合定轨可以大大提高定轨精度．     

星载测高仪交叉点数据定轨研究

从海洋测高卫星的发展状况入手,介绍了用测高仪数据进行定轨的意义.然后,分析了测高仪数据的误差修正模型及交叉点数据计算方法,详细介绍了JASON-1测高数据文件中所涉及的误差修正量和采用的模型.最后,通过仿真数据和实测数据的计算,分析了单用测高仪数据和交叉点数据进行轨道确定所能达到的精度,为未来的实际应用提供了有价值的参考.     

亚洲一号卫星共视法定轨及时刻同步

１９９３年中国科学院同步卫星观测网同步接收亚洲一号卫星的电视信号，进行了多台站时差法定轨及共视法时刻同步试验。试验表明：时差法定轨在经度、纬度方向精度较高，在卫星高度方向误差较大。不进行卫星位置修正时，时刻同步误差小于３μｓ；进行卫星位置修正时的时刻同步误差小于０．９μｓ，可应用于局部或大范围网络时刻同步及同步卫星轨道研究。     

亚洲一号卫星共视法定轨及时刻同步

１９９３年中国科学院同步卫星观测网同步接收亚洲一号卫星的电视信号，进行了多台站时差法定轨及共视法时刻同步试验。试验：时差法定轨在经度、纬度方向精度较高，在卫星高度方向误差较大，不进行卫星位置修正时，时刻同步误差小于３μｓ； 卫星位置修正时的时刻同步误差小于０．９μｓ，可应用于局或大范围网络时刻同步及同步卫星轨道研究。     

初轨的稳健估计算法

讨论单测角资料(数据中包含异常值)初轨的稳健估计算法．该算法既可用于提供可靠的初轨,又可用作资料预处理,剔除异常值．稳健估计算法的崩溃点为16％至25％,可剔除6σ以上的异常值．     

加速仪数据在CHAMP卫星精密定轨中的贡献

通过用几个实例计算，评估了CHAMP卫星加速仪数据对SLR数据定轨精度的贡献．结果显示用加速仪数据代替非保守力对CHAMP卫星进行精密定轨时，其SLR的残差从原来的16．5cm减少到2．7cm，与精密轨道相比，卫星位置误差由18．9cm减少到6．2cm(一天资料)．     

Meridian observations of outer planets at Bordeaux observatory

A systematic bias in right ascensions in the DE200 ephemeris for the outer planets is confirmed.     

Spatial orientations of galaxies in seven Abell clusters of BM type II

As a concluding paper in this series, we present an analysis of the spatial orientations of galaxies in seven Abell clusters of BM type II and compare the results with previous works. Four BM type II clusters (A1767, A1809, A2554, A2721) show a similar preferred alignment: spin vectors of galaxies tend to lie in the local supercluster plane and the projections of the spin vector tend to be oriented perpendicular with respect to the Virgo cluster centre. A preferred alignment of galaxies in cluster A2554 is noticed in both the two- and three-dimensional analyses. In a comparison with previous work, we noticed a systematic change (with distance, radial velocity, morphology and magnitude) in the galaxy alignments from early-type (BM I) to late-type (BM III) clusters. Possible explanations of these systematic changes will be discussed.     

Mode constraints on nonradial pulsations from polarization data

The distortion of the surface of a hot variable star undergoing nonradial pulsation produces a variable polarization. Resultant changes in the Stokes parameters are here calculated analytically in the linear regime as a function of both the pulsation mode (l, m) and the orientation of the stellar rotation axis. It is shown that all of the necessary stellar atmosphere information may be condensed into one simple integral. The analytic form permits a systematic appraisal of the mode-discriminating ability of polarization observations. A suitable systematic approach is outlined and illustrated with reference to the data for BW Vul.     

Critical diameters of detonation of explosives used in spacecraft separation systems

The article presents findings of systematic investigations of the detonation ability of advanced individual explosives, their mixtures, and materials used in spacecraft separation systems. A method for calculation and prediction of their critical detonation diameters is developed.     

Measurement of Kodaikanal white-light images

A program of digitization of the daily white-light solar images from the Kodaikanal station of the Indian Institute of Astrophysics is in progress. A similar set of white-light data from the Mount Wilson Observatory was digitized some years ago. In both cases, areas and positions of individual sunspot umbrae are measured. In this preliminary report, comparisons of these measurements from the two sites are made. It is shown that both area and position measurements are in quite good agreement. The agreement is sufficiently good that it is possible to measure motions and area changes of sunspots from one site to the next, involving time differences from about 12 hours to about 36 hours. This enables us to trace the motions of many more small sunspots than could be done from one site alone. Very small systematic differences in rotation rate between the two sites of about 0.4% are found. A portion of this discrepancy is apparently due to the difference in plate scales between the two sites. Another contributing factor in the difference is the latitude visibility of sunspots. In addition it is suggested that a small, systematic difference in the measured radii at the two sites may contribute a small amount to this discrepancy, but it has not been possible to confirm this hypothesis. It is concluded that in general, when dealing with high precision rotation results of this sort, one must be extremely careful about subtle systematic effects.Operated by the Association of Universities for Research in Astronomy, Inc., under Cooperative Agreement with the National Science Foundation.     

Coronal magnetic field patterns inferred from radio observations

A noise storm center clearly associated to an active center has been followed from January 2 to January 8, 1969. The study of the mean distribution of continuum and of type 1 bursts has shown a global and systematic displacement interpreted in terms of the coronal magnetic structure. This structure is formed by field lines connecting the active center to a stable region of inverse polarity, characterized by the existence of a coronal condensation. These results show the possibility of using this method for the systematic study of coronal magnetic structures, the knowledge of which is vital as far as the problems of particle propagation in the corona and in the interplanetary medium are concerned.     

Anomalous particle-production thresholds through systematic and non-systematic quantum-gravity effects

A growing number of studies is being devoted to the identification of plausible quantum properties of spacetime which might give rise to observably large effects. The literature on this subject is now relatively large, including studies in string theory, loop quantum gravity and non-commutative geometry. It is useful to divide the various proposals into proposals involving a systematic quantum-gravity effect (an effect that would shift the main/average prediction for a given observable quantity) and proposals involving a non-systematic quantum-gravity effect (an effect that would introduce new fundamental uncertanties in some observable quantity). The case of quantum-gravity-induced particle-production-threshold anomalies, a much studied example of potentially observable quantum-gravity effect, is here used as an example to illustrate the differences to be expected between systematic and non-systematic effects.     

Velocity Space of Galactic O-B Stars

Based on the Hipparcos proper motions and available radial velocity data of O-B stars, we have re-examined the local kinematical structure of the young disk population of-1500 O-B stars not including the Gould-belt stars. A systematic warping motion of the stars about the direction to the Galactic center has been reconfirmed. A negative K-term implying a systematic contraction of stars in the solar vicinity has been detected. Two different distance scales are used in order to find out their impact on the kinematical parameters, and we conclude that the adopted distance scale plays an important role in characterizing the kinematical parameters at the present level of the measurement uncertainty.     

Development of an Absolute Accelerometer for Extra-Solar Planet Detection

The method of stellar radial velocity variations has recently shown its capability by the first discovery of several extra-solar planets. Accuracies achieved today are in the range 3-10 m/s. The AAA (absolute astronomical accelerometer) is an instrument which aims to reach the photon noise limit for the measurement of velocity changes, with systematic errors of about 1 m/s, long term. The principle is to use a servo-controlled CCD spectrograph as a null detector, and to register always the lines of the star on the same CCD pixels. Thus, systematic errors linked to the Earth-induced large variations are cancelled. A tunable Fabry-Perot channelled spectrum is also following the star spectrum, while the FP thickness is measured by heterodyne detection of the beats between a tunable laser diode and a stabilized laser diode. A complete prototype of the instrument is operating with laboratory sources and the first results are presented. It is planned to use this system with a new spectrograph, to be coupled to the 152 cm telescope at Observatoire de Haute Provence.     

Systematic Errors of High–Precision Photometric Catalogues

A very accurate imitation of Hipparcos and Tycho Hp, B _T, andV _T magnitudes was made using W, B, V, R magnitudes from the Tien Shan photometric catalogue. The calculated magnitudes were compared to the observed ones. It is shown that there are systematic differences between calculated and observed magnitudes. The systematic errors are supposed to be bound up with the sky scanning procedure on the Hipparcos satellite. Polynomials in powers of coordinates have been proposed to take into account the systematic errors. 6558 stars have been found to be appropriate high-precision photometric standards. This revised version was published online in July 2006 with corrections to the Cover Date.