[发明专利]行星大气进入着陆器落点不确定度分析方法

说明书

技术领域

本发明属于航天器着陆与返回技术领域，涉及一种行星大气进入着陆器落点不确定度分析方法。

背景技术

在带有大气行星上完成着陆任务，需要在任务前选定预定的着陆点，但探测器在火星大气进入点处的导航控制误差，探测器的气动参数以及火星大气模型的不确定性，都会严重影响着陆器最终的着陆精度，甚至关乎任务的成败。因此，分析这些偏差以及不确定性对着陆点的影响，是一项必不可少的工作；针对带有大气的行星着陆任务，发展一种快速的落点不确定度分析方法，对降低未来火星着陆设计周期和成本，提高设计效率很有意义。目前，在处理这个问题的方法中，总的来说有三类，一是根据系统状态初值及系统方程中不确定参数的统计特性，选择足够多的采样点，进行蒙特卡洛仿真，从而得到各个时刻系统状态的统计特性；二是将系统方程进行线性化，利用线性系统理论对着陆点的统计特性进行分析；三是利用根据系统初始状态的不确定性分布将状态用Askey正交多项式逼近，然后将状态带入到系统动力学中，根据Galerkin投影法则，将表示原系统的随机微分方程转化为一个等效的高维确定性微分方程，最后利用龙格-库塔等数值积分方法，得到各时刻表示系统状态的正交多项式系数，从而得到系统状态的统计特性。

第一类方法需要较高的计算代价，利用这类方法往往需要较长的任务周期，第二类方法虽然计算效率高，但线性化地方法使得在系统初始状态偏差较大时，出现发散现象；第三类方法有完整的数学理论体系，并且计算效率比较高，具有进一步发展的潜力。参见Avinash Prabhakar，James Fisher and Raktim Bhattacharya.Polynomial Chaos-Based Analysis of Probabilistic Ucertainty in Hypersonic Flight Dynamics[J].Journal of Guidance，Control，and Dynamics.2010，33(1)：222-234.中，利用Askey正交多项式和Galerkin投影法将系统表示为等价的高阶微分方程来求解着陆器状态的统计特性，但其没有考虑用Askey正交多项式表示着陆器状态时的截断误差，从而导致在多误差源影响下，算法容易发散的问题。

发明内容

本发明针对现有的行星大气进入着陆器落点不确定性分析技术存在的计算效率低的情况，提出一种行星大气进入着陆器落点不确定度分析方法，能够准确的估计系统状态的统计特性，并且计算效率明显提高。

该行星大气进入着陆器落点不确定度分析方法：

第一步：根据系统初始状态的不确定性分布将状态用Askey正交多项式逼近，构建正交多项式基；

第二步：将系统状态和不确定参数带入到系统动力学中，将表示原系统的随机微分方程转化为一个等价的高阶确定性微分方程；

第三步：利用龙格-库塔等数值积分算法对此高阶确定性微分方程进行积分，求解确定性微分方程，同时对求得的逼近着陆器状态的正交多项式系数进行检测，如果正交多项式的非线性项系数超过预定比例，那么进入第四步，否则进入第五步；

第四步：根据此时的着陆器状态分布特性，利用施密特正交化办法构建新的正交多项式，用新的正交多项式逼近此时的着陆器状态，从新转化成等效确定性微分方程，利用龙格-库塔方法对其进行积分，并监测非线性项系数与线性项系数的比例；

第五步：利用施密特正交化办法建立新的正交多项式，以此类推，直至所需要的停止条件；

第六步：利用数学期望和数学方差的定义，结合每个时刻表示状态的正交多项式，计算此时系统状态的统计特性。

本发明的有益效果：

该发明针对大气进入类行星着陆器落点不确定度问题，能够确保对着陆器统计特性的快速准确估计，并且克服了在多误差源干扰情况下算法发散的问题。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面结合附图对本发明的实施例作详细说明：本实施例在以本发明的技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。

本部分以火星着陆落点偏差的不确定度分析问题为例，给出具体的实施方式。

火星着陆系统动力学为：

(1)