非合作目标激光测距预报的实时修正方法研究
张训方1,2, 赵雪1,2, 李荣旺1,3, 李祝莲1,3     
1. 中国科学院云南天文台, 云南 昆明 650011;
2. 中国科学院大学, 北京 100049;
3. 中国科学院空间目标与碎片观测重点实验室, 江苏 南京 210008
摘要: 非合作目标的激光测距预报一般是基于双行根数(TLE)外推出来的,往往有较大偏差,对激光测距的成功率有较大影响。结合空间目标的轨道理论和实测的数据分析,预报的偏差主要是预测模型外推的空间目标在运行轨道上的平近点角与实际平近点角存在偏差。根据非合作目标在望远镜跟踪视场中的脱靶量,利用相关算法可以找到一个最佳的时间根数偏差量修正空间目标的平近点角。经过修正,空间目标的视位置偏差得到改善,距离偏差能够从几百米减小到几十米,提升了预期回波到达时刻的准确度,可以给单光子探测器提供更高精度的距离门控,提高测距成功率。
关键词: 实时修正    时间根数偏差    非合作目标    轨道预报    
Research on Real-time Correction Method of Laser Ranging Prediction of Non-cooperative Target
Zhang Xunfang1,2, Zhao Xue1,2, Li Rongwang1,3, Li Zhulian1,3     
1. Yunnan Observatories, Chinese Academy of Sciences, Kunming 650011, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China;
3. Key Laboratory of Space Object and Debris Observation, Chinese Academy of Sciences, Nanjing 210008, China
Abstract: Prediction for laser ranging of non-cooperative target is extrapolated from the Two Line Elements (TLE) normally, which has relatively large deviation and affects the success rate of laser ranging mission. Based on the orbital theory and the measured data of space target, the deviation of the forecast is mainly due to the deviation of mean anomaly from the prediction model and the real mean anomaly of the spatial target in the orbit. According to the target's miss distance in the field of view of tracking telescope, the relative algorithm can be used to find the optimal time element deviation to correct the mean anomaly of the spatial target. After modification, the visual position deviation of the space target is improved, the distance deviation can be reduced from several hundred meters to tens of meters, which makes the echo arrival time expectation more accurate, and can provide a more precise range gate control for the single-photon detector, and improves the mission success rate of laser ranging.
Key words: Real-time correction    Time element deviation    Noncooperative target    Orbit prediction    

激光测距技术是一项综合技术,涵盖激光、电子、天文测量和卫星轨道计算等多个学科领域,在实际工作中是一个非常复杂的过程,受目标特性、望远镜指向误差、距离、接收孔径大小等多种因素的影响[1-2]。在众多影响因素中,空间目标的轨道预报准确性和稳定性直接影响测距望远镜发出的激光脉冲能否持续到达其表面,接收端的单光子探测器能否及时探测到回波光子对激光测距的成功率有很大影响[3]

卫星激光测距系统广泛采用了距离门控技术,即时间滤波技术。该技术能有效防止噪声的干扰,是在预期回波到达时刻之前一瞬间给单光子探测器送开门脉冲,这样使受到噪声干扰的概率减到最小,使淹没在众多噪声中的单光子探测成为可能[4]

非合作目标的激光测距预报,即载入跟踪望远镜的轨道预报基于双行根数(Two Line Elements, TLE),经过SGP4 (Simplified General Perturbations Version 4)轨道衍生器外推出来的,在对目标进行跟踪时产生的较大偏差主要体现在望远镜视场中的位置偏差以及测站与目标之间的距离偏差。视位置偏差过大导致激光脉冲不能击中目标;测站和目标的预报距离偏差过大,使预期回波到达时刻与实际回波到达时刻出现大的偏差,导致单光子探测器的距离门开启期间探测不到回波光子。

http://mobile.celestrak.com/NORAD/documentation/spacetrk.pdf

本文提出非合作目标激光测距预报实时修正方法,先利用双行根数原始预报对可见的非合作目标跟踪几分钟,获取其在望远镜跟踪视场中的脱靶量数据信息,根据这些信息解算时间根数偏差,并基于该时间根数偏差对轨道预报进行修正后,重新载入望远镜继续跟踪目标。激光测距卫星(例如Ajisai, Beaconc, Starlette, Stella等)具有双行根数和国际激光测距服务机构的官方预报星历(Consolidated Prediction Format, CPF)两种预报,文[5]将中低轨卫星的CPF星历与利用全球激光跟踪站的观测数据定轨出的精密轨道进行了比较,分析CPF星历精度在米级以内;文[6]分析了SGP4模型处理近地目标时,定轨精度在百米量级。为了排除望远镜指向误差等测距系统误差的影响,在理论上证明该方法的有效性,本文把中低轨道高度的激光测距卫星视为非合作目标进行研究分析,把某一时间段可见的激光卫星的CPF预报作为参考值,即模拟观测值。目标的脱靶量由同一时刻双行根数预报和CPF预报中的方位角和高度角偏差量代替,然后进行时间根数偏差的解算和轨道预报的修正。

