CN111562594A - 基于站间单差和模糊度固定的北斗精密时间传递方法 - Google Patents
基于站间单差和模糊度固定的北斗精密时间传递方法 Download PDFInfo
- Publication number
- CN111562594A CN111562594A CN202010383638.3A CN202010383638A CN111562594A CN 111562594 A CN111562594 A CN 111562594A CN 202010383638 A CN202010383638 A CN 202010383638A CN 111562594 A CN111562594 A CN 111562594A
- Authority
- CN
- China
- Prior art keywords
- difference
- ambiguity
- data
- satellite
- model
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000011159 matrix material Substances 0.000 claims abstract description 8
- 239000005436 troposphere Substances 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 5
- 101100522111 Oryza sativa subsp. japonica PHT1-11 gene Proteins 0.000 description 2
- 238000005295 random walk Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/05—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing aiding data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/256—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to timing, e.g. time of week, code phase, timing offset
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了基于站间单差和模糊度固定的北斗精密时间传递方法,属于卫星时间传递技术领域。首先,基于时间主站和用户站的北斗观测数据,组成站间单差观测方程并约束测站坐标进行解算,获得初始的单差模糊度浮点解和对应的方差阵;其次,选定参考卫星,对求解的单差模糊度进行星间差分,进行双差模糊度的归整和固定,并将该参考卫星的单差浮点解作为基准,恢复模糊度固定后其他的单差模糊度的浮点解;再后,将固定后得到的全部单差模糊度浮点解代入原始观测方程,求解两站之间的时差参数。本发明提出的技术方法不需要外部辅助条件即可以简单实施,同时也适用于广播星历,方便实时应用。
Description
技术领域
本发明属于卫星时间传递技术领域,具体涉及基于站间单差和模糊度固定的北斗精密时间传递方法。
背景技术
基于卫星导航系统的时间传递方法因成本低、连续性好、全天候等特点成为精密时间传递的重要手段之一。卫星时间传递中常规方法有共视法(CV),全视法(AV)和载波相位(CP)。CP方法因为使用了高精度载波相位观测值,时间传递精度高成为当前的主流技术。但是,当前的CP时间传递方法存在两个局限,一是采用非差模式,观测值中很多误差难以高效处理,影响时间传递的精度;二是常规的CP时间传递解算中载波相位模糊度都是浮点解,难以高效准确固定,对收敛速度和精度也有影响。
如何高效修正观测误差并实现CP中相位模糊度的固定,进一步提高CP时间传递的精度和收敛速度,对时间服务具有重要应用价值。
发明内容
针对现有的技术问题,本发明提供基于站间单差和模糊度固定的北斗精密时间传递方法,首先,基于时间主站和用户站的北斗观测数据,组成站间单差观测方程并约束测站坐标进行解算,获得初始的单差模糊度浮点解和对应的方差阵;其次,选定参考卫星,对求解的单差模糊度进行星间差分,进行双差模糊度的归整和固定,并将该参考卫星的单差浮点解作为基准,恢复模糊度固定后其他的单差模糊度的浮点解;再后,将固定后得到的全部单差模糊度浮点解代入原始观测方程,求解两站之间的时差参数。
本发明的技术方案是:基于站间单差和模糊度固定的北斗精密时间传递方法,主要包括以下步骤:
第一步:数据获取
获得主站和用户站上北斗卫星的双频伪距、相位观测数据以及数据处理需要的辅助数据;
第二步:数据预处理
首先,对全部数据进行质量检查、粗差剔除,删除无卫星星历、观测值不完整的数据以及删除不共同观测的卫星数据,得到干净可用的数据,并进行周跳探测,给出周跳探测结果;
第三步:误差模型改正
对预处理后干净的数据进行相对论、潮汐、天线相位中心、对流层和地球自转误差的修正;
第四步:建立单差观测模型
对于全部公共卫星,进行站间差分,形成单差观测方程,并确定随机模型,通过组建站间单差观测方程,可以很好的消除或削弱非模型化误差,提高观测值误差的修正水平,进而提高时间传递精度;
第五步:单差数据解算
基于第四步中得到的单差观测方程和随机模型,进行最小二乘求解,估计得到参数的初值以及对应的方差-协方差阵;
第六步:模糊度固定
基于第五步中的单差模糊度初值和对应的方差-协方差阵,选定参考卫星,转化为具有整周特性的双差模糊度,采用Lambda算法进行固定,并将参考卫星的站间单差模糊度浮点解作为基准,恢复模糊度固定后全部单差模糊度的浮点解,通过采用模糊度固定技术,实现求解参数的快速收敛和精准估计,提高时间传递性能;
第七步:时差参数确定
将第六步中求解的模糊度固定后全部单差模糊度浮点解代入到第四步的单差观测模型中,再次求解得到时差参数,即完成时间传递工作。
进一步地,第一步中数据处理需要的辅助数据,包括精密星历、天线相位中心参数、地球自转参数、测站准确坐标、潮汐参数,用来进行误差修正。
进一步地,第三步中相对论和潮汐误差使用IERS Conventions 2010中指定的模型改正,天线相位中心误差采用igs14.atx模型改正,对流层误差采用Saastamoinen模型改正,地球自转误差使用IERS EOP C04模型改正。
进一步地,第四步中形成的单差观测方程为:
LP-U-LP-R=ΔρP+A·Δx+F·Δtrop+C·Δdt+ΔmodelP+ΔεP
LΦ-U-LΦ-R=ΔρΦ+A·Δx+F·Δtrop+C·Δdt+λ·ΔN+ΔmodelΦ+ΔεΦ其中,P和Φ代表伪距和相位,L为无电离层组合观测值,A为站星间单位方向矢量,x为三维坐标改正量,Δ表示站间差分,R和U代表主站和用户站,trop代表对流层残差,F为投影系数,ρ为站星几何距离,C为光速,dt为时差,model为各种模型化误差总和,ε代表观测噪声。
进一步地,第四步中形成的随机模型如下:
更进一步地,第五步中测站坐标进行强约束估计,对流层残差进行随机游走估计,单差模糊度连续无周跳情况下当作常数估计,在有周跳时重新初始化,时差参数进行白噪声估计,从而得到各参数的估计值。
本发明的有益效果是:本发明提供一种基于站间单差和模糊度固定的北斗精密时间传递方法,具备以下的优点:
(1)本发明将主站和用户站的观测组成站间单差观测方程,进行参数估计获取单差模糊度初值和对应的方差-协方差阵,并选定参考卫星,实现双差模糊度归整和固定,同时确定参考星单差模糊度浮点解为基准,恢复模糊度固定后的全部单差模糊度浮点解,进而代入单差模型求解得到时差参数,实现精密时间传递。
(2)站间单差模式,可以消除或削弱轨道、钟差等共性误差,提高观测值误差修正水平,提高时间传递精度,相对于常规CP方法中非差观测值误差难以高效修正,本发明通过组建站间单差观测方程,可以很好的消除或削弱非模型化误差,提高观测值误差的修正水平,进而提高时间传递精度。
(3)模糊度采用固定技术,实现求解参数的快速收敛和精准估计,提高时间传递性能,相对于常规CP时间传递方法采用模糊度浮点解,其参数收敛过程很慢,本方法采用模糊度固定模式,不仅加快了收敛速度,同时也提高了时间传递精度。
(4)方法实施简单,同时也适用于广播星历,方便实时应用,不需要外部辅助条件即可以简单实施,方便操作,同时,因为采用站间单差模式,其卫星轨道、钟差、大气等误差相关性很高,可以有效消除,因此本方法也实用于广播星历,方便实时应用。
附图说明
图1为本发明的基于站间单差和模糊度固定的北斗精密时间传递技术流程图;
图2为本发明的PT11测站获取的时差序列。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
实施例:如图1、2所示的基于站间单差和模糊度固定的北斗精密时间传递方法,主要包括以下步骤:
第一步:数据获取
获得主站PTBB和用户站PT11上北斗卫星的双频伪距、相位观测数据以及数据处理需要的辅助数据,其中,数据处理需要的辅助数据包括精密星历,天线相位中心参数、地球自转参数、测站准确坐标、潮汐参数,用来进行误差修正,2020年001天单天数据为例,数据采样间隔为30s;
第二步:数据预处理
首先,对全部数据进行质量检查、粗差剔除,删除无卫星星历、观测值不完整的数据以及删除不共同观测的卫星数据,得到干净可用的数据,并进行周跳探测,给出周跳探测结果;
第三步:误差模型改正
对预处理后干净的数据进行相对论、潮汐、天线相位中心、对流层和地球自转误差的修正,其中相对论和潮汐误差使用IERS Conventions 2010中指定的模型改正,天线相位中心误差采用igs14.atx模型改正,对流层误差采用Saastamoinen模型改正,地球自转误差使用IERS EOP C04模型改正;
第四步:建立单差观测模型
对于全部公共卫星,进行站间差分,形成单差观测方程,单差观测方程具体为:
LP-U-LP-R=ΔρP+A·Δx+F·Δtrop+C·Δdt+ΔmodelP+ΔεP
LΦ-U-LΦ-R=ΔρΦ+A·Δx+F·Δtrop+C·Δdt+λ·ΔN+ΔmodelΦ+ΔεΦ
(1)
其中,P和Φ代表伪距和相位,L为无电离层组合观测值,A为站星间单位方向矢量,x为三维坐标改正量,Δ表示站间差分,R和U代表主站和用户站,trop代表对流层残差,F为投影系数,ρ为站星几何距离,C为光速,dt为时差,model为各种模型化误差总和,ε代表观测噪声;
并确定随机模型,随机模型如下:
其中,a为观测值精度,伪距观测值一般设为0.2-0.3m,相位观测值一般设为0.002-0.003m,为两测站卫星高度角平均值,单位为弧度,通过组建站间单差观测方程,可以很好的消除或削弱非模型化误差,提高观测值误差的修正水平,进而提高时间传递精度;
第五步:单差数据解算
基于第四步中得到的单差观测方程和随机模型,进行最小二乘求解,估计得到参数的初值以及对应的方差-协方差阵,测站坐标进行强约束估计,对流层残差进行随机游走估计,单差模糊度连续无周跳情况下当作常数估计,在有周跳时重新初始化,时差参数进行白噪声估计;
第六步:模糊度固定
基于第五步中的单差模糊度初值和对应的方差-协方差阵,选定参考卫星,转化为具有整周特性的双差模糊度,采用Lambda算法进行固定,并将参考卫星的站间单差模糊度浮点解作为基准,恢复模糊度固定后全部单差模糊度的浮点解,通过采用模糊度固定技术,实现求解参数的快速收敛和精准估计,提高时间传递性能;
第七步:时差参数确定
将第六步中求解的模糊度固定后全部单差模糊度浮点解代入到第四步的单差观测模型中,再次求解得到时差参数,即完成时间传递工作。
Claims (5)
1.基于站间单差和模糊度固定的北斗精密时间传递方法,其特征在于,主要包括以下步骤:
第一步:数据获取
获得主站和用户站上北斗卫星的双频伪距、相位观测数据以及数据处理需要的辅助数据;
第二步:数据预处理
首先,对全部数据进行质量检查、粗差剔除,删除无卫星星历、观测值不完整的数据以及删除不共同观测的卫星数据,得到干净可用的数据,并进行周跳探测,给出周跳探测结果;
第三步:误差模型改正
对预处理后干净的数据进行相对论、潮汐、天线相位中心、对流层和地球自转误差的修正;
第四步:建立单差观测模型
对于全部公共卫星,进行站间差分,形成单差观测方程,并确定随机模型;
第五步:单差数据解算
基于第四步中得到的单差观测方程和随机模型,进行最小二乘求解,估计得到参数的初值以及对应的方差-协方差阵;
第六步:模糊度固定
基于第五步中的单差模糊度初值和对应的方差-协方差阵,选定参考卫星,转化为具有整周特性的双差模糊度,采用Lambda算法进行固定,并将参考卫星的站间单差模糊度浮点解作为基准,恢复模糊度固定后全部单差模糊度的浮点解;
第七步:时差参数确定
将第六步中求解的模糊度固定后全部单差模糊度浮点解代入到第四步的单差观测模型中,再次求解得到时差参数,即完成时间传递工作。
2.根据权利要求1所述的基于站间单差和模糊度固定的北斗精密时间传递方法,其特征在于,第一步中数据处理需要的辅助数据包括精密星历、天线相位中心参数、地球自转参数、测站准确坐标、潮汐参数。
3.根据权利要求1所述的基于站间单差和模糊度固定的北斗精密时间传递方法,其特征在于,第三步中相对论和潮汐误差使用IERS Conventions 2010中指定的模型改正,天线相位中心误差采用igs14.atx模型改正,对流层误差采用Saastamoinen模型改正,地球自转误差使用IERS EOP C04模型改正。
4.根据权利要求1所述的基于站间单差和模糊度固定的北斗精密时间传递方法,其特征在于,第四步中形成的单差观测方程为:
LP-U-LP-R=ΔρP+A·Δx+F·Δtrop+C·Δdt+ΔmodelP+ΔεP
LΦ-U-LΦ-R=ΔρΦ+A·Δx+F·Δtrop+C·Δdt+λ·ΔN+ΔmodelΦ+ΔεΦ
其中,P和Φ代表伪距和相位,L为无电离层组合观测值,A为站星间单位方向矢量,x为三维坐标改正量,Δ表示站间差分,R和U代表主站和用户站,trop代表对流层残差,F为投影系数,ρ为站星几何距离,C为光速,dt为时差,model为各种模型化误差总和,ε代表观测噪声。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010383638.3A CN111562594A (zh) | 2020-05-08 | 2020-05-08 | 基于站间单差和模糊度固定的北斗精密时间传递方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010383638.3A CN111562594A (zh) | 2020-05-08 | 2020-05-08 | 基于站间单差和模糊度固定的北斗精密时间传递方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111562594A true CN111562594A (zh) | 2020-08-21 |
Family
ID=72068034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010383638.3A Pending CN111562594A (zh) | 2020-05-08 | 2020-05-08 | 基于站间单差和模糊度固定的北斗精密时间传递方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111562594A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114995098A (zh) * | 2022-04-22 | 2022-09-02 | 中国科学院国家授时中心 | 基于相时延测量的皮秒级精度的卫星干涉时间传递方法 |
CN115225245A (zh) * | 2022-09-20 | 2022-10-21 | 中国科学院国家授时中心 | 一种非差非组合ppp-rtk长基线实时时间传递方法 |
CN116660958A (zh) * | 2023-07-31 | 2023-08-29 | 武汉大学 | 整周模糊度并行滤波解算方法、装置、设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106569242A (zh) * | 2016-11-04 | 2017-04-19 | 中国科学院测量与地球物理研究所 | 固定参考星的gnss单差处理方法 |
CN108445518A (zh) * | 2018-03-16 | 2018-08-24 | 中国科学院数学与系统科学研究院 | 一种基于双差模糊度固定解约束的gnss精密时间传递方法 |
CN110146906A (zh) * | 2019-05-23 | 2019-08-20 | 中国科学院国家授时中心 | 基于单差载波相位观测值的远程时间传递方法 |
CN111045034A (zh) * | 2019-12-13 | 2020-04-21 | 北京航空航天大学 | 基于广播星历的gnss多系统实时精密时间传递方法及系统 |
-
2020
- 2020-05-08 CN CN202010383638.3A patent/CN111562594A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106569242A (zh) * | 2016-11-04 | 2017-04-19 | 中国科学院测量与地球物理研究所 | 固定参考星的gnss单差处理方法 |
CN108445518A (zh) * | 2018-03-16 | 2018-08-24 | 中国科学院数学与系统科学研究院 | 一种基于双差模糊度固定解约束的gnss精密时间传递方法 |
CN110146906A (zh) * | 2019-05-23 | 2019-08-20 | 中国科学院国家授时中心 | 基于单差载波相位观测值的远程时间传递方法 |
CN111045034A (zh) * | 2019-12-13 | 2020-04-21 | 北京航空航天大学 | 基于广播星历的gnss多系统实时精密时间传递方法及系统 |
Non-Patent Citations (1)
Title |
---|
黄令勇: "三频GNSS精密定位理论与方法研究", 《中国优秀博硕士学位论文全文数据库(博士)基础科学辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114995098A (zh) * | 2022-04-22 | 2022-09-02 | 中国科学院国家授时中心 | 基于相时延测量的皮秒级精度的卫星干涉时间传递方法 |
CN114995098B (zh) * | 2022-04-22 | 2024-05-03 | 中国科学院国家授时中心 | 基于相时延测量的皮秒级精度的卫星干涉时间传递方法 |
CN115225245A (zh) * | 2022-09-20 | 2022-10-21 | 中国科学院国家授时中心 | 一种非差非组合ppp-rtk长基线实时时间传递方法 |
CN116660958A (zh) * | 2023-07-31 | 2023-08-29 | 武汉大学 | 整周模糊度并行滤波解算方法、装置、设备及存储介质 |
CN116660958B (zh) * | 2023-07-31 | 2023-10-20 | 武汉大学 | 整周模糊度并行滤波解算方法、装置、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
CN111562594A (zh) | 基于站间单差和模糊度固定的北斗精密时间传递方法 | |
CN109085628B (zh) | 一种整周模糊度的固定方法及系统 | |
CN105301617B (zh) | 一种卫星导航系统中的整周模糊度有效性检验方法 | |
CN111580132B (zh) | 基于时间实验室增强信息的北斗局域精密时间传递方法 | |
CN112835082A (zh) | 一种gnss区域增强电离层和对流层大气产品质量指标计算方法、电子设备、存储介质 | |
CN111308504A (zh) | 一种基于相位观测值实时探测北斗卫星轨道机动的方法 | |
CN111694030A (zh) | 一种基于格网虚拟观测值的bds局域差分方法及系统 | |
CN113253314B (zh) | 一种低轨卫星间时间同步方法及系统 | |
CN109884679B (zh) | 一种单模gnss系统跨频点混合双差rtk解算方法 | |
CN113848572A (zh) | 一种基于大气误差增强的多频ppp序贯单历元定位方法 | |
CN112146557A (zh) | 一种基于gnss的实时桥梁变形监测系统及方法 | |
CN113325446B (zh) | 一种多模共频gnss载波相位时间传递方法及系统 | |
CN110208835A (zh) | 一种基于消电离层组合的跨系统紧组合差分定位方法 | |
CN111290004A (zh) | 伪距差分定位方法、装置、电子设备及存储介质 | |
CN116540303A (zh) | 一种历元间差分观测方程建立方法和地震位移计算方法 | |
CN114928428A (zh) | 一种基于增强手段的精密时间传递方法和系统 | |
CN116026226A (zh) | 一种半遮挡环境下的水闸变形监测方法及系统 | |
CN109613582B (zh) | 一种车载实时单频米级伪距定位方法 | |
CN112198540B (zh) | 一种基于动态网络基站的多模多频载波相位定位方法 | |
CN111650612A (zh) | 一种基于双差模式的实时精密时间传递方法 | |
CN115933356B (zh) | 一种虚拟原子钟的高精度时间同步系统和方法 | |
CN115144878B (zh) | 基于ppp的短距离大高差nrtk对流层延迟改正方法 | |
CN114779301B (zh) | 一种基于广播星历的卫星导航实时精密单点定位方法 | |
CN114355420B (zh) | 一种分布式北斗位置服务中心ppp产品定位方法及装置 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200821 |