CN114839656A - 一种部分模糊度固定的时间传递快速重收敛方法 - Google Patents
一种部分模糊度固定的时间传递快速重收敛方法 Download PDFInfo
- Publication number
- CN114839656A CN114839656A CN202210304325.3A CN202210304325A CN114839656A CN 114839656 A CN114839656 A CN 114839656A CN 202210304325 A CN202210304325 A CN 202210304325A CN 114839656 A CN114839656 A CN 114839656A
- Authority
- CN
- China
- Prior art keywords
- ambiguity
- satellite
- epoch
- information
- time
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000012546 transfer Methods 0.000 title claims abstract description 42
- 238000003908 quality control method Methods 0.000 claims abstract description 6
- 239000005433 ionosphere Substances 0.000 claims description 38
- 238000003089 Pariser Parr Pople method Methods 0.000 claims description 25
- 239000005436 troposphere Substances 0.000 claims description 24
- 230000008859 change Effects 0.000 claims description 19
- 238000012360 testing method Methods 0.000 claims description 15
- 238000007667 floating Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 9
- 238000007689 inspection Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 7
- 238000011084 recovery Methods 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000007812 deficiency Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000011160 research 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/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
-
- 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/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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
本发明公开了一种部分模糊度固定的时间传递快速重收敛方法，基于历元间差分模型获取重收敛时刻模糊度信息，并对模糊度子集信息进行质量控制，形成更可靠的快速重收敛PPP时间传递方法以获取更稳定、连续、精确的GNSS(Global Navigation Satellite System)时间传递结果。本发明具有良好的时效性与灵活性，算法简单易行，在实时与事后模式中均可实现，且可拓展至多频多系统的时间比对应用。
Description
技术领域
本发明属于卫星导航技术领域，尤其涉及导航信号中断时，一种基于部分模糊度固定的快速重收敛技术实现更稳定、连续、精确的时间传递方法。
背景技术
GNSS是一种重要的定位、导航以及授时(PNT)技术，能够全天候、自动化、连续的提供高精度位置和时间服务信息。近年来，GNSS时频传递技术以其低成本、易于实现、服务范围广等优势成为近年来的研究热点。根据原始观测方程的形式以及观测值类型，GNSS大地测量时间传递技术可分为共视法(Common-View,CV)、全视法(All-In-View,AV)以及PPP(Precise Point Positioning)时间传递方法。相比于共视法，全视法以及PPP时间传递依靠整体星座参考的时间基准进行站间时间比对，进而实现全球范围内的时间传递。此外，PPP时间传递方法采用载波相位观测值，能够实现更高的时间传递精度以及短期频率传递稳定度。随着广域差分产品精度以及完备性的不断完善，PPP时间传递成为了性能最优的GNSS时间传递方法。
类似于定位应用，PPP时间传递也存在初始化时收敛时间较长的问题。一般而言，其收敛时间受卫星端群延迟、接收机端群延迟以及卫星钟差产品的初始钟偏差等因素的影响较大。当采用IGS分析中心提供的事后精密产品时，其首次收敛至1ns精度的初始化时长为20～30分钟。为了缩短收敛时长，有学者尝试引入对流层信息，形成附加约束条件的PPP时间传递模型。引入外部约束条件后可有效改善收敛期间时间传递性能，但其模糊度信息仍需要5～10分钟才能完全收敛，进而实现亚纳秒量级的时间传递。在进行长时间实时时间比对的过程中，不可避免的会出现短时间内GNSS信号中断或受到干扰等情况。当导航信号的重新捕获时，相位观测值的模糊度信息将会发生变化，从而导致PPP时间传递的重收敛现象。在重收敛过程中，接收机钟差参数长达5～30分钟的初始化时间制约了时间传递的连续性与可靠性。因此，快速重收敛成为实现长期、稳定PPP时间比对的关键问题，寻找一种准确、可靠的针对PPP时间传递的快速重收敛方法十分必要且迫切。
发明内容
为了解决上述已有技术存在的不足，本发明提出一种部分模糊度固定的时间传递快速重收敛方法，具体技术方案如下：
一种部分模糊度固定的时间传递快速重收敛方法，包括如下步骤：
S1：根据无几何距离组合观测值即GF组合观测值、无几何无电离层组合观测值即MW组合观测值，判断当前历元是否为周跳历元：
若当前历元未发生周跳，则根据GF组合观测值的变化量获取信号在电离层穿刺点天顶方向电子总含量的变化率，记录当前历元的对流层延迟及模糊度信息，更新信号在电离层穿刺点天顶方向电子总含量变化率、测站天顶方向的对流层延迟及模糊度信息；并实现非组合模式下的精密单点定位即PPP时间传递，获取时间传递结果，完成当前历元的数据处理；
若当前历元为周跳历元，则执行步骤S2-步骤S4，实现周跳修复处理，进而完成非组合模式下PPP时间传递；
S2：根据电离层穿刺点天顶方向电子总含量的变化率及测站天顶方向的对流层延迟的最新信息还原信号中断时间内的大气延迟变化量，并联合坐标信息代入历元间差分观测方程，获取重收敛时刻模糊度信息的浮点解集及相应的方差集；
S3：采用相位观测值，构造最小二乘模糊度降相关平差法即LAMBDA算法整周模糊度的搜索准则，选取方差最小的浮点解模糊度作为基准模糊度
实现搜索准则中模糊度参数与接收机钟差参数的分离；再根据ratio检验对模糊度集进行质量控制，进行部分模糊度固定检验，获取重收敛时刻模糊度集的信息；
S4：根据步骤S3中获取的重收敛时刻模糊度集的信息，结合周跳前最后一个历元确定的模糊度参数，对重收敛时刻的非差模糊度集进行恢复，然后将恢复的非差模糊度信息代入非组合观测方程，获取PPP时间传递结果，完成当前历元的数据处理。
进一步地，所述步骤S1的具体过程为：
S1-1：获取GF组合观测值及MW组合观测值；
基于原始GNSS数据的非差非组合观测方程表示为：
式中，φ及P分别为载波观测与伪距观测值，i为频点编号，j为卫星编号，
和Pi j分别为j卫星对应的i频点信号的相位和伪距观测值，λi表示i频点信号的波长，ρj为j卫星的星地距，c为真空中的光速，tr与tj分别表示接收机及j卫星的钟差参数，fi为i频点信号的频率，
与
分别为j号卫星对应的电离层与对流层投影函数，T为测站天顶方向的对流层延迟，I为卫星信号在电离层穿刺点天顶方向的电子总含量，
为j卫星对应的i频点的模糊度参数，
和
分别为载波观测值噪声与伪距观测值噪声；
根据式(1)，获取历元间差分的GF组合观测值：
其中，Δ为历元间差分算子，ΔGFj为j卫星历元间差分的GF观测值；
根据式(1)，获取卫星j的MW组合观测值
以及相应的宽巷模糊度
式中，λWL为宽巷模糊度对应的波长；
S1-2：判断步骤S1-1计算的各卫星的GF观测值以及宽巷模糊度，若不能满足式(6)，则认为发生周跳：
式中，k表示历元号，
表示j卫星在k历元的宽巷模糊度，ΔGFj(k)为j卫星在k历元的历元间差分GF观测值，
为k-1及以前历元的宽巷模糊度均值，σWL(k)为k历元宽巷模糊度的方差，μGF为给定的探测阈值；
S1-3：对于未发生周跳的历元，通过式(3)和式(7)计算出I的变化率Iv，更新信号在电离层穿刺点天顶方向的电子总含量变化率：
Iv＝ΔI/Δt (7)
其中，Δt为历元间的时间间隔；
S1-4：若未发生周跳的卫星数目不小于4颗，则根据式(1)进行非组合模式下的PPP计算：由实时差分改正数或事后精密星历提供卫星钟差参数及轨道参数，由广播星历电离层模型或外部电离层精密改正数提供电离层先验信息，由先验信息提供坐标参数，待估参数列表为：
由式(8)可知，tr为接收机钟差，表征的是接收机钟相对于卫星钟差时间基准tref的差异；I为卫星信号在电离层穿刺点天顶方向的电子总含量，T为测站天顶方向对流层延迟，
为各卫星的模糊度参数组成的模糊度向量；
根据式(8)计算的测站天顶方向对流层延迟T和模糊度向量
更新测站天顶对流层延迟和模糊度信息；
根据两站的接收机钟差，获取A、B两站的时间传递结果Tcompare(A,B)：
Tcompare(A,B)＝tA-tB (9)
其中，tA和tB分别为A、B两站的本地钟时间；
S1-5：若未发生周跳的卫星数目小于4颗，则标记为周跳历元，则继续执行步骤S2-步骤S4，实现周跳修复处理，进而完成非组合模式下PPP时间传递。
进一步地，所述步骤S2的具体过程如下：
S2-1：对于信号中断时间τ内信号在电离层穿刺点天顶方向的电子总含量变化量，根据步骤S1得到的Iv预测：
ΔIτ＝Iv·τ (10)
其中，ΔIτ为信号中断时段τ内信号在电离层穿刺点天顶方向的电子总含量的变化量；
S2-2：对流层延迟在信号中断时段τ内认为变化量为0，即
ΔT＝0 (11)
S2-3：由于时间传递为静态场景，星地距Δρ及卫星钟差变化量Δtj视为已知量，将式(10)、式(11)代入式(2)，得到历元间差分观测方程：
其中，
为载波历元间差分观测方程的验前残差，
为伪距历元间差分观测方程的验前残差，按照式(12)估计出各卫星模糊度信息的浮点解
进一步地，所述步骤S3的具体过程为：
S3-1：根据式(12)中的相位观测值确定LAMBDA算法中的搜索准则为：
式中，Ek为m阶单位矩阵，Λλ为不同频点信号波长组成的2阶对角阵，
为克罗内克算子，B为接收机钟差参数Δtr与模糊度参数
组合的系数矩阵，y为Δtr与
组成的参数矩阵，
按照协方差依次升序排列，Q为式(13)中对应的相位观测值协方差矩阵；
S3-2：由于系数阵B秩亏，需要选取协方差最小的
取整，作为模糊度搜索过程中的参考模糊度
基于式(13)、(14)及参考模糊度
采用LAMBDA算法进行整周模糊度搜索并完成ratio检验；
若ratio检验通过，则认为得到模糊度集的固定解
执行步骤S4，进而实现非差模糊度信息的恢复；
若检验未通过，则去除模糊度集中协方差最大的模糊度，并重复步骤S2进行整周模糊度集的搜索以及检验；
若模糊度子集维数小于4时仍未通过ratio检验，则模糊度集采用浮点解，执行步骤S4，进而实现非差模糊度信息的恢复。
本发明的有益效果在于：
1.本发明具有良好的可靠性，考虑了观测值质量等客观因素对模糊度固定的影响，对模糊度子集实现了质量控制，进一步提高了周跳修复后时间传递结果的可靠性。
2.本发明具有良好的时效性，涉及到的精密钟差以及轨道产品可由实时差分改正数获取，也可由事后精密星历获取；此外，算法本身简单、易行，能够在实时模式和事后模式中实现。
3.本发明具有良好的灵活性，针对不同导航系统(BDS/GPS/GALILEO)的短时间内周跳现象，本发明均可实现重收敛时刻模糊度的快速修复；且可拓展至多频率多系统时间传递下的PPP时间传递快速重收敛应用。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，通过参考附图会更加清楚的理解本发明的特征和优点，附图是示意性的而不应理解为对本发明进行任何限制，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，可以根据这些附图获得其他的附图。其中：
图1是本发明的方法流程图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本发明的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
如图1所示，本发明首先基于GF以及MW组合观测值对时频终端的双频原始观测值进行周跳探测，判断当前历元是否为周跳历元。若当前历元未发生周跳现象，则按照非差非组合的PPP算法，进行精密时频传递；同时，更新信号在电离层穿刺点天顶方向的电子总含量变化率、测站天顶方向对流层延迟以及模糊度信息。若当前历元为周跳历元，则根据历元间差分观测模型，采用最新的电离层、对流层信息并联合坐标信息，获取收敛时刻模糊度浮点解。然后根据模糊度方差选取参考模糊度，采用LAMBDA方法对整周模糊度进行搜索，获取整周模糊度集。此外，本发明结合ratio探测以及模糊度子集降维方法对搜索获得的模糊度子集进行质量控制，获取更可靠的重收敛时刻模糊度信息。最后，根据搜索的模糊度信息还原重收敛时刻的非差模糊度，并实现快速收敛后的模糊度固定时间传递。
具体地，一种部分模糊度固定的时间传递快速重收敛方法，包括如下步骤：
S1：根据无几何距离组合观测值即GF组合观测值、无几何无电离层组合观测值即MW组合观测值，判断当前历元是否为周跳历元：
若当前历元未发生周跳，则根据GF组合观测值的变化量获取信号在电离层穿刺点天顶方向电子总含量的变化率，记录当前历元的对流层延迟及模糊度信息，更新信号在电离层穿刺点天顶方向电子总含量变化率、测站天顶方向的对流层延迟及模糊度信息；并实现非组合模式下的精密单点定位即PPP时间传递，获取时间传递结果，完成当前历元的数据处理；
若当前历元为周跳历元，则执行步骤S2-步骤S4，实现周跳修复处理，进而完成非组合模式下PPP时间传递；
步骤S1的具体过程为：
S1-1：获取GF组合观测值及MW组合观测值；
基于原始GNSS数据的非差非组合观测方程表示为：
式中，φ及P分别为载波观测与伪距观测值，i为频点编号，j为卫星编号，
和Pi j分别为j卫星对应的i频点信号的相位和伪距观测值，λi表示i频点信号的波长，ρj为j卫星的星地距，c为真空中的光速，tr与tj分别表示接收机及j卫星的钟差参数，fi为i频点信号的频率，
与
分别为j号卫星对应的电离层与对流层投影函数，T为测站天顶方向的对流层延迟，I为卫星信号在电离层穿刺点天顶方向的电子总含量，
为j卫星对应的i频点的模糊度参数，
和
分别为载波观测值噪声与伪距观测值噪声；
根据式(1)，获取历元间差分的GF组合观测值：
其中，Δ为历元间差分算子，ΔGFj为j卫星历元间差分的GF观测值；
根据式(1)，获取卫星j的MW组合观测值
以及相应的宽巷模糊度
式中，λWL为宽巷模糊度对应的波长；
S1-2：判断步骤S1-1计算的各卫星的GF观测值以及宽巷模糊度，若不能满足式(6)，则认为发生周跳：
式中，k表示历元号，
表示j卫星在k历元的宽巷模糊度，ΔGFj(k)为j卫星在k历元的历元间差分GF观测值，
为k-1及以前历元的宽巷模糊度均值，σWL(k)为k历元宽巷模糊度的方差，μGF为给定的探测阈值；
S1-3：对于未发生周跳的历元，通过式(3)和式(7)计算出I的变化率Iv，更新信号在电离层穿刺点天顶方向的电子总含量变化率：
Iv＝ΔI/Δt (7)
其中，Δt为历元间的时间间隔；
S1-4：若未发生周跳的卫星数目不小于4颗，则根据式(1)进行非组合模式下的PPP计算：由实时差分改正数或事后精密星历提供卫星钟差参数及轨道参数，由广播星历电离层模型或外部电离层精密改正数提供电离层先验信息，由先验信息提供坐标参数，待估参数列表为：
由式(8)可知，tr为接收机钟差，表征的是接收机钟相对于卫星钟差时间基准tref的差异；I为卫星信号在电离层穿刺点天顶方向的电子总含量，T为测站天顶方向对流层延迟，
为各卫星的模糊度参数组成的模糊度向量；
根据式(8)计算的测站天顶方向对流层延迟T和模糊度向量
更新测站天顶对流层延迟和模糊度信息；
根据两站的接收机钟差，获取A、B两站的时间传递结果Tcompare(A,B)：
Tcompare(A,B)＝tA-tB (9)
其中，tA和tB分别为A、B两站的本地钟时间；
S1-5：若未发生周跳的卫星数目小于4颗，则标记为周跳历元，则继续执行步骤S2-步骤S4，实现周跳修复处理，进而完成非组合模式下PPP时间传递。
S2：根据电离层穿刺点天顶方向电子总含量的变化率及测站天顶方向的对流层延迟的最新信息还原信号中断时间内的大气延迟变化量，并联合坐标信息代入历元间差分观测方程，获取重收敛时刻模糊度信息的浮点解集及相应的方差集；
步骤S2的具体过程如下：
S2-1：对于信号中断时间τ内信号在电离层穿刺点天顶方向的电子总含量变化量，根据步骤S1得到的Iv预测：
ΔIτ＝Iv·τ (10)
其中，ΔIτ为信号中断时段τ内信号在电离层穿刺点天顶方向的电子总含量的变化量；
S2-2：对流层延迟在信号中断时段τ内认为变化量为0，即
ΔT＝0 (11)
S2-3：由于时间传递为静态场景，星地距Δρ及卫星钟差变化量Δtj视为已知量，将式(10)、式(11)代入式(2)，得到历元间差分观测方程：
其中，
为载波历元间差分观测方程的验前残差，
为伪距历元间差分观测方程的验前残差，按照式(12)估计出各卫星模糊度信息的浮点解
S3：采用相位观测值，构造最小二乘模糊度降相关平差法即LAMBDA算法整周模糊度的搜索准则，选取方差最小的浮点解模糊度作为基准模糊度
实现搜索准则中模糊度参数与接收机钟差参数的分离；再根据ratio检验对模糊度集进行质量控制，进行部分模糊度固定检验，获取重收敛时刻模糊度集的信息；
步骤S3的具体过程为：
S3-1：根据式(12)中的相位观测值确定LAMBDA算法中的搜索准则为：
式中，Ek为m阶单位矩阵，Λλ为不同频点信号波长组成的2阶对角阵，
为克罗内克算子，B为接收机钟差参数Δtr与模糊度参数
组合的系数矩阵，y为Δtr与
组成的参数矩阵，
按照协方差依次升序排列，Q为式(13)中对应的相位观测值协方差矩阵；
S3-2：由于系数阵B秩亏，需要选取协方差最小的
取整，作为模糊度搜索过程中的参考模糊度
基于式(13)、(14)及参考模糊度
采用LAMBDA算法进行整周模糊度搜索并完成ratio检验；
若ratio检验通过，则认为得到模糊度集的固定解
执行步骤S4，进而实现非差模糊度信息的恢复；
若检验未通过，则去除模糊度集中协方差最大的模糊度，并重复步骤S2进行整周模糊度集的搜索以及检验；
若模糊度子集维数小于4时仍未通过ratio检验，则模糊度集采用浮点解，执行步骤S4，进而实现非差模糊度信息的恢复。
S4：根据步骤S3中获取的重收敛时刻模糊度集的信息，结合周跳前最后一个历元确定的模糊度参数，对重收敛时刻的非差模糊度集进行恢复，然后将恢复的非差模糊度信息代入非组合观测方程，获取PPP时间传递结果，完成当前历元的数据处理。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种部分模糊度固定的时间传递快速重收敛方法，其特征在于，包括如下步骤：
S1：根据无几何距离组合观测值即GF组合观测值、无几何无电离层组合观测值即MW组合观测值，判断当前历元是否为周跳历元：
若当前历元未发生周跳，则根据GF组合观测值的变化量获取信号在电离层穿刺点天顶方向电子总含量的变化率，记录当前历元的对流层延迟及模糊度信息，更新信号在电离层穿刺点天顶方向电子总含量变化率、测站天顶方向的对流层延迟及模糊度信息；并实现非组合模式下的精密单点定位即PPP时间传递，获取时间传递结果，完成当前历元的数据处理；
若当前历元为周跳历元，则执行步骤S2-步骤S4，实现周跳修复处理，进而完成非组合模式下PPP时间传递；
S2：根据电离层穿刺点天顶方向电子总含量的变化率及测站天顶方向的对流层延迟的最新信息还原信号中断时间内的大气延迟变化量，并联合坐标信息代入历元间差分观测方程，获取重收敛时刻模糊度信息的浮点解集及相应的方差集；
S3：采用相位观测值，构造最小二乘模糊度降相关平差法即LAMBDA算法整周模糊度的搜索准则，选取方差最小的浮点解模糊度作为基准模糊度
实现搜索准则中模糊度参数与接收机钟差参数的分离；再根据ratio检验对模糊度集进行质量控制，进行部分模糊度固定检验，获取重收敛时刻模糊度集的信息；
S4：根据步骤S3中获取的重收敛时刻模糊度集的信息，结合周跳前最后一个历元确定的模糊度参数，对重收敛时刻的非差模糊度集进行恢复，然后将恢复的非差模糊度信息代入非组合观测方程，获取PPP时间传递结果，完成当前历元的数据处理。
2.根据权利要求1所述的方法，其特征在于，所述步骤S1的具体过程为：
S1-1：获取GF组合观测值及MW组合观测值；
基于原始GNSS数据的非差非组合观测方程表示为：
式中，φ及P分别为载波观测与伪距观测值，i为频点编号，j为卫星编号，
和Pi j分别为j卫星对应的i频点信号的相位和伪距观测值，λi表示i频点信号的波长，ρj为j卫星的星地距，c为真空中的光速，tr与tj分别表示接收机及j卫星的钟差参数，fi为i频点信号的频率，VI j与
分别为j号卫星对应的电离层与对流层投影函数，T为测站天顶方向的对流层延迟，I为卫星信号在电离层穿刺点天顶方向的电子总含量，
为j卫星对应的i频点的模糊度参数，
和
分别为载波观测值噪声与伪距观测值噪声；
根据式(1)，获取历元间差分的GF组合观测值：
其中，Δ为历元间差分算子，ΔGFj为j卫星历元间差分的GF观测值；
根据式(1)，获取卫星j的MW组合观测值
以及相应的宽巷模糊度
式中，λWL为宽巷模糊度对应的波长；
S1-2：判断步骤S1-1计算的各卫星的GF观测值以及宽巷模糊度，若不能满足式(6)，则认为发生周跳：
式中，k表示历元号，
表示j卫星在k历元的宽巷模糊度，ΔGFj(k)为j卫星在k历元的历元间差分GF观测值，
为k-1及以前历元的宽巷模糊度均值，σWL(k)为k历元宽巷模糊度的方差，μGF为给定的探测阈值；
S1-3：对于未发生周跳的历元，通过式(3)和式(7)计算出I的变化率Iv，更新信号在电离层穿刺点天顶方向的电子总含量变化率：
Iv＝ΔI/Δt (7)
其中，Δt为历元间的时间间隔；
S1-4：若未发生周跳的卫星数目不小于4颗，则根据式(1)进行非组合模式下的PPP计算：由实时差分改正数或事后精密星历提供卫星钟差参数及轨道参数，由广播星历电离层模型或外部电离层精密改正数提供电离层先验信息，由先验信息提供坐标参数，待估参数列表为：
由式(8)可知，tr为接收机钟差，表征的是接收机钟相对于卫星钟差时间基准tref的差异；I为卫星信号在电离层穿刺点天顶方向的电子总含量，T为测站天顶方向对流层延迟，
为各卫星的模糊度参数组成的模糊度向量；
根据式(8)计算的测站天顶方向对流层延迟T和模糊度向量
更新测站天顶对流层延迟和模糊度信息；
根据两站的接收机钟差，获取A、B两站的时间传递结果Tcompare(A,B)：
Tcompare(A,B)＝tA-tB (9)
其中，tA和tB分别为A、B两站的本地钟时间；
S1-5：若未发生周跳的卫星数目小于4颗，则标记为周跳历元，则继续执行步骤S2-步骤S4，实现周跳修复处理，进而完成非组合模式下PPP时间传递。
3.根据权利要求1或2所述的方法，其特征在于，所述步骤S2的具体过程如下：
S2-1：对于信号中断时间τ内信号在电离层穿刺点天顶方向的电子总含量变化量，根据步骤S1得到的Iv预测：
ΔIτ＝Iv·τ (10)
其中，ΔIτ为信号中断时段τ内信号在电离层穿刺点天顶方向的电子总含量的变化量；
S2-2：对流层延迟在信号中断时段τ内认为变化量为0，即
ΔT＝0 (11)
S2-3：由于时间传递为静态场景，星地距Δρ及卫星钟差变化量Δtj视为已知量，将式(10)、式(11)代入式(2)，得到历元间差分观测方程：
其中，
为载波历元间差分观测方程的验前残差，
为伪距历元间差分观测方程的验前残差，按照式(12)估计出各卫星模糊度信息的浮点解
4.根据权利要求3所述的方法，其特征在于，所述步骤S3的具体过程为：
S3-1：根据式(12)中的相位观测值确定LAMBDA算法中的搜索准则为：
式中，Ek为m阶单位矩阵，Λλ为不同频点信号波长组成的2阶对角阵，
为克罗内克算子，B为接收机钟差参数Δtr与模糊度参数
组合的系数矩阵，y为Δtr与
组成的参数矩阵，
按照协方差依次升序排列，Q为式(13)中对应的相位观测值协方差矩阵；
S3-2：由于系数阵B秩亏，需要选取协方差最小的
取整，作为模糊度搜索过程中的参考模糊度
基于式(13)、(14)及参考模糊度
采用LAMBDA算法进行整周模糊度搜索并完成ratio检验；
若ratio检验通过，则认为得到模糊度集的固定解
执行步骤S4，进而实现非差模糊度信息的恢复；
若检验未通过，则去除模糊度集中协方差最大的模糊度，并重复步骤S2进行整周模糊度集的搜索以及检验；
若模糊度子集维数小于4时仍未通过ratio检验，则模糊度集采用浮点解，执行步骤S4，进而实现非差模糊度信息的恢复。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210304325.3A CN114839656B (zh) | 2022-03-25 | 2022-03-25 | 一种部分模糊度固定的时间传递快速重收敛方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210304325.3A CN114839656B (zh) | 2022-03-25 | 2022-03-25 | 一种部分模糊度固定的时间传递快速重收敛方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114839656A true CN114839656A (zh) | 2022-08-02 |
| CN114839656B CN114839656B (zh) | 2024-06-04 |
Family
ID=82562824
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210304325.3A Active CN114839656B (zh) | 2022-03-25 | 2022-03-25 | 一种部分模糊度固定的时间传递快速重收敛方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114839656B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115933356A (zh) * | 2023-01-09 | 2023-04-07 | 北京航空航天大学 | 一种虚拟原子钟的高精度时间同步系统和方法 |
| CN116338752A (zh) * | 2023-01-29 | 2023-06-27 | 北京航空航天大学 | 一种用于载波相位单差时间传递的精确模糊度固定方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160077214A1 (en) * | 2014-09-15 | 2016-03-17 | Fugro N.V. | Precise gnss positioning system with improved ambiguity estimation |
| US20160077213A1 (en) * | 2014-09-15 | 2016-03-17 | Fugro N.V. | Integer ambiguity-fixed precise point positioning method and system |
| CN107132558A (zh) * | 2017-06-13 | 2017-09-05 | 武汉大学 | 惯性辅助的多频多模gnss周跳修复方法及系统 |
| CN108445518A (zh) * | 2018-03-16 | 2018-08-24 | 中国科学院数学与系统科学研究院 | 一种基于双差模糊度固定解约束的gnss精密时间传递方法 |
| WO2019144528A1 (zh) * | 2018-01-29 | 2019-08-01 | 东南大学 | 基于模糊度紧约束的多星座基准站间模糊度快速解算方法及其应用 |
| CN112014860A (zh) * | 2020-07-20 | 2020-12-01 | 中国科学院空天信息创新研究院 | 一种基于北斗ppp-rtk的低轨卫星时空基准建立方法 |
| US20210072406A1 (en) * | 2018-06-01 | 2021-03-11 | Zhejiang Yat Electrical Appliance Co., Ltd | Gnss-rtk-based positioning method |
-
2022
- 2022-03-25 CN CN202210304325.3A patent/CN114839656B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160077214A1 (en) * | 2014-09-15 | 2016-03-17 | Fugro N.V. | Precise gnss positioning system with improved ambiguity estimation |
| US20160077213A1 (en) * | 2014-09-15 | 2016-03-17 | Fugro N.V. | Integer ambiguity-fixed precise point positioning method and system |
| CN107132558A (zh) * | 2017-06-13 | 2017-09-05 | 武汉大学 | 惯性辅助的多频多模gnss周跳修复方法及系统 |
| WO2019144528A1 (zh) * | 2018-01-29 | 2019-08-01 | 东南大学 | 基于模糊度紧约束的多星座基准站间模糊度快速解算方法及其应用 |
| CN108445518A (zh) * | 2018-03-16 | 2018-08-24 | 中国科学院数学与系统科学研究院 | 一种基于双差模糊度固定解约束的gnss精密时间传递方法 |
| US20210072406A1 (en) * | 2018-06-01 | 2021-03-11 | Zhejiang Yat Electrical Appliance Co., Ltd | Gnss-rtk-based positioning method |
| CN112014860A (zh) * | 2020-07-20 | 2020-12-01 | 中国科学院空天信息创新研究院 | 一种基于北斗ppp-rtk的低轨卫星时空基准建立方法 |
Non-Patent Citations (5)
| Title |
|---|
| 张小红;胡家欢;任晓东;: "PPP/PPP-RTK新进展与北斗/GNSS PPP定位性能比较", 测绘学报, no. 09, 15 September 2020 (2020-09-15) * |
| 苏珂;金双根;: "BDS/Galileo四频精密单点定位模型性能分析与比较", 测绘学报, no. 09, 15 September 2020 (2020-09-15) * |
| 赵兴旺;张翠英;: "精密单点定位部分模糊度固定方法", 中国空间科学技术, no. 01, 25 February 2013 (2013-02-25) * |
| 邹璇;陈亮;李宗楠;唐卫明;施闯;: "基于PPP-AR模糊度快速重收敛技术提供可靠的URTK服务", 武汉大学学报(信息科学版), no. 10, 5 October 2016 (2016-10-05) * |
| 郑艳丽;刘经南;宋伟伟;孙华;: "基于区域增强信息的PPP的快速收敛算法", 大地测量与地球动力学, no. 04, 15 August 2012 (2012-08-15) * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115933356A (zh) * | 2023-01-09 | 2023-04-07 | 北京航空航天大学 | 一种虚拟原子钟的高精度时间同步系统和方法 |
| CN115933356B (zh) * | 2023-01-09 | 2023-08-22 | 北京航空航天大学 | 一种虚拟原子钟的高精度时间同步系统和方法 |
| CN116338752A (zh) * | 2023-01-29 | 2023-06-27 | 北京航空航天大学 | 一种用于载波相位单差时间传递的精确模糊度固定方法 |
| CN116338752B (zh) * | 2023-01-29 | 2023-11-07 | 北京航空航天大学 | 一种用于载波相位单差时间传递的精确模糊度固定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114839656B (zh) | 2024-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109001786B (zh) | 一种基于导航卫星和低轨增强卫星的定位方法和系统 | |
| CN110231037B (zh) | 一种附加钟差模型约束的gnss机动卫星定轨方法 | |
| CN114839656B (zh) | 一种部分模糊度固定的时间传递快速重收敛方法 | |
| US8046169B2 (en) | System and method for determining the geographic location of a device | |
| EP3430436A1 (en) | Rapid recovery of precise position after temporary signal loss | |
| CN110007320B (zh) | 网络rtk解算方法 | |
| CN109709579B (zh) | 一种基于用户测距误差实时估计的gnss卫星星历故障检测方法 | |
| EP3430429B1 (en) | Satellite navigation receiver with improved ambiguity resolution | |
| US20090177382A1 (en) | Calibration of a Navigation System | |
| CN107942346B (zh) | 一种高精度gnss电离层tec观测值提取方法 | |
| CA2780675A1 (en) | Detection and correction of anomalous measurements and ambiguity resolution in a global navigation satellite system receiver | |
| CN102998681A (zh) | 一种卫星导航系统的高频钟差估计方法 | |
| CN104865587A (zh) | 基于geo星座粗定位的北斗接收机快速定位方法 | |
| CN109752739B (zh) | 一种观测数据处理方法、装置、终端设备及存储介质 | |
| CN111913201A (zh) | 一种gnss差分定位方法、装置及计算机可读存储介质 | |
| Dai et al. | Innovative algorithms to improve long range RTK reliability and availability | |
| CN114562992A (zh) | 一种基于因子图及场景约束的多径环境组合导航方法 | |
| CN109212563A (zh) | 北斗/gps三频周跳探测与修复方法 | |
| CN109946722B (zh) | 一种多系统多频段定位方法及系统 | |
| CN104898144A (zh) | 一种北斗系统参考站整周模糊度单历元确定方法 | |
| CN109143286B (zh) | 一种顾及非模型化误差的卫星导航定位方法 | |
| Li et al. | Robust Kalman filter for RTK positioning under signal-degraded scenarios | |
| Yin et al. | A novel cycle slips detection model for the high precision positioning | |
| CN115308781B (zh) | 基于bdgim辅助的相位平滑伪距高精度时间传递方法 | |
| CN115201880A (zh) | 一种精密单点定位方法、装置和计算机可读存储介质 |
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 |