考虑随机模型精化的精密GPS动态定位新方法

GPS动态定位要求建立函数模型和随机模型。函数模型描述的是观测值和待估参数之间的物理和几何关系，随机模型描述了GPS观测值的统计特征，并通过观测值的方差协方差给定了每个观测值对最后的定位结果的贡献。正确给定函数模型和随机模型对于GPS定位结果的估计和观测值的粗差探测均至关重要。由于有各种误差存在于伪距和载波相位观测值中，一般GPS动态定位模型均采用双差观测值来构建函数模型。有时候，仔细地使用单差观测值，较之双差观测值有更多的优点，给出了选用单差观测值的理由。但是单差观测值给函数模型带来了接收机钟差，如果直接使用单差观测方程，设计矩阵是奇异的。为了解决这个问题，将伪距观测值中接收机钟差项和接收机延迟项合并为一个新的未知参数。至于载波相位观测值，首先选定一个参考卫星，然后在观测方程的右端同时增加一正一负的参考卫星单差整周模糊度，将正项与接收机钟差项和接收机延迟项合并为一个新的未知参数，将负项和原观测方程中的单差整周模糊度项合并为双差整周模糊度，而参考卫星观测方程的模糊度项则为零，这样无须组建双差观测值，软件实现较容易，也可以直接使用LAMBDA法求整周模糊度，最终也解决了观测方程奇异的问题。准确理解观测值的统计特征是建立GPS随机模型的基础，长期以来GPS商业软件均采用简化模型。关于GPS随机模型的研究远没有函数模型那样受到广泛关注，静态GPS定位可以采用方差协方差分量估计等严密的方法，而动态定位无法承担方差协方差分量估计的计算负担。GPS观测值的信噪比(SNR)是GPS接收机观测过程中的副产品，影响SNR值的因素，如大气层、多路径、接收机内部电路等，也正好是GPS观测值的误差源，因此GPS观测值的方差与SNR存在一定的对应关系。利用这个对应关系来精化GPS随机模型。为了验证本文采用的函数模型的正确性和随机模型的有效性，我们对1999年的一次实测数据(包括零基线和短基线)进行了试算。与零基线的真值和GPSurvey 2．35处理的短基线静态结果比较，表明使用的函数模型是正确的。简化随机模型和精化随机模型处理的结果比较说明精化模型提高了基线处理的精度，同时说明了研究GPS随机模型精化的必要性。     

一种整数相位钟法精密单点定位模糊度固定模型

星间单差法是常用的精密单点定位PPP模糊度固定方法,但是要面临基准星转换的问题。为此,提出了一种逐级模糊度固定模型,采用法国CNES发布的整数相位钟差产品,在PPP非差观测模型基础上逐一选取两颗卫星进行模糊度固定;得到多组单差模糊度固定解后,再以此构成约束条件进行滤波得到其他参数。实验选取了8个IGS站共48个观测时段进行模糊度固定实验。结果表明,模糊度成功固定后,位置三维误差平均值由5.60 cm减小到2.72 cm;位置误差标准差由3.64 cm减小到1.50 cm。仿动态条件下,模糊度固定后位置误差由6.02 cm降至4.75 cm。     

一种GPS整周模糊度单历元解算方法

仅利用单历元的载波相位观测值进行整周模糊度解算,观测方程秩亏,给单历元模糊度解算带来很大困难.因此,本文提出一种单历元确定GPS整周模糊度的方法.利用单历元测码伪距观测值和双频载波相位观测值组成双差观测方程,根据方差矩阵对宽巷模糊度进行分组,采用基于LABMDA方法的逐步解算方法来确定双差相位观测值的宽巷模糊度.确定宽...     

长距离网络RTK基准站间整周模糊度单历元确定方法

提出一种长距离(100～200km)网络RTK基准站间的整周模糊度单历元确定方法。该方法首先利用载波相位模糊度间的线性约束关系对双差宽巷模糊度进行搜索。为了减小非弥散误差残差对载波相位模糊度解算的影响,采用了一种新的根据高度角重新选择基准卫星的方法。然后根据双差宽巷模糊度选取双频载波相位模糊度的备选组合,利用基准站间非弥散误差残差的计算值对双差载波相位模糊度进行搜索和确定。经试验算例的验证,该方法快速、稳定,不受周跳影响,只需一个历元的观测数据即可确定长距离基准站间的双差整周模糊度。     

