CN114995098B - 基于相时延测量的皮秒级精度的卫星干涉时间传递方法 - Google Patents
基于相时延测量的皮秒级精度的卫星干涉时间传递方法 Download PDFInfo
- Publication number
- CN114995098B CN114995098B CN202210427908.5A CN202210427908A CN114995098B CN 114995098 B CN114995098 B CN 114995098B CN 202210427908 A CN202210427908 A CN 202210427908A CN 114995098 B CN114995098 B CN 114995098B
- Authority
- CN
- China
- Prior art keywords
- delay
- stations
- time
- satellite
- interference
- 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
- 238000012546 transfer Methods 0.000 title claims abstract description 44
- 238000005259 measurement Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000012937 correction Methods 0.000 claims abstract description 29
- 230000005540 biological transmission Effects 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 12
- 238000012805 post-processing Methods 0.000 claims abstract description 11
- 238000005314 correlation function Methods 0.000 claims description 13
- 238000005305 interferometry Methods 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 8
- 230000010355 oscillation Effects 0.000 claims description 4
- 101100167360 Drosophila melanogaster chb gene Proteins 0.000 description 20
- 238000005516 engineering process Methods 0.000 description 9
- 239000005433 ionosphere Substances 0.000 description 8
- 238000012360 testing method Methods 0.000 description 7
- 239000005436 troposphere Substances 0.000 description 6
- 239000013307 optical fiber Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- MRMBZHPJVKCOMA-YJFSRANCSA-N biapenem Chemical compound C1N2C=NC=[N+]2CC1SC([C@@H]1C)=C(C([O-])=O)N2[C@H]1[C@@H]([C@H](O)C)C2=O MRMBZHPJVKCOMA-YJFSRANCSA-N 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/04—Tuning or receiving; 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
本发明提供的一种基于相时延测量的皮秒级精度的卫星干涉时间传递方法,通过精密单点定位(PPP)时间传递得到钟差初值;根据卫星坐标,测站坐标、各项修正量和钟差初值,计算理论时延;对观测数据依次进行相关处理和相关后处理,得到干涉群时延结果以及干涉相时延结果;根据干涉群时延结果,解算出模糊度;根据理论时延、干涉相时延结果以及解算出的模糊度,计算两个测站间的最终钟差。本发明中地面测站不主动发射信号,被动接收轨道已知的卫星的任意下行信号,无需解调信号也无需在卫星上搭载原子钟载荷,通过群时延辅助进行干涉相时延测量,利用已知的轨道信息和测站坐标解算站间钟差实现时间传递,精度可达亚皮秒级。
Description
技术领域
本发明属于卫星授时技术领域,具体涉及一种基于相时延测量的皮秒级精度的卫星干涉时间传递方法。
背景技术
目前常用的时间传递技术主要可分为星基系统和陆基系统两大类。星基系统使用卫星开展时间传递,主要使用的有国际导航卫星系统(GNSS)、区域导航卫星系统(RNSS),其中我国的系统包括有北斗导航卫星系统 (BDS)和中国区域导航试验系统(CAPS)。陆基授时系统又可分为无线电授时和有线授时两种方式,其中无线电授时系统主要包括BPL长波授时系统、BPM短波授时系统以及BPC低频时码授时系统等,有线授时系统主要有NTP互联网授时系统、PSTN公网电话授时系统等,另外新研发的基于光纤链路开展的光纤时间传递也在不断扩展应用。各授时技术的授时精度如下表1所示。
表1各类授时技术的精度对比
对于长基线(百公里-千公里)应用,卫星共视和卫星授时可提供纳秒级别的时间传递精度,这也是国际权度局(BIPM)组织国际时间传递所使用的主要方法。
目前亚纳秒级的时间传递技术主要有基于GNSS的载波共视技术、精密单点定位(PPP)时间传递技术和基于光纤的时间传递技术。基于GNSS 的PPP时间传递和载波共视技术精度可达亚纳秒级别,无法满足皮秒、亚皮秒级别的高精度的时间传递需求。光纤时间传递精度高,但作为一种路基有线技术,必须依赖铺设的光缆,以光缆作为传输介质,应用范围受到局限。目前尚缺少一种适用于短基线(十公里级以内)的亚皮秒级无线时间传递方法。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于相时延测量的皮秒级精度的卫星干涉时间传递方法。本发明要解决的技术问题通过以下技术方案实现:
第一方面,本发明提供的一种基于相时延测量的皮秒级精度的卫星干涉时间传递方法包括:
获取两个测站共同观测同一卫星的观测数据;
根据GNSS数据,使用PPP解算出钟差初值;
根据卫星坐标,测站坐标、各项修正量和钟差初值,计算理论时延;
对观测数据依次进行相关处理和相关后处理,得到干涉群时延结果以及干涉相时延结果;
根据干涉群时延结果,解算出模糊度;
根据解算出的模糊度、干涉相时延结果以及理论时延,计算两个测站间的最终钟差,实现时间传递。
可选的,理论时延表示为:
其中:Δt′AB指AB两站的钟差初值,是通过PPP时间传递获得的较低精度的纳秒级钟差产品;c是真空中的光速,dAS、dBS指卫星S分别到A、 B站传播距离,指A、B站观测卫星S的对流层时延,/>指A、 B站观测卫星S的电离层时延,/>指A、B站观测卫星S的引力时延,/>指A、B站的系统硬件时延,oth指其它改正。
可选的,对观测数据进行相关处理包括:
对两个测站的观测数据进行理论时延补偿以及相位补偿;
对理论时延补偿和相位补偿后的观测数据进行傅里叶变换;
将傅里叶变换后的两测站的观测数据进行共轭相乘,得到互相关函数;
对互相关函数进行转化,得到关于频率的相位函数。
可选的,对观测数据进行相关后处理包括:
根据相位函数,计算干涉群时延结果以及干涉相时延结果。
可选的,理论时延补偿以及相位补偿后的两个测站的观测数据表示为:
wA(t)=vA(t+δ)e-iθ
wB(t)=vB(t)
θ=2πfLOδ
其中,wA(t)、wB(t)分别表示A、B两个测站补偿后的信号,fLO为变频器的本振频率;
傅里叶变换后的信号表示为:
互相关函数表示为:
相位函数表示为:
可选的,根据解算出的模糊度、干涉相时延结果以及理论时延,计算两个测站间的最终钟差,实现时间传递包括:
根据解算出的模糊度、干涉相时延结果以及理论时延,计算出干涉测量的时差值;
根据干涉测量的时差值,引入观测方程,计算测站间的最终钟差,实现时间传递。
其中,干涉群时延结果表示为:
干涉相时延结果表示为:
模糊度表示为:
干涉测量的时差值表示为:
Δt=δ+τphase;
其中,Δt指干涉测量的时差值。
其中,观测方程表示为:
其中:ΔtAB指AB两站的最终钟差,是用于时间传递的修正量。
第二方面,本发明提供的一种基于相时延测量的皮秒级精度的卫星干涉时间传递系统包括:两个测站,两个测站中包括一个主站;主站,用于:
获取两个测站共同观测同一卫星的观测数据;
根据GNSS数据,使用PPP解算出钟差初值;
将卫星坐标,测站坐标、各项修正量和钟差初值,计算理论时延;
对观测数据依次进行相关处理和相关后处理,得到干涉群时延结果以及干涉相时延结果;
根据干涉群时延结果,解算出模糊度;
根据解算出的模糊度、干涉相时延结果以及理论时延,计算两个测站间的最终钟差,实现时间传递。
其中,理论时延表示为:
其中:Δt′AB指AB两站的钟差初值,是通过PPP时间传递获得的较低精度的纳秒级钟差产品;c是真空中的光速,dAS、dBS指卫星S分别到A、 B站传播距离,指A、B站观测卫星S的对流层时延,/>指A、B站观测卫星S的电离层时延,/>指A、B站观测卫星S的引力时延,/>指A、B站的系统硬件时延,oth指其它改正。
本发明提供的一种基于相时延测量的皮秒级精度的卫星干涉时间传递方法,通过获取两个测站共同观测同一卫星的观测数据;基于GNSS数据解算出钟差初值;根据卫星坐标,测站坐标、各项修正量和钟差初值,计算理论时延;对观测数据依次进行相关处理和相关后处理,得到干涉群时延结果以及干涉相时延结果;根据干涉群时延结果,解算出模糊度;根据理论时延、干涉相时延结果以及解算出的模糊度,计算两个测站间的钟差,实现时间传递。本发明中地面测站不主动发射信号,被动接收轨道已知的卫星的任意下行信号,无需解调信号也无需在卫星上搭载原子钟载荷,通过群时延辅助进行干涉相时延测量,利用已知的轨道信息和测站坐标解算站间钟差,实现时间传递,精度可达亚皮秒级。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明实施例提供的一种基于相时延测量的皮秒级精度的卫星干涉时间传递方法的流程示意图;
图2是本发明实施例提供的基于相时延测量的皮秒级精度的卫星干涉时间传递原理图;
图3是本发明实施例提供的干涉相时延时间传递结果图;
图4a是本发明实施例提供的干涉相时延时间传递结果的频率稳定度示意图;
图4b是本发明实施例提供的干涉相时延时间传递结果的时间稳定度示意图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。
如图1所示,本发明提供的一种基于相时延测量的皮秒级精度的卫星干涉时间传递方法包括:
S1,获取两个测站共同观测同一卫星的观测数据;
S2,根据GNSS数据,使用PPP解算出钟差初值;
本步骤可以根据两个测站的观测数据,使用PPP时间传递解算获得钟差初值。
传播时延表示为:
S3,根据卫星坐标,测站坐标、各项修正量和钟差初值,计算理论时延;
理论时延表示为:
其中:Δt′AB指AB两站的钟差初值,是通过PPP时间传递获得的较低精度的纳秒级钟差产品;c是真空中的光速,dAS、dBS指卫星S分别到A、 B站传播距离,指A、B站观测卫星S的对流层时延,/>指A、 B站观测卫星S的电离层时延,/>指A、B站观测卫星S的引力时延,/>指A、B站的系统硬件时延,oth指其它改正。
参考图2,测站A、B同时观测卫星S,由于测站坐标和卫星轨道已知,可以计算信号从卫星到达测定轨主站的理论传播时延。
值得说明的是:引力时延、质心修正、相位中心修正等均可精确修正,这里不做详细描述。其他影响较大的误差项分析如下。
(1)传播路径误差
卫星到测站的传播几何距离可通过卫星和站的坐标计算,以站A为例,令A站坐标为(xA,yA,zA),卫星S坐标为(xS,yS,zS):
其中,站坐标和卫星轨道均可能存在误差。
(2)站坐标误差
在静态接收机授时过程中,一般通过已知的精确坐标代入求解计算。实际上,站坐标会受到固体引力潮汐、大陆板块迁移、极移、海潮等各种因素的影响。陆板块迁移的年变化为厘米级,尚不需考虑。可是,地球固体潮汐和海洋潮汐较大,一般通过经验模型精确改正。
(3)轨道误差
本发明以地球同步轨道(GEO)卫星为例。目前可提供精密轨道的GEO 卫星主要为北斗的GEO卫星,轨道精度可达米级。国家授时中心的转发测定轨技术可为任意通信GEO卫星提供米级轨道,其观测的卫星也可作为本专利观测的卫星。
由观测方程可知,卫星与两测站的距离差经过一次差分,可部分消除卫星轨道引入的误差。卫星轨道引入的误差可由下式计算。
(eB-eA)·ΔS/c(4)
其中eA、eB分别是A站、B站指向卫星的单位向量,ΔS为轨道误差。由于GEO卫星轨道高度约为三万六千公里,卫星轨道高度远大于基线长度时,为便于计算可粗略认为轨道引起的误差为:
(dAB/轨道高度)·ΔS/c(5)
其中,dAB为A、B之间的距离。
要求100公里的基线轨道引入的误差优于1ps,则粗略估计GEO卫星轨道精度应优于0.1米,MEO卫星轨道精度应优于6.6厘米,1000公里LEO 卫星轨道精度应优于3毫米,可暂不考虑。
(4)对流层时延
对流层模型可通过模型修正,以国际GNSS服务(IGS)使用的GMF 模型为例。
测站观测卫星的对流层时延与卫星高度角有关。
其中τtro为对流层延迟,即公式1中的或/>ES为卫星高度角,/>表示干燥大气引起的对流层偏差中的干分量;/>表示水汽分子的偶极距所引起的对流层偏差的湿分量;映射函数也相应的分为干、湿映射函数,记作Mdry(ES)和Mwet(ES)。
其中M(ES)指Mdry(ES)或Mwet(ES),a、b、c是远小于1的常数,且分为干气部分和湿气部分,分别对应于干分量投影函数和湿分量投影函数。
当基线较短时,两站观测卫星的路径基本一致,卫星高度角十分接近,使得两站观测卫星的对流层时延有很大的相关性,可通过差分消除部分影响。若基线南北分布,则其观测卫星的高度角差异最大。
最差情况下令测站与卫星位于同一经度,易知南北分布的100公里基线的纬度差为0.9°,观测GEO卫星的高度角之差在赤道处达到最大约1.06°,由于PPP解算的对流层产品精度优于1cm,则高度角大于41°时经过差分后,对流层误差的残余影响优于1ps。若测站为东西分布,则易知其观测GEO 卫星的高度角之差在赤道处最大为0.017°,高度角大于17°时经过差分后,对流层误差的残余影响优于0.1ps。
(5)电离层时延
电离层模型可通过双频观测改正,也可通过模型改正。下文以MSLM (modifiedsingle-layer model)模型为例。
电离层延迟可通过下式计算,
STEC=MF(z)·VTEC(9)
MF(z)=1/cos(arcsin(R/(R+H)·cos(a(90°-z))))(10)
其中τion为电离层延迟,即公式1中的或/>freq为观测频点, STEC指倾斜电子浓度,VTEC指垂直电子浓度,MF为投影函数,z为卫星高度角,R为地球半径,H为电离层高度,a为常数0.9782。
电离层延迟可通过双频观测修正,在没有双频观测的时候可采用使用 IGS或国际GNSS检测评估系统(iGMAS)提供的电离层网格产品修正。电离层网格产品的精度为2-8TECU(总电子数单位),假设单频信号为C 波段(4-8GHz),电离层延迟比在GPS L1状况下缩小6至25倍,单站电离层延迟的修正精度为0.17ns至0.67ns。
与对流层类似,短基线情况下可通过站间差分消除部分电离层影响。测站南北分布时,100公里基线下电离层的映射函数差异小于0.04,则电离层的修正的残余影响8ps至34ps。测站东西分布时,100公里基线下电离层的映射函数差异小于0.0008,则电离层的修正的残余影响0.1ps至0.5ps。测站南北分布时,建议通过双频观测的方法,利用电离层时延的平方反比特性消除电离层影响。
(6)系统时延
该系统存在设备引起的硬件时延,可通过移动站进行连线干涉的方式标定。即使用一个移动站C,在A站使用其时频信号开展相时延干涉时差测量;再将该移动站C移至B站,使用其时频信号开展相时延干涉时差测量;两次测量结果互差后可消除移动站C引入的误差,得到A-B基线的系统硬件时延,具体解算过程不再详述。
令两测站A、B接收到的信号为vA(t)、vB(t),且均为同一信号源下发的信号,不妨令B站为基准,则存在时差Δt,令vA(t+Δt)=vB(t)。相关处理即是求解两测站钟差Δt的过程。由于Δt较大时相关处理难以获得准确的结果,所以要先用理论时延δ进行修正,使vA(t+δ)与vB(t)较为接近,则可转为测量较小的残余时延τ=Δt-δ。
相关处理首先需要计算理论时延进行补偿,补偿的理论时延δ是根据卫星坐标、各项修正量和精度较差的钟差初值,即初步钟差计算的。
S4,对观测数据依次进行相关处理和相关后处理,得到干涉群时延结果以及干涉相时延结果;
本步骤可以对两个测站的观测数据进行理论时延补偿以及相位补偿;对理论时延补偿和相位补偿后的观测数据进行傅里叶变换;将傅里叶变换后的两测站的观测数据进行共轭相乘,得到互相关函数;对互相关函数进行转化,得到关于频率的相位函数;根据相位函数,计算干涉群时延结果以及干涉相时延结果。
值得说明的是:对观测数据进行理论时延补偿和相位补偿(条纹旋转),令补偿后的信号为wA(t)、wB(t),其中fLO为变频器的本振频率。wA(t)、wB(t) 分别表示A、B两个测站补偿后的信号;
wA(t)=vA(t+δ)e-iθ(11)
wB(t)=vB(t)(12)
θ=2πfLOδ(13)
对补偿后的信号进傅里叶变换后的信号表示为:
互相关函数为:
据互相关函数获得相位值即相位函数,/>是关于频率的函数。
相位函数表示为:
则可以获得残余时延τ的群时延结果τgroup和相时延结果τphase。
S5,根据干涉群时延结果,解算出模糊度;
其中,N为模糊度,在一定时间内为一个常整数,可根据群时延计算。
S6,根据解算出的模糊度、干涉相时延结果以及理论时延,计算两个测站间的钟差,实现时间传递。
本发明可以根据解算出的模糊度、干涉相时延结果以及理论时延,计算出干涉测量的时差值;根据干涉测量的时差值,引入观测方程,计算测站间的最终钟差,实现时间传递。
其中,观测方程表示为:
其中:ΔtAB指AB两站的最终钟差,是用于时间传递的修正量,Δt指干涉测量的时差值,是通过相关处理获得的测量值,c是真空中的光速,dAs指卫星S分别到A、B站传播距离,指A、B站观测卫星S的对流层时延,/>指A、B站观测卫星S的电离层时延,/>指A、 B站观测卫星S的引力时延,/>指A、B站的系统硬件时延,oth指其它改正,包括天线相位中心改正等。
则两站的时差值Δt则可根据相时延和理论时延获得,两个测站间的时差值表示为:
Δt=δ+τphase(22)
则两站的最终钟差值ΔtAB则可根据观测方程获得,两个测站间的钟差表示为:
值得说明的是:两个测站通过抛物面天线观测同一颗卫星,测站基于本地时间采集卫星信号,通过相关处理可获得卫星信号到达两个测站天线的时间差,通过卫星轨道和测站坐标可计算卫星到测站的实际几何距离差,二者差分后可获得两个测站的钟差,从而实现时间传递。
本发明提供的一种基于相时延测量的皮秒级精度的卫星干涉时间传递系统包括:两个测站,两个测站中包括一个主站;主站,用于:
获取两个测站共同观测同一卫星的观测数据;
根据GNSS数据,使用PPP解算出钟差初值;
将卫星坐标,测站坐标、各项修正量和钟差初值,计算理论时延;
对观测数据依次进行相关处理和相关后处理,得到干涉群时延结果以及干涉相时延结果;
根据干涉群时延结果,解算出模糊度;
根据解算出的模糊度、干涉相时延结果以及理论时延,计算两个测站间的最终钟差,实现时间传递。
下面结合实验验证本发明的最终效果。
为验证本方法,2021年11月17日-18日组织了验证试验。由于该方法现处于预研阶段,试验系统中使用了同一个原子钟,因此未使用PPP时间传递结果作为初值。
试验使用了两面3.7米口径的抛物面天线,基线长度为19.93米;观测卫星为北斗C01星,星下点位于144.5°E;观测频段为C波段;使用DiFX 相关处理软件系统和AIPS相关后处理软件进行相关处理和相关后处理。
试验中的轨道数据为iGMAS提供的事后精密轨道。试验没有进行大气延迟修正和设备时延修正。试验结果如下图3所示,存在的起伏可能是设备和大气时延引入的,在实际使用中可进一步修正。
时差结果的修正阿伦方差和时间方差如图4a和图4b所示,频率稳定度达到了7.35e-13(1s),时间稳定度达到了4.25e-13秒(1s),达到了亚皮秒级水平。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (2)
1.一种基于相时延测量的皮秒级精度的卫星干涉时间传递方法,其特征在于,包括:
获取两个测站共同观测同一卫星的观测数据;
根据GNSS数据,使用PPP解算出钟差初值;
根据卫星坐标、测站坐标、各项修正量和所述钟差初值,计算理论时延;所述理论时延表示为:
其中:Δt′AB指AB两站的钟差初值,是通过PPP时间传递获得的较低精度的纳秒级钟差产品;c是真空中的光速,dAS、dBS指卫星S分别到A、B站传播距离,指A、B站观测卫星S的对流层时延,/>指A、B站观测卫星S的电离层时延,/>指A、B站的系统硬件时延,oth指其它改正;
对所述观测数据依次进行相关处理和相关后处理,得到干涉群时延结果以及干涉相时延结果;
对所述观测数据进行相关处理包括:
对两个测站的观测数据进行理论时延补偿以及相位补偿;理论时延补偿以及相位补偿后的两个测站的观测数据表示为:
wA(t)=vA(t+δ)e-iθ;
wB(t)=vB(t);
θ=2πfLOδ;
其中,wA(t)、wB(t)分别表示A、B两个测站补偿后的信号,fLO为变频器的本振频率;vA(t)、vB(t)表示A、B两个测站接收到的信号;
对理论时延补偿和相位补偿后的观测数据进行傅里叶变换;傅里叶变换后的信号表示为:
将傅里叶变换后的两测站的观测数据进行共轭相乘,得到互相关函数;互相关函数表示为:
相位函数表示为:
对所述互相关函数进行转化,得到关于频率的相位函数;
对观测数据进行相关后处理包括:
根据相位函数,计算干涉群时延结果以及干涉相时延结果;干涉群时延结果表示为:
干涉相时延结果表示为:
根据干涉群时延结果,解算出模糊度;模糊度表示为:
根据解算出的模糊度、所述干涉相时延结果以及所述理论时延,计算两个测站间的最终钟差,实现时间传递;
所述根据解算出的模糊度、所述干涉相时延结果以及所述理论时延,计算两个测站间的最终钟差,实现时间传递包括:
根据解算出的模糊度、所述干涉相时延结果以及所述理论时延,计算出干涉测量的时差值;干涉测量的时差值表示为:
Δt=δ+τphase;
其中,Δt指干涉测量的时差值;
根据干涉测量的时差值,引入观测方程,计算测站间的最终钟差,实现时间传递,观测方程表示为:
其中:ΔtAB指AB两站的最终钟差,是用于时间传递的修正量。
2.一种基于相时延测量的皮秒级精度的卫星干涉时间传递系统,其特征在于,包括:两个测站,两个测站中包括一个主站;所述主站,用于:
获取两个测站共同观测同一卫星的观测数据;
根据GNSS数据,使用PPP解算出钟差初值;
将卫星坐标、测站坐标、各项修正量和所述钟差初值,计算理论时延;所述理论时延表示为:
其中:Δt′AB指AB两站的钟差初值,是通过PPP时间传递获得的较低精度的纳秒级钟差产品;c是真空中的光速,dAS、dBS指卫星S分别到A、B站传播距离,指A、B站观测卫星S的对流层时延,/>指A、B站观测卫星S的电离层时延,/>指A、B站的系统硬件时延,oth指其它改正;
对所述观测数据依次进行相关处理和相关后处理,得到干涉群时延结果以及干涉相时延结果;
对所述观测数据进行相关处理包括:
对两个测站的观测数据进行理论时延补偿以及相位补偿;理论时延补偿以及相位补偿后的两个测站的观测数据表示为:
wA(t)=vA(t+δ)e-iθ;
wB(t)=vB(t);
θ=2πfLOδ;
其中,wA(t)、wB(t)分别表示A、B两个测站补偿后的信号,fLO为变频器的本振频率;vA(t)、vB(t)表示A、B两个测站接收到的信号;
对理论时延补偿和相位补偿后的观测数据进行傅里叶变换;傅里叶变换后的信号表示为:
将傅里叶变换后的两测站的观测数据进行共轭相乘,得到互相关函数;互相关函数表示为:
相位函数表示为:
对所述互相关函数进行转化,得到关于频率的相位函数;
对观测数据进行相关后处理包括:
根据相位函数,计算干涉群时延结果以及干涉相时延结果;干涉群时延结果表示为:
干涉相时延结果表示为:
根据干涉群时延结果,解算出模糊度;模糊度表示为:
根据解算出的模糊度、所述干涉相时延结果以及所述理论时延,计算两个测站间的最终钟差,实现时间传递;
所述根据解算出的模糊度、所述干涉相时延结果以及所述理论时延,计算两个测站间的最终钟差,实现时间传递包括:
根据解算出的模糊度、所述干涉相时延结果以及所述理论时延,计算出干涉测量的时差值;干涉测量的时差值表示为:
Δt=δ+τphase;
其中,Δt指干涉测量的时差值;
根据干涉测量的时差值,引入观测方程,计算测站间的最终钟差,实现时间传递,观测方程表示为:
其中:ΔtAB指AB两站的最终钟差,是用于时间传递的修正量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210427908.5A CN114995098B (zh) | 2022-04-22 | 2022-04-22 | 基于相时延测量的皮秒级精度的卫星干涉时间传递方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210427908.5A CN114995098B (zh) | 2022-04-22 | 2022-04-22 | 基于相时延测量的皮秒级精度的卫星干涉时间传递方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114995098A CN114995098A (zh) | 2022-09-02 |
CN114995098B true CN114995098B (zh) | 2024-05-03 |
Family
ID=83025460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210427908.5A Active CN114995098B (zh) | 2022-04-22 | 2022-04-22 | 基于相时延测量的皮秒级精度的卫星干涉时间传递方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114995098B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104181550A (zh) * | 2014-08-19 | 2014-12-03 | 北京无线电计量测试研究所 | 一种基于北斗卫星导航系统的共视时间频率传递方法 |
CN107579794A (zh) * | 2017-08-21 | 2018-01-12 | 中国科学院国家授时中心 | 一种基于北斗geo导航卫星的精密共视时间频率传递方法 |
CN107766291A (zh) * | 2017-09-15 | 2018-03-06 | 中国人民解放军63920部队 | 一种获取甚长基线干涉测量中残余时延的方法 |
CN108254760A (zh) * | 2017-12-28 | 2018-07-06 | 中国科学技术大学 | 一种基于三颗量子卫星的定位与导航方法与系统 |
CN108919634A (zh) * | 2018-08-13 | 2018-11-30 | 中国科学院国家授时中心 | 一种北斗三频非差非组合观测值时间传递系统及方法 |
CN110018505A (zh) * | 2019-04-22 | 2019-07-16 | 中国科学院国家授时中心 | 一种利用通信卫星干涉测量实现舰船定位的方法 |
CN110146906A (zh) * | 2019-05-23 | 2019-08-20 | 中国科学院国家授时中心 | 基于单差载波相位观测值的远程时间传递方法 |
KR102082975B1 (ko) * | 2018-12-05 | 2020-02-28 | 한국항공우주연구원 | 위성항법장치를 이용한 신관 동작 방법 및 장치 |
CN111562594A (zh) * | 2020-05-08 | 2020-08-21 | 中国科学院国家授时中心 | 基于站间单差和模糊度固定的北斗精密时间传递方法 |
CN111983650A (zh) * | 2020-07-07 | 2020-11-24 | 北京自动化控制设备研究所 | 一种基于gnss的高精度时间传递方法 |
CN113589679A (zh) * | 2021-06-17 | 2021-11-02 | 中国科学院国家授时中心 | 一种基于干涉测量的卫星精密时间传递方法 |
CN113671543A (zh) * | 2021-05-21 | 2021-11-19 | 中国科学院国家授时中心 | 一种基于三频模式的高精度星地时间比对方法及系统 |
-
2022
- 2022-04-22 CN CN202210427908.5A patent/CN114995098B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104181550A (zh) * | 2014-08-19 | 2014-12-03 | 北京无线电计量测试研究所 | 一种基于北斗卫星导航系统的共视时间频率传递方法 |
CN107579794A (zh) * | 2017-08-21 | 2018-01-12 | 中国科学院国家授时中心 | 一种基于北斗geo导航卫星的精密共视时间频率传递方法 |
CN107766291A (zh) * | 2017-09-15 | 2018-03-06 | 中国人民解放军63920部队 | 一种获取甚长基线干涉测量中残余时延的方法 |
CN108254760A (zh) * | 2017-12-28 | 2018-07-06 | 中国科学技术大学 | 一种基于三颗量子卫星的定位与导航方法与系统 |
CN108919634A (zh) * | 2018-08-13 | 2018-11-30 | 中国科学院国家授时中心 | 一种北斗三频非差非组合观测值时间传递系统及方法 |
KR102082975B1 (ko) * | 2018-12-05 | 2020-02-28 | 한국항공우주연구원 | 위성항법장치를 이용한 신관 동작 방법 및 장치 |
CN110018505A (zh) * | 2019-04-22 | 2019-07-16 | 中国科学院国家授时中心 | 一种利用通信卫星干涉测量实现舰船定位的方法 |
CN110146906A (zh) * | 2019-05-23 | 2019-08-20 | 中国科学院国家授时中心 | 基于单差载波相位观测值的远程时间传递方法 |
CN111562594A (zh) * | 2020-05-08 | 2020-08-21 | 中国科学院国家授时中心 | 基于站间单差和模糊度固定的北斗精密时间传递方法 |
CN111983650A (zh) * | 2020-07-07 | 2020-11-24 | 北京自动化控制设备研究所 | 一种基于gnss的高精度时间传递方法 |
CN113671543A (zh) * | 2021-05-21 | 2021-11-19 | 中国科学院国家授时中心 | 一种基于三频模式的高精度星地时间比对方法及系统 |
CN113589679A (zh) * | 2021-06-17 | 2021-11-02 | 中国科学院国家授时中心 | 一种基于干涉测量的卫星精密时间传递方法 |
Non-Patent Citations (3)
Title |
---|
GNSS空间信号法时差监测方法与结果分析;章洁君;陈俊平;张益泽;;测绘通报;20170125(01);第26-29, 92页 * |
Jacqueline Walker ; Marco Genova.Simulation Study for Commercial Time Transfer Service over Geostationary Satellite.《2013 Joint European Frequency and Time Forum & International Frequency Control Symposium (EFTF/IFC)》.2013,第1-4页. * |
基于Turbo 码的空间激光脉冲时间传递方法;付靖祺等;《激光杂志》;20220221;第43卷(第3期);第20-24页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114995098A (zh) | 2022-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109358487B (zh) | 一种基于gnss精密授时的伪卫星系统及方法 | |
CN110031881B (zh) | 高精度星间激光测距辅助精密单点定位的方法 | |
RU2318222C2 (ru) | Способ и система навигации в реальном масштабе времени, использующие три несущих радиосигнала, передаваемых спутником, и ионосферные коррекции | |
CN108196272A (zh) | 一种基于实时精密单点定位的卫星导航定位装置及方法 | |
CN113671543B (zh) | 一种基于三频模式的高精度星地时间比对方法及系统 | |
US20140002302A1 (en) | Ground Location Inertial Navigation Geopositioning System (Groundlings) | |
JPH06213992A (ja) | 衛星をベースとするナビゲーションシステムにおいて位置推定の精度を改善する方法及び装置 | |
CN113671537B (zh) | 一种三频信标信号电离层信道仿真方法 | |
Liu et al. | Relative position determination of a lunar rover using high-accuracy multi-frequency same-beam VLBI | |
CN104330806A (zh) | 基于Ka测距模式的星间系统差标定方法 | |
CN104730551A (zh) | 一种星地双基地差分干涉基线坐标以及形变量测量方法 | |
CN113985455A (zh) | 一种融合多源电离层观测数据的高精度电离层建模方法 | |
Pinell | Receiver architectures for positioning with low Earth orbit satellite signals | |
Zhang et al. | Initial performance evaluation of precise point positioning with triple-frequency observations from BDS-2 and BDS-3 satellites | |
CN114935767A (zh) | 基于干涉时差测量的卫星无源定位授时方法及系统 | |
CN103543454A (zh) | 一种嵌入在移动通讯网中的卫星定轨系统 | |
CN105158781A (zh) | 一种地球gnss卫星实现外层空间用户定位的方法 | |
CN114995098B (zh) | 基于相时延测量的皮秒级精度的卫星干涉时间传递方法 | |
CN115877415B (zh) | 一种超高精度空间航天器钟差估计方法及系统 | |
Lei et al. | Geostationary orbit determination using SATRE | |
Tamazin | High resolution signal processing techniques for enhancing GPS receiver performance | |
Martínez-Belda et al. | On the potential of Galileo E5 for time transfer | |
Hwang et al. | Clock synchronization algorithm for pseudolite | |
Cheung et al. | Differencing Methods for 3D Positioning of Spacecraft | |
Tongkasem et al. | The comparison of Klobuchar model with GPS TEC model at the low geomagnetic latitude station, Thailand |
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 |