Mapping mean lake surface from satellite altimetry and GPS kinematic surveys

Lake water height is a key variable in water cycle and climate change studies, which is achievable using satellite altimetry constellation. A method based on data processing of altimetry from several satellites has been developed to interpolate mean lake surface (MLS) over a set of 22 big lakes distributed on the Earth. It has been applied on nadir radar altimeters in Low Resolution Mode (LRM: Jason-3, Saral/AltiKa, CryoSat-2) in Synthetic Aperture Radar (SAR) mode (Sentinel-3A), and in SAR interferometric (SARin) mode (CryoSat-2), and on laser altimetry (ICESat). Validation of the method has been performed using a set of kinematic GPS height profiles from 18 field campaigns over the lake Issykkul, by comparison of altimetry’s height at crossover points for the other lakes and using the laser altimetry on ICESat-2 mission. The precision reached ranges from 3 to 7 cm RMS (Root Mean Square) depending on the lakes. Currently, lake water level inferred from satellite altimetry is provided with respect to an ellipsoid. Ellipsoidal heights are converted into orthométric heights using geoid models interpolated along the satellite tracks. These global geoid models were inferred from geodetic satellite missions coupled with absolute and regional anomaly gravity data sets spread over the Earth. However, the spatial resolution of the current geoid models does not allow capturing short wavelength undulations that may reach decimeters in mountaineering regions or for rift lakes (Baikal, Issykkul, Malawi, Tanganika). We interpolate in this work the geoid height anomalies with three recent geoid models, the EGM2008, XGM2016 and EIGEN-6C4d, and compare them with the Mean Surface of 22 lakes calculated using satellite altimetry. Assuming that MLS mimics the local undulations of the geoid, our study shows that over a large set of lakes (in East Africa, Andean mountain and Central Asia), short wavelength undulations of the geoid in poorly sampled areas can be derived using satellite altimetry. The models used in this study present very similar geographical patterns when compared to MLS. The precision of the models largely depends on the location of the lakes and is about 18 cm, in average over the Earth. MLS can serve as a validation dataset for any future geoid model. It will also be useful for validation of the future mission SWOT (Surface Water and Ocean Topography) which will measure and map water heights over the lakes with a high horizontal resolution of 250 by 250 m.     

椭圆轨道卫星对地全球覆盖电推进控制

针对椭圆轨道卫星近/远地点的星下点对全球或特定纬度区域的访问问题，提出一种连续小推力下的对地覆盖控制策略。首先，推导了自然摄动对卫星拱线变化的影响，并探讨了进行小推力覆盖控制的必要性。然后，针对燃料消耗的优化问题，将控制方程展开成含傅里叶级数的形式，用以获得便于星上计算的解析形式的次优解，同时探讨了截取阶数与优化程度的关系。在进行拱线控制的同时，通过合理设置约束，对椭圆轨道的近地点高度进行保护，确保卫星安全运行。仿真结果表明，提出的方法能够以适当的燃料消耗代价实现椭圆轨道的近/远地点的全球覆盖控制或特定纬度区域的反复推扫，且控制力在可接受的范围内。     

Substitution of satellite-based land surface temperature defective data using GSP method

Land surface temperature (LST) as an important environmental variable provides valuable information for earth environmental system modelling. Currently, LST is obtained through satellite thermal sensors at various spatial and temporal resolutions. Although spatially continuous satellite-based LST measurements are intended to overcome the shortcomings of sparse ground-based LST measurements, LST images often contain anomalous values due to the existence of clouds or sensor malfunctioning. The problem becomes more serious where the users deal with high spatial resolution characterized by low temporal resolution. This study examines the capability of a newly developed graph signal processing (GSP) method using two-dimensional single-date thermal data. For this purpose, four Landsat/TIRS datasets are analyzed. The data of five elliptical regions on thermal images are eliminated and then reconstructed through the GSP method and using the LST values of the enclosing rectangles containing the ellipsoids. The results indicate that the temperature variation determined by the GSP method generally conforms to the original image LST values. According to a correlation test conducted on the original image LST and those obtained through the GSP method, the values vary from 58% to 95%, which is an above-the-average rate (RMSE from 0.69 to 2.27). The statistical analysis of the original image LST in both the elliptical regions and the enclosing rectangles containing the ellipsoids indicates that an increase in the variance of LST data causes an increased error in the calculation of temperature by the GSP method, and vice versa. The results of the analysis of variance (ANOVA) and Duncan test indicated that an increase in the number of the non-zero spectral bins would result in increased RMSE values for all the dates and the regions. Moreover, the model errors were significant at the 0.05 level across all the image date and five elliptical study regions. Based on the results, the use of this method is recommended for the reconstruction of LST missing values, where dissimilarity of atmospheric conditions limits the use of other methods that depend on the time series data of various dates and a great deal of data calculation.     

利用微波探测仪(ATMS)对在轨微波辐射计观测精度的模拟分析

