多模GNSS融合PPP系统间偏差特性分析

多模全球导航卫星系统（Global Navigation Satellite System,GNSS）精密单点定位（precise point positioning,PPP）存在系统间偏差（inter-system bias,ISB）,构建了顾及系统间偏差的多模GNSS融合PPP算法,对多星座实验（the multi-GNSS experiment,MGEX）监测网中的7个测站观测数据进行静态解算,获得Galileo、GLONASS、北斗与全球定位系统之间的ISB值。分析结果表明,四系统PPP融合定位在水平分量和高程分量的精度分别为8.9 mm、5.3 mm和10.9 mm,体现出较高的融合定位精度。不同系统ISB值在单天内的稳定性较好,均优于0.12 ns。从多天ISB序列看,ISB存在不规律跳变,变化幅度可达近20 ns。不同类型接收机ISB存在一定差异,同一类型接收机结果相近。综合来看,Galileo ISB值最稳定且结果最优,北斗与GLONASS结果相当。     

Characterization of between-receiver GPS-Galileo inter-system biases and their effect on mixed ambiguity resolution

The Global Positioning System (GPS) and Galileo will transmit signals on similar frequencies, that is, the L1–E1 and L5–E5a frequencies. This will be beneficial for mixed GPS and Galileo applications in which the integer carrier phase ambiguities need to be resolved, in order to estimate the positioning unknowns with centimeter accuracy or better. In this contribution, we derive the mixed GPS + Galileo model that is based on “inter-system” double differencing, that is, differencing the Galileo phase and code observations relative to those corresponding to the reference or pivot satellite of GPS. As a consequence of this, additional between-receiver inter-system bias (ISB) parameters need to be solved as well for both phase and code data. We investigate the size and variability of these between-receiver ISBs, estimated from L1 and L5 observations of GPS, as well as E1 and E5a observations of the two experimental Galileo In-Orbit Validation Element (GIOVE) satellites. The data were collected using high-grade multi-GNSS receivers of different manufacturers for several zero- and short-baseline setups in Australia and the USA. From this analysis, it follows that differential ISBs are only significant for receivers of different types and manufacturers; for baselines formed by identical receiver types, no differential ISBs have shown up; thus, implying that the GPS and GIOVE data are then fully interoperable. Fortunately, in case of different receiver types, our analysis also indicates that the phase and code ISBs may be calibrated, since their estimates, based on several datasets separated in time, are shown to be very stable. When the single-frequency (E1) GIOVE phase and code data of different receiver types are a priori corrected for the differential ISBs, the short-baseline instantaneous ambiguity success rate increases significantly and becomes comparable to the success rate of mixed GPS + GIOVE ambiguity resolution based on identical receiver types.     

Combined GPS and BDS for single-frequency continuous RTK positioning through real-time estimation of differential inter-system biases

Double-differenced (DD) ambiguities between overlapping frequencies from different GNSS constellations can be fixed to integers if the associated differential inter-system biases (DISBs) are well known. In this case, only one common pivot satellite is sufficient for inter-system ambiguity resolution. This will be beneficial to ambiguity resolution (AR) and real-time kinematic (RTK) positioning especially when only a few satellites are observed. However, for GPS and current operational BDS-2, there are no overlapping frequencies. Due to the influence of different frequencies, the inter-system DD ambiguities still cannot be fixed to integers even if the DISBs are precisely known. In this contribution, we present an inter-system differencing model for combined GPS and BDS single-frequency RTK positioning through real-time estimation of DISBs. The stability of GPS L1 and BDS B1 DISBs is analyzed with different receiver types. Along with parameterization and using the short-term stability of DISBs, the DD ambiguities between GPS and BDS pivot satellites and the between-receiver single-difference ambiguity of the GPS pivot satellite can be estimable jointly with the differential phase DISB term from epoch to epoch. Then the inter-system differencing model can benefit from the near time-constant DISB parameters and thus has better multi-epoch positioning performance than the classical intra-system differencing model. The combined GPS and BDS single-frequency RTK positioning performance is evaluated with various simulated satellite visibilities. It will be shown that compared with the classical intra-system differencing model, the proposed model can effectively improve the positioning accuracy and reliability, especially for severely obstructed situations with only a few satellites observed.     

多GNSS系统精密定轨ISB/IFB估计及特性分析

