基于准流形方法日地系L3点转移轨道设计

针对深空探测中轨道转移时间长且能量消耗较大的问题,提出基于准流形实现从地球停泊轨道到日地系L3点转移轨道的设计方法。在日地限制性三体问题模型下,在L1点或L2点Halo轨道上施加扰动推力,构造准流形,利用其非线性三体动力学特性,通过霍曼转移轨道与近地轨道进行拼接,使航天器进入准流形后能够无动力滑行到L3点附近区域。在准流形与L3点周期轨道交点,施加速度脉冲,使航天器进入相应周期轨道,从而完成轨道转移。仿真结果表明,利用该方法所得结果与基于不变流形的转移轨道相比,能将速度增量从4398m/s减少为4014m/s,并将转移时间从9年以上缩短到7.3年以内,有效地提高了航天器的工作效率。     

日地系L2点与日火系L1点转移轨道设计

将太阳-地球-火星-飞行器组成的四体问题分解成由太阳-地球-飞行器和太阳-火星-飞行器两个共面圆形限制性三体问题,设计日地系L2点与日火系L1点Lyapunov轨道之间的转移轨道,该转移轨道可以作为探测火星时的低能中间转移轨道.采用Richardson三阶近似解作为初始值,运用微分修正方法分别得到两个不同三体系统下拉格朗日点的精确Lyapunov轨道.基于Lyapunov轨道不变流形以及微分修正方法,设计了日地系L2点与日火系L1点间的转移轨道.将所得结果与基于Halo轨道不变流形设计的转移轨道进行了对比.结论表明:利用Lyapunov轨道不变形设计探火中间转移轨道相较于利用Halo轨道不变流形设计探火中间转移轨道在能量消耗以及飞行时间上都存在优势.     

基于电推进卫星飞往Halo轨道的转移轨道研究

利用电推进及轨道力学的特性实现节能优化,将限制性三体问题中的稳定不变流形与小推力轨道优化相结合,研究全电推进卫星从地球停泊轨道飞向日地拉格朗日L2点Halo轨道的低消耗转移轨道.航天器的转移轨道分为逃逸段、拼接段与无动力滑行段.在逃逸段卫星沿速度方向加速脱离地球引力,拼接段采用Radau伪谱法进行优化,使航天器以最短时间到达目标Halo轨道的稳定不变流形上,随后航天器电推进系统关机,沿稳定不变流形无动力滑行至目标轨道.基于雅克比积分常数给出拼接段轨道初始猜测值,以先提高切向方向航天器能量避免了全程优化离散点过多难以求解的问题.仿真结果表明,该方法收敛速度较快,对平动点工程任务的初期轨道特性计算具有实际意义.      

一种适用于平动点周期轨道初值计算的简化路径搜索修正法

在限制性三体问题中，路径搜索修正法是一种基于平动点周期轨道垂直穿越Poincare截面的几何对称性计算平面及空间平动点周期轨道近似初值的方法.采用路径搜索修正法的一种简化形式，在圆形限制性三体模型中，对地月系中几种典型的平面及空间周期轨道近似初值进行了计算.结果表明，该简化方法得到的周期轨道近似初值不唯一，由近似初值经微分修正得到的精确结果中通常同时存在Halo轨道和大幅值逆行轨道（DRO）.进一步分析表明，在某些临界初值下，精确结果中Halo轨道将消失，同时可能出现平面Lyapunov轨道及Vertical轨道.上述计算中，搜索初值与结果中轨道类型的对应关系值得进一步研究.      

基于不变流形的地-月L2点转移轨道优化设计

在地–月L2点月球中继卫星轨道转移设计中,采用高比冲、小推力的电推进器可以大大增加卫星的有效载荷比,但会增加轨道设计的难度。基于地球GEO轨道为初始轨道,地–月L2点的halo轨道为目标轨道,通过最优控制中的混合法及平动点轨道的不变流形,研究了作为拓展任务的利用地月系统不变流形的小推力变轨方案,可以有效简化转移的轨道设计。仿真结果表明:得到了任意推力情况下最节省推进剂燃料的推力方向控制方案,对月球中继卫星的轨道设计及其平动点轨道设计具有工程意义。     

