Galileo/GPS载波相位组合观测值可匹配的研究

详细讨论了Galileo系统的4种载波与GPS L2载波的相位组合观测值的可匹配问题.在模糊度保持为整数的前提下,给出了Galileo/GPS载波相位组合观测值的定义,并对包括系统噪声和观测噪声在内的有关误差影响加以分析,最后给出了组合观测值可匹配的定义和判断可匹配的充要条件.      

精密进近阶段的多系统GNSS组合RAIM可用性算法及分析

给出了多系统全球卫星导航系统（GNSS）组合接收机自主完好性监测（ReceiverAutonomousIntegrityMonitoring,RAIM）可用性计算方法,在此基础上利用GPS、GLONASS实测数据与BDS、Galileo全星座仿真数据,分析了BDS、GPS、GLONASS和Galileo不同组合在精密进近阶段的RAIM可用性。通过试验分析发现,BDS的5颗地球同步轨道卫星和3颗倾斜地球同步轨道卫星对亚洲、非洲和欧洲大部分地区的RAIM可用性有很大的贡献。这些地区站星间几何观测结构得到改善,使得RAIM可用性相对于其他地区有很大幅度的提升。在亚太地区APV-I阶段单系统导航情况下,北斗导航系统RAIM可用性达到99.5%,高于其他三个导航系统。在精密进近阶段（APV-I、APV-II和CAT-I）,BDS与其他导航系统（GPS、GLONASS和Galileo）的组合导航可以满足全球大部分区域的RAIM可用性需求,大多可达到100%。     

SINS/GPS组合系统姿态角误差可观测性研究

首先建立起捷联惯导系统（ＳＩＮＳ）与全球定位系统（ＧＰＳ）相结合系统的数学模型,从研究ＳＩＮＳ／ＧＰＳ组合系统的可观测性出发,通过理论分析得出了组合系统姿态角误差可观测性很差的结论,并给出使可观测阵满秩的条件。对３种飞行条件所做仿真分析验证了理论分析的正确性。     

Galileo卫星导航系统广播Nequick模型及其参数拟合

提出了基于IGRF模型的Galileo广播Nequick模型及其参数拟合算法, 解决了Galileo信号仿真中地理场景映射与地磁坐标下的电离层延时修正参数拟合问题. 应用IGRF模型, 可计算出任意给定位置和时间点的地磁参数以及E层、 F₁层、F₂层的电子密度, 从而计算出Galileo电离层修正参数. 仿真结果表明, 该算法拟合的全球电离层延时与IGS提供的实际观测值基本一致, 仿真精度高于一般的经验电离层模型, 实现了Galileo卫星信号的电离层延时修正参数的精确仿真.      

GPS/SINS超紧组合导航系统可观测性分析

系统状态的可观测性决定了系统的滤波精度.目前,分析系统状态可观测性的主要方法是针对线性系统的,其中奇异值分解(SVD,Singular Value Decomposition)法可以对系统进行事前分析,并且计算简单,成为广泛使用的方法之一.然而集中滤波结构的全球定位系统/捷联惯性导航系统(GPS/SINS,Global Positioning System/Strapdown Inertial Navigation System)超紧组合模型的量测方程是非线性的,对非线性系统进行可观测性分析是非常困难的.针对这个问题,推导出了超紧组合模型线性的量测方程,并证明了该模型满足分段定常系统(PWCS,Piece-Wise Constant System)定理;使用奇异值分解法对其做可观测性仿真,并对实验结果做了全面分析.结果表明:系统状态的可观测性受载体运动状态影响,载体机动性越强,可观测的状态越多,可观测性越高.     

北斗三频相位观测值线性组合研究

以北斗三频数据为基础,给出了北斗三频组合观测值的数学模型并系统分析了组合后的波长、电离层误差以及观测噪声误差;在保持整周模糊度整数特性的前提下,以整数线性变换法为基础,通过Matlab编程实现组合系数的优化选取。通过将不同组合系数应用于伪距相位法三频数据周跳探测中,发现组合系数之和不为0的组合周跳探测值均在0.5周以上,且组合系数之和越大,波动越大;通过将不同组合系数应用于无几何CIR法三频数据模糊度解算中,发现系数之和不为0,但波长较长的组合求得的模糊度残差值大致在0.6周,通过合适的模糊度搜索方法可以得到正确的模糊度固定解,在短基线模糊度解算中可以考虑使用。     

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

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

