CN116299585B - 一种顾及历元间差分信息的gnss载波相位时间传递方法 - Google Patents

一种顾及历元间差分信息的gnss载波相位时间传递方法 Download PDF

Info

Publication number
CN116299585B
CN116299585B CN202310542390.4A CN202310542390A CN116299585B CN 116299585 B CN116299585 B CN 116299585B CN 202310542390 A CN202310542390 A CN 202310542390A CN 116299585 B CN116299585 B CN 116299585B
Authority
CN
China
Prior art keywords
carrier phase
time transfer
representing
epoch
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.)
Active
Application number
CN202310542390.4A
Other languages
English (en)
Other versions
CN116299585A (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.)
National Time Service Center of CAS
Original Assignee
National Time Service Center of CAS
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 National Time Service Center of CAS filed Critical National Time Service Center of CAS
Priority to CN202310542390.4A priority Critical patent/CN116299585B/zh
Publication of CN116299585A publication Critical patent/CN116299585A/zh
Application granted granted Critical
Publication of CN116299585B publication Critical patent/CN116299585B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/258—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 the satellite constellation, e.g. almanac, ephemeris data, lists of satellites in view
    • 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
    • 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/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/37—Hardware or software details of the signal processing chain
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90—Details of database functions independent of the retrieved data types
    • G06F16/906—Clustering; Classification
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • 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)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明涉及载波相位时间传递领域，具体涉及一种顾及历元间差分信息的GNSS载波相位时间传递方法。实时获取两个时间传递站的伪距观测值、载波相位观测值以及精密卫星轨道和钟差产品数据；对伪距观测值以及载波相位观测值进行预处理，并构建载波相位时间传递模型；建立历元间的差分解算模型，基于载波相位时间传递模型以及历元间的差分解算模型进行钟差解算，获取当前历元两个时间传递站的GNSS接收机钟差；根据当前历元两个时间传递站的GNSS接收机钟差获取时间传递量。本发明通过历元间差分信息的获取对时间传递模型进行有效增强，提升了接收机钟差的解算强度，并进一步增强了GNSS载波相位时间传递的性能。

Description

