两步法快速解算编队卫星GPS模糊度

为克服卫星编队飞行实时相对定位中双频模糊度解算速度慢的缺点,结合扩展Kalman滤波(EKF,Extended Kalman Filter),首先采用少数个历元(如10个)相位平滑伪距相对定位结果与 L₆的平均值对滤波初始化,再根据两步法解算双频模糊度,即先解算并正确固定宽巷模糊度,获得较准确的基线分量估值,然后采用选权拟合方法,将基线分量作为约束条件解算并固定双频模糊度.仿真算例计算结果表明,当宽巷模糊度正确固定后,编队卫星间相对定位误差在5cm以内,两步法可以在较短时间(约3min)内固定双频模糊度,为精确解算编队卫星的相对状态提供保障.     

基于固定模糊度PPP的低轨卫星实时精密定轨

传统动力学定轨法受制于动力学模型精度，传统几何定轨法精度受限，只能达到亚米级，而基于精密单点定位（PPP）模式的几何定轨法一般采用浮点解，定轨精度及可靠性较基于双差模式的相对定位较差。为提高PPP模式低轨定轨的定位性能，利用中国区域内外的IGS测站计算出当前所有卫星的宽巷和窄巷相位小数偏差产品，对经过中国大陆区域上空的国产低轨卫星海洋二号（HY-2）和资源三号 (ZY-3) 卫星进行固定模糊度PPP的定轨解算，与事后精密轨道结果进行比较，分析其外符合精度。结果表明：仅利用约10min弧段的HY-2和ZY-3卫星数据，切向与径向的定轨精度可达2cm左右，法向为5cm左右，较浮点解定轨精度大幅提升。基于固定模糊度PPP的定轨方法能够满足厘米级的实时精密定轨。     

基于MLAMBDA的BDS/GPS单频RTK性能分析

GNSS RTK技术以其高精度、高效率、实时性的优点，被广泛应用于航空航天等领域.目前双频RTK技术已非常成熟并且应用较广.相比于双频，单频GNSS RTK在数据质量控制、定位误差处理等方面存在难点.因此单频RTK服务精度可能会受到限制，其定位性能有待研究.本文基于扩展卡尔曼滤波模型，通过MLAMBDA模糊度搜索方法和Ratio检验法，结合实测数据，对比分析BDS，GPS，BDS/GPS三种模式下的单频RTK定位性能.实验证明在静态场景下，三种模式的单频RTK定位精度都在厘米级，可满足高精度定位需求；动态场景下三种模式的模糊度固定率都在70%以上，可满足日常定位需求.在静态及动态应用场景下，北斗的模糊度固定率最高，模糊度解算所用时间短，能实现快速RTK定位.      

Kalman滤波估算电离层延迟的一种优化方法

频间偏差（Inter Frequency Bias，IFB）通常会给电离层延迟的解算带来误差.目前从电离层延迟中消除频间偏差的方法是基于GPS双频观测数据建立垂直电离层模型，利用卡尔曼滤波实时估算电离层模型系数和频间偏差.然而滤波过程中的测量噪声协方差矩阵没有考虑系统观测量之间的相关性，这会导致滤波模型不准确，进而影响最后求解的电离层延迟的准确性.本文选取了美国19个参考站的GPS双频观测数据，利用卡尔曼滤波实时估算电离层模型系数以及频间偏差.在滤波过程中，通过将先验频间偏差的估计方差引入测量噪声方差，实现对测量噪声协方差矩阵的优化.计算结果表明，优化后得到的卫星频间偏差与欧洲定轨中心（Center for Orbit Determination in Europe，CODE）得到的频间偏差更接近.将优化后的电离层延迟代入伪距解算，得到的位置误差的标准差在东向和天顶向分别下降了12.5%和15.4%，天顶向误差平均值下降了17.6%，定位精度得到提高.      

基于实时SSR数据的LEO卫星精密单点定位研究