软件GNSS信号模拟器架构与中频信号生成

软件GNSS(Global Navigation Satellite Systems)信号模拟器对于GNSS接收机的高效研发将做出重要贡献,因其结构灵活、开放性以及低成本.以GPS/Galileo组合系统为例,讨论了软件GNSS中频信号模拟器的架构,主要功能模块包括卫星星座仿真、接收机轨迹生成、传播通道特性仿真(包括电离层模型、对流层模型、多径模型等)、数字中频信号生成.在此基础上,着重阐述了数字中频信号生成模块的实现,功率谱图及分析结果验证了所生成的信号,包括GPS L1 C/A,Galileo E1 CBOC(Composite Binary Offset Carrier),Galileo E5a和E5b信号.     

卫星不足情况下低成本MEMS-INS/GPS伪松组合导航

为了解决GPS可观测卫星不足情况下低成本微电子机械-惯性导航系统/全球定位系统（MEMS-INS/GPS）组合导航精度维持问题,提出基于灰色模型和自适应卡尔曼滤波的MEMS-INS/GPS伪松组合导航方法。以MEMS-INS/GPS松组合导航模式为框架,建立了伪松组合导航系统的状态空间模型。基于MEMS-INS/GPS的历史观测数据,使用灰色模型对MEMSINS/GPS观测差值进行预测,称为系统伪观测量。当GPS可观测卫星充分时,使用噪声自适应估计卡尔曼滤波对MEMS-INS/GPS进行松组合导航;当GPS可观测卫星不足时,使用噪声自适应估计卡尔曼滤波依据系统伪观测量,将MEMS-INS/GPS进行伪松组合导航。以车载低成本MEMSINS/GPS组合导航系统为例进行仿真和实验验证,结果表明：当GPS可观测卫星不足时,传统的MEMS-INS/GPS松组合导航精度迅速下降并发散,而MEMS-INS/GPS伪松组合导航精度与GPS正常工作时的导航精度相差不大,维持了较高精度的导航状态。     

利用GNSS反射信号载波测量湖面高度变化

给出一种利用导航卫星反射信号的载波相位测量光滑湖面高度变化的方法.首先建立载波相位观测方程,通过在直射与反射信号的观测量之间求单差,再在历元间求双差的方法,消除了大量误差并得到了直射与反射信号的几何路径差;然后根据各颗卫星的观测值和高度角等信息计算出湖面高度变化,并把各卫星得到的湖面高度度变化值进行加权平均以提高测量系统的精度和可靠性.仿真验证了方法在-15dB信噪比下达到了毫米级的测量精度.     

The impact on tsunami detection from space using GNSS-reflectometry when combining GPS with GLONASS and Galileo

The devastating Sumatra tsunami in 2004 demonstrated the need for a tsunami early warning system in the Indian Ocean. Such a system has been installed within the German-Indonesian Tsunami Early Warning System (GITEWS) project. Tsunamis are a global phenomenon and for global observations satellites are predestined. Within the GITEWS project a feasibility study on a future tsunami detection system from space has therefore been carried out. The Global Navigation Satellite System Reflectometry (GNSS-R) is an innovative way of using GNSS signals for remote sensing. It uses ocean reflected GNSS signals for sea surface altimetry. With a dedicated Low Earth Orbit (LEO) constellation of satellites equipped with GNSS-R receivers, densely spaced sea surface height measurements could be established to detect tsunamis. Some general considerations on the geometry between LEO and GNSS are made in this simulation study. It exemplary analyzes the detection performance of a GNSS-R constellation at 900 km altitude and 60° inclination angle when applied to the Sumatra tsunami as it occurred in 2004. GPS is assumed as signal source and the combination with GLONASS and Galileo signals is investigated. It can be demonstrated, that the combination of GPS and Galileo is advantageous for constellations with few satellites while the combination with GLONASS is preferable for constellations with many satellites. If all three GNSS are combined, the best detection performance can be expected for all scenarios considered. In this case an 18 satellite constellation will detect the Sumatra tsunami within 17 min with certainty, while it takes 53 min if only GPS is considered.     

