方差分量估计方法在ERS—2卫星精密定轨中的应用

将方差分量估计（VCE）方法应用于ERS－2卫星的精密定轨,用SLR和PRARE资料计算了1998年前3个月的23个长度为5天的弧段（除了调整轨道的时段外,相邻弧段有两天的重叠）,从观测值残差分析和重叠弧段比较两个方面,考察VCE方法对定轨计算的影响,并给出了各组观测值的平均验后均方差,对观测值残差的分析表明,使用VCE方法明显地改善了观测值的拟合程度,但从阿卑（Abbey)标准对观测值残差的检验结果来看,使用VCE方法不能消除轨道中由力学模型和几何模型误差引入的系统差,重叠弧段比较的结果表明：（1）使用VCE方法缩小了重叠弧段的平均距离差,并改善了一部分权段明显不合理的偏离,使最后得到的轨道具有更均匀的精度,（2）相比较而言,VCE方法使相邻弧段靠拢的趋势在轨道切向体现得较为明显,由各组观测值的平均验后方差可见,说单个标准点观测值而言,部分SLR台站的观测资料在定轨计算中占有比其他观测资料更重的地位,纵观全文,使用VCE后得到的观测值的平衡验后均方差来确定资料的双重将比使用均方差更为合理。     

Application of Accelerometer Data in Precise Orbit Determination of GRACE -A and -B

We investigate how well the GRACE satellite orbits can be determined using the onboard GPS data combined with the accelerometer data.The preprocessing of the accelerometer data and the methods and models used in the orbit determination are presented.In order to assess the orbit accuracy,a number of tests are made,including external orbit comparison,and through Satellite Laser Ranging (SLR) residuals and K-band ranging (KBR) residuals.It is shown that the standard deviations of the position differences between the so-called precise science orbits (PSO) produced by GFZ,and the single-difference (SD) and zero-difference (ZD) dynamic orbits are about 7 cm and 6 cm,respectively.The independent SLR validation indicates that the overall root-mean-squared (RMS) errors of the SD solution for days 309-329 of 2002 are about 4.93cm and 5.22cm,for GRACE-A and B respectively; theoverall RMS errors of the ZD solution are about 4.25 cm and 4.71 cm,respectively.The relative accuracy between the two GRACE satellites is validated by the KBR data to be on a level of 1.29 cm for the SD,and 1.03 cm for the ZD solution.     

TLE-Aided Orbit Determination Using Single-Station SLR Data

It is difficult to use the single-station satellite laser ranging (SLR) data for orbit determination, due to the singular geometrical distribution of the observations. The single-station data produced by performing the diffuse- reflection SLR on the earth-orbiting space debris are therefore ineffective for orbit improvement. To solve this problem, we propose an orbit determination method by using single-station SLR data in aid of the two-line element set (TLE). For verifying its feasibility, this method is implemented and applied to the orbit determination of the satellite Ajisai, using the single-station SLR data of five passes in one day and the corresponding TLE. And on this basis, the five-day orbit prediction is generated, the result indicates that the errors of predicted positions are less than 40 m. In addition, the potential application of this method in the orbit improvement of space debris is discussed.     

Determination of the orbit of a natural satellite with hyperbolic flybys

We show that, when a natural satellite like Titan is invisible (e.g., due to an opaque atmosphere) its planetary orbit and its mass can be determined by tracking a spacecraft in close flybys. This is an important problem in the Cassini mission to the Saturnian system, which will be greatly improved by a good astrometric model for all its main components; in particular, an accuracy of a few hundred meters for the orbit of Titan is necessary to allow a measurement of its moment of inertia. The orbit of the spacecraft is the union of elliptical arcs, joined by short hyperbolic transitions: a problem of singular perturbation theory, whose solution leads to a matching condition between the inner hyperbolic orbit and the elliptical orbital elements. Since the inner elements are given in terms of the relative position and velocity of the spacecraft, accurate Doppler measurements in both regions can provide a satisfactory determination of Titan's position and velocity, hence of its Keplerian elements. The errors in this determination are discussed on the basis of the expected Allan deviation of the Doppler method; it is found that the driving errors are those in the elliptical arcs; the fractional errors in Titan's orbital elements are expected to be 10^–7. It is also possible to measure the mass of the satellite; however, when the eccentricity e of the flybys is large, the mass and a scaling transformation are highly correlated and the fractional error in the mass is expected to be e times worse.     

