观测资料和卫星钟差对PPP参数收敛影响分析

在GPS 单点定位中, 参数解算的收敛时间和收敛稳定性是重要的研究内容之一, 影响收敛时间和收敛稳定性的因素很多, 本文主要就观测资料的不同采样间隔、卫星钟差资料的不同采样间隔、不同的定位精度要求对精密单点定位中参数收敛时间的影响进行了深入的分析探讨, 以中国上海GPS综合应用网中的12个测站资料为例, 分析了采样间隔、定位精度要求与收敛时间的关系, 并对不同采样间隔的收敛时间进行了统计分析, 得出一些初步结论.     

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

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

Near real-time estimation of tropospheric water vapour content from ground based GNSS data and its potential contribution to weather now-casting in Austria

The importance of high resolution meteorological analysis of the atmosphere increased over the past years. A detailed analysis of the humidity field is an important precondition for a better monitoring of local and regional extreme precipitation events and for forecasts with improved spatial resolution. For this reason, the Austrian Meteorological Agency (ZAMG) is operating the spatial and temporal high resolution INCA system (Integrated Now-casting through Comprehensive Analysis) since begin of 2005. Errors in this analysis occur mainly in the areas of rapidly changing and hard to predict weather conditions or rugged topography with extreme differences in height such as the alpine area of Austria. The aim of this work is to provide GNSS based measurements of the tropospheric water vapour content with a temporal resolution of 1 h and a temporal delay of less than 1 h to assimilate these estimates into the INCA system. Additional requirement is an accuracy of better than 1 mm of the precipitable water (PW) estimates.     

Displacement monitoring performance of relative positioning and Precise Point Positioning (PPP) methods using simulation apparatus

Besides the classical geodetic methods, GPS (Global Positioning System) based positioning methods are widely used for monitoring crustal, structural, ground etc., deformations in recent years. Currently, two main GPS positioning methods are used: Relative and Precise Point Positioning (PPP) methods. It is crucial to know which amount of displacement can be detected with these two methods in order to inform their usability according to the types of deformation. Therefore, this study conducted to investigate horizontal and vertical displacement monitoring performance and capability of determining the direction of displacements of both methods using a developed displacement simulator apparatus. For this purpose, 20 simulated displacement tests were handled. Besides the 24?h data sets, 12?h, 8?h, 4?h and 2?h subsets were considered to examine the influence of short time spans. Each data sets were processed using GAMIT/GLOBK and GIPSY/OASIS scientific software for relative and PPP applications respectively and derived displacements were compared to the simulated (true) displacements. Then statistical significance test was applied. Results of the experiment show that using 24?h data sets, relative method can determine up to 6.0?mm horizontal displacement and 12.3?mm vertical displacement, while PPP method can detect 8.1?mm and 19.2?mm displacements in horizontal and vertical directions respectively. Minimum detected displacements are found to grow larger as time spans are shortened.     

The realization and convergence analysis of combined PPP based on raw observation

In order to speed up Precise Point Positioning (PPP)’s convergence, a combined PPP method with GPS and GLONASS which is based on using raw observations is proposed, and the positioning results and convergence time have been compared with that of single system. The ionospheric delays and receiver’s Differential Code Bias (DCB) corrections are estimated as unknown parameters in this method. The numerical results show that the combined PPP has not caused significant impacts on the final solutions, but it greatly improved Position Dilution of Precision (PDOP) and convergence speed and enhanced the reliability of the solution. Meanwhile, the convergence speed is greatly influenced by the receiver’s DCB, positioning results in horizontal which are better than 10 cm can be realized within 10 min. In addition, the ionosphere and DCB products can be provided with high precision.     

An improved approach to model ionospheric delays for single-frequency Precise Point Positioning

PPP with low-cost, single-frequency receivers has been receiving increasing interest in recent years because of its large amount of possible users. One crucial issue in single-frequency PPP is the mitigation of ionospheric delays which cannot be removed by combining observations on different frequencies. For this purpose, several approaches have been developed, such as, the approach using ionospheric model corrections with proper weight, the GRAPHIC (Group and Phase Ionosphere Calibration) approach, and the method to model ionospheric delays over a station with a low polynomial or stochastic process. From our investigation on the stochastic characteristics of the ionospheric delay over a station, it cannot be precisely represented by either a deterministic model in the form of a low-order polynomial or a stochastic process for each satellite, because of its strong irregular spatial and temporal variations. Therefore, a novel approach is developed accordingly in which the deterministic representation is further refined by a stochastic process for each satellite with an empirical model for its power density. Furthermore, ionospheric delay corrections from a constructed model using GNSS data are also included as pseudo-observations for a better solution. A large data set collected from about 200 IGS stations over one month in 2010 is processed with the new approach and several commonly adopted approaches for validation. The results show its significant improvements in terms of positioning accuracy and convergence time with a negligible extra processing time, which is also demonstrated by data collected with a low-cost, handheld, single-frequency receiver.     

