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 ⎧=-⎪=⎨⎪=⎩；。

2013年“认证杯”数学中国数学建模网络挑战赛
第一阶段
A题护岸框架
在江河中，堤岸、江心洲的迎水区域被水流长期冲刷侵蚀。

在河道整治工程中，需要在受侵蚀严重的部位设置一些人工设施，以减弱水流的冲刷，促进该处泥沙的淤积，以保护河岸形态的稳定。

现在常用的设施包括四面六边透水框架1等。

这是一种由钢筋混凝土框杆相互焊接而成的正四面体结构，常见的尺寸为边长约1m，框杆截面约0.1×0.1m。

将一定数量的框架投入水中，在水中形成框架群，可以使水流消能减速，达到减弱冲击，防冲促淤的效果。

对四面六边透水框架群，框架尺寸、架空率和长度2都直接或间接地影响着其消能减速的效果。

当前人们通过水槽实验等方法，已经积累了一些实验数据，也见诸各类文献当中。

但由于数据量和数据精度仍有各种不足之处，还没有形成足够完善的经验公式。

请你建立合理的数学模型，给出这三个参数与其减速效果之间的关系。

1“护岸防洪防冲四面六边透水框架”是水利部西北水利科学研究所拥有的实用新型专利，
专利号98233050。

2此处指顺水流方向的铺设长度。

1
2。

城市土壤地质环境重金属污染摘要：针对问题一，结合数据运用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记录采样点的位置。

应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。

单独政策下的人口变化摘要人口的数量和结构是影响经济社会发展的重要因素。

从20世纪70年代后期以来，我国开始实施计划生育政策，该政策的实施，有效的控制了我国人口的过快增长，对经济发展和人民生活的改善做出了积极的贡献。

但另一方面，其负面影响也开始显现。

劳动人口绝对数量开始步入下降通道，人口老龄化，负担过重。

党的十八届三中全会提出了开放单独二孩，今年以来许多省、市、自治区相继出台了具体的政策。

政策出台前后各方面人士对开放“单独二孩”的效应有过大量的研究和评论。

本文通过建模来探究单独政策下的人口结构变化。

针对问题一，采用短期内精确度较高的灰色模型，进行编程求解。

以2012年以前的十年人口普查结果作为原始数据，在实行单独政策前预测2013年-2022年的人口总数变化。

得到结果，人口呈缓慢增长。

再进一步进行了为何实行单独政策的验证，对影响人口变化的各因素进行探讨研究，得到结果，我国在实施计划生育政策以来，出生率低于超低出生率，自然增长率更是停滞不前，人口老龄化加重等，研究发现造成这些问题主要原因是生育率的下降，而生育率下降又跟计划生育国策有着紧密联系。

针对问题二，预测“单独二孩”政策对未来人口结构的影响。

选择一个地区，利用该地区全国第六次人口普查数据，通过年龄推算方法得出2014-2017年分年龄的女性人数和已婚妇女数，再根据《2013年生育意愿调查数据报告》估算出符合“单独二孩”生育政策要求的潜在人群和在“单独二孩”政策下的新增出生人数，以及这些新生婴儿会什么时候出生，最后分析新增加的出生人口会导致总和生育率出现多大程度的反弹，从而研究对未来人口结构的影响。

关键词：单独二胎灰色模型年龄推算法生育率一、问题重述1.1问题背景人口的数量和结构是影响经济社会发展的重要因素。

从20世纪70年代后期以来，我国鼓励晚婚晚育，提倡一对夫妻生育一个孩子。

该政策实施30多年来，有效地控制了我国人口的过快增长，对经济发展和人民生活的改善做出了积极的贡献。

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的条件下给出了巡逻效果最好的方案。

2014年数学建模模拟训练4【A题】举世瞩目的2014年世界杯决赛阶段的比赛2014年6月12日至7月13日在南美洲国家巴西举行。

巴西世界杯共有32支球队参赛。

除去东道主巴西自动获得参赛资格以外，其他31个国家需通过参加2011年6月开始的预选赛获得参赛资格。

