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

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

BDS/GPS对流层延迟估计及对精密单点定位的影响

针对对流层湿延迟参数估计及水平梯度处理策略对精密单点定位影响的问题, 以9个多系统实验网实测数据为基础, 设计了6种对流层延迟处理方案并进行实验. 实验结果表明, 观测时长在8 h时段内, BDS/GPS定位模式下, 对流层延迟的水平梯度处理不宜采用随机游走模型; 采用随机游走估算对流层湿延迟参数, 且不估算或使用分段常数处理水平梯度的策略相对占优. 动静态精密单点定位结果显示, BDS/GPS采用分段常数或随机游走进行对流层湿延迟参数估算和分段常数处理水平梯度的策略相对占优, 采用随机游走处理水平梯度的策略稍弱. 总体而言, BDS较GPS估算的对流层延迟及定位精度稍弱; 海洋周边6种方案解算结果差异较内陆站大. 静态精密单点定位模式下, 随着观测时长的增加, 6种方案定位误差趋势趋于一致, 误差大小基本在2 cm内, 观测时长大于12 h时, 随机游走估算对流层湿延迟参数和不考虑水平梯度的处理策略整体较优.     

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

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

基于相位偏差的实时PPP模糊度固定

采用CNES发布的实时相位偏差数据,实现包含模糊度固定的实时精密单点定位.对全球10个IGS测站10天观测数据进行RTPPP解算,分别统计模糊度首次固定时间和定位精度,结果显示利用实时相位偏差数据能在平均30min内实现模糊度首次固定,模糊度固定时水平位置误差由6cm迅速降低至2cm左右,三维位置误差由10cm迅速降低至5cm左右,同时RTPPP模糊度固定在3h观测内可保持水平3cm、三维5cm左右的定位精度.通过分析得出,基于相位偏差的RTPPP模糊度固定技术具有较高的定位精度和定位稳定性,能够快速实现cm级定位.     

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

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

地磁暴期间北半球高纬度地区电离层变化特征及对精密定位的影响

基于加拿大地区高纬度电离层观测网的电离层闪烁观测数据,分析了2018年8月26日地磁暴事件引发的北半球高纬度地区电离层总电子含量(TEC)异常变化、TEC变化率指数(ROTI)及电离层相位闪烁的变化特征.结果表明:加拿大地区最大异常值约6 TECU,磁暴引发全球电离层TEC异常峰值高达20 TECU;加拿大地区电离层相...     

超精密微牛级微波离子推进技术及其在卫星超精度控制中的应用研究

空间技术的快速发展使得利用空间卫星的编队飞行构建大型空间星座成为可能,在引力波探测、射电望远镜编队、星座组网等任务方面具有重要作用。超精度控制是实现卫星高精度编队飞行的关键技术。推进系统是实现卫星编队长期高度稳定飞行的保证,从而实现内部科学装置的正确运行。不同于常规的推进系统,卫星精密编队超精度控制对推进系统的推力可调范围、分辨率、响应时间、推力的一致性等有着极高的要求。根据卫星精密编队任务需求,对微牛级推进系统的功能及技术要求进行了分析,提出了基于M2微波离子推力器的卫星超精度控制推进系统。阐述了M2超精密微牛级推进系统的关键技术和研究进展,为后续M2推力器在无拖曳控制方面的应用奠定了基础。     

基于实时解空间的三维声源定位

提出了一种基于实时解空间库的快速三维声源定位算法。该算法使用来自四个声传感器 ,一个风速传感器和一个风向传感器的信号能够准确快速地计算出声源的位置。     

起始点位置偏差对精密GPS差分定位影响的研究

从理论上讨论了起始点位置误差对ＧＰＳ基线解算精度的影响,并对实测资料进行了模拟计算,得到了与理论计算值相吻合的结果。指出：ＧＰＳ网中起始点的位置偏差不影响基线解的内部精度,但使ＧＰＳ基线解产生系统性偏差,偏差的大小与起始点位置偏差的大小成正比,与基线平均纬度的正切函数值成正比。文章最后,对在“ＳＡ”和“ＡＳ”技术影响下,进行高精度ＧＰＳ差分定位时,如何测定网中具有较高精度的ＷＧＳ８４坐标起始点问题,提出了建议。     

