CN111965673B - 基于多gnss的单频精密单点定位算法的时间频率传递方法 - Google Patents

基于多gnss的单频精密单点定位算法的时间频率传递方法 Download PDF

Info

Publication number
CN111965673B
CN111965673B CN202010586214.7A CN202010586214A CN111965673B CN 111965673 B CN111965673 B CN 111965673B CN 202010586214 A CN202010586214 A CN 202010586214A CN 111965673 B CN111965673 B CN 111965673B
Authority
CN
China
Prior art keywords
receiver
frequency
clock
satellite
representing
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
Application number
CN202010586214.7A
Other languages
English (en)
Other versions
CN111965673A (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.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
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 Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202010586214.7A priority Critical patent/CN111965673B/zh
Publication of CN111965673A publication Critical patent/CN111965673A/zh
Application granted granted Critical
Publication of CN111965673B publication Critical patent/CN111965673B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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

Abstract

基于多GNSS的单频精密单点定位算法的时间频率传递方法，包括下列步骤：接收机A接有原子钟A；接收机A通过天线接收到多个卫星系统的单频观测数据；进行钟差参数估计，解算得到本地的钟差参数ClockA‑参考时间；接收机B接有原子钟B；接收机B通过天线接收到多个卫星系统的单频观测数据；基于单频观测数据与精密产品数据，进行钟差参数估计；对两地得到的钟差参数ClockA、钟差参数ClockB进行统一，对两地的接收机进行仪器测量的方式或者使用接收机给出的标校值，完成两地的时间频率传递；ClockA‑ClockB＝(ClockA‑参考时间)‑(ClockB‑参考时间)。该方法能够实现将多GNSS精密单点定位算法应用于高精度的时间频率传递，可保证高精度时间频率传递的鲁棒性和连续性。

Description

基于多GNSS的单频精密单点定位算法的时间频率传递方法
技术领域
本发明涉及基于导航卫星系统单频数据的高精度时间频率传递方法，具体涉及一种基于多个全球定位导航系统(GNSS，包含GPS、BDS、Galileo、GLONASS、QZSS)的单频数据精密单点定位算法的时间频率传递方法。
背景技术
当前高精度的时间频率传递主要有共视、全视、精密单点定位(PPP)方法与卫星双向比对(TWSTT)方法。PPP方法，传统上使用消电离层组合的方法，需要双频数据，而由于单频接收机成本低廉，在当前卫星导航定位，时间传递领域拥有大量用户，因此，使用廉价的单频接收机，去实现高精度单机时间频率传递，恰好能满足低成本高精度单机用户的要求。在复杂多变的外部环境无法保证双频数据的准确性和连续性的条件下，单频算法能够使双频接收机维持使用。而保证时间传递的连续性与准确性是时间频率传递在实际应用中的急需解决的瓶颈问题。开展基于多GNSS单频伪距与载波相位观测值的精密单点定位算法用于时间频率传递领域的研究，将显著提升基于导航卫星的高精度时间频率传递的性能与鲁棒性，可更好地服务于国防与民生等应用领域。
发明内容
为克服现有技术的不足，本发明提供一种基于多GNSS的单频伪距载波相位观测值的两地时间频率传递模型，该模型能够实现将多GNSS精密单点定位算法应用于高精度的时间频率传递，可保证高精度时间频率传递的鲁棒性和连续性。
本发明基于多GNSS的单频精密单点定位算法的时间频率传递方法，包括下列步骤：
步骤1：第一接收机A接有高精度的第一外接原子钟A，第一外接原子钟A为第一接收机A提供1pps/10M的高精度时间频率信号；
步骤2：第一接收机A通过天线接收到多个卫星系统的单频点伪距、载波相位观测数据，将二者统称为第一单频观测数据；
步骤3：基于第一单频观测数据与精密产品数据，进行钟差参数估计，解算得到本地的第一钟差参数ClockA-参考时间，也就是本地的第一接收机A的时钟与使用的精密产品参考时间尺度的差值；
步骤4：第二接收机B接有高精度的第二外接原子钟B，第二外接原子钟B为第二接收机B提供1pps/10M的高精度时间频率信号；
步骤5：第二外接接收机B通过天线接收到多个卫星系统的单频点伪距、载波相位观测数据，将二者统称为第二单频观测数据；
步骤6：基于第二单频观测数据与精密产品数据，进行钟差参数估计，解算得到本地的第二钟差参数ClockB-参考时间，也就是本地的第二接收机B的时钟与使用的精密产品参考时间尺度的差值；
步骤7：对两地得到的第一钟差参数ClockA、第二钟差参数ClockB进行统一，同时对两地的接收机进行仪器测量的方式或者使用接收机给出的标校值，进而完成两地的时间频率传递，得到相应的应用；
ClockA-ClockB＝(ClockA-参考时间)-(ClockB-参考时间)。
在本发明的一个实施例中，钟差参数估计具体包括下列步骤：
步骤1：首先从外接有高精度原子钟的GNSS接收机获取到单频点伪距、载波相位观测数据；
步骤2：根据是否需要电离层产品选择具体的处理策略，如果选择需要电离层产品，那么根据实际精度的要求，从外部获取电离层产品；
步骤3：如果不需要获取电离层产品，或者需要获取，在获取电离层产品成功后，同时获取相应时间的精密轨道、精密钟差、差分码偏差DCB产品；
步骤4：上述所有数据准备好后，开始逐一历元进行处理，将上述数据代入以下载波相位与伪距观测方程中，如果有电离层产品，将电离层产品作为已知参数代入；
以下方程为基础，通过模型修正、外部产品以及参数估计手段，对与卫星有关的误差，包括卫星钟差、卫星轨道误差、相对论效应、与接收机和测站有关的误差，包括接收机钟差、接收机天线相位误差、地球潮汐、地球自转、与信号传播有关的误差，包括对流层延迟误差、电离层延迟误差和多路径效应，进行处理；
Figure GDA0004236099920000021
其中，P与Φ分别代表伪距与载波相位观测值，所有的上角标s代表卫星号，所有的下角标r代表接收机，所有的下角标j代表频点；所以
Figure GDA0004236099920000022
与/>
Figure GDA0004236099920000023
分别代表j频点数据，在卫星s与接收机r之间的伪距与载波相位观测值；
Figure GDA0004236099920000024
代表接收机r与卫星s的站星几何距离，c代表光速，dtr代表接收机钟差，dts代表卫星s的钟差，Zw代表湿的天顶方向对流层延迟，w代表“湿”，MFw(e)代表投影函数，本公式中
Figure GDA0004236099920000031
代表j频点的基于1频点的电离层延迟系数，/>
Figure GDA0004236099920000032
代表第一频点卫星S到接收机的电离层延迟，/>
Figure GDA0004236099920000033
与/>
Figure GDA0004236099920000034
分别代表j频点伪距在卫星与接收机端的硬件延迟，/>
Figure GDA0004236099920000035
与/>
Figure GDA0004236099920000036
分别代表j频点载波相位在卫星与接收机端的硬件延迟，/>
Figure GDA0004236099920000037
代表j频点信号，卫星s到r接收机的模糊度参数，
Figure GDA0004236099920000038
与/>
Figure GDA0004236099920000039
分别代表j频点信号，卫星s到r接收机伪距观测值P与载波相位观测值Φ的残差，也即函数模型对所有误差没有量化的部分；
步骤5：如果电离层
Figure GDA00042360999200000310
参数未知，无法从外部获取已知的电离层值/>
Figure GDA00042360999200000311
使用半和模型进行消电离层参数，由公式(6)推导得到如下公式：
Figure GDA00042360999200000312
其中，这里的L对应公式(1)的Φ，
Figure GDA00042360999200000313
代表j频点数据，在卫星s与接收机r之间的载波相位观测值，/>
Figure GDA00042360999200000314
为线性化方程得到的参数，而由于使用的是基于P1/P2码得到的卫星钟差产品，所以/>
Figure GDA00042360999200000315
代表的是对于差分码偏差进行修正项，其中/>
Figure GDA00042360999200000316
为DCB修正系数，使用不同的频点以及码组合将需要不同的修正系数，/>
Figure GDA00042360999200000317
代表基于P1/P2码的差分码偏差产品；X代表线性化待估参数，其中包含钟差参数；/>
Figure GDA00042360999200000318
代表消电离层组合的卫星钟差产品；
Figure GDA00042360999200000319
的方程为整合得到的函数观测方程；/>
Figure GDA00042360999200000320
为整合得到的延迟参数,/>
Figure GDA00042360999200000321
为待估计参数；
通过半和模型，对电离层参数
Figure GDA00042360999200000322
进行消除，除了电离层的处理，再获取相应的精密轨道、精密钟差、DCB的产品代入公式(7)；
步骤6：对准备好单个历元数据，进行基本的质量控制，包括对粗差的剔除、周跳的探测；
步骤7：进行扩展卡尔曼滤波处理，对钟差参数dtr进行估计，得到相应的观测数据残差
Figure GDA00042360999200000323
估计参数公式如下所示：
时间更新方程：
Figure GDA0004236099920000041
其中n维向量
Figure GDA0004236099920000042
为k时刻的系统状态变量，n维向量/>
Figure GDA0004236099920000043
是k-1时刻的系统状态变量；A是状态转移矩阵或者过程增益矩阵，是n×n阶方阵，其上标T为“转置”，它将k-1时刻状态和当前的k时刻状态联系起来；B是可选的控制输入的增益，uk-1为驱动噪声，在大多数实际情况下并没有控制增益，所以Buk-1这一项为零；/>
Figure GDA0004236099920000044
为X参数的协方差矩阵，Q为状态方程噪声；
状态更新方程：
Figure GDA0004236099920000045
其中Kk表示卡尔曼增益，是滤波的中间结果，H表示量测矩阵，是m×n阶矩阵，它把m维测量值转换到n维与状态变量相对应；Zk表示测量值，是m阶向量；R为观测噪声矩阵，
Figure GDA0004236099920000046
为更新的待估计参数，Pk为更新的待估计参数协方差矩阵；
步骤8：对应卫星残差
Figure GDA0004236099920000047
过大，则剔除该卫星，回到步骤7进行重新估计，直到所有数据满足阈值条件为止；阈值的使用根据实际处理过程经验值；
步骤9：保存相关结果钟差参数dtr，估计的同时，得到坐标、模糊度、系统间偏差参数，保存该历元结果的同时，回到步骤6，进行下一个历元结果的处理。
本发明的有益效果是：
第一，可以对多GNSS与单频伪距与载波相位用于高精度时间频率传递数据处理模型进行统一。本发明将多GNSS和单频两种数据处理模型进行有机统一，在复杂多变的外部环境和需求情况下，无法实现双频解算时，可更好地实现基于GNSS单频数据获得高度鲁棒性与稳定性。基于统一的数据处理模型，可实现多系统和单频两种算法的有机统一。
第二，为高精度时间频率传递用户提供选择多样性。对于电离层的处理，实现多种模式切换，以适应复杂多变的环境，更好地适应用户需求。
附图说明
图1是基于导航卫星的高精度时间频率传递方法实施框图；
图2是基于GNSS单频伪距与载波相位的钟差解算方法流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。
本发明将多GNSS和单频高精度时间频率传递方法进行有机融合，建立高鲁棒性的统一模型，用户可以获取高精度时间频率的稳定性和连续性。本发明可将北斗卫星导航系统(BDS)应用于时间频率行业，更好地服务于国防、军事、电力、通信、银行、金融等行业。
本发明的技术方案主要包括三个核心技术：
(1)各项误差的处理
通过模型修正、外部产品以及参数估计等手段，对与卫星有关的误差(卫星钟差、卫星轨道误差、相对论效应)、与接收机和测站有关的误差(接收机钟差、接收机天线相位误差、地球潮汐、地球自转等)、与信号传播有关的误差(对流层延迟误差、电离层延迟误差和多路径效应)进行处理。这些处理均为本领域技术人员熟知，不再累述。
(2)单频点电离层误差的处理
第1种方法是采用电离层模型改正(参见各个卫星导航系统的导航文件ICD接口文件)，不同系统的广播星历里面包含不同的电离层参数模型；第2种方法是根据电离层引起的群延迟和相延迟大小相等、方向相反这一特点，采用码和相位的半合组合进行修正(参见高扬论文Gao Y,Shen X.A New Method for Carrier-Phase-Based Precise PointPositioning[J].Navigation,2002,49(2):p.109-116)；第3种方法是将电离层延迟当作未知参数，类似与本专利要估计的钟差参数同时进行估计。
(3)建立函数模型进行钟差参数估计
根据选择的电离层处理方式不一样，建立不同的函数模型，以下对具体涉及的函数模型进行解释。首先对于接收机观测到的伪距与载波相位可以直接给出伪距与载波相位观测方程，如下所示：
Figure GDA0004236099920000051
其中，P与Φ分别代表伪距与载波相位观测值，所有的上角标s代表卫星号，所有的下角标r代表接收机，所有的下角标j代表频点。所以
Figure GDA0004236099920000061
与/>
Figure GDA0004236099920000062
分别代表j频点数据，在卫星s与接收机r之间的伪距与载波相位观测值。
Figure GDA0004236099920000063
代表接收机r与卫星s的站星几何距离，c代表光速，dtr代表接收机钟差，dts代表卫星s的钟差，Zw代表湿的天顶方向对流层延迟，w代表湿(wet)的意思，MFw(e)代表投影函数，/>
Figure GDA0004236099920000064
代表j频点的基于1频点的电离层延迟系数，/>
Figure GDA0004236099920000065
代表第一频点卫星S到接收机的电离层延迟，/>
Figure GDA0004236099920000066
与/>
Figure GDA0004236099920000067
分别代表j频点伪距在卫星与接收机端的硬件延迟，/>
Figure GDA0004236099920000068
与/>
Figure GDA0004236099920000069
分别代表j频点载波相位在卫星与接收机端的硬件延迟，/>
Figure GDA00042360999200000610
代表j频点信号，卫星s到r接收机的模糊度参数，/>
Figure GDA00042360999200000611
与/>
Figure GDA00042360999200000612
分别代表j频点信号，卫星s到r接收机伪距观测值P与载波相位观测值Φ的残差，也即函数模型(公式)对所有误差没有量化的部分。
由于选择单频数据，这里只以GNSS的第一频点数据为例，由于卡尔曼滤波只能是对线性方程进行，所以对上两式进行线性化处理得到：
Figure GDA00042360999200000613
其中，这里的L对应上面公式(1)的Φ，
Figure GDA00042360999200000620
代表j频点数据，在卫星s与接收机r之间的载波相位观测值，/>
Figure GDA00042360999200000614
为线性化方程得到的参数，而由于我们使用的是基于P1/P2码得到的卫星钟差产品，所以这里/>
Figure GDA00042360999200000615
代表的是对于差分码偏差进行修正项，其中/>
Figure GDA00042360999200000616
为DCB修正系数，使用不同的频点以及码组合将需要不同的修正系数，/>
Figure GDA00042360999200000617
代表基于P1/P2码的差分码偏差产品。X代表线性化待估参数，其中包含钟差参数。/>
Figure GDA00042360999200000618
代表消电离层组合的卫星钟差产品。
如果没有电离层产品，对电离层使用半和模型，半和模型为基于伪距观测值与载波相位观测值在穿越电离层时延迟大小相等，符号相反特性进行整合，整合上述两个公式得到如下所示：
Figure GDA00042360999200000619
其中，
Figure GDA0004236099920000071
的方程为整合得到的函数观测方程。/>
Figure GDA0004236099920000072
为整合得到的延迟参数,/>
Figure GDA0004236099920000073
为待估计参数。
如果是多系统(BDS,GPS,Galileo,GLONASS,QZSS)，则多估计一个ISB(系统间偏差)参数，得到待估参数，有电离层估计与无电离层估计参数序列分别为：
X1＝[x,dtr,ISB,Zw，I,Nb,Ns] (4)
X2＝[x,dtr,ISB,Zw,Nb,Ns] (5)
其中，x代表待估计位置参数，dtr代表待估计参数接收机钟差，ISB为多系统条件下系统参数，Zw为湿的天顶对流层延迟，I为电离层延迟，其中b与s分别代表基准系统与增加的系统，Nb为基准系统模糊度，Ns为增加的系统模糊度。
最后，根据(2)或(3)公式，利用扩展卡尔曼滤波估算得到本地钟差数据dtr，同时对两地的接收机硬件延迟进行标校，进而完成对两地的钟差进行统一与比对，实现高精度的时间频率传递。这个过程可以是事后，也可以是实时模式，取决于用户的具体需求。
图1描述本发明基于多GNSS的单频精密单点定位算法的时间频率传递方法的总体实施方案流程，具体如下：
步骤1：接收机A接有高精度的外接原子钟A，原子钟A为接收机A提供1pps/10M的高精度时间频率信号。
步骤2：接收机A通过天线接收到多个卫星系统的单频点伪距，载波相位观测数据。
步骤3：基于单频观测数据与精密产品数据，进行钟差参数估计，详细处理方法如图2所示，精密产品为例如来源于欧洲定轨中心CODE、德国地学研究中心GFZ等的精密轨道、精密钟差(本领域技术人员熟知)，解算得到本地的钟差参数ClockA-参考时间，也就是本地的接收机钟与使用的精密产品参考时间尺度的差值。
步骤4：接收机B接有高精度的外接原子钟B，原子钟B为接收机B提供1pps/10M的高精度时间频率信号。
步骤5：接收机B通过天线接收到多个卫星系统的单频点伪距，载波相位观测数据。
步骤6：基于单频观测数据与精密产品数据，进行钟差参数估计，详细处理方法如图2所示，单频观测数据与精密产品，解算得到本地的钟差参数ClockB-参考时间，也就是本地的接收机钟B与使用的精密产品参考时间尺度的差值。
步骤7：对两地得到的钟差参数(ClockA,ClockB)进行统一，同时对两地的接收机进行仪器测量的方式或者使用接收机给出的标校值，进而完成两地的时间频率传递，得到相应的应用。
ClockA-ClockB＝(ClockA-参考时间)-(ClockB-参考时间)
图2描述图1中子流程，描述了其中钟差参数估计的详细算法步骤。
步骤1：首先从GNSS接收机(外接有高精度原子钟)获取到单频点的伪距，载波相位数据。
步骤2：根据是否需要电离层产品选择具体的处理策略，如果选择需要电离层产品，那么从外部获取，根据实际精度的要求，有很多种方式。例如国际上CODE等机构获取ion产品精度较高，也可以从广播星历中提取参数进行模型改正，精度较低。
步骤3：如果不需要获取电离层产品，或者需要获取，获取电离层产品成功后，同时获取相应时的精密轨道，精密钟差，差分码偏差(DCB)等产品，如果上述过程是实时模式，就实时获取，如果是事后模式，可采用事后下载方式获取。
步骤4：将所有数据准备好以后，开始逐一历元进行处理，例如1:00:00估计一次，1:00:05估计一次，每隔5秒估计一次，逐一历元进行估计，实际估计间隔根据观测数据的采样间隔。可以代入以下载波相位与伪距观测方程中，如果有电离层产品，将电离层产品作为已知参数代入。
以下方程为基础，通过模型修正、外部产品以及参数估计等手段，对与卫星有关的误差(卫星钟差、卫星轨道误差、相对论效应)、与接收机和测站有关的误差(接收机钟差、接收机天线相位误差、地球潮汐、地球自转等)、与信号传播有关的误差(对流层延迟误差、电离层延迟误差和多路径效应)进行处理。这些处理均为本领域技术人员熟知，不再累述。
Figure GDA0004236099920000081
步骤5：如果电离层
Figure GDA0004236099920000082
参数未知，无法从外部获取已知的电离层值/>
Figure GDA0004236099920000083
可以使用半和模型(参见高扬论文Gao Y,Shen X.A New Method for Carrier-Phase-Based PrecisePoint Positioning[J].Navigation,2002,49(2):p.109-116)进行消电离层参数，由公式(6)推导得到如下公式：
Figure GDA0004236099920000084
通过半和模型，对电离层参数
Figure GDA0004236099920000099
进行了消除，除了电离层的处理，再获取相应的精密轨道，精密钟差，DCB等产品代入公式(7)。此时所有公式中已知数据准备好，可以进下一步。
步骤6：对准备好单个历元数据，进行基本的质量控制，包括对粗差的剔除，周跳的探测等等。该步骤的具体实施方法为本领域技术人员熟知，不再累述。
步骤7：进行扩展卡尔曼滤波处理，对钟差参数dtr进行估计，得到相应的观测数据残差
Figure GDA0004236099920000091
估计参数公式如下所示：
时间更新方程：
Figure GDA0004236099920000092
其中n维向量
Figure GDA0004236099920000093
为k时刻的系统状态变量，n维向量/>
Figure GDA0004236099920000094
是k-1时刻的系统状态变量。A是状态转移矩阵或者过程增益矩阵，是n×n阶方阵，其上标T为“转置”，它将k-1时刻状态和当前的k时刻状态联系起来。B是可选的控制输入的增益，uk-1为驱动噪声，在大多数实际情况下并没有控制增益，所以Buk-1这一项为零。/>
Figure GDA0004236099920000095
为X参数的协方差矩阵，Q为状态方程噪声。
状态更新方程：
Figure GDA0004236099920000096
其中Kk表示卡尔曼增益，是滤波的中间结果，H表示量测矩阵，是m×n阶矩阵，它把m维测量值转换到n维与状态变量相对应。Zk表示测量值，是m阶向量。R为观测噪声矩阵，
Figure GDA0004236099920000097
为更新的待估计参数，Pk为更新的待估计参数协方差矩阵。
步骤8：对应卫星残差
Figure GDA0004236099920000098
过大，则剔除该卫星，回到步骤7进行重新估计，直到所有数据满足阈值条件为止。阈值的使用根据实际处理过程经验值。
步骤9，保存相关结果钟差参数dtr，估计的同时，我们会得到坐标，模糊度，系统间偏差等参数，在此专利中我们不予使用，保存该历元结果的同时，我们回到步骤6，进行下一个历元结果的处理。
自全球卫星导航系统(GNSS)问世以来,由于其具有全球性、全天候、高精度、高效率、保密性强等一系列特点,现已被广泛应用于卫星导航、测量定位、形变监测、大气探测，高精度时间频率传递等领域。目前,GNSS精密导航与定位与授时，一般都采用双频接收机,其主要原因是为了消除电离层误差的影响。在利用GNSS技术开展区域变形监测、大气探测、时间频率传递、时间同步等研究中,由于需要采集高时空分辨率的信息,不得不布设大量的GNSS接收机。如全部采用GNSS双频接收机进行施测,其成本无疑将非常昂贵,这必然会极大限制GNSS技术在这些领域未来的发展和应用。正因为如此,如何消除电离层误差的影响,利用相对廉价的单频接收机实现大范围的高精度时间频率传递测试与使用，正是本发明解决的问题。现有技术尚未出现使用过单频方法进行时间传递的解决方案，本发明针对单频解决方案，在有无外部辅助电离层产品的情况下，提出对应解决方案。

Claims (1)

1.基于多GNSS的单频精密单点定位算法的时间频率传递方法，其特征在于，包括下列步骤：
步骤1：第一接收机A接有高精度的第一外接原子钟A，第一外接原子钟A为第一接收机A提供1pps/10M的高精度时间频率信号；
步骤2：第一接收机A通过天线接收到多个卫星系统的单频点伪距、载波相位观测数据，将二者统称为第一单频观测数据；
步骤3：基于第一单频观测数据与精密产品数据，进行钟差参数估计，解算得到本地的第一钟差参数ClockA-参考时间，也就是本地的第一接收机A的时钟与使用的精密产品参考时间尺度的差值；
其中钟差参数估计具体包括下列步骤：
步骤31：首先从外接有高精度原子钟的GNSS接收机获取到单频点伪距、载波相位观测数据；
步骤32：根据是否需要电离层产品选择具体的处理策略，如果选择需要电离层产品，那么根据实际精度的要求，从外部获取电离层产品；
步骤33：如果不需要获取电离层产品，或者需要获取，在获取电离层产品成功后，同时获取相应时的精密轨道、精密钟差、差分码偏差DCB产品；
步骤34：上述所有数据准备好后，开始逐一历元进行处理，将上述数据代入以下载波相位与伪距观测方程中，如果有电离层产品，将电离层产品作为已知参数代入；
以下方程为基础，通过模型修正、外部产品以及参数估计手段，对与卫星有关的误差，包括卫星钟差、卫星轨道误差、相对论效应、与接收机和测站有关的误差，包括接收机钟差、接收机天线相位误差、地球潮汐、地球自转、与信号传播有关的误差，包括对流层延迟误差、电离层延迟误差和多路径效应，进行处理；
Figure FDA0004236099910000011
其中，P与Φ分别代表伪距与载波相位观测值，所有的上角标s代表卫星号，所有的下角标r代表接收机，所有的下角标j代表频点；所以
Figure FDA0004236099910000012
与/>
Figure FDA0004236099910000013
分别代表j频点数据，在卫星s与接收机r之间的伪距与载波相位观测值；
Figure FDA0004236099910000014
代表接收机r与卫星s的站星几何距离，c代表光速，dtr代表接收机钟差，dts代表卫星s的钟差，Zw代表湿的天顶方向对流层延迟，w代表湿，MFw(e)代表投影函数，/>
Figure FDA0004236099910000021
代表j频点的基于1频点的电离层延迟系数，/>
Figure FDA0004236099910000022
代表第一频点卫星S到接收机的电离层延迟，/>
Figure FDA0004236099910000023
与/>
Figure FDA0004236099910000024
分别代表j频点伪距在卫星与接收机端的硬件延迟，/>
Figure FDA0004236099910000025
与/>
Figure FDA0004236099910000026
分别代表j频点载波相位在卫星与接收机端的硬件延迟，/>
Figure FDA0004236099910000027
代表j频点信号，卫星s到r接收机的模糊度参数，/>
Figure FDA0004236099910000028
与
Figure FDA0004236099910000029
分别代表j频点信号，卫星s到r接收机伪距观测值P与载波相位观测值Φ的残差，也即函数模型对所有误差没有量化的部分；
步骤35：如果电离层
Figure FDA00042360999100000210
参数未知，无法从外部获取已知的电离层值/>
Figure FDA00042360999100000211
使用半和模型进行消电离层参数，由公式(6)推导得到如下公式：
Figure FDA00042360999100000212
其中，这里的L对应公式(1)的Φ，
Figure FDA00042360999100000213
代表j频点数据，在卫星s与接收机r之间的载波相位观测值，/>
Figure FDA00042360999100000214
为线性化方程得到的参数，而由于使用的是基于P1/P2码得到的卫星钟差产品，所以/>
Figure FDA00042360999100000215
代表的是对于差分码偏差进行修正项，其中/>
Figure FDA00042360999100000216
为DCB修正系数，使用不同的频点以及码组合将需要不同的修正系数，/>
Figure FDA00042360999100000217
代表基于P1/P2码的差分码偏差产品；X代表线性化待估参数，其中包含钟差参数；/>
Figure FDA00042360999100000218
代表消电离层组合的卫星钟差产品；
Figure FDA00042360999100000219
的方程为整合得到的函数观测方程；/>
Figure FDA00042360999100000220
为整合得到的延迟参数,/>
Figure FDA00042360999100000221
为待估计参数；
通过半和模型，对电离层参数
Figure FDA00042360999100000222
进行消除，除了电离层的处理，再获取相应的精密轨道、精密钟差、DCB的产品代入公式(7)；
步骤36：对准备好单个历元数据，进行基本的质量控制，包括对粗差的剔除、周跳的探测；
步骤37：进行扩展卡尔曼滤波处理，对钟差参数dtr进行估计，得到相应的观测数据残差
Figure FDA00042360999100000223
估计参数公式如下所示：
时间更新方程：
Figure FDA0004236099910000031
其中，其中n维向量
Figure FDA0004236099910000032
为k时刻的系统状态变量，n维向量/>
Figure FDA0004236099910000033
是k-1时刻的系统状态变量；A是状态转移矩阵或者过程增益矩阵，是n×n阶方阵，其上标T为转置，它将k-1时刻状态和当前的k时刻状态联系起来；B是控制输入的增益，uk-1为驱动噪声，并没有控制增益，所以Buk-1这一项为零；/>
Figure FDA0004236099910000034
为X参数的协方差矩阵，Q为状态方程噪声；
状态更新方程：
Figure FDA0004236099910000035
其中Kk表示卡尔曼增益，是滤波的中间结果，H表示量测矩阵，是m×n阶矩阵，它把m维测量值转换到n维与状态变量相对应；Zk表示测量值，是m阶向量；R为观测噪声矩阵，
Figure FDA0004236099910000036
为更新的待估计参数，Pk为更新的待估计参数协方差矩阵；
步骤38：对应卫星残差
Figure FDA0004236099910000037
过大，则剔除该卫星，回到步骤7进行重新估计，直到所有数据满足阈值条件为止；阈值的使用根据实际处理过程经验值；
步骤39：保存相关结果钟差参数dtr，估计的同时，得到坐标、模糊度、系统间偏差参数，保存该历元结果的同时，回到步骤6，进行下一个历元结果的处理；
步骤4：第二接收机B接有高精度的第二外接原子钟B，第二外接原子钟B为第二接收机B提供1pps/10M的高精度时间频率信号；
步骤5：第二外接接收机B通过天线接收到多个卫星系统的单频点伪距、载波相位观测数据，将二者统称为第二单频观测数据；
步骤6：基于第二单频观测数据与精密产品数据，进行钟差参数估计，解算得到本地的第二钟差参数ClockB-参考时间，也就是本地的第二接收机B的时钟与使用的精密产品参考时间尺度的差值；
步骤7：对两地得到的第一钟差参数ClockA、第二钟差参数ClockB进行统一，同时对两地的接收机进行仪器测量的方式或者使用接收机给出的标校值，进而完成两地的时间频率传递，得到相应的应用；
ClockA-ClockB＝(ClockA-参考时间)-(ClockB-参考时间)。
CN202010586214.7A 2020-06-24 2020-06-24 基于多gnss的单频精密单点定位算法的时间频率传递方法 Active CN111965673B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010586214.7A CN111965673B (zh) 2020-06-24 2020-06-24 基于多gnss的单频精密单点定位算法的时间频率传递方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010586214.7A CN111965673B (zh) 2020-06-24 2020-06-24 基于多gnss的单频精密单点定位算法的时间频率传递方法

Publications (2)

Publication Number Publication Date
CN111965673A CN111965673A (zh) 2020-11-20
CN111965673B true CN111965673B (zh) 2023-06-20

Family

ID=73360636

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010586214.7A Active CN111965673B (zh) 2020-06-24 2020-06-24 基于多gnss的单频精密单点定位算法的时间频率传递方法

Country Status (1)

Country Link
CN (1) CN111965673B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112924992B (zh) * 2021-01-25 2022-11-04 中国科学院国家授时中心 一种geo轨道精度评估方法、装置、电子设备及存储介质
CN113325446B (zh) * 2021-06-01 2023-08-08 中国科学院国家授时中心 一种多模共频gnss载波相位时间传递方法及系统
CN114923519A (zh) * 2022-05-16 2022-08-19 中咨数据有限公司 工程地质监测系统及多期航拍测绘的三维地质建模方法
CN115407371B (zh) * 2022-09-02 2023-08-15 中国人民解放军国防科技大学 基于PPP-B2b的实时高精度时间传递方法及装置
CN116299617B (zh) * 2023-02-16 2023-10-31 北京航空航天大学 一种gnss星地原子钟实时动态模型的建立方法
CN116318509B (zh) * 2023-03-01 2024-04-09 中国人民解放军国防科技大学 基于模糊度固定残差后验定权的ppp时间频率传递方法
CN116299618B (zh) * 2023-03-24 2024-03-19 中国科学院精密测量科学与技术创新研究院 基于ppp解算参数的载波相位卫星共视时间传递方法
CN116540282B (zh) * 2023-05-08 2023-09-29 北京交通大学 一种基于惯性传感器辅助的gnss的移动载体精准时间传递方法
CN116299585B (zh) * 2023-05-15 2023-09-08 中国科学院国家授时中心 一种顾及历元间差分信息的gnss载波相位时间传递方法
CN117270004B (zh) * 2023-11-22 2024-01-26 中交路桥科技有限公司 基于卫星定位的工程测绘方法、设备、系统及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101403790A (zh) * 2008-11-13 2009-04-08 浙江师范大学 单频gps接收机的精密单点定位方法
CN104199058A (zh) * 2014-09-18 2014-12-10 中国人民解放军国防科学技术大学 基于Kalman滤波器实时预测值调整时间尺度算法
CN107356947A (zh) * 2017-05-31 2017-11-17 中国科学院测量与地球物理研究所 基于单频导航卫星数据确定卫星差分伪距偏差的方法
CN108169774A (zh) * 2017-12-26 2018-06-15 北方信息控制研究院集团有限公司 支持rtppp和rtk的多模gnss单频周跳探测与修复方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101403790A (zh) * 2008-11-13 2009-04-08 浙江师范大学 单频gps接收机的精密单点定位方法
CN104199058A (zh) * 2014-09-18 2014-12-10 中国人民解放军国防科学技术大学 基于Kalman滤波器实时预测值调整时间尺度算法
CN107356947A (zh) * 2017-05-31 2017-11-17 中国科学院测量与地球物理研究所 基于单频导航卫星数据确定卫星差分伪距偏差的方法
CN108169774A (zh) * 2017-12-26 2018-06-15 北方信息控制研究院集团有限公司 支持rtppp和rtk的多模gnss单频周跳探测与修复方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
基于初始运行阶段Galileo PPP国际时间传递;欧阳明俊;杨旭海;孙保琪;葛玉龙;;时间频率学报(第02期);80-86 *
基于无电离层模型的动态精密单点定位;李黎;苏强;叶和军;;四川兵工学报(第01期);51-54 *
多系统钟差预报方法研究及对实时精密单点定位与时间传递的影响;李成钢;史小雨;欧阳明俊;左翔;陈伟;;测绘(第05期);195-200,205 *
涂锐 ; 黄观文 ; 张勤 ; 王利 ; .利用单基准站改正信息和电离层参数估计的单频PPP算法.武汉大学学报(信息科学版).2012,37(2),170-173,182. *

Also Published As

Publication number Publication date
CN111965673A (zh) 2020-11-20

Similar Documents

Publication Publication Date Title
CN111965673B (zh) 基于多gnss的单频精密单点定位算法的时间频率传递方法
CN109709591B (zh) 一种面向智能终端的gnss高精度定位方法
CN109581452B (zh) 一种gnss参考站载波相位整周模糊度解算方法
EP2064567B1 (en) Relative positioning
CN1864078B (zh) 使用三个gps频率消除载波相位整数模糊度的方法
CN111045034B (zh) 基于广播星历的gnss多系统实时精密时间传递方法及系统
CN107710017A (zh) 用于在实时运动模式和相对定位模式之间切换的卫星导航接收器及方法
EP3502747B1 (en) Advanced navigation satellite system positioning method and system using seeding information
CN107607971B (zh) 基于gnss共视时间比对算法的时间频率传递方法及接收机
Abd Rabbou et al. Performance analysis of precise point positioning using multi-constellation GNSS: GPS, GLONASS, Galileo and BeiDou
CN113325446B (zh) 一种多模共频gnss载波相位时间传递方法及系统
CN114966760B (zh) 一种电离层加权的非差非组合ppp-rtk技术实现方法
CN117836669A (zh) 使用参考站载波相位测量的dgnss
Li et al. Review of PPP–RTK: Achievements, challenges, and opportunities
CN112505735B (zh) 对终端进行定位的方法、装置和存储介质
US10830898B2 (en) Method and apparatus applicable to positioning in NLOS environment
CN114935770B (zh) 一种多历元加快精密单点定位收敛速度的方法及装置
CN110068848B (zh) 一种高性能rtk处理技术方法
WO2001022111A1 (en) Global navigation satellite systems and methods
CN111983641A (zh) 一种用于实时生成北斗星基增强系统完好性参数的方法
CN116540303B (zh) 一种历元间差分观测方程建立方法和地震位移计算方法
JP4723801B2 (ja) 相対測位装置
CN104991265A (zh) 一种北斗卫星导航系统用户统一性定位方法
CN113109840A (zh) 一种基于gnss接收机的电离层tec实时测量方法
Hong et al. Analysis of dual-frequency solution method for single-frequency precise point positioning based on SEID model for GPS and BDS

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