CN115267843B - 多频多模gnss高频精密卫星钟差的实时非差估计方法 - Google Patents

多频多模gnss高频精密卫星钟差的实时非差估计方法 Download PDF

Info

Publication number
CN115267843B
CN115267843B CN202210718533.8A CN202210718533A CN115267843B CN 115267843 B CN115267843 B CN 115267843B CN 202210718533 A CN202210718533 A CN 202210718533A CN 115267843 B CN115267843 B CN 115267843B
Authority
CN
China
Prior art keywords
frequency
estimation
satellite
difference
real
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
CN202210718533.8A
Other languages
English (en)
Other versions
CN115267843A (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.)
Institute of Precision Measurement Science and Technology Innovation of CAS
Original Assignee
Institute of Precision Measurement Science and Technology Innovation 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 Institute of Precision Measurement Science and Technology Innovation of CAS filed Critical Institute of Precision Measurement Science and Technology Innovation of CAS
Priority to CN202210718533.8A priority Critical patent/CN115267843B/zh
Publication of CN115267843A publication Critical patent/CN115267843A/zh
Application granted granted Critical
Publication of CN115267843B publication Critical patent/CN115267843B/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/27—Acquisition or tracking or demodulation of signals transmitted by the system creating, predicting or correcting ephemeris or almanac data within the receiver
    • 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/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • G01S19/072—Ionosphere corrections
    • 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/32—Multimode operation in a single same satellite system, e.g. GPS L1/L2
    • 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • G01S19/426—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Power Engineering (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种多频多模GNSS高频精密卫星钟差的实时非差估计方法。基于低频线程和高频线程两组非差并行处理的参数估计模式：低频线程利用伪距和相位观测值信息，构建低频满秩模型进行实时钟差估计，实现相位模糊度、天顶对流层延迟、斜电离层延迟和伪距偏差联合估计；高频线程建立高频估计模型，在实时钟差估计中，对高频相位观测值添加相位模糊度改正，同时将天顶对流层延迟作为滤波的先验信息实现实时高频卫星钟差估计。在低频线程和高频线程实时处理模式中，各线程独立并行进行参数估计，优化参数解算效率，以满足实时钟差估计的时效性需求。

Description

多频多模GNSS高频精密卫星钟差的实时非差估计方法
技术领域
本发明涉及高精度定位领域，尤其涉及一种多频多模GNSS高频精密卫星钟差的实时非差估计方法。
背景技术
高精度GNSS定位服务依赖于增强信息，精密卫星钟差作为增强信息中的重要内容，在精密定位应用中起着关键作用。随着GNSS实时高精度应用服务范围的不断拓展，其对GNSS实时卫星钟差估计提出了更迫切的要求。实时卫星钟差估计通常采用历元间单差或混合差分模式进行，此方式虽然可以消除相位模糊度带来的大量计算负荷影响，保证时效性需求，却牺牲了观测值中大量的有用信息。此外，由于当前的实时钟差更新频率基本在5s或者10s，无法满足部分高频实时PPP的应用。
发明内容
本发明主要目的在于：提供一种多频多模GNSS高频精密卫星钟差的实时非差估计方法，有助于实现实时高频数据处理。
本发明所采用的技术方案是：一种多频多模GNSS高频精密卫星钟差的实时非差估计方法，本方法包括：
基于低频线程和高频线程并行的参数估计模式；
所述的低频线程利用伪距和相位观测值信息，构建低频满秩模型，基于非差模式进行实时钟差估计，实现相位模糊度、天顶对流层延迟、斜电离层延迟和伪距偏差联合估计；
所述的高频线程建立高频估计模型，基于非差模式进行实时钟差估计，对高频相位观测值添加相位模糊度改正，同时将天顶对流层延迟作为滤波的先验信息实现实时高频卫星钟差估计；
在低频线程和高频线程实时处理模式中，GPS、GLOANSS、BDS、Galileo四系统独立并行进行参数估计，优化参数解算效率，以满足实时钟差时效性需求。
按上述方法，所述的低频满秩模型的构建过程如下：
GNSS非组合观测值原始伪距和载波相位方程表示为：
Figure BDA0003694111860000021
Figure BDA0003694111860000022
式中，
Figure BDA0003694111860000023
和/>
Figure BDA0003694111860000024
分别为伪距和相位的观测值与计算值的差值；T表示的是系统类型；S表示的是卫星号；r为对应的接收机编号；f为频率号，f＝1,2,3...；c为真空光速；tr和tS分别为接收机钟差和卫星钟差；br,f和/>
Figure BDA0003694111860000025
为频率相关的接收机伪距硬件延迟与卫星伪距硬件延迟；/>
Figure BDA0003694111860000026
为天顶对流层延迟；/>
Figure BDA0003694111860000027
为电离层斜延迟的系数；/>
Figure BDA0003694111860000028
为第一频点的斜电离层延迟；λf为频率f对应的载波相位波长；/>
Figure BDA0003694111860000029
为载波相位模糊度，其吸收接收机端和卫星端相位硬件延迟；/>
Figure BDA00036941118600000210
为伪距观测噪声及其他未模型化误差之和；/>
Figure BDA00036941118600000211
为载波相位观测噪声及其他未模型化误差之和；
构建多个频点统一的接收机钟差参数，将接收机钟差吸收接收机端伪距硬件延迟的消电离层组合，则组合后的接收机钟差参数表示为：
Figure BDA00036941118600000212
其中
Figure BDA00036941118600000213
构建多个频点统一的卫星钟差参数，将卫星钟差吸收接卫星端伪距硬件延迟的消电离层组合，则组合后的卫星钟差参数表示为：
Figure BDA00036941118600000214
其中
Figure BDA00036941118600000215
统一伪距和载波相位观测方程中对应的电离层信息，则参数重组的电离层吸收接收机端和卫星端之间的伪距硬件延迟线性组合，第1频点的电离层参数的组合形式表示为：
Figure BDA00036941118600000216
其他频点的电离层延迟基于μf和/>
Figure BDA00036941118600000217
线性表示；其中
Figure BDA00036941118600000218
Figure BDA00036941118600000219
基于以上的参数重组统一卫星钟差、接收机钟差、电离层参数后，相位模糊度参数因此也吸收接收机端和卫星端之间的伪距偏差的线性组合：
Figure BDA00036941118600000220
将参数重组后的参数代入原始的GNSS非组合模型中，则双频非组合卫星钟差估计模型表示为：
Figure BDA00036941118600000221
考虑到三频信号和多频信号与L1、L2对应的伪距偏差差异，因此在其对应的伪距观测方程中添加对应的接收机伪距偏差
Figure BDA00036941118600000222
和卫星伪距偏差/>
Figure BDA00036941118600000223
参数，改正后的多频多模GNSS卫星钟差估计的低频满秩模型表示为：
Figure BDA0003694111860000031
对应的伪距偏差的形式表示为：
Figure BDA0003694111860000032
按上述方法，在低频满秩模型构建过程中，接收机钟差与卫星钟差之间相关性引起的秩亏，各系统通过所有卫星钟差构建零均值基准
Figure BDA0003694111860000033
进行约束，m为各系统对应的全星座信息中最大卫星编号；
频率编号大于2时，接收机伪距偏差和卫星伪距偏差之间的相关性引起的秩亏问题，基于零均值基准
Figure BDA0003694111860000034
进行约束。
按上述方法，所述的改正后的多频多模GNSS卫星钟差估计的低频满秩模型，采用扩展的卡尔曼滤波实现参数估计；
扩展的卡尔曼滤波参数估计对应的观测方程和状态方程形式表示为
Figure BDA0003694111860000035
其中Xj为j时刻的状态向量，Aj和/>
Figure BDA0003694111860000036
为实际观测值Lj和预报参数Xp,j的设计矩阵，其可经过观测方程线性化得到；Φj为状态转移矩阵，此处为单位阵，Wj为系统噪声；预报参数Xp,j的协方差表示为/>
Figure BDA0003694111860000037
此处为单位阵。经量测更新对应的状态向量Xj通过/>
Figure BDA0003694111860000038
对应的状态向量的方差表示为
Figure BDA0003694111860000039
按上述方法，高频估计模型引入对应的相位模糊度改正，基于L1和L2双频消电离层组合构建如下模型：
Figure BDA00036941118600000310
其中对应的消电离层组合相位模糊度改正信息基于低频模型顾及的L1和L2模糊度线性组合表示为
Figure BDA00036941118600000311
且对应的相位观测值方差信息为/>
Figure BDA00036941118600000312
在此过程中，高频观测值/>
Figure BDA00036941118600000313
在高频第j历元，引入低频第i历元的相位模糊度信息/>
Figure BDA00036941118600000314
进行改正，并且在下一个低频历元i+1完成解算之前，高频历元仍然沿用相同的相位模糊度信息。
按上述方法，所述高频估计模型中，接收机钟差和卫星钟差参数之间相关性导致的秩亏问题，采用各系统内所有卫星钟差构建零均值基准
Figure BDA0003694111860000041
进行约束。
一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述多频多模GNSS高频精密卫星钟差的实时非差估计方法的步骤。
本发明产生的有益效果是：
1、本发明采用高、低频双线程并行处理，高频模式就近采用低频相位模糊度改正可减小参数解算负荷，提升实时卫星钟差估计效率，同时分系统并行处理的方式，能够降低法方程维度，便于参数高效求解，有助于实现实时高频数据处理；低频线程和高频线程在进行实时钟差估计时均基于非差模式进行，能够保留观测信息，具备较高的钟差估计精度和可靠性，同时可以实现相位模糊度固定，促进估计的卫星钟差产品的高精度定位应用。
2、本发明使用的函数模型能兼容双频、三频及以上多频GNSS卫星钟差估计。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1为完整的GNSS实时卫星钟差估计算法及其流程图。
图2为基于高低频并行的非差非组合多频多模实时卫星钟差估计模块流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
本发明所涉及的数据处理，是基于实时数据流和SSR改正信息进行，结合图1，实际处理过程如下：
Ntrip用户向Ntrip Caster发送请求，接收其提供的实时数据流和SSR信息，根据具体需求将其存入计算机共享内存，然后同步缓存观测数据和SSR信息，避免由于时延造成的不稳定性，然后对实时数据质量控制和实时精密卫星轨道产品恢复等数据预处理工作。在用户解码SSR信息恢复轨道信息中，获得参考时刻tr的轨道(δOr,δOa,δOc)、
Figure BDA0003694111860000042
信息后，基于三维方向向量可获得对应的ECEF坐标系下的轨道改正数信息
Figure BDA0003694111860000051
然后基于广播星历还原得到卫星位置/>
Figure BDA0003694111860000052
本发明提供一种多频多模GNSS高频精密卫星钟差的实时非差估计方法，基于预处理后的实时数据和SSR信息进行GNSS实时卫星钟差解算，实时钟差解算通过高低频并行的参数估计模式进行。所述的低频线程利用伪距和相位观测值信息，构建低频满秩模型，基于非差模式进行实时钟差估计，实现相位模糊度、天顶对流层延迟、斜电离层延迟和伪距偏差联合估计；所述的高频线程建立高频估计模型，基于非差模式进行实时钟差估计，对高频相位观测值添加相位模糊度改正，同时将天顶对流层延迟作为滤波的先验信息实现实时高频卫星钟差估计；在低频线程和高频线程实时处理模式中，各系统独立并行进行参数估计，优化参数解算效率，以满足实时钟差时效性需求。所述的系统为GPS、GLONASS、BDS和Galileo4个子模块。
基于RTCM协议将对应的实时高频卫星钟差编码成SSR信息，经过Ntrip Caster对外播发。编码过程中，广播星历中对应参考时刻tr的卫星钟差参考信息为(a0,a1,a2)，对应的广播钟差可表示为
Figure BDA0003694111860000053
c为真空光速。实时精密钟差表示为多项式的形式dts＝[c0+c1·(t-tr)+c2·(t-tr)2]/c。则钟差改正数(Δc0,Δc1,Δc2)可表示为/>
Figure BDA0003694111860000054
发送至SSR的钟差改正信息。实时高精度用户基于接收的SSR改正信息进行解码，开展高精度定位应用。
如图2所示，具体的高低频处理过程如下：
GNSS非组合观测值原始伪距和载波相位方程表示为：
Figure BDA0003694111860000055
Figure BDA0003694111860000056
式中，
Figure BDA0003694111860000057
和/>
Figure BDA0003694111860000058
分别为伪距和相位的观测值与计算值的差值；T表示的是系统类型；S表示的是卫星号；r为对应的接收机编号；f为频率号，f＝1,2,3...；c为真空光速；tr和tS分别为接收机钟差和卫星钟差；br,f和/>
Figure BDA0003694111860000059
为频率相关的接收机伪距硬件延迟与卫星伪距硬件延迟；/>
Figure BDA00036941118600000510
为天顶对流层延迟；/>
Figure BDA00036941118600000511
为电离层斜延迟系数；/>
Figure BDA00036941118600000512
为第1频点的斜电离层延迟；λf为频率f对应的载波相位波长；/>
Figure BDA00036941118600000513
为载波相位模糊度，其吸收接收机端和卫星端相位硬件延迟；/>
Figure BDA0003694111860000061
为伪距观测噪声及其他未模型化误差之和；/>
Figure BDA0003694111860000062
为载波相位观测噪声及其他未模型化误差之和；/>
Figure BDA0003694111860000063
为频率f对应的载波相位波长；
构建多个频点统一的接收机钟差参数，将接收机钟差吸收接收机端伪距硬件延迟的消电离层组合，则组合后的接收机钟差参数表示为：
Figure BDA0003694111860000064
其中
Figure BDA0003694111860000065
构建多个频点统一的卫星钟差参数，将卫星钟差吸收接卫星端伪距硬件延迟的消电离层组合，则组合后的卫星钟差参数表示为：
Figure BDA0003694111860000066
其中
Figure BDA0003694111860000067
统一伪距和载波相位观测方程中对应的电离层信息，则参数重组的电离层吸收接收机端和卫星端之间的伪距硬件延迟线性组合，第1频点的电离层参数的组合形式表示为：
Figure BDA0003694111860000068
其他频点的电离层延迟基于μf和/>
Figure BDA0003694111860000069
线性表示；其中
Figure BDA00036941118600000610
Figure BDA00036941118600000611
基于以上的参数重组统一卫星钟差、接收机钟差、电离层参数后，相位模糊度参数因此也吸收接收机端和卫星端之间的伪距偏差的线性组合：
Figure BDA00036941118600000612
将参数重组后的参数代入原始的GNSS非组合模型中，则双频非组合卫星钟差估计模型表示为：
Figure BDA00036941118600000613
考虑到三频信号和多频信号与L1、L2对应的伪距偏差差异，因此在其对应的伪距观测方程中添加对应的接收机伪距偏差
Figure BDA00036941118600000614
和卫星伪距偏差/>
Figure BDA00036941118600000615
参数，改正后的多频多模GNSS卫星钟差估计的低频满秩模型表示为：
Figure BDA00036941118600000616
对应的伪距偏差的形式表示为：
Figure BDA00036941118600000617
/>
Figure BDA00036941118600000618
和/>
Figure BDA00036941118600000619
分别为当频率编号大于2时对应的接收机伪距偏差和卫星伪距偏差信息。
所述的改正后的多频多模GNSS卫星钟差估计的低频满秩模型，采用扩展的卡尔曼滤波实现参数估计。扩展的卡尔曼滤波参数估计对应的观测方程和状态方程形式表示为
Figure BDA0003694111860000071
其中Xj为j时刻的状态向量，Aj和/>
Figure BDA0003694111860000072
为实际观测值Lj和预报参数Xp,j的设计矩阵，其可经过观测方程线性化得到；Φj为状态转移矩阵，此处为单位阵，Wj为系统噪声；预报参数Xp,j的协方差表示为/>
Figure BDA0003694111860000073
经量测更新对应的状态向量Xj通过
Figure BDA0003694111860000074
对应的状态向量的方差表示为
Figure BDA0003694111860000075
基于扩展的卡尔曼滤波，针对不同的待估参数添加了虚拟观测值，使得参数处理过程更加灵活高效。为高低频处理中不同待估参数的灵活改正和估计提供了合理的解决方案。
在低频满秩模型构建过程中，接收机钟差与卫星钟差之间相关性引起的秩亏，各系统通过所有卫星钟差构建零均值基准
Figure BDA0003694111860000076
进行约束，m为各系统对应的全星座信息中最大卫星编号；频率编号大于2时，接收机伪距偏差和卫星伪距偏差之间的相关性引起的秩亏问题，基于零均值基准/>
Figure BDA0003694111860000077
进行约束。
在此参数估计的过程中，形成GPS、GLONASS、BDS和Galileo 4个子模块独立运行，对应求解的对流层参数可作为高频处理中的先验信息，相位模糊度参数实现对高频处理中的参数改正，低频和高频两者之间通过共享内存建立联系。
高频估计模型引入对应的相位模糊度改正，基于L1和L2双频消电离层组合构建如下模型：
Figure BDA0003694111860000078
此过程依赖于双频相位观测值经过消电离层组合实现高频钟差估计。而在此过程中，对应的消电离层组合相位模糊度改正信息基于低频模型顾及的L1和L2相位模糊度线性组合表示为
Figure BDA0003694111860000079
重组后对高频相位观测值进行改正，且对应的相位观测值方差信息为/>
Figure BDA00036941118600000710
所述高频估计模型中，接收机钟差和卫星钟差参数之间相关性导致的秩亏问题，采用各系统内所有卫星钟差构建零均值基准
Figure BDA00036941118600000711
进行约束。
扩展卡尔曼滤波参数估计中，针对高低频处理过程中对应的对流层参数和低频中的相位模糊度参数构建虚拟观测值信息，有助于处理实时参数估计中卫星升降和周跳处理等问题；同时对于高低频处理时效性的差异，在高低频并行处理中，低频解算完成后及时对高频中的相位模糊度和对流层信息进行更新，在低频运行时间段内继续保持原有的相位模糊度信息进行改正，高效的保证高频估计参数的高精度。
本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述多频多模GNSS高频精密卫星钟差实时非差估计方法的步骤。
应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (6)

1.一种多频多模GNSS高频精密卫星钟差的实时非差估计方法，其特征在于，本方法包括：
基于低频线程和高频线程并行的参数估计模式；
所述的低频线程利用伪距和相位观测值信息构建低频满秩模型，基于非差模式进行实时钟差估计，实现相位模糊度、天顶对流层延迟、斜电离层延迟和伪距偏差联合估计；
所述的高频线程建立高频估计模型，基于非差模式进行实时钟差估计，对高频相位观测值添加相位模糊度改正，同时将天顶对流层延迟作为滤波的先验信息实现实时高频卫星钟差估计；
在低频线程和高频线程实时处理模式中，GPS、GLOANSS、BDS、Galileo四系统独立并行进行参数估计，优化参数解算效率，以满足实时钟差时效性需求；
所述的低频满秩模型的构建过程如下：
GNSS非组合观测值原始伪距和载波相位方程表示为：
Figure FDA0004189868440000011
Figure FDA0004189868440000012
式中，
Figure FDA0004189868440000013
和/>
Figure FDA0004189868440000014
分别为伪距和相位的观测值与计算值的差值；T表示的是系统类型；S表示的是卫星号；r为对应的接收机编号；f为频率号，f＝1,2,3...；c为真空光速；tr和tS分别为接收机钟差和卫星钟差；br,f和/>
Figure FDA0004189868440000015
为频率相关的接收机伪距硬件延迟与卫星伪距硬件延迟；Zr为天顶对流层延迟；/>
Figure FDA0004189868440000016
为电离层斜延迟的系数；/>
Figure FDA0004189868440000017
为第1频点的斜电离层延迟；λf为频率f对应的载波相位波长；/>
Figure FDA0004189868440000018
为相位模糊度，其吸收接收机端和卫星端相位硬件延迟，因此具有浮点解特性；/>
Figure FDA0004189868440000019
为伪距观测噪声及其他未模型化误差之和；/>
Figure FDA00041898684400000110
为载波相位观测噪声及其他未模型化误差之和；
构建多个频点统一的接收机钟差参数，将接收机钟差吸收第1、2两个频点接收机端伪距硬件延迟的消电离层组合，则组合后的接收机钟差参数表示为：
Figure FDA0004189868440000021
其中
Figure FDA0004189868440000022
构建多个频点统一的卫星钟差参数，将卫星钟差吸收第1、2两个频点卫星端伪距硬件延迟的消电离层组合，则组合后的卫星钟差参数表示为：
Figure FDA0004189868440000023
其中
Figure FDA0004189868440000024
统一伪距和载波相位观测方程中对应的电离层信息，则经过参数重组后的电离层吸收接收机端和卫星端之间的伪距硬件延迟线性组合，第1个频点的电离层参数的组合形式表示为：
Figure FDA0004189868440000025
其他频点的电离层延迟基于μf和/>
Figure FDA0004189868440000026
线性表示；其中
Figure FDA0004189868440000027
基于以上的参数重组统一卫星钟差、接收机钟差、电离层参数后，相位模糊度参数因此也吸收接收机端和卫星端之间的伪距偏差的线性组合：
Figure FDA0004189868440000028
将参数重组后的参数代入原始的GNSS非组合模型中，则双频非组合卫星钟差估计模型表示为：
Figure FDA0004189868440000029
考虑到三频信号和多频信号与L1、L2对应的伪距偏差差异，因此在其对应的伪距观测方程中添加对应的接收机伪距偏差
Figure FDA00041898684400000210
和卫星伪距偏差/>
Figure FDA00041898684400000211
参数，改正后的多频多模GNSS卫星钟差估计的低频满秩模型表示为：
Figure FDA00041898684400000212
对应的伪距偏差的形式表示为：
Figure FDA00041898684400000213
2.根据权利要求1所述的多频多模GNSS高频精密卫星钟差的实时非差估计方法，其特征在于，在低频满秩模型构建过程中，接收机钟差与卫星钟差之间相关性引起的秩亏，各系统通过所有卫星钟差构建零均值基准
Figure FDA0004189868440000031
进行约束，m为各系统对应的全星座中最大卫星编号；
频率编号大于2时，接收机伪距偏差和卫星伪距偏差之间的相关性引起的秩亏问题，基于零均值基准
Figure FDA0004189868440000032
进行约束。
3.根据权利要求1所述的多频多模GNSS高频精密卫星钟差的实时非差估计方法，其特征在于，所述的改正后的多频多模GNSS卫星钟差估计的低频满秩模型，采用扩展的卡尔曼滤波实现参数估计；
扩展的卡尔曼滤波参数估计对应的观测方程和状态方程形式表示为
Figure FDA0004189868440000033
其中Xj为j时刻的状态向量，Aj和/>
Figure FDA0004189868440000034
为实际观测值Lj和预报参数Xp,j的设计矩阵，经过观测方程线性化得到；Φj为状态转移矩阵，此处为单位阵，Wj为系统噪声；预报参数Xp,j的协方差表示为/>
Figure FDA0004189868440000035
经量测更新对应的状态向量Xj通过
Figure FDA0004189868440000036
对应的状态向量的方差表示为
Figure FDA0004189868440000037
4.根据权利要求1至3中任意一项所述的多频多模GNSS高频精密卫星钟差的实时非差估计方法，其特征在于，高频估计模型引入对应的相位模糊度改正，基于第1和2频点的双频消电离层组合构建如下模型：
Figure FDA0004189868440000038
其中对应的消电离层组合相位模糊度改正信息基于低频模型估计的第1和第2频点相位模糊度线性组合表示
Figure FDA0004189868440000039
且对应的相位观测值方差信息为：/>
Figure FDA00041898684400000310
在此过程中，高频观测值/>
Figure FDA00041898684400000311
在高频第j历元，引入低频第i历元的相位模糊度信息/>
Figure FDA00041898684400000312
进行改正，并且在下一个低频历元i+1完成解算之前，高频历元仍然沿用相同的相位模糊度信息。
5.根据权利要求4所述的多频多模GNSS高频精密卫星钟差的实时非差估计方法，其特征在于，所述高频估计模型中，接收机钟差和卫星钟差参数之间相关性导致的秩亏问题，采用各系统内所有卫星钟差构建零均值基准
Figure FDA0004189868440000041
进行约束。
6.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至5中任意一项所述多频多模GNSS高频精密卫星钟差的实时非差估计方法的步骤。
CN202210718533.8A 2022-06-14 2022-06-14 多频多模gnss高频精密卫星钟差的实时非差估计方法 Active CN115267843B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210718533.8A CN115267843B (zh) 2022-06-14 2022-06-14 多频多模gnss高频精密卫星钟差的实时非差估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210718533.8A CN115267843B (zh) 2022-06-14 2022-06-14 多频多模gnss高频精密卫星钟差的实时非差估计方法

Publications (2)

Publication Number Publication Date
CN115267843A CN115267843A (zh) 2022-11-01
CN115267843B true CN115267843B (zh) 2023-05-30

Family

ID=83762118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210718533.8A Active CN115267843B (zh) 2022-06-14 2022-06-14 多频多模gnss高频精密卫星钟差的实时非差估计方法

Country Status (1)

Country Link
CN (1) CN115267843B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116338754B (zh) * 2023-05-25 2023-09-12 中国科学院精密测量科学与技术创新研究院 一种固定模糊度的非差非组合gnss时间频率传递方法
CN116540280B (zh) * 2023-07-07 2023-09-15 中国科学院空天信息创新研究院 多频卫星导航数据状态域改正信息综合处理方法和系统
CN116626723B (zh) * 2023-07-21 2023-10-17 中国科学院空天信息创新研究院 GNSS的LcUPD和TcUPD的估计方法及装置
CN117031507A (zh) * 2023-10-08 2023-11-10 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) 适用于BDS-3 B1CB2a双频信号的精密单点定位方法
CN117169941B (zh) * 2023-11-01 2024-01-16 中国船舶集团有限公司第七〇七研究所 一种基于非差osr改正数增强的rappp定位方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109358487A (zh) * 2018-10-10 2019-02-19 武汉大学 一种基于gnss精密授时的伪卫星系统及方法
CN109597105A (zh) * 2018-12-13 2019-04-09 东南大学 一种顾及载波系统间偏差的gps/glonass紧组合定位方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2985631B1 (en) * 2014-08-14 2019-08-07 Trimble Inc. Navigation satellite system based positioning involving the generation of receiver-specific or receiver-type-specific correction information
WO2017181221A1 (en) * 2016-04-18 2017-10-26 Spatial Information Systems Research Ltd A method of analysing multiple signals transmitted by discrete global navigation satellite systems
CN107356947B (zh) * 2017-05-31 2019-06-18 中国科学院测量与地球物理研究所 基于单频导航卫星数据确定卫星差分伪距偏差的方法
CN108387912B (zh) * 2018-01-31 2020-05-08 中国科学院测量与地球物理研究所 一种Multi-GNSS精密单点定位的解算方法
CN111025346A (zh) * 2019-11-18 2020-04-17 广州南方卫星导航仪器有限公司 一种快速估计gnss精密卫星钟差的方法及存储介质
CN111190203B (zh) * 2020-01-14 2021-10-08 中国科学院测量与地球物理研究所 基于非差观测模型和消秩亏的ppp-rtk定位方法
CN112014860B (zh) * 2020-07-20 2023-07-14 中国科学院空天信息创新研究院 一种基于北斗ppp-rtk的低轨卫星时空基准建立方法
CN112505736B (zh) * 2020-11-06 2022-08-30 中国科学院精密测量科学与技术创新研究院 基于多频多系统的系统内和系统间偏差建模的rtk定位方法
CN113325453B (zh) * 2021-06-22 2022-11-15 中国科学院精密测量科学与技术创新研究院 基于参数约束的gnss非差模糊度确定方法及快速定位方法
CN113885055A (zh) * 2021-09-07 2022-01-04 中国科学院上海天文台 星地双向钟差约束的导航卫星多星定轨系统及方法
CN113848572A (zh) * 2021-09-16 2021-12-28 东南大学 一种基于大气误差增强的多频ppp序贯单历元定位方法
CN114019550A (zh) * 2021-10-25 2022-02-08 中国电子科技集团公司第五十四研究所 模糊度加权约束的gnss非差非组合观测方程解算方法
CN114019551A (zh) * 2021-10-26 2022-02-08 中国电子科技集团公司第五十四研究所 一种gnss观测站网原始观测方程解算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109358487A (zh) * 2018-10-10 2019-02-19 武汉大学 一种基于gnss精密授时的伪卫星系统及方法
CN109597105A (zh) * 2018-12-13 2019-04-09 东南大学 一种顾及载波系统间偏差的gps/glonass紧组合定位方法

