一类禁飞区后方安全撤离轨迹的设计方法研究

对交会对接过程中目标航天器后方的一类安全撤离轨迹进行研究,针对长方形禁飞区提出两类撤离模式,并分析与禁飞区最可能相交的点的特征,通过两个定理给出这两类撤离模式自由漂移轨迹安全的充分必要条件．最后针对实际问题提出撤离方案,给出两组仿真算例,仿真结果验证了撤离方法的正确性和撤离方案的可行性．     

参数不确定空间机械臂系统的变结构鲁棒控制

本文讨论载体位置与姿态均无控的空间机械臂的控制问题．基于增广静态反馈控制模型,提出当载荷参数存在不确定性时空间机械臂追踪惯性空间期望轨迹的变结构鲁棒控制方法．通过仿真运算,证实了方法的有效性．     

月面飞行器动力下降段多目标轨迹规划

针对未来月面着陆动力下降段轨迹规划需综合考虑多性能指标的问题,提出一种对飞行轨迹先优化后决策的多目标轨迹规划方法.在多目标进化算法MOEA/D-AWA(multi-objective evolutionary algorithm based on decomposition with adaptive weight adjustment)的框架下对轨迹规划的多个指标进行分解,得到若干个单指标的子问题.将凸优化算法作为求解单目标轨迹优化子问题的底层算法,嵌套在MOEA/D-AWA的框架中,经过迭代优化获得一组动力下降段飞行轨迹,其构成多目标轨迹规划问题的帕累托最优解集.根据模糊决策理论对各个帕累托最优解对应的多个轨迹指标逐步降阶并进行综合评估,经过决策得到多指标约束下的飞行轨迹.仿真实验表明,该轨迹规划方法能够在综合多目标的情况下,优化获得一组动力下降轨迹集合,且能够根据不同任务要求从中决策出最优的动力下降段轨迹,可有效解决月面飞行器的多目标轨迹规划问题.     

一种月球表面飞跃转移轨迹设计方法

针对飞跃器在月球表面飞跃转移轨迹设计问题,提出了基于凸优化方法的整个飞跃过程燃料最优轨迹设计方法.与经典凸优化方法对轨迹分段求解后再拼接得到全轨迹设计的方法不同,在假设垂直上升、着陆时间固定条件下,根据实际工程需要对轨迹进行了分段设计约束,利用黄金分割法搜索上升着陆时间,通过将原问题转化为求解一个二阶锥问题得到了全飞行...     

基于伪谱法的小天体最优下降轨迹优化方法

针对小天体着陆探测下降阶段的多约束轨迹优化问题,基于Gauss伪谱法进行了燃耗最优下降轨迹优化设计,得出了燃耗最优下降轨迹。建立了小天体下降轨迹优化问题的最优控制问题模型,采用Gauss伪谱法进行离散化,转化为非线性规划问题进行了求解。数学仿真结果显示:优化结果符合各项约束条件,以零速度到达了目标着陆点,且符合燃耗最优的优化目标。利用Gauss伪谱法进行小天体最优下降轨迹优化,计算速度快,求解精度高。     

基于高斯伪谱法的空天飞机上升段最优轨迹设计

针对复杂多约束条件下空天飞机上升段燃料最优轨迹优化问题,提出一种基于高斯伪谱法的上升段轨迹优化策略.依据发动机的推力特性将上升轨迹合理分段,使原最优控制问题转化为多段最优控制问题后,采用高斯伪谱法进行并行优化计算.数值仿真结果表明采用这种轨迹优化策略能够满足组合动力系统工作模态转换时对飞行状态的约束条件,可以在较短的时间内完成高精度的上升段轨迹优化任务,从而验证了该方法的有效性.     

运载火箭推力故障下轨迹在线重规划方法

针对运载火箭非入轨飞行段出现动力系统故障的问题,提出了一种基于能量最优的在线轨迹重规划方法。通过参考轨迹进行引力项的凸化以及根据故障时刻和故障程度的不同设置不同的飞行时长,同时采用无损凸化、线性化和梯形离散等方法,将原始最优控制问题转化为可以采用原对偶内点法迭代求解的二阶锥规划问题,重新规划其飞行轨迹。仿真结果表明该算法能有效提高运载火箭应对突发故障的能力,实现轨迹在线快速收敛,满足动力系统故障后飞行轨迹重构需求。     

