CN111884748B - 应用于卫星共视法时间同步中实时数据交互的编码方法 - Google Patents

应用于卫星共视法时间同步中实时数据交互的编码方法 Download PDF

Info

Publication number
CN111884748B
CN111884748B CN202010721494.8A CN202010721494A CN111884748B CN 111884748 B CN111884748 B CN 111884748B CN 202010721494 A CN202010721494 A CN 202010721494A CN 111884748 B CN111884748 B CN 111884748B
Authority
CN
China
Prior art keywords
time
data
seconds
clock difference
delta
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010721494.8A
Other languages
English (en)
Other versions
CN111884748A (zh
Inventor
张思德
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Navcompass Science & Technology Co ltd
Original Assignee
Beijing Navcompass Science & Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Navcompass Science & Technology Co ltd filed Critical Beijing Navcompass Science & Technology Co ltd
Priority to CN202010721494.8A priority Critical patent/CN111884748B/zh
Publication of CN111884748A publication Critical patent/CN111884748A/zh
Application granted granted Critical
Publication of CN111884748B publication Critical patent/CN111884748B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J3/00Time-division multiplex systems
    • H04J3/02Details
    • H04J3/06Synchronising arrangements
    • H04J3/0602Systems characterised by the synchronising information used
    • H04J3/0605Special codes used as synchronising signal
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/185Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851Systems using a satellite or space-based relay
    • H04B7/18519Operations control, administration or maintenance

Abstract

本申请公开了一种应用于卫星共视法时间同步中实时数据交互的编码方法,包括:设定起始跟踪时刻;计算每秒单个或者组合伪距对应的钟差;在第i个起始跟踪时刻计算钟差;输出所有时刻不同间隔的钟差数据。本申请方法兼容了卫星共视统一的国际标准格式CGGTTS,又扩展了不同时间间隔的数据,可有效应用在更低精度的晶振和铷钟的远程驯服;提高用户在卫星共视中的可靠性。有利于时间频率领域的远程校准和计量的推广;本申请方法和卫星导航差分数据格式具有兼容性,可同时应用于高精度定位及高精度时钟源的系统中。

Description

应用于卫星共视法时间同步中实时数据交互的编码方法
技术领域
本发明属于卫星导航技术领域,具体涉及一种应用于卫星共视法时间同步中实时数据交互的编码方法。
背景技术
目前GPS或北斗同步时钟大部分采用卫星导航单向授时,为用户提供准确的时间频率。随着工业领域、5G通信领域、时间频率校准计量领域等对时间频率的可溯源性、安全性、高精度均有较大的需求。采用卫星共视法时间同步可解决这些问题,但这种双向授时法没有统一的数据交互格式,且当前采用的国际统一的卫星共视数据交互格式,CGGTTS格式采用16分钟间隔进行拟合,不能满足利用卫星共视法对铷钟、恒温晶振的实时驯服。
发明内容
为解决现有技术中的不足,本申请提出一种应用于卫星共视法时间同步中实时数据交互的编码方法,实现了统一一种卫星共视中实时的交互格式,并和现有的卫星导航差分格式有一定的兼容性,可实现差分定位、卫星共视驯服时钟同步实现。
一种应用于卫星共视法时间同步中实时数据交互的编码方法,包括如下:
设定起始跟踪时刻;
计算每秒单个或者组合伪距对应的钟差;
在第i个起始跟踪时刻计算钟差;
输出所有时刻钟差。
所述设定起始跟踪时刻,包括:
设定参考起始跟踪时刻为每天的0点,记作Time_Ref(0);
则第i个起始跟踪时刻Time(i)=Time_Ref(0)+i*△T
其中i=0、1、2、…;△T=300s、600s、…,△T为30秒的整数倍。
所述每秒单个或者组合伪距对应的钟差,计算方法按照国际计量局BIPM 制定的格式规范《CGGTTS-Version 2E:an extended standard for GNSS Time Transfer》中约定的单频用户和双频用户对钟差的计算。
当第i个起始跟踪时刻Time(i)小于等于60秒时,所述第i个起始跟踪时刻钟差,直接取瞬时钟差值。
当第i个起始跟踪时刻Time(i)大于60秒时,所述第i个起始跟踪时刻钟差计算方法如下:
将第i个起始跟踪时刻分为m组,每组间隔为n秒,则m×n=i;
每n秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的钟差数据作为中间值数据△ti(i=0,1…,m-1);
对中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
所述第i个起始跟踪时刻,i=10时,300秒间隔数据的钟差值:
每15秒的数据做一次最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti(i=0,1,…,19);
对20组中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
所述第i个起始跟踪时刻,i=20时,600秒间隔数据的钟差值:
每15秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti(i=0,1,…,39);
对40组中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
所述第i个起始跟踪时刻,i=32时,960秒间隔数据的钟差值:
剔除前120秒、最后60秒的数据,剩余780秒的数据;
对780秒的数据进行分组,每15秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti(i=0,1,…,51);
再对52组中间值数据△ti进行最小二乘法一次直线拟合,然后取中点时刻对应的数值△t作为钟差值。
本申请所达到的有益效果:
1、本申请方法兼容了卫星共视统一的国际标准格式CGGTTS,又扩展了不同时间间隔的数据,可有效应用在更低精度的晶振和铷钟的远程驯服中。经过实物验证,对铷钟的驯服已取得较好的结果,时间准确度可达到10ns(95%,溯源至UTC(k))。
2、本申请方法兼容北斗、GPS、GLONASS、Galileo等卫星系统,所有频点的信号,可提高用户在卫星共视中的可靠性。
3、本申请方法包含的共视数据的生成频率最高达1秒,使时间频率的远程溯源和RTK差分一样简单方便,有利于时间频率领域的远程校准和计量的推广;
4、本申请方法和卫星导航差分数据格式具有兼容性,可同时应用于高精度定位及高精度时钟源的系统中。
附图说明
图1本申请实施例的一种应用于卫星共视法时间同步中实时数据交互的编码方法;
图2本申请实施例的第i个起始跟踪时刻钟差计算方法流程图。
具体实施方式
下面结合附图对本申请作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本申请的保护范围。
一种应用于卫星共视法时间同步中实时数据交互的编码方法,如图1所示,包括如下:
步骤S101:设定起始跟踪时刻;
步骤S102:计算每秒单个或者组合伪距对应的钟差;
步骤S103:在第i个起始跟踪时刻计算钟差;
步骤S104:输出所有时刻钟差。
所述设定起始跟踪时刻,包括:
设定参考起始跟踪时刻为每天的0点,记作Time_Ref(0);
则第i个起始跟踪时刻Time(i)=Time_Ref(0)+i*△T
其中i=0、1、2、…;△T=300s、600s、…,△T为30秒的整数倍。
所述每秒单个或者组合伪距对应的钟差,计算方法按照国际计量局BIPM 制定的格式规范《CGGTTS-Version 2E:an extended standard for GNSS Time Transfer》中约定的单频用户和双频用户对钟差的计算。
当第i个起始跟踪时刻Time(i)小于等于60秒时,所述第i个起始跟踪时刻钟差,直接取瞬时钟差值。
当第i个起始跟踪时刻Time(i)大于60秒时,所述第i个起始跟踪时刻钟差计算方法如下:
步骤S201:将第i个起始跟踪时刻分为m组,每组间隔为n秒,则m×n=i;
步骤S202:每n秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti(i=0,1…,m-1);
步骤S203:对中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
所述第i个起始跟踪时刻,i=10时,300秒间隔数据的钟差值:
每15秒的数据做一次最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti(i=0,1,…,19);
对20组中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
所述第i个起始跟踪时刻,i=20时,600秒间隔数据的钟差值:
每15秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti(i=0,1,…,39);
对40组中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
所述第i个起始跟踪时刻,i=32时,960秒间隔数据的钟差值:
剔除前120秒、最后60秒的数据,剩余780秒的数据;
对780秒的数据进行分组,每15秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti(i=0,1,…,51);
再对52组中间值数据△ti进行最小二乘法一次直线拟合,然后取中点时刻对应的数值△t作为钟差值。
本实施例具体应用见如下表格:
表1为数据类型定义;
表1数据类型的定义
Figure BDA0002600179820000041
Figure BDA0002600179820000051
共视法实时数据交互格式:
表2:信息类型的格式
Figure BDA0002600179820000052
表3:信息类型列表
Figure BDA0002600179820000053
Figure BDA0002600179820000061
信息类型4010~4095中,每个数据域用下述的一个数据类型DFxxx表示。
表4数据域信息表
Figure BDA0002600179820000062
Figure BDA0002600179820000071
Figure BDA0002600179820000081
表5 4010-4049信息类型的信息头列表
Figure BDA0002600179820000082
表6 4010~4089信息类型数据体列表
Figure BDA0002600179820000083
说明:GNSS码伪距对应的信息,先按伪距类型排列,然后按卫星编号循环排列。如L1CA(PRN01)、L1P(PRN01)、L2P(PRN01)、L3P(PRN01);…… L1CA(PRN32)、L1P(PRN32)、L2P(PRN32)、L3P(PRN32);
表7 4090-4095信息类型的信息头列表
Figure BDA0002600179820000091
表8 4090-4095信息类型的信息体列表
Figure BDA0002600179820000092
本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述,但是本领域技术人员应该理解,以上实施示例仅为本发明的优选实施方案,详尽的说明只是为了帮助读者更好地理解本发明精神,而并非对本发明保护范围的限制,相反,任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。

Claims (4)

1.一种应用于卫星共视法时间同步中实时数据交互的编码方法,其特征在于,包括如下:
设定起始跟踪时刻;
计算每秒单个或者组合伪距对应的钟差;
在第i个起始跟踪时刻计算钟差;
输出所有时刻不同间隔的钟差;
所述设定起始跟踪时刻,包括:
设定参考起始跟踪时刻为每天的0点,记作Time_Ref(0);
则第i个起始跟踪时刻Time(i)=Time_Ref(0)+i*△T
其中i=0、1、2、…;△T=300s、600s、…,△T为300秒的整数倍;
所述每秒单个或者组合伪距对应的钟差,计算方法按照国际计量局BIPM制定的格式规范中约定的单频用户和双频用户对钟差的计算;
当所述第i个起始跟踪时刻Time(i)小于等于60秒时,所述第i个起始跟踪时刻钟差,直接取瞬时钟差值;
当第i个起始跟踪时刻Time(i)大于60秒时,所述第i个起始跟踪时刻钟差计算方法如下:
将第i个起始跟踪时刻分为m组,每组间隔为n秒,则m×n=i;
每n秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti,i=0,1…,m-1;
对中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
2.根据权利要求1所述的应用于卫星共视法时间同步中实时数据交互的编码方法,其特征在于:所述第i个起始跟踪时刻,i=10时,300秒间隔数据的钟差值:
每15秒的数据做一次最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti,i=0,1,…,19;
对20组中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
3.根据权利要求1所述的应用于卫星共视法时间同步中实时数据交互的编码方法,其特征在于:所述第i个起始跟踪时刻,i=20时,600秒间隔数据的钟差值:
每15秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti,i=0,1,…,39;
对40组中间值数据△ti进行最小二乘法一次直线拟合,取中点时刻对应的数值△t作为钟差值。
4.根据权利要求1所述的应用于卫星共视法时间同步中实时数据交互的编码方法,其特征在于:所述第i个起始跟踪时刻,i=32时,960秒间隔数据的钟差值:
剔除前120秒、最后60秒的数据,剩余780秒的数据;
对780秒的数据进行分组,每15秒的数据做最小二乘法二次曲线拟合,获取中点时刻对应的数据作为中间值数据△ti,i=0,1,…,51;
再对52组中间值数据△ti进行最小二乘法一次直线拟合,然后取中点时刻对应的数值△t作为钟差值。
CN202010721494.8A 2020-07-24 2020-07-24 应用于卫星共视法时间同步中实时数据交互的编码方法 Active CN111884748B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010721494.8A CN111884748B (zh) 2020-07-24 2020-07-24 应用于卫星共视法时间同步中实时数据交互的编码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010721494.8A CN111884748B (zh) 2020-07-24 2020-07-24 应用于卫星共视法时间同步中实时数据交互的编码方法

Publications (2)

Publication Number Publication Date
CN111884748A CN111884748A (zh) 2020-11-03
CN111884748B true CN111884748B (zh) 2023-02-03

Family

ID=73200244

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010721494.8A Active CN111884748B (zh) 2020-07-24 2020-07-24 应用于卫星共视法时间同步中实时数据交互的编码方法

Country Status (1)

Country Link
CN (1) CN111884748B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022192870A1 (en) * 2021-03-08 2022-09-15 Microchip Technology Incorporated Scalable common view time transfer and related apparatuses and methods

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181550A (zh) * 2014-08-19 2014-12-03 北京无线电计量测试研究所 一种基于北斗卫星导航系统的共视时间频率传递方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105607089B (zh) * 2016-01-18 2017-12-22 西安测绘研究所 一种基于伪距辅助载波相位的北斗共视时间传递的方法
CN111030774A (zh) * 2019-12-19 2020-04-17 北京无线电计量测试研究所 一种基于北斗卫星导航系统的实时共视数据处理方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181550A (zh) * 2014-08-19 2014-12-03 北京无线电计量测试研究所 一种基于北斗卫星导航系统的共视时间频率传递方法

Also Published As

Publication number Publication date
CN111884748A (zh) 2020-11-03

Similar Documents

Publication Publication Date Title
Gurtner et al. Rinex-the receiver independent exchange format-version 3.00
EP2331984B1 (en) Method and apparatus for a global navigation satellite system receiver coupled to a host computer system
CN109633701B (zh) 基于守时实验室资源的gnss定时接收机系统延迟校准方法
Odijk et al. PPP-RTK by means of S-system theory: Australian network and user demonstration
Steigenberger et al. Performance evaluation of the early CNAV navigation message
CN110196419B (zh) 用于gnss信号采集回放设备的伪距精度校准方法及系统
CN114994727B (zh) 一种实现高精度时间校准及卫星定位的设备
CN102004441B (zh) 自适应晶振频率守时方法
CN111884748B (zh) 应用于卫星共视法时间同步中实时数据交互的编码方法
CN103269262A (zh) 一种时间同步装置的守时方法
CN111191293B (zh) 基于区块链的北斗数据处理方法及装置
CN113110016A (zh) 共视数据生成方法、接收机和时间校准系统
CN116679329A (zh) 基于状态域空间的北斗时空基准精确表征方法和系统
Boriskin et al. The RTCM multiple signal messages: A new step in GNSS data standardization
CN107505637B (zh) 一种基于相对伪距的热启动定位方法
CN202533743U (zh) 时间分析仪测量基准系统
CN201266923Y (zh) Gps组合时间频率仪
CN211180593U (zh) 时频检测设备
CN111766614B (zh) 一种守时实验室联合组网的亚纳秒实时时间比对方法
Gurtner et al. The Receiver Independent Exchange Format
Gogoi et al. NavIC receiver clock offsets estimation with common view master clock method
WO2021077588A1 (zh) 接收机的冷启动控制方法、装置、终端设备及介质
Gogoi et al. NavIC Time Transfer With Satellite Common View: An Analysis on Choosing Different Reference Times
CN114397678A (zh) 一种基于实时精密单点定位技术的小型高精度频率源装置
Levine GPS and the legal traceability of time

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant