基于单频星载GPS数据的低轨卫星精密定轨

为满足搭载单频GPS接收机低轨卫星的精密定轨需求以及深化单频定轨研究,文中解决了单频星载GPS数据的周跳探测问题,并利用“海洋二号”（HY-2A）卫星及“资源三号”（ZY-3）卫星的单频星载GPS实测数据采用两种方法确定了二者的简化动力学轨道,并通过观测值残差分析、与双频精密轨道比较、激光测卫数据检核等方法对所得轨道精度进行评定。结果表明,在不考虑电离层延迟影响的情况下,HY-2A卫星定轨精度为2~3dm,ZY-3卫星为1m左右;而采用半和改正组合消除电离层延迟一阶项影响后,二者定轨精度均显著提高,HY-2A卫星三维精度提高至1dm左右,ZY-3卫星提高至1~2dm。文章的研究成果表明,搭载单频GPS接收机的低轨卫星也可获得厘米级的定轨精度。     

星载GPS低轨卫星几何法精密定轨研究

本文讨论了星载ＧＰＳ接收机为单频情形的代轨卫星几何法定轨,包括载波一相对定轨法和动态网定轨法,并利用Ｔｏｐｅｘ／Ｐｏｓｅｉｄｏｎ卫星星载ＧＰＳ实测数据中Ｌ１载波相位观测值进行验证,结果表明,载波相位相对定轨精度与地面基准站的观测质量有关,其三轨道位置精度为分米级;载波相位动态风定轨精度介于各基准站皮相位相对定轨之间,它相当于在各基准站相对定轨之间加权均衡,而且提高了定轨的可靠性。     

附加先验轨道约束的LEO几何法实时精密定轨验证分析

低轨卫星的实时精密定轨能够极大拓展其完成复杂科学任务的能力,例如实时环境监测、机动控制和卫星自主导航等.本文根据几何法实时精密定轨模型,提出了附加LEO先验轨道约束从而改善实时定轨的精度、收敛速度和稳定性的构想.分别采用广播星历、超快速星历预报部分和实时精密星历,设计了6种实时定轨方案,并利用Swarm-A,B,C星7天的观测数据进行方案验证与分析.结果表明,使用广播星历、IGU和IGC星历的方案精度递增,附加先验轨道约束能够进一步提升精度.使用IGC星历并附加标准差为1m的先验轨道约束后,在径向、切向和法向的定轨精度分别达到6.12cm,5.55cm和4.98cm.此外,附加先验轨道约束能够显著提升收敛速度,使用IGC星历平均收敛时间约为31min,附加标准差为1m的先验轨道约束后收敛仅需约4min.      

卫星与测站钟差支持条件下的GEO卫星精密定轨

在基于伪距的GEO卫星精密定轨中, GEO卫星的静地特性导致定轨解算无法对星地组合钟差进行有效估计, 需要独立的时间同步支持. 本文讨论了卫星和测站钟差支持条件下的GEO卫星定轨原理, 利用仿真数据系统地分析了中国区域网跟踪条件下GEO卫星的定轨精度, 从定性和定量角度分析了钟差二次项、星地时间同步精度、站间时间同步精度及系统差等因素对定轨精度的影响.      

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

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

基于双频GPS观测的简化动力学最小二乘批处理精密定轨

星载双频GPS载波相位和伪距观测量已成为低轨卫星获取精确三维位置和速度信息的主要方式. 本文以非差消电离载波相位和伪距组合作为观测量,应用简化动力学最小二乘批处理方法进行地球低轨卫星的精密定轨,并给出完整定轨流程. 采用逐段常量的经验加速度对动力学模型误差进行补偿,描述了经验加速度敏感矩阵及稀疏带状矩阵求逆的有效计算方法. 利用GRACE-A卫星GPS观测数据对定轨位置精度进行分析,结果显示,三维位置定轨精度优于5cm,经验加速度在径向、切向和法向上的补偿水平不超过40nm·s-2,大气阻力系数和辐射光压系数的估计值符合物理实际,星载接收机钟差大致呈线性并具有短周期小波动.      

星载GPS几何法实时定轨有关问题的研究

首先讨论了星载GPS几何法实时定轨的绝对定位方法和各种差分技术。由于伪距差技术能克服GPS卫星的星历误差、卫星钟误差，特别是SA误差的影响，而且实现难度不大，所以应用它来实时定轨。实测数据的处理表明，它能明显提高定轨的精度。然后分析了星载GPS所受扰动影响的情况，对应用抗差估计削弱GPS卫星信号扰动的影响进行了试验，试验的结果说明抗差估计能进一步提高星载GPS几何法定轨的精度。     

基于自适应滤波的编队卫星实时相对定轨

提出基于自适应滤波的编队卫星实时相对定轨算法,利用2005-12-09—10两颗GRACE（Gravity Recovery and Climate Experiment）卫星的GPS（Global Positioning System）实测数据进行实时相对定轨试验计算,采用JPL（Jet Propulsion Laboratory）轨道对试验结果外部检核,结果表明:①自适应滤波相对定轨通过自适应因子,可以较好地平衡编队卫星的观测信息和相对动力学信息,其相对定轨结果精度优于Kalman滤波相对定轨结果;②自适应滤波相对定轨结果随着星间基线缩短而精度提高;③两颗GRACE卫星采用单频伪距和广播星历进行自适应滤波相对定轨,可以得到精度优于6cm的星间基线。     

