基于可控域的定点返回轨道全局敏感性分析.pdf

第 卷第 期 年 月系统工程与电子技术 文章编号：（）网址：收稿日期：；修回日期：；网络优先出版日期：。网络优先出版地址：基金项目：国家自然科学基金（）；湖南省自然科学基金（）；载人航天工程科技创新团队项目资助课题通讯作者引用格式：王奇，陆林，李海阳，等基于可控域的定点返回轨道全局敏感性分析系统工程与电子技术，（）：犚犲 犳 犲 狉 犲 狀 犮 犲犳 狅 狉犿犪 狋：，（）：基于可控域的定点返回轨道全局敏感性分析王奇，陆林，李海阳，杨路易，（国防科技大学空天科学学院，湖南 长沙 ；空天任务智能规划与仿真湖南省重点实验室，湖南 长沙 ）摘要：针对载人月球探测定点返回轨道在多变量、多约束的复杂参数空间内参数之间影响关系不确定的问题，对载人探月定点返回轨道的参数敏感性问题进行了研究。基于近月点伪参数，建立了定点返回轨道的可控域计算模型，并提出了一种基于可控域的全局参数敏感性分析方法，以计算得到的可控域集合点构建敏感性分析矩阵，通过基于方差的方法对各参数的敏感性进行求解。在仿真分析中，得到了在给定约束条件下定点返回轨道可控域范围的一般分布规律，并给出了可控域范围内近月点伪参数对终端参数的敏感性指标。相关研究结论可为未来的载人月球探测轨道任务分析与设计提供重要参考。关键词：载人月球探测；定点返回轨道；可控域；全局敏感性分析中图分类号：文献标志码：犇犗犐：犌 犾 狅 犫 犪 犾狊 犲 狀 狊 犻 狋 犻 狏 犻 狋 狔犪 狀 犪 犾 狔 狊 犻 狊狅 犳狆 狅 犻 狀 狋狉 犲 狋 狌 狉 狀狅 狉 犫 犻 狋犫 犪 狊 犲 犱狅 狀犮 狅 狀 狋 狉 狅 犾 犾 犪 犫 犾 犲犱 狅犿犪 犻 狀，（犆狅 犾 犾 犲 犵 犲狅 犳犃犲 狉 狅 狊 狆犪 犮 犲犛 犮 犻 犲 狀 犮 犲犪 狀犱犈狀犵 犻 狀 犲 犲 狉 犻 狀犵，犖犪 狋 犻 狅 狀 犪 犾犝狀 犻 狏 犲 狉 狊 犻 狋 狔狅 犳犇犲 犳 犲 狀 狊 犲犜犲 犮 犺 狀 狅 犾 狅 犵狔，犆犺 犪 狀犵 狊 犺 犪 ，犆犺 犻 狀 犪；犎狌 狀 犪 狀犓犲 狔犔犪 犫 狅 狉 犪 狋 狅 狉 狔狅 犳犐 狀 狋 犲 犾 犾 犻 犵 犲 狀 狋犘 犾 犪 狀 狀 犻 狀犵犪 狀犱犛 犻犿狌 犾 犪 狋 犻 狅 狀犳 狅 狉犃犲 狉 狅 狊 狆犪 犮 犲犕犻 狊 狊 犻 狅 狀 狊，犆犺 犪 狀犵 狊 犺 犪 ，犆犺 犻 狀 犪）犃犫 狊 狋 狉 犪 犮 狋：，犓犲 狔狑狅 狉 犱 狊：；引言月球作为距离地球最近的天体，是人类进行深空探测的重要中转站。近年来，随着航天探测技术的不断进步，月球的应用价值和开发潜力逐步显现，世界各主要航天国家纷纷将下一步的航天目标投向月球，其中包括“阿尔忒弥斯”计划、“嫦娥工程”计划等。新一轮探月热潮将主要目标集中在载人月球探测，并最终进行月球基地的建设。第 期王奇等：基于可控域的定点返回轨道全局敏感性分析 在任务过程中，对于从月球返回地球的航天器，不仅需要确保航天员可以安全返回地球，同时为了最小化回收资源，以便对航天员进行快速搜救、开展航天器回收工作，客观上要求航天器可以准确地返回指定的着陆场。由于定点返回轨道涉及到较多的约束参数，并且各参数之间的影响程度各不相同，这对定点返回轨道的进一步认识和理解带来了困难。因此，针对定点返回轨道的相关参数进行深入的参数敏感性分析对于未来的月球探测任务具有重要意义。参数敏感性分析是一种基于数学统计来研究数学模型的方法，通过对系统数学模型进行研究，可用于确定系统的各个不确定性输入对模型不确定性输出的影响程度 ，广泛运用于水文过程模拟、材料结构分析、风险评估 等领域。运用参数敏感性分析方法可以进一步对定点返回轨道的参数敏感性特性进行分析，确定对返回轨道影响较高的参数，从而为后续的轨道优化设计提供参考。目前，已有诸多学者对月地返回轨道问题进行了研究，文献 基于双二体假设动力学模型，对月地返回轨道特性进行了一定程度的分析。文献 利用非线性规划算法优化求解了不同机动方案的月球逃逸轨道。文献 结合三体问题模型，运用兰伯特算法求解了航天器在环月飞行阶段的返回问题。文献 针对月地返回飞行任务，提出了一种双重优化算法的两层串行求解策略。文献 提出了一种利用单脉冲和三脉冲机动方案构造月地返回轨道的方法。文献 针对返回近地空间站的月地转移轨道优化问题，提出了一种基于改进的多圆锥曲线的高精度轨道设计方法。文献 在考虑满足给定约束条件的情况下，通过一种返回轨迹的构建算法构建了从月球返回地球大气层再入点的单脉冲和三脉冲轨道机动方案。上述研究大多基于一般月地转移轨道，在定点返回轨道的研究方面。文献 采用多起点算法对单脉冲定点返回轨道进行了优化求解，但搜索近乎遍历，效率不高。文献 针对载人月球极地探测任务，对定点返回轨道优化设计问题进行了研究。文献 提出了一种解析方法来快速生成单脉冲定点返回轨道的初始解。文献 针对三脉冲返回轨道方案，提出了一种从初步计算到精确计算的串行求解策略。文献 基于单脉冲月地返回轨道的可达域分析，推导了返回轨道存在性判据，并给出了单脉冲定点返回轨道的多层快速迭代设计方法。上述研究针对定点返回轨道的设计方法和轨道参数相关特性进行了分析，但尚未对定点返回轨道的可控域问题进行深入的研究。定点返回轨道可控域表征了可行解的空间范围，可作为可达域问题 的逆问题进行求解，相关结果可直接用于相关的空间任务分析和设计。对于地月空间的可达域问题，文献 提出了精确可控域数值延拓方法，系统研究了载人登月各阶段轨道的可达域问题。文献 对从地月南向近直线晕轨道出发的载人奔月轨道到达月面的可达域特性进行了分析。文献 基于近月伪参数对一般月地返回轨道的可达域进行了分析，但未进一步分析定点返回轨道的可达域。可以看出，目前关于可达域问题的研究较多，而较少有针对可控域问题进行研究。由于近月点参数的可控域直接决定了定点回轨道的可行解范围，因此直接基于定点返回轨道的可控域进行参数敏感性分析，可以得出在可行的空间范围内各参数之间的影响程度。同时，目前关于区间的全局敏感性分析方法大多基于空间采样进行计算，而定点返回轨道的可控域为在约束条件下的有限区间范围，直接采用空间采样方法会对结果造成较大的误差，因此需要进一步发展一种针对可控域的敏感性分析方法，从而更准确的反映出可控域范围内的参数敏感性特性。基于上述考虑，本文针对给定约束条件下的定点返回轨道，具体分析了近月点参数的可控域分布范围。并在定点返回轨道可控域分析的基础上，以参数可控域集合点为初始计算点，提出了一种基于可控域的参数敏感性分析方法。在仿真分析中，给出了定点返回轨道近月点参数的可控域分布规律，并基于提出的可控域敏感性分析方法，详细分析了在可控域空间范围内各主要参数间的影响程度，给出了各参数间的敏感性大小顺序。本文的相关研究结果可帮助工程人员进一步了解定点返回轨道的相关特性，同时也可为进一步的定点返回轨道的优化设计提供参考。问题描述月地返回轨道一般指航天器从环月轨道施加脉冲进行加速转移，在地球和月球引力的作用下，最终进入地球大气层的转移轨道。对于定点返回轨道，航天器需要以合适的状态到达地球大气层，之后通过大气减速降落到地面并进入到预定的着陆场，本文主要针对航天器从近月出发点到地球再入点的飞行轨道进行研究。在月地返回轨道设计中，选取近月点伪经度犿、近月点伪纬度犿、近月点轨道偏心率犲犿、近月点速度倾角犻犿作为优化设计变量，当给定近月点时刻狋犿和近月点高度犺犿时，可采用圆锥曲线拼接法进行月地返回轨道的设计，得到转移轨道的终端参数。对于定点返回轨道，还需要同时满足工程约束和轨道约束，其中月地返回轨道的终端特征点如图所示，图中为返回轨道升交点，为再入点，为真空近地点，为月球反垂点，为着陆点。图月地返回轨道终端特征点 系统工程与电子技术第 卷轨道约束主要包括出发环月轨道约束、地球终端约束和返回窗口约束，环月轨道约束主要为环月轨道高度约束。犺犿犺（）式中：犺为初始环月轨道高度，通常由工程任务给定。返回窗口约束指在返回时刻，月球反垂点的赤纬犿满足：犞犿犛（）式中：犞和犛分别为真空近地点和着陆点的赤纬；犞可由几何关系求得 （）犻）（）式中：为升交点到再入点地心夹角；为再入点到真空近地点的地心夹角，可通过迭代求解获得 犲犲 犲犲 （）式中：为再入角；犲犲为返回轨道偏心率。地球终端约束需要满足着陆点位置要求：犺犺 烅烄烆（）式中：犺为航天器月地转移轨道真空近地点的高度；、为航天器到达再入点时的经度和纬度。约束条件中的再入点的经纬度、可以由着陆点的经纬度计算，真空近地点高度的约束犺可由给定的大气层高度和再入角求得。此外，不等式约束还包括工程约束中的飞行时间约束，再入升降轨约束，出发升降轨约束。犜犜 犲，；犲，犿，；犿，烅烄烆 （）式中：狉 为再入点地心距；犿、犲分别为出发和到达时轨道的相位角；和 分别表示升段和降段。基于可控域的全局参数敏感性分析方法 可控域计算模型可达域通常指在如式（）的动力学系统中，狋狋，狋犳系统连续，给定狋时刻初始状态狓（狋）狀，存在控制量狌（狋）犝犿，使得狋犳时刻终端状态狓（狋犳）狀，则称集合狀为初始状态集合狀对应的可达域。狓（狋）犳（狋，狓（狋），狌（狋）狔（狋）犮狋，狓（狋）狓（狋）狀，狌（狋）犝犿；狔（狋犳）犢犽犳（）：狀犝犿狀犮（）：狀犢烅烄烆犽（）式中：（）为对时间进行微分；狔（狋）为狋时刻的状态。在有些情况下，终端状态可达域并不直观，也不是所有参数都需要关注，转而求解对应的可达域犢犽。在本文所研究的定点返回轨道，其可控域是在给定返回轨道终端约束的情况下，由初始时刻的近月点轨道参数组成的多维参数集，实际上是上述可达域问题的逆问题。飞行器在狋犳时刻的约束为狀狋犳狓（狋犳）：狋狋犳，犺犺；牔 牔 ；牔犲，牔犛犻烅烄烆 （）则飞行器在狋时刻的参数可控域可以用一般化的数学描述为狀狋狓（狋）：狋狋犳，犺犺；犳，牔（狋犳狋）犜烅烄烆（）在实际中对定点返回轨道进行分析设计时，主要关注近月点的设计参数在初始时刻的分布范围，其对应的可控域可以表示为犢犽狋（犿，犿，犲犿，犻犿）：狓（狋）狀狋（）在可控域计算中，首先给定返回时刻以及约束条件，然后对返回窗口进行判断，若当前时刻满足定点返回轨道窗口约束，则遍历近月点伪参数，通过圆锥曲线拼接法对定点返回轨道进行求解，得到转移轨道的终端参数，然后判断该轨道是否满足当前时刻的约束条件，若满足定点返回轨道要求，则输出近月点参数以及对应的转移轨道参数，若不满足则选择下一组近月点参数进行计算，直至近月点参数遍历完成。计算可控域的流程图如图所示。图可控域计算流程图 全局参数敏感性计算方法对于本文的定点返回轨道问题，可将计算模型在数学上描述为一个函数映射：犢犢（犡，犡，犡犱）犢（犡）（）第 期王奇等：基于可控域的定点返回轨道全局敏感性分析 式中：犡犡，犡，犡犱为用向量表示的模型输入参数；犱为输入参数数目；犢为模型的输出参数，可表示为犢犢，犢，犢犿（）式中：犿为输出参数的数目。在文本研究的定点返回轨道问题中，犡为近月点伪参数，犢则为转移轨道的终端参数。目前，常用的敏感性分析方法通常分为局部敏感性分析方法 和全局敏感性分析方法，对于本文所研究的可控域敏感性问题，局部敏感性只能反映局部点处的敏感性信息，无法在区域范围内对各参数的敏感性进行综合分析。因此，需要使用全局敏感性分析方法对可控域内的参数集合点进行敏感性分析。在基于方差的全局敏感性分析基础上提出了一种蒙特卡罗计算方法，该方法只需要对样本点进行单层抽样，假设采样次数为犖，参数数目为犱。该方法首先在参数空间范围内进行采样生成两个初始分析矩阵犃和犅，然后通过矩阵间的变换总计得到犖犱个计算输入矩阵，再将上述矩阵代入到蒙特卡罗计算公式中，可得到各参数的敏感性结果。但通常计算得到的定点返回轨道的可控域结果为限制在空间某一范围中的结果，若采用上述方法在整个空间范围内进行采样会产生很多非可控域内的参数集合点。以三参数的可控域范围为例，如图（）所示，可以看出在使用采样方法采样时，会在空间立方体中生成一系列采样点。该立方体中除了可控域集合中的状态点，还包含了许多非可控域内的状态点，在敏感性计算中仍使用空间中生成的序列采样点进行计算会对结果产生较大的误差。由于在计算可控域的过程中已经计算出一组参数集合点，该集合点在该区域内为近似均匀分布，同时在计算定点返回轨道中的等式约束一般是限制在允许的误差范围内，在可控域集合点附近加上微小扰动仍存在满足定点返回轨道要求的集合点。因此，可以在原有的可控域集合点施加小扰动变量生成另一组可控域集合点，如图（）所示。这样将这两组可控域集合点作为初始分析矩阵进行全局敏感性分析得到的结果仍然是限制在可控域范围内的敏感性结果，可以较为准确的反映在可行解范围内参数的敏感性。图空间采样点与可控域集合点关系图 基于上述考虑，针对定点返回轨道的可控域问题，结合 的蒙特卡罗方法，提出一种基于可控域的全局敏感性分析方法，该方法的基本计算步骤如下。（）首先通过模型计算得到满足要求的近月点参数可控域集合犃：犃犡（）犡（）犡（）犡（）犡（）犡（）犡（）犡（）犡（犖）犡（犖）犡（犖）犡（犖）熿燀燄燅（）式中：犡，犡，犡，犡犿，犿，犲犿，犻犿表示近月点伪参数；犖为计算得到的可控域集合总数。接着，定义另一矩阵犅，该矩阵中各参数的值为矩阵犃加上各参数对应的比例区间的结果：犅犡（）犽犡（）犽犡（）犽犡（）犽犡（）犽犡（）犽犡（）犽犡（）犽犡（犖）犽犡（犖）犽犡（犖）犽犡（犖）犽熿燀燄燅（）式中：、分别为参数犿、犿、犲犿、犻犿对应的可控域区间范围；犽为比例因子，当给定合适的犽值大小时，可以得到另一组接近原可控域集合点矩阵犅。（）通过矩阵犃和犅构建新的矩阵犃犻犅，该矩阵中第犻列来自矩阵犅，其余犱列来自矩阵犃：犃犻犅犡（）犡（）犻犽犻犡（）犡（）犡（）犻犽犻犡（）犡（犖）犡（犖）犻犽犻犡（犖）熿燀燄燅（）式中：犻，犱，通过上述变换可得到犱个犃犻犅矩阵，总计得到犱个分析样本矩阵。（）将上述犱个样本矩阵代入计算模型中，可得到犱个输出矩阵：系统工程与电子技术第 卷犢犃犢（犃（）犢（犃（）犢（犃（犖）熿燀燄燅），犢犅犢（犅（）犢（）犢（犅（犖）熿燀燄燅），犢犃犻犅犢（犃犻犅（）犢（犃犻犅（）犢（犃犻犅（犖）熿燀燄燅）（）式中：犃（犼）（犼，犖）表示矩阵犃的第犼行，其余各项定义类似。（）选择输出参数犢犾（犾，犿），将上述结果代入到蒙特卡罗计算公式中计算各阶方差：犡犻犡犻（犢狘犡犻）犖犖犼犢（犼）犅（犾）（犢（犼）犃犻犅（犾）犢（犼）犃（犾）（）犡犻犞犡犻（犢狘犡犻）犖犖犼（犢（犼）犃（犾）犢（犼）犃犻犅（犾）（）犡犻 犼犞犡犻犡犼（犢狘犡犻 犼）犖犖狑（犢（狑）犃犻犅（犾）犢（狑）犃犼犅（犾）（）式中：犢犃（犾）表示矩阵犢犃中对应犢犾的输出结果向量；犡犻表示除了参数犡犻之外的其他参数；犡犻 犼表示除了犡犻和犡犼之外的其他参数。总方差为犞（）犖犖犻（犢（犻）犃（犾）（犖犖犻犢（犻）犃（犾）（）然后通过敏感性指标可以计算出各阶敏感性系数：犛犻犞犡犻犡犻（犢犾犡犻）犞（）犛犻 犼犡犻 犼犡犻犡犼（犢犾犡犻 犼）犞（）犛犜 犻犡犻犡犻（犢犾犡犻）犞（）式中：（）表示方差；（）表示期望；犛犻为一阶敏感性指数，表示单一模型参数对模型输出的影响；犛犻 犼为二阶敏感性指数，表示两个模型参数的组合对模型输出的影响；犛犜 犻为全局敏感性指数，表征的是所有包括第犻个模型参数在内的参数组合对于模型输出的影响。上述敏感性系数在计算时基于已有的可控域集合点，不需要进行额外的采样点生成，因此敏感性结果可以较为准确的反映在定点返回轨道可行解范围内各参数之间的敏感性信息。仿真分析 可控域计算在仿真分析中，为统计不同返回时刻下的定点返回轨道可控域分布的一般规律，分别选取近月点时刻为犜：、犜：、犜：，着陆场位置为（，），近月点高度为 ，返回再入角，返回时间不超过 天，返回轨道采用降轨出发，升轨到达方式。在计算中遍历近月点伪经纬度，最终计算得到的可控域结果。为验证本文方法计算得到的可控域的准确性，首先选择犜时刻，然后将本文计算得到的可控域与采用多圆锥截线法 计算得到的可控域进行对比，结果如图所示。图不同计算方法得到的可控域范围 图中，图（）为两个可控域范围的对比图，图（）表示由多圆锥截线法计算得到的可控域范围，图（）表示由圆锥曲线拼接法得到的可控域范围。从两个图的比较分析中可以看出，由两种方法计算得到的可控域范围基本一致，本文的计算方法可较为准确地反映定点返回轨道的可控域范围，因此继续采用本文方法进行后续的计算。进一步比较不同时刻下定点返回轨道的可控域分布范围，结果如图所示。第 期王奇等：基于可控域的定点返回轨道全局敏感性分析 图近月点伪经纬度可控域 从图中结果可以看出，对于不同返回时刻下的定点返回轨道，满足定点返回条件的近月点伪经纬度可控域分布均较为连续，且可控域分布的范围和形状均比较相似，伪经度范围约在 ，伪纬度范围在 。由于不同定点返回轨道间近月点参数可控域的相似性，因此在本文后续的分析中，只选取时刻犜下的定点返回轨道进行分析。图近月点偏心率与速度倾角范围 图所示为犜时刻下近月点偏心率和速度倾角与近月点伪经纬度间的关系。从图可以看出，犜时刻下的定点返回轨道的伪参数近似分布在一连续的区域范围内，具体分布范围如表所示。表中，表示参数区间范围下限，表示参数区间范围上限，为参数区间范围。偏心率和倾角的范围均和伪经纬度的位置有关，偏心率的大小分布在，且偏心率的大小大致随着近月点伪经度的增加而减小。近月点速度倾角的范围分布在 ，与偏心率不同，速度倾角的大小受近月点伪纬度的影响较大，其大小近似随着伪纬度的增加而减小。表近月点参数可控域区间范围犜 犪 犫 犾 犲犐 狀 狋 犲 狉 狏 犪 犾狉 犪 狀 犵 犲狅 犳犮 狅 狀 狋 狉 狅 犾 犾 犪 犫 犾 犲犱 狅犿犪 犻 狀狅 犳狆 犲 狉 犻 犾 狌 狀 犲狆 犪 狉 犪犿犲 狋 犲 狉 狊参数（）（）犲犻（）参数敏感性分析由上述可控域分析已知在给定约束条件下的定点返回轨道伪参数的可控域范围区间，为进一步判断近月点参数对定点返回轨道终端参数具体的影响程度，在此基础上使 系统工程与电子技术第 卷用本文提出的基于可控域的参数敏感性分析方法对定点返回轨道的参数进行敏感性分析。在参数敏感性分析中，输入参数为个伪参数设计变量，输出参数选择工程中的主要关注参数，分别为再入点经度、再入点纬度、返回时间、返回轨道倾角以及真空近地点高度。在进行参数敏感性分析前，首先进行误差收敛性判断以确定犽的取值大小。不同的犽值下，计算得到再入点经度误差限的结果如图所示，本文收敛误差限取值为，图中所示结果即为计算得到的约束误差除以收敛误差限后的结果，当误差限小于即认为满足约束要求。从图中可以看出，当犽值取值较大时，有较多的集合点不满足约束参数的收敛要求，随着犽取值的减小，落在误差限内的集合点数目逐渐增加。进一步统计满足定点返回轨道所有约束的集合点个数，如图所示，当犽 时，满足约束的集合点个数为 （集合点总数为 ），达到了 ，此时可认为矩阵犅中的集合点已基本满足定点返回轨道约束要求。图不同犽值下集合点收敛情况 犽 图符合约束要求的集合点总数随犽值变化情况 犽 进一步分析在不同犽值下参数敏感性分析结果的变化情况，如图所示。图中结果为近月点经度对再入点经度的敏感性结果，从图中可以看出，当犽取值较大时，参数敏感性计算结果不稳定，出现较大的变化幅度。而随着犽取值减小，参数敏感性的计算结果逐渐稳定，当犽取值为 时，参数敏感性计算结果为 ，相比于犽取值为 的计算结果 变化幅度仅为 ，当犽值小于 时，的取值对参数敏感性结果影响很小，因此在本文中取犽 进行后续的分析计算。图不同犽值下的敏感性系数结果 犽 近月点伪参数的可控域范围通过上述计算已知，然后在计算得到的伪参数可控域矩阵犃基础上通过比例因子犽生成矩阵犅，接着通过上述两个矩阵进行参数敏感性分析可得到近月点参数对转移轨道终端参数的敏感性指标。由于在可控域范围内，每个可行状态点的变化往往伴随着多个参数的改变，需要考虑各参数在整个可控域内与各参数相互作用下对结果的综合影响，因此主要统计各参数的全局敏感性系数即犛犜犻，得到的结果如图 所示。第 期王奇等：基于可控域的定点返回轨道全局敏感性分析 图 全局参数敏感性分析结果柱状图 从图 的敏感性分析结果可以看出，通过敏感性分析可以具体地区别出各参数对输出结果的敏感性大小，从而判断不同参数间的影响程度。对于不同的输出参数，近月点参数的敏感性大小也各不相同，具体敏感性结果如表所示。从表中可以看出，对于所选的输出参数来说，近月点偏心率均为最敏感的参数，说明在个伪参数设计变量中偏心率对定点返回轨道的转移终端参数最为敏感，因此在对定点返回轨道进行进一步优化和调整时需要主要考虑偏心率对结果的影响。此外，从表中也可以看出，在输出参数的低敏感性参数中，除返回轨道倾角的低敏感性参数为近月点速度倾角外，其余参数的低敏感性参数均为近月点伪纬度。因此，从整体上看近月点伪纬度对定点返回轨道的终端参数影响较低，在对定点返回轨道进行进一步优化和调整时可降低对近月点伪纬度调整的优先级。表参数敏感性排序犜 犪 犫 犾 犲犛 狅 狉 狋 犻 狀 犵狅 犳狆 犪 狉 犪犿犲 狋 狉 犻 犮狊 犲 狀 狊 犻 狋 犻 狏 犻 狋 狔狉 犲 狊 狌 犾 狋 狊参数敏感性高低再入点纬度犲犿犿犻犿犿再入点经度犲犿犻犿犿犿返回时间犲犿犿犻犿犿返回轨道倾角犲犿犿犿犻犿真空近地点高度犲犿犿犻犿犿进一步分析参数间相互作用下对结果的影响。如图 所示，对再入点纬度影响最大的参数组合为犻犿犲犿，对再入点经度影响最大的参数组合为犻犿犲犿，对返回时间影响最大的参数组合为犿犲犿，对返回轨道倾角影响最大的参数组合为犻犿犲犿，对真空近地点高度影响最大的参数组合为犿犲犿。从结果可以看出，对个转移轨道终端参数影响最大的参数间组合中均包含近月点偏心率犲犿，这也与近月点偏心率始终为高敏感性参数有关。除返回时间外和真空近地点高度外，对结果影响最大的均为近月点速度倾角犻犿近月点偏心率犲犿的组合。因此，在考虑参数间相互作用对结果的影响时，可根据上述的敏感性分析结果进行进一步的确定。系统工程与电子技术第 卷图 参数组合的敏感性结果 结论本文针对月球探测任务中的定点返回轨道，对给定约束条件下的近月点参数可控域进行了分析，并在此基础上提出了一种基于可控域的参数敏感性分析方法。在可控域分析中，给出了不同时刻下的定点返回轨道的近月点参数可控域范围，以此给出了定点返回轨道可控域分布范围的一般规律。在计算得到近月点参数可控域的基础上，具体分析了在定点返回轨道可行解范围内的近月点参数对轨道终端参数的影响程度，给出了参数的全局敏感性指标以及参数间组合的敏感性指标。在实际的轨道设计任务中，可根据求得的参数敏感性结果，对于高敏感性参数需要精确求解范围或者提高权重比，对于低敏感性参数则可缩小其搜索范围或者固定为常量。同时，上述结果为二体条件下的计算结果，因此相关敏感性分析可作为进一步高精度求解的初始判断条件，从而提高求解效率。本文的相关研究结果可为未来的载人月球探测任务提供重要参考。参考文献，：，（）：，（）：，（）：，（），：，（，），（）：，（）：，（），（）：，：，（）：，：，：，（）：，：张磊，于登云，张説月地转移轨道快速设计与特性分析中国空间科学技术，（）：，（）：，（）：，：陆林，杨路易，李海阳，等载人飞船月地返回轨道优化设计与特性分析系统工程与电子技术，（）：，（）：，（）：杨路易，李海阳，张进，等基于改进多圆锥截线的月地返回第 期王奇等：基于可控域的定点返回轨道全局敏感性分析 轨道快速设计方法系统工程与电子技术，（）：，（）：，（）：，（）：陆林，李海阳，刘将辉，等载人月球极地探测定点返回轨道设计中国空间科学技术，（）：，（）：，：，：周晚萌载人探月序列任务有限推力轨道逆动力学设计方法研究长沙：国防科技大学，：，（）：贺波勇载人登月轨道精确可控域数值延拓分析方法长沙：国防科技大学，：，（）：，：，：，（）：，（）：作者简介王奇（），男，硕士研究生，主要研究方向为载人月球探测轨道设计与分析。陆林（），男，博士研究生，主要研究方向为载人月球探测轨道设计研究。李海阳（），男，教授，博士，主要研究方向为载人航天系统分析与仿真。杨路易（），男，博士研究生，主要研究方向为载人月球探测任务智能规划。