中学生数学杂志2014年全国大学生数学建模A题论文

2020年11月19日发(作者：瞿颖)

2012高教社杯全国大学生数学建模竞赛

承 诺 书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白，在竞赛开始后参赛队员 不能以任何方式（包括电话、电子邮件、网
上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛 题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的< br>资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参
考文献中明 确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规
则 的行为，我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式 进行公开展
示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B
我们的参赛报名号为（如果赛区设置报名号的话）：
所属学校（请填写完整的全名）： 大学
参赛队员 (打印并签名) ：1.
2.
3.
指导教师或指导教师组负责人 (打印并签名)：

日期： 2014 年 9 月 15 日



赛区评阅编号（由赛区组委会评阅前进行编号）：



2012高教社杯全国大学生数学建模竞赛

编 号 专 用 页




评
阅
人


评
分


备
注


赛区评阅编号（由赛区组委会评阅前进行编号）：



赛区评阅记录（可供赛区评阅时使用）：







全国统一编号（由赛区组委会送交全国前编号）：





全国评阅编号（由全国组委会评阅前进行编号）：





嫦娥三号软着陆轨道设计与控制策略
一、摘要

本 文针对嫦娥三号软着陆轨道设计与控制策略的实际问题，以理论力学（万有引力、开
普勒定律、万能守恒 定律等）和卫星力学知识为理论基础，结合微分方程和微元法，借
助MATLAB软件解决了题目所要求 解的问题。
针对问题（1），在合理的假设基础上，利用物理理论知识、解析几何知识和微元法，分
析并求解出近月点和远月点的位置，即139.1097 。再运用能量守恒定律和相关数< br>据，计算出速度1v（近月点的速度）=1750.78/ms,2v（远月点的速度）=1669.77 /ms，，
最后利用曲线的切线方程，代入点（近月点与远月点）的坐标求值，计算出方向余弦即
为相应的速度方向。
针对问题（2）


关键词：模糊评判，聚类分析，流体交通量，排队论，多元非线性回归

二、问题重述

嫦娥三号于2013年12月2日1时30分成功发射，12月6日抵 达月球轨道。嫦娥
三号在着陆准备轨道上的运行质量为2.4t，其安装在下部的主减速发动机能够产生
1500N到7500N的可调节推力，其比冲（即单位质量的推进剂产生的推力）为2940m/s，
可以满足调整速度的控制要求。在四周安装有姿态调整发动机，在给定主减速发动机的
推力方向 后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。嫦娥三号
的预定着陆点为19.51 W，44.12N，海拔为-2641m（见附件1）。
嫦娥三号在高速飞行的情况下，要保证准确 地在月球预定区域内实现软着陆，关
键问题是着陆轨道与控制策略的设计。其着陆轨道设计的基本要求： 着陆准备轨道为近
月点15km，远月点100km的椭圆形轨道；着陆轨道为从近月点至着陆点，其软 着陆过程
共分为6个阶段（见附件2），要求满足每个阶段在关键点所处的状态；尽量减少软着陆
过程的燃料消耗。
根据上述的基本要求，请你们建立数学模型解决下面的问题：
（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与
方向。
（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。
（3）对于你们设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。


1