Network based real-time precise point positioning

Network based real-time precise point positioning system includes two stages, i.e. real-time estimation of satellite clocks based on a reference network and real-time precise point positioning thereafter. In this paper, a satellite- and epoch-differenced approach, adopted from what is introduced by Han et al. (2001), is presented for the determination of satellite clocks and for the precise point positioning. One important refinement of our approach is the implementation of the robust clock estimation. A prototype software system is developed, and data from the European Reference Frame Permanent Network on September 19, 2009 is used to evaluate the approach. Results show that our approach is 3 times and 90 times faster than the epoch-difference approach and the zero-difference approach, respectively, which demonstrates a significant improvement in the computation efficiency. The RMS of the estimated clocks is at the level of 0.1 ns (3 cm) compared to the IGS final clocks. The clocks estimates are then applied to the precise point positioning in both kinematic and static mode. In static mode, the 2-h estimated coordinates have a mean accuracy of 3.08, 5.79, 6.32 cm in the North, East and Up directions. In kinematic mode, the mean kinematic coordinates accuracy is of 4.63, 5.82, 9.20 cm.     

北斗卫星导航系统实时定轨与钟差处理策略

目前鲜有对北斗卫星导航系统（BeiDouNavigationSatelliteSystem,BDS）实时精密定轨与钟差确定的研究,文章提出了BDS实时轨道与实时钟差处理策略,包括了观测与动力学模型、实时轨道与实时钟差处理流程与评估方法。尤其对于实时钟差,为了提高计算效率,联合使用两个独立并行的线程估计非差绝对钟差和历元间相对钟差。利用多模全球卫星导航系统试验(MGEX)与全球连续检测评估系统（iGMAS）实测数据进行了北斗实时轨道与钟差解算,BDS实时轨道径向平均精度对于GEO卫星优于20cm,对于IGSO与MEO一般优于10cm;钟差精度对于GEO卫星为0.5~4.5ns,对于IGSO／MEO为0.2~2.0ns。基于目前的轨道与钟差结果,实时精密单点定位（PrecisePointPositioning,PPP）结果可以达到分米量级。     

The GLONASS-M satellite yaw-attitude model

The proper modeling of the satellites’ yaw-attitude is a prerequisite for high-precision Global Navigation Satellite System (GNSS) positioning and poses a particular challenge during periods when the satellite orbital planes are partially eclipsed. Whereas a lot of effort has been put in to examine the yaw-attitude control of GPS satellites that are in eclipsing orbits, hardly anything is known about the yaw-attitude behavior of eclipsing GLONASS-M satellites. However, systematic variations of the carrier phase observation residuals in the vicinity of the orbit’s noon and midnight points of up to ±27 cm indicate significant attitude-related modeling issues. In order to explore the GLONASS-M attitude laws during eclipse seasons, we have studied the evolution of the horizontal satellite antenna offset estimates during orbit noon and orbit midnight using a technique that we refer to as “reverse kinematic precise point positioning”. In this approach, we keep all relevant global geodetic parameters fixed and estimate the satellite clock and antenna phase center positions epoch-by-epoch using 30-second observation and clock data from a global multi-GNSS ground station network. The estimated horizontal antenna phase center offsets implicitly provide the spacecraft’s yaw-attitude. The insights gained from studying the yaw angle behavior have led to the development of the very first yaw-attitude model for eclipsing GLONASS-M satellites. The derived yaw-attitude model proves to be much better than the nominal yaw-attitude model commonly being used by today’s GLONASS-capable GNSS software packages as it reduces the observation residuals of eclipsing satellites down to the normal level of non-eclipsing satellites and thereby prevents a multitude of measurements from being incorrectly identified as outliers. It facilitates continuous satellite clock estimation during eclipse and improves in particular the results of kinematic precise point positioning of ground-based receivers.     

Satellite availability and positioning performance of uncombined precise point positioning using BeiDou-2 and BeiDou-3 multi-frequency signals