考虑随机干扰的高超声速滑翔飞行器轨迹优化

为解决高超声速滑翔式飞行器具有干扰不确定性的再入轨迹优化问题,提出了一种基于广义正交多项式方法的随机轨迹优化数值求解方法。针对随机干扰,基于广义正交多项式方式对其分布进行采样,形成采样空间,而后将每个采样值代入确定性轨迹优化问题中进行反复迭代求解,得出观测值的样本空间,并计算其期望、方差和协方差,估计输出值。以最大纵程为代价函数对具有随机干扰的高超声速滑翔式飞行器最优轨迹进行了数字仿真。仿真结果表明,基于广义正交多项式的随机轨迹优化方法能够有效处理随机干扰对轨迹优化问题的影响,与蒙特卡洛法相比计算效率大大提高。     

综合热力学模型的平流层飞艇上升轨迹优化

针对平流层飞艇的上升轨迹优化问题,综合热力学模型进行了研究.主要分析了飞艇基本热力学行为,研究了蒙皮及内部气体的能量方程并建立了详细的飞艇动力学和运动学模型.在热力学、动力学和运动学分析的基础上,建立了以飞行时间为优化目标的平流层飞艇的轨迹优化模型.利用直接配点法将轨迹优化问题转化为非线性优化问题,再通过非线性求解器SNOPT(Sparse Nonlinear Optimizer)对不同场景的问题进行最优化轨迹求解.优化结果表明:热力学效应对优化轨迹有较大影响,在上升过程中,太阳能辐射为主要影响因素,另外风场也对换热量有一定影响.     

变构型多操纵面RLV进场着陆轨迹优化设计

针对变构型多操纵面可重复使用飞行器(RLV)的进场着陆问题,提出了一种进场着陆轨迹设计方法.将进场着陆轨迹分为深下滑着陆轨迹和拉平着陆轨迹;考虑了气动舵面可操纵偏转限制及深下滑拟平衡约束条件,对深下滑着陆轨迹进行了优化设计;考虑了RLV气动舵面调节余量和构型变化过程,以起落架放下时间、拉平法向过载及接地状态为约束,基于RLV动力学方程,采用轨迹推演的方法对拉平着陆轨迹进行了优化设计.深下滑着陆轨迹和拉平着陆轨迹组成完整的进场着陆轨迹,所设计的着陆轨迹综合考虑了多种约束条件,提高了着陆过程平稳性和安全性.      

First estimations of the aerosol optical thickness using Langley Method at Southern Brazil (29.4°S, 53.8°W)

The first estimations of the aerosol optical thickness (AOT) using Langley Method at Southern Space Observatory (SSO) at Southern Brazil (29.4°S, 53.8°W) are presented. In addition to ozone and sulphur dioxide columns, AOT can be obtained using Brewer Spectrophotometer at specific wavelengths: 306.3, 310.1, 313.5, 316.8 and 320.1 nm. The AOT was obtained for the period from November/2002 to May/2003. Very low AOT averages were obtained, whose values were about 0.21 ± 0.03 at 306.3 nm, 0.21 ± 0.02 at 310.1 nm, 0.19 ± 0.02 at 313.5 nm, 0.20 ± 0.02 at 316.8 nm and 0.20 ± 0.02 at 320.0 nm for all period analysed. Different behaviour of AOT were obtained at two daily specific periods of aerosol accumulation, one in the afternoons from November/2002 to February/2003, caused mainly by a mild biomass burning season’s in the region and other in the mornings from March to May/2003, due the high relative humidity presented in the region studied.     

Simultaneous determination of aerosol- and surface characteristics from top-of-atmosphere reflectance using MERIS on board of ENVISAT