对目前低轨卫星实时定位的方法进行了研究,现在通常采用GPS定位,使用广播星历和普通晶振,实时定位精度一般在10m以内,不能满足高精度实时定位的需求。IGS组织在全球范围内对GPS跟踪分析,生成精密星历、精密钟差产品、按SSR格式的广播星历和钟差修正产品并在网上发布。对这些IGS产品进行了调查,提出在现有测控支持情况下,可以通过高密度上注SSR信息流实现在轨高精度定位。以某型号低轨微小卫星在轨导航增强载荷为应用背景,用IGS03产品中的1057和1058数据对双频GPS接收机的星历和钟差进行修正,采用递推最小二乘估计和LAMDA模糊度固定过对载波相位和伪距信息进行处理,在短时间内获得亚米级定位结果。     

新的GPS单频单历元定姿算法

LAMBDA算法依赖于初始的模糊度浮点解,但仅用载波相位观测方程需要多个历元才能获得浮点解,将导致初始化时间过长.针对这一问题,对GPS(Global Positioning System)单历元的载波相位单差方程进行特殊变换,将未知的整周模糊度看成噪声,从而构造出新的观测方程,和原始的观测方程进行组合求解,克服了仅用载波相位双差观测方程因为亏秩而无法在单历元获得浮点解的缺点,解决了初始化时间的问题.通过深入研究浮点解和固定解之间的关系,提出一种将低精度浮点解映射到固定解的方法,降低了LAMBDA算法对高精度浮点解的依赖性,避免了用多个历元获取浮点解的高精度,从而实现了单频、单历元的整周模糊度估计.通过实际测试,该算法成功率高于97%,能够有效地用于实时动态姿态解算.      

基于增强型单频GPS的高精度星间相对定位样条方法

为了实现高精度的分布式SAR星间相对定位,提出了一种基于增强型单频GPS的相对定位样条方法,即在单频GPS测量信息的基础上,增加星间距离观测信息,并结合相对位置参数的连续特性,建立了相对定位的样条模型,最后利用最小二乘法进行参数估计.仿真结果表明,新方法不仅能大大提高相对定位精度,而且还能有效地减少固定整周模糊度所需的历元.最后的理论分析证明了仿真的正确性.      

精密单点定位部分模糊度固定方法

提出一种新的模糊度固定方法,并把算法应用到精密单点定位(Precise Point Positioning,PPP)中,将模糊度进行选择固定,以实现恶劣环境下快速、有效获得较高PPP固定解的目的。该方法以模糊度固定成功率和Ratio因子为约束条件,基于最小方差优先固定的模糊度子集选取方法实现PPP中模糊度的部分固定。算例从模糊度固定可靠性、收敛速度以及PPP固定解精度等方面对部分模糊度固定方法进行了分析。结果表明,该方法使得PPP模糊度固定的平均可靠性比传统固定方法提高了近一倍、参数解算的平均收敛时间缩短了1/4左右。因此,更适合在初始阶段观测信息较少或残余误差较大而难以固定全部模糊度的时候,快速获得较高精度的PPP固定解。     

可容忍信标误差的三维传感网节点定位方法

针对信标位置存在误差情况下的三维无线传感器网络节点定位问题,提出一种基于正交回归的多跳定位方法.同时考虑到自变量误差和因变量偏差对节点坐标估计的影响,基于约束加权正交回归参数估计准则,建立可容忍信标位置误差的三维多跳定位模型,解决了信标位置和距离估计两方面的误差并存时的节点自定位问题,并给出求解节点坐标最优值的数值方法;推导出相应的坐标估计精度评估标准3D-MCRB (3D Multi-hop Cramér-Rao Bound).仿真结果表明:此方法对信标位置误差和距离估计误差都具有较好的抑制能力,在大多数实验条件下,能将定位精度提高10%以上.     

北斗三频相位观测值线性组合研究

