BM1308-52型GPS/BDS时间传递接收机研制与性能分析CSCD

随着北斗卫星导航系统(BDS)的发展与完善,基于BDS的时间传递应用需求越来越迫切。简要介绍了为开展北斗时间传递研究自研的多通道多频GPS/BDS时间传递接收机BM1308-52。接收机可同时接收GPS、BDS的码信息和载波相位信息,输出GPS、BDS的CGGTTS标准共视文件和Rinex观测文件,观测时间、处理方法及数据输出格式符合国际规范。最后,利用实测数据测试了BM1308-52的性能,测试结果表明,GPS单向时间比对和零基线共视比对不确定度优于2ns,BD单向时间比对不确定度优于3ns,与国际水平相当。BM1308-52的系统稳定可靠,观测精度高,可以更好地为时间频率传递服务。     

BM1308-52型GPS/BDS时间传递接收机研制与性能分析

随着北斗卫星导航系统(BDS)的发展与完善,基于BDS的时间传递应用需求越来越迫切。简要介绍了为开展北斗时间传递研究自研的多通道多频GPS/BDS时间传递接收机BM1308-52。接收机可同时接收GPS、BDS的码信息和载波相位信息,输出GPS、BDS的CGGTTS标准共视文件和Rinex观测文件,观测时间、处理方法及数据输出格式符合国际规范。最后,利用实测数据测试了BM1308-52的性能,测试结果表明,GPS单向时间比对和零基线共视比对不确定度优于2ns,BD单向时间比对不确定度优于3ns,与国际水平相当。BM1308-52的系统稳定可靠,观测精度高,可以更好地为时间频率传递服务。     

多通道GPS共视法时频传递接收机的研制

GPS共视法是国际上流行的远距离时间频率传递技术,核心是共视法接收机。我们成功研制了多通道GPS共视法时频传递接收机系统,硬件部分主要由自主研制的高精度时间间隔计数器和Motorola生产的VPONCORE GPS引擎组成,软件符合时间频率咨询委员会(CCTF)发布的GPS共视法数据处理软件标准化指南的要求,与单通道GPS定时接收机相比,界面更友好,操作更方便,具有很强的分析处理数据功能。经测试证明多通道GPS接收机零基线共钟共视时间比对的不确定度小于4 ns(仰角40°),与国外报道基本相同。     

GPS和类GPS测距技术在双星编队星座状态确定中的联合应用

为提高状态确定的精度,借鉴AFF技术引入一种辅助手段——类似GPS的伪距和载波相位测距技术(简称类GPS测距技术)．它的伪码码元和载波相位波长比GPS的更短,因而可获得更高精度的相对测量信息．针对对地观测的双星编队星座的状态确定任务,建立了联合GPS和类GPS测距技术进行编队星座状态整体确定的数学模型,并快速同时解算出GPS星间单差模糊度和类GPS星内单差模糊度,最后进行了数学仿真．仿真结果表明,编队状态的精度有明显提高,其中相对位置精度为10~(-3)m,姿态角精度为10~(-4)rad．仿真证明该方法有效．     

适用于单频GPS用户的周跳探测方法

针对单频GPS(Global Position System)用户的载波相位周跳探测问题,在传统的载波相位和伪距组合方法的基础上提出了一种新的周跳探测方法.该方法将载波相位和码伪距的组合观测量在历元间做差形成差分序列,以此构造实时的周跳探测量,由于该探测量易受观测噪声的影响,因此使用仰角指数模型作为经验模型对于测距噪声进行估计,基于假设检验的方法,通过对该探测量是否超过相应的检测门限来判断是否存在周跳.在实测的GPS观测数据基础上,对这一方法进行了验证,结果表明:相对于传统的周跳探测方法,该方法可以更加及时、准确地发现较小的周跳,对于单频的GPS用户具有较好的适应性.      

高动态环境扩频系统伪码延时的精确估计方法