二、问题分析
2.1问题（1）的分析
首先根据问题的假设、题目中所提供的数据及图片分析，可以知道嫦
娥三号绕月球的轨道是由 圆形轨道变为椭圆形轨道，借助开普勒定律、能量守恒定律求
解出近月点的速度。
为了确定 近月点和元月点的精确位置及相应的速度方向，我们建立以赤道（月球的
赤道）平面为xoy平面、月心 为原点、月心与零度经线和零度纬线交线的交点的连线为
坐标轴的坐标系和赤道（月球的赤道）平面为x oy平面，为极轴（月球的极轴）为z轴
建立空间直角坐标系，x轴与极坐标系的轴相重合。
)的空
间直角坐标。
其次利用两点间的距离公式，并借助MATLAB软件求解出近月点与着陆点最短距离。
=。
最后利用卫星的轨迹是以月心为其中一个焦点，以近月点与远月点的距离为长轴的
椭圆，从而求解出卫星 的轨迹方程，再运用隐函数求导的应用的知识，求解出在近月点
和远月点的方向导数，进而求解近月点和 远月点方向余即为近月点和远月点的速度的方
向。
2.2问题（2）的分析
首先在根据题意，将嫦娥三号软着陆问题，分为6个阶段依次为主减速、
快速调整、粗避障、精避障 、缓慢下降、自由下降，我们先将6个阶段分为4个阶
段，依次为第一阶段（主减速和快速调整）、第二 阶段（粗避障）第三阶段(精避障),第
四阶段（缓慢下降和自由下降）。
其次在第一阶段
粗避障阶段，嫦娥三号悬停在月球表面约2400米上方，对星下月表进行二维和三维
成像，利 用遗传算法的思想，从图像中先随机选取部分点，能直接从三维图像中得知该
点的海拔高度，再分别扫描 这些点附近的地貌，找出一些地势平坦的区域，我们用区域
内所有点与中心点海拔的均方差作为地势判断 依据之一，保留这些坐标， 并进行重新
组合，并改变某些坐标以便能获得其他新区域的坐标，再次搜索 地势平坦的区域，重复
进行多次搜索，直到没有出现崎岖地势的时候， 我们将此时地势最平坦的地方作为全
局最优降落地点。

三、模型假设及符号说明

2


1、不考虑空间飞行器上各点因燃料消耗而产生的位移；
2、在对卫星和空间飞行器进行轨道估计时，认为作用于其上的所有外力都通过其质
心；
3、卫星和空间飞行器的运动是在真空中进行的；
4、卫星只受重力影响，空间飞行器除自身推力外只受重力影响； 5、卫星的观测
图片及数据精准；
四、变量与符号说明
C
一条车道的基本通行能力 L
连续车流的车头间距 C n 条车道的基本通行能力 y
排队长度 1x 车流量
2x 横断面通行能力系数车流量 3x
持续时间
五、模型建立与求解
5.1 问题(1)的分析、模型建立与求解 5.1.1建模准备
（1）开普勒定律
开普勒第一定律开普勒第一定律开普勒第 一定律，也称椭圆定律：每一个行星都沿
各自的椭圆轨道环绕太阳，而太阳则处在椭圆的一个焦点中。开 普勒第二定律开普勒定
律开普勒第二定律，也称面积定律：在相等时间内，太阳和运动着的行星的连线所 扫过
的面积都是相等的。 这一定律实际揭示了行星绕太阳公转的角动量守恒。用公式表示为
开 普勒定律开普勒第三定律开普勒定律开普勒第三定律，也称调和定律：各个行星绕太
阳公转周期的平方和 它们的椭圆轨道的半长轴的立方成正比。由这一定律不难导出：行
星与太阳之间的引力与半径的平方成反 比。这是牛顿的万有引力定
律的一个重要基础。用公式表示为3
2aKT

这里，是行星公转轨道半长轴，是行星公转周期，是常数 。
（2）万有引力
万有引力：任意两个质点有通过连心线方向上的力相互吸引。该引力大小与它们
质量的乘积成正比与它们 距离的平方成反比，与两物体的化学组成和其间介质种类无
关。即：
12
2MMFG
r 其中12MM，为两物体的质量，11
226.6710..GNmkg
克）
5.1.2 模型的建立
根据 以上的分析，建立以月球赤道平面为xOy平面,月心为原点O、Ox为月心与零
度经线和零度纬线交线 的交点的连线，Oz为极轴（月球的极轴），Oy与Ox和Oz满足
右手标架，建立空间直角坐标系(如 图5-1所示）。
由于着陆点在球面上且近月点与远月点是由月球的经度、纬度及高度唯一确定，在

3


此为了便于计算 将极坐标转化为空间直角坐标，并代数题中相关数据，反解出经度
极坐标转化为空间直角坐标。




4






5





6





7





8