【关键词】软着陆；节省燃油；坐标反推；MTALAB
1.问题重述
2.问题分析
这个问题属于航天着陆轨道设计。
（1）轨道面确定：飞行器可以视为一个小型天体，它在椭圆绕月轨道上需要服从开普勒第一定律、为节省燃油轨道面过预定着陆点、选定典型极地轨道，三个条件决定着陆准备轨道和月球19.51W经线圈位于同一平面。
（2）着陆阶段水平位移确定：软着陆6大阶段中主减速阶段（月面上空15km至3km）的水平位移占据总水平位移中的绝大部分，其余5大阶段多产生竖向位移或历时较短，所以仅需要考虑主减速阶段的水平位移；在设定飞行器初末俯仰角后，不断增大俯角，并以满推力7500N制动。
（3）主减速阶段推力控制策略：因比冲衡定，飞行器推力越大从而消耗燃油越快；为使得嫦娥三号减少制动耗能，应在高空尽可能多地使用燃油抵消机械能；进而在主减速阶段可以采用喷射器满推力7500N工作的控制策略。
（4）坐标反推：在计算主减速阶段水平位移后，依据主减速阶段平均高度（距离月心）折算维度变化，参考预定着陆点坐标得到近月点坐标。
3.模型假设
（1）在嫦娥三号动力下降过程中，忽略地球及其他引力对嫦娥三号的影响，只考虑月球对于飞行器的影响作用。
（2）因月球自转速度相对缓慢，可忽略自转影响与科里奥利力影响。
（3）考虑燃料的消耗对于嫦娥三号的质量的影响。
（5）月球表面没有大气层，几乎无空气阻力，嫦娥三号机械能守恒。
（6）在0.1s内，嫦娥三号的速度，加速度的大小和方向以及质量保持不变。
4.模型建立
4.1建立动坐标系
坐标系X轴以嫦娥三号进入主减速阶段为起始点，距离月球表面15km的圆轨道的切向方向为坐标轴原点，以起始点沿月球径向远离球心的方向Y轴的正方向。即坐标原点与嫦娥三号在同一月球半径上。
4.2力学模型与运动模型
嫦娥三号在运动过程中受到月球的引力与重力成正比，月球的引力指向月球圆心，分析嫦娥三号的运动状态，在椭圆的任一点处，嫦娥三号受到的月球引力均为F；第t时刻，飞行器的水平速度的大小为Vxt，竖直方向的速度大小为Vyt：
一方面，嫦娥三号在极地椭圆轨道的平面内做曲线运动，满足法向力的平衡：引力、惯性力、重力合力为零， 另一方面，嫦娥三号在曲线运动时切向满足牛顿第二定律。
质量变化由燃油损耗引起，由于比冲恒定等于2940m/s，所以单位时间消耗的燃料数量恒定，为2.55kg/s。
嫦娥三号主减速阶段的边界条件已知：高度由15公里降至3公里，初速度为1.695km/s，方向为水平（即椭圆轨道近月点处切向）。
4.3水平位移的计算
4.4近月点位置的确定
所以降落圈的最南端为C（19.41°W，31.37°N），最北端为D（19.41°W，56.87°N）。
选择由南向北的典型极地绕月轨道，近月点位置即可确定：（19.41°W，31.37°N），月面上空15km处，即可进入主减速阶段。
4.5远月点位置及近远月点速度
远月点由月面坐标的中心对称可得，为（340.59°E，31.37°N）；由于月球球心处于椭圆轨道的一个焦点上，近、远月点的速度可由角动量守恒、机械能守恒两大定律联合解得：近月点速度为1.695km/s，远月点速度为1.617km/s。
在求得近、远月点位置与速度后，若要模拟实际着落点，应以此计算结果为初始条件，分主减速、快速调整、粗避障、细避障等6阶段进行计算。
5.误差分析
5.1月球自转影响
在大约750s的软着陆过程中，月球自转经过的角度为0.11°；由月球自转方向自西向东，故近月点应向东修正0.11°。
5.2科氏力影响
科氏力作用下飞行器偏向223m，折算为经度偏差即0.01°，方向与自转相反。
所以近月点坐标修正为：（19.31°W，31.37°N）。
6.模型推广
此模型考虑在着陆过程中，飞行器在变质量的情况下制动的控制策略，可应用于一般的飞行器着陆控制问题；近月点计算模型可适用于一般的着陆准备轨道选取。