多全球导航卫星系统（Global Navigation Satellite System,GNSS）系统联合精密定轨需要考虑系统间及频率间偏差的影响。推导了多GNSS定轨系统间偏差（inter system bias,ISB）/频率间偏差（inter frequency bias,IFB）解算模型,以GPS系统硬件延迟为基准,给出了一种消除ISB/IFB秩亏的约束方法。试验数据结果表明,各系统ISB/IFB均表现出良好的稳定性及同一系统各卫星时间序列的一致性,BDS ISB的标准差为0.36 ns,Galileo ISB的标准差为0.18 ns,GLONASS IFB的标准差为0.51 ns;在接收机类型相同的情况下,不同跟踪站的ISB比较接近,但仍可达到ns级差异;GLONASS IFB在同一跟踪站相同频道号的卫星及不同跟踪站相同频道号卫星均表现出了良好的一致性。     

GPS,Galileo, QZSS and IRNSS differential ISBs: estimation and application

Knowledge of inter-system biases (ISBs) is essential to combine observations of multiple global and regional navigation satellite systems (GNSS/RNSS) in an optimal way. Earlier studies based on GPS, Galileo, BDS and QZSS have demonstrated that the performance of multi-GNSS real-time kinematic positioning is improved when the differential ISBs (DISBs) corresponding to signals of different constellations but transmitted at identical frequencies can be calibrated, such that only one common pivot satellite is sufficient for inter-system ambiguity resolution at that particular frequency. Recently, many new GNSS satellites have been launched. At the beginning of 2016, there were 12 Galileo IOV/FOC satellites and 12 GPS Block IIF satellites in orbit, while the Indian Regional Navigation Satellite System (IRNSS) had five satellites launched of which four are operational. More launches are scheduled for the coming years. As a continuation of the earlier studies, we analyze the magnitude and stability of the DISBs corresponding to these new satellites. For IRNSS this article presents for the first time DISBs with respect to the L5/E5a signals of GPS, Galileo and QZSS for a mixed-receiver baseline. It is furthermore demonstrated that single-frequency (L5/E5a) ambiguity resolution is tremendously improved when the multi-GNSS observations are all differenced with respect to a common pivot satellite, compared to classical differencing for which a pivot satellite is selected for each constellation.     

BDS–GPS inter-system bias of code observation and its preliminary analysis

Navigation applications will benefit significantly from the improved reliability, availability, and accuracy offered by combining BeiDou Navigation Satellite System (BDS) and Global Positioning System (GPS). In the BDS/GPS navigation data fusion model, the effect of inter-system bias (ISB) must be considered. We present a detailed analysis of the pseudorange measurements for BDS and GPS and demonstrate the existence of code ISB in BDS/GPS measurements. The ISB mainly consists of the time system offset, the coordinate system difference, and the inter-system hardware delay bias. A method based on statistical hypothesis testing is proposed to assess the stability and difference of the BDS–GPS ISB. Real data from 18 stations equipped with six types of receivers are used to compute the ISB. The results illustrate that (a) receiver-dependent ISBs are evident and comparatively consistent, with the maximum ISB observed in our experiments being ?1516 ns, (b) these receiver-dependent ISBs exhibit great stability in terms of their standard deviation and intra-day variation, and (c) the estimated ISBs for each BDS satellite type with respect to GPS are consistent.     

Estimation of satellite position,clock and phase bias corrections

Precise point positioning with integer ambiguity resolution requires precise knowledge of satellite position, clock and phase bias corrections. In this paper, a method for the estimation of these parameters with a global network of reference stations is presented. The method processes uncombined and undifferenced measurements of an arbitrary number of frequencies such that the obtained satellite position, clock and bias corrections can be used for any type of differenced and/or combined measurements. We perform a clustering of reference stations. The clustering enables a common satellite visibility within each cluster and an efficient fixing of the double difference ambiguities within each cluster. Additionally, the double difference ambiguities between the reference stations of different clusters are fixed. We use an integer decorrelation for ambiguity fixing in dense global networks. The performance of the proposed method is analysed with both simulated Galileo measurements on E1 and E5a and real GPS measurements of the IGS network. We defined 16 clusters and obtained satellite position, clock and phase bias corrections with a precision of better than 2 cm.     

一种无须变换参考星的GNSS单基线卡尔曼滤波算法

