On Lambert’s problem and the elliptic time of flight equation: A simple semi-analytical inversion method

In the two-body model, time of flight between two positions can be expressed as a single-variable function and a variety of formulations exist. Lambert’s problem can be solved by inverting such a function. In this article, a method which inverts Lagrange’s flight time equation and supports the problematic 180°$180°$ transfer is proposed. This method relies on a Householder algorithm of variable order. However, unlike other iterative methods, it is semi-analytical in the sense that flight time functions are derived analytically to second order vs. first order finite differences. The author investigated the profile of Lagrange’s elliptic flight time equation and its derivatives with a special focus on their significance to the behaviour of the proposed method and the stated goal of guaranteed convergence. Possible numerical deficiencies were identified and dealt with. As a test, 28 scenarios of variable difficulty were designed to cover a wide variety of geometries. The context of this research being the orbit determination of artificial satellites and debris, the scenarios are representative of typical such objects in Low-Earth, Geostationary and Geostationary Transfer Orbits. An analysis of the computational impact of the quality of the initial guess vs. that of the order of the method was also done, providing clues for further research and optimisations (e.g. asteroids, long period comets, multi-revolution cases). The results indicate fast to very fast convergence in all test cases, they validate the numerical safeguards and also give a quantitative assessment of the importance of the initial guess.     

基于核部分非负矩阵分解的亚像元级地物光谱分析

为了进一步提高亚像元级地物的光谱分析精度,提出了一种基于核部分非负矩阵分解(Kernel Protection Non-negative Matrix Factorization,KPNMF)的非线性解混算法.首先通过基于凸面几何理论的端元提取方法提取纯像元端元候选像素集合,然后根据候选像素的空间纯度指数判断纯像元端元.在纯像元端元信息已知的条件下,利用核方法对部分非负矩阵分解(Protection Non-negative Matrix Factorization,PNMF)进行推广,构造相应的目标函数,推导迭代求解过程,分解求得亚像元端元光谱和所有端元的丰度.试验结果表明,提出的解混算法具有良好的非线性分解能力,解混结果优于线性解混算法.     

无反作用力矩空间机器人轨迹跟踪控制

针对目前关节驱动方式的不足,提出一种采用无反作用力矩驱动方式的空间机器人设计概念。系统平台与各节机械臂上均安装一组金字塔构型的控制力矩陀螺（CMGs）作为力矩执行机构,各臂间由自由球铰连接。采用Rodriguez参数描述平台姿态和机械臂角位移,利用Kane方程建立系统动力学模型。在此基础上,设计了渐近稳定的轨迹跟踪控制律,使得平台姿态和机械臂角位移跟踪期望运动轨迹;并设计了CMGs操纵律,使之准确输出期望控制力矩。此外,研究了机械臂工作空间到关节空间的轨迹规划算法,使得所设计的控制律也可应用于工作空间的轨迹跟踪控制。由三关节系统的仿真结果,验证了无反作用力矩设计概念的可行性和轨迹规划算法的有效性。     

基于HLA的航天飞行任务联合仿真系统设计

简要说明了HLA的基本特点,描述了一个已经实现的航天测控仿真系统,并将其结构思想与HLA进行了对照。在此基础上,提出了一个依据HLA标准,对航天员座舱仿真、火箭和飞船仿真、测控网仿真等进行集成的航天飞行任务联合仿真系统框架,分析了各个仿真联邦成员的功能需求和信息接口,并着重介绍了针对任务要求的时间管理、数据通信等关键RTI机制的设计和实现方法,以及参照FEDEP模型和借助商业成品软件工具的系统开发策略。     

基于小波变换的单脉冲雷达目标RCS测量方法研究

随着现代信号处理技术的发展,对非平稳信号分析和处理的小波分析技术已成功应用于雷达目标特性分析领域,大功率单脉冲雷达作为我国航天测控网的主干设备,具有一定的目标特性识别能力。本文主要针对脉冲雷达RCS测量原理,讨论了基于小波变换的单脉冲雷达空间目标RCS测量方法,提出应发挥窄带低分辨率雷达在未来空间目标识别中的作用。     

空间碎片材料检测技术研究