The determination of aerosol optical thickness (AOT) from nadir scanning multi-spectral radiometers, like SeaWiFS, MERIS or MODIS, requires the separation of spectral atmospheric and surface properties. Since SeaWiFS and MERIS do not provide information at 2.1 μm, like MODIS, the estimation of the surface reflectance cannot be made by the cross correlation approach described by Kaufman et al., 1997. The BAER approach (Bremen AErosol Retrieval), von Hoyningen-Huene et al., 2003, uses a linear mixing model of spectra for ‘green vegetation’ and ‘bare soil’, tuned by the NDVI, determining an apparent surface to enable this separation of aerosol and surface properties from VIS and NIR channels. Thus AOT can be derived over a wide range of land surfaces for wavelengths <0.67 μm. Using MERIS L1 data over Europe, the AOT retrieved is comparable with ground-based observations, provided by AERONET. Regional variation of AOT can be observed, showing the atmospheric variability for clear sky conditions by: large scale variation of aerosol turbidity, regional pollution, urban regions, effects of contrails and cases of aerosol-cloud interaction. Simultaneously with the spectral AOT also spectral surface reflectance is obtained, where all atmospheric influences have been considered (molecules, aerosols and absorbing gases (O₃)) for channels with wavelengths <0.67 μm. The AOT is extrapolated by Angström power law to NIR channels and the atmospheric correction for land surface properties is performed, enabling the further investigation of land use and spectral land properties.     

Different mechanisms controlling interannual phytoplankton variation in the South China Sea and the western North Pacific subtropical gyre: A satellite study

A time series of remotely-sensed chlorophyll a (chl a) in 1997–2010 was evaluated to determine mechanisms of phytoplankton variation in recent decade in the South China Sea (SCS) and the western North Pacific subtropical gyre (WNPSG). Satellite-derived sea surface temperature (SST) and aerosol optical thickness (AOT) were used as proxies for vertical nutrient supply and atmospheric aerosol, respectively. Chl a in the WNPSG was not significantly correlated with SST (r = 0.18, p > 0.05), but was with AOT (r = 0.31, p < 0.05), indicating the chl a was influenced by atmospheric deposition. Chl a in the SCS was negatively correlated with SST (r = −0.60, p < 0.05) and was positively with AOT (r = 0.20, p < 0.05). The correlation between AOT and chl a in the SCS does not reflect a major contribution from atmospheric deposition to chl a; instead, the relationship resulted from concurrence of the peaks of AOT and wind speed, which drive water mixing and nutrient supply. Consequently, chl a in the SCS would be regulated primarily by the nutrient supply from deep waters. Because SST was controlled by the ENSO teleconnection in the SCS, the chl a was coupled with ENSO events. The present study demonstrated that interannual phytoplankton variation could be controlled by different factors even in neighboring oligotrophic regions.     

Study of atmospheric air pollutants during the partial solar eclipse on 15 January 2010 over Udaipur: A semi-arid location in Western India

The paper describes behavior of surface ozone, its precursor gases, BC along with TOCC, TWVC, AOT1020 nm as well as UV and IR radiation intensities observed during the partial solar eclipse of 15th January, 2010 over Udaipur, where 52% solar disc is obscured due to the moon’s shadow. During the beginning to main eclipse phase, the deviation values of several air pollutants concentrations from eclipse to control day values vary in a small range from −9 to −2 ppb in case of surface ozone and −180 to −80 ppb for CO. The corresponding change in the values of BC observed from −3.3 to −.5 μg/m³. No significant change is found in NO₂, NO or in ratio of NO₂/NO values during the partial eclipse time. TOCC values decrease from 3 to 5 DU along with a reduction in UV radiation intensity from 20 to 35% from starting to the main eclipse phase. The AOT1020 nm values are found to increase from .2 to 1.0 along with a reduction in IR radiation intensity order of 50%. However, TWVC values decrease from .22 to .1 cm during the eclipse hours. The low level of dilution in surface ozone in eclipse period may be attributed with change in local atmospheric boundary layer dynamic conditions or limited air pollutants dispersion, in term of decreases in planetary boundary layer height, wind speed and hence ventilation coefficient in the same eclipse hours. Thus, present studies support the argument for the leading roles of photochemical reactions with its precursor gases under presence of solar radiation in surface ozone variability. Other possible controlling factors are advection of air pollutants from the polluted region as evident from backward wind trajectories and altering the local meteorological conditions.     

委托实现的决策问题研究

应用代理理论的分析手段,对委托实现的决策问题进行了分析.首先对该问题的背景作了与代理理论相一致的假设,然后在此基础上建立了问题的一般模型;接着,对这个一般模型进行了分解,并据此对委托实现决策问题的实质进行了剖析;然后,对分解出的纯实现问题按单个代理人作了进一步的简化,分析了简化实施的条件和简化后解的变化情况.最后,指出了本文的不足和进一步研究的方向.     