网络RTK流动站整周模糊度的单历元解算

传统网络RTK模糊度解算方法需要多个历元的观测数据,并且要进行周跳的探测和修复,影响模糊度解算的效率。本文提出一种单历元确定网络RTK双差整周模糊度的新方法。首先利用测码伪距观测值和载波相位观测值的单历元数据组成双差联合观测方程,采用改进LABMDA算法进行两步搜索确定GPS双差相位观测值的宽巷模糊度。确定宽巷模糊度后,再用宽巷模糊度值和载波相位观测值组成新的联合观测方程,大大改善了方程的状态,可以准确解算出GPS双差整周模糊度,显著提高了网络RTK整周模糊度固定的效率。     

GPS精密卫星钟差估计研究

随着精密单点定位技术的发展,对于精确的卫星坐标以及卫星钟差改正精度的要求越来越高,精密卫星星历以及精密卫星钟差的求解成为制约精密单点定位技术发展的瓶颈。本文基于修复周跳的载波相位观测值与相位平滑伪距观测值,采用无电离层延迟星间单差精密卫星钟差估计模型,在先估计出整周模糊度后,进行了精密卫星钟差的估计,并采用与IGS事后精密钟差作二次差的方法进行精度分析,这对于提高精密单点定位精度具有一定的意义。     

BDS网络RTK中距离参考站整周模糊度单历元解算方法

提出了一种BDS网络RTK中距离(50~100 km)参考站间的双频载波相位整周模糊度单历元解算方法。该方法首先利用B1、B2载波相位整周模糊度间的线性关系选取B1、B2载波相位整周模糊度备选值。利用双频载波相位整周模糊度备选值计算双差电离层延迟误差,根据参考站各卫星电离层延迟误差间的空间关系,使用双差电离层延迟误差构建双差电离层延迟误差的线性计算模型。通过双差电离层延迟误差线性计算模型的建立搜索和确定B1、B2载波相位的整周模糊度。经CORS网实测数据试验算例的验证,该方法只需一个历元的观测数据即可确定参考站间双差B1、B2载波相位整周模糊度,且不受周跳影响。     

零基线约束的参考站间模糊度固定方法研究

论文提出一种GNSS零基线约束的参考站间模糊度快速固定方法,该方法基于“一天线+双接收机”设计进行参考站间模糊度固定。首先,利用两个参考站四台接收机观测值形成四对参考站间无电离层组合双差观测值。其次,结合无电离层组合观测值模糊度和宽巷整周模糊度进行卡尔曼滤波和N1模糊度实数解计算。然后在模糊度域内进行线性变换,单历元固定零基线模糊度;利用零基线模糊度约束解算四组参考站间模糊度中的一组模糊度。最后,利用固定的模糊度还原所有原始双差整周模糊度。实验证明:与传统参考站间模糊度固定方法比,充分利用了零基线模糊度固定的约束条件进行参考站间模糊度固定,加快了模糊度收敛速度,极大地提高了模糊度固定的速度和成功率。      

RDSS动基线整周模糊度求解方法及仿真

提出了一种利用动基线求解卫星无线电测定业务RDSS(Radio Determination Satellite Service)定向技术中单差整周模糊度的方法.接收机首先在基线处于某一初始状态下测量两颗卫星的载波相位单差,然后旋转基线两个较小角度,并分别测量不同状态下两颗卫星的载波相位单差.通过两次旋转获得的各基线状态之间的关系获得旋转平面,进而求得初始状态下的基线矢量.将求得的基线矢量初值代入单差观测方程可求得整周模糊度.整周模糊度的求解精度与基线长度、基线初始状态方位角以及基线转角等多种因素有关.仿真分析表明,当基线长度为3 m,基线两次旋转角度分别为30°,基线初始状态方位角在(0°,120°)和(180°,300°)范围内时,都可以正确求解整周模糊度.     

