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.     

一种地月转移轨道中途修正的制导算法

针对地月转移轨道中途修正问题，提出了一种求解修正速度增量的制导算法。该算法由初值设计和精确解求解两部分组成。首先，利用伪状态理论，通过简单迭代设计中途修正的初值，并通过Vinti预报方法修正了地球扁率的影响。然后，在求解精确解时，提出了一种基于伪状态理论的状态转移矩阵解析算法。该算法通过设计高精度的初值，降低了求解地月转移轨道中途修正问题的难度，而且避免了传统数值方法计算状态转移矩阵的复杂性。数值仿真结果表明，该算法可有效求解中途修正问题。     

复杂动态环境下无人飞行器动态避障近似最优轨迹规划

针对复杂动态环境下无人飞行器的动态障碍规避问题，基于合理假设建立了无人飞行器和动态障碍的运动学模型，并综合考虑无人飞行器飞行过程中的终端约束、控制输入约束、安全避障约束等，以能量最少为性能指标构建动态避障问题数学描述。之后，针对终端约束和控制输入约束，依据优化模型预测静态规划算法(OMPSP)生成初始轨迹；针对动态避障问题的不等式约束，引入松弛变量并结合滑模变结构控制方法设计松弛变量动力学，实现对一个、多个或同时多个动态障碍的安全规避；最后，依据有限时间微分动态规划(RHDDP)算法进行轨迹优化，获得满足上述各种约束并能规避动态障碍的近似最优轨迹。     

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

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

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

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

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

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

基于有限Fourier级数的航天器转移轨迹设计

针对航天器小推力转移轨迹的初始设计问题，利用基于三阶Fourier级数的设计方法实现了航天器小推力的多圈转移。同时，基于有限Fourier级数的形状法，对具有多个约束条件的小推力多圈转移轨迹进行了优化设计。选取了共面同轴同偏心率的初始和末端轨道位置，对所提出的方法进行了仿真验证。结果表明：与改进逆五阶多项式形状法相比，所提出的方法虽然增加了转移时间，但当转移圈数为5圈、有限Fourier级数的项数为10时，可减少将近75%的转移速度增量，同时大大减小了所需的最大推力加速度的值。     

反推力装置运动学与动力学仿真

为了研究叶栅式反推力装置各部件在工作过程中的运动学与动力学特性,根据机构运动原理对反推力装置进行简化并建立了运动学与动力学数学模型。以滑动整流罩位移与阻流门所受气动负荷为输入进行运动学与动力学仿真,得到了反推力装置各部件的位移、速度及受力特性曲线,并对比分析了在不同尺寸参数下各部件特征点的运动轨迹和反推力装置负载力变化。结果表明:运动学及动力学仿真结果与工程实际相符；反推力装置机构参数选择不合理时,各设计点会发生干涉现象并导致机构无法运动;机构参数变化对负载力最大值影响尤为突出,在阻流门AC段长度值增大6%,阻流门CB段长度值减小9%的情况下,负载力正向最大值将增大19.53%,负向最大值增大12.67%。研究方法及研究结果可为反推力装置运动学及动力学分析,以及为反推力装置机构优化设计提供参考。