针对在高动态环境直接扩频系统伪码延时测量所遇到的问题,分析了多普勒频移对伪码延时的影响.利用扩展卡尔曼算法(EKF)对高动态直接扩频信号的载波相位和频率进行了估计,并利用载波辅助技术测量载波相位的变化值来校正伪码延时环,减小多普勒频移对伪码延时的影响,得到了精确的延时估计值,提高了伪码延时锁定环的动态跟踪性能.     

基于相位条纹的高精度GPS码相位测量方法

GPS接收机在测量卫星到接收机的传播距离时,通常能得到码相位和载波相位2个基本测量值。虽然载波相位测量值比码相位测量值精度高,但存在整周模糊度的问题,在实际应用中比采用码相位的技术付出的代价高很多。因此,基于相位条纹技术,提出了一种高精度的码相位测量方法。在传统码跟踪环的基础上,通过提取互功率谱相位条纹的频率,得到高精度的码相位测量值,从而组装出高精度的码伪距。仿真实验结果表明:在信噪比为-15 dB的情况下,码相位测量误差均方差约0.37 m,优于传统延迟锁定环在相同条件下约1.82 m的跟踪精度。得到了比传统码跟踪环更高的码相位测量精度的同时,不需要解算载波相位的整周模糊度,对提高GPS定位精度具有研究意义和应用价值。      

GPS载波相位时间传递软件的实现

研究了基于精密单点定位原理的GPS载波相位时间传递方法的数学模型和数据处理方法,设计了相应的程序流程,并基于Visual C++和Fortran混合编程的方式编写了相应的计算程序GPSCPTT V1.0。采用实际观测数据进行实验验证,结果表明,软件可以实现高精度的时间传递。     

基于载波相位的多通道相位一致性测量方法CSCD

对于卫星导航系统,由导航天线进行波束合成的多路并行导航信号需要具有严格的相位一致性,本文给出基于专用地检设备的多通道相位一致性测量方法。该地检设备采用载波相位辅助伪码测距技术,通过载波相位和伪码相位的联合解算,实现了伪码相位对载波相位解整周期模糊,最终测距精度达到1ps。采用该测试方法对系统进行配相,实现了各通道无周期模糊的相位一致性。     

基于信道复用方法的GNSS共视信号模拟

全球导航卫星系统（GNSS）共视（CV）技术应用中需要对GNSS共视信号进行模拟仿真,可以降低对共视接收机和共视算法进行测试过程中的成本。为此,提出了一种基于信道复用方法的GNSS共视信号的双路信号模拟方法。首先,对GNSS共视技术原理进行了分析。然后,根据GNSS直射信号的模拟思路,设计了基于GNSS直射信号模拟器的GNSS共视信号模拟方法,对共视信号传播过程中可能产生的误差进行了分析。最后,对零基线、短基线、长基线3种场景下仿真的共视信号,以及实场采集的试验数据进行了验证分析。验证的结果表明,仿真的GNSS共视信号定位准确,定位精度在米级;共视比对结果均方根值（RMS）精度优于12 ns,可以进行共视法时间传递,证明了提出的共视信号模拟方法能够有效地用于GNSS共视信号生成。对GNSS共视信号模拟器、共视接收机的研制和共视算法的研究具有一定的理论参考意义和实际应用价值。      

Combining GPS and GLONASS for time and frequency transfer

Geodetic time and frequency transfer (TFT) consists in a comprehensive modeling of code and carrier phase observations from Global Navigation Satellite System (GNSS) in order to determine the synchronization errors between two remote clocks connected to GNSS receivers. Using either common view (CV), or Precise Point Positioning (PPP), current GNSS time transfer uses only GPS measurements. This study combines GPS and GLONASS observations in geodetic TFT in order to determine the added value of the GLONASS data in the results. Using the software Atomium, we demonstrate on one hand that using both constellations improves the solution for both CV and PPP results when analysing short data batches. In that case, there are not enough GPS code data to calibrate the solution, and additional GLONASS code data allows us to retrieve a correct absolute value for the solution. On the other hand, the CV results of frequency transfer are not significantly affected by adding GLONASS data, while in PPP the combination with GLONASS modifies the frequency transfer results, and in particular the daily frequency offset, with maximum differences of 150 ps between the TFT solutions obtained with GPS-only or GPS + GLONASS.     

