2014年全国数学建模a题解析
2014高教社杯全国大学生数学建模竞赛A题论文答辩
70.9 48.8 29.9 91.3 2.588 1.056 2.498
75.7 37.4 33.3 90.8 1.838 1.168 1.702
总计
1.347 2.437 2.984 3.784 2.763
求解参数N与P的关系为
N (P 3) 3
P值太大，反而会影响计算效率，因此，取
P 30 为宜。
rpGM 1.6139 103 m / s ra a
沿运动轨迹切线方向
第2页，共15页。
1.问题一：着陆准备轨道近月点和远月点的位置
加速度为：
d 2Z dt 2
e i
d 2r dt 2
r d
dt
2
i
r
d 2
dt 2
2 dr dt
d
dt
对嫦娥三号进行受力分析，由牛顿第二定律得：
mMG ei
2014年高教社杯全国大学生数学建模竞赛
A题： 嫦娥三号软着陆轨道设计
与控制策略
第1页，共15页。
1. 问题一：嫦娥三号速度的大小和方向
vp
(1 e )
(1 e )a
(1 e )
va (1 e )a
联立上式可得近月点（近拱点），远月点（远拱点）的速度：
vp
va
raGM 1.6922 103 m / s rp a
当 rp 1752.013 103 m 时，解得 cos ，则-1 ； 180
当 ra 1837.013 103 m 时，解得 cos，则1 。 0
则在近月点的位置是 （180，1752.013 103 )
远月点的位置是 （0，1837.013 103 )
第4页，共15页。
2014数学建模国赛A题
模型一：软着陆轨道运动学与力学模型问题分析由于月球没有大气,探测器着陆时无法利用大气制动,只能利用制动发动机来减制了探测器所能携带有效载荷的质量。
探测器在月面着陆可以分为硬着陆和软着陆。
硬着陆对月速度不受限制,探测器撞上月球后设备将损坏,只能在接近月球的过程中传回月面信息;软着陆对月速度比较小,探测器着陆后可继续在月面进行考察,因此相比于硬着陆,软着陆更具有实用意义模型的建立与求解：1、建立坐标系：建立坐标系如图所示, 月心惯性坐标系m O XYZ :原点位于月球中心, m O Z 轴指向动力下降起始点, m O X 轴位于环月 轨道平面内且指向前进方向, m O Y 轴按右手定则确定, 着陆器在m O XYZ 系下的位置用极坐标(,,)r a 来表示, r 为月心到着陆器的距离矢量( r 表示大小) , A 和B 分别表示经度和纬度。
轨道坐标系Oxyz 原点位于着陆器的质心, Oz 轴为月心指向着陆器质心的方向, Ox 轴位于当地水平面内指向着陆器运动方向, Oy 轴按照右手定则确定。
制动推力F 的方向与着陆器本体轴重合, 推力方位角W 和推力仰角H 描述了制动推力F 与轨道坐标系之间的位置关系, 推力方位角W 绕正Oz 轴逆时针旋转为正, 推力仰角H 绕正Oy 轴顺时针旋转为正。
2、确定卫星在椭圆轨迹方程式：嫦娥三号在轨道上高速飞行时，设嫦娥三号卫星的质量为m ，卫星在轨道上任意点速度为v ，设月球的质量为M ;卫星与月亮之距为r ；卫星—月球系统的总能量为E ；由能量守恒定律可得：212GMm mv E r -=得：v =对于月球—卫星系统，当行星在椭圆轨迹上运动时,在卫星轨迹上有存在一点p ，月球中心和p 点的矢径为p r 该点的卫星速度为p v ，p r 与p v 之间的夹角为p β。
如果月球中心与卫星运动方向之间的垂直距离为b 。
sin p p b r β=；这个b 是卫星轨道为椭圆的短半轴；根据角动量守恒定律：sin sin p p p mvr mv r ββ=因为 sin p p b r β= 所以 sin p vr v b β= 即：sin p v b vr β=结合上述式子可以推出：sin β= （1） 式（1）为卫星椭圆轨道表达式;3、确定近月点的速度：开普勒第二定律：根据开普勒第二定律知：行星和太阳的连线在相等时间内扫过相等的面积；设近月点A 与远月点B 距离月球的距离为,A B L a c L a c =-=+，在近月点与远月点两点分别取极短的相等的时间，故有A B S S ∆=∆ 代入得：B B a c v v a c-=+ 卫星运动的总机械能等于其动能和引力势能之和，故当卫星分别经过A 、B 时的机械能为：222211()2211()22A A A A B B B B GMm GMm E mv mv L a cGMm GMm E mv mv L a c =+=--=+=--由于卫星在椭圆轨道上只受万有引力作用，所以遵循机械能守恒：A B E E =最后由椭圆方程可以求出：A B b v a c v =-=由于B v 大小为在100km 轨道上的速度，可以根据万有引力求出，即,A B v v 可以求出，其中1.73/A v km s =;它的方向为轨迹的切线方向；4、制动过程的力学与运动学分析：忽略其他星球对卫星的引力影响，则可以把嫦娥卫星的制动过程看成是一个类平抛运动；其中v 为平抛的初速度；设',F F 为万有引力和卫星的推力，S 为主减速区的竖直高度 由物理关系可以得出下列等式：'2'1()2()B BA mv F F s F F t mv l v t →→→→=++== （图）解得：728.4l km =再根据数学几何关系可以求出近月点与处在主减速区的着陆点的直线距离为746.46L km =设近着陆点与月点的坐标分别为'''(,,),(,,)x y z x y z 则可以列出下列等式L = （2）5、坐标系的转换设纬度α，经度β，海拔为h （米）月球上任意一点(,,)h αβ表示三维体系中的点(,,)x y z ，则：东经：(1737.01/1000)(cos )(sin )x h αβ≈+西经：(1737.01/1000)(cos )(sin(180))x h αβ≈++(1737.01/1000)(cos )(sin())y h αβ≈+北纬：(1737.01/1000)sin z h α≈+ 南纬：(1737.01/1000)sin(90)z h α≈++ 海拔计算时单位是米，,,x y z 单位是千米则着陆点(,,)x y z 为（(1138.08,1173.58,579.23)）近月点'''(,,)x y z 为(1752.013cos()sin(180),1752.013cos sin ,1752.013sin )αβαβα+ 综上所述，联立上述的式（1）和式（2），并且将所有已知的条件带入公式中，得到近月点的位置坐标为：'''1387.28470.281752.01x y z ⎧=⎪=⎨⎪=⎩所以容易求出远日点的位置坐标为： ''''''1452.49492.391834.37x y z ⎧=-⎪=⎨⎪=⎩；。
2014数学建模国赛A题教程
承诺书
我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参 赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下 载）。
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网 上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
在模型优化中，我们考虑了在桌面上均匀分布的力的情况，通过建立空间力系的平
衡模型，在临界条件下（桌子支撑腿受到指向桌内的摩擦力取最大值），由理论力学知
识推导出桌面上均匀分布的力 F 与 角、钢筋位置之间的函数式。计算得出桌子的稳定
性与钢筋位置无关，桌子在这种受力情况下的稳定性只与支撑腿与竖直方向的夹角有
2. 提出问题
（1）. 给定长方形平板尺寸为 120 cm × 50 cm × 3 cm，每根木条宽 2.5 cm， 连接桌腿木条的钢筋固定在桌腿最外侧木条的中心位置，折叠后桌子的高度为 53 cm。 试建立模型描述此折叠桌的动态变化过程，在此基础上给出此折叠桌的设计加工参数 （例如，桌腿木条开槽的长度等）和桌脚边缘线的数学描述。
我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展 示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。
我们参赛选择的题号是（从 A/B/C/D 中选择一项填写）：
B
我们的报名参赛队号为（8 位数字组成的编号）：
27006025
所属学校（请填写完整的全名）：
长安大学
二、 问题分析
(1).折叠桌以铰链连接，外形由直纹曲面构成。通过反复研究折叠桌的动态视频， 分析出折叠桌的运动特性，我们采用几何投影法，化三维运动为二维运动，简化模型。 同时，为了便于分析几何关系，我们仅对单组木条中最长与最短两根木条进行探究。并 通过 Solidwoks 软件绘画其几何关系图。根据各木条之间的连动原理推导出所有木条间 的关系，建立曲线参数方程表示折叠桌整体的动态变化过程。最后计算出折叠桌的设计 加工参数，并通过函数式和三维曲线图描述桌角边缘线。
2014高教社杯全国大学生数学建模竞赛(A)题目
2014高教社杯全国大学生数学建模竞赛题目（请先阅读“全国大学生数学建模竞赛论文格式规范”）A题嫦娥三号软着陆轨道设计与控制策略嫦娥三号于2013年12月2日1时30分成功发射，12月6日抵达月球轨道。
嫦娥三号在着陆准备轨道上的运行质量为 2.4t，其安装在下部的主减速发动机能够产生1500N到7500N的可调节推力，其比冲（即单位质量的推进剂产生的推力）为2940m/s，可以满足调整速度的控制要求。
在四周安装有姿态调整发动机，在给定主减速发动机的推力方向后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。
嫦娥三号的预定着陆点为19.51W，44.12N，海拔为-2641m（见附件1）。
嫦娥三号在高速飞行的情况下，要保证准确地在月球预定区域内实现软着陆，关键问题是着陆轨道与控制策略的设计。
其着陆轨道设计的基本要求：着陆准备轨道为近月点15km，远月点100km的椭圆形轨道；着陆轨道为从近月点至着陆点，其软着陆过程共分为6个阶段（见附件2），要求满足每个阶段在关键点所处的状态；尽量减少软着陆过程的燃料消耗。
根据上述的基本要求，请你们建立数学模型解决下面的问题：（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。
（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。
（3）对于你们设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。
附件1：问题的背景与参考资料；附件2：嫦娥三号着陆过程的六个阶段及其状态要求；附件3：距月面2400m处的数字高程图；附件4：距月面100m处的数字高程图。
附件1：问题A的背景与参考资料1．中新网12月12日电(记者姚培硕)根据计划，嫦娥三号将在北京时间12月14号在月球表面实施软着陆。
嫦娥三号如何实现软着陆以及能否成功成为外界关注焦点。
目前，全球仅有美国、前苏联成功实施了13次无人月球表面软着陆。
北京时间12月10日晚，嫦娥三号已经成功降轨进入预定的月面着陆准备轨道，这是嫦娥三号“落月”前最后一次轨道调整。
2013-2014年全国数模竞赛a题讲解
2013-2014年全国数模竞赛a题讲解2013-2014年全国数模竞赛A题是一道涉及建模和优化等数学概念的综合性问题。
本文将对该题进行详细的解析和讲解，帮助读者理解题目的要求，并提供一些解题思路和方法。
第一部分：理解题目该题目的题面由多个部分组成，涉及到原问题、目标、约束条件等内容。
在进行解题之前，我们首先需要完全理解题目的要求。
原问题是一个货车经过N个城市，每个城市都有相应的货物量，目标是使得货车的路径长度最短。
同时，题目要求我们设计一个数据模型，来描述这个问题。
第二部分：建立数学模型为了更好地解决问题，我们需要建立一个数学模型来描述货车的路径以及货物量的分配。
在本部分，我们将详细讲解如何建立这个模型。
假设有N个城市，每个城市的货物量分别为w1, w2, ..., wN。
我们可以将货车的路径表示为一个N*N的矩阵D，其中D[i][j]表示从第i个城市到第j个城市的距离。
同时，我们引入一个N维的向量x，其中x[i]表示从第i个城市运送的货物量。
我们的目标是最小化路径长度，即最小化下式：Minimize ∑∑D[i][j]*x[i]*x[j] (i从1到N, j从1到N)同时，我们有一些约束条件需要满足：1. 每个城市必须运送货物：∑x[i] = W，其中W是总的货物量。
2. 每个城市的货物量不能超过其容量：x[i] <= C，其中C是城市i的容量。
第三部分：优化求解在第二部分中，我们已经建立了数学模型，现在我们需要找到一种优化方法来求解这个模型。
在现实生活中，这类问题通常是NP难问题，因此我们需要采用一些启发式搜索算法。
在本部分，我们将介绍一种常用的优化方法，即遗传算法。
遗传算法模拟了自然界中的进化过程，通过不断筛选和演化来得到最优解。
遗传算法的优化步骤如下：1. 初始化种群：随机生成一组初始解，也就是一组路径和货物分配方案。
2. 评估适应度：根据路径长度和货物量是否满足约束条件，计算每个解的适应度。
2014年数学建模全国赛A题
城市土壤地质环境重金属污染摘要：针对问题一，结合数据运用MATLAB画出该区域城区分布图和8种主要重金属元素在该城区的空间分布图，接着用单因子指数得出重金属单因子污染指数，最后用内梅罗污染指数法得出重金属在功能区的污染程度由大到小依次为：工业区、交通区、生活区、公园路地区、山区。
针对问题二，根据问题一求出的单因子污染指数，用excel分别作出各重金属在各区域所占比重的饼图，并据此分析重金属污染的原因。
针对问题三、运用SPSS软件，采用因子分析法得出重金属污染物的传播特征，并在其基础上建立灰色关联分析模型，确定重金属污染源的为（2383,3692，7）、（2708,2295，22）、（4777,4891,8）、（1647,2728,6）、（18134,10046，41）、（13797,9621,18）、（21439,11383,45）、（13694,2357,33）8个位置。
针对问题四，为了更好地反映城市地质环境变化的演变模式，考虑时间、土壤的饱和度、PH、水分等因素对污染物浓度的影响，建立三维传播模型，分析演变过程。
关键词：单因子指数法、综合指数法、因子分析、灰色关联分析一、问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。
现对某城市城区土壤地质环境进行调查。
为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。
应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。
2014深圳数学建模A题
单独政策下的人口变化摘要人口的数量和结构是影响经济社会发展的重要因素。
从20世纪70年代后期以来，我国开始实施计划生育政策，该政策的实施，有效的控制了我国人口的过快增长，对经济发展和人民生活的改善做出了积极的贡献。
但另一方面，其负面影响也开始显现。
劳动人口绝对数量开始步入下降通道，人口老龄化，负担过重。
党的十八届三中全会提出了开放单独二孩，今年以来许多省、市、自治区相继出台了具体的政策。
政策出台前后各方面人士对开放“单独二孩”的效应有过大量的研究和评论。
本文通过建模来探究单独政策下的人口结构变化。
针对问题一，采用短期内精确度较高的灰色模型，进行编程求解。
以2012年以前的十年人口普查结果作为原始数据，在实行单独政策前预测2013年-2022年的人口总数变化。
得到结果，人口呈缓慢增长。
再进一步进行了为何实行单独政策的验证，对影响人口变化的各因素进行探讨研究，得到结果，我国在实施计划生育政策以来，出生率低于超低出生率，自然增长率更是停滞不前，人口老龄化加重等，研究发现造成这些问题主要原因是生育率的下降，而生育率下降又跟计划生育国策有着紧密联系。
针对问题二，预测“单独二孩”政策对未来人口结构的影响。
选择一个地区，利用该地区全国第六次人口普查数据，通过年龄推算方法得出2014-2017年分年龄的女性人数和已婚妇女数，再根据《2013年生育意愿调查数据报告》估算出符合“单独二孩”生育政策要求的潜在人群和在“单独二孩”政策下的新增出生人数，以及这些新生婴儿会什么时候出生，最后分析新增加的出生人口会导致总和生育率出现多大程度的反弹，从而研究对未来人口结构的影响。
关键词：单独二胎灰色模型年龄推算法生育率一、问题重述1.1问题背景人口的数量和结构是影响经济社会发展的重要因素。
从20世纪70年代后期以来，我国鼓励晚婚晚育，提倡一对夫妻生育一个孩子。
该政策实施30多年来，有效地控制了我国人口的过快增长，对经济发展和人民生活的改善做出了积极的贡献。
2014全国大学生数学建模竞赛A题题目及参考答案_
2014全国大学生数学建模竞赛A题题目及参考答案_ 2011高教社杯全国大学生数学建模竞赛题目,请先阅读“全国大学生数学建模竞赛论文格式规范”,A题城市表层土壤重金属污染分析随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。
现对某城市城区土壤地质环境进行调查。
为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土(0~10 厘米深度)进行取样、编号，并用GPS记录采样点的位置。
应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。
另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。
附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。
现要求你们通过数学建模来完成以下任务:(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。
(2) 通过数据分析，说明重金属污染的主要原因。
(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。
(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息,有了这些信息，如何建立模型解决问题,DJHFSJKDHFKDSJKFHSJKDFHJKDSHFDJKSFHJKDSHFJKDSHFJK题目 A题城市表层土壤重金属污染分析摘要,本文研究的是某城区警车配置及巡逻方案的制定问题，建立了求解警车巡逻方案的模型，并在满足D1的条件下给出了巡逻效果最好的方案。
2013-2014年全国数模竞赛a题讲解
2013-2014年全国数模竞赛a题讲解摘要：一、全国数模竞赛简介1.竞赛背景与历史2.竞赛级别与影响力3.对参赛者的意义与价值二、2013-2014年数模竞赛A题解析1.题目概述与背景2.题目难点与关键点3.解题思路与步骤4.答案与解析三、数模竞赛对数学教育的启示1.培养数学建模思维2.提高实际问题解决能力3.团队协作与沟通能力4.对未来数学研究的影响正文：一、全国数模竞赛简介全国数模竞赛，全名为全国大学生数学建模竞赛，是由中国数学会主办的一项面向全国大学生的数学竞赛活动。
自1992年首次举办以来，已经发展成为具有广泛影响力的国家级竞赛。
竞赛旨在激发大学生学习数学的兴趣，培养学生的创新意识和团队协作精神，提高学生解决实际问题的能力。
数模竞赛对于参赛者来说，既是一次锻炼自己的机会，也是与其他优秀学生交流学习的平台。
二、2013-2014年数模竞赛A题解析2013-2014年全国数模竞赛A题是一道具有较高难度的数学建模题目。
题目背景涉及到生物学、物理学等多个领域，要求参赛者具有较强的知识储备和综合分析能力。
在解题过程中，关键点在于如何将复杂问题抽象为数学模型，并运用合适的数学方法求解。
通过分析题目，我们可以将问题划分为以下几个部分：1.题目概述与背景：题目描述了一种生物学现象，要求参赛者基于这一现象建立数学模型，并分析其动力学性质。
2.题目难点与关键点：难点主要在于如何将生物学现象抽象为数学模型，以及如何运用数学方法分析模型的动力学性质。
解决这一问题的关键在于对题目背景知识的掌握和对数学建模方法的理解。
3.解题思路与步骤：首先，需要深入理解题目背景，提取关键信息；其次，根据题目要求建立数学模型；最后，运用数学方法分析模型的性质。
4.答案与解析：根据解题思路，参赛者可以得到最终答案，并对答案进行解析，阐述答案的合理性和正确性。
三、数模竞赛对数学教育的启示全国数模竞赛对于数学教育具有重要的启示作用。
首先，竞赛培养了学生的数学建模思维，使他们能够将现实问题抽象为数学问题，运用数学方法解决实际问题。
2014高教社杯全国大学生数学建模竞赛A题_共26页
2014 高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号（由赛区组委会评阅前进行编号）：
赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注
全国统一编号（由赛区组委会送交全国前编号）：
全国评阅编号（由全国组委会评阅前进行编号）：
嫦娥三号软着陆轨道设计与控制策略 摘要
本文针对嫦娥三号软着陆轨道设计与控制策略问题，通过提取题目中的信 息，利用拱点的概念、B 样条函数逼近的统计定位方法、非线性规划问题及哈 密尔顿函数为理论基础进行了完整的建模工作。首先，通过建立坐标系结合物 理学运动公式求解出了近月点与远月点的位置及相应的速度；在此基础上，利 用 B 样条函数逼近的方法确定了嫦娥三号的着陆轨；最后通过分解着陆过程并 利用非线性规划问题及哈密尔顿函数确定着陆阶段的最优控制策。
参赛队员 (打印并签名) ：1.
2.
3.
指导教师或指导教师组负责人 (打印并签名)：
（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上
内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖
资格。）
日期： 2014 年 9 月 15 日
赛区评阅编号（由赛区组委会评阅前进行编号）：
我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开 展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。
我们参赛选择的题号是（从 A/B/C/D 中选择一项填写）：
A
我们的报名参赛队号为（8 位数字组成的编号）：
07033001
所属学校（请填写完整的全名）：吉林师范大学博达学院
针对问题二，采用 B 样条函数逼近的运动学统计定位方法确定了在着陆弧 段上任意时刻的位置方程，从而刻画出了嫦娥三号的着陆轨道，并用 matlab 对轨 迹进行了模拟。在 6 个阶段的最优控制策略上，先通过直角坐标系得出质心的运 动方程，再通过对 6 个阶段初始条件和终端状态的分解，利用非线性规划问题 求解哈密尔顿函数，得出性能指标（耗燃量）的最小值为：382.6531kg，从而确 定了最优控制策略。
2014数学建模A题资料
（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。
（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。
（3）对设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。
二、
本问题的一个关键点在于求解出嫦娥三号卫星的绕月运行轨道。由于已知卫星的运行轨迹为近月点15km，远月点100km的椭圆形轨道。
嫦娥三号软着陆轨道设计与控制策略
摘要：
嫦娥三号卫星采用的是软着陆方式登陆月球，在卫星高速飞行的情况下，我们要精确地在月球预定区域内实现软着陆，需要对其运行轨道进行设计并制定相应控制策略。由于天体的运动均满足开普勒三大定律以及总能量守恒定律，我们据此建立一系列的方程，最终求得卫星在近月点处的速度大小
求得卫星在远月点处的速度大小 。其速度方向均为当前运动轨道的切方向。
我们将嫦娥三号卫星绕月运行的椭圆轨道抽象出如下图所示的一个简单的几何图形，月球的月心位于椭圆的一个焦点F上，椭圆的半长轴为 ，半短轴为 ，半焦距为 。A点为近月点，速度为 ，B点为远月点，速度为 。易知A、B两点距月心的距离 ， 。在一个极小的时间段 内，卫星与月心连线扫过的面积分别为 ， 。由开普勒第二定律可知，卫星与中心天体连线在单位时间内扫过的面积相等，所以 ，代入化简后可得公式(1) ；由于卫星运动的总机械能等于其动能和引力势能之和，所以在A点，卫星的总机械能 （公式2），同理B点的总机械能 （公式3）。卫星在运行过程中只有动能和引力势能之间的转化，机械能守恒，所以 （公式4）。[1]
对于第二问，要设计一种方案使得在软着陆的过程中耗能最少，并达到预设的各项指标。我们需要不断的去设计、计算、调整，在不断的尝试摸索中寻找出一个比较不错的软着陆方案。
2014研究生数学建模竞赛A题
2014年全国研究生数学建模竞赛A题小鼠视觉感受区电位信号(LFP)与视觉刺激之间的关系研究人类脑计划（Human Brain Project, HBP）是继人类基因组计划之后，又一国际性科研计划，其核心是神经信息学(Neuroinformatics）。
该研究旨在努力探究数十亿个神经元的信息，以期对知觉、行动以及意识等有更进一步的了解。
科学家们预期这是一条开发新技术的好途径，由此可能进一步认识像老年痴呆和帕金森综合症等疾病，有望为各种精神疾病研究出新的治疗方法。
此外，该计划还可以更好地为人工智能服务。
目前该计划已经取得一定进展，例如在2014年足球世界杯巴西开幕式上，脊髓损伤患者开球就是利用了美国杜克大学的神经生物学家、大脑-计算机界面研究的先驱、巴西人Miguel Nicolelis开发的一种由大脑脑电波控制的外骨骼系统。
它既是人类脑计划的结晶，也是2014年Nature十大科学展望领域之一神经科学的代表性成就。
1.脑电波介绍脑是支配人和高级动物活动的司令部和信息中心，神经系统承担着感受外界刺激，产生、处理、传导和整合信号，实现各种认知活动（如知觉、学习、记忆、情绪、语言、意思和思维等），以及运动控制等众多功能。
神经系统的基本结构单元是神经元，其放电活动涉及复杂的物理化学过程，表现出丰富的非线性动力学行为。
神经系统整体可视为由数目众多的神经元组成的庞大而复杂的信息网络，通过对信息的处理、编码、整合，转变为传出冲动，从而联络和调节机体的各系统和器官的功能。
神经元对信息的处理和加工是神经元集群共同完成的，而神经元集群的同步形成较强的电信号就是脑电波。
脑电波（Electroencephalogram,EEG）是大脑在活动时，脑皮质细胞群之间形成电位差，从而在大脑皮质的细胞外产生电流。
它们是脑神经细胞的电生理活动在大脑皮层或头皮表面的总体反映，将大脑活动时这种电波的变化记录下来就得到脑电图。
2.脑电波的应用脑电波或脑电图是一种比较敏感的客观指标，不仅用于脑科学的基础理论研究，而且更重要的在于临床实践的应用。
2014国赛a题讲评
2014国赛a题讲评摘要：一、2014 国赛A 题概述二、试题解析1.题目要求2.考察重点3.解题思路三、参赛作品亮点分析1.创新思维2.团队协作3.实际应用价值四、讲评总结正文：一、2014 国赛A 题概述2014 年全国大学生数学建模竞赛（简称“国赛”）A 题主要考察参赛选手的数学建模能力、解决实际问题的能力和创新思维。
题目具有一定的难度，要求参赛选手在规定时间内完成对题目的解析、建模和求解。
二、试题解析1.题目要求2014 国赛A 题要求参赛选手针对某一具体问题进行数学建模，并完成模型的求解和验证。
题目要求参赛选手在论文中详细阐述所建立的数学模型，并给出求解过程和结果。
2.考察重点此题主要考察参赛选手的数学建模能力、解决实际问题的能力和创新思维。
要求参赛选手灵活运用所学的数学知识，对题目进行深入分析，并建立合适的数学模型。
3.解题思路针对这类题目，解题思路通常包括以下几个步骤：（1）认真阅读题目，理解题意，明确题目要求。
（2）对题目进行分析，找出问题的关键点，确定建模方向。
（3）建立数学模型，根据题目要求选择合适的数学方法。
（4）求解模型，得出结果，并对结果进行分析和讨论。
（5）撰写论文，详细介绍模型的建立和求解过程，阐述结果的意义和影响。
三、参赛作品亮点分析1.创新思维在2014 国赛A 题的参赛作品中，许多队伍都展现出了很强的创新思维。
他们能够灵活运用所学的数学知识，对题目进行深入分析，并建立合适的数学模型。
在求解过程中，他们也能够尝试不同的方法，寻找最优解。
2.团队协作在此次竞赛中，很多队伍都表现出了很好的团队协作精神。
他们能够充分发挥团队成员的优势，共同分析问题，共同完成模型的建立和求解。
这种团队协作精神对于解决复杂问题具有重要意义。
3.实际应用价值许多参赛作品都具有一定的实际应用价值。
他们所建立的数学模型不仅可以解决竞赛题目，还能够应用于实际生活中的问题。
这充分体现了数学建模竞赛的宗旨，即通过解决实际问题，推动数学知识的应用和发展。
2014全国大学生数学建模竞赛A题论文解析
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则．我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题．我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料)，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出．我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性．如有违反竞赛规则的行为，将受到严肃处理．我们参赛选择的题号是（从A/B/C/D中选择一项填写）赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：嫦娥三号软着陆轨道设计与控制策略摘要本文针对嫦娥三号软着陆轨道设计与控制策略的实际问题，以理论力学（万有引力、开普勒定律、万能守恒定律等）和卫星力学知识为理论基础，结合微分方程和微元法，借助MATLAB软件解决了题目所要求解的问题。
针对问题（1），在合理的假设基础上，利用物理理论知识、解析几何知识和微元法，分析并求解出近月点和远月点的位置，即139.1097 。
再运用能量守恒定律和相关数据，计算出速度v（近月点的速度）1=1750.78/v（远月点的速度）=1669.77/m s，，最后利用曲线的切线m s,2方程，代入点（近月点与远月点）的坐标求值，计算出方向余弦即为相应的速度方向。
针对问题（2）关键词：模糊评判，聚类分析，流体交通量，排队论，多元非线性回归一、问题重述嫦娥三号于2013年12月2日1时30分成功发射，12月6日抵达月球轨道。
嫦娥三号在着陆准备轨道上的运行质量为2.4t，其安装在下部的主减速发动机能够产生1500N到7500N的可调节推力，其比冲（即单位质量的推进剂产生的推力）为2940m/s，可以满足调整速度的控制要求。
2014全国数学建模A题一等奖论文
v2 = 526.94m / s 。即远月点的速度为 526.94 m / s .方向为水平方向。
图 3 着陆准备轨道环绕模拟
由于近月点和远月点分别是椭圆轨道的长半轴的两个端点， 且两点的连线经过月心 (图 3)，因此由对称性可知远月点的位置为(19.51E，32.31S)，高度为 100000 米。
360 = 30.301千米 / 度 2πR P 为纬度改变 1o 水平距离的改变量。 p=
根据能量守恒定律可知：
1 2 1 2 mv1 + mg ′h近 = mv 2 + mg ′h远 2 2 其中： v1 为近月点的速度； v2 为远月点的速度。
⑵模型的求解 在本题中由于我们无法确定任意时刻减速动力以及速度的大小及方向， 因此我们通 过假设简化模型，从而对问题进行求解。由于发动机推力主要是用于减少飞行器的横向 速度，同时克服由月球引力引起的径向速度，我们假设了嫦娥三号可以通过自身调节机 制使得自己在运动过程中竖直方向受恒力作用，方向向下，水平方向也受恒力作用，方 向与水平速度方向相反，初速度为 1700m/s。 因此我们可以将抛物线下降的过程分解成竖直方向匀加速，水平方向匀减速的运 动。（如图 1）由附件 2 可知，嫦娥三号在 3000m 时已经基本位于目标上方，所以我们 认为在 3000 米处水平速度近似为 0，57 m / s 为其竖直方向速度。
§3 模型的假设
1．由给出的附件月球的形状扁率为1/963.7256，数量级较小，假设月球为一个球体。 2．由于从近月点100km左右的高度降落到地球表面的时间比较对短，假设嫦娥三号不受 非球项、日月引力摄动等影响因素的影响。 3．假设月球引力场为平行定常引力场，嫦娥三号着陆轨道不受月球自转的影响。 4．假设月球表面海拔为零的球面势能为0。 5．假设嫦娥三号水平移动的距离近似为着陆划过月球表面弧度长度。 6．假设月球的重力加速度恒定，为 1 / 6 g 。
2014年全国数学建模大赛A题
2014高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。
如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的报名参赛队号为（8位数字组成的编号）：25001113所属学校（请填写完整的全名）：云南大学参赛队员(打印并签名) ：1. 林博文2. 张竞文3. 方春晖指导教师或指导教师组负责人(打印并签名)：李海燕（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。
以上内容请仔细核对，提交后将不再允许做任何修改。
如填写错误，论文可能被取消评奖资格。
）日期：2014年9月15日赛区评阅编号（由赛区组委会评阅前进行编号）：2014高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：嫦娥三号软着陆轨道设计与控制策略优化摘 要 嫦娥三号是中国国家航天局嫦娥工程第二阶段的登月探测器，包括着陆器和玉兔号月球车。
嫦娥三号在高速飞行的情况下，要保证准确地在月球预定区域内实现软着陆，关键问题是着陆轨道与控制策略的设计。
2014全国大学生数学建模a题
2014高教社杯全国大学生数学建模竞赛a题摘要2013年嫦娥三号成功发射，标志着我国航天事业上的又一个里程碑，针对嫦娥三号软着陆问题，分别建立着陆前轨道准备模型和软着陆轨道模型，建立动力学方程，以燃料最省为目标进行求解。
问题一：在软着陆前准备轨道上利用开普勒定律、能量守恒定律以及卫星轨道的相关知识，利用牛顿迭代法分别确定了近月点和远月点的速度分别为 1.6925km/s、1.6142km/s，位置分别为（19.91W，20.96N），（160.49E,69.31S）。
问题二：在较为复杂的软着陆阶段，因为相对于月球的半径，嫦娥三号到月球的表面的距离太小，如果以月球中心建立坐标系会造成比较大的误差，因此选择在月球表面建立直角坐标系，在主减速阶段的类平抛面上建立相应的动力学模型，求出关键点的状态和并设计出相应的轨道，接下来通过利用灰度值阀值分割方法和螺旋搜索法对粗避障阶段和精避障阶段的地面地形进行相应的分析，找出安全点，然后调整嫦娥三号的方向以便安全降落，最后在落地时通过姿态发动机调整探测器的姿态，使之可以平稳的落到安全点上，在以上的各个阶段都可以以燃料最省为最优指标，从而建立非线性的最优规划的动力学模型，并基于该动力学模型可以对各个阶段的制导率进行优化设计由此就可以得到各个阶段的最优控制策略，问题三：最后针对所设计的轨道和各个阶段的控制策略进行了误差分析和灵敏度分析。
对系统误差和偶然误差都做了解释；通过灵敏度分析发现，嫦娥三号在近月点的位置对结果的影响最大。
关键字牛顿迭代法，灰度值阀值分割，螺旋搜索法，灵敏度分析一、问题重述嫦娥三号于2013年12月2日1时30分成功发射，12月6日抵达月球轨道。
嫦娥三号在着陆准备轨道上的运行质量为 2.4t，其安装在下部的主减速发动机能够产生1500N到7500N的可调节推力，其比冲（即单位质量的推进剂产生的推力）为2940m/s，可以满足调整速度的控制要求。
在四周安装有姿态调整发动机，在给定主减速发动机的推力方向后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。
2014年数学建模题目
2014年数学建模题目
2014年全国研究生数学建模竞赛A题是关于小鼠视觉感受区电位信号（LFP）与视觉刺激之间的关系研究。
题目要求探究数十亿个神经元的信息，以期对知觉、行动以及意识等有更进一步的了解，并可能为各种精神疾病研究出新的治疗方法。
具体而言，需要解决的问题包括：
1. 脑电波介绍：脑是支配人和高级动物活动的司令部和信息中心，神经系统承担着感受外界刺激，产生、处理、传导和整合信号，实现各种认知活动（如知觉、学习、记忆、情绪、语言、意思和思维等），以及运动控制等众多功能。
2. 脑电波的应用：脑电波可以用来研究大脑的活动和功能，包括神经元的信息传递和整合。
脑电波的应用范围广泛，包括神经科学、心理学、生理学、医学等领域。
对于这个题目，需要运用数学建模的方法，结合脑电波数据和视觉刺激数据，分析小鼠视觉感受区的电位信号与视觉刺激之间的关系，探究大脑的信息处理机制。
具体而言，需要解决的问题包括：
1. 数据预处理：对脑电波数据进行预处理，包括滤波、去噪、归一化等操作，以提高数据的质量和可靠性。
2. 特征提取：从脑电波数据中提取出与视觉刺激相关的特征，包括时域和频域特征。
3. 模型建立与优化：根据提取的特征，建立数学模型，并不断优化模型参数和结构，以提高预测准确率和稳定性。
4. 结果解释与讨论：对模型的结果进行解释和讨论，探究小鼠视觉感受区的电位信号与视觉刺激之间的关系，以及大脑的信息处理机制。
总之，2014年全国研究生数学建模竞赛A题是一个具有挑战性和重要意义的题目，需要运用数学建模的方法，结合脑电波数据和视觉刺激数据，分析小鼠视觉感受区的电位信号与视觉刺激之间的关系，探究大脑的信息处理机制。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。
如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。
（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。
以上内容请仔细核对，提交后将不再允许做任何修改。
如填写错误，论文可能被取消评奖资格。
）赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：嫦娥三号软着陆轨道设计与控制策略摘要嫦娥三号卫星着陆器实现了我国首次地外天体软着陆任务。
要保证准确的在月球预定区域内实现软着陆轨道与控制策略的设计。
问题一运用活力公式[1]来建立速度模型，利用matlab软件代入数值计算出。
所求速度33⨯⨯(=1.692210m/s,=1.613910m/s)v v远近采用轨道六根数[2]来建立近月点，远月点位置的模型。
轨道根数是六个确定椭圆轨道的物理量，也是联系赤道直角坐标与轨道极坐标重要夹角的关系。
通过着陆点的位置求出轨道根数各个值的数据，从而确定近月点，远月点的位置，坐标分别为(19.51W 27.88N 15KM),(160.49 27.885S 100KM)E。
问题二“嫦娥三号”软着陆过程中需要经历6个不同的阶段,对于主减速阶段，在极坐标系下建立其运动方程。
结合Pontryagin极大值原理[3]和哈密顿函数[4]，化简出燃料最省的软着陆轨道方程，得出最优控制变量的变化规律。
对于其它各阶段，将其简化为加速度不同的线性运动模型，利用动能定理得出相应轨道方程和控制策略。
问题三对第二问中求出的“嫦娥三号”推力和速度切线方向夹角ϕ,给ϕ增加或减小一个角度ϕ,分别求出各个对应的近月点坐标'y。
之后求各个坐标与其原始值之间的变化量'y并求其平均值'y，得到其敏感性因数，敏感性系数越大，说明该属性对模型的影响越大。
关键字：活力公式轨道六根数 Pontryagin极大值原理燃料最省一、问题重述嫦娥三号于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）对于你们设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。
对于误差因数分析,通过计算着陆轨道与策略的理论值与实际值之间的变化量,并求其平均值,得出平均值与实际值的比值,其中比值越大说明其误差越大,越不可行。
二、模型假设1.月球可看做一质量均匀、形状标准的球体；2.反向推力大小为常定值；3.飞行器为一质点，不考虑飞行器的姿态对轨道的影响，也不考虑飞行器姿态；4.忽略重力而只考虑空气阻力的作用；忽略地球曲率的影响，在在入轨道是直线轨道；5.不考虑地球等其他天体的影响；三、符号说明G------万有引力常量；M------月球的质量；F ------发动机的推力；ϕ------推力的方向角，即推力和切向速度的夹角；r ------嫦娥三号卫星的极半径；θ------极角；r v ------径向速度； v θ------切向速度；m ------ 任意时刻嫦娥三号卫星的质量；m ------发动机单位时间消耗的燃料质量；m ------嫦娥三号在着陆轨道上的质量。
四、模型的分析、建立与求解4.1问题一的建立与求解 4.1.1近月点，远月点的速度近月点，远月点均在椭圆轨道上，建立以月心为原点，椭圆轨道长半轴为x 轴，短半轴为y 轴的平面直角坐标系。
运用活力公式建立速度模型并求解数值。
活力公式，又叫轨道能量.这个公式是二体问题的一个积分。
是反映了天体的位置、速度和轨道半长径之间的相互关系。
平面运动的面积定律: 二体问题中作用于“嫦娥三号”卫星上的力总是指向地心，结果是轨道是是始终保持在固定平面上。
因为力总是与位置矢量相反，没有垂直于轨道平面上的加速度，所以卫星不可能脱离轨道平面。
卫星加速度••r 可由牛顿万有引力得出：3GMr ••=-r r(1)作为这一事实的数学描述，式（1）两边叉乘位置矢量r ，则3+r GM••⨯=-r r r r （）=0 （2）上面方程右边为0，因为一个矢量本身叉乘为0，方程左边可展开为td=+=d •••••••⨯⨯⨯⨯r r r r r r r r （）（3）v •=r因为•⨯r r 对时间的导数等于0，因此•⨯r r 本身必须为常数，也就是： =st con •⨯=r r h （4）两个矢量叉乘所产生的矢量几何上垂直于这两个矢量。
因此，位置矢量r 和速度矢量•r 总是垂直于h ，换句话说，运行轨道在一个平面。
矢量h 为单位质量的角动量或者说是特殊角动量，它和角动量l 关联，有=l mh ，其中m 是卫星质量。
给（1）式两边叉乘矢量h ，可以发现轨道的其他特性:()GM r ••⨯=-rh r (5)关于开普勒运动的能量积分定律，它涉及卫星和地心距的关系。
为此，将式(5)两边平方，得22222222+2(12cos )(6)=(2(12cos )(1))GM GMrGM e e GM e e •⨯=- =-++ -+--rh r v h v v （）（）（）（） 因为矢量h 和•r 互相垂直，所以上式左边的值22h v ,其中表示卫星速度。
代入半长轴的导数221(1)GM e a h -=,利用圆锥截面方程，得任意开普勒轨道(椭圆曲线轨道),活力公式的表达式为221=()()v G M m r a +- (7)在此，因为卫星的质量相对于月球的质量来说太小，我们计算时忽略卫星的质量，得到简化的活力公式表达式：221()v GM r a =- (8) v ------表示两天体间的相对速度r ------表示两天体间的相对距离a ------表示半长轴(椭圆:0a >;抛物线:a =∞或10a =;双曲线:a <∞)G ------表示万有引力常数M ,m ------表示两天体的质量在matlab环境下，编程求解速度分别为3=1.692210v⨯近，3=1.613910v⨯远，速度方向为轨道切线方向。
4.1.2 近月点，远月点的位置（1）轨道根数:轨道根数[1](或称轨道要素或轨道参数)是对选定的两个质点,在牛顿运动定律和平方反比定律的重力吸引下,确认特定轨道所必须要的参数。
1.轨道半长轴a:既为平均轨道半径,但是不是长轴与短轴的算术平均数。
2.轨道偏心率e:为椭圆扁平程度的一种量度,定义是椭圆两焦点间的距离与长轴长度的比值就是cea=。
3.轨道倾角i:行星轨道面对黄道面的倾角或在升交点处从黄道面逆时针方向量到行星轨道的角度。
4.升交点黄道经度Ω:行星轨道升交点的黄道经度。
5.近月点幅角ω:从升交点沿行星运功轨道逆时针量到近日点的角度。
6.指定历元的平近点角0M:行星对应0t时刻的平近点角在使用以上的轨道根数,可找出天体按开普勒轨道(即二体问题中的轨道)运行位置,但在实际问题中,若天体所受的其他作用力不可忽略,便需加这些摄动(因素)项来修正其位置（2）建立坐标系建立月心赤道坐标系，它与月球自转轴和赤道方向对齐。
原点是月心，z轴是指向北极，赤道平面组成了x—y参考平面。
月球的自转和公转是一样的时间，所以就只能看见一面，所以x轴指向月球始终面对地球的那面中心点。
如图所示：图一赤道坐标系中心点的位置可以通过三维直角坐标(,,)x y z 或者极坐标（,,αδγ）来表示。
两种坐标转换如下：cos cos cos sin sin x y r z δαδαδ ⎛⎫⎛⎫ ⎪ ⎪== ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭r （9）图二在轨道根数这六个物理量里，a 和 e 确定了轨道形状，M 确定了沿轨道的位置， ωi Ω这三个根数则是确定了轨道在空间的定向，即就是与赤道直角坐标的角度的联系。
（3）建模轨道根数的计算模型：由（4）式可知角动量矢量：(10)y z z y z x x z x y y x •••••••⎛⎫- ⎪ ⎪=⨯=- ⎪ ⎪- ⎪⎝⎭h r r和角动量的模||h =h 。
/sin sin sin cos /(11)cos /x x y y z z h h w i i h h w i h h w ⎛⎫⎛⎫++Ω⎛⎫⎪ ⎪ ⎪Ω=-=- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪++⎝⎭⎝⎭⎝⎭式中/h =w h 作为i 和Ω的函数，i 和Ω由（9）式得到。
因此，倾角和交升点赤经可得到如下公式：arctan z i =⎝⎭ (12)cos /z i h h=arctan()arctan()y x xyh hh hΩ==--w w （13）由瞬时角动量可导出半通径与椭圆的基本特性半通径式相等：22(1)h p a e GM ==- （14）有（14）式可以推出轨道偏心率e =（15） 由轨道极坐标方程根据椭圆基本特性可以得到椭圆上各点的极径r 与真近点角θ公式：1cos pr e θ=- （16）·0,[0,180]0,[180,360]·r v r v θθ⎧>⎪ ⎨<⎪⎩。