2014高教社杯全国大学生数学建模竞赛A题

2014高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。

如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。

我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A
我们的报名参赛队号为（8位数字组成的编号）：07033001 所属学校（请填写完整的全名）：吉林师范大学博达学院
参赛队员(打印并签名) ：1.
2.
3.
指导教师或指导教师组负责人(打印并签名)：
（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。

以上内容请仔细核对，提交后将不再允许做任何修改。

如填写错误，论文可能被取消评奖资格。

）
日期： 2014 年 9 月 15 日赛区评阅编号（由赛区组委会评阅前进行编号）：
2014高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号（由赛区组委会评阅前进行编号）：
全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：
嫦娥三号软着陆轨道设计与控制策略
摘要
本文针对嫦娥三号软着陆轨道设计与控制策略问题，通过提取题目中的信息，利用拱点的概念、B 样条函数逼近的统计定位方法、非线性规划问题及哈密尔顿函数为理论基础进行了完整的建模工作。

首先，通过建立坐标系结合物理学运动公式求解出了近月点与远月点的位置及相应的速度；在此基础上，利用B 样条函数逼近的方法确定了嫦娥三号的着陆轨；最后通过分解着陆过程并利用非线性规划问题及哈密尔顿函数确定着陆阶段的最优控制策。

针对问题一，利用拱点的概念及物理学中天体的运动的方程先求解出了嫦娥三号在近月点和远月点得速度大小分别为：1692.2m/s 、1613.9m/s 及方向为运行轨道方向的切线方向，并利用复平面上建立坐标系的方法建立了极坐标，运用椭圆轨道方程进而得出了近月点与远月点的位置分别为：
)10013.17520183⨯︒，（，)10013.183703⨯︒，（
针对问题二，采用B 样条函数逼近的运动学统计定位方法确定了在着陆弧段
上任意时刻的位置方程，从而刻画出了嫦娥三号的着陆轨道，
并用matlab 对轨迹进行了模拟。

在6个阶段的最优控制策略上，先通过直角坐标系得出质心的运动方程，再通过对6个阶段初始条件和终端状态的分解，利用非线性规划问题求解哈密尔顿函数，得出性能指标（耗燃量）的最小值为：382.6531kg ，从而确定了最优控制策略。

针对问题三，对于设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。

通过对月球半径、陨石坑对嫦娥三号着陆所需燃料及轨道的影响、着陆过程中消耗燃料对嫦娥三号运行质量的影响三个方面进行了误差分析，合理解释了所求数据与实际数据之间的误差。

又通过对问题二所采用的B 样条函数逼近方法中的参数进行敏感性分析，得出模拟着陆轨道的最佳节点数。

关键词：
拱点复平面建立坐标系 B 样条函数非线性规划哈密尔顿函数
一、问题重述
1.1引言
航空航天技术对一个国家的综合国力发展具有重要的战略意义，一直是世界各国重点建设的内容。

中国的航空航天技术创建于1956年，50年来，在中国人民的努力下取得了伟大成就，为国家的科技发展、经济建设和国家安全作出了巨大贡献。

嫦娥三号是我国探月工程第二期的第二颗月球探测器，在此之前，我国已成功实施了探月的一期工程，发射了嫦娥一号，并发射了探月工程二期的嫦娥二号探测器，为嫦娥三号的成功发射打下了坚实的基础。

嫦娥三号的软着陆，与嫦娥一号和嫦娥二号不同的是，嫦娥三号使用了先进大推力变推力发动机。

嫦娥三号使用了专门研制的1500～7500N变推力发动机作为推进系统的主发动机，这台发动机不仅用于完成嫦娥三号直接奔月轨道飞行期间的中途修正、近月制动和变轨任务，还用于嫦娥三号的软着陆。

可见，嫦娥三号的软着陆过程非常复杂，变推力主发动机在其中起到了至关重要的作用，也是嫦娥三号研制的主要难点之一。

1.2问题的提出
要保证嫦娥三号能够准确的实现软着陆，其关键问题是着陆轨道与控制策略的设计。

其着陆轨道设计的基本要求：着陆轨道为近月点15km，远月点100km 的椭圆形轨道；着陆轨道从近月点至着陆点，其过程共经历6个阶段，要求满足每个阶段所在关键点的状态，并尽量减少着陆过程的燃料消耗。

