高纬度地区GPT2w模型的适应性分析

GPT2w模型是现有精度最高的天顶对流层模型，但是应用在高纬度地区时存在较大误差.为更好地保障卫星导航定位系统在高纬度地区的高精度应用，评定了GPT2w模型在高纬度地区的精度，获取天顶对流层湿延迟、干延迟和总延迟，探讨了GPT2w模型改正对精密单点定位的影响.试验结果表明: GPT2w模型在高纬度地区的精度为厘米级，优于其在中低纬度地区的精度；南北极地区天顶对流层呈现明显季节变化特征和区域一致性特征，夏季天顶对流层总延迟高于冬季，北极地区天顶对流层湿延迟明显高于南极地区，北极地区天顶对流层随季节的变化幅度大于南极地区.PPP试验结果表明，GPT2w模型能够有效改善定位精度，适应高纬度地区的高精度定位需求.      

Kalman滤波估算电离层延迟的一种优化方法

频间偏差（Inter Frequency Bias,IFB）通常会给电离层延迟的解算带来误差.目前从电离层延迟中消除频间偏差的方法是基于GPS双频观测数据建立垂直电离层模型,利用卡尔曼滤波实时估算电离层模型系数和频间偏差.然而滤波过程中的测量噪声协方差矩阵没有考虑系统观测量之间的相关性,这会导致滤波模型不准确,进而影响最后求解的电离层延迟的准确性.本文选取了美国19个参考站的GPS双频观测数据,利用卡尔曼滤波实时估算电离层模型系数以及频间偏差.在滤波过程中,通过将先验频间偏差的估计方差引入测量噪声方差,实现对测量噪声协方差矩阵的优化.计算结果表明,优化后得到的卫星频间偏差与欧洲定轨中心（Center for Orbit Determination in Europe,CODE）得到的频间偏差更接近.将优化后的电离层延迟代入伪距解算,得到的位置误差的标准差在东向和天顶向分别下降了12.5%和15.4%,天顶向误差平均值下降了17.6%,定位精度得到提高.      

基于N-BEATS的单站对流层天顶总延迟预报

高精度对流层延迟先验值有助于加速精密单点定位的快速收敛。基于高精度高分辨率气象数据库，采用深度学习N BEATS算法，进行了单站对流层天顶总延迟的预报试验。试验选取了9个IGS跟踪站，试验弧段从2002年1月至2019年6月共185a。首先基于N BEATS算法，设计了3种预报策略，然后基于前175a针对不同预报策略进行模型训练，并对最后365d的对流层天顶总延迟进行预报。试验结果表明，以该气象数据库为基准，12h以内预报弧段的预报残差均值量级大多可达亚毫米，2h、4h、6h的预报残差的标准差分别约为5mm、9mm、13mm。     

基于多模型最优融合的双星定位系统一体化精密定轨方法

针对卫星动力学模型复杂且不准,考虑到卫星定轨中待估参数在时间和空间上的相关性,提出了一种基于多模型最优融合的双星定位系统一体化参数建模的近地卫星精密定轨新方法.利用节点自由分布B样条描述卫星运动,实现了对卫星粗略动力学模型的抑制作用;同时结合双星观测模型,使该方法转化为关于求解卫星轨道样条表示参数和定轨系统误差的多模型融合的非线性优化问题;通过引入模型结构确定最优融合权值的选取准则,在最小二乘准则下,采用非线性最优化方法搜寻样条的最优节点分布,得到了待估参数的最优估计,完成了近地卫星的精密定轨.理论分析和仿真计算表明,该方法确实有效,不仅提高了卫星的定轨精度,而且使状态估计的结构更加稳定.      

附加先验轨道约束的LEO几何法实时精密定轨验证分析

低轨卫星的实时精密定轨能够极大拓展其完成复杂科学任务的能力,例如实时环境监测、机动控制和卫星自主导航等.本文根据几何法实时精密定轨模型,提出了附加LEO先验轨道约束从而改善实时定轨的精度、收敛速度和稳定性的构想.分别采用广播星历、超快速星历预报部分和实时精密星历,设计了6种实时定轨方案,并利用Swarm-A,B,C星7天的观测数据进行方案验证与分析.结果表明,使用广播星历、IGU和IGC星历的方案精度递增,附加先验轨道约束能够进一步提升精度.使用IGC星历并附加标准差为1m的先验轨道约束后,在径向、切向和法向的定轨精度分别达到6.12cm,5.55cm和4.98cm.此外,附加先验轨道约束能够显著提升收敛速度,使用IGC星历平均收敛时间约为31min,附加标准差为1m的先验轨道约束后收敛仅需约4min.      

基于Bayes统计模型的卫星融合定轨方法研究

根据目前天基导航系统现状,结合中国对低、中轨卫星精密定轨的要求,给出了天地基信息融合定轨的原理;结合卫星待估融合参数的先验信息,提出了基于Bayes统计模型的卫星精密定轨方法;在卫星观测的线性化融合模型中引入观测噪声,利用概率估计融合模型,根据Bayes理论进行卫星状态改进量的最大后验估计,并分析了Bayes估计方法的定轨精度;依据期望融合的待估改进量方差最小规则建立了相应的参数求解算法;最后以导航融合测控系统中测距和测速数据的融合定轨为例进行了仿真实验,表明该融合方法能够得到很好的定轨效果。     

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

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

