-
卫星激光测距(satellite laser ranging,SLR)作为空间大地测量的重要技术之一,在卫星轨道的确定与评估以及地球参考架的建立与维持等方面起着重要作用[1-14],如在国际地球参考架ITRF2014(International Terrestrial Reference Frame 2014)中,SLR是定义参考架原点的唯一手段,并与甚长基线干涉测量(very long baseline interferometry,VLBI)技术共同定义了参考架的尺度因子[10]。目前参加国际SLR网联测的中国SLR站有上海站、北京站、武汉站、长春站和昆明站5个固定站[15-16],为弥补固定站的不足,我国自主研制了流动卫星激光测距[17-18]。2016-05-23至2016-06-16,直径为35 cm的流动SLR测距仪在西安进行了试验性观测,本文采用全球Lageos-1卫星数据计算了西安流动SLR站坐标。
在利用全球SLR数据解算站坐标时,不仅需要可靠的函数模型(包括动力学模型、观测值模型),同样也需要正确的随机模型[19],以便根据各台站数据的实际观测精度,合理确定观测资料的贡献,提高坐标解的精度及可靠性。本文从权值确定和粗差编辑两方面,研究不同随机模型对西安流动SLR站坐标解算的影响,并从坐标解算精度、坐标解算稳定性、残差加权均方差以及观测资料利用率4个方面进行了分析。
-
假设k组相互独立的观测值Li(i=1,2…k)的权阵分别为Pi,其误差方程为:
$$ V = AX - L $$ (1) 式中,
$$ \begin{array}{l} V = \left[ {\begin{array}{*{20}{c}} {{V_1}}\\ {{V_2}}\\ \vdots \\ {{V_k}} \end{array}} \right];A = \left[ {\begin{array}{*{20}{c}} {{A_1}}\\ {{A_2}}\\ \vdots \\ {{A_k}} \end{array}{\rm{}}\begin{array}{*{20}{c}} B\\ 0\\ \vdots \\ 0 \end{array}} \right];X = \left[ {\begin{array}{*{20}{c}} {{{\bar X}_1} - {{\bar X}_1}^{\rm{*}}}\\ {{{\bar X}_2} - {{\bar X}_2}^{\rm{*}}} \end{array}} \right];\\ L = \left[ {\begin{array}{*{20}{c}} {{L_{1O}} - {L_{1C}}}\\ {{L_{2O}} - {L_{2C}}}\\ \vdots \\ {{L_{kO}} - {L_{kC}}} \end{array}} \right];P = \left[ {\begin{array}{*{20}{c}} {{P_1}}&{}&{}&0\\ 0&{{P_2}}&{}&0\\ \vdots &{}& \ddots & \vdots \\ 0&{}&{}&{{P_k}} \end{array}} \right]。 \end{array} $$ 其中,下标O、C分别代表观测值和计算值;X 1、X 1*分别为卫星相关参数(包括卫星的位置、速度及动力学参数)的真值和参考解;X 2、X 2*分别为西安站参数(包括站坐标及测站距离偏差)的真值和参考解;Ai为观测值对卫星相关参数的偏导数矩阵;B为西安SLR观测值对西安站参数的偏导数矩阵;L1O为西安站SLR观测值。式(1)的解为:
$$ \hat X = {({A^{\rm{T}}}PA)^{ - 1}}{A^{\rm{T}}}PL $$ (2) $$ {\mathit{\pmb{\text{ }\!\!\Sigma\!\!\text{ }}}_{X}}={{({{A}^{\text{T}}}PA)}^{-1}}\hat{\sigma }_{0}^{2} $$ (3) $$ \hat \sigma _0^2 = {V^{\rm{T}}}PV/\left( {m - n} \right) $$ (4) 式中,m为观测个数;n为待估参数个数; $ {\Sigma _\mathit{X}}$为X的协方差矩阵;$ \hat \sigma _0^2$为单位权方差。
借鉴抗差估计将粗差归入随机模型处理的思想[20-22],观测值权值取值如下:
$$ p_{ij}^{} = \left\{ {\begin{array}{*{20}{c}} {{p_i}, {\rm{}}\left| {{u_{ij}}} \right| \le {k_0}}\\ {0, {\rm{}}\left| {{u_{ij}}} \right| > {k_0}} \end{array}} \right. $$ (5) 1)当pi=1,$ {u_{ij}} = {v_{ij}}/{\hat \sigma _0}$,k0=3时,vij为第i台站的第j个观测值残差,${\hat \sigma _0} $为所有台站观测资料的均方差。此时解为等权最小二乘解,粗差编辑标准为所有站残差3倍中误差。
2)当${p_i} = 1/\hat \sigma _i^2 $,$ {u_{ij}} = {v_{ij}}/{\hat \sigma _0}$,k0=3时,vij和$ {\hat \sigma _0}$意义同上,$ {\hat \sigma _i}$为第i台站观测资料的均方差。此时解为基于方差分量估计的最小二乘解,粗差编辑标准为所有站残差3倍中误差。
3)当$ {p_i} = 1/\hat \sigma _i^2$,${u_{ij}} = {v_{ij}}/\hat \sigma _i^2 $,k0=3时,vij、${\hat \sigma _0} $和$ {\hat \sigma _i}$意义同上。此时解为基于方差分量估计的最小二乘解,粗差编辑标准为各台站残差3倍中误差。
观测值权值采用抗差估计等价权函数时[22]:
$$ p_{ij}^{} = \left\{ {\begin{array}{*{20}{l}} {{p_i}, \left| {{u_{ij}}} \right| \le {k_0}}\\ {\frac{{{p_i}{k_0}}}{{\left| {{u_{ij}}} \right|}}{{\left( {\frac{{{k_1} - \left| {{u_{ij}}} \right|}}{{{k_1} - {k_0}}}} \right)}^2}, {k_0} < \left| {{u_{ij}}} \right| \le {k_1}}\\ {0, \left| {{u_{ij}}} \right| > {k_1}} \end{array}} \right. $$ (6) 当$ {p_i} = 1/\hat \sigma _i^2$,${u_{ij}} = {v_{ij}}/\hat \sigma _i^2 $,k1 > k0时,${v_{ij}} $和$ {\hat \sigma _i}$意义同上,解为基于抗差方差分量估计的最小二乘解。
-
2016-05-23至2016-06-16,西安流动SLR站共观测到Lageos-1卫星5 d 36个标准点数据,其中5月23日标准点10个,5月24日标准点8个,6月3日标准点5个,6月15日标准点9个,6月16日标准点4个。标准点即选取固定时间段(Lageos-1卫星为2 min),按一定技术把该时间段内所有观测数据压缩成的一个数据点。
西安流动SLR站坐标是与Lageos-1卫星轨道一起求解的,采用全球Lageos-1观测资料,来自美国国家宇航局(National Aeronautics and Space Administration,NASA),数据格式为CRD(consolidated laser ranging data)。这种数据格式可以满足高重复率激光测距和远距离行星际激光测距的数据要求。为了提高参数求解精度,本文考虑了尽可能完善的力学模型和测量模型, SLR站坐标采用卫星激光测距框架2014台站坐标序列。解算参数为初始历元的6个轨道根数,即类阻力系数、太阳光压系数、T和N方向的经验加速度、西安站坐标及距离偏差。其中类阻力系数和太阳光压系数3 d估计一组,T和N方向的经验加速度每天估计一组,其他参数30 d估计一组。
本文计算站坐标采用的3个弧段分别为:①2016-05-01至2016-05-30的30 d观测数据作为弧段1,全球SLR测站总观测数为6 422,其中西安流动站观测数据个数为18;②2016-05-31至2016-06-29的30 d观测数据作为弧段2,全球SLR测站总观测数为5 473,其中西安流动站观测数据个数为18;③2016-05-20至2016-06-18的30 d观测数据作为弧段3,全球SLR测站总观测数为5 740,其中西安流动站观测数据个数为36。根据采用的权值及粗差编辑标准不同,本文采用下述4种方案计算。
方案1 等权最小二乘,粗差编辑标准为所有站残差3倍中误差。
方案2 基于方差分量估计的最小二乘,粗差编辑标准为所有站残差3倍中误差。
方案3 基于方差分量估计的最小二乘,粗差编辑标准为各台站残差3倍中误差。
方案4 基于抗差方差分量估计的最小二乘解,k0和k1分别取1.5和6.0。
-
按照上述4种方案分别计算了西安流动SLR站的三维坐标,表 1给出了4种计算方案采用的有效观测个数及残差加权均方差(weighted root mean square,WRMS)。
从表 1可以看出, 4种方案的残差WRMS都小于3 cm,方案1采用等权最小二乘计算,没有考虑到各SLR台站数据质量相差较大的实际情况,在4种方案中轨道拟合的精度最低,残差WRMS最大。方案2采用方差分量估计,虽然在轨道计算过程中考虑了各台站数据质量的不同,提高了轨道拟合精度,但是粗差剔除没有考虑各台站数据质量的不同,因此删除的观测资料个数最多,残差WRMS较方案1小。方案3在轨道计算和粗差剔除的过程中都考虑了各台站数据质量的不同,既提高了轨道拟合精度,又增加了采用的观测资料个数,残差WRMS较方案1、2都小。方案4采用的观测资料最多,残差WRMS最小。这是由于方案4采用抗差方差分量估计,不仅赋予了各个台站的观测资料合理的权值,而且合理有效地利用了观测资料,削弱了观测资料中粗差的影响。
表 1 SLR 3个弧段4种方案数据处理有效观测数、剔除观测数及残差WRMS情况比对
Table 1. Comparison of the Number of SLR Data Used and Rejected and WRMS Between Three Arcs and Four Projects
弧段 观测总数 方案 有效观测数 剔除观测数 WRMS/cm 1 6 422 1 5 717 705 2.70 2 5 675 747 1.75 3 6 031 391 1.61 4 6 084 338 1.11 2 5 473 1 4 958 515 2.89 2 4 881 592 2.03 3 5 213 260 1.79 4 5 273 200 1.43 3 5 740 1 5 137 603 2.77 2 5 061 679 1.98 3 5 453 287 1.94 4 5 518 222 1.55 图 1给出了西安SLR流动站观测资料的验后残差。从图 1可以看出,西安SLR流动站的观测资料达到了厘米级。
表 2给出了4种方案计算的西安流动SLR站坐标值及其均方差,均方差mX、mY、mZ是X、Y、Z的内符合精度,并不能完全反映站坐标的解算精度。为客观反映站坐标精度和稳定性,本文计算了每种方案的站坐标平均值及其均方差,结果如表 3所示。
表 2 西安流动SLR站坐标及其均方差
Table 2. Coordinates and Root Mean Squares of Xi'an Mobile SLR Station
弧段 方案 X/m Y /m Z/m mX/cm mY/cm mZ/cm 1 1 -1 718 607.239 0 4 995 239.455 5 3 563 076.560 5 2.76 5.93 6.85 2 -1 718 607.242 7 4 995 239.469 0 3 563 076.585 4 0.93 1.92 2.24 3 -1 718 607.257 0 4 995 239.477 9 3 563 076.600 5 0.90 1.83 2.14 4 -1 718 607.271 6 4 995 239.507 4 3 563 076.629 4 0.68 1.48 1.71 2 1 -1 718 607.288 8 4 995 239.472 7 3 563 076.630 6 2.95 6.02 5.05 2 -1 718 607.295 4 4 995 239.502 9 3 563 076.656 5 1.13 2.26 1.94 3 -1 718 607.296 0 4 995 239.494 2 3 563 076.655 9 1.16 2.33 2.00 4 -1 718 607.305 5 4 995 239.486 2 3 563 076.666 0 0.77 1.51 1.33 3 1 -1 718 607.278 0 4 995 239.489 0 3 563 076.651 0 1.81 3.62 3.17 2 -1 718 607.285 4 4 995 239.505 2 3 563 076.664 5 0.88 1.76 1.55 3 -1 718 607.285 1 4 995 239.504 2 3 563 076.663 9 0.88 1.75 1.54 4 -1 718 607.290 6 4 995 239.507 6 3 563 076.668 1 0.68 1.33 1.18 表 3 西安流动SLR站坐标平均值及其均方差
Table 3. The Mean Value and the Root Mean Square of Xi'an Mobile SLR Station Coordinate
方案 X/m Y /m Z/m mX/cm mY/cm mZ/cm 1 -1 718 607.268 6 4 995 239.472 4 3 563 076.614 0 2.62 1.68 4.75 2 -1 718 607.274 5 4 995 239.492 4 3 563 076.635 5 2.80 2.03 4.35 3 -1 718 607.279 4 4 995 239.492 1 3 563 076.640 1 2.01 1.33 3.45 4 -1 718 607.289 2 4 995 239.500 4 3 563 076.654 5 1.70 1.23 2.18 从表 3可以看出,4种方案计算的站坐标值比较接近,不同方案和弧段之间X、Y、Z轴坐标的最大差异分别为6.7 cm、5.2 cm、10.8 cm。对于同一计算方案而言,弧段3的坐标估计精度是最高的,这是由于弧段3采用的西安站SLR数据最多。对于同一弧段而言,方案4计算的坐标均方差精度都是最高的,这是由于方案4采用抗差方差分量估计,能够合理利用各台站的观测资料。
从表 3可以看出,在4种方案中,方案4计算的坐标重复性最好,坐标解算的稳定性最强, 这是由于方案4可以合理利用各台站的观测资料,使参数解算更符合观测数据的实际观测情况。不同方案、不同弧段的计算结果表明,西安流动SLR站的坐标解算精度可以达到厘米级。
由分析可得出结论,本文通过优化所有观测数据的随机模型,提高了轨道拟合精度,改善了西安流动站坐标的解算精度和稳定性。
-
卫星激光测距是空间大地测量的重要技术之一,在卫星轨道的确定与评估以及地球参考架的建立与维持等方面[23-24]起着重要作用。在利用全球SLR数据解算站坐标时,正确的随机模型能够合理确定观测资料的贡献,提高坐标解的精度。本文从权函数和粗差编辑两方面研究了不同随机模型对西安SLR站坐标解算的影响,采用全球Lageos-1卫星观测数据,利用4种方案计算了西安流动SLR站坐标,并对各方案结果的精度进行了评价。结果表明,西安流动SLR站的观测精度和坐标解算精度均达到厘米级。同时,随机模型直接影响SLR站坐标的解算结果及可靠性;对于相同的计算弧段,抗差方差分量估计得到的站坐标精度最高、结果最稳定,残差加权均方差最小,观测资料利用率也最高;对于相同的计算方案,采用的SLR数据越多,坐标估计精度越高。下一步,将在利用SLR资料计算极移的研究中分析随机模型的影响。
Influence of Stochastic Model on the Coordinate of Xi'an Mobile SLR Station Determination
-
摘要: 合理的参数估计及精度评定不仅需要可靠的函数模型,而且需要正确的随机模型。从权函数和粗差编辑两方面,研究了不同随机模型对西安流动卫星激光测距(satellite laser ranging,SLR)站坐标解算的影响,采用全球Lageos-1卫星观测数据计算了西安流动SLR站坐标。计算结果表明:①西安流动SLR站的观测精度和坐标解算精度均达到厘米级。②随机模型直接影响SLR站坐标的解算结果及可靠性;对于相同的计算弧段,抗差方差分量估计得到的站坐标精度最高、结果最稳定,残差加权均方差最小,观测资料利用率也最高;对于相同的计算方案,采用的SLR数据越多,坐标估计精度越高。
-
关键词:
- 卫星激光测距 /
- 随机模型 /
- 抗差方差分量估计 /
- Lageos-1卫星 /
- 流动SLR站
Abstract: Reasonable parameter estimation and precision evaluation not only need reliable function model but also need right stochastic model. The influence of stochastic model on the coordinate of Xi'an mobile satellite laser ranging (SLR) station determination is researched in the ways of weight function and outlying observations edit. The coordinate of Xi'an mobile SLR station is determined by using the data of Lageos-1 of all SLR observation stations in the world, which observed from May 1 to June 29, 2016. When the coordinates of Xi'an mobile SLR station are computed, SLRF2014 station positions and velocities are employed. The coordinates of Xi'an mobile SLR station are computed with the orbit of satellite Lageos-1 determination. Three thirty-day orbital arcs and four computation plans are used to compute the coordinates of Xi'an mobile SLR station. There are sufficient numbers of SLR observations to compute reliable orbits, and the values of post-fit range residual root mean square (RMS) is about 2 cm. The difference of the results of the three thirty-day orbital arcs and the four computation plans are compared. By comparison and analysis, we find that the precision of Xi'an mobile SLR observation and the precision of the coordinate computed both reached the level of centimeter. Furthermore, the stochastic model has direct effects on the coordinate of SLR station determination result and reliability. The plan of robust variance component estimation can provide the most precision coordinate, the most stable result, the smallest weighted root mean square (WRMS) and the largest number of SLR data used, when using same computation arc. The more SLR data is used, the higher precision of coordinate is gotten, when using same computation plan. -
表 1 SLR 3个弧段4种方案数据处理有效观测数、剔除观测数及残差WRMS情况比对
Table 1. Comparison of the Number of SLR Data Used and Rejected and WRMS Between Three Arcs and Four Projects
弧段 观测总数 方案 有效观测数 剔除观测数 WRMS/cm 1 6 422 1 5 717 705 2.70 2 5 675 747 1.75 3 6 031 391 1.61 4 6 084 338 1.11 2 5 473 1 4 958 515 2.89 2 4 881 592 2.03 3 5 213 260 1.79 4 5 273 200 1.43 3 5 740 1 5 137 603 2.77 2 5 061 679 1.98 3 5 453 287 1.94 4 5 518 222 1.55 表 2 西安流动SLR站坐标及其均方差
Table 2. Coordinates and Root Mean Squares of Xi'an Mobile SLR Station
弧段 方案 X/m Y /m Z/m mX/cm mY/cm mZ/cm 1 1 -1 718 607.239 0 4 995 239.455 5 3 563 076.560 5 2.76 5.93 6.85 2 -1 718 607.242 7 4 995 239.469 0 3 563 076.585 4 0.93 1.92 2.24 3 -1 718 607.257 0 4 995 239.477 9 3 563 076.600 5 0.90 1.83 2.14 4 -1 718 607.271 6 4 995 239.507 4 3 563 076.629 4 0.68 1.48 1.71 2 1 -1 718 607.288 8 4 995 239.472 7 3 563 076.630 6 2.95 6.02 5.05 2 -1 718 607.295 4 4 995 239.502 9 3 563 076.656 5 1.13 2.26 1.94 3 -1 718 607.296 0 4 995 239.494 2 3 563 076.655 9 1.16 2.33 2.00 4 -1 718 607.305 5 4 995 239.486 2 3 563 076.666 0 0.77 1.51 1.33 3 1 -1 718 607.278 0 4 995 239.489 0 3 563 076.651 0 1.81 3.62 3.17 2 -1 718 607.285 4 4 995 239.505 2 3 563 076.664 5 0.88 1.76 1.55 3 -1 718 607.285 1 4 995 239.504 2 3 563 076.663 9 0.88 1.75 1.54 4 -1 718 607.290 6 4 995 239.507 6 3 563 076.668 1 0.68 1.33 1.18 表 3 西安流动SLR站坐标平均值及其均方差
Table 3. The Mean Value and the Root Mean Square of Xi'an Mobile SLR Station Coordinate
方案 X/m Y /m Z/m mX/cm mY/cm mZ/cm 1 -1 718 607.268 6 4 995 239.472 4 3 563 076.614 0 2.62 1.68 4.75 2 -1 718 607.274 5 4 995 239.492 4 3 563 076.635 5 2.80 2.03 4.35 3 -1 718 607.279 4 4 995 239.492 1 3 563 076.640 1 2.01 1.33 3.45 4 -1 718 607.289 2 4 995 239.500 4 3 563 076.654 5 1.70 1.23 2.18 -
[1] Pearlman M R, Degnan J J, Bosworth J M. The International Laser Ranging Service[J]. Advances in Space Research, 2002, 30(2): 135-143 doi: 10.1016/S0273-1177(02)00277-6 [2] 黄诚.利用Lageos激光测距资料精密测定地球自转参数[D].上海: 中国科学院上海天文台, 1985 Huang Cheng. Accurately Determining the Earth Orientation Parameters with the Data of Lageos Satellites Laser Ranging[D]. Shanghai: Shanghai Astronomical Observatory, Chinese Academy of Science, 1985 [3] 赵罡, 王小亚, 吴斌.不同测站卫星质心不同改正对卫星激光测距定轨精度的影响分析[J].测绘学报, 2012, 41(2):165-170 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=chxb201202003 Zhao Gang, Wang Xiaoya, Wu Bin. Effect Analysis of System Dependent Center-of-Mass Correction on Precision of SLR Orbit Determination[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(2):165-170 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=chxb201202003 [4] 孔垚, 孙保琪, 杨旭海, 等.利用SLR数据进行北斗卫星广播星历精度分析[J].武汉大学学报·信息科学版, 2017, 42(6):831-837 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201706019 Kong Yao, Sun Baoqi, Yang Xuhai, et al. Precision Analysis of BeiDou Broadcast Ephemeris by Using SLR Data[J]. Geomatics and Information Science of Wuhan University, 2017, 42(6):831-837 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201706019 [5] 周旭华, 王晓慧, 赵罡, 等. HY2A卫星的GPS/DORIS/SLR数据精密定轨[J].武汉大学学报·信息科学版, 2015, 40(8):1 000-1 005, 1 011 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whchkjdxxb201508002 Zhou Xuhua, Wang Xiaohui, Zhao Gang, et al. The Precise Orbit Determination for HY2A Satellite Using GPS, DORIS and SLR Data[J]. Geomatics and Information Science of Wuhan University, 2015, 40(8):1 000-1 005, 1 011 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whchkjdxxb201508002 [6] 秦显平, 焦文海, 程芦颖, 等.利用SLR检核CHAMP卫星轨道[J].武汉大学学报·信息科学版, 2005, 30(1):38-41 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200501010 Qin Xianping, Jiao Wenhai, Cheng Luying, et al. Evaluation of CHAMP Satellite Orbit with SLR Measurements[J]. Geomatics and Information Science of Wuhan University, 2005, 30(1):38-41 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200501010 [7] 彭汉兵, 杨元喜, 王刚, 等.星蚀期北斗卫星轨道性能分析——SLR检核结果[J].测绘学报, 2016, 41(6):639-645 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=chxb201606002 Peng Hanbing, Yang Yuanxi, Wang Gang, et al. Performance Analysis of BDS Satellite Orbits During Eclipse Periods: Results of Satellite Laser Ranging Validation[J]. Acta Geodaetica et Cartographica Sinica, 2016, 41(6):639-645 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=chxb201606002 [8] 魏二虎, 万丽华, 金双根, 等.联合GNSS和SLR观测对地球自转参数的解算与分析[J].武汉大学学报·信息科学版, 2014, 39(5):581-585 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201405015 Wei Erhu, Wan Lihua, Jin Shuanggen, et al. Estimation of ERP with Combined Observations of GNSS and SLR[J]. Geomatics and Information Science of Wuhan University, 2014, 39(5):581-585 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201405015 [9] 秦显平, 杨元喜, 王刚, 等.SLR资料精密测定GLONASS卫星轨道[J].武汉大学学报·信息科学版, 2003, 28(4):440-443 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200304016 Qin Xianping, Yang Yuanxi, Wang Gang, et al. GLONASS Orbit Determination by Using SLR Data[J]. Geomatics and Information Science of Wuhan University, 2003, 28(4):440-443 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200304016 [10] Altamimi Z, Rebischung P, Métivier L. Analysis and Results of ITRF2014 (IERS Technical Note No. 38)[M]. Westerstede, Germany: MegaDruck.de Produktions-und Vertriebs GmbH, 2017 [11] 耿涛, 赵齐乐.基于全球激光观测资料的COMPASS-M1轨道质量评定[J].武汉大学学报·信息科学版, 2009, 34(11): 1 290-1 292 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whchkjdxxb200911007 Geng Tao, Zhao Qile. Determining Orbit of COMPASS-M1 Using International Laser Ranging Data[J]. Geomatics and Information Science of Wuhan University, 2009, 34(11): 1 290-1 292 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whchkjdxxb200911007 [12] 李建文, 陈刘成.单颗导航卫星精密定轨与激光精度评估[J].武汉大学学报·信息科学版, 2009, 34(11): 1 275-1 279 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whchkjdxxb200911004 Li Jianwen, Chen Liucheng. Single Navigation Satellite Orbit Determination and Evaluation by SLR[J]. Geomatics and Information Science of Wuhan University, 2009, 34(11):1 275-1 279 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=whchkjdxxb200911004 [13] 路晓峰, 贾小林, 杨志强.利用SLR和GPS双频相位平滑伪距资料测定导航卫星钟差[J].武汉大学学报·信息科学版, 2008, 33(3):237-240 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200803005 Lu Xiaofeng, Jia Xiaolin, Yang Zhiqiang. Determination of Navigation Satellite Clock Bias Using SLR and GPS Dual Frequency Phase-Smoothed Pseudo-Range Data[J]. Geomatics and Information Science of Wuhan University, 2008, 33(3):237-240 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200803005 [14] 毛悦, 贾小林, 秦显平.利用SLR及伪距数据进行卫星预报钟差精度验证[J].武汉大学学报·信息科学版, 2007, 32(12):1 160-1 163 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200712018 Mao Yue, Jia Xiaolin, Qin Xianping. Validating Clock Bias of Satellite Using SLR and Pseudorange Data[J]. Geomatics and Information Science of Wuhan University, 2007, 32(12):1 160-1 163 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200712018 [15] 赵春梅, 瞿锋, 程鹏飞, 等.阿根廷圣胡安激光测距系统的SLR数据质量分析[J].测绘学报, 2008, 37(3):226-232 http://d.old.wanfangdata.com.cn/Periodical/chxb200803012 Zhao Chunmei, Qu Feng, Cheng Pengfei, et al.Data Quality Analysis of Argentina San Juan Laser Ranging System[J]. Acta Geodaetica et Cartographica Sinica, 2008, 37(3):226-232 http://d.old.wanfangdata.com.cn/Periodical/chxb200803012 [16] 朱元兰, 冯初刚, 张飞鹏.用中国卫星激光测距资料解算地球定向参数[J].天文学报, 2006, 47(4):441-449 doi: 10.3321/j.issn:0001-5245.2006.04.011 Zhu Yuanlan, Feng Chugang, Zhang Feipeng. Earth Orientation Parameter Solved by Lageos Chinese SLR Data[J]. Acta Astronomica Sinica, 2006, 47(4):441-449 doi: 10.3321/j.issn:0001-5245.2006.04.011 [17] 冯初刚, 蔡五三, 朱元兰, 等.中国卫星激光测距流动站地心坐标的归算[J].科学通报, 2003, 48(13): 1 402-1 404 http://d.old.wanfangdata.com.cn/Periodical/kxtb200313008 Feng Chugang, Cai Wusan, Zhu Yuanlan, et al. The Coordinate of Chinese Mobile SLR Station Determination[J]. Chinese Science Bulletin, 2003, 48(13): 1 402-1 404 http://d.old.wanfangdata.com.cn/Periodical/kxtb200313008 [18] 郭唐永, 李欣, 谭业春, 等.中国新一代流动卫星激光测距系统研究[J].武汉大学学报·信息科学版, 2006, 31(1): 31-34 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200601009 Guo Tangyong, Li Xin, Tan Yechun, et al. New Generation Mobile SLR System TROS-I in China[J]. Geomatics and Information Science of Wuhan University, 2006, 31(1): 31-34 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb200601009 [19] 李博峰, 沈云中, 楼立志.基于等效残差的方差-协方差分量估计[J].测绘学报, 2010, 39(4): 349-354 http://d.old.wanfangdata.com.cn/Periodical/chxb201004004 Li Bofeng, Shen Yunzhong, Lou Lizhi. Variance-Covariance Component Estimation Based on the Equivalent Residuals[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(4): 349-354 http://d.old.wanfangdata.com.cn/Periodical/chxb201004004 [20] Yang Y X. Robust Bayesian Estimation[J]. Bull Geod, 1991, 65(2): 145-150 http://d.old.wanfangdata.com.cn/Periodical/dianzixb201309005 [21] Yang Y X. Robust Estimation for Dependent Observations[J]. Manuscripta Geodaetica, 1994, 19:10-17 doi: 10.1214-aos-1176347266/ [22] Yang Y X, Cheng M K, Shum C K, et al. Robust Estimation of Systematic Errors of Satellite Laser Range[J]. Journal of Geodesy, 1999, 73(7):345-349 doi: 10.1007/s001900050252 [23] 李晓杰, 郭睿, 吴杉, 等.不同姿态控制模式下的北斗卫星定轨策略研究[J].武汉大学学报·信息科学版, 2019, 44(10): 1 465-1 471 http://www.cnki.com.cn/Article/CJFDTotal-WHCH201910007.htm Li Xiaojie, Guo Rui, Wu Shan, et al. Orbit Determining Strategy Analysis for BeiDou Satellite in Different Attitude Control Modes[J]. Geomatics and Information Science of Wuhan University, 2019, 44(10): 1 465-1 471 http://www.cnki.com.cn/Article/CJFDTotal-WHCH201910007.htm [24] 程鹏飞, 文汉江, 刘焕玲, 等.卫星大地测量学的研究现状及发展趋势[J].武汉大学学报·信息科学版, 2019, 44(1): 48-54 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201901005 Cheng Pengfei, Wen Hanjiang, Liu Huanling, et al. Research Situation and Future Development of Satellite Geodesy[J]. Geomatics and Information Science of Wuhan University, 2019, 44(1): 48-54 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201901005 -