其着陆点大致位置：19.51W，44.12N，海拔-2641m。

为了更好的找到嫦娥三号软着陆的轨道设计与控制策略，本文依次提出以下的几个问题：
（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小和方向。

（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。

（3）对设计的着陆轨道和控制策略做出相应的误差分析和敏感性分析。

二、模型假设
1.月球引力场均匀。

2.月球自转对嫦娥三号的着陆无影响。

3.嫦娥三号运行的轨道是以月球球心为一个焦点的椭圆。

4、月球的形状扁率对嫦娥三号的着陆轨道无影响。

三、 符号说明
在近月点时的速度
a v 在远月点时的速度
a 半主轴长
*
e 椭圆轨道离心率
μ 月球标准重力参数值
a r
远月点到月球球心的距离 b r 近月点到月球球心的距离
G 万有引力常数
M 月球质量
m 嫦娥三号运行质量
四、 问题分析
问题一：
问题要求确定近月点与远月点的位置，以及嫦娥三号相应速度大小和方向。

首先，可以确定嫦娥三号近月的运动轨迹是一个椭圆，其中把月球球心作为椭圆长轴负半轴的焦点，则距离焦点较近的长轴交点为近月点，距离焦点较远的长轴交点为远月点。

利用椭圆运动的天体方程，可求出嫦娥三号在近月点及远月点的相应速度。

再以月球球心为坐标原点建立极坐标，通过建立常微分方程并通过matlab 软件进行求解，确定了嫦娥三号在近月点及远月点的相应位置。

问题二：
由问题一解出的准备着陆的位置和着陆的6个阶段的运动特点，可以推算出这6个阶段关键点的运动状态及大致位置，再利用B 样本函数逼近的统计定位方法得出着陆轨 迹上任意时刻探测器的位置，从而模拟出了嫦娥三号的着陆轨道。

在最优控制策略上，本文采用非线性规划求解月球软着陆最优控制中的两点边值问题，构造性能指标函数，由哈密尔顿函数取极小值的思想得出最优着陆曲线，进而对函数 进行积分求得最终消耗燃料，即最少燃料。

问题三：
通过对月球半径、陨石坑对嫦娥三号着陆所需燃料及轨道的影响、着陆过程中消耗燃料对嫦娥三号运行质量的影响三个方面进行了误差分析，合理解释了所求数据与实际数据之间的误差。

同时，通过对问题二所采用的B 样条函数逼近方p v
法中的参数进行敏感性分析，得出模拟着陆轨道的最佳节点数。

五、 模型的建立与求解
5.1问题一的建模与解答
5.1.1 嫦娥三号在近月点及远月点相应速度的大小和方向
Ⅰ 模型的建立
通常，可以把月球作为椭圆的焦点，把与焦点（月球）距离较近的交点叫做近月点，把与焦点（月球）距离较远的交点叫做远月点，如图1所示：
1.拱点
拱点在天文学上的意义是在椭圆轨道上运行的天体最接近或最远离它的引力中心，通常也就是系统的质量中心的点。

最靠近引力中心的点称为近拱点（也就是本文中的近月点）而距离最远的点就称为远拱点（也就是本文中的远月点）连接近拱点和远拱点的直线称为拱点线，是椭圆的长轴，也是椭圆内最长的线段。