Influence of ionospheric perturbations in GPS time and frequency transfer

The stability of GPS time and frequency transfer is limited by the fact that GPS signals travel through the ionosphere. In high precision geodetic time transfer (i.e. based on precise modeling of code and carrier phase GPS data), the so-called ionosphere-free combination of the code and carrier phase measurements made on the two frequencies is used to remove the first-order ionospheric effect. In this paper, we investigate the impact of residual second- and third-order ionospheric effects on geodetic time transfer solutions i.e. remote atomic clock comparisons based on GPS measurements, using the ATOMIUM software developed at the Royal Observatory of Belgium (ROB). The impact of third-order ionospheric effects was shown to be negligible, while for second-order effects, the tests performed on different time links and at different epochs show a small impact of the order of some picoseconds, on a quiet day, and up to more than 10 picoseconds in case of high ionospheric activity. The geomagnetic storm of the 30th October 2003 is used to illustrate how space weather products are relevant to understand perturbations in geodetic time and frequency transfer.     

GLONASS-based precise point positioning and performance analysis

Current precise point positioning (PPP) techniques are mainly based on GPS which has been extensively investigated. With the increase of available GLONASS satellites during its revitalization, GLONASS observations were increasingly integrated into GPS-based PPP. Now that GLONASS has reached its full constellation, there will be a wide interest in PPP systems based on only GLONASS since it provides a PPP implementation independent of GPS. An investigation of GLONASS-based PPP will also help the development of GPS and GLONASS combined PPP techniques for improved precision and reliability. This paper presents an observation model for GLONASS-based PPP in which the GLONASS hardware delay biases are addressed. In view of frequently changed frequency channel number (FCN) for GLONASS satellites, an algorithm has been developed to compute the FCN for GLONASS satellites using code and phase observations, which avoids the need to provide the GLONASS frequency channel information during data processing. The observation residuals from GLONASS-based PPP are analyzed and compared to those from GPS-based PPP. The performance of GLONASS-based PPP is assessed using data from 15 globally distributed stations.     

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.     

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

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

GNSS geodetic techniques for time and frequency transfer applications

We performed an initial analysis of the pseudorange data of the GIOVE-B satellite, one of the two experimental Galileo satellites currently in operation, for time transfer.¹ For this specific aim, software was developed to process the GIOVE-B raw pseudoranges and broadcast navigation messages collected by the Galileo Experimental Sensor Stations (GESS) tracking network, yielding station clock phase errors with respect to the Experimental Galileo System Time (EGST). The software also allows processing the Global Positioning System (GPS) P1 and P2 pseudorange data with broadcast navigation message collected at the same stations to obtain the station clock phase errors with respect to the GPS system time (GPST). Differencing these solutions between stations provides two independent means of GNSS time transfer. We compared these time transfer results with Precise Point Positioning (PPP) method applied to GPS data in combined carrier-phase and pseudorange mode as well as in pseudorange-only mode to show their relative merits. The PPP solutions in combined carrier-phase and pseudorange mode showed the least instability of the methods tested herein at all scales, at few parts in 10¹⁵ at 1 day for the stations processed, following a tau^−½ interval dependency. Conversely, the PPP solutions in pseudorange-only mode are an order of magnitude worst (few parts in 10¹⁴ at 1 day for the stations processed) following a tau⁻¹ power-law, but slightly better than the single-satellite raw GPS time transfer solutions obtained using the developed software, since the PPP least-squares solution effectively averages the pseudorange noise. The pseudorange noise levels estimated from PPP pseudorange residuals and from clock solution comparisons are largely consistent, providing a validation of our software operation. The raw GIOVE-B time transfer, as implemented in this work, proves to be slightly better than single-satellite raw GPS satellite time transfer, at least in the medium term. However, one of the processed stations shows a combined GPS P1 and P2 pseudorange noise level at 2 m, a factor 2 worst than usually seen for geodetic receivers, so the GPS time transfer results may not be at their best for the cases processed. Over the short term, the GPS single-satellite time transfer instability outperforms the GIOVE-B by an order of magnitude at 1 s interval, which would be due to the different characteristics of the tracking loop filters for GPS P1 and P2 on one hand and the GIOVE-B signals on the other. Even at this preliminary stage and using an experimental satellite system, results show that the GIOVE-B (and hence Galileo) signals offer interesting perspectives for high precision time transfer between metrological laboratories.     