1 激光测距预报的生成 1.1 时间系统的转换

协调世界时(UTC)是一种修正过的原子时,空间目标的状态以及观测数据都使用UTC作为输入、输出的时间参数。格林尼治平恒星时(Greenwich sidereal time, GST) θG是坐标系转换的一个关键参数,要由协调世界时计算出θG,详见文[7]。

1.2 TEME坐标系与站心地平坐标系之间的转换

利用预测模型外推出来的空间目标位置信息是在真赤道平春分点坐标系(True Equator Mean Equinox Frame, TEME)中表达的,而激光测距望远镜使用的轨道预报是在站心地平坐标系中表达的,所以要实现由真赤道平春分点坐标系向站心地平坐标系的转换,如图 1

$ {\rm{图 1中, }}{R_{\rm{z}}}\left( {{{\bar \theta }_{\rm{G}}}} \right) = \left( {\begin{array}{*{20}{c}} {{\rm{cos}}{{\bar \theta }_{\rm{G}}}}&{{\rm{sin}}{{\bar \theta }_{\rm{G}}}}&0\\ {-{\rm{sin}}{{\bar \theta }_{\rm{G}}}}&{{\rm{cos}}{{\bar \theta }_{\rm{G}}}}&0\\ 0&0&1 \end{array}} \right); $ (1)
图 1 TEME坐标系向站心地平坐标系转换示意图 Fig. 1 Schematic diagram of the conversion from the TEME system to the station local system

$\overrightarrow R $为空间目标在真赤道平春分点坐标系中的位置矢量;Rz(θG)为转换矩阵;θG为格林尼治平恒星时;$\overrightarrow r $为空间目标在地固坐标系中的位置矢量;$\overrightarrow ρ $为空间目标在站心地平坐标系中的位置矢量${\overrightarrow r ^{{\rm{sta}}}}$为激光测距望远镜在地固坐标系中的位置矢量;(MLT)为转换矩阵:

$ (MLT) = {R_{\rm{x}}}\left( {\frac{{\rm{ \mathsf{ π} }}}{2}-\phi } \right){R_{\rm{z}}}\left( {\frac{{\rm{ \mathsf{ π} }}}{2} + \lambda } \right), $ (2)

其中,

$ {R_{\rm{x}}}\left( {\frac{{\rm{ \mathsf{ π} }}}{2}-\phi } \right) = \left( {\begin{array}{*{20}{c}} 1&0&0\\ 0&{{\rm{sin}}\phi }&{{\rm{cos}}\phi }\\ 0&{-{\rm{cos}}\phi }&{{\rm{sin}}\phi } \end{array}} \right), $ (3)
$ {R_{\rm{z}}}\left( {\frac{{{\rm{ \mathsf{ π} }}}}{2} + \lambda } \right) = \left( {\begin{array}{*{20}{c}} {-{\rm{sin}}\lambda }&{{\rm{cos}}\lambda }&0\\ {-{\rm{cos}}\lambda }&{-{\rm{sin}}\lambda }&0\\ 0&0&1 \end{array}} \right), $ (4)

ϕλ分别为激光测距望远镜所在位置的纬度和经度;${R_{\rm{x}}}\left({\frac{{\rm{ \mathsf{ π} }}}{2}-\phi } \right)$${R_{\rm{z}}}\left({\frac{{\rm{ \mathsf{ π} }}}{2} + \lambda } \right)$均为转换矩阵。

1.3 光行差改正

当对中低轨空间目标进行激光测距时,测站和目标之间往往具有几千千米的距离,在激光脉冲的上行和下行期间,空间目标相对测站发生了较大变化,因此对于激光测距一般用发射时刻作为时标生成跟踪文件,即望远镜的实时指向有个前指量。针对本文的应用,由于下行光行差的原因,目标可见对应的预报应该滞后,即将激光接收方向作为模拟观测值。CPF预报中,前指量的算法如下:

