Fast solar radiation pressure modelling with ray tracing and multiple reflections

Physics based SRP (Solar Radiation Pressure) models using ray tracing methods are powerful tools when modelling the forces on complex real world space vehicles. Currently high resolution (1?mm) ray tracing with secondary intersections is done on high performance computers at UCL (University College London). This study introduces the BVH (Bounding Volume Hierarchy) into the ray tracing approach for physics based SRP modelling and makes it possible to run high resolution analysis on personal computers. The ray tracer is both general and efficient enough to cope with the complex shape of satellites and multiple reflections (three or more, with no upper limit). In this study, the traditional ray tracing technique is introduced in the first place and then the BVH is integrated into the ray tracing. Four aspects of the ray tracer were tested for investigating the performance including runtime, accuracy, the effects of multiple reflections and the effects of pixel array resolution.Test results in runtime on GPS IIR and Galileo IOV (In Orbit Validation) satellites show that the BVH can make the force model computation 30–50 times faster. The ray tracer has an absolute accuracy of several nanonewtons by comparing the test results for spheres and planes with the analytical computations. The multiple reflection effects are investigated both in the intersection number and acceleration on GPS IIR, Galileo IOV and Sentinel-1 spacecraft. Considering the number of intersections, the 3rd reflection can capture $99.12\%\text{,}99.14\%$, and $91.34\%$ of the total reflections for GPS IIR, Galileo IOV satellite bus and the Sentinel-1 spacecraft respectively. In terms of the multiple reflection effects on the acceleration, the secondary reflection effect for Galileo IOV satellite and Sentinel-1 can reach 0.2?$\text{nm}/{\text{s}}^2$ and 0.4?$\text{nm}/{\text{s}}^2$ respectively. The error percentage in the accelerations magnitude results show that the 3rd reflection should be considered in order to make it less than $0.035\%$. The pixel array resolution tests show that the dimensions of the components have to be considered when choosing the spacing of the pixel in order not to miss some components of the satellite in ray tracing. This paper presents the first systematic and quantitative study of the secondary and higher order intersection effects. It shows conclusively the effect is non-negligible for certain classes of misson.     

基于GPS和里程计的快速行进间粗对准方法

针对车载捷联惯导系统(SINS),提出一种快速行进间粗对准方法。将捷联矩阵分解为3个矩阵相乘,利用GPS提供的载体在导航坐标系的速度和里程计提供的载体在载体坐标系的速度构建不共线向量,求解初始载体坐标系相对惯性坐标系的常值转换矩阵,进而求得初始姿态矩阵。该对准方法对载车的唯一要求是在对准过程中做一个转弯机动。与现有对准方法相比,该方法没有用到加速度计信息。仿真结果表明,该方法能在1 min之内完成粗对准,采用零偏为0.1(°)/h的低精度陀螺,对准误差小于0.3°。     

机载对地观测飞行轨迹设计与对比分析

针对机载对地观测载机长时间匀速直线飞行时分布式位置姿态测量系统(POS)姿态误差随时间积累的问题,基于载机有效机动能够提高分布式POS系统可观测度进而提高系统估计精度的思想,设计并对比多种机动方式下分布式POS系统的估计精度,并对机动后分布式POS系统进入测区前的滤波估计稳定时间、直线飞行成像时长和全球定位系统(GPS)基线长度进行了测试。仿真结果表明,设计的飞行轨迹能够提高成像段运动参数的测量精度,可为机载对地观测最优飞行轨迹的选择和设计提供理论指导。     

一种基于UKF滤波器的卫星定位算法

针对基于通用最小二乘迭代法的卫星定位解算方法在定位精度和鲁棒性等方面的不足,提出了一种基于UKF(Unscented Kalman Filter)滤波器的卫星定位解算方法。该方法直接利用GPS接收机测得的伪距和多普勒频移作为观测量,对接收机的位置和速度进行估计。经过工程上的实际验证表明:与基于通用的最小二乘迭代法的卫星定位解算方法相比,所提出的基于UKF滤波器的卫星定位解算方法的定位结果有更高的精度和更强的鲁棒性。     