Contribution analysis of QZSS to single-frequency PPP of GPS/BDS/GLONASS/Galileo

The Quasi-Zenith Satellite System (QZSS) established by the Japan Aerospace Exploration Agency mainly serves the Asia-Pacific region and its surrounding areas. Currently, four in-orbit satellites provide services. Most users of GNSS in the mass market use single-frequency (SF) receivers owing to the low cost. Therefore, it is meaningful to analyze and evaluate the contribution of the QZSS to SF precise point positioning (PPP) of GPS/BDS/GLONASS/Galileo systems with the emergence of GNSS and QZSS. This study compares the performances of three SF PPP models, namely the GRoup and PHase Ionospheric Correction (GRAPHIC) model, GRAPHIC with code observation model, and an ionosphere-constrained model, and evaluated the contribution of the QZSS to the SF PPP of GPS/BDS/GLONASS/Galileo systems. Moreover, the influence of code bias on the SF PPP of the BDS system is also analyzed. A two-week dataset (DOY 013–026, 2019) from 10 stations of the MGEX network is selected for validation, and the results show that: (1) For cut-off elevation angles of 15, 20, and 25°, the convergence times for the static SF PPP of GLONASS + QZSS are reduced by 4.3, 30.8, and 12.7%, respectively, and the positioning accuracy is similar compared with that of the GLONASS system. Compared with the BDS single system, the convergence times for the static SF PPP of BDS + QZSS under 15 and 25° are reduced by 37.6 and 39.2%, the horizontal positioning accuracies are improved by 18.6 and 14.1%, and the vertical components are improved by 13.9 and 21.4%, respectively. At cut-off elevation angles of 15, 20, and 25°, the positioning accuracy and precision of GPS/BDS/GLONASS/Galileo + QZSS is similar to that of GPS/BDS/GLONASS/Galileo. And the convergence times are reduced by 7.4 and 4.3% at cut-off elevation angles of 20 and 25°, respectively. In imitating dynamic PPP, the QZSS significantly improves the positioning accuracy of BDS and GLONASS. However, QZSS has little effect on the GPS-only, Galileo-only and GPS/BDS/GLONASS/Galileo. (2) The code bias of BDS IGSO and MEO cannot be ignored in SF PPP. In static SF PPP, taking the frequency band of B1I whose multipath combination is the largest among the frequency bands as an example, the vertical component has a systematic bias of approximately 0.4–1.0 m. After correcting the code bias, the positioning error in the vertical component is lower than 0.2 m, and the positioning accuracy in the horizontal component are improved accordingly. (3) The SF PPP model with ionosphere constraints has a better convergence speed, while the positioning accuracy of the three models is nearly equal. Therefore the GRAPHIC model can be used to get good positioning accuracy in the absence of external ionosphere products, but its convergence speed is slower.     

GNSS实时动态授时精度分析

基于GNSS载波相位观测值的实时动态授时,可有效规避PPP授时对实时精密轨道和钟差产品的依赖问题,对短距离动、静态高精度时间用户具有重要意义。为了更好地验证GNSS实时动态授时性能,基于中国科学院国家授时中心时间频率资源和三个GNSS跟踪站长达2个月的观测数据,以GPS系统为例开展了授时试验。与事后PPP时间传递相比,实时动态授时结果差异STD优于0.15ns;与光纤双向时间传递结果相比,实时动态授时结果差异STD优于0.5ns。试验表明,GNSS实时动态授时精度能够达到亚纳秒量级,可为下一步推广应用提供重要参考。     

