基于UKF的姿态无关GPS双差整周模糊度确定

基于运动的GPS双差整周模糊度确定方法需要预先知道载体的姿态,为了消除对姿态的依赖,研究了一种新的整周模糊度确定方法。首先将多天线载波双差相位矢量观测方程转化为标量观测方程,消去了对载体姿态矩阵的依赖。经过短时离线观测,以中心估计算法和协方差白化算法估计和修正双差整周模糊度的近似解及协方差矩阵。再以该近似解和协方差矩阵为初值,由无迹卡尔曼滤波(UKF)实时估计双差整周模糊度的精确解。最后对该方法进行了仿真,仿真结果表明经过离线批处理算法的预处理UKF可准确收敛于真实解,证明了方法的正确和有效性。     

脉冲星导航的整周模糊度解算方法研究

脉冲相位是X射线脉冲星导航方法的基本观测量,在测量的脉冲相位与预报相位比对的过程中存在整周期模糊度问题．借鉴GPS载波相位模糊度的解算方法,提出了单差搜索法、最小二乘搜索法和模糊度函数法3种空间搜索方法用以解算脉冲星导航系统的整周模糊度,并针对地球卫星完全迷失的情况进行了搜索仿真验证．结果表明,这3种方法都能够准确地解算出整周模糊度,均可用于脉冲星导航系统的在轨解算．其中最小二乘搜索法的性能稳定,搜索速度快;模糊度函数法搜索速度慢,但适用于装配单X射线探测器的脉冲星导航系统．     

基于DSP＋LAMBDA算法的GPS实时定位技术研究

主要分析了GPS载波相位整周模糊度LAMBDA求解算法,通过数据模拟测试来验证该算法在DSP上的工作状况。仿真计算证明,在DSP上实现LAMBDA算法可以满足GPS实时动态定位的要求。     

载波相位DGPS快速定位算法

介绍ＤＧＰＳ系统原理，分析讨论了载波相位ＤＧＰＳ定位中主要残差在中、低动态下的影响，提出了载波相位ＤＧＰＳ接收机中一种不需要解算整周模糊度而快速解算定位的一种算法。     

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

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

单频GPS动态相对定位的模糊度逼近/搜索解法

针对目前实时动态定位(RTK)中常用搜索解法可靠度不高的缺点,利用各历元诸双差方程所形成的法方程,对其进行解耦处理以便于实时计算.在此基础上综合了逼近、搜索两种解法各自的局部优势,构成了逼近/搜索的联合解法.实际测试表明,该算法在单频全球定位系统(GPS)载波相位差分的动态相对定位过程中,用2～3min的时间即可可靠地在航解算出整周模糊度.     

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

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

GPS/Galileo组合导航的模糊度搜索空间

为了组合导航的载波相位模糊度固定,将目前在GPS中常用的模糊度固定方法--最小二乘降相关平差（LAMBDA）法直接应用于GPS/Galileo组合模糊度固定,发现其搜索空间的确定方法并不能很好地适应GPS/Galileo组合中模糊度维数较高的情况。通过对常规LAMBDA搜索空间确定方法的分析比较,在传统方法的基础上提出了一种专门针对高维模糊度固定的搜索空间确定方法--修正法确定模糊度搜索空间。通过对修正法进行仿真试验,证明该方法能保证在GPS/Galileo组合定位模式下实际备选模糊度个数基本与预先设定的备选模糊度个数一致,进而能在不降低模糊度固定成功率的基础上有效提高LAMBDA模糊度固定的搜索效率,其性能优于传统的模糊度搜索空间确定方法。     

联合CDGPS技术和星间相对测量进行编队星座状态确定

以空间圆3星编队星座为对象,建立了联合GPS载波相位差分(Carrier phase Differential GPS,CDGPS)和星间相对测量进行编队星座状态确定的数学模型;利用高精度的星间相对测量信息给星间公里级基线提供厘米级约束,极大地缩小了星间单差模糊度的搜索空间,进而在卫星无需机动的情况下采用Bayes最小二乘法快速解算出星间GPS载波相位单差整周模糊度;最后数学仿真证明了方法的有效性,结果表明卫星间相对位詈确定精度达10^-2m．卫星姿态确定精度达10^-3rad．     

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

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

Using only observation station data for PPP ambiguity resolution by UPD estimation