基于经验模分解的PPP时间比对数据消噪方法

精密单点定位(PPP)时间比对数据会受到观测噪声的影响,因此对时间比对数据进行消噪是一项重要工作。提出一种基于经验模分解(EMD)的PPP时间比对数据消噪方法,并将该方法与Vondrak滤波方法的消噪效果进行对比。结果表明,两种方法的消噪效果相当,均能有效滤除PPP时间比对数据中的随机噪声,明显改善时间比对的稳定度。     

Evaluation of COMPASS ionospheric model in GNSS positioning

As important products of GNSS navigation message, ionospheric delay model parameters are broadcasted for single-frequency users to improve their positioning accuracy. GPS provides daily Klobuchar ionospheric model parameters based on geomagnetic reference frame, while the regional satellite navigation system of China’s COMPASS broadcasts an eight-parameter ionospheric model, COMPASS Ionospheric Model(CIM), which was generated by processing data from continuous monitoring stations, with updating the parameters every 2 h. To evaluate its performance, CIM predictions are compared to ionospheric delay measurements, along with GPS positioning accuracy comparisons. Real observed data analysis indicates that CIM provides higher correction precision in middle-latitude regions, but relatively lower correction precision for low-latitude regions where the ionosphere has much higher variability. CIM errors for some users show a common bias for in-coming COMPASS signals from different satellites, and hence ionospheric model errors are somehow translated into the receivers’ clock error estimation. In addition, the CIM from the China regional monitoring network are further evaluated for global ionospheric corrections. Results show that in the Northern Hemisphere areas including Asia, Europe and North America, the three-dimensional positioning accuracy using the CIM for ionospheric delay corrections is improved by 7.8%–35.3% when compared to GPS single-frequency positioning ionospheric delay corrections using the Klobuchar model. However, the positioning accuracy in the Southern Hemisphere is degraded due apparently to the lack of monitoring stations there.     

基于经验模分解的PPP时间比对数据消噪方法

精密单点定位(PPP)时间比对数据会受到观测噪声的影响,因此对时间比对数据进行消噪是一项重要工作。提出一种基于经验模分解(EMD)的PPP时间比对数据消噪方法,并将该方法与Vondrak滤波方法的消噪效果进行对比。结果表明,两种方法的消噪效果相当,均能有效滤除PPP时间比对数据中的随机噪声,明显改善时间比对的稳定度。     

Inertial aided BDS triple-frequency integer ambiguity rounding method

The integer ambiguity resolution (AR) of carrier phase is significant for Global Navigation Satellite System (GNSS) precise positioning. However, in kinematic case, single-epoch AR methods based on alone GNSS are usually not reliable due to the instable pseudorange accuracy. Moreover, the computation of classical AR method Least Squares Ambiguity Decorrelation Adjustment (LAMBDA) is large. Thus, the inertial measurement unit (IMU) is introduced, a new inertial-aided AR method that directly rounds the float ambiguity of BeiDou triple-frequency combined observations, which is characterized by long wavelength, low carrier-phase noise and ionospheric delay, is proposed. The mathematical model of the new method is derived first. Then the impacts of the carrier-phase noise, ionospheric delay and inertial navigation system (INS) position error on the AR success ratio of combined observation are analyzed through probabilistic approach. Based on above investigation, the combinations (0, ?1, 1), (1, 4, ?5) and (4, ?2, ?3) are selected to resolve the original ambiguity. A vehicular integrated navigation test is performed to demonstrate the proposed method. The results show that the average AR success ratios of the three selected combinations, whose float ambiguity errors are 0.041, 0.146, 0.279 cycle respectively, are above 97.25% without regard to low-elevation C05. With respect to positioning accuracy based on our AR method when compared with IE software, the east, north, up error RMS of position are 0.042, 0.024, 0.069 m, respectively. In terms of the AR recover after the BeiDou signals outage, as long as 62 s BeiDou signal complete outage, all the ambiguities of all satellites could be re-fixed immediately. Besides, during the 90 s signals partial outage, the AR is not influenced by the position error, since the float ambiguity errors are all below half-cycle. The research of this contribution demonstrates the effectiveness of the proposed new method, which indicates it is applicable to kinematic positioning, even in BDS degraded and denied environments.     