巴西世界杯期间，总共在巴西境内举办共计64场比赛角逐出冠军。

假如你是中国体育彩票中心研究员，请根据赛制、赛程安排、分组形势及各自的实力，请建立数学模型进行分析，并给出：1.中国体育彩票中心设计若干世界杯竞猜游戏，并分析各种奖项出现的可能性，奖项和奖金额设置对彩民的吸引力等各因素评价游戏的合理性。

例如：给出本次世界杯32强的各级（32进16，16进8，8进4，4进2，夺冠）赔率。

2.给足球彩民写一篇短文，供买彩票参考。

【B题】众所周知，吸烟不仅危害自身健康，而且由此引起的被动吸烟更是危害公众身心健康的主要原因。

为此，如何帮助相关人士摆脱烟瘾的困扰也就成为一个重要的研究课题。

本文研究数据涉及234人，他们都自愿表示戒烟但还未戒烟。

在他们戒烟的这一天，测量了每个人的CO（一氧化碳）水平并记下他们抽最后一支烟到CO 测定时间.。

CO的水平提供了一个他们先前抽烟数量的客观指标，但其值也受到抽最后一支烟的时间的影响, 因此抽最后一支烟的时间可以用来调整CO的水平。

记录下研究对象的性别、年龄及自述每日抽烟支数。

这个调查跟踪1年, 考察他们一直保持戒烟的天数, 由此估计这些人中再次吸烟的累加发病率, 也就是原吸烟者戒烟一段时间后又再吸烟的比例. 其中假设原烟民戒烟的可信度是很低的（更恰当地说多数是再犯者）戒烟天数是从0到他（她）退出戒烟或研究截止时间（1 年）的天数。

假定他们全部没有人中途退出研究。

请回答下列问题：1）试分析上述234人中再次吸烟的累加发病率分布情况（如不同年龄段、不同性别等因素下的累加发病率分布情况）。

2）你认为年龄、性别、每日抽烟支数及调整的CO浓度等因素会影响戒烟时间（天数）长短吗？如果影响请利用附录中的数据，分别给出戒烟时间与上述你认为有影响的因素之间的定量分析结果。

2014年全国研究生数学建模竞赛A题小鼠视觉感受区电位信号(LFP)与视觉刺激之间的关系研究人类脑计划（Human Brain Project, HBP）是继人类基因组计划之后，又一国际性科研计划，其核心是神经信息学(Neuroinformatics）。

该研究旨在努力探究数十亿个神经元的信息，以期对知觉、行动以及意识等有更进一步的了解。

科学家们预期这是一条开发新技术的好途径，由此可能进一步认识像老年痴呆和帕金森综合症等疾病，有望为各种精神疾病研究出新的治疗方法。

此外，该计划还可以更好地为人工智能服务。

目前该计划已经取得一定进展，例如在2014年足球世界杯巴西开幕式上，脊髓损伤患者开球就是利用了美国杜克大学的神经生物学家、大脑-计算机界面研究的先驱、巴西人Miguel Nicolelis开发的一种由大脑脑电波控制的外骨骼系统。

它既是人类脑计划的结晶，也是2014年Nature十大科学展望领域之一神经科学的代表性成就。

1.脑电波介绍脑是支配人和高级动物活动的司令部和信息中心，神经系统承担着感受外界刺激，产生、处理、传导和整合信号，实现各种认知活动（如知觉、学习、记忆、情绪、语言、意思和思维等），以及运动控制等众多功能。

神经系统的基本结构单元是神经元，其放电活动涉及复杂的物理化学过程，表现出丰富的非线性动力学行为。

神经系统整体可视为由数目众多的神经元组成的庞大而复杂的信息网络，通过对信息的处理、编码、整合，转变为传出冲动，从而联络和调节机体的各系统和器官的功能。

神经元对信息的处理和加工是神经元集群共同完成的，而神经元集群的同步形成较强的电信号就是脑电波。

脑电波（Electroencephalogram,EEG）是大脑在活动时，脑皮质细胞群之间形成电位差，从而在大脑皮质的细胞外产生电流。

它们是脑神经细胞的电生理活动在大脑皮层或头皮表面的总体反映，将大脑活动时这种电波的变化记录下来就得到脑电图。

2.脑电波的应用脑电波或脑电图是一种比较敏感的客观指标，不仅用于脑科学的基础理论研究，而且更重要的在于临床实践的应用。

承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则．我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题．我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料)，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出．我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性．如有违反竞赛规则的行为，将受到严肃处理．我们参赛选择的题号是（从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，可以满足调整速度的控制要求。