基于光学测角数据的风云四号同步轨道卫星精密定轨

针对风云四号同步卫星的精密定轨和精度评估需求,首先利用地面光学测角数据对FY-4A卫星进行精密定轨,定轨后方位角和高度角的残差rms分别为0.25＂和0.45＂。与基于测距数据的轨道相比,位置精度在有测角数据的弧段内小于50m。进一步联合测角数据和测距数据对FY-4A卫星进行联合定轨,定轨后轨道重叠精度优于15m。利用联合定轨结果评估了基于测距数据的实时轨道产品精度,可以明显发现轨道精度随着测距数据的积累而逐步提高。     

载波相位/捷联惯导组合对高轨飞行器的导航

高轨飞行器可用卫星数目较少,信号空间链路损耗大,使用伪距进行测量的精度较低.提出基于GNSS(Global Navigation Satellite System)卫星载波相位与捷联惯导紧组合方法对高轨机动飞行器进行自主导航.该方法将连续跟踪的卫星初次可用时的整周模糊度的浮点解作为状态变量,通过平方根UKF建立了组合导航非线性滤波模型,提出了基于整周浮点解交集的滤波器故障检测方法.研究表明,提出的组合导航方法充分利用了载波相位高测量精度和系统性误差缓变的特点,提高了系统的可靠性和精度.     

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.     

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 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.     

基于经验加速度的低轨卫星轨道预报新方法

研究将定轨过程中的经验加速度应用于地球低轨卫星轨道预报的新方法. 利用GPS伪距观测数据和简化动力学最小二乘批处理方法对地球低轨卫星定 轨, 其中卫星位置、速度及大气阻力系数和辐射光压系数可以直接用于轨道预报. 作为简化动力学最重要特征的经验加速度呈现准周期、余弦曲线特点, 可通过 傅里叶级数拟合建模. 确定性动力学模型与补偿大气阻力模型误差的切向经验 加速度级数拟合模型组成增强型动力学模型用于提高轨道预报精度. 应用 GRACE-A星载GPS伪距观测数据和IGS超快星历定轨并进行轨道预报, 结果表明 轨道预报初值位置精度达到0.2m, 速度精度达到1×10-4m·s-1, 预报3天位置精度优于60m, 比只利用确定性动力学模型进行预报精度平 均提高2.3倍. 先定轨后预报的模式可用在星上自主精确导航系统中.      

Orbit determination of BeiDou navigation satellite system maneuvered satellites based on instantaneous velocity pulses parameters

The BeiDou navigation satellite system (BDS) comprises geostationary earth orbit (GEO) satellites as well as inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO) satellites. Owing to their special orbital characteristics, GEO satellites require frequent orbital maneuvers to ensure that they operate in a specific orbital window. The availability of the entire system is affected during the maneuver period because service cannot be provided before the ephemeris is restored. In this study, based on the conventional dynamic orbit determination method for navigation satellites, multiple sets of instantaneous velocity pulses parameters which belong to one of pseudo-stochastic parameters were used to simulate the orbital maneuver process in the orbital maneuver arc and establish the observed and predicted orbits of the maneuvered and non-maneuvered satellites of BeiDou regional navigation satellite system (BDS-2) and BeiDou global navigation satellite system (BDS-3). Finally, the single point positioning (SPP) technology was used to verify the accuracy of the observed and predicted orbits. The orbit determination accuracy of maneuvered satellites can be greatly improved by using the orbit determination method proposed in this paper. The overlapping orbit determination accuracy of maneuvered GEO satellites of BDS-2 and BDS-3 can improve 2–3 orders of magnitude. Among them, the radial orbit determination accuracy of each maneuvered satellite is basically better than 1 m. simultaneously, the combined orbit determination of the maneuvered and non-maneuvered satellites does not have a great impact on the orbit determination accuracy of the non-maneuvered satellites. Compared with the multi GNSS products (indicated by GBM) from the German Research Centre for Geosciences (GFZ), the impact of adding the maneuvered satellites on the orbit determination accuracy of BDS-2 satellites is less than 9 %. Furthermore, the orbital recovery time and the service availability period are significantly improved. When the node of the predicted orbit is traversed approximately 3 h after the maneuver, the accuracy of the predicted orbit of the maneuvered satellite can reach that of the observed orbit. The SPP results for the BDS reached a normal level when the node of the predicted orbit was 2 h after the maneuver.     

Analysis of GNSS clock prediction performance with different interrupt intervals and application to real-time kinematic precise point positioning