对于某一时刻T,空间目标和测站在地固坐标系中的位置矢量为$\overrightarrow r $(X, Y, Z)和${\overrightarrow r ^{{\rm{sta}}}}$(X0Y0Z0),则测站和空间目标之间的相对距离$D = \left| {\overrightarrow r-{{\overrightarrow r }^{{\rm{sta}}}}} \right|$。在不考虑光行差的情况下,激光脉冲发射方向为T时刻空间目标的预报位置,激光脉冲在上行光路传播的时间Tout=D/c (c为光速),空间目标在站心地平坐标系中的位置矢量$\overrightarrow ρ $(ρx, ρy, ρz)。在考虑光行差的情况下,激光脉冲发射方向为TB=T+Tout时刻空间目标的预报位置。空间目标在站心地平坐标系中的位置矢量${\overrightarrow \rho _{\rm{B}}}$(ρBx, ρBy, ρBz)。此处需要特别注意的是地固坐标系是一个旋转坐标系,计算中需要进行相应改正。这样就得到上行光路的光行差引起的激光发射前指量:

$ {\rm{d}}A = A-{A_{\rm{B}}} = {\rm{ta}}{{\rm{n}}^{-1}}\left( {-\frac{{{\rho _y}}}{{{\rho _x}}}} \right) - {\rm{ta}}{{\rm{n}}^{ - 1}}\left( { - \frac{{{\rho _{{\rm{B}}y}}}}{{{\rho _{{\rm{B}}x}}}}} \right), $ (5)
$ {\rm{d}}E = E-{E_{\rm{B}}} = {\rm{ta}}{{\rm{n}}^{-1}}\frac{{{\rho _z}}}{{\sqrt {{\rho _x}^2 + {\rho _y}^2} }}-{\rm{ta}}{{\rm{n}}^{ - 1}}\frac{{{\rho _{{\rm{B}}z}}}}{{\sqrt {{\rho _{{\rm{B}}x}}^2 + {\rho _{{\rm{B}}y}}^2} }}. $ (6)

利用空间几何相关知识,下行光路的光行差引起的激光接收滞后量为dA,dE

2 时间根数偏差的解算以及预报的修正 2.1 理论表述

图 2中,阴影部分为地心天球坐标系O-XYZ的基本平面(XY平面),接近J2000.0平赤道面,X轴指向接近J2000.0平春分点,Z轴指向地球的北极。r为空间目标的轨道半径,$\overrightarrow h $为空间目标轨道平面的法向量,$\overrightarrow N $为通过轨道升交点的矢量,$\overrightarrow e $为通过轨道近地点的矢量。Ω为从X轴方向起量的轨道升交点方向$\overrightarrow N $的经度;i为轨道面与基本平面的夹角;ω为近地点幅角;θ为极坐标变量。空间目标的轨道极坐标方程[8]

$ r = \frac{{a\left( {1-{e^2}} \right)}}{{1 + e{\rm{cos}}\left( {\theta-\omega } \right)}}, $ (7)
图 2 环绕地球运行的空间目标的轨道示意图 Fig. 2 Schematic diagram representing the orbit of earth-circling spatial target

其中,a为半长径;e为偏心率;近地点幅角ω和极坐标变量θ都是从升交点方向起量的。当有摄动时,椭圆轨道随时间变化,升交点方向也在变化,ω从变化的升交点方向起量,而极坐标变量θ仍从一个定义的不变方向起量。

平近点角M定义为

$ M = n\left( {t-\tau } \right), $ (8)

从近地点开始起量,其中,n为平运动角速度。空间目标的平近点角是一个与时间相关的参数,能够确定对应时刻目标在椭圆轨道上的位置。当t=τ时,M=0,相应地,空间目标在近地点处,r达到极小值,τ是空间目标过近地点的时刻。

空间目标在轨运行时会受到多种摄动力的影响,包括地球的非球形引力、日月引力、大气阻力、太阳辐射压力等,对处于中低轨道的空间目标,受到大气阻力的影响比较大。双行根数考虑了地球非球形引力、日月引力的长期和周期摄动影响,以及大气阻力模型产生的引力共振和轨道衰退。双行根数是“平均”根数,它用特定的方法去掉了周期扰动项,预测模型SGP4用同样的方法重构了周期扰动项。在观测中,空间目标往往在沿迹方向与望远镜的指向中心有提前或者滞后的偏差。该偏差的主要原因是SGP4模型采用的大气模型简化为B*一个阻力参数,沿迹方向的偏差量可以简化归因于时间根数τ的偏差。