微波辐射计的观测精度及其对数值模式同化应用的影响评估是微波辐射计观测指标设计的重要参考。基于微波探测仪(ATMS)资料,利用三维变分同化系统模拟分析在轨微波辐射计的观测精度指标。针对ATMS观测误差特征,在其观测基础上增加均值为零、标准偏差分别为0.5,1.0,1.5,2.0K的正态随机扰动,进而获得不同精度的观测模拟值序列,然后利用Harris和Kelley的辐射资料偏差订正经验方法订正不同精度的观测资料。偏差订正后,利用三维变分同化模式(WRFDA)直接同化ATMS资料。通过2016年6月6h预报场的同化试验,评估了不同观测精度的模拟资料对数值模式的同化影响。     

上升段飞行器目标的视频图像跟踪

粒子滤波是一种基于贝叶斯估计理论和蒙特卡罗理论的实时目标跟踪方法，具有较为灵活的并行化跟踪方式，能够较好地维持跟踪目标的假设状态，具有较好的跟踪效果和鲁棒性。上升段飞行器目标飞行视频图像跟踪是火箭等目标飞行监控的重要阶段，但现阶段对飞行器上升段的视频图像跟踪主要依靠人工手动操作云台控制器，实现视频图像中的飞行器跟踪，跟踪图像存在跟踪滞后、画面抖动等现象，跟踪效果受人为因素影响较大。本文提出一种基于粒子滤波方法的上升段飞行器目标视频图像跟踪方法，建立飞行器目标粒子滤波跟踪模型实现对飞行器目标的识别和跟踪，在识别和跟踪的基础上建立云台控制模型，通过对云台的智能控制获得飞行器上升段的高质量图像。采用火箭发射的视频图像作为模型验证的实验数据，检验飞行器目标的跟踪效果。     

类IXV飞行器初期再入制导与姿态控制方法研究

摘要： 针对类IXV飞行器无翼式升力体构型及采用RCS/尾襟翼组合控制特点，研究其初期再入段的高精度制导与姿态控制方法.设计带过程约束的数值预测 校正制导律以提高制导系统的鲁棒性及精确性；根据其执行机构配置特点，设计基于RCS的航向控制律及基于RCS/气动舵的纵、横向复合控制律，并采用鲁棒伺服LQR技术进行控制参数快速设计.通过蒙特卡洛打靶仿真来验证算法的精度及鲁棒性.     

涡轮基组合循环发动机分布式控制系统通信网络混合拓扑结构优化方法

涡轮机组合循环(Turbine based combined cycle,TBCC)发动机控制系统通信网络拓扑结构是其分布式控制系统方案设计的重要部分,优化网络拓扑结构可提高发动机推重比和控制系统可靠性。本文基于智能优化算法提出TBCC分布式控制系统网络拓扑结构优化方法。基于图论建立TBCC几何模型和网格模型,以重量和可靠性为优化性能指标,同时考虑发动机表面高温区域以及控制节点的工作可靠性,分别采用粒子群算法和遗传算法优化星形结构中智能中央节点位置、中央节点的环形拓扑结构,获得星形-环形混合拓扑结构。仿真实例表明,基于本文方法优化所得的混合拓扑结构相较于星形集中式控制结构,系统重量降低了51.9%。     

多臂空间机器人的视觉伺服与协调控制

针对多臂空间机器人自主目标抓捕任务，首先建立多臂空间机器人的运动模型和其与目标的相对运动模型，采用Kane方法建立多臂空间机器人的动力学模型；其次，研究基于视觉伺服的机械臂在线轨迹规划算法，并引入零反作用机动，消除机械臂运动对平台姿态的扰动；再次，在不使用零反作用机动功能时，分别使用基于角动量前馈补偿的协调控制算法和逆动力学方法设计了协调控制器，在机械臂运动时保持平台姿态和相对目标的位置。最后，开发了基于Matlab的仿真软件MASS（多臂空间机器人仿真），仿真结果校验了上述方法的有效性。     

一种飞机大仰角飞行的模糊控制方法

提出了一种飞机大仰角飞行的动态逆模糊集成控制方法。首先使用动态逆方法求得系统的逆，使系统反馈线性化，然后，用具有在线规则自调整的模糊控制方法来确保系统的最终特性。该方法应用在高性能收音机俯仰机动指令系统的设计中，仿真结果表明，所设计的控制器满足操纵品质的要求，并且，在大仰角飞行时具有较好的鲁棒性。     

多臂协调操作自由飞行空间机器人的位置和内力混合控制

基于多臂自由飞行空间机器人多臂协调操使负载沿着期望的轨迹运动并且控制负载所受的内力为目的，在多臂自由飞行空间机器人系统协调操作的动力学方程基础上，推导了各个机械臂协调操作负载按期望轨迹运动时各个机械臂关节的驱动力矩的计算方法；给出了作用在负载的内力的定义，根据关节力矩计算方法和PID反馈控制原理，建立了多臂自由飞行空间机器人协调操作负载时的位置和内力的控制算法；讨论了所提出的控制算法的稳定性问题，得到了负载的位置误差和内力误差的约束条件。通过仿真实验证明该控制算法能够使负载的实际位姿和内力收敛到期望的轨迹和内力。