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.     

基于双滤波器的高精度目标跟踪估计方法

针对连续推力的合作航天器，采用双重无迹卡尔曼滤波(DUKF)算法估计其状态和加速度。通过状态滤波器和参数滤波器的配合，提升滤波精度，完成运动状态和参数的估计，从而实现合作目标的运动轨迹跟踪。与合作航天器相比，非合作航天器存在大小未知、发生时刻未知的机动，无法获得加速度，且信息获取和运动状态的估计难度大。针对非合作航天器，基于简化的相对运动方程，结合天基平台获得目标的观测信息，采用两个扩展卡尔曼滤波(EKF)及基于半正焦弦的机动检测策略实现多未知脉冲机动的运动状态的估计。仿真结果表明：相比于无迹卡尔曼滤波（UKF），DUKF在对合作航天器的状态和加速度估计方面具有更快的收敛速度和更高的滤波精度；对于存在未知机动的非合作航天器，通过对比验证机动检测策略与滤波器切换策略相结合的方法的有效性，该方法能够检测到多次机动并且减少误判。     

速率模式飞轮姿态控制系统飞轮组合平稳切换方法

对速率模式飞轮姿态控制系统的能控性、稳定性和飞轮组合切换过程进行了研究，揭示了在飞轮组合切换过程中引起卫星姿态抖动的原因，提出了通过对残余角动量进行卸载并将卸载函数同时前馈到系统输入的方法，以保持卫星姿态的稳定度。从理论上证明了可引入残余角动量进行卸载，将卸载函数同时前馈到系统输入，并且不改变原系统的能控性和稳定性。给出了某卫星飞轮切换应用实例，说明了飞轮组合平稳切换方法的有效性。     

Modelling ignition probability with pairwise mixing-reaction model for flame particle tracking

Reduced order models for ignition analysis can offer insights into ignition processes and facilitate the combustor optimization. In this study, a Pairwise Mixing-Reaction (PMR) model is formulated to model the interaction between the flame particle and the surrounding cell mixture during Lagrangian flame particle tracking. Specifically, the model accounts for the two-way coupling of mass and energy between the flame particle and the surrounding shell layer by modelling the corresponding turbulent mixing, chemical reaction and evaporation process if present. The state of a flame particle, e.g., burnt, hot gas or extinguished, is determined based on particle temperature. This model can properly describe the ignition process with a spark kernel being initiated in a nonflammable region, which is of practical importance in certain turbine engines and has not been rigorously accounted for by the existing models based on the estimation of local Karlovitz number. The model is integrated into an ignition probability analysis platform and is demonstrated for a methane/air bluff-body flame with the flow and fuel/air mixing characteristics being extracted from a non-reacting simulation. The results show that for the spark location being at the extreme fuel-lean outer shear layer of the recirculation zone, PMR can yield ignition events with a significant number of active flame particles. The mechanisms for the survival of the initial flame particles and the entrainment of the survived flame particles into the recirculation zone are analyzed. The results also show that the ignition probability map from PMR agrees well with the experimental observation: a high ignition probability in the shear layer of the recirculation zone near the mean stoichiometric surface, and low ignition probabilities inside the recirculation zone and the top stagnation region of the recirculation zone. The parametric study shows that the predicted shape of the ignition progress factor and ignition probability is in general insensitive to the model parameters and the model is adequate for quantifying the regions with high ignition probabilities.     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

Implementation of tactical maneuvers with maneuver libraries

This study creates and combines the general maneuver libraries for fixed-wing aircraft to implement tactical maneuvers. First, the generalized maneuver libraries are established by analyzing the characteristics of tactical maneuvers required in battlefields. The 7th order polynomial is applied to both the creation of the maneuver libraries and the generation of the trajectories or flight paths for modal inputs. To track the desired trajectory, we design the Attitude Command Attitude Hold (ACAH) system and the Rate Command Rate Hold (RCRH) system using Model Following Controller (MFC). Moreover, the Line-of-Sight (LOS) guidance law is designed. In particular, the CONDUIT® is employed to optimize the gains so that the control systems meet the aircraft Handling Qualities (HQ) criteria. Finally, flight simulations are performed for the longitudinal loop, immelmann-turn, and climb-slalom-descent maneuvers to verify that tactical aggressive maneuvers are realizable via the combination of maneuver libraries. This study can contribute to the development of flight techniques for aircraft tactical maneuvers and to the revision of air force operational manuals.     

增广误差模型算法在目标跟踪中的应用

针对目标机动运行过程中,滤波模型与机动状态模型失配的问题,提出了一种新的增广状态误差滤波模型。不同于现有增广方案,该模型从模型失配所致状态滤波误差的角度出发,将状态估计误差增广为一状态量,通过滤波估计后用其校正原状态量。算法分析表明,该增广滤波模型具有自适应调节多重渐消因子的等效特性,增强了对目标的跟踪能力。基于该增广状态误差滤波模型,给出了滤波算法设计并进行了仿真实验。实验结果表明,基于该模型的滤波算法在对机动目标进行跟踪时具有更强的鲁棒性。     

一种小天体旋转参数估计方法

小天体旋转参数是科学数据，也是小天体测绘，着陆导航的基础数据.研究一种在小天体探测接近段过程中使用的基于图像上特征点跟踪和扩展卡尔曼滤波器的小天体旋转参数估计方法.该方法首先建立小天体旋转关系模型，表示小天体在相机坐标系中的姿态变化；然后定义小天体旋转的状态方程，推到了扩展卡尔曼滤波器的计算过程.通过对观测图像序列上的特征点跟踪，利用扩展卡尔曼滤波器方法得到小天体旋转轴指向和自转角速度估计.实验中，在仿真相机与小天体相机100 km距离上，分析了相机坐标系小天体坐标系之间的四元数初值，图像上特征点跟踪精度，相机的观测指向等因素对小天体旋转参数估计精度的影响.实验结果表明，提出的基于图像特征点跟踪和扩展卡尔曼滤波器的小天体旋转参数估计方法能够得到具有较高精度的估计结果.     

折叠式变体飞行器轨迹优化及控制分析

针对临近空间高超声速飞行器存在的问题，设计了一种折叠翼飞行器，可以通过折叠机翼来适应各种飞行状态，保持最优的气动特性。并针对临近空间滑翔式高超声速的特点，采用高斯伪谱法对固定翼飞行器和折叠翼飞行器的轨迹优化，通过将折叠翼飞行器与传统固定翼飞行器在射程能力、规避热流能力方面进行对比，提出了一种综合目标的轨迹优化思想。设计的折叠翼飞行器相比传统固定翼飞行器性能更加优越，更适合临近空间环境，提高了17.67%的航程，减少了热流率峰值的35.72%，并通过控制系统的设计和仿真加以验证，仿真结果表明变体飞行器机动能力相比固定翼飞行器有了显著的提高。     

航天器姿态预设性能控制方法综述

针对存在不确定惯量矩阵及其它未知不确定性(如执行器故障、航天器结构发生突变等)下的航天器姿态控制问题，综述了预设性能控制及其应用的研究进展。首先阐释了预设性能控制方法的基本内涵及其关键步骤；然后分析总结了现有航天器姿态控制以及预设性能控制研究的基本概况与发展趋势；最后面向未来复杂化、智能化的空间任务对航天器控制系统提出的新需求，提出航天器预设性能控制值得研究的问题和方向。