Multivariate bootstrapped relative positioning of spacecraft using GPS L1/Galileo E1 signals

GNSS-based precise relative positioning between spacecraft normally requires dual frequency observations, whereas attitude determination of the spacecraft, mainly due to the stronger model given by the a priori knowledge of the length and geometry of the baselines, can be performed precisely using only single frequency observations. When the Galileo signals will come available, the number of observations at the L1 frequency will increase as we will have a GPS and Galileo multi-constellation. Moreover the L1 observations of the Galileo system and modernized GPS are more precise than legacy GPS and this, combined with the increased number of observations, will result in a stronger model for single frequency relative positioning. In this contribution we will develop an even stronger model by combining the attitude determination problem with relative positioning. The attitude determination problem will be solved by the recently developed Multivariate Constrained (MC-) LAMBDA method. We will do this for each spacecraft and use the outcome for an ambiguity constrained solution on the baseline between the spacecraft. In this way the solution for the unconstrained baseline is bootstrapped from the MC-LAMBDA solutions of each spacecraft in what is called: multivariate bootstrapped relative positioning. The developed approach will be compared in simulations with relative positioning using a single antenna at each spacecraft (standard LAMBDA) and a vectorial bootstrapping approach. In the simulations we will analyze single epoch, single frequency success rates as the most challenging application. The difference in performance for the approaches for single epoch solutions, is a good indication of the strength of the underlying models. As the multivariate bootstrapping approach has a stronger model by applying information on the geometry of the constrained baselines, for applications with large observation noise and limited number of observations this will result in a better performance compared to the vectorial bootstrapping approach. Compared with standard LAMBDA, it can reach a 59% higher success rate for ambiguity resolution. The higher success rate on the unconstrained baseline between the platforms comes without extra computational load as the constrained baseline(s) problem has to be solved for attitude determination and this information can be applied for relative positioning.     

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

Cycle slip detection using multi-frequency GPS carrier phase observations: A simulation study

The detection and repair of the cycle slip or gross error is a key step for high precision global positioning system (GPS) carrier phase navigation and positioning due to interruption or unlocking of GPS signal. A number of methods have been developed to detect and repair cycle slips in the last two decades through cycle slip linear combinations of available GPS observations, but such approaches are subject to the changing GPS sampling and complex algorithms. Furthermore, the small cycle slip and gross error cannot be completely repaired or detected if the sampling is quite longer under some special observation conditions, such as Real Time Kinematic (RTK) positioning. With the development of the GPS modernization or Galileo system with three frequencies signals, it may be able to better detect and repair the cycle slip and gross error in the future. In this paper, the cycle slip and gross error of GPS carrier phase data are detected and repaired by using a new combination of the simulated multi-frequency GPS carrier phase data in different conditions. Results show that various real-time cycle slips are completely repaired with a gross error of up to 0.2 cycles.     

Calibration errors in determining slant Total Electron Content (TEC) from multi-GNSS data

The global navigation satellite system (GNSS) is presently a powerful tool for sensing the Earth's ionosphere. For this purpose, the ionospheric measurements (IMs), which are by definition slant total electron content biased by satellite and receiver differential code biases (DCBs), need to be first extracted from GNSS data and then used as inputs for further ionospheric representations such as tomography. By using the customary phase-to-code leveling procedure, this research comparatively evaluates the calibration errors on experimental IMs obtained from three GNSS, namely the US Global Positioning System (GPS), the Chinese BeiDou Navigation Satellite System (BDS), and the European Galileo. On the basis of ten days of dual-frequency, triple-GNSS observations collected from eight co-located ground receivers that independently form short-baselines and zero-baselines, the IMs are determined for each receiver for all tracked satellites and then for each satellite differenced for each baseline to evaluate their calibration errors. As first derived from the short-baseline analysis, the effects of calibration errors on IMs range, in total electron content units, from 1.58 to 2.16, 0.70 to 1.87, and 1.13 to 1.56 for GPS, Galileo, and BDS, respectively. Additionally, for short-baseline experiment, it is shown that the code multipath effect accounts for their main budget. Sidereal periodicity is found in single-differenced (SD) IMs for GPS and BDS geostationary satellites, and the correlation of SD IMs over two consecutive days achieves the maximum value when the time tag is around 4?min. Moreover, as byproducts of zero-baseline analysis, daily between-receiver DCBs for GPS are subject to more significant intra-day variations than those for BDS and Galileo.     