一种顾及历元间差分信息的GNSS载波相位时间传递方法
技术领域
本发明涉及载波相位时间传递领域，具体涉及一种顾及历元间差分信息的GNSS载波相位时间传递方法。
背景技术
利用全球卫星导航系统（GNSS）的远程时间传递技术作为一种集应用效率高、成本低，设备体积小，机动灵活、维护简单等优势为一体的时间空间时频传递技术，已广泛应用于卫星导航，电力、通信、交通、国防军工、科学研究等领域。随着高精度的时间频率标准的不断发展，尤其是远程时间传递技术作为与国家标准时间频率建立联系的重要手段，对于时间传递性能也提出了更高的要求。
然而，传统的GNSS载波相位时间传递技术仅仅利用了当前历元的伪距和载波相位观测量，受观测环境和接收机硬件设备的影响，在观测量中随着传递时长的增加会产生异常，造成时间传递的精度衰减，究其原因是传统的GNSS载波相位时间传递方法并未考虑到两类观测量在时间域内的相关性，损失了大量有益的观测信息，因此如何有效融合GNSS观测量在时间域内的有效信息，进一步提升GNSS时间传递的可靠性是时间频率领域关注的热点问题。
发明内容
为了解决现有技术中在观测量中随着传递时长的增加会产生异常，造成时间传递的精度衰减的问题，本发明提供一种顾及历元间差分信息的GNSS载波相位时间传递方法，该方案包括：实时获取两个时间传递站的伪距观测值、载波相位观测值以及精密卫星轨道和钟差产品数据；对伪距观测值以及载波相位观测值进行预处理，并构建载波相位时间传递模型；建立历元间的差分解算模型，基于载波相位时间传递模型以及历元间的差分解算模型进行钟差解算，获取当前历元两个时间传递站的GNSS接收机钟差；根据当前历元两个时间传递站的GNSS接收机钟差获取时间传递量。本发明通过历元间差分信息的获取对时间传递模型进行有效增强，提升了接收机钟差的解算强度，并进一步增强了GNSS载波相位时间传递的性能。
本发明采用如下技术方案，一种顾及历元间差分信息的GNSS载波相位时间传递方法，包括：
实时获取两个时间传递站每个历元的GNSS双频伪距观测值、载波相位观测值，以及对应的GNSS精密卫星轨道和钟差产品数据；
对所述双频伪距观测值以及载波相位观测值进行预处理，并根据预处理后的双频伪距观测值以及载波相位观测值构建载波相位时间传递模型；
当未发生周跳时，根据上一个历元与当前历元的双频伪距观测值以及载波相位观测值建立历元间的差分解算模型；
基于载波相位时间传递模型以及历元间的差分解算模型进行钟差解算，获取当前历元两个时间传递站的GNSS接收机钟差；
根据当前历元两个时间传递站的GNSS接收机钟差获取时间传递量。
进一步的，对所述双频伪距观测值以及载波相位观测值进行预处理的方法为：
对获取到的GNSS双频伪距观测值以及载波相位观测值进行数据分类、野值剔除以及潮汐和天线相位中心误差的修正。
进一步的，对所述双频伪距观测值以及载波相位观测值进行预处理，还包括：
对获取到的GNSS双频伪距观测值以及载波相位观测值进行周跳探测，并构建双频无电离层伪距和载波相位组合量，表示为：
；
其中，表示伪距的无电离层组合观测量，表示载波相位的无电离层组合观 测量，和表示两个双频伪距观测量，和为两个双频载波相位观测量，和分 别为两个频率的频点。
进一步的，建立历元间的差分解算模型为：
；
其中，为历元间一次差分运算符，k表示第k个历元，表示第k个历元的下 一个历元，i为卫星标识，r为接收机标识，表示卫星i与接收机r之间的距离，表示接 收机钟差，表示卫星钟差，表示伪距噪声，表示载波相位噪声，c为光速，表 示伪距的无电离层组合观测量，表示载波相位的无电离层组合观测量，表示电离层延 迟。
进一步的，基于载波相位时间传递模型以及历元间的差分解算模型进行钟差解算的方法为：
基于载波相位时间传递模型以及历元间的差分解算模型进行融合，构建历元间差分信息的接收机钟差解算模型；
所述接收机钟差解算模型包括函数模型和随机模型，函数模型为：
；
其中，表示伪距的无电离层组合观测量，表示载波相位的无电离层组合观测 量，表示卫星i与接收机r之间的距离，表示接收机钟差，表示卫星钟差，表示伪 距噪声，表示载波相位噪声，c为光速，为对流层延迟，为载波相位模糊度，为 历元间一次差分运算符，为第k个历元与第个历元之间的钟差差分量。
随机模型为：
；
其中，为时间传递模型的权重分配，/>为时间传递模型中对应的中误差，/>为历元间的差分解算模型的权重分配，/>为历元间的差分解算模型中对应的中误差。
进一步的，进行钟差解算时还包括：
对接收机钟差解算模型中的参数进行修正，包括：
对流层延迟按照模型改正并结合随机游走的随机过程进行模拟和估计；
载波相位模糊度在连续无周跳情况下作为常数估计；
接收机钟差作为高斯白噪声逐历元估计。
进一步的，根据当前历元两个时间传递站的GNSS接收机钟差获取时间传递量的方法为：
获取两个时间传递站当前历元的接收机钟差量之间的差值，并根据该差值与时间传递链路中固有延迟之和，得到对应两个时间传递站之间的时间传递量。
进一步的，根据预处理后的双频伪距观测值以及载波相位观测值构建载波相位时间传递模型，该模型具体为：
；
其中，表示伪距的无电离层组合观测量，/>表示载波相位的无电离层组合观测量，i为卫星标识，r为接收机标识，c为光速，/>表示接收机钟差，/>表示卫星钟差，/>表示伪距噪声，/>表示载波相位噪声，/>为对流层延迟，/>为载波相位模糊度。
本发明的有益效果是：本发明建立基于历元间差分的钟差解算模型，能够有效获取历元间的差分信息，同时通过对基于历元间的差分信息模型化，并与载波相位时间传递模型进行融合，将GNSS观测量在时间域内的相关信息融入到GNSS载波相位时间传递模型中，能够考虑到观测量在时间域内的相关性，利用这些相关性提升了接收机钟速参数的估计精度，为进一步提升GNSS载波相位时间传递性能奠定了基础，同时为提升时间传递的可靠性提供了理论支撑。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例的一种顾及历元间差分信息的GNSS载波相位时间传递方法流程示意图；
图2为本发明实施例中一种构建接收机钟差解算模型的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
首先利用GNSS观测数据，构建历元间差分模型，获得接收机钟差参数在历元间的变化率，然后结合传统GNSS载波相位时间传递模型，构建基于历元间差分增强的时间传递模型，增强接收机钟差参数的解算强度，进一步达到提升GNSS时间传递性能的目的。如图1所示，给出了本发明实施例的一种顾及历元间差分信息的GNSS载波相位时间传递方法流程示意图，包括：
101.实时获取两个时间传递站每个历元的GNSS双频伪距观测值、载波相位观测值，以及对应的GNSS精密卫星轨道和钟差产品数据；
102.对所述双频伪距观测值以及载波相位观测值进行预处理，并根据预处理后的双频伪距观测值以及载波相位观测值构建载波相位时间传递模型；
本发明对两个时间传递站的GNSS观测数据分别处理，包括测量完整性检查、异常值探测与剔除及周跳探测与标记，得到干净的观测数据。紧接着，通过对电离层误差、对流层、卫星及接收机相关的偏差，以及地球潮汐、海潮等误差的改正，完成模型观测方程生成前的各种准备。
具体的，对所述双频伪距观测值以及载波相位观测值进行预处理的方法为：
对获取到的GNSS双频伪距观测值以及载波相位观测值进行数据分类、野值剔除以及潮汐和天线相位中心误差的修正。
对所述双频伪距观测值以及载波相位观测值进行预处理，还包括：
对获取到的GNSS双频伪距观测值以及载波相位观测值进行周跳探测，并构建双频无电离层伪距和载波相位组合量，表示为：
；
其中，表示伪距的无电离层组合观测量，表示载波相位的无电离层组合观测 量，和表示两个双频伪距观测量，和为两个双频载波相位观测量，和分别为 两个频率的频点。
根据预处理后的双频伪距观测值以及载波相位观测值构建载波相位时间传递模型，该模型具体为：
；
其中，表示伪距的无电离层组合观测量，/>表示载波相位的无电离层组合观测量，i为卫星标识，r为接收机标识，c为光速，/>表示接收机钟差，/>表示卫星钟差，/>表示伪距噪声，/>表示载波相位噪声，/>为对流层延迟，/>为载波相位模糊度。
103.当未发生周跳时，根据上一个历元与当前历元的双频伪距观测值以及载波相位观测值建立历元间的差分解算模型；
建立历元间的差分解算模型为：
；
其中，为历元间一次差分运算符，k表示第k个历元，表示第k个历元的下一 个历元，i为卫星标识，r为接收机标识，表示卫星i与接收机r之间的距离，表示接收机 钟差，表示卫星钟差，表示伪距噪声，表示载波相位噪声，c为光速，表示伪 距的无电离层组合观测量，表示载波相位的无电离层组合观测量，表示电离层延迟。
在一个具体实施例中，若当前历元发生周跳，则在发生周跳的情况下通过标记周跳历元，并重新选取一个计算弧段，在该弧段的第二个历元处进行基于历元间差分的钟差解算模型的构建，从而进行历元间信息的提取和模型化。
104.基于载波相位时间传递模型以及历元间的差分解算模型进行钟差解算，获取当前历元两个时间传递站的GNSS接收机钟差；
基于载波相位时间传递模型以及历元间的差分解算模型进行钟差解算的方法为：
基于载波相位时间传递模型以及历元间的差分解算模型进行融合，构建历元间差分信息的接收机钟差解算模型，构建该接收机钟差解算模型的流程示意图如图2所示。
所述接收机钟差解算模型包括函数模型和随机模型，函数模型为：
；
其中，表示伪距的无电离层组合观测量，表示载波相位的无电离层组合观测 量，表示卫星i与接收机r之间的距离，表示接收机钟差，表示卫星钟差，表示伪 距噪声，表示载波相位噪声，c为光速，为对流层延迟，为载波相位模糊度，为历 元间一次差分运算符，为第k个历元与第个历元之间的钟差差分量。
随机模型为：
；
其中，为时间传递模型的权重分配，/>为时间传递模型中对应的中误差，/>为历元间的差分解算模型的权重分配，/>为历元间的差分解算模型中对应的中误差，在一个具体实施例中，通常按照2:1的经验值直接给出。
进行钟差解算时还包括对接收机钟差解算模型中的参数进行修正，包括：对流层延迟按照模型改正并结合随机游走的随机过程进行模拟和估计；载波相位模糊度在连续无周跳情况下作为常数估计；接收机钟差作为高斯白噪声逐历元估计。
105.根据当前历元两个时间传递站的GNSS接收机钟差获取时间传递量。
根据当前历元两个时间传递站的GNSS接收机钟差获取时间传递量的方法为：
获取两个时间传递站当前历元的接收机钟差量之间的差值，并根据该差值与时间传递链路中固有延迟之和，得到对应两个时间传递站之间的时间传递量，具体表达式为：
；
其中，为时间传递量，A和B为时间传递站的标识，/>为时间传递链路的固有延迟。
本发明建立基于历元间差分的钟差解算模型，能够有效获取历元间的差分信息，同时通过对基于历元间的差分信息模型化，并与载波相位时间传递模型进行融合，将GNSS观测量在时间域内的相关信息融入到GNSS载波相位时间传递模型中，能够考虑到观测量在时间域内的相关性，利用这些相关性提升了接收机钟速参数的估计精度，为进一步提升GNSS载波相位时间传递性能奠定了基础，同时为提升时间传递的可靠性提供了理论支撑。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种顾及历元间差分信息的GNSS载波相位时间传递方法，其特征在于，包括：
实时获取两个时间传递站每个历元的GNSS双频伪距观测值、载波相位观测值，以及对应的GNSS精密卫星轨道和钟差产品数据；
对所述双频伪距观测值以及载波相位观测值进行预处理，并根据预处理后的双频伪距观测值以及载波相位观测值构建载波相位时间传递模型；
当未发生周跳时，根据上一个历元与当前历元的双频伪距观测值以及载波相位观测值建立历元间的差分解算模型；
基于载波相位时间传递模型以及历元间的差分解算模型进行钟差解算，获取当前历元两个时间传递站的GNSS接收机钟差；
基于载波相位时间传递模型以及历元间的差分解算模型进行融合，构建历元间差分信息的接收机钟差解算模型，所述接收机钟差解算模型包括函数模型和随机模型，函数模型为：
；
其中，表示伪距的无电离层组合观测量，/>表示载波相位的无电离层组合观测量，表示卫星i与接收机r之间的距离，/>表示接收机钟差，/>表示卫星钟差，/>表示伪距噪声，/>表示载波相位噪声，/>为光速，/>为对流层延迟，/>为载波相位模糊度，/>为历元间一次差分运算符，/>为第k个历元与第/>个历元之间的钟差差分量；
随机模型为：
；
其中，为时间传递模型的权重分配，/>为时间传递模型中对应的中误差，/>为历元间的差分解算模型的权重分配，/>为历元间的差分解算模型中对应的中误差；
根据当前历元两个时间传递站的GNSS接收机钟差获取时间传递量。
2.根据权利要求1所述的一种顾及历元间差分信息的GNSS载波相位时间传递方法，其特征在于，对所述双频伪距观测值以及载波相位观测值进行预处理的方法为：
对获取到的GNSS双频伪距观测值以及载波相位观测值进行数据分类、野值剔除以及潮汐和天线相位中心误差的修正。
3.根据权利要求1所述的一种顾及历元间差分信息的GNSS载波相位时间传递方法，其特征在于，对所述双频伪距观测值以及载波相位观测值进行预处理，还包括：
对获取到的GNSS双频伪距观测值以及载波相位观测值进行周跳探测，并构建双频无电离层伪距和载波相位组合量，表示为：
；
其中，表示伪距的无电离层组合观测量，/>表示载波相位的无电离层组合观测量，和/>表示两个双频伪距观测量，/>和/>为两个双频载波相位观测量，/>和/>分别为两个频率的频点。
4.根据权利要求1所述的一种顾及历元间差分信息的GNSS载波相位时间传递方法，其特征在于，建立历元间的差分解算模型为：
；
其中，为历元间一次差分运算符，k表示第k个历元，/>表示第k个历元的下一个历元，i为卫星标识，r为接收机标识，/>表示卫星i与接收机r之间的距离，/>表示接收机钟差，/>表示卫星钟差，/>表示伪距噪声，/>表示载波相位噪声，/>为光速，/>表示伪距的无电离层组合观测量，/>表示载波相位的无电离层组合观测量，/>表示电离层延迟。
5.根据权利要求1所述的一种顾及历元间差分信息的GNSS载波相位时间传递方法，其特征在于，进行钟差解算时还包括：
对接收机钟差解算模型中的参数进行修正，包括：
对流层延迟按照模型改正并结合随机游走的随机过程进行模拟和估计；
载波相位模糊度在连续无周跳情况下作为常数估计；
接收机钟差作为高斯白噪声逐历元估计。
6.根据权利要求1所述的一种顾及历元间差分信息的GNSS载波相位时间传递方法，其特征在于：根据当前历元两个时间传递站的GNSS接收机钟差获取时间传递量的方法为：
获取两个时间传递站当前历元的接收机钟差量之间的差值，并根据该差值与时间传递链路中固有延迟之和，得到对应两个时间传递站之间的时间传递量。
7.根据权利要求1所述的一种顾及历元间差分信息的GNSS载波相位时间传递方法，其特征在于，根据预处理后的双频伪距观测值以及载波相位观测值构建载波相位时间传递模型，该模型具体为：
；
其中，表示伪距的无电离层组合观测量，/>表示载波相位的无电离层组合观测量，i为卫星标识，r为接收机标识，c为光速，/>表示接收机钟差，/>表示卫星钟差，/>表示伪距噪声，/>表示载波相位噪声，/>为对流层延迟，/>为载波相位模糊度。
CN202310542390.4A 2023-05-15 2023-05-15 一种顾及历元间差分信息的gnss载波相位时间传递方法 Active CN116299585B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310542390.4A CN116299585B (zh) 2023-05-15 2023-05-15 一种顾及历元间差分信息的gnss载波相位时间传递方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310542390.4A CN116299585B (zh) 2023-05-15 2023-05-15 一种顾及历元间差分信息的gnss载波相位时间传递方法