Also Published As

Publication number Publication date
CN115267843A (zh) 2022-11-01

Similar Documents

Publication Publication Date Title
CN115267843B (zh) 多频多模gnss高频精密卫星钟差的实时非差估计方法
JP7122023B2 (ja) 高速精密測位方法及びシステム
Xu et al. Performance evaluation of BDS-3 PPP-B2b precise point positioning service
CN112014860B (zh) 一种基于北斗ppp-rtk的低轨卫星时空基准建立方法
CN111190203B (zh) 基于非差观测模型和消秩亏的ppp-rtk定位方法
Leandro et al. RTX positioning: The next generation of cm-accurate real-time GNSS positioning
Elsobeiey Precise point positioning using triple-frequency GPS measurements
WO2011034615A2 (en) Gnss signal processing with rover ambiguity fixing
Li et al. Review of PPP–RTK: Achievements, challenges, and opportunities
CN111638535B (zh) 一种用于gnss实时精密单点定位的混合模糊度固定方法
CN111273327A (zh) 一种基于组合和非组合混合观测模型的精密单点定位方法
CN115933356B (zh) 一种虚拟原子钟的高精度时间同步系统和方法
CN112526564A (zh) 一种精密单点定位重新收敛方法
CN104316943A (zh) 一种伪距离和多普勒组合差分定位系统及方法
Liu et al. An efficient undifferenced method for estimating multi-GNSS high-rate clock corrections with data streams in real time
CN114879239B (zh) 一种增强瞬时ppp固定解的区域三频整数钟差估计方法
CN114966760A (zh) 一种电离层加权的非差非组合ppp-rtk技术实现方法
WO2023197487A1 (zh) 一种应用于全球导航卫星系统gnss的周跳探测与修复方法
Cao et al. BDS-3/GNSS multi-frequency precise point positioning ambiguity resolution using observable-specific signal bias
CN114935770A (zh) 一种多历元加快精密单点定位收敛速度的方法及装置
Gu et al. Real-time precise point positioning based on BDS-3 global short message communication
Choy et al. An evaluation of various ionospheric error mitigation methods used in single frequency PPP
CN116679329A (zh) 基于状态域空间的北斗时空基准精确表征方法和系统
Bisnath et al. Innovation: Examining precise point positioning now and in the future
CN110095796B (zh) 一种跨系统mw紧组合实时动态精密导航定位方法

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