如果在跟踪目标时,适当调整τ的值,可以减小空间目标与望远镜指向中心的偏移量,修正空间目标与测站的距离量,进而提高跟踪的稳定性。

时间根数τ偏差的解算算法如下,记空间目标在站心地平坐标系中的方位角为A,高度角为E,离测站距离值为s,实测的序列记为(ti, Ai, Ei, si), i=1, 2, 3, …, n,在计算预报时适当调整τ的值,令τ + Δτ取代τ,Δτ是调整量,空间目标在真赤道平春分点坐标系中相应时刻点的位置${{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\rightharpoonup$}} \over R} }_i}$(Xi, Yi, Zi, ti + Δτ)代替${{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\rightharpoonup$}} \over R} }_i}$ (Xi, Yi, Zi, ti):

$ \overrightarrow r \left( t \right) = {R_{\rm{z}}}(\overline \theta {_{\rm{G}}}\left( t \right))\overrightarrow R {\rm{ }}(t + \Delta \tau ), \overrightarrow \rho \left( t \right) = \left( {MLT} \right)[\overrightarrow r \left( t \right)-{\overrightarrow r ^{{\rm{sta}}}}]; $

由此可得预报序列(ti, A′i, E′i, s′i; Δτ),i=1, 2, 3, …, n。记ΔAi=A′i-Ai,ΔEi=E′i-Ei,当Δτ的选值使

$ RM{S^*} = \sqrt {\frac{{{\sum _i}{{(\Delta {A_i}{\rm{cos}}{E_i})}^2} + {\sum _i}\Delta {E_i}^2}}{n}} $ (9)

为跟踪时间段上的极小值时,Δτ为预报的时间根数偏差。上述问题是一个复杂的非线性问题,直接求解比较困难,考虑到时间根数偏差精确到毫秒,因此一个简单直接的方法就是在一定范围内进行搜索。

当空间目标的指向信息得到修正时,距离信息也得到了相应的修正,进而能够提高单光子探测器的距离门控的精度,增大接收回波光子的概率。

2.2 实例分析

Ajisai, Beaconc, Starlette, Stella等激光测距卫星位于近地轨道,与大部分的非合作目标的轨道高度相仿,空间环境类似,能够很好地替代仅有双行根数预报的非合作目标进行研究分析。

本文以Stella激光测距卫星为例进行分析。Stella卫星的轨道高度为792.3 × 803.1 km,处于非合作目标的密集区。当对空间目标进行激光测距时,测站和目标之间往往具有上千千米的距离,在激光脉冲的上行和下行期间,空间目标相对测站发生了较大变化,指向滞后量如图 3,横轴为日积秒,纵轴分别为方位角偏差和高度角偏差。从图 3可以看出,滞后量波动幅度比较大,最高时可达十几角秒,是解算时间根数偏差时不可忽视的因素。

图 3 激光卫星Stella某一过境弧段中望远镜接收回波光子应有的指向滞后量 Fig. 3 The pointing lag of the telescope during receiving the echo photon in certain passing segment of Stella

图 4显示了在Stella的某次过境观测弧段,分别用前1 min,2 min,3 min以及全弧段跟踪数据进行时间根数搜索的结果,图中横轴为时间根数偏差Δτ,纵轴为对应的RMS*,时间根数搜寻范围为[-100, 200]。由图 4可见,无论用前1 min,2 min,3 min还是全弧段跟踪数据进行搜寻,RMS*都有一个极小值,对应时间根数偏差分别为26 ms,27 ms,29 ms,32 ms并无显著差异。对于其它的激光测距卫星,根据不同时长的跟踪弧段,搜索的时间根数偏差结果如表 1

图 4 激光测距卫星Stella分别以过境弧段的前1、2、3分钟和整个弧段的跟踪数据搜索出的时间根数偏差结果 Fig. 4 The time element deviations of Stella calculated by using the tracking data of previous 1, 2, 3 minutes and the entire passing segment
表 1 若干激光测距卫星根据不同时长的跟踪数据解算出的时间根数偏差 Table 1 The time element deviations of some Laser Ranging Stallites by using different duration of the tracking data
激光卫星 轨道高度/km 过境时长 前一分钟/ms 前两分钟/ms 前三分钟/ms 全弧段/ms
Ajisai 1 479 × 1 497 12 min 20 s 35 35 35 34
Beaconc 928.6 × 1 305 8 min 30 s 30 28 26 23
Starlette 805.5 × 1 107 6 min 9 s 36 38 40 43