基于单频星载GPS数据的低轨卫星精密定轨

为满足搭载单频GPS接收机低轨卫星的精密定轨需求以及深化单频定轨研究,文中解决了单频星载GPS数据的周跳探测问题,并利用“海洋二号”（HY-2A）卫星及“资源三号”（ZY-3）卫星的单频星载GPS实测数据采用两种方法确定了二者的简化动力学轨道,并通过观测值残差分析、与双频精密轨道比较、激光测卫数据检核等方法对所得轨道精度进行评定。结果表明,在不考虑电离层延迟影响的情况下,HY-2A卫星定轨精度为2~3dm,ZY-3卫星为1m左右;而采用半和改正组合消除电离层延迟一阶项影响后,二者定轨精度均显著提高,HY-2A卫星三维精度提高至1dm左右,ZY-3卫星提高至1~2dm。文章的研究成果表明,搭载单频GPS接收机的低轨卫星也可获得厘米级的定轨精度。     

基于双星系统定轨的数学方法

通过对双星系统定轨原理的分析,结合卫星的动力学模型,采用联合多时刻观测系统的数据融合方法,得到了一种解算卫星轨道状态的解析方法,同时还给出了误差传递方式。该方法不需要标称轨道,计算方便。理论分析和仿真计算表明,方法确实可行,且定轨精度有所提高。     

基于移动站的转发式地面站设备时延标校方法

在转发式卫星测定轨系统中,基于伪码测距原理的星地距离测量是实现卫星精密定轨和高精度时间比对的基础。为获得高精度的星地距离,需要将地面站设备时延从伪码测距值中精确扣除。在转发式卫星测轨原理的基础上,提出了基于移动站的转发式地面站设备时延标校方法,实现了对转发式地面站设备时延的标校,标校精度能够优于0.5ns,对提高转发式卫星定轨精度和卫星双向时间比对精度具有重要作用。     

Orbit determination of BDS-3 satellite based on regional ground tracking station and inter-satellite link observations

The BeiDou global navigation satellite system (BDS-3) has established the Ka-band inter-satellite link (ISL) to realize a two-way ranging function between satellites, which provides a new observation technology for the orbit determination of BDS-3 satellites. Therefore, this study presents a BDS satellite orbit determination model based on ground tracking station (GTS) observations and ISL ranging observations firstly to analyze the impact of the ISL ranging observations on the orbit determination of BDS-3 satellites. Subsequently, considering the data fusion processing, the variance component estimation (VCE) algorithm is applied to the parameter estimation process of the satellite orbit determination. Finally, using the measured data from China’s regional GTS observations and BDS-3 ISL ranging observations, the effects of ISL ranging observations on the orbit determination accuracy of BDS-3 satellites are analyzed. Moreover, the impact of the VCE algorithm on the fusion data processing is evaluated from the aspects of orbit determination accuracy, Ka-band hardware delay parameter stability, and ISL ranging observation residuals. The results show that for China’s regional GTSs, the addition of BDS-3 ISL ranging observations can significantly improve the orbit determination accuracy of BDS-3 satellites. The observed orbit determination accuracy of satellite radial component is improved from 48 cm to 4.1 cm. In addition, when the initial weight ratio between GTS observations and ISL ranging observations is not appropriate, the various indicators which include orbit determination accuracy, ISL hardware delay, and ISL observation residuals were observed to have improved after the adjustment of the VCE algorithm. These results validate the effectiveness of the VCE algorithm for the fusion data processing of the GTS observations and ISL ranging observations.     

Modeling and analysis of Earth’s gravity field measurement performance by inner-formation flying system

The Earth’s gravity field can be measured with high precision by constructing the purely gravitational orbit of the inner-satellite in Inner-formation Flying System (IFS), which is independently proposed by Chinese scholars and offers a new way to carry out gravity field measurement by satellite without accelerometers. In IFS, for the purpose of quickly evaluating the highest degree of recovered gravity field model and geoid error as well as analyzing the influence of system parameters on gravity field measurement, an analytical formula was established by spectral analysis method. The formula can reflect the analytical relationship between gravity field measurement performance and system parameters such as orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and total measurement time. This analytical formula was then corrected by four factors introduced from numerical simulation of IFS gravity field measurement. By comparing computation results from corrected analytical formula and the actual gravity field measurement performance by CHAMP, the correctness and rationality of this analytical formula were verified. Based on this analytical formula, the influences of system parameters on IFS gravity field measurement were analyzed. It is known that gravity field measurement performance is a monotone decreasing function of orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and the reciprocal of total measurement time. There is a match relationship between the inner-satellite orbit determination error and residual disturbances, in other words, the change rate of gravity field measurement performance with one of them is seriously restricted by their relative size. The analytical formula can be used to quantitatively evaluate gravity field measurement performance fast and design IFS parameters optimally. It is noted that the analytical formula and corresponding conclusions are applied to any gravity satellite which measures gravity field by satellite perturbation orbit.     