Publications (2)

Publication Number Publication Date
CN116299585A CN116299585A (zh) 2023-06-23
CN116299585B true CN116299585B (zh) 2023-09-08

Family

ID=86803443

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310542390.4A Active CN116299585B (zh) 2023-05-15 2023-05-15 一种顾及历元间差分信息的gnss载波相位时间传递方法

Country Status (1)

Country Link
CN (1) CN116299585B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117705099B (zh) * 2023-11-17 2024-05-28 武汉大学 基于变参考历元tdcp/mems imu的实时动态定位方法及系统

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5828336A (en) * 1996-03-29 1998-10-27 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Robust real-time wide-area differential GPS navigation
JPH11281728A (ja) * 1998-03-31 1999-10-15 Japan Aviation Electronics Ind Ltd キネマテック法ｇｐｓと慣性航法装置併用測位方法
CN106324629A (zh) * 2016-11-07 2017-01-11 中国电子科技集团公司第二十研究所 一种bds_gps_glonass融合精密单点定位方法
CN107607973A (zh) * 2017-09-12 2018-01-19 武汉大学 Gnss网络rtk参考站间模糊度快速固定方法及系统
CN107966722A (zh) * 2018-01-20 2018-04-27 中国人民解放军61540部队 一种gnss钟差解算方法
CN108196267A (zh) * 2017-12-20 2018-06-22 中国科学院国家授时中心 一种基于gnss cp技术的不间断时间传递方法
CN109683186A (zh) * 2018-12-20 2019-04-26 中国科学院国家授时中心 一种消除多卫星导航系统载波相位时间传递天跳变的方法
CN111323796A (zh) * 2020-03-18 2020-06-23 中国科学院国家空间科学中心 一种gnss接收机高采样钟差解算方法
WO2020133711A1 (zh) * 2018-12-28 2020-07-02 上海海积信息科技股份有限公司 一种卫星定轨方法、装置及电子设备
CN111522032A (zh) * 2020-05-25 2020-08-11 北京航空航天大学 一种北斗三号系统用户完好性处理的优化方法及优化装置
CN111751853A (zh) * 2020-06-20 2020-10-09 北京华龙通科技有限公司 一种gnss双频载波相位整周模糊度解算方法
CN111948682A (zh) * 2020-08-20 2020-11-17 山东科技大学 一种基于bds三频组合的伪距和载波相位随机模型构建方法
CN111965673A (zh) * 2020-06-24 2020-11-20 中山大学 基于多gnss的单频精密单点定位算法的时间频率传递方法
CN113253314A (zh) * 2021-05-11 2021-08-13 中国科学院国家授时中心 一种低轨卫星间时间同步方法及系统
CN113325446A (zh) * 2021-06-01 2021-08-31 中国科学院国家授时中心 一种多模共频gnss载波相位时间传递方法及系统
CN114002727A (zh) * 2021-10-14 2022-02-01 武汉理工大学 非组合载波相位历元间差分测速方法、系统、装置及介质
CN115201870A (zh) * 2022-07-04 2022-10-18 中国科学院精密测量科学与技术创新研究院 附先验约束的多频多模gnss非差非组合时间传递方法
CN115877415A (zh) * 2022-12-07 2023-03-31 中国科学院国家授时中心 一种超高精度空间航天器钟差估计方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050203702A1 (en) * 2004-03-12 2005-09-15 Sharpe Richard T. Method for backup dual-frequency navigation during brief periods when measurement data is unavailable on one of two frequencies
ES2426865T3 (es) * 2009-10-08 2013-10-25 Deutsches Zentrum für Luft- und Raumfahrt e.V. Resolución parcial de la ambigüedad para la estimación del retardo ionosférico multifrecuencia
US9405012B2 (en) * 2012-04-12 2016-08-02 Trimble Navigation Limited Advanced global navigation satellite systems (GNSS) positioning using precise satellite information
US10690777B2 (en) * 2014-07-31 2020-06-23 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. Multi-antenna-GNSS receiver-system to raise the probability of line of sight
EP3109673B1 (en) * 2015-06-24 2020-01-08 Centre National d'Etudes Spatiales Gnss receiver with an improved capability to resolve sub-carrier tracking ambiguities

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5828336A (en) * 1996-03-29 1998-10-27 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Robust real-time wide-area differential GPS navigation
JPH11281728A (ja) * 1998-03-31 1999-10-15 Japan Aviation Electronics Ind Ltd キネマテック法ｇｐｓと慣性航法装置併用測位方法
CN106324629A (zh) * 2016-11-07 2017-01-11 中国电子科技集团公司第二十研究所 一种bds_gps_glonass融合精密单点定位方法
CN107607973A (zh) * 2017-09-12 2018-01-19 武汉大学 Gnss网络rtk参考站间模糊度快速固定方法及系统
CN108196267A (zh) * 2017-12-20 2018-06-22 中国科学院国家授时中心 一种基于gnss cp技术的不间断时间传递方法
CN107966722A (zh) * 2018-01-20 2018-04-27 中国人民解放军61540部队 一种gnss钟差解算方法
CN109683186A (zh) * 2018-12-20 2019-04-26 中国科学院国家授时中心 一种消除多卫星导航系统载波相位时间传递天跳变的方法
WO2020133711A1 (zh) * 2018-12-28 2020-07-02 上海海积信息科技股份有限公司 一种卫星定轨方法、装置及电子设备
CN111323796A (zh) * 2020-03-18 2020-06-23 中国科学院国家空间科学中心 一种gnss接收机高采样钟差解算方法
CN111522032A (zh) * 2020-05-25 2020-08-11 北京航空航天大学 一种北斗三号系统用户完好性处理的优化方法及优化装置
CN111751853A (zh) * 2020-06-20 2020-10-09 北京华龙通科技有限公司 一种gnss双频载波相位整周模糊度解算方法
CN111965673A (zh) * 2020-06-24 2020-11-20 中山大学 基于多gnss的单频精密单点定位算法的时间频率传递方法
CN111948682A (zh) * 2020-08-20 2020-11-17 山东科技大学 一种基于bds三频组合的伪距和载波相位随机模型构建方法
CN113253314A (zh) * 2021-05-11 2021-08-13 中国科学院国家授时中心 一种低轨卫星间时间同步方法及系统
CN113325446A (zh) * 2021-06-01 2021-08-31 中国科学院国家授时中心 一种多模共频gnss载波相位时间传递方法及系统
CN114002727A (zh) * 2021-10-14 2022-02-01 武汉理工大学 非组合载波相位历元间差分测速方法、系统、装置及介质
CN115201870A (zh) * 2022-07-04 2022-10-18 中国科学院精密测量科学与技术创新研究院 附先验约束的多频多模gnss非差非组合时间传递方法
CN115877415A (zh) * 2022-12-07 2023-03-31 中国科学院国家授时中心 一种超高精度空间航天器钟差估计方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵坤娟等.基于北斗GEO 卫星双频观测值连续监测电离层延迟.宇航计测技术.2018,第38卷(第4期),67-72. *

