具体实施方式

下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。

针对搭载超高精度原子钟组(一般优于[email protected])的空间航天器与地面站的时间比对问题，为了获得更高的时间比对精度，需要更加稳定的测量链路和高精度星地时间比对系统及方法的支撑。本发明的目的在于提出一种基于同频模式的高精度星地双向时间比对系统及方法，通过星载设备特殊的天线设置和采用同一频点收发信号，在无需考虑大气误差的情况下，利用适当的参数输入及合理的误差修正技术，该比对系统及方法能够实现ps级的时间比对精度。

本发明所述的一种基于同频模式的高精度星地双向时间比对系统，该系统包括：

地面设备。其向空间航天器的星载设备发射频点为f₁的微波信号，并接收空间航天器星载设备发射的频点为f₁微波信号。

星载设备。其接收地面设备发射的频点为f₁的微波信号，并向地面发射频点为f₁的微波信号。

在上述基于同频模式的高精度星地双向时间比对系统中，所述的地面设备包括：

产生所述f₁频点微波信号的地面微波时频载荷；

与所述地面微波时频载荷连接的微波时频传递地面天线，其一方面将所述微波信号发射至空间航天器，另一方面接收从空间航天器上发射的所述不同频段的微波信号；

地面数据传输系统，其接收所述微波测距信号下传的伪距观测数据；

与所述地面数据传输系统连接的星地时频链路控制中心，其计算得到所述空间航天器钟和地面站钟的相对钟差。

在上述基于同频模式的高精度星地双向时间比对系统中，所述星载设备包括：

产生所述f₁频点微波信号的空间时频载荷；

与所述空间微波时频载荷连接的微波时频传递对地天线A，其将所述微波信号发射至地面站；

与所述空间微波时频载荷连接的微波时频传递对地天线B，其接收从地面站发射的所述不同频段的微波信号；

与所述空间微波时频载荷连接的精密定轨接收机，其接收来自卫星的信号，并进行空间航天器的定轨。

空间数据传输系统，其向地面发射所述的微波测距信号下传的观测数据、轨道数据和钟差数据等。

本发明所述的一种基于同频模式的高精度星地双向时间比对方法，其特征在于方法步骤如下：

步骤1，通过上述时间比对系统的地面设备时频传递天线和空间航天器星载设备的发射天线A和分别向对方发射同一频点的微波测距信号实现地面站和空间航天器之间的伪距测量，获得f₁频点的下行观测伪距和上行观测伪距。

步骤2，将步骤1获得的f₁频点下行伪距和上行伪距配对组成双向单程测量伪距组合，然后计算获得地面站和空间航天器时频设备之间的星地双向相对钟差。

步骤3，计算并扣除双向相对钟差计算时的各项误差，实现星地高精度时间同步；其中，各项误差包括：设备时延误差、运动时延误差(空间距离不对称误差和钟差非同时误差)、周期性相对论效应误差和引力时延误差。

本发明的进一步改进在于，步骤1具体包括：

在同频模式下的一次测量过程中，通过地面站和空间航天器双方相互发出测距信号供对方测量，然后交换测量数据，进行解算、配对，从而得到双向单程测量伪距组合；

其中，双向单程测量伪距组合为地面站和空间航天器相互同时收或发的两个单向测量的组合；

本发明的进一步改进在于，由于上下行伪距频点和路径几乎一致，双向时间同步方法可消除大部分的大气时延误差，因此可不做大气时延误差修正。

本发明的进一步改进在于，由于精密定轨处理估计的是空间航天器质心位置，使用的空间航天器的精密轨道也为质心轨道，而空间航天器的伪距相位测量距离均为相位中心，因此在数据处理时需对空间站的质心坐标进行相位中心修正，利用提前标定的相位中心结果将精密定轨处理结果分别修正到天线A和B的相位中心坐标值，表达如下：

其中，R_S是空间航天器的质心坐标，R_A和R_B分别为空间站对地天线A和B的相位中心位置矢量，F₁和F₂分别代表空间航天器质心修正到对地天线相位中心A和B转换表达。

本发明的进一步改进在于，步骤1中，进行时间比对的双方分别为地面站G和空间航天器S，为降低收发信号的同频干扰，空间航天器S采用不同的微波时频天线A 和B对地收发信号，因此在同频模式下，双向单程测量伪距组合表达为：

式(1)中，c为真空光速，t^Snd和t^Rcv分别为发射和接收时刻，和分别为f₁频点的下行(AG代表传播路径，星载设备天线A发射信号，地面站G接收信号)测量伪距和f₁频点的上行(GB代表传播路径，地面站G发射信号，星载设备天线B接收信号)测量伪距；R_A和R_B分别为空间航天器对地天线A和B的相位中心位置矢量，x_S和x_G分别为空间航天器和地面站钟差，δ^Snd和δ^Rcv为信号发射通道和接收通道的硬件传输时延，δ^rel为周期性相对论效应等效时延，δ^OF为相位中心偏移等效时延，δ^g是引力时延，ε为测距噪声。