处理单基线全球导航卫星系统观测数据可获取位置、时间、大气延迟等信息,其应用包括相对定位、时频传递等。为实现实时性,常采用卡尔曼滤波递归地估计各类参数;为确保可靠性,还需形成一组独立的双差模糊度,并将其正确地固定为整数。实践中,滤波函数模型较常采用双差观测方程(即双差滤波模型)。若在当前历元原先的参考星不再可视时,双差滤波模型则需要定义新的参考星,并"映射"双差模糊度预报值以确保滤波连续。此外,双差滤波模型所计算的接收机相位钟差估值吸收了对应于参考星的站间单差模糊度,因此当参考星变换后可能会发生"整周跳跃"。在仍将双差模糊度作为一类可估参数的前提下,本文推导出以站间单差观测方程为滤波函数模型的算法(单差滤波模型),并证明了其与双差滤波模型具备理论上的等价性和实施上的差异性。与双差滤波模型相比,单差滤波模型不再需要"映射"双差模糊度预报值等运算,从而具备了更高的计算效率和灵活性;单差滤波模型所提供的接收机相位钟差估值也不受"整周跳跃"的影响,因此特别有利于频率传递应用。     

A cycle slip fixing method with GPS + GLONASS observations in real-time kinematic PPP

A key limitation of precise point positioning (PPP) is the long convergence time, which requires about 30 min under normal conditions. Frequent cycle slips or data gaps in real-time operation force repeated re-convergence. Repairing cycle slips with GPS data alone in severely blocked environments is difficult. Adding GLONASS data can supply redundant observations, but adds the difficulty of having to deal with differing wavelengths. We propose a single-difference between epoch (SDBE) method to integrate GPS and GLONASS for cycle slip fixing. The inter-system bias can be eliminated by SDBE, thus only one receiver clock parameter is needed for both systems. The inter-frequency bias of GLONASS satellites also cancels in the SDBE, so cycle slips are preserved as integers, and the LAMBDA method is adopted to search for cycle slips. Data from 7 days of 20 globally distributed IGS sites were selected to test the proposed cycle slip fixing procedure with artificial blocking of the signal; cycle slips were introduced for all un-blocked satellites at each epoch. For a 30-s sampling interval, the average success rate of fixing can be improved from 73 to 98 % by adding GLONASS. Even for a 180-s sampling interval, GPS + GLONASS can achieve a success rate of 81 %. A real-time kinematic PPP experiment was also performed, and the results show that using GPS + GLONASS can achieve continuous high-accuracy real-time PPP without re-convergence.     

长距离GPS/BDS参考站网多频载波相位整周模糊度解算方法

长距离网络RTK是实现GPS/BDS高精度实时定位的主要手段之一,其核心是长距离参考站网GPS/BDS整周模糊度的快速准确确定。本文提出了一种长距离GPS/BDS参考站网载波相位整周模糊度解算方法,首先利用GPS双频观测数据计算和确定宽巷整周模糊度,同时利用BDS的B2、B3频率观测值确定超宽巷整周模糊度。然后建立GPS载波相位整周模糊度和大气延迟误差的参数估计模型,附加双差宽巷整周模糊度的约束,解算双差载波相位整周模糊度,并建立参考站网大气延迟误差的空间相关模型。根据B2、B3频率的超宽巷整周模糊度建立包含大气误差参数的载波相位整周模糊度解算模型,利用大气延迟误差空间相关模型约束BDS双差载波相位整周模糊度的解算。克服了传统的使用无电离层组合值解算整周模糊度的不利影响。采用实测长距离CORS网GPS、BDS多频观测数据进行算法验证,试验结果证明该方法可实现长距离参考站网GPS/BDS载波相位整周模糊度的准确固定。     

SPODS软件GPS/GNSS网解的模糊度解算方法

