CN106383438A - 一种基于滑动窗口时间扩展高精度驯钟方法 - Google Patents

一种基于滑动窗口时间扩展高精度驯钟方法 Download PDF

Info

Publication number
CN106383438A
CN106383438A CN201610998692.2A CN201610998692A CN106383438A CN 106383438 A CN106383438 A CN 106383438A CN 201610998692 A CN201610998692 A CN 201610998692A CN 106383438 A CN106383438 A CN 106383438A
Authority
CN
China
Prior art keywords
frequency
time slot
clock
fpga
digital
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.)
Granted
Application number
CN201610998692.2A
Other languages
English (en)
Other versions
CN106383438B (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.)
NANJING YINSHI SOFTWARE Co
Original Assignee
NANJING YINSHI SOFTWARE Co
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 NANJING YINSHI SOFTWARE Co filed Critical NANJING YINSHI SOFTWARE Co
Priority to CN201610998692.2A priority Critical patent/CN106383438B/zh
Publication of CN106383438A publication Critical patent/CN106383438A/zh
Application granted granted Critical
Publication of CN106383438B publication Critical patent/CN106383438B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G04HOROLOGY
    • G04RRADIO-CONTROLLED TIME-PIECES
    • G04R20/00Setting the time according to the time information carried or implied by the radio signal
    • G04R20/02Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明基于滑动窗口时间扩展高精度驯钟方法,设时间扩展窗口大小为N,参考时钟频率f0为FPGA工作时钟,以1PPS信号上升沿为基准,在FPGA中生成周期性的时隙计数器TSC,计数值从0累加到N‑1,然后循环;每个时隙关联一个64位计数器,在每一个f0的上升沿累加1,从本帧时隙开始计数,到下帧相应时隙开始结束计数;MCU每秒取出相应时隙的计数值并存入到队列中,队列首元素去除,新数据从尾部加入,窗口下滑;利用最小二乘法做数值滤波并获得中间均值,计算得到卫星模块与本地时钟频率差;在FPGA中设计数字DDS电路,根据MCU给出的频差值,生成准确的频率相位数据送入数模转换DA,最终生成驯服频率f1

Description