近拱点（近月点）：在最短距离处有最快的速度即
a
e e v p )1()1(**-+=μ （4） 远拱点（远月点）：在最远距离上有最慢的速度即
a e e v a )1()1(**+-=μ
（5）
其中是标准重力参数值是半主轴，μa 。

2.轨道离心率
121+-=+-=*p a p a p
a r r r r r r e （6）
在太空动力学上，一个天体的标准重力参数是万有引力常数和它质量之积即
GM =μ （7）
其中：a r 远月点到月球中心的距离；p r 近月点到月球中心的距离。

Ⅱ模型的求解
联立（4）（6）（7）式可得近月点（近拱点）的速度：
a r GM
r v p a p =
（8）
联立（5）（6）（7）式可得远月点（远拱点）的速度：
a r GM
r v a p a = （9）
通过matlab 计算（程序详见附录）可得：
s m v p /106922.13⨯=
s m v a /106139.13⨯=
于是，可得知嫦娥三号在近月点的速度为s m v p /106922.13⨯=，方向是沿运动的切线方向；嫦娥三号在远月点的速度为 s m v a /106139.13⨯=，方向是沿运动切线方向。

5.1.2着陆准备轨道近月点与远月点的位置
Ⅰ 模型的建立
嫦娥三号经过轨道修正、减速制动、近月制动，从而进入环月轨道。

环月轨道可以确定为是个椭圆轨道（如图2所示）。

我们要确定近月点与远月点的位置，就必须研究嫦娥三号软着陆过程中的主减速过程。

在此，我们要确定嫦娥三号着陆准备轨道近月点与远月点的位置，就可以先确定嫦娥三号运动的轨迹方程。

通过轨迹方程的求出，从而确定近月点与远月点的具体位置，进而得到嫦娥三号着陆准备轨道近月点与远月点的位置。

设月球中心所在的位置为复平面的原点O ，在时刻t ，嫦娥三号位于θi re t Z =)(所表示的点P 。

这里)(),(t t r r θθ==均为t 的函数，分别表示)(t Z 的模和辐角]1[。

于是嫦娥三号的速度为： )(dt d ir dt dr e dt d ire e dt dr dt dZ i i i θθθθθ+=+= 加速度为：
⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛-=dt d dt dr dt d r i dt d r dt r d e dt Z d i θθθθ22222222
（10） 月球对嫦娥三号的引力依据物理学中的万有引力定律，大小为
2
r mMG F =引 方向指向月球的中心O ，故有，
θi e r
mMG F 2-=引 其中M 为月球质量，m 为嫦娥三号质量。

依据牛顿第二定律我们可得
222dt Z d m e r mMG i =-θ
（11） 将公式（10）带入（11），然后比较其实部与虚部就得到：
22222
220()d dr d r dt dt dt d r d MG
r dt dt r θθθ⎧+=⎪⎪⎨⎪-=-⎪⎩
这是两个未知函数的二阶微分方程组。

在确定某一卫星轨道时，需要加上定解条件。

假设当t=0时，卫星正处于远月点，而远月点位于正实轴上，距原点O 为0r 卫星的速度为0v 那么就有初值条件：
⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧========0000000|0|0||r v dt d dt dr r r t t t t θθ
因此问题转化为求解带初值问题的微分方程组为：
⎪⎪⎪⎪⎪⎪⎪⎪⎪⎩
⎪⎪⎪⎪⎪⎪⎪⎪⎪⎨⎧====-=-=+====00000002
22222|0|0||)(02r v dt d dt dr r r r MG dt d r dt r d dt d dt dr dt d r t t t t θθθθθ 又将0222=+dt d dt dr dt
d r θθ两边同乘以r 即得： 12c dt d r =θ
（12）
其中001v r c =。

这样有向线段p o 在时间t ∆内扫过的面积等于
2
2112t c dt dt d r t
t t ∆=⎰∆+θ 这正是开普勒第二定律，从月球指向嫦娥三号的线段在单位时间内扫过的面积相等。

将（12）式代入
222)(r
GM dt d r dt r d -=-θ 得：
2321221r
MG r c dt r d -=- 于是我们可以得到嫦娥三号运行较为简单形式的数学模型：
⎪⎪⎪⎪⎪⎪⎪⎩
⎪⎪⎪⎪⎪⎪⎪⎨⎧====-=-===0|0||00002
1232122t t t dt dr r r r c dt d r MG r c dt r d θθ 对上式进行化简得到：
)cos(10θθ--=*e p
r
其中 A 和0θ是待定参数；Ap e =*；MG c p 21=。

上式即为嫦娥三号的轨道方程，是一条平面二次曲线。

由于嫦娥三号绕月球运行，故必有0〈*e 〈1.由于r 在t=0时取到最大值0r （远月点），意味着此时 函数)cos(0θθ-取到最大值1，于是有：
01,0r p
e -
==*θ （13） 从而轨迹方程为：
（14） 对于嫦娥三号而言，又知嫦娥三号的近月点到月球的距离为：
**+=
-=
e p
e p r m 1cos 1π
Ⅱ 模型的求解
由（13）式 化简可得：
0)1(r e P *-=
利用matlab 求解得（程序详见附录）：
6107935.1⨯=P
将所得数值代入（14）式中， 即可得到θcos 与P 之间的关系：
θ
cos 0237.01107935.16
⨯-⨯=r
又因为，a r 表示远月点到月球中心的距离，p r 表示近月点到月球中心的距离，那么即可求得：
当m r a 3
10013.1837⨯=时，解得1cos =θ，则︒=0θ；
当m r 3
p 10013.7521⨯=时，解得1-cos =θ，则︒=180θ。

则在远月点的位置是)10013.183703⨯︒，（；近月点的位置是)10013.17520183⨯︒，（。