为了尽可能多地固定模糊度,需要定义一组"最容易被固定"的独立双差模糊度,当前最优的方法(传统方法)是分基线层和网层对候选双差模糊度进行独立性检验,其中候选模糊度按照其综合固定成功概率进行排列。考虑到通常的网解中,测站数量远多于卫星数量,提出了先星座层再网层的分层独立模糊度选取方法以及基于更新协方差阵上三角平方根的序贯模糊度固定方法,并应用于西安测绘研究所自主设计和开发的SPODS软件中。采用包含64个IGS监测站的实测GPS数据进行单天解算试验,验证了该方法的正确性。试验结果表明,该方法和传统方法得到的单天GPS轨道解与IGS最终综合轨道比较的1DRMS都约为0.012m,独立双差模糊度成功固定的比例约为92%,两种方法仅有非常细微的差异。不同测站数量的进一步试验表明,该方法和传统方法在网层需要采用Grams-Schmidt方法进行独立性检验的候选双差模糊度数量的比约等于卫星和测站的数量比,这与理论分析结论一致。对于实际应用中更为普遍的测站数量大于卫星数量的情况,该方法将减少独立双差模糊度选取的计算时间,且测站数越多,其优势越明显。     

BDS-3和Galileo组合的RTK定位性能分析

北斗三号卫星导航系统（BeiDou-3 navigation satellite system, BDS-3）已全面建成并向全球用户提供可靠的定位、导航和授时（positioning, navigation and timing, PNT）服务。为了实现与其他全球卫星导航系统（global navigation satellite system, GNSS）的兼容性和互操作性,BDS-3在BDS-2的基础上调制了B1C和B2a两个新信号,与伽利略系统（Galileo）的E1和E5a实现了频率的复用。系统间偏差（inter-system bias, ISB）对于实现不同GNSS之间的融合处理至关重要,为此提出了基于单差模型的ISB估计与应用算法,并对BDS-3与Galileo重叠频率之间的ISB进行了分析。基于可跟踪BDS-3新信号的几类接收机,揭示了BDS-3和Galileo之间的ISB的特性,在此基础上分析了BDS-3和Galileo组合的实时动态（real-time kinematic, RTK）定位性能。结果表明,基于相同类型的接收机B1C-E1和B2a-E5a之间是不存在ISB...     

Performance analysis of integrated GPS/GLONASS carrier phase-based positioning

Due to the different signal frequencies for the GLONASS satellites, the commonly-used double-differencing procedure for carrier phase data processing can not be implemented in its straightforward form, as in the case of GPS. In this paper a novel data processing strategy, involving a three-step procedure, for integrated GPS/GLONASS positioning is proposed. The first is pseudo-range-based positioning, that uses double-differenced (DD) GPS pseudo-range and single-differenced (SD) GLONASS pseudo-range measurements to derive the initial position and receiver clock bias. The second is forming DD measurements (expressed in cycles) in order to estimate the ambiguities, by using the receiver clock bias estimated in the above step. The third is to form DD measurements (expressed in metric units) with the unknown SD integer ambiguity for the GLONASS reference satellite as the only parameter (which is constant before a cycle slip occurs for this satellite). A real-time stochastic model estimated by residual series over previous epochs is proposed for integrated GPS/GLONASS carrier phase and pseudo-range data processing. Other associated issues, such as cycle slip detection, validation criteria and adaptive procedure(s) for ambiguity resolution, is also discussed. The performance of this data processing strategy will be demonstrated through case study examples of rapid static positioning and kinematic positioning. From four experiments carried out to date, the results indicate that rapid static positioning requires 1 minute of single frequency GPS/GLONASS data for 100% positioning success rate. The single epoch positioning solution for kinematic positioning can achieve 94.6% success rate over short baselines (<6 km).     

PERFORMANCE ANALYSIS OF INTEGRATED GPS/GLONASS CARRIER PHASE-BASED POSITIONING

1　ＩｎｔｒｏｄｕｃｔｉｏｎＲｅａｌ＿ｔｉｍｅｋｉｎｅｍａｔｉｃＧＰＳｐｒｅｃｉｓｅｐｏｓｉｔｉｏｎｉｎｇｈａｓｂｅｅｎｐｌａｙｉｎｇａｎｉｎｃｒｅａｓｉｎｇｒｏｌｅｉｎｂｏｔｈｓｕｒｖｅｙｉｎｇａｎｄｎａｖｉｇａｔｉｏｎ ,ａｎｄｈａｓｂｅｃｏｍｅａｎｅｓｓｅｎｔｉａｌｔｏｏｌｆｏｒｐｒｅｃｉｓｅｒｅｌａｔｉｖｅｐｏｓｉｔｉｏｎｉｎｇ .Ｈｏｗｅｖｅｒ,ｒｅｌｉａｂｌｅａｎｄｃｏｒｒｅｃｔａｍｂｉｇｕｉｔｙｒｅｓｏｌｕｔｉｏｎｄｅｐｅｎｄｓｏｎｏｂｓｅｒｖａ ｔｉｏｎｓｕｐｏｎａｌａｒｇｅｎｕｍｂｅ…     