基于三频数据的北斗卫星导航系统DCB参数精度评估方法

差分码偏差（Differential Code Biases,DCB）参数作为导航电文中重要的一项,是影响用户PNT服务的主要误差源之一。北斗卫星导航系统（后文简称“北斗系统”）发射三个频点的导航信号,在导航电文中需要发播卫星的2个TGD（Timing Group Delay）参数。文章首先介绍了北斗系统卫星DCB参数最小二乘解算与形式误差评估;其次根据北斗系统三频特点,提出了不同频点组合计算垂直方向电离层电子总含量（VTEC）互差的DCB精度定量评估方法,并与IGS(International GNSS Service)提供的GPS卫星DCB精度进行比较;最后,详细分析了DCB参数精度对用户等效距离误差（UERE）计算和定位计算的影响,分别采用卫星出场标定DCB参数和经过解算DCB参数进行评估。实测数据分析结果表明,北斗系统卫星DCB参数解算形式误差与IGS解算GPS卫星DCB参数形式误差相当,但受卫星类型和解算测站的几何分布限制,北斗系统卫星DCB参数解算不确定度相比IGS略差,估计精度优于0.5ns,不同频率组合计算VTEC互差绝对值均值优于0.6TECU。相比采用卫星出场标定值,采用系统解算DCB参数后,双频用户三维位置误差改善13.80%～47.42%。     

Evaluation on real-time dynamic performance of BDS in PPP,RTK, and INS tightly aided modes

Since China’s BeiDou satellite navigation system (BDS) began to provide regional navigation service for Asia-Pacific region after 2012, more new generation BDS satellites have been launched to further expand BDS’s coverage to be global. In this contribution, precise positioning models based on BDS and the corresponding mathematical algorithms are presented in detail. Then, an evaluation on BDS’s real-time dynamic positioning and navigation performance is presented in Precise Point Positioning (PPP), Real-time Kinematic (RTK), Inertial Navigation System (INS) tightly aided PPP and RTK modes by processing a set of land-borne vehicle experiment data. Results indicate that BDS positioning Root Mean Square (RMS) in north, east, and vertical components are 2.0, 2.7, and 7.6?cm in RTK mode and 7.8, 14.7, and 24.8?cm in PPP mode, which are close to GPS positioning accuracy. Meanwhile, with the help of INS, about 38.8%, 67.5%, and 66.5% improvements can be obtained by using PPP/INS tight-integration mode. Such enhancements in RTK/INS tight-integration mode are 14.1%, 34.0%, and 41.9%. Moreover, the accuracy of velocimetry and attitude determination can be improved to be better than 1?cm/s and 0.1°, respectively. Besides, the continuity and reliability of BDS in both PPP and RTK modes can also be ameliorated significantly by INS during satellite signal missing periods.     

A new stochastic model considering satellite clock interpolation errors in precise point positioning

Precise clock products are typically interpolated based on the sampling interval of the observational data when they are used for in precise point positioning. However, due to the occurrence of white noise in atomic clocks, a residual component of such noise will inevitable reside within the observations when clock errors are interpolated, and such noise will affect the resolution of the positioning results. In this paper, which is based on a twenty-one-week analysis of the atomic clock noise characteristics of numerous satellites, a new stochastic observation model that considers satellite clock interpolation errors is proposed. First, the systematic error of each satellite in the IGR clock product was extracted using a wavelet de-noising method to obtain the empirical characteristics of atomic clock noise within each clock product. Then, based on those empirical characteristics, a stochastic observation model was structured that considered the satellite clock interpolation errors. Subsequently, the IGR and IGS clock products at different time intervals were used for experimental validation. A verification using 179 stations worldwide from the IGS showed that, compared with the conventional model, the convergence times using the stochastic model proposed in this study were respectively shortened by 4.8% and 4.0% when the IGR and IGS 300-s-interval clock products were used and by 19.1% and 19.4% when the 900-s-interval clock products were used. Furthermore, the disturbances during the initial phase of the calculation were also effectively improved.     

Comparison of convergence time and positioning accuracy among BDS,GPS and BDS/GPS precise point positioning with ambiguity resolution