实现海洋测高卫星厘米级定轨精度的策略分析

基于对多种飞机的机载GPS测量实践,提出我国海洋二号（HY-2）卫星实现厘米级精度的星载GPS定轨测量的基本要求：1）选择适合天线,捕获多颗在视GPS卫星;2）注重天线安装位置,减弱多路径效应影响;3）选用适合的GPS信号接收机,确保星载GPS测量数据优质。     

ISO的GPS基本概念与体系

概要介绍了ISO的产品几何技术规范产生的原由、基本概念与体系的构成、主要内容和今后的发展方向 ,以及我国采用产品几何技术规范的情况。     

GLONASS-based precise point positioning and performance analysis

Current precise point positioning (PPP) techniques are mainly based on GPS which has been extensively investigated. With the increase of available GLONASS satellites during its revitalization, GLONASS observations were increasingly integrated into GPS-based PPP. Now that GLONASS has reached its full constellation, there will be a wide interest in PPP systems based on only GLONASS since it provides a PPP implementation independent of GPS. An investigation of GLONASS-based PPP will also help the development of GPS and GLONASS combined PPP techniques for improved precision and reliability. This paper presents an observation model for GLONASS-based PPP in which the GLONASS hardware delay biases are addressed. In view of frequently changed frequency channel number (FCN) for GLONASS satellites, an algorithm has been developed to compute the FCN for GLONASS satellites using code and phase observations, which avoids the need to provide the GLONASS frequency channel information during data processing. The observation residuals from GLONASS-based PPP are analyzed and compared to those from GPS-based PPP. The performance of GLONASS-based PPP is assessed using data from 15 globally distributed stations.     

Calibration of Envisat radar altimeter over Lake Issykkul

This study presents the results of calibration/validation (C/V) of Envisat satellite radar altimeter over Lake Issykkul located in Kyrgyzstan, which was chosen as a dedicated radar altimetry C/V site in 2004. The objectives are to estimate the absolute altimeter bias of Envisat and its orbit based on cross-over analysis with TOPEX/Poseidon (T/P), Jason-1 and Jason-2 over the ocean. We have used a new method of GPS data processing in a kinematic mode, developed at the Groupe de Recherche de Geodesie Spatiale (GRGS), which allows us to calculate the position of the GPS antenna without needing a GPS reference station. The C/V is conducted using various equipments: a local GPS network, a moving GPS antenna along the satellites tracks over Lake Issykkul, In Situ level gauges and weather stations. The absolute bias obtained for Envisat from field campaigns conducted in 2009 and 2010 is between 62.1 and 63.4 ± 3.7 cm, using the Ice-1 retracking algorithm, and between 46.9 and 51.2 cm with the ocean retracking algorithm. These results differ by about 10 cm from previous studies, principally due to improvement of the C/V procedure. Apart from the new algorithm for GPS data processing and the orbit error reduction, more attention has been paid to the GPS antenna height calculation, and we have reduced the errors induced by seiche over Lake Issykkul. This has been assured using cruise data along the Envisat satellite track at the exact date of the pass of the satellite for the two campaigns. The calculation of the Envisat radar altimeter bias with respect to the GPS levelling is essential to allow the continuity of multi-mission data on the same orbit, with the expected launch of SARAL/Altika mission in 2012. Implications for hydrology in particular, will be to produce long term homogeneous and reliable time series of lake levels worldwide.     

GPS载波相位时间传递软件的实现

研究了基于精密单点定位原理的GPS载波相位时间传递方法的数学模型和数据处理方法,设计了相应的程序流程,并基于Visual C++和Fortran混合编程的方式编写了相应的计算程序GPSCPTT V1.0。采用实际观测数据进行实验验证,结果表明,软件可以实现高精度的时间传递。     

Galileo系统与GPS卫星定位系统相位组合观测值的模型研究

在介绍Galileo系统空间信号的基础上,以模糊度保持整数为前提,给出了Galileo系统的4个频率载波与GPS L2载波的组合相位观测值的定义,并对有关误差影响加以分析,最后根据一定的组合标准论述了具有相应特性的组合观测值,并给出一些典型的组合.