Trajectory optimization for asteroid landing with two-phase free final time

This work develops an autonomous trajectory planning algorithm for 6-DOF asteroid landing. The trajectory planning problem is formulated as a nonconvex time-optimal optimization problem with two-phase free final time, while the cost is regularized by augmenting a fuel consumption penalty. The nonconvex optimization problem is solved in successive solution method, and successive convexification is used to convert the original nonconvex problem into a sequence of convex subproblems, where each subproblem is obtained by linearizing the nonconvex dynamics and state constraints and using the velocity increment to give a convex expression of the fuel consumption penalty in cost function. Specifically, in the linearization, we divide the flight time interval into two parts and normalize each part using a time dilation coefficient to solve the problem that both the final times for the two flight phases are unknown, so that the original free final time problem turns to a fixed-time problem by minimizing the sum of the two time dilation coefficients and fuel consumption penalty. Besides, trust regions and virtual control are used to increase robustness of the algorithm. A convergence analysis is presented which indicates the successive solution will recover the local optimality of the original problem. Then the validity of the proposed algorithm and effects of different factors on flight time and fuel consumption are examined by simulations of landing on an irregular asteroid.     

并行加工系统一类排序问题的最优解

讨论了将多个零件分派给并行加工系统加工的排序问题.假设同一零件分配给不同的设备,其加工时间不同,分析了使所有零件的总花费时间(加工时间与等待时间之和)最小的排序方法.首先建立了该类问题的数学模型,然后将其转化为指派问题,通过匈牙利算法可以得到最优解.所得算法的时间复杂性是多项式界的.最后给出了一个数值例子说明求解过程.     

限制性三体问题共线平动点轨道近似解析解

随着深空探测成为航天领域的研究热点,与其密切相关的三体问题基础研究也日益重要,尤其是在深空探测任务设计中处于基础地位的共线平动点附近运动的研究,更是具有重要的工程应用价值。在圆型限制性三体问题下,对共线平动点附近运动近似解析解的研究已经较为全面,但在更接近真实情况、更具一般性的椭圆型限制性三体问题下,相应的研究却相对较少。针对此背景,参考借鉴圆型限制性三体问题的研究方法,首先根据平动点的特性计算出平动点的位置,然后将非线性三体动力学模型在共线平动点处线性化,最后结合线性系统理论,获得了椭圆型限制性三体问题下共线平动点附近运动的近似解析解,并将其与经典的圆型限制性三体问题下的近似解析解进行对比分析,仿真结果证明了方法的有效性,同时也表明所推导的椭圆型限制性三体问题解析解相比圆型限制性三体问题解析解具有更高的精度。     

空间交会中多圈Lambert变轨算法研究

用两个位置和飞行时间来决定轨道的问题 ,即Gauss问题 ,是航天动力学中的一个基本问题。而Lambert交会问题是Gauss问题在航天器空间交会领域的具体化。文中提出了一种新的算法 ,解决了航天器多圈Lambert变轨的求解问题 ,由此来寻求在飞行时间较长的情况下 ,航天器的Lambert交会燃料最优轨道 ;并由算例验证了这种算法的正确性和鲁棒性     

A faster optimization method based on support vector regression for aerodynamic problems

In this paper, a new strategy for optimal design of complex aerodynamic configuration with a reasonable low computational effort is proposed. In order to solve the formulated aerodynamic optimization problem with heavy computation complexity, two steps are taken: (1) a sequential approximation method based on support vector regression (SVR) and hybrid cross validation strategy, is proposed to predict aerodynamic coefficients, and thus approximates the objective function and constraint conditions of the originally formulated optimization problem with given limited sample points; (2) a sequential optimization algorithm is proposed to ensure the obtained optimal solution by solving the approximation optimization problem in step (1) is very close to the optimal solution of the originally formulated optimization problem. In the end, we adopt a complex aerodynamic design problem, that is optimal aerodynamic design of a flight vehicle with grid fins, to demonstrate our proposed optimization methods, and numerical results show that better results can be obtained with a significantly lower computational effort than using classical optimization techniques.