This article proposes a new method for uncalibrated phase delay (UPD) estimation to improve the accuracy of precise point positioning (PPP), which uses only observation station data. This means that the station used to generate the UPDs is the same station to which they are applied. First, dual-frequency observation equations based on a raw PPP model are developed. Then, the UPDs are calculated from integer linear combinations of float ambiguities. Third, with the UPD corrections, the least-squares ambiguity decorrelation adjustment (LAMBDA) method is utilized to obtain the integer ambiguities. Since only observation station data are used for UPD estimation, the partial ambiguity resolution (PAR) method is adopted to increase the possibility of finding a subset of integer ambiguities. The UPD estimation and ambiguity resolution are performed in each epoch. To obtain the correct integer ambiguity, the ratio test and success rate (bootstrapping) are used to evaluate the estimated integer ambiguity. Finally, by treating the integer ambiguities as constants, fixed solutions can be obtained. Quality control is also applied throughout the entire data processing procedure to obtain high quality float and fixed solutions. Data from 22 stations of the International Global Navigation Satellite System (GNSS) Service (IGS) in East Asia on day of year (DOY) 206, 2017, are used to verify the feasibility of this method. The experimental results show that compared with the float solution, the proposed method can significantly improve the accuracy in the east, north and up directions by 24%, 21% and 18% for static PPP and 36%, 18% and 34% for dynamic PPP, respectively. However, the accuracy of the proposed method is still lower than that of the fixed solutions obtained by the PRIDE-PPPAR software, in which the fractional cycle bias is computed based on reference network data. These findings sufficiently show that the proposed method can offer better solution accuracy than the float solution. However, the quality of the UPDs estimated only from observation station data is not as good as that of the estimates obtained based on reference network data.     

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

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

Accuracy assessment of single and double difference models for the single epoch GPS compass

The single epoch GPS compass is an important field of study, since it is a valuable technique for the orientation estimation of vehicles and it can guarantee a total independence from carrier phase slips in practical applications. To achieve highly accurate angular estimates, the unknown integer ambiguities of the carrier phase observables need to be resolved. Past researches focus on the ambiguity resolution for single epoch; however, accuracy is another significant problem for many challenging applications. In this contribution, the accuracy is evaluated for the non-common clock scheme of the receivers and the common clock scheme of the receivers, respectively. We focus on three scenarios for either scheme: single difference model vs. double difference model, single frequency model vs. multiple frequency model and optimal linear combinations vs. traditional triple-frequency least squares. We deduce the short baseline precision for a number of different available models and analyze the difference in accuracy for those models. Compared with the single or double difference model of the non-common clock scheme, the single difference model of the common clock scheme can greatly reduce the vertical component error of baseline vector, which results in higher elevation accuracy. The least squares estimator can also reduce the error of fixed baseline vector with the aid of the multi-frequency observation, thereby improving the attitude accuracy. In essence, the “accuracy improvement” is attributed to the difference in accuracy for different models, not a real improvement for any specific model. If all noise levels of GPS triple frequency carrier phase are assumed the same in unit of cycles, it can be proved that the optimal linear combination approach is equivalent to the traditional triple-frequency least squares, no matter which scheme is utilized. Both simulations and actual experiments have been performed to verify the correctness of theoretical analysis.     

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.     

Predicting atmospheric delays for rapid ambiguity resolution in precise point positioning

Integer ambiguity resolution in precise point positioning (PPP) can shorten the initialization and re-initialization time, and ambiguity-fixed PPP solutions are also more reliable and accurate than ambiguity-float PPP solutions. However, signal interruptions are unavoidable in practical applications, particularly while operating in urban areas. Such signal interruptions can cause discontinuity of carrier phase arc, which introduces new integer ambiguities. Usually it will take approximately 15 min of continuous tracking to a reasonable number of satellites to fix new integer ambiguities. In many applications, it is impractical for a PPP user to wait for such a long time for the re-initialization. In this paper, a method for rapid ambiguity fixing in PPP is developed to avoid such a long re-initialization time. Firstly, the atmospheric delays were estimated epoch by epoch from ambiguity-fixed PPP solutions before the data gap or cycle slip occurs. A random walk procedure is then applied to predict the atmospheric delays accurately over a short time span. The predicted atmospheric delays then can be used to correct the observations which suffer from signal interruptions. Finally, the new ambiguities can be fixed with a distinct WL-LX-L3 (here LX denotes either of L1, L2) cascade ambiguity resolution strategy. Comprehensive experiments have demonstrated that the proposed method and strategy can fix zero-difference integer ambiguities successfully with only a single-epoch observation immediately after a short data gap. This technique works even when all satellites are interrupted at the same time. The duration of data gap bridged by this technique could be possibly extended if a more precise atmospheric delay prediction is found or on-the-fly (OTF) technology is applied. Based on the proposed method, real-time PPP with integer ambiguity fixing becomes more feasible in practice.     