To realize the smooth transition from regional BeiDou Navigation Satellite System (BDS-2) to the global one (BDS-3), the integration of BDS-2 and BDS-3 is important for providing continuous, stable and reliable positioning, navigation and timing (PNT) services for global users. This work used 154 globally distributed multi-GNSS (Global Navigation Satellite System) experiment stations spanning 30 days to analyze the satellite availability and positioning performance of uncombined precise point positioning (UC-PPP) under current BDS-2 and BDS-3 constellations. We focused on three issues: the influence of BDS-3 receiver tracking abilities, the positioning performance among different areas, and the benefit of multi-frequency observations. The results show that the elliptical zone caused by poor BDS-2 satellite visibility is disappeared when the evenly distributed BDS-3 medium earth orbit satellites are introduced. When BDS-3 are integrated with BDS-2, the area with the Position Dilution of Precision (PDOP) less than 2 can be expanded to 75° S-75° N and 30° E-150° W. The positioning performance of BDS-3 and BDS-2/BDS-3 UC-PPP are seriously affected by the receiver tracking abilities of BDS-3 signals. When the maximum pseudo-random noise sequences (PRNs) of BDS-3 satellites tracked by stations are within 30 or 37, the positioning accuracy of static UC-PPP can be improved by 22.94% or 8.27% due to the integration of BDS-2 and BDS-3. Besides, the most improvement of BDS-2 and BDS-3 integration is achieved in Asia-Pacific regions, especially for the kinematic UC-PPP or the poor receiver tracking abilities of BDS-3. Similar to the multi-frequency BDS-2 UC-PPP, the benefit of multi-frequency signals for BDS-3 or BDS-2/BDS-3 UC-PPP is also non-vital. The three-dimensional positioning accuracy of BDS-2/BDS-3 multi-frequency UC-PPP in static mode and kinematic mode are 2.24 cm and 5.39 cm, while the corresponding convergence time are 49.62 min and 73.80 min, respectively. Compared with BDS-2, both the positioning accuracy and the convergence time of BDS-2/BDS-3 joint UC-PPP are improved by approximately over 50%, which indicates that BDS-3 has a great potential to provide high-quality PNT services as other global navigation satellite systems.     

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.     

Comparative analysis of four different single-frequency PPP models on positioning performance and atmosphere delay retrieval

Single-frequency precise point positioning (SF-PPP) has attracted increasing attention due to its high precision and cost effectiveness. With various strategies to handle the dominant error, i.e., ionosphere delay, the ionosphere-float (IF), ionosphere-free-half (IFH), ionosphere-corrected (IC), and ionosphere-weighted (IW) SF-PPP models are certain to possess different characteristics and performance levels. This study is dedicated to assessing and comparing the four models from model characteristics, positioning performance, and atmosphere delay retrieval. The model comparison shows that IC and IW models are full-rank while IF and IFH models have a rank deficiency of size one that will result in biased estimations, which means the better solvability of IC and IW models. The experiments are carried out based on the 7-day Global Positioning System (GPS) observations collected at 57 global Multi-GNSS Experiment (MGEX) stations and Global Ionosphere Map (GIM) products. The results indicate that the IW model can accelerate SF-PPP convergence and achieve higher positioning accuracy compared to the other three SF-PPP models, especially in kinematic mode. With convergence criteria of 0.25 m in horizontal and 0.5 m in vertical, the east/north/up convergence times of IW model are 0.5/15.0/25.0 min and 0.5/16.0/36.5 min for static and kinematic modes, respectively. The IW model is able to achieve an instantaneous positioning accuracy of 0.28/0.35/0.75 m. In addition, a real kinematic test also demonstrates the best positioning solutions of IW model. Regarding troposphere delay retrieval, the IF, IFH, and IW models obtain a comparable daily accuracy of 3.0 cm on average, while the IC model achieves the worst accuracy of 8.0 cm. For precise ionosphere delay estimation, IW model only needs an average initialization time of 34.3 min, but a longer initialization time of 51.6 min is required for IF model. The daily precision of ionosphere delay estimation for IW model can reach up to 10.8 cm. At the present accuracy of GIM products, it is suggested that the IW model should be adopted for SF-PPP first due to its superior performance in positioning and atmosphere delay retrieval.