为了对空间碎片材料的成分进行快速检测,分析了空间碎片的来源和构成,并对空间碎片主要成分的国内外检测技术手段进行了详细对比,采用了一种基于LIBS的网络化空间碎片检测方案。该方案具有分析速度快、测量精度和灵敏度高、可多元素同时测量、空间适应性好且无需样品预处理等特点,同时可实现检测流程全自动化和多线程化。     

空间太阳能电站微波能量传输验证方案设计

微波能量传输技术作为空间太阳能电站（SpaceSolarPowerStation,SSPS）的关键技术之一,目前的研究和验证工作均集中在各单项技术的突破和验证,缺乏针对SSPS系统特点的全面优化设计。文章根据SSPS的工作模式给出了全面验证空间太阳能电站微波能量传输的验证系统方案设计,对收发天线进行了一体化设计,利用了幅度近似高斯分布的发射阵列场分布设计和低反射的接收整流阵列设计,以高精度来波方向测量和高精度移相控制为波束指向控制的技术途径。对验证系统的波束收集效率进行了分析,收集效率可达94.2%,比传统均匀分布系统高出17.6%。验证系统可从系统规模缩比、波束扫描范围、发射天线口径场分布、整流天线处功率密度、反向波束控制方法等方面模拟SSPS微波能量传输工作模式,推动SSPS系统技术的发展。     

微波能量传输系统设计及试验

面向空间太阳能电站应用,进行了固态体制微波能量传输技术研究。针对能量传输波束扩散导致收集效率低的问题,研究了基于人工媒质理论设计的完美匹配层的能量接收整流表面,通过调节人工媒质单元的结构参数实现天线输出阻抗与整流电路输入阻抗的共轭匹配,同时抑制整流电路高次谐波,省去原有匹配及低通滤波器,简化电路结构、实现高效微波能量吸收与转换。以空间太阳能电站规定的波束中心传输微波功率密度限制作为能量接收整流表面设计的约束条件,设计能量接收整流表面,结合固态体制微波能量发射端,搭建5.8GHz小规模微波能量传输系统开展了地面试验验证,实测结果显示整流表面能量转换效率最高为57.7%。此次试验验证了先进的固态能量传输试验系统,为空间太阳能地面缩比试验及未来空间太阳能电站的建设提供了技术支撑。     

Yarkovsky-Schach effect on space debris motion

The Yarkovsky-Schach effect is a small perturbation affecting Earth satellites and space debris illuminated by the Sun. It was first applied to the orbit of LAGEOS satellites as an explanation of the residuals in orbital elements. In this work, we carry out several numerical integration tests taking into consideration various orbit and rotation parameters, in order to analyse this effect in a broader context. The semi-major axis variations remain small and depend on the spin axis attitude with respect to the Sun. We show that the force amplitude is maximised for orbits inclined with $i\approx$?20–30°. We also observe the influence on other orbital elements, notably on the orbit inclination. However, these effects are clearly observed only on long timescales; in our simulations, we propagated the orbits for 200?y. The Yarkovsky-Schach effect is thus confirmed to have a minuscule magnitude. It should be taken into account in studies requiring high-precision orbit determination, or on expanded timescales.     

Optimal Earth’s reentry disposal of the Galileo constellation

Nowadays there is international consensus that space activities must be managed to minimize debris generation and risk. The paper presents a method for the end-of-life (EoL) disposal of spacecraft in Medium Earth Orbit (MEO). The problem is formulated as a multiobjective optimisation one, which is solved with an evolutionary algorithm. An impulsive manoeuvre is optimised to reenter the spacecraft in Earth’s atmosphere within 100?years. Pareto optimal solutions are obtained using the manoeuvre $\mathrm{\Delta }v$ and the time-to-reentry as objective functions to be minimised. To explore at the best the search space a semi-analytical orbit propagator, which can propagate an orbit for 100?years in few seconds, is adopted. An in-depth analysis of the results is carried out to understand the conditions leading to a fast reentry with minimum propellant. For this aim a new way of representing the disposal solutions is introduced. With a single 2D plot we are able to fully describe the time evolution of all the relevant orbital parameters as well as identify the conditions that enables the eccentricity build-up. The EoL disposal of the Galileo constellation is used as test case.