专利名称：行星大气进入着陆器落点不确定度分析方法
技术领域：
本发明属于航天器着陆与返回技术领域，涉及一种行星大气进入着陆器落点不确定度分析方法。
背景技术：
在带有大气行星上完成着陆任务，需要在任务前选定预定的着陆点，但探测器在火星大气进入点处的导航控制误差，探测器的气动参数以及火星大气模型的不确定性，都会严重影响着陆器最终的着陆精度，甚至关乎任务的成败。因此，分析这些偏差以及不确定性对着陆点的影响，是一项必不可少的工作；针对带有大气的行星着陆任务，发展一种快速的落点不确定度分析方法，对降低未来火星着陆设计周期和成本，提高设计效率很有意义。
目前，在处理这个问题的方法中，总的来说有三类，一是根据系统状态初值及系统方程中不确定参数的统计特性，选择足够多的采样点，进行蒙特卡洛仿真，从而得到各个时刻系统状态的统计特性；二是将系统方程进行线性化，利用线性系统理论对着陆点的统计特性进行分析；三是利用根据系统初始状态的不确定性分布将状态用Askey正交多项式逼近，然后将状态带入到系统动力学中，根据Galerkin投影法则，将表示原系统的随机微分方程转化为一个等效的高维确定性微分方程，最后利用龙格-库塔等数值积分方法，得到各时刻表示系统状态的正交多项式系数，从而得到系统状态的统计特性。第一类方法需要较高的计算代价，利用这类方法往往需要较长的任务周期，第二类方法虽然计算效率高，但线性化地方法使得在系统初始状态偏差较大时，出现发散现象；第三类方法有完整的数学理论体系，并且计算效率比较高，具有进一步发展的潜力。参见Avinash Prabhakar, James Fisher and Raktim Bhattacharya. Polynomial Chaos-BasedAnalysis of Probabilistic Ucertainty in Hypersonic Flight Dynamics[J]. Journalof Guidance, Control, and Dynamics. 2010, 33 (I) :222-234.中，利用 Askey 正交多项式和Galerkin投影法将系统表示为等价的高阶微分方程来求解着陆器状态的统计特性，但其没有考虑用Askey正交多项式表示着陆器状态时的截断误差，从而导致在多误差源影响下，算法容易发散的问题。
发明内容
本发明针对现有的行星大气进入着陆器落点不确定性分析技术存在的计算效率低的情况，提出一种行星大气进入着陆器落点不确定度分析方法，能够准确的估计系统状态的统计特性，并且计算效率明显提高。该行星大气进入着陆器落点不确定度分析方法第一步根据系统初始状态的不确定性分布将状态用Askey正交多项式逼近，构建正交多项式基；第二步将系统状态和不确定参数带入到系统动力学中，将表示原系统的随机微分方程转化为一个等价的高阶确定性微分方程；
第三步利用龙格-库塔等数值积分算法对此高阶确定性微分方程进行积分，求解确定性微分方程，同时对求得的逼近着陆器状态的正交多项式系数进行检测，如果正交多项式的非线性项系数超过预定比例，那么进入第四步，否则进入第五步；第四步根据此时的着陆器状态分布特性，利用施密特正交化办法构建新的正交多项式，用新的正交多项式逼近此时的着陆器状态，从新转化成等效确定性微分方程，利用龙格-库塔方法对其进行积分，并监测非线性项系数与线性项系数的比例；第五步利用施密特正交化办法建立新的正交多项式，以此类推，直至所需要的停止条件；第六步利用数学期望和数学方差的定义，结合每个时刻表示状态的正交多项式，计算此时系统状态的统计特性。本发明的有益效果·该发明针对大气进入类行星着陆器落点不确定度问题，能够确保对着陆器统计特性的快速准确估计，并且克服了在多误差源干扰情况下算法发散的问题。
具体实施例方式为使本发明的目的、技术方案和优点更加清楚，下面结合附图
对本发明的实施例作详细说明本实施例在以本发明的技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。本部分以火星着陆落点偏差的不确定度分析问题为例，给出具体的实施方式。火星着陆系统动力学为
秦 V sin γ^ -λρν2 j2Β - μ ηγj^Rm +h)2 (I)务V^Rm + hf+ vcosy/^l^ + h)其中，h表示着陆器距离火星表面的距离，V表示着陆器速度的大小，Y表示航迹角，μ表示火星引力系数，Rm表示火星半径，B表示着陆器的弹道系数，k表示着陆器的升阻比，Φ表示倾侧角，λ表示大气模型不确定性因子，P表示火星大气密度,其与着陆器距离火星表面高度的关系如式(2)所示，它是根据NASA开发的火星大气模型MarsGram所生成的数据进行最小二乘拟合得到的。T = I. 4X l(T13h3-8. 85 X l(T9h2-I. 245Xl(T3h+205. 3645(2)P = 559. 351005946503e_cl·00010511p = P/188. 95110711075T假设系统初始状态及不确定性参数的标称状态及不确定性如下表所示
权利要求
1.行星大气进入着陆器落点不确定度分析方法，其特征在于包括以下步骤 第一步根据系统初始状态的不确定性分布将状态用Askey正交多项式逼近，构建正交多项式基作为第三步检测的基础； 第二步将系统状态和不确定参数带入到系统动力学中，将表示原系统的随机微分方程转化为一个等价的高阶确定性微分方程作为第三步积分的基础； 第三步根据第一步和第二步的结果利用龙格-库塔等数值积分算法对此高阶确定性微分方程进行积分，求解确定性微分方程，同时对求得的逼近着陆器状态的正交多项式系数进行检测，如果正交多项式的非线性项系数超过预定比例，那么进入第四步，否则进入第五步； 第四步根据此时的着陆器状态分布特性，利用施密特正交化办法构建新的正交多项式，用新的正交多项式逼近此时的着陆器状态，从新转化成等效确定性微分方程，利用龙格-库塔方法对其进行积分，并监测非线性项系数与线性项系数的比例； 第五步利用施密特正交化办法建立新的正交多项式，以此类推，直至所需要的停止条件; 第六步利用数学期望和数学方差的定义，结合每个时刻表示状态的正交多项式，计算此时系统状态的统计特性。
2.如权利要求I所述的行星大气进入着陆器落点不确定度分析方法，其特征在于利用Galerkin投影法则将表示原系统的随机微分方程转化为一个等价的高阶确定性微分方程。
全文摘要
本发明属于航天器着陆与返回技术领域，涉及一种行星大气进入着陆器落点不确定度分析方法。首先根据系统初始状态的不确定性分布将状态用Askey正交多项式逼近，然后将状态带入到系统动力学中，根据Galerkin投影法则，将表示原系统的随机微分方程转化为一个等效的高维确定性微分方程，最后利用龙格-库塔等数值积分方法，得到各时刻表示系统状态的正交多项式系数，从而得到系统状态的统计特性，并且在整个过程中根据着陆器状态的统计特性自适应调整正交多项式基底，克服截断误差带来的影响。该发明能够准确的估计系统状态的统计特性，并且计算效率明显提高。
文档编号G06F19/00GK102890743SQ20111020259
公开日2013年1月23日 申请日期2011年7月19日 优先权日2011年7月19日
发明者徐瑞, 崔平远, 朱圣英, 崔祜涛, 任高峰 申请人:北京理工大学