通过表 1可以看出,上述3个激光测距卫星在过境可见的时间段内,用前1 min,2 min,3 min以及整个弧段解算的时间根数偏差无显著差异。根据搜索时间根数偏差的思想可知,利用的跟踪数据段越长,解算的结果越精确,利用全弧段跟踪数据解算的结果在理论上最精确。考虑到中低轨目标过境时间短等问题,选择前1 min数据解算的时间根数偏差适用于整个弧段。

激光卫星Stella的轨道预报经过相应的修正后,视位置偏差明显减小,且波动幅度也很小。图 5 (a)中,实线表示Stella的双行根数原始预报中的方位角与CPF参考值的残差随观测时间的变化,下方的虚线分别表示双行根数预报经过26 ms,27 ms,29 ms,32 ms修正的方位角与CPF参考值的残差随观测时间的变化。类似地,图 5 (b)中,实线表示Stella的双行根数原始预报中的高度角与CPF参考值的残差随观测时间的变化,中间的虚线分别表示双行根数预报经过26 ms,27 ms,29 ms,32 ms修正的高度角与CPF参考值的残差随观测时间的变化。图 6是以ΔAcosE为横坐标,ΔE为纵坐标画出了整个跟踪弧段以0 ms (未修正)和26 ms(前一分钟)分别修正后的视位置偏差量的变化。综合以上可以看出,双行根数预报经过时间根数偏差改正后,视位置偏差量明显减小,且在某一固定值附近小范围波动,波动范围在几角秒之内。

图 5 Stella的轨道预报分别以0 ms (未修正),26 ms,27 ms,29 ms,32 ms为时间根数偏差修正后与参考轨道在高度角(a)和方位角(b)的偏差量 Fig. 5 The deviations of Stella′s orbit prediction corrected by 0ms, 26ms, 27ms, 29ms, 32ms compared with the reference orbit in azimuth angle and elevation angle
图 6 Stella的轨道预报分别以0 ms (未修正),26 ms为时间根数偏差修正后与参考轨道的偏差量 Fig. 6 The deviations of Stella′s orbit prediction corrected by 0ms, 26ms compared with the reference orbit

表 2展示了处于近地轨道的4颗激光卫星的双行根数预报在未修正及基于前一分钟,整个弧段解算出的时间根数偏差分别修正后与CPF参考轨道在方位角、高度角的偏差均值以及相应的位置偏差RMS*$\overline {\Delta {A_0}{\rm{cos}}{E_0}}, \overline {\Delta {E_0}} $是双行根数原始预报与参考轨道的方位角、高度角在整个过境弧段的偏差均值。$\overline {\Delta {A_1}{\rm{cos}}{E_1}}, \overline {\Delta {E_1}}, \overline {\Delta {A_2}{\rm{cos}}{E_2}}, \overline {\Delta {E_2}} $分别是基于前一分钟、整个弧段解算的时间根数偏差进行修正后的轨道预报与CPF参考轨道在方位角、高度角的偏差均值。可以看出后两者的位置偏差RMS*比双行根数原始预报的要小,即更接近实测运行轨迹。后两者在方位角、高度角的偏差均值也很接近,更进一步说明可以用前一分钟的跟踪数据解算出的时间根数偏差代替整个弧段的跟踪数据解算的时间根数偏差进行轨道预报修正。

表 2 若干激光测距卫星修正后的轨道预报在方位角、高度角的偏差均值以及相应的RMS* Table 2 The mean deviations and corresponding RMS in azimuth angle and elevation angle of the corrected orbital predictions of some laser ranging satellites
激光星 $\overline {\Delta {A_0}{\rm{cos}}{E_0}} $ $\overline {\Delta {E_0}} $ $\overline {\Delta {A_1}{\rm{cos}}{E_1}} $ $\overline {\Delta {E_1}} $ $\overline {\Delta {A_2}{\rm{cos}}{E_2}} $ $\overline {\Delta {E_2}} $ RMS0* RMS1* RMS2*
Ajisai -14.730 1 15.160 2 -20.168 6 4.003 7 -20.012 4 4.322 0 19.206 6 13.580 2 13.582 6
Beaconc -13.451 4 4.064 5 -9.397 2 -8.887 3 -10.350 1 -5.871 8 12.174 0 8.769 7 4.160 2
Starlette -69.325 8 15.458 1 -62.306 6 -0.363 6 -60.936 5 -3.439 9 89.690 7 84.207 6 87.890 5
Stella 13.286 7 18.568 2 4.530 2 5.710 9 2.503 2 2.746 0 13.693 9 1.721 3 2.403 6