一种基于滑动窗口时间扩展高精度驯钟方法
技术领域
本发明属于电子测控领域,能够应用于卫星导航定位系统、电力高精度时间同步系统等等。
背景技术
对于基于全球定位系统GPS/北斗的导航定位和高精度授时装置来说,本身晶体钟/原子钟的频率准确度和稳定度对其性能影响很大。如果装置本地频率漂移造成3ns偏差,卫星定位误差就会超过1米。而且频差的影响结果是持续性的,系统实际性能影响更大。
整个卫星空间系统的运行是基于一个受地面控制管理的铯原子钟群,由于铯原子钟频率的稳定性,整个卫星星座系统的频率基准源非常稳定,这是GPS/北斗全球导航系统精确稳定运行的基础。为了达到较好的导航定位性能,地面的卫星接收装置通常采用驯服技术使本地晶体钟/原子钟与GPS/北斗频率一致。
现有的本地晶体钟/原子钟驯服技术主要采用频差测量压控驯服方法,如图1所示。整个本地钟驯服子系统由三个子模块组成:卫星接收模块、频率测量模块和可调校本地钟模块。卫星TOD串口报文每秒会输出准确的时间信息和定位信息,输出的秒脉冲1PPS信号含有卫星频率信息。频率测量模块能够测量出本地钟与卫星星座的频差,然后根据这些信息产生调整电压信号,修正本地钟的频率,从而实现本地钟与卫星星座铯钟频率的一致性。
上述本地钟驯服方法优点是实现了负反馈式的频率跟随,具有一定的精度和频率稳定度,提高了应用系统的定位精度和授时准确度,在大量工业装置和系统中得到应用。但是目前使用该方法也具有如下缺点:
(1)本地钟驯服精度不高,频率稳定度低。由于该方案采用压控方法,估测的频差数据转换成调控电压模拟量,而本地钟压控频率调整曲线不一定是线性的。实际应用中调整频率通常采用目标尝试法,调大了就往小里调一些,调小了就往大里调一些,反复调整修正,使频差逐渐趋向稳定。这种频率驯服机制的时间较长,精度不高。
(2)数字频差测算精度不高。这种基于计数法的频率测量精度受本地钟频率限制。利用一个较高的本地频率对卫星模块的1PPS进行计数,每周期的计数值之间的差异能够反映本地频率与卫星星座频率之间的频差。由于数字电路最高工作频率是有限制的,目前一般为200MHz左右,测量分辨率只有5ns。频差测量精度不高导致最后时钟的驯服精度也不高。
发明内容
发明目的:在现有架构的基础上,提出新方法,提高目前导航定位设备和时间同步装置中本地时钟的驯服方法的精度,从而大幅度提高时间同步装置的授时精度和守时的准确度和稳定性,提高车载导航设备的定位精度,对于工业应用具有重要意义。
本发明要解决的技术问题主要有如下两个:
(1)解决本地时钟和卫星星座时钟之间频差测量精度问题。通常卫星模块输出的信号主要有两个:每秒1个脉冲的1PPS和每秒1帧数据的时间报文TOD。1PPS信号的上升沿反映了每一秒的开始准时刻,而TOD中的年月日时分秒等信息反映了该时刻对应的准确时间。本地配置一个时钟频率源,通常为恒温晶振或原子钟。在现有条件如何精确地测量本地时钟和卫星时钟之间的频差是一个需要解决的关键技术问题。目前文献中提出的基于模拟时间扩展计数法、基于AD变换器的模拟时间-幅度转换法、基于延迟线的时间-数字变换器(TDC)法和基于冲击振荡器的频率游标法等等在这里有的不适用,有的测量精度难以达到预期水平。易于实现的计数法计数分辨率受限于当前数字电路工作频率的极限,在实际应用中难以达到设计精度。本发明利用硬件和软件相结合,采用基于滑动窗口时间扩展测量方法,有效提高了测量精度,较好地解决了这个问题。
(2)解决高精度频率调整问题。目前工业上广泛采用的恒温晶振、铷原子钟大都采用压控技术来进行频率微调。它的输入除了电源,还有一个电压信号输入端子,用来控制微调输出频率。但是这种电压与频率的对应关系不是线性的,对于一个给定的测量出的频差值,换算成一个准确的微调电压值是困难的。因此经常是调过了头反过来再调,造成频率调节不准确、收敛速度慢、频率稳定度差等不良结果。本发明提出了精确频率再生方法,直接合成准确频率,有效的解决了频率调节不准确、收敛速度慢、频率稳定度差的问题。
为了达到上述发明目的,本发明采用的技术方案为:一种基于滑动窗口时间扩展高精度驯钟方法,整个本地时钟驯服子系统是由卫星接收模块、本地时钟、FPGA、MCU+Memory、数模转换DA+滤波组成,本地时钟提供参考时钟频率信号为f0,最终驯服输出的时钟频率信号为f1,设时间扩展窗口大小为N,其特征在于:
参考时钟频率信号f0为FPGA工作主时钟,以1PPS信号上升沿为基准,在FPGA中生成周期性的时隙计数器TSC,计数值从0累加到N-1,然后循环;
每个时隙都关联一个64位计数器,它们在每一个参考时钟频率信号f0的上升沿累加1;
时隙0计数器C64TS0在时隙0开始时开始计数,到下一时隙0结束计数,周而复始;
时隙1计数器C64TS1在时隙1开始计数,在下一个时隙1结束计数,周而复始;
其它时隙计数器工作以此类推;
MCU每秒取出相应时隙的计数值,存入到队列中,队列首元素去除,新数据从尾部加入,窗口下滑1个;对队列窗口中的所有数值利用最小二乘法获得中间均值;这个均值就是测量得到的频差数据;
在FPGA中设计数字DDS电路,根据MCU给出的频差值,生成准确的频率相位数据,这个数据送入数模转换DA,最终生成驯服频率f1。
有益效果:在现有的卫星导航定位系统、高精度授时设备中,普遍采用基于数字频差测量、压控频率微调的本地时钟驯服方法。这种方法的特点是在硬件和软件上易于实现,成本低。但是该方法同时具有频率驯服精度不高(最高仅能达到10-9Hz量级)、反复调整频率收敛速度慢、频率稳定度差的不足,制约着设备的定位和授时精度以及稳定度。
针对现有卫星星座原子钟与本地时钟频差数字测量精度不高的不足,本发明提出基于滑动窗口时间扩展测量方法,在现有硬件和软件技术条件下有效提高了频差测量精度,为提高本地时钟频率驯服精度提供了有力支撑。针现有方案对频率调整过程收敛速度慢、调整精度受限、频率稳定性不高的缺点,本发明提出了基于FPGA的数字DDS方法,直接合成准确频率,频率调节速度快,调整分辨率细精度高,使得整个频率调节过程从负反馈式的频率跟随到频率直接生成,实现了本地时钟驯服的快速、准确和有效。
附图说明
图1为现有技术中的导航/授时装置本地钟驯服方法示意图。
图2为本发明基于滑动窗口时间扩展高精度驯钟方法的系统框图。
图3为本发明实施例的 FPGA频差测量图。
图4为本发明实施例的在FPGA中实现数字DDS示意图。
图5为本发明实施例的MCU中程序流程图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。
本发明公开了北斗/GPS卫星导航定位和高精度授时系统中的一种基于滑动窗口时间扩展高精度驯钟方法。整个本地时钟驯服子系统是由卫星接收模块、本地时钟、FPGA、MCU+Memory、数模转换DA+滤波组成。本地时钟提供参考时钟频率f0,最终驯服输出的时钟频率为f1。设时间扩展窗口大小为N,参考时钟频率f0为FPGA工作时钟,以卫星模块输出的1PPS信号上升沿为基准,在FPGA中生成周期性的时隙计数器TSC,计数值从0累加到N-1,然后循环;每个时隙都关联一个64位计数器,它们在每一个f0的上升沿累加1,从本帧时隙开始计数,到下帧相应时隙开始结束计数;MCU每秒取出相应时隙的计数值并存入到队列中,队列首元素去除,新数据从尾部加入,窗口下滑;利用最小二乘法做数值滤波并获得中间均值,计算得到卫星模块与本地时钟频率差;在FPGA中设计数字DDS电路,根据MCU给出的频差值,生成准确的频率相位数据,这个数据送入数模转换DA,最终生成驯服频率f1。本发明时钟驯服精度相对于现有方法能够实现数量级式提升,而且驯服时间短,输出频率稳定,具有重要的应用价值。
本实施例的系统框图如图2所示。整个本地时钟驯服子系统是由卫星接收模块、本地时钟、FPGA、MCU+Memory、数模转换DA+滤波等几部分组成。本地时钟提供参考时钟频率信号为f0,最终驯服输出的时钟频率信号为f1。设时间扩展窗口大小为N:
1.参考f0的上升沿,以1PPS信号上升沿为基准,在FPGA中生成周期性的时隙计数器TSC,计数值从0累加到N-1,然后循环。
2.每个时隙都关联一个64位计数器,它们在每一个f0的上升沿累加1。
3.时隙0计数器C64TS0在时隙0开始时开始计数,到下一时隙0结束计数,周而复始。
4.时隙1计数器C64TS1在时隙1开始计数,在下一个时隙1结束计数,周而复始。其它时隙计数器工作类似。
5.MCU每秒取出相应时隙的计数值,存入到队列中,队列首元素去除,新数据从尾部加入,窗口下滑1个。对队列窗口中的所有数值利用最小二乘法获得中间均值。这个均值就是测量得到的频差数据。
6.在FPGA中设计数字DDS电路,根据MCU给出的频差值,生成准确的频率相位数据,这个数据送入DA,最终生成驯服频率f1。
本发明的技术特征有:
1.频差测量精度与窗口大小N相关。测量精度提高的量级约为logN。
2.频差测量速度为每秒1个数值。
3.在FPGA中生成数字DDS电路,根据MCU送入的频差数据,瞬间生成准确的输出频率,对本地时钟实现了快速、准确、稳定的驯服。
1.频差快速循环测量,如图3所示。
利用本地时钟f0,对卫星模块输出的1PPS信号进行计数测量,测量区间为N秒。TS0测量区间从本0时隙开头到下一个0时隙开头结束,TS1测量区间从本1时隙开头到下一个1时隙开头结束,以此类推,如图3所示。每秒都有一个新的测量值,测量速度为每秒1个数值。
2.频差计算(f0=100MHz)。
设测量计算得到的计数均值为P,则每秒相位差为(P/N-100000000),频率差为
3.FPGA中数字DDS实现,如图4所示,在FPGA中实现数字DDS示意图。
使用64位频率控制字FTW,在每个f0的上升沿,64位相位累加器ACC累加FTW,取ACC的高8位查正弦表,相应正弦数据输出到DA。
4.图5为本发明的MCU中程序流程图。
初始化变量和数据结构;
1PPS中断是否到来;
读取时隙寄存器;
读取相应时隙测量寄存器;
利用最小二乘法做数值滤波并获得中间均值;
计算频差值;
把频差数据写入FPGA中的数字DDS电路,生成驯服频率;
继续查询下一个1PPS信号,周而复始。
本实施例针对卫星导航和授时设备应用所公开的快速高精度本地时钟驯服方法,能够有效提高导航和授时设备的准确度和稳定度,具有重要的实用价值。
以上所述仅为本发明之较佳实施例而已,并非以此限制本发明的实施范围,凡熟悉此项技术者,运用本发明的原则及技术特征,所作的各种变更及装饰,若无显著技术改变,皆应涵盖于本权利要求书所界定的保护范畴之内。