Coronal radio-sounding detection of a CME during the 1997 Galileo solar conjunction

Frequency fluctuations of the Galileo S-band radio signal were recorded nearly continuously during the spacecraft’s solar conjunction from December 1996 to February 1997. A strong propagating disturbance, most probably associated with a coronal mass ejection (CME), was detected on 7 February when the radio ray path proximate point was on the west solar limb at about 54 solar radii from the Sun. The CME passage through the line of sight is characterized by a significant increase in the fluctuation intensity of the recorded frequency and by an increase in the plasma speed from about 234 km s⁻¹ up to about 755 km s⁻¹. These velocity estimates are obtained from a correlation analysis of frequency fluctuations recorded simultaneously at two widely-separated ground stations. The density turbulence power spectrum is found to steepen behind the CME front. The Galileo radio-sounding data are compared with SOHO/LASCO observations of the CME in the corona and with WIND spacecraft data near the Earth’s orbit.     

The Planetary Laboratory for Image Analysis (PLIA)

The observation, characterization, and understanding of planetary atmospheres are key components of Solar System exploration. Their study requires, among other observations, the use of image data acquired from spacecrafts. The basic tasks that are generally used to work with planetary image data are: image navigation, projection, image processing, limb-darkening correction, composition of multiwavelength images and mosaicing. We illustrate the mathematical techniques required for the study of planetary images obtained from spacecrafts using as an example our own software package (PLIA). We present examples of analysis of Galileo SSI, Cassini ISS and Venus Express VIRTIS observations covering a large range of image peculiarities. The paper is organized as a guide for anyone willing to replicate one or several of the capabilities incorporated in the PLIA software. An open source version of PLIA for astronomical image processing is available.     

Estimating inter-system biases for tightly combined Galileo/BDS/GPS RTK

The overlapping carrier frequencies L1/E1, L5/E5a and B2/E5b from GPS/Galileo/BDS allow inter-system double-differencing of observations, which shows a clear advantage over differencing of the observations of each constellation independently. However, the inter-system biases destroy the integer nature of the inter-system double-differencing ambiguities. Two methods of direct rounding and parameter estimation are used to determine the ISB value. By analyzing data collected from Curtin University from 2015 to 2018, the phase and code inter-system bias (ISB) are related to the receiver type, firmware version and the selected overlapping frequency. Upgrade of receiver firmware version results in changes of ISB values. For example, the upgrade of Javad firmware in Dec, 15, 2017 causes the difference of 0.5 cycles ISB between BDS GEO and non-GEO satellites. By comparing the three dynamic models which include white noise process, random walk process, and random constant in the parameter estimation method, the ISB determined by the random constant model is consistent with the value obtained by the direct rounding method. After the calibration of ISBs, the performances of tightly combined positioning are assessed. The success rate of ambiguity resolution and accuracy of positioning for the tight combination (TC) are significantly improved in comparison with that for the loose combination (LC) over short baselines. For L5/E5a, on which only few satellites can be observed, the maximum increase in success rates of ambiguity resolution can reach 31.7%, i.e., from 54.9% of LC to larger than 86.6% of TC, and the positioning accuracies can even be increased by 0.13 m, i.e., from 0.208 of LC to 0.074 m of TC in East direction for the mix-receiver TRIMBLE NETR9-SEPT POLARX4 in 2018.