以北斗三频数据为基础,给出了北斗三频组合观测值的数学模型并系统分析了组合后的波长、电离层误差以及观测噪声误差;在保持整周模糊度整数特性的前提下,以整数线性变换法为基础,通过Matlab编程实现组合系数的优化选取。通过将不同组合系数应用于伪距相位法三频数据周跳探测中,发现组合系数之和不为0的组合周跳探测值均在0.5周以上,且组合系数之和越大,波动越大;通过将不同组合系数应用于无几何CIR法三频数据模糊度解算中,发现系数之和不为0,但波长较长的组合求得的模糊度残差值大致在0.6周,通过合适的模糊度搜索方法可以得到正确的模糊度固定解,在短基线模糊度解算中可以考虑使用。     

Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products

The state-space representation (SSR) product of satellite orbit and clock is one of the most essential corrections for real-time precise point positioning (RTPPP). When it comes to PPP ambiguity resolution (PPP-AR), the fractional cycle bias (FCB) matters. The Japan Aerospace Exploration Agency (JAXA) has developed a multi-GNSS (i.e., global navigation satellite system) advanced demonstration tool for orbit and clock analysis (MADOCA), providing free and precise orbit and clock products. Because of the shortage of relevant studies on performance evaluation, this paper focuses on the performance assessment of RTPPP and PPP-AR by real-time and offline MADOCA products. To begin with, the real-time MADOCA products are evaluated by comparing orbit and clock with JAXA final products, which gives an objective impression of the correction. Second, PPP tests in static and simulated kinematic mode are conducted to further verify the quality of real-time MADOCA products. Finally, the offline MADOCA products are assessed by PPP and PPP-AR comparisons. The results are as follows: (1) Orbit comparisons produced an average error of about 0.04–0.13 m for the global positioning system (GPS), 0.14–0.16 m for the global navigation satellite system (GLONASS), and 0.07–0.08 m for the quasi-zenith satellite system (QZSS). The G15 satellite had the most accurate orbit, with a difference of 0.04 m between the JAXA orbit products and MADOCA’s counterpart, while the R07 satellite had the least accurate orbit with a difference of 0.16 m. Clock products had an accuracy of 0.4–1.3 ns for GPS, 1.4–1.6 ns for GLONASS, and 0.7–0.8 ns for QZSS in general. The G15 satellite had the most accurate clock with a difference of only 0.40 ns between the JAXA clock products and MADOCA products, and the R07 satellite had the least accurate clock with a difference of 1.55 ns. The orbit and clock products for GLONASS performed worse than those of GPS and QZSS. (2) After convergence, the positioning accuracy was 3.0–8.1 cm for static PPP and 8.1–13.7 cm for kinematic PPP when using multi-GNSS observations and precise orbit and clock products. The PFRR station performed the good performance both in static and kinematic mode with an accuracy of 2.99 cm and 8.08 cm, respectively, whereas the CPNM station produced the worst static performance with an error of 8.09 cm, and the ANMG station produced the worst kinematic performance with a counterpart of 13.69 cm. (3) The PPP-AR solution was superior to the PPP solution, given that, with respect to PPP, post-processing PPP-AR improved the positioning accuracy and convergence time by 13–32 % (3–89 %) in GPS-only mode by 2–15 % (5–60 %) in GPS/QZSS mode. Thus, we conclude that the current MADOCA products can provide SSR corrections and FCB products with positioning accuracy at the decimeter or even centimeter level, which could meet the demands of the RTPPP and PPP-AR solutions.     

The GFZ real-time GNSS precise positioning service system and its adaption for COMPASS