本发明的进一步改进在于，步骤2中，实现星地双向时间同步的双方分别为地面站G和空间航天器S，由于空间航天器S的上下行信号收发天线相位中心不一致，且信号的空间传输延迟和卫星的高速运动，星地双向伪距测量值中包含了不同时刻的卫星位置与钟差信息，为此需要将双向伪距中包含的时空信息归算至同一时刻t₀。假设地面站和空间航天器之间的钟差为Δt，则t₀时刻的空间航天器与地面站的相对钟差表达式为：

式(2)中，Δd为空间距离不一致误差修正量，dx(t,t₀)为t时刻到t₀时刻的钟差非同时误差修正量。空间距离非对称造成的双向空间距离修正量由空间航天器的精密轨道计算，钟差非同时造成的钟差修正量由空间航天器的预报钟差计算，硬件传输时延定期标定，大气延迟、周期性相对论效应、引力时延以及相位中心偏移可依托数学模型计算。

本发明的进一步改进在于，设备时延误差通过精确标定的方法进行补偿。

本发明的进一步改进在于，由于空间航天器和地面站之间高速的相对运动导致了上下行信号传播路径和钟差并不一致，因此需对这种运动时延误差进行补偿，补偿方法包括：利用精密星历(航天器位置精度优于10cm，速度精度优于1mm/s)和精密钟差文件，对不同时刻和天线位置的同频双向单程观测伪距进行时标归算和运动时延误差修正，距离补偿和钟差补偿计算分别如下式(3)和(4)：

本发明的进一步改进在于，对于星载原子钟的相对论效应包括标称频率偏差(NFO)和周期性相对论效应，其中标称频率偏差可通过在航天器上天前进行提前校准消除，而周期性相对论效应可利用模型消除。由于星地双向时间比对中的周期性相对论效应，不具有路径对称性，必须依赖星历进行补偿，而补偿误差最终将会直接传递到时间比对结果中，因此在考虑ps级时间精度时需充分考虑J₂项及更高阶项周期性相对论误差的影响，对于空间航天器和地面站之间的周期性相对论修正模型为：

其中，W₀＝6.969290134×10^-10·c²为地球大地水准面上的重力势；GM为地球引力常数；a为空间航天器运行的半长轴；R为地球非球形引力势扰动；r为原子钟所处的地心距；τ为原子钟原时。

本发明的进一步改进在于，需消除从空间航天器到地面站接收机的信号传输由于广义相对论效应的作用导致的几何传输时延和引力传输时延。由于地面站和空间航天器坐标系都是地心惯性系(ECI)，因此可忽略双向几何传输时延的不一致。然而，在皮秒级的时间传递精度下，需要考虑地球中心引力引起的引力时延，修正模型如下：

其中，ρ_G和ρ_S分别为地面测站G和航天器S到地心的距离，为地面测站G到航天器S的地心的距离，ECI下可忽略地球自转效应，但是如果处在地心地固系则需考虑地球自转效应。

本发明的进一步改进在于，结合式(2)进行时标归算和各种误差修正后，利用同一参考时刻的双向单程观测数据，可以解算出相对钟差Δt。

本发明的进一步改进在于，能够实现星地时间同步精度达到ps级。

本发明实施例的一种基于同频模式的高精度星地双向时间比对系统，请参阅图1。其中同频模式下星载设备的安装位置请参阅图2。

本发明实施例的一种基于同频模式的高精度星地双向时间比对方法，请参阅图3和图4，主要包含双向单程伪距差分解算钟差和运动时延误差修正，具体包括以下步骤：

(1)通过f₁频点的微波信号双向单程测量可实现空间航天器与地面站之间的伪距测量请参阅图5，双向单程测量的一次测量过程中，通过空间航天器和地面站双方互发测距信号供对方测量，然后交换测量数据，进行解算，配对，从而得到相对测量值。

图5中，需要实现时间同步的双方分别为空间航天器S和地面站G，航天器S的对地天线A在时刻发射频点为f₁的微波信号，该信号经过空间航天器S发射设备时延、空间传播时延、其他附加时延，并经过地面站接收设备的设备时延后，最终在地面站G的时刻被检测，得到信号下行伪距观测值同理，地面站G在时刻发射频点为f₁的微波信号，空间航天器S的对地天线B在时刻检测到该信号，得到信号上行伪距观测值则一次双向单程伪距观测方程为：

其中，c为真空光速，t^Snd和t^Rcv分别为发射和接收时刻，和分别为f₁频点的下行测量伪距和f₁频点的上行测量伪距；R_A和R_B分别为空间站对地天线A和B在地心惯性系中的相位中心位置矢量，x_S和x_G分别为空间站和地面站钟差，δ^Snd和δ^Rcv为信号发射通道和接收通道的硬件传输时延，δ^rel为周期性相对论效应等效时延，δ^OF为相位中心偏移等效时延，δ^g是引力时延，ε为测距噪声。