Differential code bias estimation based on uncombined PPP with LEO onboard GPS observations

Differential Code Bias (DCB) is an essential correction that must be provided to the Global Navigation Satellite System (GNSS) users for precise position determination. With the continuous deployment of Low Earth Orbit (LEO) satellites, DCB estimation using observations from GNSS receivers onboard the LEO satellites is drawing increasing interests in order to meet the growing demands on high-quality DCB products from LEO-based applications, such as LEO-based GNSS signal augmentation and space weather research. Previous studies on LEO-based DCB estimation are usually using the geometry-free combination of GNSS observations, and it may suffer from significant leveling errors due to non-zero mean of multipath errors and short-term variations of receiver code and phase biases. In this study, we utilize the uncombined Precise Point Positioning (PPP) model for LEO DCB estimation. The models for uncombined PPP-based LEO DCB estimation are presented and GPS observations acquired from receivers onboard three identical Swarm satellites from February 1 to 28, 2019 are used for the validation. The results show that the average Root Mean Square errors (RMS) of the GPS satellite DCBs estimated with onboard data from each of the three Swarm satellites using the uncombined PPP model are less than 0.18 ns when compared to the GPS satellite DCBs obtained from IGS final daily Global Ionospheric Map (GIM) products. Meanwhile, the corresponding average RMS of GPS satellite DCBs estimated with the conventional geometry-free model are 0.290, 0.210, 0.281 ns, respectively, which are significantly larger than those obtained with the uncombined PPP model. It is also noted that the estimated GPS satellite DCBs by Swarm A and C satellites are highly correlated, likely attributed to their similar orbit type and space environment. On the other hand, the Swarm receiver DCBs estimated with uncombined PPP model, with Standard Deviation (STD) of 0.065, 0.037 and 0.071 ns, are more stable than those obtained from the official Swarm Level 2 products with corresponding STD values of 0.115, 0.101, and 0.109 ns, respectively. The above indicates that high-quality DCB products can be estimated based on uncombined PPP with LEO onboard observations.     

Performance assessment of real-time multi-GNSS integrated PPP with uncombined and ionospheric-free combined observables

For precise position services, the real-time precise point positioning (PPP) is a promising technology. The real-time PPP performance is expected to be improved by multi-system combination. The performance of real-time multi-system PPP needs to be periodically investigated, with the increasing number of available satellites and the continuously improved quality of real-time precise products of satellite clocks and orbits. In this study, a comprehensive performance assessment is conducted for the four-system integrated real-time PPP (FSIRT-PPP) with GPS, BDS, Galileo and GLONASS in both static and kinematic modes. The datasets from 118 stations spanning approximately a month are used for analysis, and the real-time stream CLK93 is employed. The superior performance of FSIRT-PPP is validated by comparing with the results of GPS/BDS, GPS/Galileo, GPS/GLONASS, GPS-only, BDS-only, Galileo-only and GLONASS-only cases. The FSIRT-PPP using ionospheric-free (IF) combined observables can achieve a convergence time of 10.9, 4.8 and 11.8 min and a positioning accuracy of 0.4, 0.5 and 0.7 cm in the static mode in the east, north and up directions, respectively, while the derived statistic is 15.4, 7.0 and 16.4 min, and 1.6, 1.2 and 3.4 cm in the kinematic mode in the three directions, respectively. Moreover, we also compare the position solutions of real-time PPP adopting IF combined and uncombined (UC) observables, and prove the mathematical equivalence between the two PPP models in the converged stage, provided that there are no external ionospheric corrections or constraints given to the estimated ionospheric delays in the UC model. The difference between the fully converged positioning accuracy of IF-based and UC-based real-time PPP is marginal, but the UC-based real-time PPP has longer convergence time due to the influence of the significant unmodeled time-varying errors in the real-time precise products as well as the different parameterization between them. For completeness, the real-time kinematic PPP results in harsh environments and the post-processed PPP results are also presented.