Precise point positioning (PPP) usually takes about 30?min to obtain centimetre-level accuracy, which greatly limits its application. To address the drawbacks of convergence speed and positioning accuracy, we develop a PPP model with integrated GPS and BDS observations. Based on the method, stations with global coverage are selected to estimate the fractional cycle bias (FCB) of GPS and BDS. The short-term and long-term time series of wide-lane (WL) FCB, and the single day change of narrow-lane (NL) FCB are analysed. It is found that the range of GPS and BDS non-GEO (IGSO and MEO) WL FCB is stable at up to a 30-day-time frame. At times frame of up to 60?days, the stability is reduced a lot. Whether for short-term or long-term, the changes in the BDS GEO WL FCB are large. Moreover, BDS FCB sometimes undergoes a sudden jump. Besides, 17 and 10 stations were used respectively to investigate the convergence speed and positioning errors with six strategies: BDS ambiguity-float PPP (Bfloat), GPS ambiguity-float PPP (Gfloat), BDS/GPS ambiguity-float PPP (BGfloat), BDS ambiguity-fixed PPP (Bfix), GPS ambiguity-fixed (Gfix), and BDS/GPS ambiguity-fixed (BGfix). The average convergence speed of the ambiguity-fixed solution is greatly improved compared with the ambiguity-float solution. In terms of the average convergence time, the Bfloat is the longest and the BGfix is the shortest among these six strategies. Whether for ambiguity-float PPP or ambiguity-fixed PPP, the convergence reduction time in three directions for the combined system is the largest compared with the single BDS. The average RMS value of the Bfix in three directions (easting (E), northing (N), and up (U)) are 2.0?cm, 1.5?cm, and 5.9?cm respectively, while those of the Gfix are 0.8?cm, 0.5?cm, and 1.7?cm. Compared with single system, the BDS/GPS combined ambiguity-fixed system (BGfix) has the fastest convergence speed and the highest accuracy, with average RMS as 0.7?cm, 0.5?cm, and 1.9?cm for the E, N, U components, respectively.     

GNSS global real-time augmentation positioning: Real-time precise satellite clock estimation,prototype system construction and performance analysis

Lots of ambiguities in un-differenced (UD) model lead to lower calculation efficiency, which isn’t appropriate for the high-frequency real-time GNSS clock estimation, like 1 Hz. Mixed differenced model fusing UD pseudo-range and epoch-differenced (ED) phase observations has been introduced into real-time clock estimation. In this contribution, we extend the mixed differenced model for realizing multi-GNSS real-time clock high-frequency updating and a rigorous comparison and analysis on same conditions are performed to achieve the best real-time clock estimation performance taking the efficiency, accuracy, consistency and reliability into consideration. Based on the multi-GNSS real-time data streams provided by multi-GNSS Experiment (MGEX) and Wuhan University, GPS + BeiDou + Galileo global real-time augmentation positioning prototype system is designed and constructed, including real-time precise orbit determination, real-time precise clock estimation, real-time Precise Point Positioning (RT-PPP) and real-time Standard Point Positioning (RT-SPP). The statistical analysis of the 6 h-predicted real-time orbits shows that the root mean square (RMS) in radial direction is about 1–5 cm for GPS, Beidou MEO and Galileo satellites and about 10 cm for Beidou GEO and IGSO satellites. Using the mixed differenced estimation model, the prototype system can realize high-efficient real-time satellite absolute clock estimation with no constant clock-bias and can be used for high-frequency augmentation message updating (such as 1 Hz). The real-time augmentation message signal-in-space ranging error (SISRE), a comprehensive accuracy of orbit and clock and effecting the users’ actual positioning performance, is introduced to evaluate and analyze the performance of GPS + BeiDou + Galileo global real-time augmentation positioning system. The statistical analysis of real-time augmentation message SISRE is about 4–7 cm for GPS, whlile 10 cm for Beidou IGSO/MEO, Galileo and about 30 cm for BeiDou GEO satellites. The real-time positioning results prove that the GPS + BeiDou + Galileo RT-PPP comparing to GPS-only can effectively accelerate convergence time by about 60%, improve the positioning accuracy by about 30% and obtain averaged RMS 4 cm in horizontal and 6 cm in vertical; additionally RT-SPP accuracy in the prototype system can realize positioning accuracy with about averaged RMS 1 m in horizontal and 1.5–2 m in vertical, which are improved by 60% and 70% to SPP based on broadcast ephemeris, respectively.