CN110865532B - 星地双向时间频率同步方法 - Google Patents
星地双向时间频率同步方法 Download PDFInfo
- Publication number
- CN110865532B CN110865532B CN201911165460.9A CN201911165460A CN110865532B CN 110865532 B CN110865532 B CN 110865532B CN 201911165460 A CN201911165460 A CN 201911165460A CN 110865532 B CN110865532 B CN 110865532B
- Authority
- CN
- China
- Prior art keywords
- satellite
- time
- ground station
- calculating
- frequency
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting 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
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting 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
- G04R20/06—Decoding time data; Circuits therefor
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
本发明公开了一种星地双向时间频率同步方法,包括:卫星的上行载波频率f1,下行载波频率为和f2和f3;计算f1、f2和f3的周数;计算电离层电子浓度TEC;消除f3伪码测量噪声;计算卫星和地面站之间的时差;计算卫星和地面站之间的频率差。相比于现有技术,本发明提出的技术方案基于宽带扩频理论与三频载波相位测量机制,克服了双向法中载波相位只能进行频率同步的问题,实现利用载波相位进行时间同步的目的,时差同步水平达到皮秒量级,可以满足当前氢原子钟和铯原子钟等高性能原子钟远程时间频率传递。
Description
技术领域
本发明涉及时间同步领域,特别是涉及一种卫星与地面站之间双向时间频率同步方法。
背景技术
双向法时间频率传递是目前国际公认的最高精度的时间频率传递方法,广泛应用于高精度时间频率量值远程比对,在时间频率量传和溯源方法中的地位无可替代。卫星与地面站之间双向的时间频率传递,利用信号扩频调制技术对原子钟时间频率信息进行调制形成上行信号。通过星地链路将两地面站的上行调制信号进行实时转发,两地面站同时对各自下行比对信号进行快速捕获、精密跟踪和精确解算,得到两个信号传播时延,可以精确的获得两地面站间的时间差信息。
传统双向时间频率同步方法主要包括两种,分别是基于码相位测量的时间频率同步方法和基于载波相位测量的频率同步方法。基于码相位测量的时间频率同步方法测量参数为码相位,受到码速率的限制,在典型值2.5MHz码速率情况下,时间比对精度仅为1ns,短期频率稳定度只能达到2×10-10/s。基于载波相位的频率同步方法虽然在频率同步方面有巨大提升,短期稳定度比基于码相位测量的方法高一个量级,但是由于存在载波相位模糊问题,只能进行频率同步,不能进行时间同步。
鉴于此,本发明的目的在于提供一种星地双向时间频率同步方法,基于宽带扩频理论与三频载波相位测量机制,以缓解现有技术中存在的同步精度受限问题。
发明内容
为了解决上述技术问题,本发明的目的在于提供一种星地双向时间频率同步方法,以缓解现有技术存在的问题。
一种星地双向时间频率同步方法,包括:卫星的上行载波频率f1,下行载波频率为和f2和f3;计算f1、f2和f3的周数;计算电离层电子浓度TEC;消除f3伪码测量噪声;计算卫星和地面站之间的时差;计算卫星和地面站之间的频率差。
进一步地,计算f1周数的方法,包括:地面站从t1时刻向卫星发射载波频率为f1测距信号,卫星在t2时刻接收到测距信号,获得载波相位观测量,计算f1上行整周模糊度的公式为N1为f1整周数,为卫星t2时刻接收的f1信号的载波相位,Δτs(τs(t2))为卫星t2时刻接收的f1信号的码相位,c为光速,TEC为电离层电子浓度。
进一步地,计算f2周数的方法,包括:卫星从t3时刻向地面站发射载波频率为f2测距信号,地面站在t4时刻接收到测距信号,获得载波相位观测量,计算f2上行整周模糊度的公式为N2为f2整周数,为地面站t4时刻接收的f2信号的载波相位,Δτg(τg(t4))为地面站t4时刻接收的f2信号的码相位,c为光速,TEC为电离层电子浓度。
进一步地,计算f3周数的方法,包括:卫星从t5时刻向地面站发射载波频率为f3测距信号,地面站在t6时刻接收到测距信号,获得载波相位观测量,计算f3上行整周模糊度的公式为N3为f3整周数,为地面站t6时刻接收的f3信号的载波相位,Δτg(τg(t6))为地面站t6时刻接收的f3信号的码相位,c为光速,TEC为电离层电子浓度。
进一步地,消除f3伪码测量噪声的方法,包括:对N3和TEC的计算公式修正为
进一步地,计算卫星和地面站之间时差的方法,包括:时差的计算公式为
进一步地,计算卫星和地面站之间频率差的方法,包括:对时差数据进行Allan方差计算。
进一步地,f3选择比f1和f2码速率低的伪码扩频信号。
本发明的有益效果如下:
本发明提供的技术方案可以包括以下有益效果:本发明提出的技术方案基于宽带扩频理论与三频载波相位测量机制,克服了双向法中载波相位只能进行频率同步的问题,实现利用载波相位进行时间同步的目的,时差同步水平达到皮秒量级,可以满足当前氢原子钟和铯原子钟等高性能原子钟远程时间频率传递。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本公开。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一种实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例的一种星地双向时间频率同步方法流程示意图;
图2是本发明实施例的一种星地双向时间频率同步方法各信号传送关系示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,所描述的实施例是本发明一部分实施例,而不是全部的实施例。
图1是本发明实施例的一种星地双向时间频率同步方法流程示意图,如图1所示,该方法包括如下五个步骤。
步骤S101:计算f1、f2和f3的周数。如图2所示,为了计算f1、f2和f3的周数,地面站与卫星之间进行三个频率为f1、f2和f3的信号在两者之间传递。需要进行说明的是,f3选择比f1和f2码速率低的伪码扩频信号。上行f1载波和下行f2载波主要用于星地双向高精度时差测量,下行f3载波主要用于与下行f2载波进行双频实时消除电离层延迟。f3频点的码相位观测量可通过f1和f2频点联合解算消除,因此f3频点可选择码速率较低伪码扩频信号。
在一个可选的实施例中,f1载波频率为14GHz,其上调制扩频码信号带宽为200MHz;采用f2载波频率为12GHz,其上调制扩频码信号带宽为200MHz;采用f3载波频率为4GHz,其上调制扩频码信号带宽为2MHz。
详细地,计算f1周数的方法,包括:地面站从t1时刻向卫星发射载波频率为f1测距信号,卫星在t2时刻接收到测距信号,获得载波相位观测量,计算f1上行整周模糊度的公式为N1为f1整周数,为卫星t2时刻接收的f1信号的载波相位,Δτs(τs(t2))为卫星t2时刻接收的f1信号的码相位,c为光速,TEC为电离层电子浓度。
计算f2周数的方法,包括:卫星从t3时刻向地面站发射载波频率为f2测距信号,地面站在t4时刻接收到测距信号,获得载波相位观测量,计算f2上行整周模糊度的公式为N2为f2整周数,为地面站t4时刻接收的f2信号的载波相位,Δτg(τg(t4))为地面站t4时刻接收的f2信号的码相位,c为光速,TEC为电离层电子浓度。
计算f3周数的方法,包括:卫星从t5时刻向地面站发射载波频率为f3测距信号,地面站在t6时刻接收到测距信号,获得载波相位观测量,计算f3上行整周模糊度的公式为N3为f3整周数,为地面站t6时刻接收的f3信号的载波相位,Δτg(τg(t6))为地面站t6时刻接收的f3信号的码相位,c为光速,TEC为电离层电子浓度。
步骤S103:消除f3伪码测量噪声。具体地,消除f3伪码测量噪声的方法,包括:
对N3和TEC的计算公式修正为
步骤S104:计算卫星和地面站之间的时差。具体地,时差的计算公式为
步骤S105:计算卫星和地面站之间的频率差。具体地,对时差数据进行Allan方差计算。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (2)
1.一种星地双向时间频率同步方法,其特征在于,包括:
卫星的上行载波频率f1,下行载波频率为和f2和f3;
计算所述f1、f2和f3的周数;
计算电离层电子浓度TEC;
消除所述f3伪码测量噪声;
计算所述卫星和地面站之间的时差;
计算所述卫星和地面站之间的频率差;
其中,计算所述f1周数的方法,包括:
所述地面站从t1时刻向卫星发射载波频率为f1测距信号,所述卫星在t2时刻接收到所述测距信号,获得载波相位观测量,计算所述f1上行整周模糊度的公式为所述N1为f1整周数,为所述卫星t2时刻接收的f1信号的载波相位,Δτs(τs(t2))为所述卫星t2时刻接收的f1信号的码相位,c为光速,TEC为电离层电子浓度;
计算所述f2周数的方法,包括:
所述卫星从t3时刻向地面站发射载波频率为f2测距信号,所述地面站在t4时刻接收到所述测距信号,获得载波相位观测量,计算所述f2上行整周模糊度的公式为所述N2为f2整周数,为所述地面站t4时刻接收的f2信号的载波相位,Δτg(τg(t4))为所述地面站t4时刻接收的f2信号的码相位,c为光速,TEC为电离层电子浓度;
计算所述f3周数的方法,包括:
所述卫星从t5时刻向地面站发射载波频率为f3测距信号,所述地面站在t6时刻接收到所述测距信号,获得载波相位观测量,计算所述f3上行整周模糊度的公式为所述N3为f3整周数,为所述地面站t6时刻接收的f3信号的载波相位,Δτg(τg(t6))为所述地面站t6时刻接收的f3信号的码相位,c为光速,TEC为电离层电子浓度;
所述计算电离层电子浓度TEC的方法,包括:所述电离层电子浓度TEC的计算公式为
所述消除f3伪码测量噪声的方法,包括:
对所述N3和TEC的计算公式修正为
所述计算卫星和地面站之间时差的方法,包括:所述时差的计算公式为
所述计算卫星和地面站之间频率差的方法,包括:对所述时差数据进行Allan方差计算。
2.根据权利要求1所述的方法,其特征在于,所述f3选择比f1和f2码速率低的伪码扩频信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911165460.9A CN110865532B (zh) | 2019-11-25 | 2019-11-25 | 星地双向时间频率同步方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911165460.9A CN110865532B (zh) | 2019-11-25 | 2019-11-25 | 星地双向时间频率同步方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110865532A CN110865532A (zh) | 2020-03-06 |
CN110865532B true CN110865532B (zh) | 2021-04-02 |
Family
ID=69656245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911165460.9A Active CN110865532B (zh) | 2019-11-25 | 2019-11-25 | 星地双向时间频率同步方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110865532B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113098582B (zh) * | 2021-03-25 | 2022-06-21 | 中国科学院国家授时中心 | 一种卫星双向时间频率传递方法 |
CN113489529B (zh) * | 2021-07-26 | 2023-01-20 | 中电防务科技有限公司 | 低轨星座卫星通信***中站间同步方法 |
US11864140B2 (en) * | 2022-01-19 | 2024-01-02 | Intelligent Fusion Technology, Inc. | Methods and systems for time synchronization among unmanned aerial systems |
CN115776312B (zh) * | 2023-02-13 | 2023-05-05 | 中国人民解放军国防科技大学 | 一种面向星地高精度微波时频传递的信号体制设计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104407326A (zh) * | 2014-10-31 | 2015-03-11 | 中国运载火箭技术研究院 | 一种适用于集群目标间的射频自主相对姿态测量*** |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101251594B (zh) * | 2008-04-03 | 2011-04-06 | 北京航空航天大学 | 双向测距与时间比对处理终端 |
CN101702030B (zh) * | 2009-10-29 | 2012-10-24 | 中国科学院国家授时中心 | 一种站间钟差辅助的导航卫星精密定轨方法 |
CN102545993B (zh) * | 2011-12-20 | 2014-04-02 | 中国科学院国家授时中心 | 一种基于载波相位的卫星双向时间传递方法 |
CA2887071A1 (en) * | 2012-10-04 | 2014-04-10 | Thomas Ronald Tessier | A method of routing of data messages from mobile devices through satellite and terrestrial communication networks |
CN104298860B (zh) * | 2014-09-22 | 2017-04-12 | 中国科学院国家授时中心 | 一种geo卫星转发信号码和载波相位偏差的计算方法 |
CN108134625B (zh) * | 2017-12-11 | 2020-09-01 | 北京无线电计量测试研究所 | 一种卫星双向时间频率传递方法 |
CN108845338B (zh) * | 2018-07-26 | 2022-04-01 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 星地联合载波平滑伪距测距方法 |
CN109507690A (zh) * | 2018-11-09 | 2019-03-22 | 中国科学院国家授时中心 | 基于gnss的国家标准时间亚纳秒级授时方法 |
-
2019
- 2019-11-25 CN CN201911165460.9A patent/CN110865532B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104407326A (zh) * | 2014-10-31 | 2015-03-11 | 中国运载火箭技术研究院 | 一种适用于集群目标间的射频自主相对姿态测量*** |
Also Published As
Publication number | Publication date |
---|---|
CN110865532A (zh) | 2020-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110865532B (zh) | 星地双向时间频率同步方法 | |
Levine | A review of time and frequency transfer methods | |
Fujieda et al. | Carrier-phase two-way satellite frequency transfer over a very long baseline | |
US11221418B2 (en) | Method for correcting a pseudorange in a receiver for satellite navigation | |
US20190271758A1 (en) | Radio receiver for determining location of a signal source | |
US8456353B2 (en) | Method and system for determining clock corrections | |
CN108134625B (zh) | 一种卫星双向时间频率传递方法 | |
US20150029058A1 (en) | Accurate timing distribution by high-frequency radio | |
US20170299692A1 (en) | Interference mitigation for positioning systems | |
CN103516457A (zh) | 一种高精度远程时间同步方法 | |
CN103576166A (zh) | 基于单颗静止轨道卫星对干扰源进行定位的方法及*** | |
CN113777598A (zh) | 深空x频段测距和干涉一体化信标装置、测量方法和*** | |
CN105425262B (zh) | 一种实现卫星转发导航***载波相位精密测量的方法 | |
US20160179748A1 (en) | Method and apparatus for estimating waveform onset time | |
CN116299603B (zh) | 应用于地面站的导航上行伪距修正量的获取方法和装置 | |
CN116961786A (zh) | 一种长波发射机的发射通道时延自校准装置及方法 | |
CN116626720A (zh) | 一种低轨卫星的终端定位方法、设备和介质 | |
Cannon | Mitigation and exploitation of the ionosphere: A military perspective | |
CN105959092B (zh) | 基于直发与转发信号的区域高精度授时与校频方法 | |
RU166018U1 (ru) | Устройство для синхронизации часов | |
Moon et al. | Position determination using the DTV segment sync signal | |
Vilnrotter et al. | Carrier-phase aided pseudo-noise range estimation at rf frequencies | |
Motella et al. | An interference impact assessment model for GNSS signals | |
Yang et al. | Improved two-way satellite time and frequency transfer with Multi-GEO in BeiDou navigation system | |
Christensen et al. | Solution for wireless time synchronization using sub-Nyquist sampling rates |
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 |