基于GPS载波相位的微卫星姿态算法研究

GPS姿态系统是利用GPS载波相位测量来确定载体的航向和姿态角。本文对微卫星的GPS姿态系统进行了研究。重点解决短基线的GPS状态算法，讨论了姿态价格函数的了小化方法。测试结果表明，所提的算法对小于1m的基线是有效的。     

利用GPS载波相位实时测定动态飞行器姿态

本文阐述了利用GPS载波相位信号进行实时飞行器姿态测定的基本思想,建立了整个系统实时确定载体姿态的数学模型;给出了动态整周模糊单差的求解过程,阐明了其物理意义．     

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.     

使用三天线的联合定姿和欺骗检测方法

全球卫星导航系统（GNSS）的安全性已经引发了广泛关注。使用多天线的欺骗检测方法由于其独一无二的空间特性,已成为当前最有效的欺骗检测方法之一。提出一种使用三天线的联合定姿和欺骗检测方法,能够在确定天线载体姿态的同时检测出欺骗信号的存在。针对直接定姿法受限于基线向量精度的问题,使用长度约束的基线向量估计方法以获得高精度定姿结果。在姿态信息已知的前提下,根据星历信息、姿态变换矩阵及天线的几何关系可以获得载波相位单差的期望值。使用误差平方和（SSE）来评估载波相位单差观测值和期望值之间的偏差,并构建了欺骗信号的二元检验。结果表明：在无欺骗的情形下,所提方法能降低定姿的标准差76.1%以上;在有欺骗的情形,所提方法能够实现100%检测率,并降低定姿的标准差77.3%以上。     

新的GPS单频单历元定姿算法

LAMBDA算法依赖于初始的模糊度浮点解,但仅用载波相位观测方程需要多个历元才能获得浮点解,将导致初始化时间过长.针对这一问题,对GPS(Global Positioning System)单历元的载波相位单差方程进行特殊变换,将未知的整周模糊度看成噪声,从而构造出新的观测方程,和原始的观测方程进行组合求解,克服了仅用载波相位双差观测方程因为亏秩而无法在单历元获得浮点解的缺点,解决了初始化时间的问题.通过深入研究浮点解和固定解之间的关系,提出一种将低精度浮点解映射到固定解的方法,降低了LAMBDA算法对高精度浮点解的依赖性,避免了用多个历元获取浮点解的高精度,从而实现了单频、单历元的整周模糊度估计.通过实际测试,该算法成功率高于97%,能够有效地用于实时动态姿态解算.      

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.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

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.     

Carrier phase prediction method for GNSS precise positioning in challenging environment

Since the signals of global navigation satellite system (GNSS) are blocked frequently in challenging environments, the discontinuous carrier phases seriously affect the application of GNSS precise positioning. To improve the carrier phase continuity, this paper proposes a carrier phase prediction method based on carrier open-loop tracking. In the open-loop tracking mode, the carrier numerically controlled oscillator (NCO) is controlled by the predicted Doppler, but not by the loop filter output. To improve the phase prediction effective time, accurate receiver clock drift estimation is studied in the prediction method. The phase prediction performance is tested on GNSS software receiver. In the phase prediction effective time tests, open-loop processes were set for the tested channel. The test results show that, when some satellite signals are blocked in 15?s, the probability of carrier phase error less than quarter cycles is more than 94%. In the real time kinematic (RTK) positioning tests, some satellite signals are blocked in 10–15?s repeatedly. The test results show that, the carrier phase continuity is basically not affected by the signal interruption, and the RTK can almost keep continuous centimeter-level positioning accuracy without re-fixing the integer ambiguity.     

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/INS组合定姿算法

研究了一种利用GPS载波相位三差观测信息的多天线GPS/INS组合定姿算法,其中包含一个基于惯性测量信息的GPS载波相位周跳检测算法。最后,通过仿真分析验证了该组合算法可以有效提高定姿精度,同时具有较好的稳定性。     

