2014高教社杯全国大学生数学建模竞赛a题
摘要
2013年嫦娥三号成功发射，标志着我国航天事业上的又一个里程碑，针对嫦娥三号软着陆问题，分别建立着陆前轨道准备模型和软着陆轨道模型，建立动力学方程，以燃料最省为目标进行求解。
问题一：
在软着陆前准备轨道上利用开普勒定律、能量守恒定律以及卫星轨道的相关知识，利用牛顿迭代法分别确定了近月点和远月点的速度分别为1.6925km/s、1.6142km/s，位置分别为（19.91W，20.96N），（160.49E,69.31S）。
问题二：
在较为复杂的软着陆阶段，因为相对于月球的半径，嫦娥三号到月球的表面的距离太小，如果以月球中心建立坐标系会造成比较大的误差，因此选择在月球表面建立直角坐标系，在主减速阶段的类平抛面上建立相应的动力学模型，求出关键点的状态和并设计出相应的轨道，接下来通过利用灰度值阀值分割方法和螺旋搜索法对粗避障阶段和精避障阶段的地面地形进行相应的分析，找出安全点，然后调整嫦娥三号的方向以便安全降落，最后在落地时通过姿态发动机调整探测器的姿态，使之可以平稳的落到安全点上，在以上的各个阶段都可以以燃料最省为最优指标，从而建立非线性的最优规划的动力学模型，并基于该动力学模型可以对各个阶段的制导率进行优化设计由此就可以得到各个阶段的最优控制策略，
问题三：
最后针对所设计的轨道和各个阶段的控制策略进行了误差分析和灵敏度分析。对系统误差和偶然误差都做了解释；通过灵敏度分析发现，嫦娥三号在近月点的位置对结果的影响最大。
关键字 牛顿迭代法 ，灰度值阀值分割，螺旋搜索法，灵敏度分析
0
一、问题重述
嫦娥三号于2013年12月2日1时30分成功发射，12月6日抵达月球轨道。嫦娥三号在着陆准备轨道上的运行质量为2.4t，其安装在下部的主减速发动机能够产生1500N到7500N的可调节推力，其比冲（即单位质量的推进剂产生的推力）为2940m/s，可以满足调整速度的控制要求。在四周安装有姿态调整发动机，在给定主减速发动机的推力方向后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。嫦娥三号的预定着陆点为19.51W，44.12N，海拔为-2641m。 嫦娥三号在高速飞行的情况下，要保证准确地在月球预定区域内实现软着陆，关键问题是着陆轨道与控制策略的设计。其着陆轨道设计的基本要求：着陆准备轨道为近月点15km，远月点100km的椭圆形轨道；着陆轨道为从近月点至着陆点，其软着陆过程共分为6个阶段，要求满足每个阶段在关键点所处的状态；尽量减少软着陆过程的燃料消耗。
根据上述的基本要求，请你们建立数学模型解决下面的问题：
（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应的速度大小和方向。
（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。
（3）对于所设计设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。
二、问题分析 对于问题一，本文通过查找资料，发现嫦娥三号的运行轨道正好符合开普勒定律第一定律，其次根据开普勒第二定律和能量守恒定律，我们可以得出位于近月点以及远月点的速度大小以及方向，然后利用卫星轨道的相关知识，以月球赤道为平面建立空间直角坐标系，根据嫦娥三号的绕行轨道和赤道平面的的夹角计算出近月点和远月点的位置。 对于问题二，本文分为六段建立模型考虑问题，因为嫦娥三号距离月球地面的位置相对月球半径来说太小，所以我们在月球表面建立直角坐标系，根据的要求要在主减速阶段要求到达预订着陆点上方，利用抛物线相关知识建立精确动力学模型，用最优化方法求出结果，得到相应的再该阶段的控制策略。其次在粗避障和精避障阶段，利用距离地面2400米和100米的高程图，使用图像灰度值阀值分割方法和螺旋搜索法，将图中的不同高度的地面进行分割，分两次缩小安全点的位置，然后再最后下落过程中启动小型姿态发动机来进行水平调以便整最终安全着陆。
针对问题三我们从多个方面出发回归整个建模过程，对一些误差进行了分析，得到了减少误差的方法。
1
三、基本假设
1. 假设忽略月球的自传
2. 假设嫦娥三号飞行器的降落轨道和其绕月运行轨道在同一平面内 3. 假设在每阶段推力的大小不变
4. 假设忽略其他星球对嫦娥三号的引力作用
四、模型准备
开普勒定律是指行星在宇宙空间绕太阳公转所遵循的定律。
开普勒第一定律，也称椭圆定律；也称轨道定律：每一个行星都沿各自的椭圆轨道环绕太阳，而太阳则处在椭圆的一个焦点中。
开普勒第二定律，也称面积定律：在相等时间内，太阳和运动中的行星的连线（向量半径）所扫过的面积都是相等的。 这一定律实际揭示了行星绕太阳公转的角动量守恒。
开普勒第三定律，也称调和定律，也称周期定律：是指绕以太阳为焦点的椭圆轨道运行的所有行星，其椭圆轨道半长轴的立方与周期的平方之比是一个常量。这里，a是行星公转轨道半长轴，T是行星公转周期，K是常数，其大小只与中心天体的质量有关。常用于椭圆轨道的计算。即：
a3K?2 T其中，K?GM，M为中心天体的质量。 4?五、符号说明 a椭圆轨道长半轴 b椭圆轨道短半轴 c焦点
va嫦娥三号在远月点A时的速度 vb嫦娥三号在近月点B时的速度
符号说明
n平均角速度 M平均近点角 E偏心近点角 ?真近点角 i轨道倾角
?月球自传角速度
G宇宙常量约为
6.672?10?11Ngm2/kg2
?0零时刻的升交点经度 ?s(t)t时刻经度
M月球质量7.3477?1022kg T椭圆轨道周期
2
?s(t)t时刻纬度 F着陆轨道近月点推力
Fx推力水平分力 Fy推力竖直分力
v1主减速段结束速度57m/s
?主减速段结束速度方向
h1主减速段高度
ax推力F在水平方向的加速度
F1快速调整段推力
ay推力F在竖直方向的加速度
F1x推力水平分力
g月球重力加速度1.633m/s^2 （即单位质量的推进剂产生的推I比冲
力）为2940m/s
m嫦娥三号的质量2400kg
t0主减速段所用时间
F1y推力竖直分力
a1x推力F1在水平方向的加速度 a1y推力F1在竖直方向的加速度
v0主减速开始点既近月点速度 h2快速调整段高度
六、模型的建立与分析
1、问题一：近月点和远月点的速度大小、方向以及位置表示
近月点A和远月点,B距月球的距离分别为LA?a?c，LB?a?c,在极短的时间内嫦娥三号与月球连线扫过的面积分别为
11VSA?vAgVtgLA,VSB?vBgVtgLB
22根据开普勒第二定律;嫦娥三号与月球的连线在相等的时间内扫过相等的面
a?c积，故有VSA?VSB带入得 vB?vA。其次嫦娥三号运动的总机械能等于其动
a?c能和引力势能之和，故嫦娥三号分别经过A,B两点时的机械能为
?GMm?11GMm2EA?mvA2???， ??mvA?2LB?2a?c??GMm?11GMm2EB?mvB2????mv?. ?B2L2a?c?B?由于行星在运动过程中受万有引力作用，机械能守恒EA?EB?0，故由以上公式得出
3
vA??a?c?GM?a?c?a,vB??a?c?GM?a?c?a. 对模型进行求解，带入已知条件a2?b2?c2，得
bGMbGM, vB? a?caa?cavA?利用MATLAB求解得
vA=1.6925km/s
vB=1.6142km/s
方向分别沿近月点和远月点切线的方向，沿逆时针方向。
如图（1）所示，建立空间直角坐标系OXYZ，以月球为中心，赤道平面为XOY面，
根据卫星空间定位可近似计算出远月点和近月点的经纬度，
图1 坐标系中椭圆轨道
t为嫦娥三号从远月点开始出发的时间
a3椭圆运动周期为T?2? u平均角速度n?2? T其中卫星在t0通过近月点，并以平均角速度n绕椭圆轨道移动平均近点角
4
2014全国大学生数学建模a题