根据所求位置，利用matlab 做出图像（程序详见附录）：
θ
cos 1*
-=
e p
r
图3 嫦娥三号绕月运行的轨道曲线
5.2 问题二的建模与解答
嫦娥三号的软着陆过程是从近月点到月面的全过程。

嫦娥三号的软着陆过程共分为6个阶段，前3个阶段主要是利用火箭的反推力制动减速，消除探测器较大的水平速度，通过姿态调整发动机调整探测器姿态为竖直向下，并使探测器在到达预定高度时的速度接近0；第四阶段，悬停避障段主要是根据图像敏感器对着陆区的成像选择安全着陆点；第五阶段，缓速下降段是保证探测器在离月面约4米时速度为0，从而关闭发动机；最后一阶段探测器以自由落体方式撞机月球表面着陆。

（具体阶段详见图4）
5.2.1 嫦娥三号的着陆轨道 Ⅰ 模型的建立
B 样条函数：在曲线拟合设计中，B 样条曲线主要用于要求局部交互式修改的自由曲线设计和实验数据平滑。

设）
（k n ,,21≥n P P P 为给定空间的n 个点，称下列参数曲线 1n k k i i n
1
i t t t t )t (+=≤≤∑=），（，B P P
为k 阶（k-1次）的B 样条曲线。

称折线n P P P ,,21为)(t P 的控制多变形
（如图5）。

图5
本文采用B 样条函数逼近的运动学统计定位方法对落月的轨迹进行仿真刻画，对着陆轨迹确定的策略进行了分析。

1. 计算原理：
三阶B 样条具有良好的二阶光滑度，拟合的灵活性比较强，适合拟合曲率变化大、拐曲严重的任意形状函数，而且具有较好的稳定性和收敛性。

本文选择三阶标准B 样条，函数形式如下：
⎪⎪⎪
⎩⎪⎪⎪
⎨⎧
≤≤+-+-≤+-≥=21346
11
32
21202323t t t t t t t t B 将样本数据处理区间为[]
12,-p T T ，P 为节点数 则：
B 样条曲线及其控制多变形
⎪⎪⎪
⎪⎪⎪⎪
⎩
⎪
⎪
⎪⎪⎪⎪⎪⎨
⎧
⎪⎪⎭
⎫ ⎝
⎛-=⎪
⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝
⎛-=⎪⎪⎭⎫
⎝⎛-=⎪⎪⎭⎫
⎝
⎛-=⎪⎪⎭⎫ ⎝⎛-==-+=--=∑∑∑∑∑∑=•
•
==••
==•
•
=-p
j j
j p
j j j p
j j
j p j j j p
j j
j p
j j j P h
T t B t z h T t B t z h
T t B t y h T t B t y h
T t B t x h T t B t x P
j h j T T P T T h 1111112
21)()()()()()(,2,1,
)2(3γγββα ⎪⎪⎪⎩⎪⎪⎪
⎨⎧
≤≤-+-≤-≥=21222
112232022t k t k t t t k t t B
其中：)(t sign k =。

以上可以得出：该弧段上任意时刻探测器的位置便可用B 样条函数表示。

2. 采用B 样条函数逼近方法的本质：
从统计定位的思想出发，把需要定位的弧段分成若干个小区间，在区间上用B 样条函数逼近，然后利用最小二乘法求解各节点的参数。

可调参数P 越大，样本数据区间就会分的越小，求解出的参数也会越多]2[。

Ⅱ 模型的求解
用matlab 程序作图得到嫦娥三号的着陆轨道如图6（matlab 程序详见附录）。

图6 定点直线拟合与模拟仿真模型
5.2.2 6个阶段的最优控制策略 Ⅰ 模型的建立
建立着陆坐标系如图7所示：
取月球球心为坐标原点；oy 指向着陆准备轨道的近月点；r 为探测器到球心的距离；θ为oy 与or 的夹角；ϕ(t)为推力方向与垂线的夹角；F 为制动火箭的推力大小，F 取值为[1500-7500]或0。

根据题意及运动轨迹得出嫦娥三号的质心运动方程：
⎪⎪⎪⎪⎪
⎪⎪⎪⎪
⎪⎩⎪⎪⎪
⎪
⎪
⎪⎪⎪⎪
⎪
⎨⎧
-
=+-==+-==••
•••
C F m v m F
r r r m F v v r )
2cos (1sin 2
2ωϕωωθωμ
ϕ （1）
其中，v 为探测器在r 方向上的速度；ω为探测器方位角θ的角速度；m 为探测器的质量；μ为月球的引力常数；C 为制动火箭的排气速度，是一个常值。

嫦娥三号的软着陆分为6个阶段，根据近月点的位置和着陆的轨迹，可以推算出每个阶段的初始状态和终端状态：
表1：着陆各阶段状态
在上表中，r 为月球半径；6144.3=f v （h 22月g v f =，其中2s /m 633.1g =月）探为测器到达月面是的速度。

通过matlab 编程得到着陆过程中月心距与时间的关系（程序详见附录）
图8 月心距变化曲线
0100200
300
400500600
1738
17401742174417461748
175017521754 月心距变化曲线
时间
月心距
初始条件
终端条件
)(t r
)(t v
)(t θ )(t ω )(t m )(t r
)(t v
)(t θ )(t ω )(t m
主减速 15000+r p v
1r v p
m 3000+r 57
1θ 2r v p
1m
快速调整
3000+r 57
1θ 2r v p
1m
2400+r 57
2θ 3
r v p
2m
粗避障
2400+r 57
2θ 3
r v p
2m
100+r 0
3θ
3m
精避障
100+r
3θ
3m
30+r 1v 4θ
4m
缓速下降
30+r
4θ
4m
4+r
3.6144 5θ
5m 自由落体
4+r
3.6144 5θ
0 5m
r
3.6144
2
π
5m
约束条件
变量
阶段
在计算过程中，由于状态变量的量级相差很大，在积分过程中会导致有效位数的丢失，所以本文采取归一化处理来提高计算的精准度，这样也可以使变量保持在同一量级。

因此，做以下处理：令
⎪⎪⎪⎩
⎪
⎪⎪⎨
⎧
======000
0///,//m m m v
r r v v v v r r r ωωθθμ， 及：⎩⎨⎧====02
00/,//,/r v m F F F F v
r t t t t
则质心运动方程（1）可以改写为
⎪⎪
⎪⎪
⎪
⎪⎪⎪⎩⎪⎪⎪
⎪⎪⎪⎪
⎪
⎨
⎧
-=+-==+-==••
•••
02
2)
2cos (11
sin r C F m v m F
r r r m F v v r μωϕωωθωϕ
（2） 最优耗燃问题：
运动方程（1）的状态方程[]
3为),(u x f x =•
，其中状态变量T m v r x ],,,,[ωθ=，控制变量T F u ],[ϕ=。

终端性能指标
)
(1)()0()]([f f f t m t m m t x J -=-==ϕ
（3）
若求耗然最少，则需构造哈密尔顿函数
),(),,(u x f x H T λμλ= （4）
其中T m v r ],,,,[λλλλλλωθ=，满足
x u x H ∂∂-
=•
)
,,(λλ
（5）
即：
⎪⎪⎪⎪⎪⎪⎪⎩⎪
⎪
⎪⎪⎪⎪⎪⎨
⎧-=+--==+-=----=•
•••
•ϕ
λϕλλλλωλλλωλλλωλϕλωλλλωωθωθω
ωωcos sin 22022cos 2222
223r m F
m F r v r r r v r m F r v m v r v v v r
（6） 终端约束为：
0)(/)(/)(0=⎥⎥⎥
⎦
⎤⎢⎢⎢⎣⎡--=f f f f f t v v t v r r t r G ω
由（6）式知，沿最优轨迹有0≡θλ，得出：
F C
r m m r v r r v H m v v v
r )cos sin (222
λ
ϕλϕλωλωλμ
λλωω
--+-+-= 由于ϕ不受约束，因此由
0=∂∂ϕ
H
得最优制导律 T F u ],[***ϕ= （7）
将（7）式代入（2）、（6）式积分可得最优轨道。

Ⅱ 模型的求解
根据非线性规划解两点边值问题的方法，可以把性能指标J 看作[]T
t v r 0
,,=ωλλλ与f t 的函数，即可把问题转化为：
优化目标函数：)]([f t x J ϕ=；
优化参数：[]T
t v r 0,,=ωλλλ，f t ；
约束：⎩⎨
⎧
=）积分）、（
式（620G 即:
dt t C F J f
t t ⎰=0
)(/
其中e v 为以米/秒为单位的比冲。

利用matlab 编程解得:J= 382.6531kg(程序详见附录)
5.3问题三
5.3.1 误差分析
本文在计算15km 处速度时采用的月球半径为月球的平均半径，而此位置的半径应该1735.843km ≤R≤1737.646km ，所以导致计算速度与实际速度有微小的偏差。

本文在燃料最优计算时所用的时间为750s ，而实际的时间应该比750s 还要小，如：嫦娥三号在离地面4m 时就关闭了发动机，之后就做自由落体运动。

根据：
有：（2
/633.1s m g ≈） 代入数据得：
实际用时要比750s 小。

并且随着嫦娥三号的燃料消耗，会导致其质量随时间的增加而减小，从而影响发动机的推力大小，导致计算所需的燃料与实际所用燃料的偏差。

嫦娥三号进行实施软着陆中，时间作为主要影响因素。

除此之外，还受到地形因素的影响。

所用燃料的实际值还与嫦娥三号软着陆位置的地形有关；因为软着陆的地形会影响嫦娥三号姿态调整发动机所需燃料的多少；当遇到的陨石坑面积大时，水平移动的距离就大，姿态调整发动机所需的燃料就多；当遇到的陨石坑面积小或地势平坦时，水平移动的距离就小，姿态调整发动机所需的燃料就少（见图9,、图10）。

22
1
gt h =g
h t 2=
s t 2134.2≈
图9 嫦娥三号着陆过程中距月面2400m 处月面情况变化趋势图
图10 嫦娥三号着陆过程中距月面100m 处月面情况变化趋势图
1.像素
像素是指基本原色素及其灰度的基本编码。

像素是构成数码影像的基本单元，
05001000150020002500
50
100
150
200
250
嫦娥三号着陆过程中距月面2400m 处月面情况变化趋势图
01002003004005006007008009001000
50
100
150
200
250
嫦娥三号着陆过程中距月面100m 处月面情况变化趋势图
通常以像素每英寸PPI 为单位来表示影像分辨率的大小。

通常，一个像素被视为图像的最小完整采样。

根据附录3、附录4所给图片。

我们发现图片阴影部分较深的地方为陨石坑，所以对应着色较深。

图片中成像较浅的部分，可认为月球表面且较平坦。

把两张图片所对像素点作为图片的基本单位，我们发现像素所对应单位可以做成一个
n m ⨯的矩阵记为n A 。

⎥⎥⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎢⎢⎣⎡=mn m m m n n n n a a a a a a a a a a a a a a a a A
3
2
1
3333231
2232221
1131211
然后再依据每个像素点作为单位，以每一行的像素点（地面平坦程度）变化
作为周期，做出整张图片的像素变化趋势。

5.3.2 敏感性分析]4[
在问题二中，确定嫦娥三号着陆轨道时用到了B 样条函数，而本文只利用题中涉及的7个节点确定了着陆轨迹，事实上，当节点个数增多时，会解得更多参数，从而更利于确定轨迹。

下表示为取不同节点数时得到的位置和速度的数据：
平方根。

统计学计算公式为：n
x x x n
2
2
22
1+++
从这些数据可以看出，可调参数（节点）从10增加到30是，位置误差的RMS 再逐渐减小，P 取到30是，位置误差的RMS 为80m 左右，但是，继续增加节点数，P 增大到50时，位置误差的RMS 并没有继续减小，而是增大到了90m 左右。

这是因为如果P 增加，需要求解的参数也会怎多，而求解参数N 与P 的关系为3)3(⨯+=P N ,所以，P 值太大，反而会影响计算效率，因此，取30=P 为宜。

六、 模型的评价
6.1模型的优点
1. 模型一.将嫦娥三号的运行轨道方程化[]5，可以直观的求出运行轨道上的
每一点的具体坐标。

将近拱点（近月点）与远拱点（远月点）公式化。

试用于运动轨道为椭圆的天体之间的运动。

2．模型二利用B样条差值法，处理嫦娥三号软着陆过程的拟合。

使试验结果更加逼真，准确。

3.通过使用matlab[]6，excel，Word作图，来处理相关数据，使抽象的具体
数值变成更加详尽的图像，更加直观的观察出数据变化趋势与规律。

并且，利用matlab来进行大量的数据处理。

6.2模型的缺点
1. 模型具有局限性，不适用于运动轨道为圆形的天体间的运动。

2．本文所建模型，只适用于软着陆过程
参考文献
[1] 钟玉泉复变函数论第四版北京：高等教育出版社.
[2] 昌胜骐，黄勇，宋叶志，嫦娥三号动力落月段轨迹确定策略飞行器测控学报2014,33（3）.
[3] 单永正，段广仁，吕世良月球探测器软着陆的最优控制光学精密工程第17卷第9期2009
年九月.
[4] 姜启源，谢金星，叶俊. 数学模型. 第三版. 北京：高等教育出版社，2003
[5] 东北师范大学微分方程教研室常微分方程第二版北京：高等教育出版社，2012年12月
P157-169
[6] 刘卫国MATLAB程序设计教程第二版北京：中国水利水电出版社，2005
附录：
v
1.求
p
clear
G=6.67e-011;
M=7.3477e+022;
ra=1837.013;
rp=1752.013;
a=1794.513e+003
vp=sqrt((ra*G*M)/(rp*a))
ansvp =
1.6922e+03
v
2.求
a
clear
G=6.67e-011;
M=7.3477e+022;
ra=1837.013;
rp=1752.013;
a=1794.513e+003
va=sqrt((rp*G*M)/(ra*a)) ansva =
1.6139e+03
r所对应的角度
3.求
a
clear
ra=1837013;
rp=1752013;
e=(ra-rp)/(ra+rp);
p=(1-e)*ra;
ra=1837013;
e=(ra-rp)/(ra+rp);
p=(1-e)*ra;
coso=((1/e)-(p/(ra*e))); y=acos(1)
ans y=
0.0000
r所对应的角度
4求
b
clear
ra=1837013;
rp=1752013;
e=(ra-rp)/(ra+rp);
p=(1-e)*ra;
coso=((1/e)-(p/(rp*e)));
y=acos(-1)
ans y=
3.1416
5.轨道曲线
clear;
plot(-1,0,'*');
hold on
plot(1,0,'*');
hold on
plot(-0.4,0,'*')
t=0:2*pi/30:2*pi;
x=cos(t);
y=sin(t);
plot(x,y)
title('图5 嫦娥三号绕月球运行的轨道曲线') 6.优化燃料
clear
s1=2940;
a=0;
b=750;
t=b-a;
s=1500;
y=s/s1*(b-a)
ans=
382.6531
7着陆过程中月心距与时间的关系
xi=0:100:600;
yi=[1753,1753.134,1752.932,1752.234,1748.921,1742.329,1738.143] plot(xi,yi)
grid on
title('图7 月心距变化曲线 ');
xlabel('时间');
ylabel('月心距 ')
8 球
[x,y,z]=sphere(20);
colormap(copper);
surf(x,y,z);
axis equal
grid off
9图像的导入
clear;
A=imread('C:\Users\Administrator\Desktop\1.tif');
plot(A);
title('嫦娥三号着陆过程中距月面2400m处月面情况变化趋势图')
grid on
clear;
B=imread('C:\Users\Administrator\Desktop\2.tif');
plot(B);
title('嫦娥三号着陆过程中距月面100m处月面情况变化趋势图')
grid on。