Claims (3)

1.一种基于滑动窗口时间扩展高精度驯钟方法,整个本地时钟驯服子系统是由卫星接收模块、本地时钟、FPGA、MCU+Memory、数模转换DA+滤波组成,本地时钟提供参考时钟频率信号为f0,最终驯服输出的时钟频率信号为f1,设时间扩展窗口大小为N,其特征在于:
参考时钟频率f0为FPGA工作主时钟,以1PPS信号上升沿为基准,在FPGA中生成周期性的时隙计数器TSC,计数值从0累加到N-1,然后循环;
每个时隙都关联一个64位计数器,它们在f0的上升沿累加1;
时隙0计数器C64TS0在时隙0开始时开始计数,到下一时隙0结束计数,周而复始;
时隙1计数器C64TS1在时隙1开始计数,在下一个时隙1结束计数,周而复始;
其它时隙计数器工作以此类推;
MCU每秒取出相应时隙的计数值,存入到队列中,队列首元素去除,新数据从尾部加入,窗口下滑1个;对队列窗口中的所有数值利用最小二乘法获得中间均值,这个均值就是测量得到的频差数据;计算得到卫星接收模块与本地时钟频率差;
在FPGA中设计数字DDS电路,根据MCU给出的频差值,生成准确的频率相位数据,这个数据送入数模转换DA,经过滤波最终生成驯服频率f1
2.根据权利要求1所述的基于滑动窗口时间扩展高精度驯钟方法,其特征在于,所述频差数据计算:
设测量计算得到的计数均值为P,则每秒相位差为(P/N-f0),频率差为
频差测量速度快,每秒都能计算出一个高精度的频差数据,用来生成最终的驯服频率。
3.根据权利要求1所述的基于滑动窗口时间扩展高精度驯钟方法,其特征在于,频率调整在FPGA中由数字DDS电路实现:
使用64位频率控制字FTW,在每个参考时钟频率信号f0的上升沿,64位相位累加器ACC累加频率控制字FTW,取相位累加器ACC的高8位查正弦表,相应正弦数据输出到数模转换DA。
CN201610998692.2A 2016-11-14 2016-11-14 一种基于滑动窗口时间扩展高精度驯钟方法 Expired - Fee Related CN106383438B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610998692.2A CN106383438B (zh) 2016-11-14 2016-11-14 一种基于滑动窗口时间扩展高精度驯钟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610998692.2A CN106383438B (zh) 2016-11-14 2016-11-14 一种基于滑动窗口时间扩展高精度驯钟方法