承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则．我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题．我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料)，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出．我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性．如有违反竞赛规则的行为，将受到严肃处理．我们参赛选择的题号是（从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，可以满足调整速度的控制要求。

clfclear all%build the GUI%define the plot buttonplotbutton=uicontrol('style','pushbutton',...'string','Run', ...'fontsize',12, ...'position',[100,400,50,20], ...'callback', 'run=1;');%define the stop buttonerasebutton=uicontrol('style','pushbutton',...'string','Stop', ...'fontsize',12, ...'position',[200,400,50,20], ...'callback','freeze=1;');%define the Quit buttonquitbutton=uicontrol('style','pushbutton',...'string','Quit', ...'fontsize',12, ...'position',[300,400,50,20], ...'callback','stop=1;close;');number = uicontrol('style','text', ...'string','1', ...'fontsize',12, ...'position',[20,400,50,20]);%CA setupn=100;%数据初始化z=zeros(1,n);%元胞个数z=roadstart(z,5);%道路状态初始化，路段上随机分布5辆cells=z;vmax=3;%最大速度v=speedstart(cells,vmax);%速度初始化x=1;%记录速度和车辆位置memor_cells=zeros(3600,n);memor_v=zeros(3600,n);imh=imshow(cells);%初始化图像白色有车，黑色空元胞set(imh, 'erasemode', 'none')axis equalaxis tightstop=0; %wait for a quit button pushrun=0; %wait for a drawfreeze=0; %wait for a freeze（冻结）while (stop==0)if(run==1)%边界条件处理，搜素首末车，控制进出，使用开口条件a=searchleadcar(cells);b=searchlastcar(cells);[cells,v]=border_control(cells,a,b,v,vmax); i=searchleadcar(cells);%搜索首车位置for j=1:iif i-j+1==n[z,v]=leadcarupdate(z,v);continue;else%======================================加速、减速、随机慢化 if cells(i-j+1)==0;%判断当前位置是否非空continue;else v(i-j+1)=min(v(i-j+1)+1,vmax);%加速%=================================减速k=searchfrontcar((i-j+1),cells);%搜素前方首个非空元胞位置if k==0;%确定于前车之间的元胞数d=n-(i-j+1);else d=k-(i-j+1)-1;endv(i-j+1)=min(v(i-j+1),d);%==============================%减速%随机慢化v(i-j+1)=randslow(v(i-j+1));new_v=v(i-j+1);%======================================加速、减速、随机慢化%更新车辆位置z(i-j+1)=0;z(i-j+1+new_v)=1;%更新速度v(i-j+1)=0;v(i-j+1+new_v)=new_v;endendendcells=z;memor_cells(x,:)=cells;%记录速度和车辆位置memor_v(x,:)=v;x=x+1;set(imh,'cdata',cells)%更新图像%update the step number diaplaypause(0.2);stepnumber = 1 + str2num(get(number,'string'));set(number,'string',num2str(stepnumber))endif (freeze==1)run = 0;freeze = 0;enddrawnowend///////////////////////////////////////////////////////////////////////Function[new_matrix_cells,new_v]=border_control(matrix_cells,a,b,v,vmax) %边界条件，开口边界，控制车辆出入%出口边界，若头车在道路边界，则以一定该路0.9离去n=length(matrix_cells);if a==nrand('state',sum(100*clock)*rand(1));%¶¨ÒåËæ»úÖÖ×Óp_1=rand(1);%产生随机概率if p_1<=1 %如果随机概率小于0.9，则车辆离开路段，否则不离口matrix_cells(n)=0;v(n)=0;endend%入口边界，泊松分布到达，1s内平均到达车辆数为q，t为1sif b>vmaxt=1;q=0.25;x=1;p=(q*t)^x*exp(-q*t)/prod(x);%1s内有1辆车到达的概率rand('state',sum(100*clock)*rand(1));p_2=rand(1);if p_2<=pm=min(b-vmax,vmax);matrix_cells(m)=1;v(m)=m;endendnew_matrix_cells=matrix_cells;new_v=v;///////////////////////////////////////////////////////////////////////function [new_matrix_cells,new_v]=leadcarupdate(matrix_cells,v)%第一辆车更新规则n=length(matrix_cells);if v(n)~=0matrix_cells(n)=0;v(n)=0;endnew_matrix_cells=matrix_cells;new_v=v;///////////////////////////////////////////////////////////////////////function [new_v]=randslow(v)p=0.3;%慢化概率rand('state',sum(100*clock)*rand(1));%¶¨ÒåËæ»úÖÖ×Óp_rand=rand;%产生随机概率if p_rand<=pv=max(v-1,0);endnew_v=v;///////////////////////////////////////////////////////////////////////function [matrix_cells_start]=roadstart(matrix_cells,n)%道路上的车辆初始化状态，元胞矩阵随机为0或1，matrix_cells初始矩阵，n初始车辆数k=length(matrix_cells);z=round(k*rand(1,n));for i=1:nj=z(i);if j==0matrix_cells(j)=0;elsematrix_cells(j)=1;endendmatrix_cells_start=matrix_cells;///////////////////////////////////////////////////////////////////////function[location_frontcar]=searchfrontcar(current_location,matrix_cells)i=length(matrix_cells);if current_location==ilocation_frontcar=0;elsefor j=current_location+1:iif matrix_cells(j)~=0location_frontcar=j;break;elselocation_frontcar=0;endendend///////////////////////////////////////////////////////////////////////function [location_lastcar]=searchlastcar(matrix_cells)%搜索尾车位置for i=1:length(matrix_cells)if matrix_cells(i)~=0location_lastcar=i;break;else %如果路上无车，则空元胞数设定为道路长度location_lastcar=length(matrix_cells);endend///////////////////////////////////////////////////////////////////////function [location_leadcar]=searchleadcar(matrix_cells)i=length(matrix_cells);for j=1:iif matrix_cells(i-j+1)~=0location_leadcar=i-j+1;break;elselocation_leadcar=0;endend///////////////////////////////////////////////////////////////////////function [v_matixcells]=speedstart(matrix_cells,vmax)%道路初始状态车辆速度初始化v_matixcells=zeros(1,length(matrix_cells));for i=1:length(matrix_cells)if matrix_cells(i)~=0v_matixcells(i)=round(vmax*rand(1));endend。

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

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

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

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

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

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

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的报名参赛队号为（8位数字组成的编号）： A02004029 所属学校（请填写完整的全名）：天津科技大学参赛队员 (打印并签名) ：1. 裴彰明2. 刘伟3. 安兴娟指导教师或指导教师组负责人 (打印并签名)：刘素娟（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。

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

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

）日期： 2014 年 9 月 14 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：嫦娥三号软着陆轨道设计与控制策略摘要月球软着陆探测和巡视勘察是中国月球探测工程“三步走”战略的第二期任务目标。

2013 年12 月14 日，中国的嫦娥三号探测器成功着陆在月球虹湾地区，沉寂了近40 年的月球再次迎来了地球人的问候。

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

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

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

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

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

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

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

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

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

）日期：2014年9月15日赛区评阅编号（由赛区组委会评阅前进行编号）：2014高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：嫦娥三号软着陆轨道设计与控制策略优化摘 要 嫦娥三号是中国国家航天局嫦娥工程第二阶段的登月探测器，包括着陆器和玉兔号月球车。

嫦娥三号在高速飞行的情况下，要保证准确地在月球预定区域内实现软着陆，关键问题是着陆轨道与控制策略的设计。

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，可以满足调整速度的控制要求。

在四周安装有姿态调整发动机，在给定主减速发动机的推力方向后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。