Continuous and timely real-time satellite orbit and clock products are mandatory for real-time precise point positioning (RT-PPP). Real-time high-precision satellite orbit and clock products should be predicted within a short time in case of communication delay or connection breakdown in practical applications. For prediction, historical data describing the characteristics of the real-time orbit and clock can be used as the basis for performing the prediction. When historical data are scarce, it is difficult for many existing models to perform precise predictions. In this paper, a linear regression model is used to predict clock products. Seven-day GeoForschungsZentrum (GFZ) final clock products sampled at 30 s are used to analyze the characteristics of GNSS clocks. It is shown that the linear regression model can be used as the prediction model for the satellite clock products. In addition, the accuracy of the clock prediction for different satellites are analyzed using historical data with different periods (such as 2 and 10 epochs). Experimental results show that the accuracy of the clock with the linear regression prediction model using historical data with 10 epochs is 1.0 ns within 900 s. This is higher accuracy than that achieved using historical data of 2 epochs. Finally, the performance analysis for real-time kinematic precise point positioning (PPP) is provided using GFZ final clock prediction results and state space representation (SSR) clock prediction results when communication delay or connection breakdown occur. Experimental results show that the positioning accuracy without prediction is better than that with prediction in general, whether using the final clock product or the SSR clock product. For the final clock product, the positioning accuracy in the north (N), east (E), and up (U) directions is better than 10.0 cm with all visible GNSS satellites with prediction. In comparison, the 3D positioning accuracy of N, E, and U directions with visible GNSS satellites whose prediction accuracy is better than 0.1 ns using historical data of 10 epochs is improved from 15.0 cm to 7.0 cm. For the SSR clock product, the positioning accuracy of N, E, and U directions is better than 12.0 cm with visible GNSS satellites with prediction. In comparison, the 3D positioning accuracy of N, E, and U directions with visible GNSS satellites whose prediction accuracy is better than 0.1 ns using historical data of 10 epochs is improved from 12.0 cm to 9.0 cm.     

The impact of orbital errors on the estimation of satellite clock errors and PPP

Precise satellite orbit and clocks are essential for providing high accuracy real-time PPP (Precise Point Positioning) service. However, by treating the predicted orbits as fixed, the orbital errors may be partially assimilated by the estimated satellite clock and hence impact the positioning solutions. This paper presents the impact analysis of errors in radial and tangential orbital components on the estimation of satellite clocks and PPP through theoretical study and experimental evaluation. The relationship between the compensation of the orbital errors by the satellite clocks and the satellite-station geometry is discussed in details. Based on the satellite clocks estimated with regional station networks of different sizes (∼100, ∼300, ∼500 and ∼700 km in radius), results indicated that the orbital errors compensated by the satellite clock estimates reduce as the size of the network increases. An interesting regional PPP mode based on the broadcast ephemeris and the corresponding estimated satellite clocks is proposed and evaluated through the numerical study. The impact of orbital errors in the broadcast ephemeris has shown to be negligible for PPP users in a regional network of a radius of ∼300 km, with positioning RMS of about 1.4, 1.4 and 3.7 cm for east, north and up component in the post-mission kinematic mode, comparable with 1.3, 1.3 and 3.6 cm using the precise orbits and the corresponding estimated clocks. Compared with the DGPS and RTK positioning, only the estimated satellite clocks are needed to be disseminated to PPP users for this approach. It can significantly alleviate the communication burdens and therefore can be beneficial to the real time applications.     

Comprehensive outage compensation of real-time orbit and clock corrections with broadcast ephemeris for ambiguity-fixed precise point positioning

The main challenge in real-time precise point positioning (PPP) is that the data outages or large time lags in receiving precise orbit and clock corrections greatly degrade the continuity and real-time performance of PPP positioning. To solve this problem, instead of directly predicting orbit and clock corrections in previous researches, this paper presents an alternative approach of generating combined corrections including orbit error, satellite clock and receiver-related error with broadcast ephemeris. Using ambiguities and satellite fractional-cycle biases (FCBs) of previous epoch and the short-term predicted tropospheric delay through linear extrapolation model (LEM), combined corrections at current epoch are retrieved and weighted with multiple reference stations, and further broadcast to user for continuous enhanced positioning during outages of orbit and clock corrections. To validate the proposed method, two reference station network with different inter-station distance from National Geodetic Survey (NGS) network are used for experiments with six different time lags (i.e., 5 s, 10 s, 15 s, 30 s, 45 s and 60 s), and one set of data collected by unmanned aerial vehicle (UAV) is also used. The performance of LEM is investigated, and the troposphere prediction accuracy of low elevation (e.g., 10–20degrees) satellites has been improved by 44.1% to 79.0%. The average accuracy of combined corrections before and after LEM is used is improved by 12.5% to 77.3%. Without LEM, an accuracy of 2–3 cm can be maintained only in case of small time lags, while the accuracies with LEM are all better than 2 cm in case of different time lags. The performance of simulated kinematic PPP at user end is assessed in terms of positioning accuracy and epoch fix rate. In case of different time lags, after LEM is used, the average accuracy in horizontal direction is better than 3 cm, and the accuracy in up direction is better than 5 cm. At the same time, the epoch fix rate has also increased to varying degrees. The results of the UAV data show that in real kinematic environment, the proposed method can still maintain a positioning accuracy of several centimeters in case of 20 s time lag.