2014年本科数学建模优秀论文

2014高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的报名参赛队号为（8位数字组成的编号）：10009062所属学校（请填写完整的全名）：东南大学参赛队员(打印并签名)：1.陶雷2.廖如天3.胡尊丽指导教师或指导教师组负责人(打印并签名)：（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）日期：2014年9月12日赛区评阅编号（由赛区组委会评阅前进行编号）：2014高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：1嫦娥三号软着陆轨道设计与控制策略摘要本文围绕嫦娥三号软着陆过程运行轨道以及控制方案，首先，运用由繁到简思想对问题进行简化理出思路，再由简到繁，应用Pontryagin极大值原理，求解主减速段运行轨迹。建立基于局部燃料最优的最优评价模型，同时以各阶段任务需求为基础，结合实际，考虑燃料最优和实际着陆安全性，对于结果进行一定的修正。利用MATLAB，实现对于问题的求解。对于问题一，在主减速阶段建立了以运行轨迹为抛物线为基础的燃料最优简化模型。首先，根据天体物理学知识，结合嫦娥三号着陆准备轨道为长轴为1837.013km,短轴为1752.013km，运用能量法结合开普勒第二定律，得出近月点和远月点的速度大小。然后，由繁到简，逐步分析简化，从建立立体空间的两个月球着陆极坐标系——轨道坐标系和月球惯性坐标系分析，简化到运行轨迹平面分析模型。将运行轨迹的确定问题转化为轨迹为抛物线的前提下燃料最优问题，列出动力学方程，求出运行轨迹。过程消耗燃料1.0988t,末端速度与径向夹角为97.10，历经428.9333s,水平方向移动距离435.36km。所求得的近月点位置在一个圆上，根据报道椭圆轨道过极点上方，得到近月点位置为（19.51W,29.7595N），高度为1743.013km，速度skmv/6927.11，方向沿轨迹切向向北（如图13），远月点（160.49E,29.7595S)，高度为1837.013km,速度为skmv/6144.12，方向沿轨迹切向向南（如图13）。对于问题二，建立基于局部最优的着陆轨迹和6个阶段控制策略最优评价模型。对于主减速段，根据问题一中的以抛物线为基础的燃料最优简化模型分析思路，利用Pontryagin极大值原理，实行燃料最优制导。耗时438.9s，嫦娥三号质量1301.2kg。对于快速调整段，可采用重力转弯制导，通过对发动机的开关来控制着陆器的斜向速度。耗时7.43s，1296.4。对于粗避障段，通过阈值分析将图像分割以方便计算，选择100×100的区域来计算其方差确定平坦程度，选择距离中心点水平距离较近，且相对平坦的位置作为初步的落月地点，本文中得到的初步落月地点极为图像中心点。此过程嫦娥三号做匀减速运动，实现制导。耗时64.54s，嫦娥三号1221.0kg.对于精避障段，采用中值滤波器对图像进行预处理，然后进行多灰度级阈值分割分析图像，运用障碍圆法求着陆点，选择最近的点（754,643）为圆心作为降落点。同时考虑下一阶段距离以及着陆安全性问题，必须使该阶段末端速度较小。耗时31.11s,嫦娥三号质量1201.0kg.对于缓速下降阶段，该阶段的控制策略为让发动机推力方向向下，保持匀减速运动，直到在4m处速度减小为0，实现悬停。耗时11.56s,嫦娥三号质量1191.5kg.对于自由落体阶段，耗时2.22s,嫦娥三号质量1191.5kg.整个过程的总时间为t=555.76s，燃料消耗量为m=1208.5kg。由于后面几个阶段高度较低，可以直接采用简化的动力学模型。对于问题三，结合第二问中最优轨迹和控制策略的求取过程，对各参数进行误差分析，误差主要包括：数据测量误差、截断误差、模型假设误差。然后在参数标称值的基础上加上+-10%的偏差，通过绘出各个轨迹图像，得出比冲偏差，初始速度偏差，初始质量偏差具有较强的敏感性。在建立模型的过程中，问题的简化是解决问题的一个重点，对于一个复杂的航天问题，简化为我们所熟知的物理数学模型，使得问题理解方便，给予更加清晰的思路，由繁到简，再逐步分析，结合实际情况，由简归繁，就可以解得问题的答案。在模型评价与推广中提出了评价制度更优越的基于遗传算法的多目标着陆轨迹和6个阶段控制策略最优评价模型。在建模的过程中主要运用matlab编程。关键词：燃料最优Pontryagin极大值原理局部最优阈值分析障碍圆法2一、问题重述嫦娥三号于2013年12月2日1时30分成功发射，12月6日抵达月球轨道。嫦娥三号在着陆准备轨道上的运行质量为2.4t，其安装在下部的主减速发动机能够产生1500N到7500N的可调节推力，其比冲（即单位质量的推进剂产生的推力）为2940m/s，可以满足调整速度的控制要求。在四周安装有姿态调整发动机，在给定主减速发动机的推力方向后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。嫦娥三号的预定着陆点为19.51W，44.12N，海拔为-2641m（见附件1）。嫦娥三号在高速飞行的情况下，要保证准确地在月球预定区域内实现软着陆，关键问题是着陆轨道与控制策略的设计。其着陆轨道设计的基本要求：着陆准备轨道为近月点15km，远月点100km的椭圆形轨道；着陆轨道为从近月点至着陆点，其软着陆过程共分为6个阶段（见附件2），要求满足每个阶段在关键点所处的状态；尽量减少软着陆过程的燃料消耗。根据上述的基本要求，请你们建立数学模型解决下面的问题：（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。（3）对于你们设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。附件1：问题的背景与参考资料；附件2：嫦娥三号着陆过程的六个阶段及其状态要求；附件3：距月面2400m处的数字高程图；附件4：距月面100m处的数字高程图。二、问题分析2.1问题一2.1.1软着陆过程分析根据题意，嫦娥三号的着陆点为19.51W，44.12N，海拔为-2641m，因此，可以通过计算嫦娥三号软着陆过程的轨迹，进而推出嫦娥三号软着陆的出发点的运动状态，即近月点的位置和速度。计算软着陆轨迹难点主要集中在距离月球表面15km处到距离月球表面3km处的轨迹方程根据附件一我们知道此段轨迹为抛物线，通过分析其动力方程，结合抛物线的终点位置，进而得出轨迹方程。对于远月点，由于嫦娥三号运行轨道为椭圆轨道，月球处于椭圆轨道的焦点，已知该椭圆轨迹的长轴为100km,短轴为15km,所以椭圆方程可以很快得出，进而得出远月点的位置。根据角动量守恒定理，就可以得到远月点的速度。2.2问题二对于软着陆过程各个阶段的关键点，为了方便描述，我们将其转化为坐标点如下图。图1嫦娥三号软着陆轨迹示意图32.2.1主减速阶段主减速段是软着陆过程用时最长、推进剂消耗最多的任务段。该阶段的主要是减速，消除较大的动力下降段初始水平速度(约1.7km/s)，实现到距离月面3公里处嫦娥三号的速度降到57m/s。因此推进剂消耗优化是该段制导律的主要设计目标，另外还要兼顾自主性和工程可实现性要求。我们通过大量不同模型的筛选，参考张建辉等《月球软着陆轨道优化方法比较研究》]1[等资料，制定嫦娥三号月球软着陆主减速阶段的制导方式。由于彭坤等的《基于混合法的月球软着陆轨迹优化》]2[过程复杂，其他算法需求的测量参数较多，不适合本题，因此利用Pontryagin极大值原理，选择共轭方程的初始值，得到着陆器最优着陆轨迹。2.2.2快速调整阶段为了平缓地从主减速段过渡到粗调整段，有必要快速调整着陆器的姿态和推力，保证粗调整段入口状态的需求。根据主发动机推力和着陆器姿态匀速过渡的要求，并且要兼顾最优性和安全性，采用重力转弯制导。重力转弯软着陆,是指在着陆过程中通过姿控系统使制动加速度方向与速度矢量的反方向始终保持一致。它的主要优点是着陆姿控及测量比较简单。2.2.3粗避障段粗避障段的主要任务是对着陆区成像，实现在设计着陆点上方100m处悬停，并粗步避开大陨石坑，初步确定落月地点。终端相对月面速度接近于0m/s。粗避障段制导必须能够满足制导目标的位置、速度、姿态以及初始高度和速度等多项约束。我们在该阶段考虑如何简化运算来基于粗障碍识别和速度增量综合确定的安全着陆区，避开大障碍选取较为平坦的点，又在原预定着陆点附近，实现综合最优，减少时间消耗以及燃料消耗。2.2.4精避障段精避障段的主要任务是精确避障和下降。区间是距离月面100m到30m。要求嫦娥三号悬停在距离月面100m处，对着陆点附近区域100m范围内拍摄图像，并获得三维数字高程图，得到着陆点。根据悬停段给出的安全着陆点相对位置信息，水平移动到选择的安全着陆点上方，30m处使水平速度为0m/s。精避障段制导律水平方向控制速度和位置,垂直方向控制高度、速度和加速度。2.2.5缓速下降段缓速下降阶段的区间是距离月面30m到4m，主要考虑到着陆安全性，主要控制着陆器在距离月面4m处的速度为0m/s，即实现在距离月面4m处相对月面静止。为了保证着陆月面的速度和姿态控制精度,，尽可能以较小的设定速度匀速垂直下降，消除水平速度和加速度，直到收到关机敏感器信号。考虑到推进剂的消耗和导航位置漂移，选择合适的下降速度。2.3问题三测量误差和系统参数偏差是影响软着陆制导精度的主要原因，在存在测量误差和参数偏差的情况下实现软着陆，满足终端条件是对软着陆控制系统的一个基本要求，实际中测量误差可以通过滤波消除，对闭环系统的影响不大。而系统参数如制动发动机的推力，比冲以及着陆舱的质量都是不可测的，它们是在发射以前在地面标定给出的，但在飞行过程中参数会由于一定的影响而产生一定的偏差考察这些偏差对制导过程的影响就显得十分重要。2.3.1误差分析根据问题二中提出的最优着陆轨迹和控制策略，误差来源可能有三种：一是方程中采用预测参数值或者参数的近似值；二是方案选择过程中为了简化计算，对于一些参数或者条件的忽略；三是计算过程的简化，例如采用简化模型。从中选择对于误差影响较大的因素进行分析。2.3.2敏感性分析可对参数进行一定的变化，画出轨迹，分析着陆轨迹对这些参数的敏感性。三、模型假设1、忽略月球自转2、忽略月球摄动力3、在距月球表面3000m处，基本位于着陆点上方4、忽略月球形状扁率，将月球看做球体5、月球为一个均匀引力场四、符号说明4符号意义符号意义M月球质量F推力m嫦娥三号质量0v动力下降初始速度，即近月点速度G万有引力常数0t动力下降初始时间1v近月点速度fv主减速段终端速度2v远月点速度ft主减速段终端时间ev比冲五、模型建立与求解5.1问题一5.1.1近月点和远月点速度嫦娥三号在着陆准备轨道运行情况如图2，为了更好地观察空间运动情况，忽略实际大小。根据近月点和远月点的机械能守恒和开普勒第二定