一种低轨小卫星之间的分级时钟同步方法 |
|||||||
申请号 | CN202110085543.8 | 申请日 | 2021-01-22 | 公开(公告)号 | CN112904705B | 公开(公告)日 | 2022-05-17 |
申请人 | 重庆邮电大学; | 发明人 | 赵国锋; 石东; 徐川; 邢媛; 刘俊斌; | ||||
摘要 | 本 发明 属于通信技术领域,特别涉及一种低轨 小卫星 之间的分级时钟同步方法,包括以IEEE 802.11协议作为星间链路通信协议,根据低轨小卫星群的拓扑,卫星 位置 信息、速度信息和运行状态信息,通过最优选择 算法 选出最佳主时钟卫星;以最佳主时钟卫星为根 节点 ,将卫星同步区域划分为多个时钟同步域,同时为每个时钟同步域分配唯一的同步域标识;通过双向时钟测量方法获得从时钟卫星与主时钟卫星之间的星间链路传播时延;通过时钟传递同步信息的方法对从时钟卫星进行时钟修正;本发明利用时间戳传递的方式,精确测量主从同步卫星之间的链路时延,提高主从卫星的时钟同步 精度 ,采用分级同步的思想,扩大小卫星之间的同步范围。 | ||||||
权利要求 | 1.一种低轨小卫星之间的分级时钟同步方法,其特征在于,具体包括以下步骤: |
||||||
说明书全文 | 一种低轨小卫星之间的分级时钟同步方法技术领域[0001] 本发明属于通信技术领域,特别涉及一种低轨小卫星之间的分级时钟同步方法。 背景技术[0002] 近年来,随着空间信息网络的发展,新业务尤其是时间敏感型业务的不断涌现,对信息传输与处理提出了新的需求。为保证时间敏感业务的确定性传输,时钟同步技术是关键。现有的星间时钟同步技术包括GPS定时同步技术、双向时间比对技术等,这些时钟同步方法虽然能满足星间卫星高精度时钟同步的需求,但这些方法无一例外都需要搭载高精度原子钟或需要全球导航卫星系统 (GNSS,Global Navigattion Satellite System)的支持,其造价过于昂贵、不易部署且维护难度大。对于数量多、体积较小的低轨小卫星群,更倾向于一种结构简单易于部署的时钟同步方案;且出于安全性考虑,如有特殊飞行任务的军事侦察编队小卫星群,更需要一种内部时钟同步协议。 [0003] IEEE 802.1AS作为时间敏感网络(TSN,Time Sensitive Networking)中的时钟同步协议,工作在全双工以太网的数据链路层,通过时间戳传递的方式实现高精度的时钟同步。IEEE 802.1AS协议结构简单,扩展性强,易于部署且同步精度可精确到纳秒级,为自动化工厂、汽车控制、音视频传输等众多领域提供高精度时钟同步服务。但现有的IEEE 802.1AS协议将链路传播时延测量过程和同步信息传递过程绑定在一起,针对无线信道广播的特性,直接部署到无线链路中会加大无线信道的开销,影响系统整体性能和时钟同步精度。而目前国内对于IEEE 802.1AS协议的关注更多的是停留在有线网络侧时钟同步的研究和实现上,将IEEE 802.1AS协议应用到无线网络中仍处于理论研究阶段且文献较少。 发明内容[0004] 为了解决上述问题,本发明提出一种低轨小卫星之间的分级时钟同步方法,具体包括以下步骤: [0005] S1、根据小卫星群星间链路拓扑情况、卫星的位置信息、运行速度信息和运行状态信息,选取最佳主时钟卫星,GPS/北斗授时卫星为最佳主时钟卫星提供授时服务,为低轨小卫星群提供高精度参考时钟源,并以最佳主时钟卫星为根节点,生成时钟同步生成树; [0006] S2、从时钟卫星和主时钟卫星通过双向时间测量的方式,在从时钟端获知测量所需时间戳,计算星间链路的传播时延,其中从时钟卫星包括中继卫星和普通从时钟卫星,主时钟卫星包括最佳主时钟和中继卫星;中继卫星在链路传播时延测量阶段作为从时钟卫星,在同步信息时作为主时钟卫星; [0007] S3、进行同步信息,最佳主时钟卫星向无线信道广播同步信息,从时钟卫星接收同步信息并根据时钟修正算法修正当前节点的时钟,如果从时钟节点是中继卫星,则中继卫星在修正当前节点的时钟后向下一时钟同步域传递时钟信息。 [0008] 进一步的,选取最佳主时钟卫星的过程包括:在信标帧阶段扩展信标帧结构,扩展部分包括自身ID、卫星状态信息Role、最佳主时钟权值系数Weight,最佳主时钟权值系数越大成为最佳主时钟卫星的可能性越高,最佳主时钟权值系数通过卫星位置信息、速度信息和运行状态信息进行量化。 [0009] 进一步的,通过量化卫星位置信息、速度信息和运行状态信息获取最佳主时钟权值系数的过程包括: [0010] 计算当前卫星最大传输距离范围内的邻居卫星节点与当前卫星之间的平均距离; [0011] 计算当前卫星最大传输距离范围内的邻居卫星节点与当前卫星之间的平均速度; [0013] 根据当前卫星的平均距离、平均速度以及状态模型进行加权融合获取最佳主时钟权值系数,表示为: [0014] [0015] 其中,α为位置信息的量化的权值, 为当前卫星最大传输距离范围内的邻居卫星节点与当前卫星之间的平均距离,σd为当前卫星最大传输距离范围内的邻居卫星节点与当前卫星之间的标准偏差;β为卫星速度信息的量化的权值, ni为当前卫星最大传输距离范围内的邻居卫星节点数量,vk表示当前卫星最大传输距离范围内的第k个邻居卫星节点的相对速度,vmax当前卫星最大传输距离范围内的邻居节点的最大相对速度;λ为卫星运行状态信息的量化的权值,Πh表示采用马尔科夫过程对当前卫星运行状态建模。 [0016] 进一步的,当前卫星最大传输距离表示为: [0017] [0018] 其中,Smax为当前卫星最大传输距离;λ为载波波长,PTX为发射功率,GTX和 GRX分别为发射天线增益和接收天线增益,k为玻尔兹曼常量,SPS为噪声带宽, M为每个发射符号数量。 [0019] 进一步的,采用马尔科夫过程对当前卫星运行状态建模Πh表示为: [0020] [0021] 其中,Ve表示卫星从工作状态变为故障状态的速率,Vr表示卫星故障状态时的响应速率,Vf表示卫星失效速率, 表示卫星从故障状态到失败状态的速率,λr表示卫星自愈速率,λf表示卫星从失效状态到重新开始正常工作的速率。 [0022] 进一步的, [0023] 以最佳主时钟为根节点,同步域号为0,广播层级发现数据包,其中层级发现数据包中包括数据包的标识符和自己的同步域号; [0024] 根节点的邻居节点接收层级发现数据包,首先判断自己是否有同步域号,没有则标记层级数据包中的同步域; [0025] 在标记自己的同步域号后,当前卫星节点重新封装并广播层级发现数据包,其中新层级发现数据包中的同步域号为当前节点的同步域号加一; [0026] 根据先到先标记的原则,一个卫星标记同步域号后不再接受其它卫星的标记信息。 [0027] 进一步的,步骤S2中在从时钟端获知测量所需时间戳时,基于IEEE 802.1AS 的时钟同步协议工作在数据链路层;通过IEEE 802.11的MAC层的物理层媒体访问控制状态机所在的管理实体MLME来获取高精度硬件时间戳。 [0028] 进一步的,计算星间链路的传播时延包括ClockMaster实体和ClockSlave 实体,用于完成链路传播时延测量和同步信息传递过程中数据的发起和接收处理。 [0030] 每个时钟同步域中的从时钟向无线信道广播Pdelay_request链路传播时延测量请求帧,其中在该帧中封装当前卫星的同步域号,并在发起端记录发起时刻t1; [0031] 主时钟卫星的ClockMaster实体收到并解析Pdelay_request广播帧,首先判断Pdelay_request帧中的同步域号是否与当前节点的域号相同,如果不同则丢弃,相同则在接收端记录接收时刻t2,保存Pdelay_request帧的源MAC地址; [0032] 主时钟完成Pdelay_request帧的处理后,主时钟卫星在t3时刻封装时刻t2和t3在Pdelay_response链路传播时延测量响应帧中,以步骤2中记录的源MAC地址为目的MAC地址向无线信道单播Pdelay_response帧; [0033] 从时钟卫星接收并处理Pdelay_response帧,获得时刻t2和t3,并记录接收时刻t4; [0034] 从时钟卫星计算链路传播时延,即: [0035] 其中,r为主从卫星节点的时钟频率偏移比率。 [0036] 进一步的,最佳主时钟卫星的ClockMaster实体向无线信道广播同步信息,时钟信息传递过程包括: [0037] 最佳主时钟卫星的ClockMaster实体在Sync同步信息传递帧中封装当前卫星的原始时间戳,同步域号,时钟频率偏移比率和当前卫星的时钟矫正信息 correctionField,最佳主时钟卫星的correctionField为0,时钟频率偏移比率为1,向无线信道广播Sync帧; [0038] 从时钟卫星的ClockSlave实体接收并解析Sync帧,首先判断帧中同步域号是否与当前节点的同步域号相同,不同则丢弃,相同则记录Sync帧中的时钟信息; [0039] 从时钟卫星根据Sync帧中的同步信息修正当前节点的时钟,即: [0040] Tsync=Torigin+Pdelay+correctionField(N)+Tsend+Tproce; [0041] 如果从时钟卫星是中继卫星,中继卫星在时钟同步后,更新当前卫星的矫正信息,当前卫星的矫正信息等于上一卫星节点的矫正信息加上与上一卫星节点的链路传播时延加上数据包的发送时延加上数据包的等待时延; [0042] 中继卫星向下一时钟同步域广播Sync帧,传递同步信息; [0043] 其中,Torigin为Sync帧中的原始时间戳,Pdelay为链路的传播时延, correctionField(N)为帧中的矫正信息,Tsend为数据包的发送时延,Tproce为数据包的处理时延。 [0044] 本发明的有益技术效果为: [0045] (1)系统部署简单,成本低,易于维护; [0046] (2)通过平均层数最小的最佳主时钟选取方式,减小时钟同步误差随着同步层数增加而增加,使同步性能达到最优; [0047] (3)对IEEE 802.1AS协议进行改进,适应星间同步的需要; [0048] (4)利用时间戳传递的方式,精确测量主从同步卫星之间的链路时延,同步误差,提高了主从卫星的时钟同步精度; [0050] 图1是本发明一种低轨小卫星之间的分级时钟同步方法的流程图; [0051] 图2是本发明中最佳主时钟选取流程图; [0052] 图3是本发明中链路传播时延测量流程图; [0053] 图4是本发明中同步信息传递流程图; [0054] 图5是本发明的扩展信标帧结构; [0055] 图6是本发明的低轨小卫星拓扑图; [0056] 图7是本发明的星间时钟同步生成树; [0057] 图8是本发明的主从同步卫星内部结构图; [0058] 图9是本发明的星间链路传播时延测量原理图; [0059] 图10是本发明的星间同步信息传递原理图。 具体实施方式[0060] 下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。 [0061] 本发明提出一种低轨小卫星之间的分级时钟同步方法,如图1所示,具体包括以下步骤: [0062] S1、根据小卫星群星间链路拓扑情况、卫星的位置信息、运行速度信息和运行状态信息,选取最佳主时钟卫星,GPS/北斗授时卫星为最佳主时钟卫星提供授时服务,为低轨小卫星群提供高精度参考时钟源,并以最佳主时钟卫星为根节点,生成时钟同步生成树; [0063] S2、从时钟卫星和主时钟卫星通过双向时间测量的方式,在从时钟端获知测量所需时间戳,计算星间链路的传播时延,其中从时钟卫星包括中继卫星和普通从时钟卫星,主时钟卫星包括最佳主时钟和中继卫星;中继卫星在链路传播时延测量阶段作为从时钟卫星,在同步信息时作为主时钟卫星; [0064] S3、进行同步信息,最佳主时钟卫星向无线信道广播同步信息,从时钟卫星接收同步信息并根据时钟修正算法修正当前节点的时钟,如果从时钟节点是中继卫星,则中继卫星在修正当前节点的时钟后向下一时钟同步域传递时钟信息。 [0065] 在本实施例中,将一种低轨小卫星之间的分级时钟同步方法分为最佳主时钟选取、链路传播时延测量和同步信息传递三个部分,具体包括以下步骤: [0066] (一)最佳主时钟选取 [0067] 根据当前卫星最大传输距离,将当前卫星通信范围内的卫星节点作为当前卫星的邻居卫星节点。根据小卫星的通信任务不同,星间通信距离一般在几十到上千公里不同,为满足星间长距离的通信需求,通过加大卫星发射功率和天线增益的方式来加大通信范围。本发明考虑小卫星电源、尺寸等资源的限制,设计卫星发射功率为30dbm,天线增益为 10dbi,计算可得最大通信距离为 500km,所述最大通信距离计算公式如下: [0068] [0069] 其中,λ为载波波长,PTX为发射功率,GTX和GRX分别为发射天线增益和接收天线增‑23益,k为玻尔兹曼常量等于1.381×10 J/K,SPS为噪声带宽,M为每个发射符号数量。 [0070] 为了实现最佳主时钟卫星的选取,在信标帧阶段扩展信标帧结构,如图5 所示,扩展部分包括自身ID,卫星状态信息Role,最佳主时钟权值系数Weight。其中最佳主时钟系数作为选择最佳主时钟的判断依据,最佳主时钟系数越大成为最佳主时钟卫星的可能性越高,最佳主时钟系数通过卫星位置信息、速度信息和运行状态信息进行量化。 [0071] 最佳主时钟的选取过程如图2,包括以下步骤: [0072] 卫星网络进行初始化,并进行邻居发现,获取邻居信息表,该表至少包括卫星位置信息、速度信息和运行状态信息; [0073] 根据邻居信息表中的信息计算最佳主时钟权值系数,并广播计算的权值并接收邻居的权值,判断当前卫星的权值系数是否为最大; [0075] 以选取的最佳主时钟卫星为根节点,该卫星的同步域号为0,广播层级发现数据包,接收到数据包的卫星判断是否已经标记同步域号,若没有标记则标记同步域号,若已经标记则向下一级广播层级发现数据包,并进行标记同步域号,标记完所有层级之后获得时钟同步生成树。 [0076] 在以上过程中,最佳主时钟权值系数通过量化卫星位置信息、速度信息和运行状态信息获取,具体包括以下步骤: [0077] 计算当前卫星最大传输距离范围内的邻居卫星节点与当前卫星之间的平均距离,若当前卫星的邻居卫星节点的位置集合为 则所有卫星的平均距离可表示为: [0078] [0079] 其中,ni表示第i个卫星的传输范围内有ni个邻居卫星节点,当前卫星的位置作为(xi,yi),j∈{1,2,…,ni}; [0080] 计算当前卫星最大传输距离范围内的邻居卫星节点与当前卫星之间的平均速度,若邻居节点的速度集合为为 那么平均相对速度 可表示为: [0081] [0082] 最大相对速度vmax可表示为: [0083] [0084] 将当前卫星运行状态分为工作状态、故障状态、响应状态和失败状态,所有状态构成一个有限的集合,采用马尔科夫过程对其建模,表示为: [0085] [0086] 其中,Ve表示卫星从工作状态变为故障状态的速率,Vr表示卫星故障状态时的响应速率,Vf表示卫星失效速率, 表示卫星从故障状态到失败状态的速率,λr表示卫星自愈速率,λf表示卫星从失效状态到重新开始正常工作的速率; [0087] 根据当前卫星的平均距离、平均速度以及状态模型进行加权融合获取最佳主时钟权值系数,表示为: [0088] [0089] 其中,Ps为卫星位置量化标准,Pv为卫星速度量化标准,Pz为卫星状态量化标准;α为位置信息的量化的权值, 为当前卫星最大传输距离范围内的邻居卫星节点与当前卫星之间的平均距离,σd为当前卫星最大传输距离范围内的邻居卫星节点与当前卫星之间的标准偏差;β为卫星速度信息的量化的权值,ni为当前卫星最大传输距离范围内的邻居卫星节点数量,vk表示当前卫星最大传输距离范围内的第k个邻居卫星节点的相对速度,vmax当前卫星最大传输距离范围内的邻居节点的最大相对速度;λ为卫星运行状态信息的量化的权值,Πh表示采用马尔科夫过程对当前卫星运行状态建模。 [0090] 根据最佳主时钟卫星为根节点,生成时钟同步生成树的过程如图6~7,具体包括以下步骤: [0091] 以最佳主时钟为根节点,同步域号为0,广播层级发现数据包,其中层级发现数据包中包括数据包的标识符和自己的同步域号; [0092] 根节点的邻居节点接收层级发现数据包,首先判断自己是否有同步域号,没有则标记层级数据包中的同步域; [0093] 在标记自己的同步域号后,当前卫星节点重新封装并广播层级发现数据包,其中新层级发现数据包中的同步域号为当前节点的同步域号加一; [0094] 根据先到先标记的原则,一个卫星标记同步域号后不再接受其它卫星的标记信息。 [0095] 本实施例中,基于IEEE 802.1AS的时钟同步协议工作在数据链路层,通过 IEEE 802.11的MAC层管理实体MLME来获取高精度硬件时间戳。 [0096] (二)链路传播时延测量 [0097] 如图8所示,主时钟卫星包括GPS接收模块,用于接收GPS的授时信号,保证低轨小卫星群拥有高精度参考时钟源,从时钟卫星搭载时钟晶振保证系统时钟源;ClockMaster模块和ClockSlave模块,用于完成主从时钟卫星的链路传播时延测量和同步信息传递过程中数据的发起和接收处理;Clock模块用于维护和统计实时更新的时钟信息;星间链路通信协议为IEEE 802.11。 [0098] 链路传播时延测量的流程如图3所示,图6展示在星间链路传播时延测量阶段,每个时钟同步域中从时钟卫星的ClockSlave实体发起链路传播时延测量请求,具体包括以下测量流程: [0099] 每个时钟同步域中的从时钟向无线信道广播Pdelay_request链路传播时延测量请求帧,其中在该帧中封装当前卫星的同步域号,并在发起端记录发起时刻t1; [0100] 主时钟卫星的ClockMaster实体收到并解析Pdelay_request广播帧,判断 Pdelay_request帧中的同步域号是否与当前节点的域号相同,如果不同则丢弃,相同则在接收端记录接收时刻t2,保存Pdelay_request帧的源MAC地址; [0101] 主时钟完成Pdelay_request帧的处理后,主时钟卫星在t3时刻封装时刻t2和t3在Pdelay_response链路传播时延测量响应帧中,以记录的源MAC地址为目的 MAC地址向无线信道单播Pdelay_response帧; [0102] 从时钟卫星接收并处理Pdelay_response帧,获得时刻t2和t3,并记录接收时刻t4; [0103] 从时钟卫星计算链路传播时延,所述计算公式如下: [0104] [0105] 其中,r为主从卫星节点的时钟频率偏移比率。 [0106] (三)同步信息传递 [0107] 如图7所示,在同步信息传递阶段,最佳主时钟卫星的ClockMaster实体向无线信道广播同步信息,所示时钟信息传递过程如图4,包括: [0108] 最佳主时钟卫星的ClockMaster实体在Sync同步信息传递帧中封装当前卫星的原始时间戳,同步域号,时钟频率偏移比率和当前卫星的时钟矫正信息 correctionField,最佳主时钟卫星的correctionField为0,时钟频率偏移比率为1,向无线信道广播Sync帧; [0109] 从时钟卫星的ClockSlave实体接收并解析Sync帧,首先判断帧中同步域号是否与当前节点的同步域号相同,不同则丢弃,相同则记录Sync帧中的时钟信息; [0110] 从时钟卫星根据Sync帧中的同步信息修正当前节点的时钟,所述的修正公式如下: [0111] Tsync=Torigin+Pdelay+correctionField(N)+Tsend+Tproce; [0112] 其中,Torigin为Sync帧中的原始时间戳,Pdelay为链路的传播时延, correctionField(N)为帧中的矫正信息,Tsend为数据包的发送时延,Tproce为数据包的处理时延; [0113] 如果从时钟卫星是中继卫星,中继卫星在时钟同步后,更新当前卫星的矫正信息,当前卫星的矫正信息等于上一卫星节点的矫正信息加上与上一卫星节点的链路传播时延加上数据包的发送时延加上数据包的等待时延; [0114] 中继卫星向下一时钟同步域广播Sync帧,传递同步信息。 [0115] 本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:ROM、RAM、磁盘或光盘等。 |