Publications (2)

Publication Number Publication Date
CN106383438A true CN106383438A (zh) 2017-02-08
CN106383438B CN106383438B (zh) 2018-08-24

Family

ID=57958578

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610998692.2A Expired - Fee Related CN106383438B (zh) 2016-11-14 2016-11-14 一种基于滑动窗口时间扩展高精度驯钟方法

Country Status (1)

Country Link
CN (1) CN106383438B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147874A (zh) * 2020-11-05 2020-12-29 北京航天发射技术研究所 基于卫星授时和cpt原子钟守时的时频基准生成装置和方法
CN113985719A (zh) * 2021-10-25 2022-01-28 中国科学院国家授时中心 基于滑动窗口的脉冲星时驯服铯原子钟的方法
CN114137819A (zh) * 2021-12-06 2022-03-04 上海珉嵘科技有限公司 时钟频偏调整装置、方法及卫星信号采集预处理板卡
CN115657448A (zh) * 2022-12-09 2023-01-31 电信科学技术第五研究所有限公司 用于时钟驯服的频率偏移量分析方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1996036154A1 (en) * 1995-05-09 1996-11-14 Nokia Telecommunications Oy Data transmission system with sliding-window data flow control
CN101562451A (zh) * 2009-05-27 2009-10-21 西安华伟电力电子技术有限责任公司 二级频标的精密驯服保持方法
CN102147474A (zh) * 2010-12-21 2011-08-10 西安市双合软件技术有限公司 一种基于gps/北斗系统的时间频率驯服模块
CN102436174A (zh) * 2011-10-26 2012-05-02 东莞市泰斗微电子科技有限公司 一种守时设备晶振频率驯服方法及相应装置
CN203773267U (zh) * 2013-12-27 2014-08-13 中国科学院国家授时中心 可驯频率基准
CN104570717A (zh) * 2013-10-25 2015-04-29 沈阳工业大学 基于gps/北斗和有限状态机的守时系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1996036154A1 (en) * 1995-05-09 1996-11-14 Nokia Telecommunications Oy Data transmission system with sliding-window data flow control
CN101562451A (zh) * 2009-05-27 2009-10-21 西安华伟电力电子技术有限责任公司 二级频标的精密驯服保持方法
CN102147474A (zh) * 2010-12-21 2011-08-10 西安市双合软件技术有限公司 一种基于gps/北斗系统的时间频率驯服模块
CN102436174A (zh) * 2011-10-26 2012-05-02 东莞市泰斗微电子科技有限公司 一种守时设备晶振频率驯服方法及相应装置
CN104570717A (zh) * 2013-10-25 2015-04-29 沈阳工业大学 基于gps/北斗和有限状态机的守时系统
CN203773267U (zh) * 2013-12-27 2014-08-13 中国科学院国家授时中心 可驯频率基准

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147874A (zh) * 2020-11-05 2020-12-29 北京航天发射技术研究所 基于卫星授时和cpt原子钟守时的时频基准生成装置和方法
CN112147874B (zh) * 2020-11-05 2022-04-15 北京航天发射技术研究所 基于卫星授时和cpt原子钟守时的时频基准生成装置和方法
CN113985719A (zh) * 2021-10-25 2022-01-28 中国科学院国家授时中心 基于滑动窗口的脉冲星时驯服铯原子钟的方法
CN113985719B (zh) * 2021-10-25 2022-09-16 中国科学院国家授时中心 基于滑动窗口的脉冲星时驯服铯原子钟的方法
CN114137819A (zh) * 2021-12-06 2022-03-04 上海珉嵘科技有限公司 时钟频偏调整装置、方法及卫星信号采集预处理板卡
CN114137819B (zh) * 2021-12-06 2023-11-03 上海珉嵘科技有限公司 时钟频偏调整装置、方法及卫星信号采集预处理板卡
CN115657448A (zh) * 2022-12-09 2023-01-31 电信科学技术第五研究所有限公司 用于时钟驯服的频率偏移量分析方法
CN115657448B (zh) * 2022-12-09 2023-03-28 电信科学技术第五研究所有限公司 用于时钟驯服的频率偏移量分析方法

