2018年数学建模大赛暨全国选拔赛论文.DOC

1、12018 年数学建模大赛暨全国选拔赛论文论文题目： 小组成员信息： 姓名 学 号 专业班级 联系方式 选题论文总体情况：假设的合理性：合理 比较合理 一般 不太合理建模的恰当性：恰当 比较恰当 一般 不恰当计算的准确性：准确 比较准确 一般 不准确表述的清晰度：清晰 比较清晰 一般 不太清晰综合评分： 大连交通大学理学院数学建模竞赛指导组2018 年 5 月2数学建模竞赛论文撰写模版论文题目：太阳影子定位分析模型摘要月球是距离地球最近的天体，对月球资源和环境进行科学研究和考察，是人类走出地球，探索未知世界所必需经历的重要步骤。为了求解近月点、远月点位置以及探测器位于这两点的速度。通过建立月心

2、-赤道坐标系，从预定着陆点的经纬度出发，导出近月点坐标，进一步确定远月点坐标。嫦娥三号软着陆过程中从着陆准备轨道进入至着陆轨道后，为了减少有限推力作用下月球探测器软着陆所需的燃料消耗,通过软着陆动力学方程，求解月球软着陆最优控制中的两点边值问题，并采用数值解法对其结果进行进一步的优化。-提示：摘要非常重要，请认真撰写 包括主要思想 方法和结论 不要有错别字和重复的内容 关键词：软着陆，非线性规划，数值迭代，曲面拟合，顺时针螺旋搜索3一、问题重述1.1 问题背景嫦娥三号于 2013 年 12 月 2 日 1 时 30 分成功发射， 12 月 6 日抵达月球轨道。嫦娥三号在着陆准备轨道上的运行质量

3、为 2。4t，其安装在下部的主减速发动机能够产生1500N 到 7500N 的可调节推力，其比冲（即单位质量的推进剂产生的推力）为 2940m/s，可以满足调整速度的控制要求。在四周安装有姿态调整发动机，在给定主减速发动机的推力-1.2 问题分析登月探测器首先被发射到环月停泊轨道上然后根据预先选择的着陆地点，着陆器绕环月轨道上经霍曼变轨后转入远月点 100km,近月点 15km 的椭圆形轨道，在该轨道的近月点主减速发动机点火实现所谓的远着陆，着陆器的大部分燃料消耗在主减速发动机的制动发动上，所以月球软着陆的设计主要以燃烧最优性为出发点，同时保证准确着陆到预定着陆点，通过将球形经纬曲面展开至经纬

4、平面的方法，将三维坐标中求近月点，远月-。二、模型的假设1.假设将推力 F 视为常量。2.假设飞行探测器降落过程中，将月球的重力加速度 月g视为常量。3.假设不考虑切向加速度对加速度的影响。4.将月球等效看成圆球体。5.在确定近月点、远月点的位置过程中，不考虑飞行探测器的水平位移。三、符号说明符号 含义 单位 备注v着陆器沿着 r方向上的速度M/sr月心距 Km推力与切线方向的夹角Rad飞行探测器的角速度Rad/sm飞行探测器的质量KG4四、模型建立与求解4.1 软着陆的动力学方程在着陆器减速下降的过程中，月球引力的非球项，日月引力摄动和月球旋转均可忽略，并假设着陆轨道在纵向平面内建立,其图如

5、下所示:图 6.1 月球软着陆极地坐标图轨道平面极坐标系，坐标圆点在月心 O，y 轴指向初始轨道近月点，x 轴指向登月器开始运动的方向。软着陆的动力学方程描述如下：代就能得到一个比较将确定结果，同样也减少了运算量。月球软着陆轨道是一个服从两点边值约束问题，由于起点处于或满转移轨道的近月的竖直投影的集合，可得圆的方程为： 22-lbYaX）（）（， （10）可得直线集合方程为:）（K- (11)并将两式联立可得: baXYSX),12（ 其中 a 和 b 分别为在此平面直角坐标系下着陆点的 X,Y 坐标值分别为 ：018rb (12)与 则分别为着陆点的经度和纬度。即 为 19.51W， 为 4

6、4.12N 代入公式得出即a=591.47km,b=1337.57km。并假设着陆器的运行轨道与月球经度在同一个平面上，那么，lbYaX,，即 为 0ra, 为 0rl，可得近月点的坐标为：5)98.2,51.(), NWYX当考虑远月点时，其 与 所符合的点于月球中心与着陆点两点所连直线上，则= orKS/1802E, = 02 /18)57.3)4.591( rKS（ S,及坐标为)9.,5.(W而关于近月点，远月点相应的速度。采用万有引力定律， )()(22hRvmMGF。在远月点与近月点分别代入 1h=100km, 2=15km。可得:smV/36.1远 月， smV/5.672近 月

7、方向与运行轨道切线方向一致。4.2 自主避障制导设计的精障碍识别和安全着陆区选取算法为: 1)数据预处理, 主要包括对着陆器姿态和平动速度的补偿以及将每个“脚印”的斜距信息转换成垂直距离,构建测量坐标系下据月球表面分别为2.4km和0.1km的地形三维高程图，如下图所示:图6.2 距2400m处的三维数字高程图图6.3 距月面100m处的数字高程图2)平均坡面构建, 采用最小二乘法拟合一定单元区域的平均坡面,最小二乘拟合多项式 )(xps，具有多项式形式 -6五、模型的评价关于软着陆问题的关键是找到最优飞行轨迹和推力大小与方向的时间历程。软着陆转移轨道为 100 km15 km 的椭圆轨道,从

8、近月点到月面为软着陆全过程。假设月球引力场均匀,忽略月球自转,建立的着陆坐标 5系 由于推力的方向与大小的多变性，其次方向是绝对会发生变化的，而大小可以等价的看成在降落过程中不变的，因此在软着陆过程就变得相对来说简单多了。还有就是关于这一过程的动力方程虽然从受力方面简化了，由-参考文献1周净扬，周荻，月球探测器软着陆精确建模及最优轨道设计，宇航学报，2007 年 11月2郭景录，付平。登月软着陆轨道优化算法研究。计算机仿真，2009 年 12 月3张洪华，梁俊，黄翔宇，赵宇，王立，关轶峰，程铭，李骥，王鹏洁，于洁，袁立。嫦娥三号自主避障软着陆控制技术。中国科学：技术科学，2014 年 6 月4

9、单永正，段广仁，吕世良。月球探测器软着陆的最优控制。光学 精密工程，2009 年9 月5李刘强，罗建军，谢剑锋。基于双二体模型的载人登月着陆点选取分析与仿真。科学技术与工程，2013 年 5 月6吴智友, 关于月球探测器的最优设计 http:/ 访问时间（2014 年 9 月 12 日）附 录拟合函数对比，如下表所示：表 6.1 拟合函数对比表Matlab 源代码：1) 最小二乘拟合曲面.m:clear allclc %清屏及清除记录a=imread(C:UsersAdministratorDesktopf1.tif,tif);7b=gradient(double(a);B=mat2cell(b,ones(100,1)*2300/100,ones(100,1)*2300/100);%将矩阵分成 100 块z=reshape(B1,1,23*23); %将第一部分矩阵转化为一维向量d=1; %求 xg=23;i=0;for l=1:23for n=d:gx(n)=1+i; endd=d+23;g=g+23;i=i+1;end