Precise line-of-sight vector estimation based on an inter-satellite radio frequency system

This paper proposes a precise line-of-sight (LOS) vector estimation using an inter-satellite radio frequency system. GNSS-like technology is inherited such that the ranging signals are locally generated inside the formation. However, the approach differs from the standard GNSS model usage in that the LOS vector to be of a unit length is fully explored as a priori constraint for the carrier phase integer ambiguity resolution. The constraint is lumped to the mapping process from the real-valued ambiguities to the integers by what is called validation or subset ambiguity bounding. These two approaches have the same rules of regarding the constraint as a gateway to accept or reject the ambiguity candidates, but differ by using “all-ambiguity-set” and “subset-ambiguity”. Both show remarkable improvement with up to 80% lower integer fixing failure rates than without treating the constraint. Validation provides a slightly better performance than the subset ambiguity bounding in terms of the integer fixing failure rates and the computational efficiency. The predefined tolerance regions that are critical for these two methods are analytically determined as function of the carrier noise. The paper also introduces a LOS dependent ambiguity dilution of precision (ADOP_LOS) measure that can serve as a metric to characterize the expectation of being able to successfully resolve the ambiguities. The region of ADOP_LOS lower than 0.21 is empirically summarized as the safe region where the integer fixing failure rates are less than 1%. A closed form of the ADOP_LOS is derived which is able to capture the impact of the various factors. Antenna baseline geometries and multiple frequencies in the form of an ultra-BOC signal structure are demonstrated as the most important influencing factors. With multiple properly arrayed antennas and using ultra-BOC structure, instantaneous ambiguity resolution can be achieved and the LOS accuracy can reach millimeter level.     

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.     

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

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

Particle filter-based real-time phase line bias estimation for GNSS-based attitude determination with common-clock receivers

Global navigation satellite system (GNSS)-based attitude determination has been widely adopted in a wide variety of terrestrial, sea, air, and space applications. Recently, the emergence of commercial multi-GNSS common-clock receivers has brought new opportunities for high-precision GNSS-based attitude determination with single-differenced (SD) model. However, the key requirement of using this approach is the accurate estimation of the troublesome line bias (LB) in real-time. In this contribution, we propose a particle filter-based real-time phase LB estimation approach that apply to SD model with single-system single-frequency observations from common-clock receiver. We first analyzed the relationship between the integer ambiguity ratio value and the phase LB. It is proved that the accuracy of a given phase LB value can be qualified by the related ambiguity resolution ratio value, and the normalized ratio value can therefore be used to represent the likelihood function of observations. Then, we presented the particle filter-based real-time phase LB estimation procedure, and assessed its performance using GPS L1/BDS B1I observations from two datasets collected with different types of common-clock receivers in terms of the accuracy and convergence time of phase LB estimation, the computation load, and the positioning and attitude determination accuracy with respect to the double-differenced (DD) model. Experimental results demonstrated that the phase LB could be accurately estimated with short convergence time (generally within 15 epochs). Moreover, compared with the classical DD approach, the particle filter-based SD approach delivers comparable positioning root-mean-square (RMS) errors in the North and East components but significantly smaller RMS errors in the Up component. Accordingly, the achievable yaw accuracy is comparable whereas the pitch accuracy is remarkably improved. The improvements of positioning accuracy in the Up component and pitch accuracy are approximately 35.7 % to 63.7 %, and 33.3 % to 63.1 %, respectively. Additionally, the single-epoch computation time with our particle filter-based SD approach is generally 0.08 s, which is obviously larger than the DD approach but could still meet the requirements of real-time applications below 10 Hz sampling.     

基于GNSS测量的天宫二号质心确定