Multi-GNSS precise point positioning (MGPPP) using raw observations

A joint-processing model for multi-GNSS (GPS, GLONASS, BDS and GALILEO) precise point positioning (PPP) is proposed, in which raw code and phase observations are used. In the proposed model, inter-system biases (ISBs) and GLONASS code inter-frequency biases (IFBs) are carefully considered, among which GLONASS code IFBs are modeled as a linear function of frequency numbers. To get the full rank function model, the unknowns are re-parameterized and the estimable slant ionospheric delays and ISBs/IFBs are derived and estimated simultaneously. One month of data in April, 2015 from 32 stations of the International GNSS Service (IGS) Multi-GNSS Experiment (MGEX) tracking network have been used to validate the proposed model. Preliminary results show that RMS values of the positioning errors (with respect to external double-difference solutions) for static/kinematic solutions (four systems) are 6.2 mm/2.1 cm (north), 6.0 mm/2.2 cm (east) and 9.3 mm/4.9 cm (up). One-day stabilities of the estimated ISBs described by STD values are 0.36 and 0.38 ns, for GLONASS and BDS, respectively. Significant ISB jumps are identified between adjacent days for all stations, which are caused by the different satellite clock datums in different days and for different systems. Unlike ISBs, the estimated GLONASS code IFBs are quite stable for all stations, with an average STD of 0.04 ns over a month. Single-difference experiment of short baseline shows that PPP ionospheric delays are more precise than traditional leveling ionospheric delays.     

改进的ARCE方法及其在单频 GPS快速定位中的应用

基于TIKHONOV正则化原理，设计了一种正则化矩阵的构造方法，将ARCE(ambiguity resolution using constraint equation)方法进行了改进。通过新的正则化矩阵的作用，减弱了GPS快速定位中少数历元情形下法矩阵的病态性，得到了比较准确的模糊度浮动解，大大减小了模糊度的搜索范围，利用ARCE方法固定模糊度的成功率高。并结合一个算例，验证了本文改进方法的效果。     

Multi-GNSS phase delay estimation and PPP ambiguity resolution: GPS,BDS, GLONASS,Galileo

This paper focuses on the precise point positioning (PPP) ambiguity resolution (AR) using the observations acquired from four systems: GPS, BDS, GLONASS, and Galileo (GCRE). A GCRE four-system uncalibrated phase delay (UPD) estimation model and multi-GNSS undifferenced PPP AR method were developed in order to utilize the observations from all systems. For UPD estimation, the GCRE-combined PPP solutions of the globally distributed MGEX and IGS stations are performed to obtain four-system float ambiguities and then UPDs of GCRE satellites can be precisely estimated from these ambiguities. The quality of UPD products in terms of temporal stability and residual distributions is investigated for GPS, BDS, GLONASS, and Galileo satellites, respectively. The BDS satellite-induced code biases were corrected for GEO, IGSO, and MEO satellites before the UPD estimation. The UPD results of global and regional networks were also evaluated for Galileo and BDS, respectively. As a result of the frequency-division multiple-access strategy of GLONASS, the UPD estimation was performed using a network of homogeneous receivers including three commonly used GNSS receivers (TRIMBLE NETR9, JAVAD TRE_G3TH DELTA, and LEICA). Data recorded from 140 MGEX and IGS stations for a 30-day period in January in 2017 were used to validate the proposed GCRE UPD estimation and multi-GNSS dual-frequency PPP AR. Our results show that GCRE four-system PPP AR enables the fastest time to first fix (TTFF) solutions and the highest accuracy for all three coordinate components compared to the single and dual system. An average TTFF of 9.21 min with \(7{^{\circ }}\) cutoff elevation angle can be achieved for GCRE PPP AR, which is much shorter than that of GPS (18.07 min), GR (12.10 min), GE (15.36 min) and GC (13.21 min). With observations length of 10 min, the positioning accuracy of the GCRE fixed solution is 1.84, 1.11, and 1.53 cm, while the GPS-only result is 2.25, 1.29, and 9.73 cm for the east, north, and vertical components, respectively. When the cutoff elevation angle is increased to \(30{^{\circ }}\), the GPS-only PPP AR results are very unreliable, while 13.44 min of TTFF is still achievable for GCRE four-system solutions.     