Also Published As

Publication number Publication date
CN106383438B (zh) 2018-08-24

Similar Documents

Publication Publication Date Title
CN102566410B (zh) 基于卫星授时的对本地时钟进行校准的方法和装置
CN101594128B (zh) 组合导航处理器用同步脉冲合成方法及同步脉冲合成器
CN106383438A (zh) 一种基于滑动窗口时间扩展高精度驯钟方法
CN202256483U (zh) 一种gps秒基实时自适应均匀间隔采样同步数据采集装置
CN106154299A (zh) 一种gps/sins组合导航系统时间同步方法
CN105049040A (zh) 一种利用gnss校正相干布居囚禁原子钟输出频率的方法
CN103269262B (zh) 一种时间同步装置的守时方法
CN104254761A (zh) 传感器时间同步
CN204631463U (zh) 一种辐射源监测定位的gps时钟同步系统
CN104410413A (zh) 原子频标频率修正方法、装置及原子频标
CN108008424A (zh) 一种卫星导航接收机秒脉冲的生成方法及装置
CN103605138A (zh) 一种卫星导航接收机时钟修正方法与装置
CN103941622A (zh) 基于fpga的高精度秒脉冲倍频出采样脉冲的方法
CN109765583A (zh) 一种基于gnss接收机秒脉冲的时钟同步方法
CN105629061A (zh) 一种基于高稳定度宽基准脉冲的精密频率测量装置
CN103699001A (zh) 利用恒温晶体振荡器实现低成本高精度的记时方法及系统
CN111766771A (zh) 一种基于压控晶振驯服的时间间隔测量方法及系统
CN103138749A (zh) 改进的频率校正方法和装置
CN115685725A (zh) 测量设备的时钟校准装置和测量设备
CN201540331U (zh) 多通道高精度同步测频装置
CN107976700A (zh) 一种卫星导航接收机秒脉冲的稳定输出方法
CN103684358A (zh) 采样脉冲生成方法和装置
CN110928177B (zh) 一种时钟同步系统及方法
CN108873669B (zh) 一种计算机同步时钟的utc时间计算方法
CN205720429U (zh) 一种基于高稳定度宽基准脉冲的精密频率测量装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180824

Termination date: 20211114