由于轨道机动燃料消耗,科学载荷加载、分离,以及伴飞小卫星在轨释放等原因引起天宫二号空间站质心(COM)发生位移,从而影响天宫二号的动力学质心定轨精度。针对这一问题,提出了基于全球导航卫星系统(GNSS)测量数据的简化动力学质心估计方法。燃料消耗是引起天宫二号质心发生位移的主要原因,质心在本体坐标系X轴方向位移最为显著。利用GNSS测量数据对天宫二号进行质心估计和精密定轨,在三轴对地稳定姿态下,本体坐标系X轴方向与轨道切向重合,定轨结果对本体坐标系X轴方向的质心位移并不敏感。但在连续偏航模式下,本体坐标系X轴在轨道法向上有较大分量,X轴方向的质心位移对基于GNSS测量计算的精密定轨结果有较大影响。定性和定量分析结果表明:偏航姿态模式下天宫二号本体坐标系X轴方向质心位移估计具有可行性。天宫二号实测数据计算结果表明:与未做质心估计的定轨结果进行对比,质心估计后表征轨道动力学建模误差的经验加速度补偿水平在轨道径向、切向和法向上分别降低62%、50%和65%;载波相位后验残差标准差降低0.04 cm;精密轨道与全球激光测距数据比较精度提高0.86 cm。所提方法可以应用于大型低轨航天器在轨质心估计。      

基于等效焦面的离轴遥感相机积分时间计算方法

目前在轨和在研的光学遥感卫星大都搭载的是离轴遥感相机,为了保证视轴正对星下点成像,相机需要整体俯仰一定角度补偿离轴角,这样就导致相机的焦平面无法平行于星下点水平面进行推扫成像,相机的积分时间计算不准确。针对以上问题,提出一种基于等效焦面的高精度积分时间计算方法,该方法构建与星下点水平面平行的等效焦面,通过建立严密的几何关系,求出真正电荷转移时间对应的像元尺寸,从而得到准确的积分时间。仿真试验表明,该方法可以将积分时间计算精度提高1.2%,为卫星在轨提高成像质量提供有效手段。     

Instantaneous BeiDou–GPS attitude determination: A performance analysis

The advent of modernized and new global navigation satellite systems (GNSS) has enhanced the availability of satellite based positioning, navigation, and timing (PNT) solutions. Specifically, it increases redundancy and yields operational back-up or independence in case of failure or unavailability of one system. Among existing GNSS, the Chinese BeiDou system (BDS) is being developed and will consist of geostationary (GEO) satellites, inclined geosynchronous orbit (IGSO) satellites, and medium-Earth-orbit (MEO) satellites. In this contribution, a BeiDou–GPS robustness analysis is carried out for instantaneous, unaided attitude determination.     

Virtual baseline method for Beidou attitude determination – An improved long-short baseline ambiguity resolution method

Ambiguity resolution (AR) is a critical step for successful attitude determination using carrier phase measurements of a satellite navigation system such as Beidou. This paper proposes an improved method for AR in support of Beidou attitude determination based on the concept of a “virtual baseline”. In the traditional long-short baseline method, the short baseline is limited to a length less than half of the carrier wave length of the Beidou signals. In the proposed method, a virtual short baseline is formed by differencing two collinear baselines. The AR equations for virtual short and long baselines are derived and the factors impacting the AR accuracy are analysed. Numerical simulation studies were carried out to evaluate the performance of the proposed AR method. The simulation results confirmed that the proposed method is an improvement over the traditional approach -- not only is it easier to deploy collinear antennas but also it keeps the capability of epoch-by-epoch AR, which makes it immune to cycle slips and there is no need for initialisation of ambiguity searching.     

GNSS遥感探测卫星星座设计

随着全球导航卫星系统(Global Navigation Satellite System,GNSS)掩星大气探测技术的兴起,GNSS遥感探测数据在气象数据资源中逐步占据重要地位,但是目前的掩星探测数量远不能满足数值天气预报等应用的需求,未来更需要充分利用GNSS信号资源,开展更大规模的GNSS掩星卫星星座探测.本文以...