BeiDou B2/Galileo E5b短基线紧组合相对定位模型及性能评估

随着GPS和GLONASS系统的现代化以及Galileo和BeiDou卫星导航系统的建设,GNSS正朝多频多系统的方向发展。本文对BeiDou B2/Galileo E5b短基线紧组合相对定位的模型与算法进行了研究,详细推导了BeiDou B2/Galileo E5b短基线紧组合相对定位的模型与算法,并对其定位性能进行了分析。重点分析了BeiDou B2与Galileo E5b频点的接收机间差分系统间偏差的长期稳定性,结果表明:基线两端的接收机类型(包括固件版本)相同时,差分系统间偏差接近于0;基线两端的接收机类型不同时,差分系统间偏差较大,但具有长期稳定性,因此能够事先标定并作为改正数用于后续的定位中。最后基于BeiDou B2/Galileo E5b单频单历元相对定位试验对系统间紧组合模型的定位效果进行了比较验证。结果表明,相对于传统的松组合模型,使用改正系统间偏差的紧组合模型能够显著提高模糊度固定的成功率,尤其是在遮挡比较严重、单系统可观测到的卫星数较少的情况下,模糊度固定成功率可以提高10%~25%。     

Assessment of PPP integer ambiguity resolution using GPS,GLONASS and BeiDou (IGSO,MEO) constellations

Although integer ambiguity resolution (IAR) can improve positioning accuracy considerably and shorten the convergence time of precise point positioning (PPP), it requires an initialization time of over 30 min. With the full operation of GLONASS globally and BDS in the Asia–Pacific region, it is necessary to assess the PPP–IAR performance by simultaneous fixing of GPS, GLONASS, and BDS ambiguities. This study proposed a GPS + GLONASS + BDS combined PPP–IAR strategy and processed PPP–IAR kinematically and statically using one week of data collected at 20 static stations. The undifferenced wide- and narrow-lane fractional cycle biases for GPS, GLONASS, and BDS were estimated using a regional network, and undifferenced PPP ambiguity resolution was performed to assess the contribution of multi-GNSSs. Generally, over 99% of a posteriori residuals of wide-lane ambiguities were within ±0.25 cycles for both GPS and BDS, while the value was 91.5% for GLONASS. Over 96% of narrow-lane residuals were within ±0.15 cycles for GPS, GLONASS, and BDS. For kinematic PPP with a 10-min observation time, only 16.2% of all cases could be fixed with GPS alone. However, adding GLONASS improved the percentage considerably to 75.9%, and it reached 90.0% when using GPS + GLONASS + BDS. Not all epochs could be fixed with a correct set of ambiguities; therefore, we defined the ratio of the number of epochs with correctly fixed ambiguities to the number of all fixed epochs as the correct fixing rate (CFR). Because partial ambiguity fixing was used, when more than five ambiguities were fixed correctly, we considered the epoch correctly fixed. For the small ratio criteria of 2.0, the CFR improved considerably from 51.7% for GPS alone, to 98.3% when using GPS + GLONASS + BDS combined solutions.     

BDS/GPS融合定位中系统偏差补偿模型及其性能分析

多系统的融合定位可有效提高用户导航定位的连续性、可靠性及定位精度。针对BDS、GPS观测量间存在系统间偏差的实际情况,建立了顾及系统误差的BDS/GPS融合定位模型,即在函数模型中增加附加参数来吸收系统间偏差,构造了新的顾及先验信息的融合定位模型,分析了这种新融合模型的特点及其对定位结果的影响。利用不同品牌接收机在中国不同地域对新的融合模型进行试验,试验结果表明:BDS、GPS观测量存在系统间偏差,且不同接收机的系统间偏差量值并不一样;增加系统参数的融合定位模型能较好地吸收BDS、GPS观测量的系统间偏差的影响,改善其融合导航定位性能;在观测卫星数不足、单系统不能定位的情况下,考虑先验信息的融合定位模型仍能获得较好的定位结果。