LAGEOS卫星精密定轨及残差分析

介绍LAGEOS卫星用卫星激光测距(SLR)资料精密定轨和在精密定轨基础上的残差分析处理．SLR资料的分析处理方法、采用的力学模型和解参数的多少根据目标的不同而不同．对不同的方案进行了详细分析、比较．其关键是对LAGEOS卫星进行精密定轨．得出了目前上海天文台在SLR资料的分析处理中已采用的解算模式．作为例证，该解算模式分析处理了1998年12月31日至1999年6月29日LAGEOS卫星的SLR资料，得到每3天一个弧段的精密轨道．结果显示每3天一个弧段的定轨残差rms都小于2厘米．上海天文台从1999年10月1日开始了全球LAGEOS—l和LAGEOS—2资料的快速分析服务，结果可从APSG的网址：http：／／center．shao．ac．cn／APSG／result获得．     

基于双频星载GPS数据的LEO卫星运动学定轨研究

运动学定轨是星载GPS特有的定轨方法,该方法不依赖于任何力学模型(地球重力场、大气阻力及太阳辐射压等),尤其适用于受大气阻力影响严重的低轨卫星定轨.基于双频星载GPS数据,研究了运动学定轨原理,讨论了数据预处理方法,建立了一套非差运动学定轨算法.并以GRACE (Gravity Recovery And Climate Experiment)-A、B卫星2008年2月实测数据作为试算验证了本研究方法的有效性和可靠性.GRACE 卫星实测数据计算结果表明:运动学定轨能达到5 cm精度(相对于SLR (Satellite Laser Ranging)),与动力学和简化动力学定轨精度相当.     

两种观测技术综合精密定轨的探讨

利用T／P卫星的SLR和DORIS实测资料,对两种观测技术综合精密定轨中的加权及其对定轨影响的问题作了初步的探讨。根据所提出的一种经验性的加权方法进行了综合定轨计算,结果表明：对于两种不同技术的观测,相对权选取的恰当与否将影响综合定轨的精度;综合定轨的最优加权不仅依赖于观测资料的精度,还与观测资料的多少和几何分布有关;通过选用最优加权,可使得综合轨的精度优于仅用其中一种技术的定轨精度,综合定轨能有效地提高定轨精度。     

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

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

Calculation of the intermediate perturbed orbit from three or more positions of a small body on the celestial sphere

Two new methods are described for finding the orbit of a small celestial body from three or more pairs of angular measurements and the corresponding time points. The methods are based on, first, the approach that has been developed previously by the author to the determination, from a minimum number of observations, of intermediate orbit considering most of the perturbations in the bodies’ motion and, second, Herget’s algorithmic procedure enabling the introduction of additional observations. The errors of orbital parameters calculated by the proposed methods are two orders of magnitude smaller than the corresponding errors of the traditional approach based on the construction of an unperturbed Keplerian orbit. The thus-calculated orbits of the minor planets 1566 Icarus, 2002 EC1, and 2010 TO48 are used to compare the results of Herget’s multiposition procedure and the new methods. The comparison shows that the new methods are highly effective in the study of perturbed motion. They are particularly beneficial if high-precision observational data covering short orbital arcs are available.     

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影响.     

On the Laplacian orbit determination of asteroids

A modified Laplacian technique is described for initial orbit determination of asteroids from CCD observations and its applications for orbit determination of the main belt asteroids and near Earth asteroids. The proposed modification is based on a simultaneous improvement of both the orbital elements and the derivatives of spherical coordinates in frames of Laplace's method. It provides an orbit which represents the used observations with the residuals comparable with errors of these observations. The improved values of the derivatives might be used as ephemeris parameters for identification of newly discovered objects.     

两行根数辅助的SLR单站定轨

单站测距资料定轨的困难限制了漫反射SLR(Satellite Laser Ranging)测距资料的应用.为此,提出利用两行根数模拟多站SLR测距资料作为辅助,实现单站SLR测距资料定轨的方法.该方法对卫星Ajisai单站SLR测距资料定轨并生成5 d预报轨道,误差小于40 m,实现利用单站测距资料的轨道改进,验证了方法的可行性.     

GPS卫星的激光测距和应用研究