Also Published As

Publication number Publication date
CN116299585A (zh) 2023-06-23

Similar Documents

Publication Publication Date Title
CN108919634B (zh) 一种北斗三频非差非组合观测值时间传递系统及方法
CN111045034B (zh) 基于广播星历的gnss多系统实时精密时间传递方法及系统
US20210223405A1 (en) Fast and Precise Positioning Method and System
CN105607089B (zh) 一种基于伪距辅助载波相位的北斗共视时间传递的方法
CN104483690B (zh) 一种gnss三频精密单点定位模糊度固定方法
CN103344978A (zh) 一种适用于大规模用户的区域增强精密定位服务方法
CN105301617B (zh) 一种卫星导航系统中的整周模糊度有效性检验方法
CN104483691B (zh) 一种gnss组合精密单点定位方法
CN111580132B (zh) 基于时间实验室增强信息的北斗局域精密时间传递方法
CN105842720B (zh) 一种基于载波相位的广域精密实时定位方法
CN116299585B (zh) 一种顾及历元间差分信息的gnss载波相位时间传递方法
CN104656108A (zh) 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法
CN110208836B (zh) 基于卡尔曼滤波的gnss高适应性周跳探测与修复方法
CN104483689B (zh) 一种bds参考站三频载波相位整周模糊度确定方法
CN113253314B (zh) 一种低轨卫星间时间同步方法及系统
CN113325446B (zh) 一种多模共频gnss载波相位时间传递方法及系统
CN103197335A (zh) 采用改进正则化方法抑制dgps整周模糊度病态性的方法
CN115933356B (zh) 一种虚拟原子钟的高精度时间同步系统和方法
CN109212562A (zh) 一种生成载波伪距观测量的方法
EP3901667A1 (en) Gnss receiver adapted to produce, use, and communicate software-generated satellite signal data
CN114879239B (zh) 一种增强瞬时ppp固定解的区域三频整数钟差估计方法
CN112415542A (zh) 一种基于北斗和ins结合的变形监测解算方法
CN111123295A (zh) 基于ssr的定位方法及装置、定位系统
CN116540280B (zh) 多频卫星导航数据状态域改正信息综合处理方法和系统
Yin et al. A novel cycle slips detection model for the high precision positioning

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