Motivated by the IGS real-time Pilot Project, GFZ has been developing its own real-time precise positioning service for various applications. An operational system at GFZ is now broadcasting real-time orbits, clocks, global ionospheric model, uncalibrated phase delays and regional atmospheric corrections for standard PPP, PPP with ambiguity fixing, single-frequency PPP and regional augmented PPP. To avoid developing various algorithms for different applications, we proposed a uniform algorithm and implemented it into our real-time software. In the new processing scheme, we employed un-differenced raw observations with atmospheric delays as parameters, which are properly constrained by real-time derived global ionospheric model or regional atmospheric corrections and by the empirical characteristics of the atmospheric delay variation in time and space. The positioning performance in terms of convergence time and ambiguity fixing depends mainly on the quality of the received atmospheric information and the spatial and temporal constraints. The un-differenced raw observation model can not only integrate PPP and NRTK into a seamless positioning service, but also syncretize these two techniques into a unique model and algorithm. Furthermore, it is suitable for both dual-frequency and sing-frequency receivers. Based on the real-time data streams from IGS, EUREF and SAPOS reference networks, we can provide services of global precise point positioning (PPP) with 5–10 cm accuracy, PPP with ambiguity-fixing of 2–5 cm accuracy, PPP using single-frequency receiver with accuracy of better than 50 cm and PPP with regional augmentation for instantaneous ambiguity resolution of 1–3 cm accuracy. We adapted the system for current COMPASS to provide PPP service. COMPASS observations from a regional network of nine stations are used for precise orbit determination and clock estimation in simulated real-time mode, the orbit and clock products are applied for real-time precise point positioning. The simulated real-time PPP service confirms that real-time positioning services of accuracy at dm-level and even cm-level is achievable with COMPASS only.     

Performance analysis of PPP ambiguity resolution with UPD products estimated from different scales of reference station networks

Integer ambiguity fixing with uncalibrated phase delay (UPD) products can significantly shorten the initialization time and improve the accuracy of precise point positioning (PPP). Since the tracking arcs of satellites and the behavior of atmospheric biases can be very different for the reference networks with different scales, the qualities of corresponding UPD products may be also various. The purpose of this paper is to comparatively investigate the influence of different scales of reference station networks on UPD estimation and user ambiguity resolution. Three reference station networks with global, wide-area and local scales are used to compute the UPD products and analyze their impact on the PPP-AR. The time-to-first-fix, the unfix rate and the incorrect fix rate of PPP-AR are analyzed. Moreover, in order to further shorten the convergence time for obtaining precise positioning, a modified partial ambiguity resolution (PAR) and corresponding validation strategy are presented. In this PAR method, the ambiguity subset is determined by removing the ambiguity one by one in the order of ascending elevations. Besides, for static positioning mode, a coordinate validation strategy is employed to enhance the reliability of the fixed coordinate. The experiment results show that UPD products computed by smaller station network are more accurate and lead to a better coordinate solution; the PAR method used in this paper can shorten the convergence time and the coordinate validation strategy can improve the availability of high precision positioning.     

应用iGMAS超快速星历的实时精密单点定位研究

中国主导建设的国际GNSS监测评估系统（iGMAS）相比国际上比较成熟的IGS系统在产品精度等方面存在差别，目前实时精密单点定位应用多采用IGS实时、近实时产品。为改变这一现状，针对iGMAS产品特性以及实时精密单点定位对超快速精密星历的需求，对iGMAS超快速星历的精度和稳定性方面进行评估，设计了iGMAS产品实时/事后下载应用程序，开展了基于iGMAS超快速星历的实时精密单点定位研究，并结合NovAtel OEM617双频接收机进行了GPS实时精密单点定位试验。实验结果表明，在连续观测23min后定位误差即可收敛到分米级，较接收机原始定位精度高一个量级，且稳定性好，最终在E/N/U方向定位误差均方根分别为7.2cm、6.4cm、15.2cm，与应用IGS超快速星历实时PPP试验取得相近的结果。研究实现了iGMAS数据获取、评估和实时PPP应用的一整套方案，验证了iGMAS超快速产品的性能，对推进iGMAS产品的应用提供了借鉴。     

Towards PPP-RTK: Ambiguity resolution in real-time precise point positioning