长距离网络RTK参考站间双差模糊度快速解算算法

提出了一种长距离网络RTK参考站间双差整周模糊度快速解算方法,该方法利用双频载波相位模糊度间的线性关系确定宽巷模糊度,然后选取双频载波相位的备选模糊度组合,通过计算参考站间对流层误差和轨道误差等非色散误差,对双频载波相位整周模糊度进行搜索。实验结果表明,此方法能够快速、可靠地解算长距离参考站间的双差整周模糊度。     

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参考网数据处理方法

全球导航卫星系统（GNSS）参考网多用于估计卫星轨道/钟差、监测地表形变和速度场、确定精密地球自转参数等方面。相关数据处理模式包括：双差基线解（DD）和非差精密单点定位（PPP）等。本文首先从GNSS基本观测方程出发,通过选取两组基准参数,导出了上述两模式下的列满秩观测方程,然后分析了它们的不足,例如：相位偏差在DD模式中吸收了钟差,丧失了时不变特性;模糊度在PPP模式中吸收了相位偏差,失去了整数性。基于上述分析,本文提出了一种新的参考网数据处理方案,以充分融合DD和PPP模式的优势。它的关键策略是精选基准参数,以达到消秩亏的目的,具体优点体现在：相位偏差独立可估,若合理约束为时不变参数,可充分减少参数个数,提高网解精度;待估模糊度具备整周特性,经由模糊度固定,可改善网解可靠性。     

卫星钟差解算及其星间单差模糊度固定

整数相位模糊度解算可以显著提高GNSS精密单点定位（PPP）的精度。本文提出一种解算卫星钟差的方法,通过固定星间单差模糊度恢复出能够支持单台接收机进行整数模糊度解算的卫星钟差,即所谓的“整数”钟差。为了实现星间单差模糊度固定,分别通过卫星端宽巷FCB解算和模糊度基准的选择与固定恢复出宽巷和窄巷模糊度的整数性质。为了证明本文方法的可行性,采用IGS测站的GPS数据进行卫星钟差解算试验。结果表明,在解算钟差时,星间单差模糊度固定的平均成功率为73%。得到的卫星钟差与IGS最终钟差产品相比,平均的RMS和STD分别为0.170和0.012 ns。448个IGS测站的星间单差宽巷和窄巷模糊度小数部分的分布表明本文得到的卫星钟差和FCB产品具备支持PPP用户进行模糊度固定的能力。基于以上产品开展了模拟动态PPP定位试验,结果表明模糊度固定之后,N、E、U和3D的定位精度（RMS）分别达到0.009、0.010、0.023和0.027 m,与不固定模糊度或采用IGS钟差的结果相比,分别提高了30.8%、61.5%、23.3%和37.2%。     

载波相位差分接收机自主完好性监测研究

首先提出了浮点变换完全去相关法,该方法能够在单历元动态确定整周模糊度。研究了基于载波相位测量的完好性监测方法。利用最小二乘残差构造统计检验量,对整周模糊度进行检测。分析了定位误差保护限与卫星构型、漏警概率的关系。实测数据表明,整周模糊度在单历元动态求解的成功率为100%,增加1颗卫星将使垂直定位误差保护限减少约0.2 m,统计检验量检测周跳的正确率为100%。     

星间单差模糊度固定的低轨卫星精密定轨精度分析

