分布式卫星网络智能时钟同步方法

1.本技术涉及卫星时钟同步领域，特别是涉及一种分布式卫星网络智能时钟同步方法。


背景技术：

2.卫星网络时钟同步方式主要有gnss授时、基于结构式网络拓扑的时间传递和分布式时钟同步。目前，卫星网络时钟同步系统主要还是依靠当前的gnss，每颗卫星与gnss保持时间同步，从而实现整网的时间同步，然而gnss信号微弱易受干扰，并且要实现优于纳秒级的用户授时精度，系统将变得很复杂，成本也将急剧增加。基于结构式网络拓扑的时间传递，是通过主次节点间两两比对进行的，机理简单且收敛速度快，但在网络拓扑发生变化或节点失效时，需要重新规划时间传递路径，使得系统维护成本较高，鲁棒性和可扩展性相对不足，只适用于卫星数量较少且网络拓扑变化较小的情况。
3.因此，如何建立一种时钟同步方法以满足探索高动态卫星网络拓扑场景下大规模dss的应用需求成为了一个亟待解决的技术问题。
4.上述内容仅用于辅助理解本发明的技术方案，并不代表承认上述内容是现有技术。


技术实现要素：

5.本发明的主要目的在于提供了一种分布式卫星网络智能时钟同步方法，旨在解决现有技术无法实现满足探索高动态卫星网络拓扑场景下大规模dss的应用需求的技术问题。
6.为实现上述目的，本发明提供了一种分布式卫星网络智能时钟同步方法，所述方法包括：
7.获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率；
8.根据所述相对时钟速率估计所述目标节点的时钟速率补偿量以对所述目标节点的时钟速率进行补偿；
9.在对所述时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数；
10.在所述时钟读数同步误差小于预设阈值时，判定同步成功。
11.可选地，所述获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率的步骤，包括：
12.获取目标节点i与所述目标节点的邻居节点j；
13.在所述目标节点i与所述目标节点的邻居节点j进行双向通信之后，根据考虑时延的相对速率估计算法输入τi(ti(m)),τi(ti(m-1)),以获得相对时钟速率
14.可选地，所述在所述目标节点i与所述目标节点的邻居节点j进行双向通信之后，根据考虑时延的相对速率估计算法输入τi(ti(m)),τi(ti(m-1)),1)),以获得相对时钟速率的步骤，包括
15.在所述目标节点i与所述目标节点的邻居节点j进行双向通信之后，估计节
[0016][0017]
其中，t是两节点间多次测量的周期间隔；t
m-1
，tm分别表示第m-1和m次测量的参考时间；δτ
i,2
(tm)表示利用星历对两个卫星在时刻t
m-1
到时刻tm相对运动导致的时差补偿量；
[0018]
ρm∈(0,1)是iir滤波器系数。
[0019]
可选地，所述根据所述相对时钟速率估计所述目标节点的时钟速率补偿量以对所述目标节点的时钟速率进行补偿的步骤，包括：
[0020]
将逻辑时钟速率估计转换为时钟速率修正系数λ'i的估计：
[0021][0022][0023]
其中，γ表示分布式控制器参数；ki表示比例估计器的增益参数；k
p
表示积分估计器的增益参数；ni是与节点i建链测量的临接节点集合；ωi表示pi控制器状态参数。
[0024]
可选地，所述根据所述相对时钟速率估计所述目标节点的时钟速率补偿量以对所述目标节点的时钟速率进行补偿的步骤之后，还包括：
[0025]
目标节点i的时钟速率同步误差为：
[0026][0027]
系统的时钟速率误差表示为：
[0028][0029]
根据所述系统的时钟速率误差对同步误差进行约束。
[0030]
可选地，所述在对时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数的步骤，包括：
[0031]
在对时钟速率补偿结束后，可以得到每个卫星节点时钟相位偏差的一致性估计，
如下式。
[0032][0033]
结合时钟速率同步算法，可以更新每个节点的时钟相位
[0034][0035]
可选地，所述在对时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数的步骤之前，还包括：
[0036]
建立误差模型，t时刻时钟读数同步误差为：
[0037][0038]
每个节点i的时钟同步误差为：
[0039][0040]
其中为节点i的虚拟时钟，为节点j的参考时钟读数；
[0041]
则有：
[0042]
可选地，所述获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率的步骤之前，还包括：
[0043]
建立线性时钟模型可以表示时钟的真实时间特性，目标节点i在参考时间时刻t的物理时钟读数τi(t)可以描述为：
[0044]
τi(t)＝λit+μi+o(t)
[0045]
其中λi是时钟速率，μi是时钟初始偏差，o(t)表示时钟读数的高阶小量，对于高精度时钟，在较短测量周期内o(t)可以忽略；
[0046]
对每个卫星节点i的时钟读数进行线性变换，得到修正的时钟读数，建立一个虚拟逻辑时钟，逻辑时钟与本地时钟之间呈线性关系，可得逻辑时钟的表达式如下：
[0047][0048]
其中，是逻辑时钟读数；是节点i逻辑时钟速率；是节点i逻辑时钟初始偏差；λ'i是t时刻节点i的时钟速率修正系数；μ'i是t时刻节点i的时钟初始偏差修正量。
[0049]
此外，为实现上述目的，本发明还提出一种分布式卫星网络智能时钟同步装置，所述装置包括：
[0050]
相对时钟估计模块，用于获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率；
[0051]
速率补偿模块，用于根据所述相对时钟速率估计所述目标节点的时钟速率补偿量
以对所述目标节点的时钟速率进行补偿；
[0052]
偏差补偿模块，用于在对所述时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数；
[0053]
时钟同步模块，用于在所述时钟读数同步误差小于预设阈值时，判定同步成功。
[0054]
此外，为实现上述目的，本发明还提出一种计算机设备，所述计算机设备包括：存储器，处理器以及存储在所述存储器上并可在所述处理器上运行的分布式卫星网络智能时钟同步程序，所述分布式卫星网络智能时钟同步程序配置为实现如上文所述的分布式卫星网络智能时钟同步方法。
[0055]
此外，为实现上述目的，本发明还提出一种介质，所述介质上存储有分布式卫星网络智能时钟同步程序，所述分布式卫星网络智能时钟同步程序被处理器执行时实现如上文所述的分布式卫星网络智能时钟同步方法的步骤。
[0056]
本发明根据获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率；根据所述相对时钟速率估计所述目标节点的时钟速率补偿量以对所述目标节点的时钟速率进行补偿；在对所述时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数；在所述时钟读数同步误差小于预设阈值时，判定同步成功，通过每个节点通过相对速率估计算法估计出和邻居节点间时钟速率的相对误差，作为时钟速率补偿算法的初始输入参与反馈。由于时钟速率是一个缓慢时变的变量，通过比例积分估计器动态地跟踪缓慢时变的速率参数，法能够提升时钟速率的同步精度。每个节点根据本地时钟周期向邻居节点广播当前的本地时钟读数，并周期性地获得邻居节点传来的时钟读数，通过时钟读数补偿算法完成时钟读数补偿，所有节点的时钟同步到虚拟时钟，从而实现整个网络的同步。
附图说明
[0057]
图1是本发明实施例方案涉及的硬件运行环境的分布式卫星网络智能时钟同步设备的结构示意图；
[0058]
图2为本发明分布式卫星网络智能时钟同步方法第一实施例的流程示意图；
[0059]
图3为本发明分布式卫星网络智能时钟同步方法第一实施例卫星网络智能时钟同步的系统模型框架图；
[0060]
图4为本发明分布式卫星网络智能时钟同步方法第一实施例分布式时钟同步算法的示意图；
[0061]
图5为本发明分布式卫星网络智能时钟同步方法第一实施例的分布式卫星网络节点间双单向测量过程图。
具体实施方式
[0062]
应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0063]
参照图1，图1为本发明实施例方案涉及的硬件运行环境的分布式卫星网络智能时钟同步设备结构示意图。
[0064]
如图1所示，该分布式卫星网络智能时钟同步设备可以包括：处理器1001，例如中
央处理器(central processing unit，cpu)，通信总线1002、用户接口1003，网络接口1004，存储器1005。其中，通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(display)、输入单元比如键盘(keyboard)，可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如无线保真(wireless-fidelity，wi-fi)接口)。存储器1005可以是高速的随机存取存储器(random access memory，ram)存储器，也可以是稳定的非易失性存储器(non-volatile memory，nvm)，例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
[0065]
本领域技术人员可以理解，图1中示出的结构并不构成对分布式卫星网络智能时钟同步设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
[0066]
如图1所示，作为一种存储介质的存储器1005中可以包括操作系统、数据存储模块、网络通信模块、用户接口模块以及分布式卫星网络智能时钟同步程序。
[0067]
在图1所示的分布式卫星网络智能时钟同步设备中，网络接口1004主要用于与网络服务器进行数据通信；用户接口1003主要用于与用户进行数据交互；本发明分布式卫星网络智能时钟同步设备中的处理器1001、存储器1005可以设置在分布式卫星网络智能时钟同步设备中，所述分布式卫星网络智能时钟同步设备通过处理器1001调用存储器1005中存储的分布式卫星网络智能时钟同步程序，并执行本发明实施例提供的分布式卫星网络智能时钟同步方法。
[0068]
本发明实施例提供了一种分布式卫星网络智能时钟同步方法，参照图2，图2为本发明分布式卫星网络智能时钟同步方法第一实施例的流程示意图。
[0069]
本实施例中，所述分布式卫星网络智能时钟同步方法包括以下步骤：
[0070]
步骤s10：获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率。
[0071]
需要说明的是，在分布式卫星网络中，每个卫星节点均配备一个高稳晶振或原子钟，建立线性时钟模型可以表示时钟的真实时间特性，目标节点i在参考时间时刻t的物理时钟读数τi(t)可以描述为：
[0072]
τi(t)＝λit+μi+o(t)
[0073]
其中λi是时钟速率，μi是时钟初始偏差，o(t)表示时钟读数的高阶小量，对于高精度时钟，在较短测量周期内o(t)可以忽略。
[0074]
对每个卫星节点i的时钟读数进行线性变换，得到修正的时钟读数，建立一个虚拟逻辑时钟，逻辑时钟与本地时钟之间呈线性关系，可得逻辑时钟的表达式如下：
[0075][0076]
其中，是逻辑时钟读数；
[0077]
是节点i逻辑时钟速率；
[0078]
是节点i逻辑时钟初始偏差；
[0079]
λ'i是t时刻节点i的时钟速率修正系数；
[0080]
μ'i是t时刻节点i的时钟初始偏差修正量。
[0081]
分布式卫星网络智能时钟同步算法的目标是通过某种方式修正卫星网络中所有节点的物理时钟速率和时钟初始偏差，最终实现每个节点的逻辑时钟速率和读数渐近收敛到同一个虚拟参考时钟
[0082][0083]
其中，是参考时钟速率，是参考时钟初始偏差。
[0084]
因此，分布式卫星网络智能时钟同步需要设计一种算法，使得每个节点的逻辑时钟速率和读数渐近收敛，即：
[0085][0086][0087]
需要说明的是，天基智能网联系统是一个高速动态、快速响应的智能敏捷系统，需要实现高速数据传输、动态资源共享、分布式计算与分布式信号级协同等，离不开智能弹性的高精度时钟基准。高精度智能时钟基准建立技术，支持分布式相干探测、多源协同探测、分布式敏捷波形生成等应用。在弹性分散化星座编队协作互联系统中，星座编队之间相互组成虚拟mimo系统，由于发射天线相互分散形成分布式mimo系统。不同于传统mimo技术中发送天线和接收天线都属于同一终端，分布式星间链路的多收发器属于不同节点且相互独立，协作传输需要多节点信号时钟同步来完成。各个节点上的时钟和频率不智能和时变性都可以建模为相对于理想信道的时变信道特性。实现分布式波束形成的基本挑战是，独立的时钟漂移以及独立的发射机运动会产生高度非平稳的信道，如果来自接收器的反馈延迟比动态信道保持恒定的周期长，那么分布式波束形成性能就会急剧下降。因此为实施多收发器协同波束形成，首先需要进行星间链路多收发器节点的高精度分布式时钟预同步与预测。
[0088]
可以理解的是，传统星间链路的时钟同步协议主要是结构式算法，一般分为两步：首先在两个星间邻居节点之间建立同步；然后构造一个多跳拓扑结构使得此种同步可以一层一层地传递到所有星间链路节点。然而结构性的时钟同步协议在构建和维持树形或簇形结构上负载过重，另外根节点或网关节点的故障会导致大量与之连接的节点发散故障，协议的可扩展性和鲁棒性有所欠缺。基于多智能体协同智能控制理论的分布式时钟同步算法是完全分布式的，不需要依赖于既定的拓扑结构，也不依赖于根节点或者参考节点，具有较强的鲁棒性和可扩展性。智能问题是多智能体系统协调控制的一个基本问题，所谓智能是指随着时间的演化，多个智能体通过相互作用、相互协调使得感兴趣的状态或者输出达到相同。智能问题有三个基本要素，即智能体、智能体间用于信息交换的通信网络以及智能体的作用规则(或称为智能协议)。
[0089]
在具体实施中，本实施例采用图论作为星间测量网络拓扑的描述工具，构建测量网络拓扑模型其中v、ε和分别表示卫星节点、测量链路及其权重的集合。研究星间双向测量中的建链机制，构建无向测量链路和无向测量网络拓扑图。同时，针对星簇、同轨面星座、异轨面星座、大规模星群等典型分布式卫星应用系统的构型，定量描述卫
星测量网络拓扑图的特征及规律。
[0090]
在具体实施中，分布式卫星网络智能时钟同步方法利用星间时钟时钟测量与分布式控制构建卫星网络智能时钟同步的系统模型框架，如图3所示，将各个卫星与邻近卫星时钟的相对时差测量值δτ
ij
和相对频率测量值δf
ij
等作为系统输入，利用智能同步协议进行分布式计算，估计自身时钟的时间偏差和频率偏差等参数，并基于这些估计参数对星上逻辑时钟进行时钟补偿，通过多次迭代测量控制达到卫星网络逻辑时钟一致。基于多智能体智能控制理论，将卫星网络时钟作为多智能体系统，根据卫星时钟特性，构建并优化智能体动力学模型及其参数，设计适用于卫星网络时钟同步的智能协议。
[0091]
在具体实施中，分布式时钟同步算法考虑了时钟速率参数缓慢时变，采用一阶时钟模型，即考虑到每个时钟晶振频率受到自身老化以及环境因素的影响，是一个随时间缓慢变化的变量，分布式时钟同步算法的示意图如图4所示。由于本地时钟速率是一个缓慢时变的参数，而静态协同控制算法只能对静态(时不变)状态实现智能趋同。因此这里考虑动态智能协同控制算法。在动态智能算法中，每个节点时变的参考输入通过驱动比例积分估计器，动态地跟踪时变输入量的平均值。算法主要包括下面三个部分：相对速率估计算法、时钟速率补偿算法以及时钟读数补偿算法。算法的核心思想是每个节点通过相对速率估计算法估计出和邻居节点间时钟速率的相对误差，作为时钟速率补偿算法的初始输入参与反馈。由于时钟速率是一个缓慢时变的变量，通过比例积分估计器动态地跟踪缓慢时变的速率参数，法能够提升时钟速率的同步精度。最后，每个节点根据本地时钟周期向邻居节点广播当前的本地时钟读数，并周期性地获得邻居节点传来的时钟读数，通过时钟读数补偿算法完成时钟读数补偿。
[0092]
进一步地，所述所述获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率的步骤，包括：获取目标节点i与所述目标节点的邻居节点j；在所述目标节点i与所述目标节点的邻居节点j进行双向通信之后，根据考虑时延的相对速率估计算法输入τi(ti(m)),τi(ti(m-1)),以获得相对时钟速率
[0093]
进一步地，所述所述在所述目标节点i与所述目标节点的邻居节点j进行双向通信之后，根据考虑时延的相对速率估计算法输入τi(ti(m)),τi(ti(m-1)),以获得相对时钟速率的步骤，包括
[0094]
在所述目标节点i与所述目标节点的邻居节点j进行双向通信之后，估计节
[0095][0096]
其中，t是两节点间多次测量的周期间隔；t
m-1
，tm分别表示第m-1和m次测量的参考时间；δτ
i,2
(tm)表示利用星历对两个卫星在时刻t
m-1
到时刻tm相对运动导致的时差补偿量；
[0097]
∈(0,1)是iir滤波器系数。
[0098]
步骤s20：根据所述相对时钟速率估计所述目标节点的时钟速率补偿量以对所述目标节点的时钟速率进行补偿。
[0099]
进一步地，所述根据所述相对时钟速率估计所述目标节点的时钟速率补偿量以对
所述目标节点的时钟速率进行补偿的步骤，包括：
[0100]
将逻辑时钟速率估计转换为时钟速率修正系数λ'i的估计：
[0101][0102][0103]
其中，γ表示分布式控制器参数；ki表示比例估计器的增益参数；k
p
表示积分估计器的增益参数；ni是与节点i建链测量的临接节点集合；ωi表示pi控制器状态参数。
[0104]
进一步地，所述根据所述相对时钟速率估计所述目标节点的时钟速率补偿量以对所述目标节点的时钟速率进行补偿的步骤之后，还包括：
[0105]
目标节点i的时钟速率同步误差为：
[0106][0107]
系统的时钟速率误差表示为：
[0108][0109]
根据所述系统的时钟速率误差对同步误差进行约束。
[0110]
步骤s30：在对所述时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数。
[0111]
进一步地，所述所述在对时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数的步骤，包括：
[0112]
在对时钟速率补偿结束后，可以得到每个卫星节点时钟相位偏差的一致性估计，如下式。
[0113][0114]
结合时钟速率同步算法，可以更新每个节点的时钟相位
[0115][0116]
进一步地，所述在对时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数的步骤之前，还包括：
[0117]
建立误差模型，t时刻时钟读数同步误差为：
[0118][0119]
每个节点i的时钟同步误差为：
[0120][0121]
其中为节点i的虚拟时钟，为节点j的参考时钟读数；
[0122]
则有：
[0123]
步骤s40：在所述时钟读数同步误差小于预设阈值时，判定同步成功。
[0124]
需要说明的是，基于多智能体智能的分布式星间网络时间同步理论与协议，使得各个卫星节点能够更新自身的时钟信息，从而实现星间网络分布式智能时钟同步。该方法具有很强的适用性，能够覆盖几十公里到几万公里的距离范围，以及ns-ps级精度需求，适合天基动态多任务时空频分集利用场景的分布式时钟同步。
[0125]
进一步地，所述所述获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率的步骤之前，还包括：建立线性时钟模型可以表示时钟的真实时间特性，目标节点i在参考时间时刻t的物理时钟读数τi(t)可以描述为：
[0126]
τi(t)＝λit+μi+o(t)
[0127]
其中λi是时钟速率，μi是时钟初始偏差，o(t)表示时钟读数的高阶小量，对于高精度时钟，在较短测量周期内o(t)可以忽略；对每个卫星节点i的时钟读数进行线性变换，得到修正的时钟读数，建立一个虚拟逻辑时钟，逻辑时钟与本地时钟之间呈线性关系，可得逻辑时钟的表达式如下：
[0128][0129]
其中，是逻辑时钟读数；是节点i逻辑时钟速率；是节点i逻辑时钟初始偏差；λ'i是t时刻节点i的时钟速率修正系数；μ'i是t时刻节点i的时钟初始偏差修正量。
[0130]
在具体实施中，获取星历信息及星间网络拓扑，在此基础上建立星间双向测量链路。通过星间链路周期性地进行分布式卫星网络星间时频测量。
[0131]
每条星间链路每次通过双向测量获取4个时戳值，即时间测量值τ
ij
包含了{τ
i,1
，τ
j,2
，τ
j,1
，τ
i,2
}，分别表示节点i的发送时间、节点j的接收时间、节点j的发送时间和节点i的接收时间。每两个卫星节点间的一次建链可以获取多组测量值。星间双向测量原理如下：星间双向测量可以通过微波和激光，目前微波技术相对成熟，并且在bds等系统中已得到应用。这里利用微波扩频技术，两个卫星节点基于双单向测量构建时频同步系统，节点间双单向测量过程如图5所示。节点a在t1时刻发出信号，该信号是经过伪码调制的扩频信号，其伪码相位携带了发送时刻t1的信息(发送时戳)。发射系统的电路延迟为ta，无线电的空间延迟为τa，该信号经过节点b的电路处理延迟rb，在t2时刻被检测到(接收时戳)。由于信号在t1时刻发出，在t2被检测到，考虑到卫星存在相对运动，其相对距离是时间的函数，因此测量到的延迟τa是这t1和t2两个时刻的函数τa(t1,t2)，该延迟包括了所有空间路径上的延迟。节点b基于自身的时钟测量这两个时刻的差值为tb。节点b也在时刻t3时刻发出测距信号，从而提
供双单向的时间和距离测量。设发射系统的电路延迟为tb，无线电的空间延迟为τb，该信号经过节点a的电路处理延迟ra，在t4时刻被检测到。和节点a同样的道理，空间延迟为τb在卫星存在运动的情况下，其延迟是时间函数τb(t3,t4)。节点a测量出节点b发出信号到节点a接收到信号的时刻差值为ta。假设两个节点之间钟差是δt，根据上图中的时间关系可以得出，
[0132]
ta＝δt+tb+τb(t3,t4)+ra[0133]
tb＝-δt+ta+τa(t1,t2)+rb[0134]
得到：
[0135][0136]
其中，第一项是两次单向测量值的差值，是星间时差结果组成的主要部分；第二项反应了收发系统电路延迟的差异，与设备零值有关。对于两个同样的设备，该值将是一个确定的小量，可以通过标定的方法减小或消除，得到优于亚纳秒级精度的校正；第三项包含了双单向路径传输延迟的信息，当双向传播路径一致时(处于同一轨道或相对静止)，其差值为零，抵消了传播路径延迟的影响，可以得到较高的时间同步精度。当星间存在相对运动时，两个单向测量的传播路径延迟不等，即τa(t1,t2)≠τb(t3,t4)，该项是不能忽略的。
[0137]
进行历元归算，消除星间远距离传输、高速运动等因素导致的误差。
[0138]
分布式卫星网络与地面无线传感器网络不同，节点间距离非常远且相对运动速度很大。尤其是不同轨道面上卫星的高速运动会导致星间双向测量的传输路径存在严重的不对称，并且星间双向测量并不是同时进行的。针对此问题，可以进行历元归算，即利用星历中的位置与速度信息将两次测量数据计算归化到同一时刻，以便于对两次不同时刻的单向测量值进行比对和综合处理，消除远距离测量和高速运动导致的星间测量存在的较大误差。历元归算算法大致分为基于星历数据的逐项补偿法和基于星间伪距观测数据的拟合/插值法。伪距观测直接拟合法需要积累足够的观测数据才能进行，具有一定的滞后性。逐项补偿法在算法的灵活性和适应性上优于多项式拟合法。对于低轨星群，卫星时钟频率偏差一般小于1e-10，同步周期为秒级时归算精度优于0.1ns，若定轨精度在m级，归算误差在0.01ns以内。因此，通过历元归算将星间测量时戳归算到基于参考时钟的约定发送时刻，可以保证0.1ns的精度。基于星间双向测量机制的时间同步精度，主要依赖于星间伪距测量的精度(即时戳生成精度)。由于星间伪码扩频测距的精度可达亚纳秒甚至更高，那么经过设备零值标校和历元归算后，可以实现星间动态运动场景下的高精度时间同步。
[0139]
基于修正后的星间测量值，利用低通滤波器估计两个节点间的时钟相对速率
[0140]
因为绝对参考时间t不可知，所以卫星时钟速率的绝对值是无法获得的，只能通过星间测量获得卫星时钟之间的相对速率。利用连续两次星间测量的时戳值，估计节点j相对于节点i的相对时钟速率如下：
[0141]
[0142]
其中，t是两节点间多次测量的周期间隔；
[0143]
t
m-1
，tm分别表示第m-1和m次测量的参考时间；
[0144]
δτ
i,2
(tm)表示利用星历对两个卫星在时刻t
m-1
到时刻tm相对运动导致的时差补偿量；
[0145]
ρm∈(0,1)是iir滤波器系数。
[0146]
结合分布式卫星网络拓扑，每个卫星节点时钟利用所有相邻卫星节点的时钟相对速率，通过动态一致性估计器，估计时钟速率
[0147]
由于每个卫星节点时钟的真实时钟速率并不能直接获得，并且一致性时钟同步控制的目标是实现星群网络节点逻辑时钟的一致，可以将逻辑时钟速率实现星群网络节点逻辑时钟的一致，可以将逻辑时钟速率估计转换为时钟速率修正系数λ'i的估计。
[0148]
分布式卫星网络的分布式一致性时钟速率同步算法如下：
[0149][0150][0151]
其中，γ表示分布式控制器参数；
[0152]ki
表示比例估计器的增益参数；
[0153]kp
表示积分估计器的增益参数；
[0154]
ni是与节点i建链测量的临接节点集合；
[0155]
ωi表示pi控制器状态参数。
[0156]
利用时钟速率更新时钟读数，再通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数。
[0157]
经过时钟速率估计与修正后，各个节点的逻辑时钟频率达到同步，进一步利用星间链路测量值历元归算后得到的相对时差δτ
ij
(tm)，可以得到每个卫星节点时钟相位偏差的一致性估计，如下式。
[0158][0159]
结合时钟速率同步算法，可以更新每个节点的时钟相位
[0160][0161]
至此，分布式卫星网络中每个节点利用与临接节点的相对测量，完成了一次时钟频率和相位同步的更新。
[0162]
判断同步算法是否收敛，重复进行上述步骤迭代更新时钟参数
[0163]
通过判断时钟读数同步误差是否超过阈值，判断同步算法是否收敛，若误差较大重复进行上述步骤迭代更新时钟参数，若误差在一定范围内表明同步算法收敛，可以输出时钟同步结果。
[0164]
通过多次迭代收敛，可以实现分布式卫星网络时钟的频率同步与相位同步
[0165]
本实施例根据获取目标节点与所述目标节点的邻居节点，在所述目标节点与邻居节点进行双向通信后，通过相对速率估计算法获得所述目标节点与邻居节点的相对时钟速率；根据所述相对时钟速率估计所述目标节点的时钟速率补偿量以对所述目标节点的时钟速率进行补偿；在对所述时钟速率补偿结束后，通过相邻节点时钟读数的加权平均补偿修正本节点时钟读数；在所述时钟读数同步误差小于预设阈值时，判定同步成功，通过每个节点通过相对速率估计算法估计出和邻居节点间时钟速率的相对误差，作为时钟速率补偿算法的初始输入参与反馈。由于时钟速率是一个缓慢时变的变量，通过比例积分估计器动态地跟踪缓慢时变的速率参数，法能够提升时钟速率的同步精度。每个节点根据本地时钟周期向邻居节点广播当前的本地时钟读数，并周期性地获得邻居节点传来的时钟读数，通过时钟读数补偿算法完成时钟读数补偿，所有节点的时钟同步到虚拟时钟，从而实现整个网络的同步。
[0166]
此外，本发明实施例还提出一种介质，所述介质上存储有分布式卫星网络智能时钟同步程序，所述分布式卫星网络智能时钟同步程序被处理器执行时实现如上文所述的分布式卫星网络智能时钟同步方法的步骤。
[0167]
本发明分布式卫星网络智能时钟同步装置的实施例或具体实现方式可参照上述各方法实施例，此处不再赘述。
[0168]
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
[0169]
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
[0170]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如只读存储器/随机存取存储器、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。
[0171]
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。