(2)利用双向单程观测伪距计算空间航天器与地面站之间的相对钟差

将式(1)中双向伪距中包含的时空信息归算至同一时刻t₀，然后将归算后的上下行伪距差分可以得到空间航天器和地面站的在t₀时刻的相对钟差值Δt，计算表达式如下：

式中，Δd为空间距离不一致修正量，dx(t,t₀)为t时刻到t₀时刻的钟差非同时修正量。

(3)对运动时延误差进行修正

为避免同频信号干扰，星载设备的对地发射天线A和对地接收天线B位置之间存在一个固定距离，这将导致信号传播路径和钟差不一致误差的增大，为满足ps级时间同步精度要求，因此需充分考虑运动时延误差的影响。

对于运动时延误差的修正，这里利用精密星历(航天器位置精度优于10cm，速度精度优于1mm/s)和精密钟差文件，对不同时刻和天线位置的同频双向单程观测伪距进行时标归算和运动时延误差修正，距离补偿和钟差补偿计算分别如下式(3)和(4)，具体流程请参阅图6。

(4)扣除各项误差(除运动时延误差)，实现超高精度时间同步

为了提高时间同步的精度，还需要对公式(2)中包含的各项误差项进行扣除，各项误差的扣除方法分别如下：

1)设备时延误差

不考虑路径的情况下，设备时延误差仅与设备零值有关。对于已知的两个稳定设备，该值将是一个确定的量，可以通过精确标定来减小或消除其影响。

2)相对论效应误差

由于空间航天器和地面站之间存在高速的相对运动，星载设备受到广义相对论和狭义相对论的作用导致测距信号在传播过程中存在偏差，因此需要对相对论误差进行补偿。

对于标称频率偏差(NFO)，可采用提前校准的方法消除。假设星载原子钟的基频为f₀，则在发射前将频率调整为f：

其中，W₀＝6.969290134×10^-10·c²为地球大地水准面上的重力势；GM为地球引力常数；a为空间航天器运行的轨道半长轴；R为地球非球形引力势扰动；r为原子钟所处的地心距；τ为原子钟原时。

对于周期性相对论效应可采用模型法消除，具体模型如下：

通过(2)和(6)可消除双向时间同步过程中的周期性相对论效应误差。

3)引力时延

在ps级的是时间同步精度下，仅需要考虑地球中心引力引起的引力时延，对这项误差的补偿模型如下：

其中，ρ_G和ρ_S分别为地面测站G和航天器S到地心的距离，为地面测站G到航天器S的地心的距离，ECI下可忽略地球自转效应，但是如果处在地心地固系则需考虑地球自转效应。

4)其他误差

硬件噪声等其他项误差可采用目前成熟的硬件方案或算法处理即可，不作为发明的重点内容介绍。

请参阅图3至图6，本发明实施例中，假定实现时间同步的两方分别为LEO卫星和地面站，则以LEO卫星和地面站时间的时间同步为例进行说明，在同频模式下，假定二者的上下行信号以同发或同收的体制进行测距，其LEO星载设备的发射天线A 和jieshou天线B时间的距离为d cm，则实现星地超高精度时间同步的具体步骤如下：

步骤1，LEO卫星S在时刻发射频点为f₁(～30GHz)的微波信号，该信号经过LEO卫星S发射设备时延、空间传播时延、其他附加时延，并经过地面站接收设备的设备时延后，最终在地面站G的时刻被检测，得到地面站G的伪距观测值同理，地面站G在时刻发射频点为f₁(～30GHz)的微波信号，LEO卫星S在时刻检测到该信号，得到LEO卫星S的伪距观测值则LEO卫星和地面站各自接收到的伪距观测值表示为：

其中，c为真空光速，t^Snd和t^Rcv分别为发射和接收时刻，和分别为f₁频点的下行测量伪距和f₁频点的上行测量伪距；R_A和R_B分别为空间站对地天线A和B在地心惯性系中的相位中心位置矢量，x_S和x_G分别为空间站和地面站钟差，δ^Snd和δ^Rcv为信号发射通道和接收通道的硬件传输时延，δ^rel为周期性相对论效应等效时延，δ^OF为相位中心偏移等效时延，δ^g是引力时延，ε为测距噪声。

步骤2，将上下行观测伪距归算至同一参考时刻的测距值，并进行配对求解

利用LEO卫星的事后精密星历数据，将不同时刻的伪距观测值换算到同一参考时刻的双向测距值。利用同一时刻的双向测距值就可以重新组成一组新测量方程，表达式如下：

步骤3，误差修正

利用各项误差扣除方法，运动时延误差、设备时延误差、相对论效应误差、引力时延等误差进行修正和扣除。其中，运动时延误差消除需借助事后精密星历文件和星载设备天线搭载情况构建误差修正模型。

步骤4，计算LEO和地面站之间的相对钟差Δt

通过(7)进行差分运算和误差模型修正后，解算得到LEO卫星和地面站之间的双向钟差为Δt：

以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。