高精度、高可靠性的卫星轨道是实现低轨卫星精密应用的重要前提,而模糊度固定技术是提高卫星定轨精度的关键途径。研究了基于整数钟的星间单差模糊度固定原理和方法,并利用2019年4月—5月的两颗GRACE-FO（gravity recovery and climate experiment follow on）卫星数据（GRACE-C/D）系统评估了固定解对低轨卫星简化动力学和运动学定轨的精度提升效果。结果表明,两颗卫星简化动力学和运动学定轨的宽巷模糊度固定率均达到99%,而窄巷模糊度固定率在95%左右。对于简化动力学定轨,GRACE-C/D固定解轨道的重叠轨道的3D均方根误差（root mean square error, RMSE）分别从7.1 mm和7.4 mm减小到了4.2 mm和3.6 mm;卫星激光测距（satellite laser ranging, SLR）残差标准差（standard deviation, STD）分别从15.9 mm和14.4 mm降低到了10.8 mm和11.0 mm,精度提升了32%和24%;K波段测距残差RMSE从8.0 mm减小到2.9 mm,进一步表明固定解还能有效提升低轨卫星间相对位置精度。对于运动学定轨,与精密科学轨道产品互差3D RMSE,浮点解分别为37.5 mm和36.4 mm,固定解分别为27.7 mm和25.5 mm,精度提升约28%,SLR残差STD也减小了约20%。     

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　ＩｎｔｒｏｄｕｃｔｉｏｎＲｅａｌ＿ｔｉｍｅｋｉｎｅｍａｔｉｃＧＰＳｐｒｅｃｉｓｅｐｏｓｉｔｉｏｎｉｎｇｈａｓｂｅｅｎｐｌａｙｉｎｇａｎｉｎｃｒｅａｓｉｎｇｒｏｌｅｉｎｂｏｔｈｓｕｒｖｅｙｉｎｇａｎｄｎａｖｉｇａｔｉｏｎ ,ａｎｄｈａｓｂｅｃｏｍｅａｎｅｓｓｅｎｔｉａｌｔｏｏｌｆｏｒｐｒｅｃｉｓｅｒｅｌａｔｉｖｅｐｏｓｉｔｉｏｎｉｎｇ .Ｈｏｗｅｖｅｒ,ｒｅｌｉａｂｌｅａｎｄｃｏｒｒｅｃｔａｍｂｉｇｕｉｔｙｒｅｓｏｌｕｔｉｏｎｄｅｐｅｎｄｓｏｎｏｂｓｅｒｖａ ｔｉｏｎｓｕｐｏｎａｌａｒｇｅｎｕｍｂｅ…     

The least-squares ambiguity decorrelation adjustment: its performance on short GPS baselines and short observation spans

The least-squares ambiguity decorrelation adjustment is a method for fast GPS double-difference (DD) integer ambiguity estimation. The performance of the method will be discussed, and although it is stressed that the method is generally applicable, attention is restricted to short-baseline applications in the present contribution. With reference to the size and shape of the ambiguity search space, the volume of the search space will be introduced as a measure for the number of candidate grid points, and the signature of the spectrum of conditional variances will be used to identify the difficulty one has in computing the integer DD ambiguities. It is shown that the search for the integer least-squares ambiguities performs poorly when it takes place in the space of original DD ambiguities. This poor performance is explained by means of the discontinuity in the spectrum of conditional variances. It is shown that through a decorrelation of the ambiguities, transformed ambiguities are obtained which generally have a flat and lower spectrum, thereby enabling a fast and efficient search. It is also shown how the high precision and low correlation of the transformed ambiguities can be used to scale the search space so as to avoid an abundance of unnecessary candidate grid points. Numerical results are presented on the spectra of conditional variances and on the statistics of both the original and transformed ambiguities. Apart from presenting numerical results which can typically be achieved, the contribution also emphasizes and explains the impact on the method's performance of different measurement scenarios, such as satellite redundancy, single vs dual-frequency data, the inclusion of code data and the length of the observation time span. Received: 31 October 1995 / Accepted: 21 March 1997     

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.     

On optimal filtering of GPS dual frequency observations without using orbit information

The concept of optimal filtering of observations collected with a dual frequency GPS P-code receiver is investigated in comparison to an approach for C/A-code units. The filter presented here uses only data gathered between one receiver and one satellite. The estimated state vector consists of a one-way pseudorange, ionospheric influence, and ambiguity biases. Neither orbit information nor station information is required. The independently estimated biases are used to form double differences where, in case of a P-code receiver, the wide lane integer ambiguities are usually recovered successfully except when elevation angles are very small. An elevation dependent uncertainty for pseudorange measurements was discovered for different receiver types. An exponential model for the pseudorange uncertainty was used with success in the filter gain computations.