2014年数学建模A题-省一等奖

2014高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)：（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）日期：2014年09月15日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：1嫦娥三号软着陆轨道设计与控制策略摘要月球软着陆是月球探测中的一项关键技术。求解出嫦娥三号的软着陆轨道与最优化控制的模型，让其燃料的消耗尽可能少。对我国的登月计划的发展有重要意义。根据软着陆的一些条件，建立一种基于SQP算法的最优化求解方法来求出软着陆轨道。然后再建立景象匹配的方法，算出6个阶段的最优控制策略，使其消耗燃料尽可能的少。针对问题一，为了求解嫦娥三号在着陆准备阶段的近月点与远月点的位置以及对应的速度大小与方向。我们提出了以运动力学为基础，先求出近月点速skmv/6925.11与远月点的速度skmv/6337.12，再结合软着陆初始条件和齐奥尔科夫斯基公式求出近月点以及远月点的大致位置，并求得方向，然后经过时间的变化范围对所求的位置范围进行优化，得出近月点以及远月点的优化值是kmLkm34.84327.451。针对问题二，我们建立了一种基于SQP方法的常推力月球软着陆轨道非线性优化方法，通过将常推力月球软着陆轨道离散化，利用离散点处状态连续作为约束条件，把常推力月球软着陆轨道优化问题归结为一个非线性规划问题。利用SQP方法求解此轨道优化问题。求出轨道之后，在求解6个阶段的最优控制策略的过程中，采用基于景象匹配的月球探测器精确软着陆轨道的方法，将6个阶段优化为三个阶段。分别为主减速及快速调整阶段、粗，精避障阶段和缓速下降及自由落体阶段。分别建立三个模型，将采样数据与该最佳平面在当地的采样值作差，确定阈值，求出最大粗糙度是mh3.0max。然后再分别求出这三个阶段在优化后所需要的最短时间，因燃料消耗与最短时间之和成正比，故此时6个阶段的最优控制策略可以求解。针对问题三，在问题二中我们对嫦娥三号所建立的导软着陆轨道的研究模型，提出了在景象匹配基础上通过分段控制的方法进行精确软着陆飞行，并满足不同降落段约束条件的精确软着陆控制方案。根据问题二中我们建立的模型，对月球软着陆建立误差模型，运用误差敏感系数矩阵进行分析。分析结果表明，与初始位置偏差相比，初始速度偏差对终端各状态的影响要大；位置、速度测量误差分析只对末端的影响较大。关键词：软着陆、SQP算法、轨道优化、景象匹配2一问题重述1.1问题的背景中国是继美国、前苏联之后的第三个能使卫星登上月球实现软着陆的国家。因此，嫦娥三号如何实现软着陆以及能否成功成为外界关注的焦点。北京时间12月10日晚，嫦娥三号已经成功降轨进入预定的月面着陆准备轨道，这是嫦娥三号“落月”前最后一次轨道调整。在实施软着陆之前，嫦娥三号还将在这条近月点高度约15公里、远月点高度约100公里的椭圆轨道上继续飞行。嫦娥三号着陆地点选在较为平坦的虹湾区。但由于月球地形的不确定性，最终“落月”地点的选择仍存在一定难度。但嫦娥三号的预定着陆点为19.51W，44.12N，海拔为-2641m。在大约距离月球15公里时，反推发动机就要点火工作；到离月球100米时，卫星将暂时处于悬停状态，此时它已不受地球上工程人员的控制，因卫星上携带的着陆器具有很高智能，它会自动选择一块平整的地方降下去，并在离月球表面4米的时候关闭推进器，卫星呈自由落体降落，确保软着陆成功。为了确保探测器能够成功在月球表面实现软着陆，需要认真设计降落过程中探测器的发动机的控制方案，使“嫦娥3号”能够顺利完成科研任务，得到最大化的应用。由于月球上没有大气，嫦娥三号无法依靠降落伞着陆，只能靠变推力发动机，才能完成中途修正、近月制动、动力下降、悬停段等软着陆任务。这将是中国航天器首次在地外天体的软着陆和巡视勘探，同时也是1976年后人类探测器首次的落月探测。嫦娥三号在着陆准备轨道上的运行质量为2.4t，其安装在下部的主减速发动机能够产生1500N到7500N的可调节推力。在给定主减速发动机的推力方向后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。要保证准确地在月球预定区域内实现软着陆，关键问题是着陆轨道与控制策略的设计。其着陆轨道设计的基本要求：着陆准备轨道为近月点15km，远月点100km的椭圆形轨道；着陆轨道为从近月点至着陆点，其软着陆过程共分为6个阶段，要求满足每个阶段在关键点所处的状态；尽量减少软着陆过程的燃料消耗。1.2提出问题根据上述的叙述以及基本要求，提出以下三个问题：（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。3（2）根据附件2，建立相应的数学模型，确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。（3）根据问题二中我们设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。二问题分析1、问题1需要确定的是嫦娥三号着陆轨道近月点和远月点的位置，以及相应的速度大小与方向。首先我们要确定着陆准备轨道的位置，然后再根据题中所给的一些数据，根据天体运动的一些物理公式和霍曼转移来求出近月点和远月点的速度，并确定速度的方向。然后由近月点速度和题中所给的软着陆条件，运用动能定理和齐奥尔科夫斯基公式确定着陆准备轨道近月点和远月点的位置。2、问题二是确定嫦娥三号的软着陆轨道和在6个阶段的最优控制策略。所谓软着陆[1],是指月球着陆器经地月转移到达月球附近后,在制动系统的作用下以很小的速度近乎垂直地降落到月面上。根据软着陆轨道离散化，利用离散点处状态连续作为节点约束条件把确定轨道的优化问题转化为一个多参数优化问题，然后利用逐次二次规划算法（SQP）求解出最优化的着陆轨道。求出轨道之后，确定六个阶段的最优控制策略。确定最优控制策略，需要考虑各个阶段的燃料消耗和每个阶段在关键点所满足的条件，然后在通过景象匹配方法，将六个阶段分为三个阶段，建立模型，然后求解出来最优控制策略。3、在问题二我们确定嫦娥三号的着陆轨道和6个阶段的最优控策略分析的基础上，对模型进行改进和做相应的误差分析与敏感性分析。在问题二中我们对嫦娥三号所建立的软着陆轨道的研究模型。根据问题二中我们建立的模型，对月球软着陆建立误差模型，运用误差敏感系数矩阵进行分析。分析结果表明，与初始位置偏差相比，初始速度偏差对终端各状态的影响要大；位置、速度测量误差分析只对末端的位置、速度影响较大。三问题假设假设嫦娥三号着陆不受月球自转的影响。忽略其他行星对嫦娥三号的引力影响。着陆设备相对于月面的速度接近于零。4假设问题一中软着陆过程力F的大小为恒值。假设问题一中软着陆过程的燃料消耗忽略不计。在观测过程中，一直有阳光的照射。假设在软着陆过程中，嫦娥三号做类平抛运动。四符号说明1v近月点速度2v远月点速度m嫦娥三号在着陆准备轨道上的运行质量M月球的质量G万有引力常量月球的标准重力参数r嫦娥三号至月球中心的距离a嫦娥三号绕月轨道的半长轴x表示着陆点与月心所构成的直线S表示近月点到x的距离L嫦娥三号从近月点到速度为57km的位置所运转的抛物线路程h卫星在主减速阶段运行的垂直高度差ft嫦娥三号从近月点到3000米处所用的时间F发动机的产生的推力五模型建立与求解5.1问题一的模型建立与求解5.1.1确定嫦娥三号相应速度的大小与方向首先，计算嫦娥三号近月点及远月点的速度大小。计算速度大小采用霍曼转移轨道公式。霍曼转移轨道[2]：移轨道上物体的总能等于动能与重力位能的和，而总能又等于重力位能（轨道半径为轨道半长轴a时的重力位能）的一半，由此通过霍曼转移轨道的公式求出近月和远月的速度大小。5霍曼转移轨道公式：aGMmrGMmmv2212化简得：arv12求得近月点和远月点的速度：1v2vskmv/6925.11skmv/6337.125.1.2求解近月点和远月点的位置由动能定理：mghFLmvmv22212121求得L的大致范围：kmLkm33.225627.451上式所求L值为嫦娥三号从近月点到速度为km57的位置所运转的抛物线路程，由分析得，Sh,则有三角形勾股定理可知，s和L的值基本相等。然后，对距离L进行优化。由近月点速度skmv/6925.11，月球半径轨道kmr013.17520，根据公式：001wrv得近月点的角速度：radw401065.9。再通过齐奥尔科夫斯基[3]和终端结束时间（ft）可得优化后的L的范围。齐奥尔科夫斯基公式：FmIIvvtspspff//exp100由齐奥尔科夫斯基公式可转换求得F：fspspftmIIvvF//exp100所以，当力NF1500时，算得stf2006；当力NF7500时，算得stf318.401当stf750时，NF2.4013，通过推力F的两个极值1500N和7500N，可以得出：stsf7503187.401NFN75002.4013故可得准备轨道的近月点的位置：6以着陆点处上方km15处为圆心，L的取值范围：kmLkm34.84327.451由于近月点与远月点是与月球中心位置在一条直线上，故远月点是在与近月点相对于月球球心对称的元月轨道上面。因此，远月点的位置是距着陆点下方月球直径加上远月点的高度处。范围与近月点的范围一样。近月点的位置：以(19.51W,44.12N)为圆心，以kmLkm34.84327.451为圆环曲面内；远月点的位置：以(160.09E,135.48S)为圆心，以kmLkm34.84327.451为圆环曲面内;如图1所示，上下阴影部分分别为近月点和远月点所在的位置范围。图15.1.3误差分析由于我们使用的是物理公式求解的速度大小，并没有考虑月球的自传和其他行星对月球的万有引力的影响，并且月球表面还有重力。求嫦娥三号的速度时，这些因素都没有考虑进去，因此，我们求的速度会有一定的误差。还有，我们在确定近月点和远月点的位置时，直接采用上述求得的近月点的速度来用。所以对我们求的位置范围也会有误差。我们建立的模型还学要进一步完善改进。5.2问题二的模型建立与求解确定嫦娥三号的着陆轨道和在六个阶段的最优控制策略[4]。着陆轨道是从近月点至着陆点，其软着陆过程共分为6个阶段，而每一个阶段都不是独立的，前一个阶段完成的动作，要考虑后几个阶段的技术要求。我们所研究建立的软着陆轨道就是指从霍尔曼转移的近月点开始到着陆点的着陆轨迹。在近月点，嫦娥三号开始准备着陆，发动机提7供减速动力。霍尔曼轨道如图2所示：图2求嫦娥三号的最优着陆轨道，我们采用逐次二次规划算