简述了SLR和GPS跟踪技术的发展概况。详细介绍了近几年来对GPS－35、36卫星的激光测距进展和应用研究的情况。给出了残差分析的初步结果,并指出了目前GPS卫星的稀少的激光测距资料对卫星精密定轨和站坐标的解算是有价值的。同时,也简单地评述了GPS卫星的微波和激光跟踪技术各自的优势和弱点．建议联合利用GPS-35、36卫星的SLR和GPS观测资料来开展有关的应用研究。最后,对今后的应用研究工作提出了建议．     

The basic parameters of static pointing model of telescope

In this paper, models of the static pointing errors and their correction of an altazimuth telescope are discussed, and a stable model with a small number of basic parameters is given. Compared with the common spherical harmonics model, the present model is not sensitive to the distribution of observing data and has a much higher correction accuracy. The residuals from this model follow the normal distribution in statistics. It can meet the need of precise orbit determination.     

Error analyses of resonant orbits for geodesy

An error analysis of resonant orbits for geodesy indicates that attempts to use resonance to recover high order geopotential coefficients may be seriously hampered by errors in the geopotential. This effect, plus the very high correlations (up to .99) of the resonant coefficients with each other and the orbital period in single satellite solutions, makesindividual resonant orbits of limited value for geodesy. Multiple-satellite, single-plane solutions are only a slight improvement over the single satellite case. Accurate determination of high order coefficients from low altitude resonant satellites requires multiple orbit planes and small drift-periods to reduce correlations and effects of errors of non-resonant geopotential terms. Also, the effects of gravity model errors on low-altitude resonant satellites make the use of tracking arcs exceeding two to three weeks of doubtful validity. Because high-altitude resonant orbits are less affected by non-resonant terms in the geopotential, much longer tracking arcs can be used for them.     

A new method for the determination of an intermediate orbit using three positions of a small body on the celestial sphere

We propose a new method for the determination of the preliminary orbit of a small celestial body using three pairs of its angular coordinates in three moments of time. The method is based on the use of the intermediate orbit we constructed earlier using three position vectors and the corresponding time moments. This intermediate orbit accounts for the main part of the perturbations of the motion of the body under study. We compare the results obtained by the classical Lagrange-Gauss method, Herrick-Gibbs method, generalized Herrick-Gibbs method, and the new method by the examples of the determination of the orbit of the small planet 1566 Icarus. The comparison showed that the new method is a highly efficient tool for the study of perturbed motion. It is especially efficient when applied to high-precision observational data covering short arcs of the orbit.     

A New Solution Assessment Approach and Its Application to Space Geodesy Data Analysis

The statistics of the residuals are used in this paper to perform a quality assessment of the solutions from space geodesy data analysis. With the stochastic estimation and the relatively arbitrary empirical parameters being employed to absorb unmodelled errors, it has long been noticed that different estimate combinations or analysis strategies may achieve the same level of fitting yet result in significantly different solutions. Based on the postulate that no conceivable signals should remain in the residuals, solutions of the same level of root mean square error (RMS) and variance–covariance may be differentiated in the sense that for reasonable solutions, the residuals are virtually identical with noise. While it is possible to develop complex noise models, the Gaussian white noise model simplifies the solution interpretation and implies the unmodelled errors have been smoothed out. Statistical moments of the residuals as well as the Pearson chi-square are computed in this paper to measure the discrepancies between the residuals and Gaussian white noise. Applying to both satellite laser ranging (SLR) and global positioning system (GPS) data analysis, we evaluate different parameter estimate combinations and/or different strategies that would be hardly discriminated by the level of fitting. Unlike most solution assessment methods broadly termed as external comparison, no information independent of the data analyzed is required. This makes the immediate solution assessment possible and easy to carry out. While the external comparison is the best and most convincing quality assessment of the solution, the statistics of the residuals provide important information on the solutions and, in some cases as discussed in this paper, can be supported with external comparison.     

Satellite laser ranging and its applications

Satellite laser ranging (SLR) provides an important capability for precise orbit determination and for geophysical parameter estimation to support a number of contemporary geodynamic and oceanographic investigations. The precision of the SLR measurement has improved from the early meter-level systems to the current capabilities of a few centimeters for the best systems. The accuracy of the orbits and geophysical parameter recovery have shown an associated improvement. Polar motion with accuracies of 2 mas, station coordinates better than 10 cm, and interstation baseline rates indicative of tectonic motion are determined routinely with the current set of global SLR data. This discussion reviews the SLR measurement, analysis approach, and some of the recent results derived from the current SLR data set.     

利用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数据联合定轨可以大大提高定轨精度．