图 7可以看出,激光卫星Stella的距离预报分别以26 ms,27 ms,29 ms,32 ms为时间根数偏差修正后,距离偏差由最大的二百多米减小到几十米,并且变化范围更稳定。

图 7 激光卫星Stella的轨道预报以不同的时间根数偏差修正后,距离预报值与参考值的偏差量 Fig. 7 The distance deviations of Stella′s orbit prediction corrected by different time elememt compared with reference values

表 3展示了处于近地轨道的4颗激光卫星的轨道预报在未修正及基于前一分钟、整个弧段解算出的时间根数偏差分别修正后相较于CPF参考轨道在预报距离方面的偏差均值D0, D1, D2以及相应的RMS,即RMS0, RMS1, RMS2。可以看出,预报距离偏差均值由几百米减小到了几十米,甚至几米;均方根值RMS也减小了一个量级,说明偏差稳定性显著增强。上述结果说明基于时间根数偏差修正双行根数预报能够提高单光子探测器距离门控的精度。

表 3 若干激光测距卫星修正后的轨道预报在预报距离方面的偏差均值以及相应的RMS Table 3 The mean deviations and corresponding RMS in predicted range of the corrected orbital predictions of some laser ranging satellites
激光星 D0 D1 D2 RMS0 RMS1 RMS2
Ajisai -215.89 -49.78 -54.52 78.80 9.68 11.60
Beaconc -154.19 8.72 -29.21 46.16 9.11 4.35
Starlette -214.18 -4.56 36.19 56.31 14.84 7.67
Stella -180.97 -37.51 -4.99 62.42 18.58 8.67
3 结论

本文介绍了一种非合作目标激光测距预报实时修正的方法。在搜索时间根数偏差的过程中,以RMS*取得极小值时的Δτ值为时间根数偏差,用于空间目标的双行根数预报修正。用不同时长的跟踪数据搜寻出的时间根数偏差无显著差异,用前一分钟数据解算的时间根数偏差可以适用于整个弧段。修正后的双行根数预报更接近空间目标的实测运行轨迹且预报距离偏差值均由几百米减小到了几十米,甚至几米,能够提高单光子探测器距离门控的精度。

本文选取若干个空间目标的双行根数预报经过修正后,视位置值和预报距离值的偏差得到了明显的改善。

参考文献
[1] 翟东升, 汤儒峰, 李祝莲, 等. 激光脉冲特性对漫反射激光测距系统精度的影响研究[J]. 天文研究与技术, 2016, 13(3): 326–332 DOI: 10.3969/j.issn.1672-7673.2016.03.008
[2] 唐美荣, 和丽娟, 翟东升, 等. 激光测距系统回波光子分布特性研究[J]. 天文研究与技术, 2017, 14(1): 32–38
[3] 李语强, 李荣旺, 李祝莲, 等. 空间碎片激光测距应用研究[J]. 红外与激光工程, 2015, 44(11): 3324–3329 DOI: 10.3969/j.issn.1007-2276.2015.11.027
[4] 李祝莲, 熊耀恒, 何妙婵, 等. 云南天文台人造卫星激光测距系统原理[J]. 天文研究与技术——国家天文台台刊, 2008, 5(3): 248–252
[5] 陈国平, 何冰, 张志斌, 等. CPF星历精度分析[J]. 中国科学院上海天文台年刊, 2010(31): 35–44
[6] 韦栋, 赵长印. SGP4/SDP4模型精度分析[J]. 天文学报, 2009, 50(3): 332–339 DOI: 10.3321/j.issn:0001-5245.2009.03.010
[7] 黄珹, 刘林. 参考坐标系及航天应用[M]. 北京: 电子工业出版社, 2015.
[8] 刘林, 胡松杰. 航天器定轨理论与应用[M]. 北京: 电子工业出版社, 2015.
由中国科学院国家天文台主办。
0

文章信息

张训方, 赵雪, 李荣旺, 李祝莲
Zhang Xunfang, Zhao Xue, Li Rongwang, Li Zhulian
非合作目标激光测距预报的实时修正方法研究
Research on Real-time Correction Method of Laser Ranging Prediction of Non-cooperative Target
天文研究与技术, 2019, 16(1): 25-32.
Astronomical Research and Technology, 2019, 16(1): 25-32.
收稿日期: 2018-03-03
修订日期: 2018-04-18

工作空间