地月空间的远距离逆行轨道族及其分岔研究

地月系统中存在着一类绕月逆行、高度稳定的轨道族,称为远距离逆行轨道族(DRO)。以圆型限制性三体问题(CR3BP)为动力学模型研究了DRO轨道族周边的动力系统结构。利用Broucke稳定性图寻找分叉点,判断分叉类型,基于数值延拓计算分岔后产生的一系列新轨道分支。分叉类型主要有切分叉与多倍周期分叉(从3倍周期开始),轨道维度包含平面轨道族与三维轨道族。计算新轨道族的特征,包括形状、周期、能量、稳定性、双曲流形结构等。探讨周期轨道的轨道周期与能量的关系,以几何化的方式展现分叉结构、多周期轨道的双曲流形结构等。该动力结构将为基于DRO轨道族的地月空间任务提供重要的理论支持。      

火卫一周期准卫星轨道及入轨分析

围绕火卫一的准卫星轨道（QSOs）因其具有良好的稳定性，是火卫一探测任务最为实用的轨道。在平面圆型限制性三体问题模型下，利用庞加莱截面和KAM环迭代方法探究了准卫星轨道的周期轨道族，并给出不同能量准卫星周期轨道的初始条件。针对火卫一周期准卫星轨道入轨，提出一种转移轨道设计方法:对准卫星周期轨道调整速度后进行反向积分，直至离开火卫一邻近区域，从而得到由火星环绕轨道向火卫一周期准卫星轨道的转移轨道，并调整转移轨道参数对燃料与时间消耗进行优化。研究结果表明，当周期准卫星轨道能量处于特定区间时，存在特定速度脉冲区间，可利用火卫一引力实现较少燃料消耗的轨道转移；在该速度脉冲区间中，通过选取较小的速度脉冲，可缩短转移时间。      

基于三体系统平动点轨道的环日全景任务轨道设计

在空间开展太阳观测是研究太阳活动周、太阳爆发、极端天气等事件起源的重要手段。环日全景探测计划是为实现从黄道面360°全方位观察太阳行星际空间而提出的。本文针对环日全景探测计划,构建了基于三体系统平动点低能量轨道的环日全景轨道部署方法。该方法以日–地L1/L2点Halo轨道幅值及Halo轨道离轨点为变量,以转移轨道飞行时间、入轨机动大小为评价指标,基于三体系统不变流形构建环日全景的转移轨道,并开展轨道优化设计。采用等高线图对设计变量及任务成本进行全局分析。仿真计算表明,轨道部署无法同时满足飞行时间最短与入轨机动最小的要求。设计了轨道机动约束条件下的最优飞行时间解,并给出了基于长三甲运载火箭的一箭双星发射及入轨方案。      

地-月L2点中继星月球近旁转移轨道设计

位于地月L2点周期轨道的中继星将首次为"嫦娥4号"月球背面着陆探测任务提供通信中继服务。中继星转移轨道设计是中继任务实施的关键环节。针对中继星转移轨道存在转移时间、近月点高度和halo轨道振幅等约束条件,系统研究了基于月球近旁的地月L2点转移轨道设计方法。首先基于限制性三体模型,分析了halo轨道族与着陆点可见性关系;然后将月球近旁转移轨道分为地月直接转移段和地月动平衡点附近周期轨道拟流形入轨转移段,采用带有状态约束的微分修正算法对这两段轨道进行拼接,得到了从地球附近至目标轨道族的月球近旁转移轨道;最后,针对南族halo轨道分析了halo轨道振幅和月球飞越高度对转移轨道设计的影响,以及转移轨道的入轨相位分布。仿真结果表明:月球近旁转移轨道设计方案具备工程上的可行性与优越性。该方案可以为实际工程任务和应用提供参考。     