Tropospheric wet refractivity tomography using multiplicative algebraic reconstruction technique

Algebraic reconstruction techniques (ART) have been successfully used to reconstruct the total electron content (TEC) of the ionosphere and in recent years be tentatively used in tropospheric wet refractivity and water vapor tomography in the ground-based GNSS technology. The previous research on ART used in tropospheric water vapor tomography focused on the convergence and relaxation parameters for various algebraic reconstruction techniques and rarely discussed the impact of Gaussian constraints and initial field on the iteration results. The existing accuracy evaluation parameters calculated from slant wet delay can only evaluate the resultant precision of the voxels penetrated by slant paths and cannot evaluate that of the voxels not penetrated by any slant path. The paper proposes two new statistical parameters Bias and RMS, calculated from wet refractivity of the total voxels, to improve the deficiencies of existing evaluation parameters and then discusses the effect of the Gaussian constraints and initial field on the convergence and tomography results in multiplicative algebraic reconstruction technique (MART) to reconstruct the 4D tropospheric wet refractivity field using simulation method.     

Precision orbit determination for TOPEX/Poseidon using TDRSS doppler tracking data

Precision orbit determination on the TOPEX/Poseidon (T/P) altimeter satellite is now being routinely achieved with sub-5cm radial and sub-15 cm total positioning accuracy using state-of-the-art modeling with precision tracking provided by a combination of: (a) global Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS), or (b) the Global Positioning System (GPS) Constellation which provides pseudo-range and carrier phase observations. The geostationary Tracking and Data Relay Satellite System (TDRSS) satellites are providing the operational tracking and communication support for this mission. The TDRSS Doppler data are of high precision (0.3 mm/s nominal noise levels). Unlike other satellite missions supported operationally by TDRSS, T/P has high quality independent tracking which enables absolute orbit accuracy assessments. In addition, the T/P satellite provides extensive geometry for positioning a satellite at geostationary altitude, and thus the TDRSS-T/P data provides an excellent means for determining the TDRS orbits. Arc lengths of 7 and 10 days with varying degrees of T/P spacecraft attitude complexity are studied. Sub-meter T/P total positioning error is achieved when using the TDRSS range-rate data, with radial orbit errors of 10.6 cm and 15.5 cm RMS for the two arcs studied. Current limitations in the TDRSS precision orbit determination capability include mismodeling of numerous TDRSS satellite-specific dynamic and electronic effects, and in the inadequate treatment of the propagation delay and bending arising from the wet troposphere and ionosphere.     

基于GPS的地球资源卫星精密星历外推新方法

随着技术的发展,通过星载GPS接收机直接确定卫星星历成为卫星定位的一个重要手段.GPS接收机获取的卫星星历数据是某一时刻的瞬时状态,要获取连续的卫星星历数据还需要进一步处理.常用的处理方法有几何法与动力学法.在GPS接收机给定瞬时星历频率较低的情况下,几何法的计算误差比较大,特别是只有一组瞬时星历时,无法用几何法进行轨道的外推.在分析地球资源卫星轨道特点的基础上,提出一种新的轨道缩减动力模型,该模型将卫星运动在直角坐标系中分解为简谐运动,利用模型实现了轨道外推的算法.通过试验验证,该算法可以达到较高的精度.      

北斗GEO卫星CEI相时延解算方法研究

连线端站干涉测量(connected element interferometry，CEI)是高精度测角技术，在中高轨卫星、月球及深空航天器定轨定位中有良好的应用前景。基于CEI技术特点，提出了一种新的测量方法，即在相干测距模式下利用测距音和载波信号作为信号源进行连线端站干涉测量。构建了CEI试验系统对北斗GEO卫星进行观测，利用相干测距模式下的下行信号解算群时延、相时延。利用北斗GEO卫星精密星历计算的时延理论值，对北斗GEO卫星CEI群时延和相时延结果进行评估。结果表明，相干测距模式下CEI群时延和相时延残差均值分别为0.47ns、0.08ns，标准差(3σ)分别4.2ns、0.13ns。该项研究验证了相干测距模式下CEI相时延解算的可行性，可为共位地球同步卫星精密相对定位、月球探测器CEI测量提供技术参考。     

一种新的星载GPS定轨滤波算法

针对传统EKF(Extended Kalman Filtering)算法应用于星载GPS(Global Positioning System)低轨卫星定轨时系统噪声方差初值不易确定的问题,提出了一种新的定轨滤波算法.该算法在非线性方程线性化过程中,在前一时刻滤波估值点进行线性化,从而得到扰动方程,并将该扰动方程引入到传统EKF进行滤波处理.该算法与传统EKF分别应用在星载GPS低轨卫星的定轨中,通过比较,结果表明改进的算法在一定程度上抑制了由于系统噪声方差阵选取偏差较大而引起的滤波发散现象,且对于系统噪声方差的初值选取有较强的鲁棒性.