Effect of biases in integrity monitoring for RTK positioning

In road transport, continuous high-accuracy positioning is required in real time. To ensure the proper functioning and safety of vehicular applications, integrity monitoring (IM) is needed to protect from the positioning errors under a certain alert limit (AL) with a pre-defined probability of misleading information (MI). In this study, a detailed threat model is developed for real-time kinematic (RTK) positioning application of short baselines. The model distinguishes between ambiguity-float and -fixed scenarios, and considers the influences of phase and code multipath as well as between-receiver atmospheric residuals. With the float ambiguities temporally constrained, the bias contribution that propagates with time-updated ambiguities was studied analytically for the horizontal protection level (HPL) in IM. Based on real data from both static and kinematic experiments, HPL was computed along the direction of the semi-major axis of the horizontal error ellipse. In ambiguity-float and -fixed cases, the HPL was mostly several meters and decimetres, respectively. It was found that time-propagated biases play a dominant role in the ambiguity-float HPL, and among them, phase and code multipath had in general the largest contributions. For ambiguity-fixed case, the phase multipath was found to play a dominant role in the HPL. This shows the importance of considering the biases in the RTK IM for both the ambiguity-float and -fixed scenarios. Given a horizontal alert limit (HAL) of 5 m, the availabilities of ambiguity-float solutions were low, i.e., below 50% for the static roof tests and below 5% for the kinematic road tests. For the ambiguity-fixed scenario, with HAL at 0.5 m, integrity availability was nearly 100% for the static roof tests and above 85% for the kinematic road tests.     

Testing a new multivariate GNSS carrier phase attitude determination method for remote sensing platforms

GNSS (Global Navigation Satellite Systems)-based attitude determination is an important field of study, since it is a valuable technique for the orientation estimation of remote sensing platforms. To achieve highly accurate angular estimates, the precise GNSS carrier phase observables must be employed. However, in order to take full advantage of the high precision, the unknown integer ambiguities of the carrier phase observables need to be resolved. This contribution presents a GNSS carrier phase-based attitude determination method that determines the integer ambiguities and attitude in an integral manner, thereby fully exploiting the known body geometry of the multi-antennae configuration. It is shown that this integral approach aids the ambiguity resolution process tremendously and strongly improves the capacity of fixing the correct set of integer ambiguities. In this contribution, the challenging scenario of single-epoch, single-frequency attitude determination is addressed. This guarantees a total independence from carrier phase slips and losses of lock, and it also does not require any a priori motion model for the platform. The method presented is a multivariate constrained version of the popular LAMBDA method and it is tested on data collected during an airborne remote sensing campaign.     

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

Integer phase ambiguity estimation in next-generation geodetic Very Long Baseline Interferometry

Next-generation Very Long Baseline Interferometry (VLBI) system designs are aiming at 1 mm global position accuracy. In order to achieve this, it is not only necessary to deploy improved VLBI systems, but also to develop analysis strategies that take full advantage of the observations taken. Since the new systems are expected to incorporate four independent radio frequency bands, it should be feasible to resolve phase ambiguities directly from post-correlation data, providing roughly an order of magnitude improvement in precision of the delay observable. As the unknown ambiguities are of integer nature, it is discussed here how they the can be resolved analytically using algorithms which have been developed for Global Navigation Satellite System (GNSS) applications. Furthermore, it will be shown that ionosphere contribution and source structure effects, so-called core-shifts, can be solved simultaneously with the delay, which is the main geodetic observable for follow-on analysis. In order to verify the proposed algorithm, simulated observations were created using parameters from actual design studies. It is shown that, even in the case of low signal-to-noise ratio observations, reliable phase ambiguity resolution can be achieved and it is discussed how the integer ambiguity recovery depends on the number of observations and signal-to-noise ratio.