小天体环的轨道动力学

随着天文观测的深入,发现不仅土星、木星等气态巨行星环绕有光环,某些小天体也环绕有光环。针对目前发现环的半人马型小天体女凯龙星（10199 Chariklo）,研究其引力场中粒子（视为质点）的轨道动力学,分析赤道椭率和小天体自转对环中粒子轨道运动的影响,通过庞加莱截面中的KAM环迭代,得到第一类周期轨道和1∶3共振周期轨道,并通过对比和分析得到小天体光环位置与小天体自转平运动共振的关系。研究结果表明:女凯龙星内环中的粒子最可能处在满足径向振荡幅值范围的第一类周期轨道及其附近的准周期轨道上,但不能排除处在1∶3共振周期轨道及其附近的准周期轨道上;外环中的粒子不可能处在1∶3共振周期轨道上,只可能处在第一类周期轨道及其附近的准周期轨道上。      

Four-body cislunar quasi-periodic orbits and their application to ballistic lunar transfer design

Examining the properties of quasi-periodic orbits provides insight into the Sun-perturbed environment in cislunar space. In this investigation, quasi-periodic trajectories and their properties are explored in the Sun-Earth-Moon four-body problem. Computation and the stability characteristics of families of invariant tori are detailed. Furthermore, this investigation offers a framework for construction of ballistic lunar transfer trajectories in the four-body problem. The framework leverages manifold trajectories to supply a set of initial conditions for construction of periapsis Poincaré maps. Periapsis maps reduce the dimensionality of the space and illuminate solutions of interest as a basis to produce feasible families of transfers. Through a continuation process, families of ballistic transfers and families of transfers that include powered lunar flybys are constructed. Ultimately, these solutions supply an initial guess for transition to the Sun-Earth-Moon ephemeris model.     

基于不变流形的高超声速飞行器动态面控制

针对存在未知气动参数的吸气式高超声速飞行器纵向运动控制问题,提出一种基于不变流形的自适应动态面控制方法.通过合理假设将高超声速飞行器纵向模型分解为弹道倾角回路和速度回路,分别实现对弹道倾角和速度参考指令的跟踪.弹道倾角回路被表达为严反馈形式进行控制器设计.采用基于不变流形的自适应方法实现了对未知参数的估计.所提出的自适应动态面控制方案能保证未知参数估计误差全局一致稳定和闭环系统全局有界稳定,且估计器和控制器设计不存在耦合,因此参数设计更加容易.仿真结果验证了该控制方法在参数估计方面的显著优势和良好的闭环系统性能.     

基于几何不变性和BP网络的二维目标识别算法

提出了一种基于几何不变性和BP网络的二维目标识别算法.该算法不仅能适应目标物体在旋转、缩放和平移变换(RST变换)下的不变性识别,而且能适应仿射及射影变换下的不变性识别.算法通过对目标物体边缘点进行规格化和对规格化后的边缘点进行5点不变量穷举计算解决了模型图像与目标图像的对应点选取问题;通过将不同观测方位和不同旋转角度的样本图像边缘点的5点不变量集合作为输入向量对BP网络进行训练解决了由于仿射和射影变换造成规格化边缘点间距变化对正确分类的影响.算法不仅能识别多边形而且能识别曲线轮廓的目标物体.      

Multipath Separation-Direction of Arrival (MS-DOA) with Genetic Search Algorithm for HF channels