Integer ambiguity resolution at a single station can be achieved by introducing predetermined uncalibrated phase delays (UPDs) into the float ambiguity estimates of precise point positioning (PPP). This integer resolution technique has the potential of leading to a PPP-RTK (real-time kinematic) model where PPP provides rapid convergence to a reliable centimeter-level positioning accuracy based on an RTK reference network. Nonetheless, implementing this model is technically subject to how rapidly we can fix wide-lane ambiguities, stabilize narrow-lane UPD estimates, and achieve the first ambiguity-fixed solution. To investigate these issues, we used 7 days of 1-Hz sampling GPS data at 91 stations across Europe. We find that at least 10 min of observations are required for most receiver types to reliably fix about 90% of wide-lane ambiguities corresponding to high elevations, and over 20 min to fix about 90% of those corresponding to low elevations. Moreover, several tens of minutes are usually required for a regional network before a narrow-lane UPD estimate stabilizes to an accuracy of far better than 0.1 cycles. Finally, for hourly data, ambiguity resolution can significantly improve the accuracy of epoch-wise position estimates from 13.7, 7.1 and 11.4 cm to 0.8, 0.9 and 2.5 cm for the East, North and Up components, respectively, but a few tens of minutes is required to achieve the first ambiguity-fixed solution. Therefore, from the timeliness aspect, our PPP-RTK model currently cannot satisfy the critical requirement of instantaneous precise positioning where ambiguity-fixed solutions have to be achieved within at most a few seconds. However, this model can still be potentially applied to some near-real-time remote sensing applications, such as the GPS meteorology.     

RTK model and positioning performance analysis using Galileo four-frequency observations

This paper proposes a real-time kinematic (RTK) model that uses one common reference satellite for the Galileo system with four frequency observations. In the proposed model, the double-differenced (DD) pseudorange and carrier phase biases among the different frequencies are estimated as unknown parameters to recover the integer features of the DD ambiguities among the different frequencies for ambiguity resolution and precise positioning. Analysis results show that the E5a, E5b, and E5 frequencies have virtually the same performance in terms of the positioning accuracy, observation residuals, and ratio values of ambiguity resolution. However, the E1 frequency performs worse than the E5a, E5b, and E5 frequencies. The RTK results for the combination of multiple frequencies are much better than those for a single-frequency observation, the coordinates’ standard deviation is improved about 20–30%, and the ambiguity fix time is improved about 10%.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

Performance improvement of real-time PPP ambiguity resolution using a regional integer clock

Obtaining reliable GNSS uncalibrated phase delay (UPD) or integer clock products is the key to achieving ambiguity-fixed solutions for real-time (RT) precise point positioning (PPP) users. However, due to the influence of RT orbit errors, the quality of UPD/integer clock products estimated with a globally distributed GNSS network is difficult to ensure, thereby affecting the ambiguity resolution (AR) performance of RT-PPP. In this contribution, by fully utilising the consistency of orbital errors in regional GNSS network coverage areas, a method is proposed for estimating regional integer clock products to compensate for RT orbit errors. Based on Centre National d’Études Spatiales (CNES) RT precise products, the regional GPS/BDS integer clock was estimated with a CORS network in the west of China. Results showed that the difference between the estimated regional and CNES global integer clock/bias products was generally less than 5 cm for GPS, whereas clock differences of greater than 10 cm were observed for BDS due to the large RT orbit error. Compared with PPP using global integer clock/bias products, the AR performance of PPP using the regional integer clock was obviously improved for four rover stations. For single GPS, the horizontal and vertical accuracies of ambiguity-fixed PPP solutions were improved by 56.2% and 45.3% on average, respectively, whereas improvements of 67.5% and 50.5% in the horizontal and vertical directions, respectively, were observed for the combined GPS/BDS situation. Based on a regional integer clock, the RMS error of a kinematic ambiguity-fixed PPP solution in the horizontal direction could reach 0.5 cm. In terms of initialisation time, the average time to first fix (TTFF) in combined GPS/BDS PPP was shortened from 18.2 min to 12.7 min. In view of the high AR performance realised with the use of regional integer clocks, this method can be applied to scenarios that require high positioning accuracy, such as deformation monitoring.