Direction-of-Arrival (DOA) defines the estimation of arrival angles of an electromagnetic wave impinging on a set of sensors. For dispersive and time-varying HF channels, where the propagating wave also suffers from the multipath phenomena, estimation of DOA is a very challenging problem. Multipath Separation-Direction of Arrival (MS-DOA), that is developed to estimate both the arrival angles in elevation and azimuth and the incoming signals at the output of the reference antenna with very high accuracy, proves itself as a strong alternative in DOA estimation for HF channels. In MS-DOA, a linear system of equations is formed using the coefficients of the basis vector for the array output vector, the incoming signal vector and the array manifold. The angles of arrival in elevation and azimuth are obtained as the maximizers of the sum of the magnitude squares of the projection of the signal coefficients on the column space of the array manifold. In this study, alternative Genetic Search Algorithms (GA) for the maximizers of the projection sum are investigated using simulated and experimental ionospheric channel data. It is observed that GA combined with MS-DOA is a powerful alternative in online DOA estimation and can be further developed according to the channel characteristics of a specific HF link.     

Online feedback motion planning for spacecraft obstacle avoidance using positively invariant sets

This paper investigates a novel reliable feedback motion planning algorithm for steering a spacecraft to the desired position while satisfying dynamics and environmental constraints, such as actuator constraints, unexpected state constraints, bounded disturbances, etc. The proposed algorithm, termed safeRRT1, applies positively invariant (PI) set transitions to an online sampling approach based on the asymptotically optimal rapidly-exploring random tree (RRT1). A time-varying safe corridor is generated which consists of a sequence of PI sets under the instantaneous constraints. A spacecraft is guaranteed to remain within the corridor when the feedback controllers corresponding to the PI set transitions are executed. Simulations of the spacecraft avoiding obstacles validate the feasibility and effectiveness of the proposed safeRRT1 algorithm.     

Low-altitude field-aligned currents used to diagnose plasma distribution and magnetic topology in the outer magnetosphere

Field-aligned currents (FAC) can be used to determine changes in the total plasma content (TPC) of convecting flux tubes. The observed steady-state FAC system is combined with the observed equipotential pattern to determine contours of TPC as mapped to the ionosphere. Criteria for a qualitative mapping of the FAC, equipotentials and TCP contours along magnetic field lines to the equatorial plane are set up and the result is shown in figure 3. Some interesting features are: (1) There is a considerable distortion, which is most obvious near midnight, due to the existence of FAC; (2) There is a dusk to dawn component of convection across the tail; (3) The reversal of this component in the pre-midnight quadrant produces the Harang discontinuity. A discussion of time-dependent flows suggests that both plasma depletion associated with FAC and neutral lines may be necessary in a substorm expansion. Between substorm expansions, convection is faster than FAC-produced collapse of flux tubes. Finally it is pointed out that the current and electric field are probably not parallel in the tail, requiring a rethinking of tail models.     

Limit cycles of accretion disks in the Rapid Burster and GRS 1915+105

The Rapid Burster is known to show rapidly repetitive bursts (Type II bursts). An interesting feature of the Type II burst is an approximate proportionality of the burst duration to the time to the next burst. The time sequence from a burst to the following quiescent phase can be said to be a time-scale invariant high (burst phase)–low (quiescent phase) transition. The Galactic superluminal source, GRS 1915+105 exhibits a variety of the time variation of the X-ray flux and often repeats transitions between a high-flux level and a low-flux level. In such high–low transitions, Belloni et al. (1997b) found an approximate proportionality between the duration of the low-flux phase and that of the following high-flux phase, over a wide range of the time scale. This low–high transition can again be said to be time-scale invariant. However, an interesting difference between the two time scale invariant transitions is an opposite order of the high and low-flux phase in the time-scale invariant sequence. In the case of the Rapid Burster, the high-flux (burst) phase is the first, while the low-flux phase is the first in the case of GRS1915+105. A limit cycle between an accretion disk in a state of the standard-disk and that in a state of the ADAF (advection dominated accretion flow) is discussed to explain the time-scale invariant high–low (or low–high) transition as well as the difference between the neutron star system and the black-hole system, qualitatively. A phenomenological relation of the accretion disk change with mass ejections from the central part of the disk is also discussed.