全国大学生数学建模经验总结

第一篇：全国大学生数学建模经验总结

全国大学生数学建模经验总结

我坚信:只有想不到的，没有做不到的！

我是一名渝州学院大二学生，2010年加入学院建模队，2011年9月参加了全国大学生数学建模竞赛，获得江西省二等奖，虽然接触建模时间不是很长，但建模给我带来的却很多。

建模对很多人来说是很模糊的东西，但作为一名建模人就应该担任起对建模的责任。作为一名专科生，我知道自己相对于别人起跑线低，所以进入建模队后以严格的纪律来要求自己，别人懂的我得懂，别人不懂的我也要弄懂，不为别的，只因为要做就要做最好！

我虽在有些方面较强于别人，但我上课时还是认真听取老师的讲解，世上没有相同的人，每个人的想法思路也不可能完全相同，何不把别人的借鉴过来为己所用，多一个思路就证明多一条出路，多一条出路解题时分析题目的能力自然会比别人想的全面！所以在有些时候我还是把自己的思路讲给同学们参考，一起讨论解题的最好办法！

记得高三班主任在送我上大学的时候说过“不论什么时候多从别人的角度出发，凡事不要只为自己，你有足够的能力去做好任何事，多从事物的本身出发去考虑，要不做就不要做，但做就要做最好，以后的一切你自己把握！”担任过班长、学习委员、各科课代表的我曾对大学的录取不屑一顾，是人才总有发挥的地方，不论我的选择怎样家人、老师总会支持我。人生是用来闯的，不做错事是不可能的，但做错事后我会勇敢承认，人无圣贤，孰能无过！

加入建模队后我认识到，不是任何错误都可以犯的，也许某步棋没走对，那满盘就尽输了！所以不论每次训练我都警告自己失误不是每个人都犯得起的，建模建的是我们的思维，我们分析问题的能力，处理问题的方法，文字的表达能力。没有人会去当面问你解题的思路，我们唯一能做的就是把自己的思路想法用精炼的文字表达出来，这也锻炼了我对人处事的方式，以前总会想差不多就可以了，不去追求理解事物的本质，但建模不允许我们这么做，任何问题都要求我们刨根问底，对题目不理解何来的思路可循？

单纯从快递的运输方面举例来说，如果让我们设计一种最快、最低廉的运送方案，大多数人会考虑到运输设备的调用、运输路线的选择、接受地点的设置等方面，但对于从建模走过的人来说，我们会考虑的更深一层，例如该选用什么运输方式可以让运费最低，可以结合当地交通情况，利用单文件多种运输方式结合进行递送，通过对各种运输方式了解后，在保证运送时间相当的情况下，选用最低费用的运输方式，当然得结合实际情况。

总之，建模让我从另一个方面看世界，能让我更深入的分析理解问题，以致让我用最好的方法进行处理！

要做就做最好，不论什么时候只要确立了目标就要坚持到底！

第二篇：全国大学生数学建模大赛经验总结

比赛须知：

1、比赛期间参赛队员可以使用各种图书资料、计算机和软件，在国际互联网上浏览。

2、但是不得与队外任何人员（包括在网上）讨论，否则按舞弊处理。比赛期间不能大声喧哗，不能影响别的队。

6、数学建模大赛时间安排：

发题时间：2014年9月12日至9月15日上午8点

参加数学建模竞赛的十大秘诀诚信是最重要的数学建模竞赛是考查学生研究能力和实践能力的一场综合性比赛，有很多方面的知识和能力可以考查，但其中我觉得最重要的是诚信。我感到中国在这方面的教育还远远不够，我知道有很多同学写论文并不是实事求是地去做，而是编造数据、修改结论，明明自己没法编程实现却硬说自己做出来了，还编了一些数据。这些行为也许能够骗过评委，也许可以因“此”而获奖，但是这对他们将来是很不利的,希望能够引起足够的注意。团队合作是能否获奖的关键

在三天的比赛中，团队交流所占用的时间可能会超过一半。在一个小组中，出现意见不一是非常正常的，如果一个队意见完全一致，我想他们肯定不会拿奖。出现分歧的时候应当如何解决是很关键的，甚至直接决定你是否可以获奖，我的建议是“妥协”，这似乎是个贬义词，但我的意思是说不要总认为自己的观点是正确的，多听听别人的观

点，在两者之间谋求共同点。如果三个人都是自傲类型的人，也许每个人都非常强，但一旦合作，分歧就无法解决，做出来的就是一团糟，也就是说“三个诸葛亮顶不上一个臭皮匠”。我奉劝这样的话最好别组成一队了。合作在竞赛前就应当培养，比如一块儿做模拟题什么的，充分利用每个人的优点，也可以张三准备图论，李四准备最优化方法，然后几天后大家一块交流，这些都是可以磨合团队之间的关系的。通常在比赛时，三个人的分工是明确的，一个是领军人物，主要是构建整个问题的框架并提出有创意的idea，自然其他部分比如论文写比如程序设计比如计算他也能参加，应该算是一名全能型的人物;第二个是算手，顾名思义，主司计算方面的问题，比如编程计算一个微积分或者手工计算一条最优路径等。优秀的团队算手一般会精通(是精通不是入门)一个软件的应用，比如C比如MATLAB比如LINGO;最后一个是写手，主要工作在于论文的写作和润色上。好的论文要让人一眼就明了其中的意思，所以写手的工作还是需要一定的技巧的。当然，最重要的还是三个队员之间的讨论和交流，同心协力，在整个比赛过程中形成一种良好的交流氛围。时间和体力的问题

竞赛中时间分配也很重要，分配不好可能完不成论文，所以开始时要大致做一下安排。不必分的太细，比如第一天做第一小题，第二天做第二小题，这样反而会有压力，一切顺其自然。开始阶段不忙写作，可以将一些小组讨论的要点记录下来，不要太工整，随便写一下，到第三天再开始写论文也不迟的。也不要到第三天晚上才开始。另外要

说的就是体力要跟上，三天一般睡眠只有不到10 个小时，所以没有体力是不行的，建议是赛前熬夜编程几次，既训练了自己的建模能力，也达到了训练体力的目的，赛前锻炼身体我觉得没什么用处，多熬夜就行了，但比赛前一天可不许熬。重视摘要

摘要是论文的门面，摘要写的不好评委后面就不会去看了，自然只能给个成功参赛奖。摘要首先不要写废话，也不要照抄题目的一些话，直奔主题，要写明自己怎样分析问题，用什么方法解决问题，最重要的是结论是什么要说清楚，在中国的竞赛中结论如果正确一般得奖是必然的，如果不正确的话评委可能会继续往下看，也可能会扔在一边，但不写结论的话就一定不会得奖了，这一点不比美国竞赛，所以要认真写。摘要至少需要琢磨两个小时，不要轻视了它的重要性。很有必要多看看优秀论文的摘要是如何写的，并要作为赛前准备的内容之一。论文写作要正规

论文一定要大致按照摘要、问题重述、模型假设、符号说明、问题分析、(建立、分析、求解模型)、模型检验、参考文献、附录等等的方式来写。一篇论文结构上如果失败的话，比赛也一定不会成功，一般初评会先淘汰一些结构失败的文章，如果论文没有好的结构，内容再好也没有用。论文前面的结构一般都不会变，后面可以按照实际情况来安排，省略的部分可以有结果说明、灵敏度分析、其他模型、模型扩展、优缺点分析等等，多看些优秀论文就知道还有哪些形式了。附

录可以贴一些算法流程图或比较大的结果或图表等等。分析问题要认真

一般竞赛题目自己肯定没有见过，而且我发现近些年来的赛题都不是书上哪个模型可以直接套成功的，很多根本就没有固定的模型可以参考，所以分析问题不是一个去找书本的过程，依赖书本就意味着自己的思想被束缚起来。可以完全按照自己的分析去完成，平时练习的时候学习的是一种方法，通过以前学到的方法来解决，不是套用书本来解决，没有模型套怎么办，只有靠自己去实际分析。我估计在前面说的五点也许会有三分之一的队可以做到，而且可以做的很好，但是这一点上就需要真本事了，平时多努力，比赛发挥正常，这一点做好是没有问题的。编程求解是重要手段

美国竞赛时，美国学生中的论文很多是编程数据的说明，比如99 年

A 题行星撞地球那题，他们也能够模拟出撞击后果，这对我们来说简直是不可思议的。美国学生实践能力较强，而中国学生擅长理论分析，所以我把编程放在了分析的后面是有中国特色的。数学建模竞赛特别强调计算机编程解决实际问题的能力，最近几年尤其强调，编程方面的能力不是一朝一夕可以练成的，需要长期刻苦的训练，常用的工具有MATLAB、Mathematica、C/C++ 等等，一个人只需要会一门语言就行了，但需要精通它。比如要画柱状图该怎么做，要用Floyd 算法怎么办，赛前不准备是没有办法在比赛中很好运用的，因此每个常用的算法都自己去编程实现一下。模型的假设与模型的建立

评委看完摘要后紧接着就是看模型假设了，有一个万能的方法就是可以抄题目中可以作为假设的几句话，这样会给人留下好的印象，毕竟说明你审题了。但不能全抄，要加上自己的一些假设。一般假设用文字描述就行了，最好不要太具体了，一些重要参数不要被定死只能取某些值，否则会让人感觉论文的局限性较强。模型的建立是根据你对问题分析而来的，提出的数学符号和建立模型最好要比较接近，在同一页最好，以便评委可以对照符号来看，数学公式要严谨，推导要严密，这些都反映了参赛者的数学素质和能力，即使你推导不对，别人看到你的阵势也首先会误以为你是对的。那么多的试卷，评委不可能顺着你的公式一直推下去，但你要写得显得有数学修养才行。9 图文表并貌可以增色

我听说一个不确切的信息是评委老师喜欢用MATLAB 编程的论文，不知道有没有这回事，但这说明了老师需要看一个具有图或表在其中的论文，一篇如果像政治书那样写的论文估计没有人会对它感兴趣的，尤其是科技论文。MATLAB 编程之所以受到青睐是因为MATLAB 提供的图形处理能力很强大。图表的说明性特别强，如果结论有很多数据的话，最好做成图表的形式加以说明，会令你的论文更有说服力，也更容易受到评委的好评。其他

其他内容还是有很多的，说也说不完，挑几个重要的讲。比如不要上网讨论，网上的人水平参差不齐，你不知道谁是对的，而且很多人想

得奖，不会告诉你正确的，反而骗你说相反的，有时真理往往掌握在少数人手里。还有就是论文写作中灵敏度分析不要写太多，大致说明一下就可以了，不要喧宾夺主。最后想到的就是要使用数学公式编辑器来写论文，不要用什么上下标来表示，论文字体用小四，分标题用四号黑体等等。

第三篇：2014全国大学生数学建模竞赛

嫦娥三号软着陆轨道设计与控制策略

摘要

随着月球探测任务的发展,未来月球探测考察目标将主要是 复杂地形特性的高科学价值区域。为了能够安全地在这些遍布岩石、的区域内完成高精度软着陆，这就要求导航和控制系统具有较强的自主性和实时性。本文针对最终着陆段安全、精确的需求，对月球软着陆导航与控制方法进行较深入研究，主要内容包括：

首先，提出一种基于单帧图像信息的障碍检测方法。该方法根据着陆区内障碍成像的特点，通过匹配相应的阴影区与光照区完成对岩石、弹坑的检测，利用图像灰度方差对粗糙区域进行提取：在检测出故障信息的基础上，选取安全着陆点以保证软着陆任务的成功。

其次，给出一种基于矢量观测信息的自主光学导航方法。该方法利用光学相机和激光测距仪测量值构建着陆点相对着陆器的矢量信息，结合着陆器的姿态信息确定着陆器的位置。为了消除测量噪声带来的干扰，利用扩展Kalman滤波理论设计了导航滤波器。

再次，提出一种李雅普诺夫函数障碍规避制导方法。该方法通过对状态函数、危险地形势函数的设计，以满足平移过程中减低障碍威胁与精确定点着陆器，设计PWPF（调频调宽）调节器实现定推理等效变推力控制效果。

最后，针对采用变推力主发动机的月球着陆器，提出一种垂直软着陆控制方法。该方法采用标称控制与闭环控制相结合的方式，规划标称轨迹以保证着陆器到达着陆点时其下降速度、加速度亦为零，设计闭环控制器产生附加控制量消除初始偏差、着陆器质量变化的干扰，以保证着陆器沿标称轨迹到达着陆点。

本文分别对所提出的最终着陆段导航与控制方法进行数学仿真以验证个方法的可行性。仿真结果表明，本文多给出导航方法能够达到较高的性能指标，满足在危险区域实现高精度软着陆的需要。

关键词： 月球软着陆；自主导航与控制；障碍检测；规避制导；适量测量

一、问题重述

嫦娥三号于2013年12月2日1时30分成功发射，12月6日抵达月球轨道。根据计划，嫦娥三号将在北京时间12月14号在月球表面实施软着陆。嫦娥三号如何实现软着陆以及能否成功成为外界关注焦点。嫦娥三号在着陆准备轨道上的运行质量为2.4t，其安装在下部的主减速发动机是目前中国航天器上最大推力的发动机，能够产生1500N到7500N的可调节推力，进而对嫦娥三号实现精准控制。其比冲（即单位质量的推进剂产生的推力）为2940m/s，可以满足调整速度的控制要求。在四周安装有姿态调整发动机，在给定主减速发动机的推力方向后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。嫦娥三号的预定着陆点为19.51W，44.12N，海拔为-2641m。嫦娥三号将在近月点15公里处以抛物线下降，相对速度从每秒1.7公里逐渐降为零。整个过程大概需要十几分钟的时间。在距月面100米处时，嫦娥三号要进行短暂的悬停，扫描月面地形，避开障碍物，寻找着陆点。之后，嫦娥三号在反推火箭的作用下继续慢慢下降，直到离月面4米高时再度悬停。此时，关掉反冲发动机，探测器自由下落。

嫦娥三号在高速飞行的情况下，要保证准确地在月球预定区域内实现软着陆，关键问题是着陆轨道与控制策略的设计。其着陆轨道设计的基本要求：着陆准备轨道为近月点15km，远月点100km的椭圆形轨道；着陆轨道为从近月点至着陆点，其软着陆过程共分为6个阶段，分别为着陆准备轨道、主减速段、快速调整段、粗避障段、精避障段、缓速下降阶段，要求满足每个阶段在关键点所处的状态；尽量减少软着陆过程的燃料消耗。

根据上述的基本要求，请你们建立数学模型解决下面的问题：

（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。

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

二、问题分析

对于问题一：

嫦娥三号从15公里左右的高度下降到月球表面，在这一过程中不考虑月球表面太阳风的影响，忽略月球的自转速度引起的科氏力的影响，由于下降时间比较短也不考虑太阳、地球对嫦娥三号的摄动影响，嫦娥三号水平速度要从1.692km/s降为0m/s由于3000m处时嫦娥三号已经基本位于着陆点上方，所以此时假设在3000m处的速度只存在竖直向下的速度而不存在水平分速度，因为降落减速时间比较短只有垂直于月面的方向运动才能实现，所以在确定着陆点位置和着陆轨迹时应当考虑燃料最优情况下推力最大，方向自由的方法即取F7500N建立主减速段动力学模型。

三、符号说明

四、模型假设

对于问题一：

忽略月球的自传和太阳、地球对嫦娥三号卫星的引力摄动 月球近似为一个质量均匀的标准球体 将嫦娥三号是为一个质点

主减速忽略动作调整所产生的燃料消耗段不考虑太阳风的影响

五、模型建立与求解

5.1问题一的建模与求解 解法一： 假设嫦娥三号在t时刻在远月点开始缓慢下降，在n时刻到达近月点，整个过程遵循开普勒第三定律，即

v00

在t时刻有：v12R1 R0R0R1r0 R0r1r2 其中v1：远月点速度

v2：近月点速度

R0：远月点月心距

R1：近月点月心距（已知月球的半径为1738千米）

R017381001838km

R11738151753km 在t1时刻处v2 k2R1 R0R0R1R00.512k0.488 R0R1利用能量平衡式求得近地点速度为

20.51249012（）1.692km/s（沿切线方向）v2，比当地的环境速度17531.672km/s大vk0.0196km/s，径向速度vk0。

1同理解得v11.6139km/s（沿切线方向）

vri0

解得主减速段动力学模型的建立：

根据题意，在横向飞行的水平距离远远小于月球半径的平均值，所以可以将整个减速段过程简化为水平和竖直方向运动方程，根据牛顿第二定律、速度计算公式有：

axTx maytTymTxta

1.692km/s m0Qdt0Tyadt57m/s t0mQdt0tT22xTy27500N

v22atS

运用matlab编程解得S451810.4m； 其中 ax：水平方向加速度

ay：竖直方面加速度

a：月球表面重力加速度a Tx：推力的水平方向分力

Ty：推力的竖直方向分力

t:主减速段时间

S:嫦娥三号主减速段水平位移

Q:嫦娥三号发动机燃料秒消耗率

根据已知资料得到嫦娥三号着陆过程中纬度改变，经度基本不变，月球赤纬和地球纬度一样也分为南北各90个分度,又因为月球极区半径为1735.843km，所以每一个纬度的竖直高度差为19.2871

4g 6千米。即近月点位置坐标为19.0464W,28.9989N海拔15km,远月点位置坐标为160.9536E,28.9989S海拔100km。

解法2:轨迹方程法。

众所周知,太阳系中的八大行星都在按照各自的椭圆轨道绕太阳进行公转,太阳位于椭圆的一个焦点上,行星的运动遵循开普勒三定律,笔者发现,在各类物理竞赛中,常会涉及到天体运动速度的计算,本文拟从能量和行星运动的轨迹方程两个不同的角度来探索行星在近日点和远日点的速度。

该解法的指导思想是对椭圆的轨迹方程求导,并结合一般曲线的曲率半径通式求出近日点和远日点的曲率半径表达式,然后利用万有引力提供向心力列方程求解。如图1所示,椭圆的轨迹方程为

x2y221 5 2ba将5式变形为

a2x2b2y2a2b2 6

根据隐函数的求导法则将6式对x求导有

2a2x2b2yy0 7 即

a2xy2 8

by将7式再次对x求导得

2a22b2(yyyy)0 9 将8、9两式联立得

a2b2y2a4x2 10 y-43by根据曲率半径公式有 r(1y)11 y122 将8、10、11式联立并将A点坐标A(0，a)代入可得A点的曲率半径为

b2RA 12

a根据椭圆的对称性,远日点B的曲率半径为

b2RBRA 13

a 由于在A、B两点行星运行速度方向与万有引力方向垂直,万有引力只改变速度方向,并不改变速度大小,故分别根据万有引力提供向心力得

GMmmvA 14 (ac)2RAGMmmvB 15 2(ac)RB将13至15式联立可得 22vAbGMbGM，vB acaaca

5.2问题二的建模与求解 模型一：动力学模型

典型的月球软着陆任务中,探测器一般首先发射到100km的环月停泊轨道,然后根据所选定的着陆位置,在合适的时间给着陆器一个有限脉冲,使得着陆器转入近月点(在着落位置附近)为15km,远月点为100km的月球椭圆轨道,这一阶段称为霍曼转移段。当着陆器运行到近月点时,制动发动机开始工作,其主要任务是抵消着陆器的初始动能和势能,使着陆器接触地面时,相对月面速度为零,即实现所谓的软着陆,这一阶段称为动力下降段。着陆器的大部分燃料都是消耗在此阶段,所以月球软着陆轨迹优化主要是针对动力下降段这一阶段。由于月球表面附近没有大气,所以在飞行器的动力学模型中没有大气阻力项。而且从15km左右的轨道高度软着陆到月球表面的时间比较短,一般在几百秒的范围内,所以诸如月球引力非球项、日月引力摄动等影响因素均可忽略不计,所以这一过程可以在二体模型下描述。其示意图如图1所示,其中o为月球质心,x轴方向为由月心指向着陆器的初始位置,y轴方向为初始位置着陆器速度方向。

图 1 月球软着陆极坐标系

其动力学方程如下: rv 

v(F/m)sin/rr

22 ((F/m)cos2v)/r

mF/ISP

在上式中r为着陆器与月心距离,v为着陆器径向速度,为着陆器极角,为着陆器极角角速度,为月球引力常数,F着陆器制动发动机推力,m为着陆器质量,为制动发动机推力方向角,其定义为F与当地水平方向夹角,ISP为制动发动机比冲。根据动力下降段的起点位置可以确定动力学方程初始条件,由于起点处于霍曼转移轨道的近地点,故其初始条件为: r0rp

00

v00 01rprp(2ra)rarp其中rp和ra分别为霍曼转移段的近地点半径和远地点半径。

终端条件为实现软着陆, 即

rfR

vf0

f0

其中R为月球半径,终端条件中对终端极角f及终端时间tf无约束。

优化变量为制动发动机推力方向角(t)。

优化的性能指标为在满足上述初始条件和终端条件的前提下, 使着陆过程中燃料消耗最少,即

Jm(t)dt

t0f设计主减速段制导控制律 2动力下降段燃料最优精确着陆问题描述 2.1 燃料最优精确着陆问题

着陆器运动方程：考虑采用变推力发动机情况，有

rv

.vga

(1)

aTmmaT..其中r[rhrxry]T，v[vhvxvy]T分别表示着陆器相对期望着陆点的位置和速度矢量；T为推力器提供的推力矢量，幅值为 T，对应控制加速度矢量 a；g为火星的重力加速度矢量，此处认为是常值；m为着陆器质量，对应推力器质量排除系数。指标函数：考虑燃料消耗

min(m0mf)min0fTdt

(2)边界条件：即初始条件和终端条件

r(0)r0,v(0)v0,m(0)m0,r(tf)v(tf)[000]

(3)控制约束：考虑发动机一旦启动不能关闭，存在最大和最小推力约束

0T1TT

2(4)状态约束：为避免在着陆前撞击到火星地表，需确保整个下降段位于火星地平面以上，即

rh0

（5）进一步地，若着陆区域附近表面崎岖不平，仅仅确保地表约束不能满足需求时，可以考虑下降倾角约束，即将着陆器下降轨线约束到以着陆点为顶点的圆锥体内

2.2 等效后燃料最优精确着陆问题 定义等效变换变量

Ttrx2ry2rhtanalt

（6）

uaT

m

（7）

Tmzlnm等效着陆器运动方程： .r0I3..

yv00.00z其中p[uT0r0vI030z07*70ug0AcyBc(pg4)

（8）],g4[gTT0]T

t指标函数：

min0f(t)dt

（9）

边界条件：同式(3)。

控制约束：由文献[10]可知，控制约束(4)可等效表示为

u1T1ez0[1(zz0)(zz0)2]T2ez0[1(zz0)]

（10）（11）

2状态约束：地表约束同式(5)，倾角约束(6)可等效表示为

T

Sycy0

（12）

其中

0100000S

0010000ctanalt

T000000

3.燃料最优精确着陆问题的离散化及变换 3.1 等效燃料最优精确着陆问题的离散化

首先将整个飞行时间均分成 n 段（对应 n +1 个点），每段步长为t，离散化后的着陆器运动方程为：yk1AykB(pkg4)

其中AR77,BR74分别为离散系统的系统矩阵和输入矩阵

12AetAcI3tActAc

2tt112BetsAcBcdsesAcdsBctBctBct2Bc

0026其中I3为三阶单位阵。

有系统性质可知，整个控制时域内系统状态满足 y3Ay2Bp2g4A3y0A2Bp0g4ABp1g4Bp2g4ynAyn1Bpn1g4Any0An1Bp0g4ABpn2g4Bpn1g4y1Ay0Bp0g4y2Ay1Bp1g4A2y0ABp0g4Bpn2g4Bp1g4

为表达方便，令

y0p00A0yp1111A ,pp2，2A2 Yy2nyn7n11pn4n11nA7n1700B1AB223ABn1An则(15)可等价于

000B012ABBB0002 ABB003AABB0n1AABBA2BABBn7n14n1000000Yy0pg4

分别定义如下常值矩阵：

最终可得离散化后的燃料最优化问题如下： 指标函数：式(9)可表示为

边界条件：式(3)可表示为

控制约束：式(10)和式(11)分别可表示为

状态约束：式(5)和式(12)分别可表示为

含有 p个线性约束和 q个二阶锥约束的最优化问题的标准形式为 指标函数

min(Tx)满足约束

DTxf0AxcibdinTiTi

（k=1,，n）

n*pp其中xR为待优化向量，R，线性约束参数DR,fR，二阶锥约束参数维数n(Ai,bi,ci,di)由相应约束确定

则式(17)～式(23)可最终转换为如下最优化问题： 指标函数：min(vpp)满足：

初值约束:MxΨ0pMx(Ψ0y0)A0g4r0末值约束:MxΨ0pMx(Ψ0y0)A0g4控制约束:Murkpvrkp 控制上限:(vzΨkTTTTv0T0

0

T1vr)p1vTz(Φky0Akg4)z0,z0 z0kT2e 控制下限：

4数值仿真结果与分析本节以某火星着陆器为例，计算了典型初始条件下满足各种约束的燃料最优精确着陆轨迹。其中探测器各参数分别取为：m02000kg,g[3.711400]ms2,c2kms,T11.3kN,T213kN.。着陆器初始位置矢量r0= [1500,-600, 800] m，初始速度矢量v0= [-30, 10, 40]m/s，倾角alt=86°。二阶锥优化问题可以通过大量免费的优化工具求解，如 CSDP、DSDP、OpenOpt、SeDuMi、SDPA、SDPLR等。本文选用 SDPT3 进行计算，通过执行线性搜索确定燃料最优下降时间tf为 43s，图 1 给出了相应的最优着陆轨迹、下降速度、加速度、控制推力、推力幅值以及探测器质量变化曲线。

由优化结果可以看出，探测器在给定时间飞行并软着陆到指定位置，且在整个下降过程始终与火星地表保持一定的安全距离，验证了下降倾角约束的有效性。其推力幅值曲线呈现“最大-最小-最大”的最优控制形式，不过为了保持发动机始终处于点火状态，在中间段对应最小推力约束，这与文献中的分析结论一致。此外，通过利用如 TOMLAB 等商业最优控制软件进行复核计算，也验证了此计算结果的燃料最优性能。

*

图 1 给定初始条件下火星着陆器动力下降段燃料最优计算结果

需要注意到，此燃料最优轨迹的获取对着陆器的实时在线计算性能提出了较高的要求，经测试，无论使用何种优化工具，计算给定飞行任务时间的最优轨迹均需数秒，而全局最优则需要数十秒甚至更长，这在实际任务中是不允许的。因此，可行的方案是通过在地面计算大量的燃料最优轨迹，并寻找规律，选取关键路径点状态存储到着陆器计算机中，通过在线查表或者在利用对计算量要求较小的反馈制导律完成安全着陆任务。

因此，为了研究探测器燃料最优轨迹特性，选取相同的探测器参数，暂不考虑推力器最小幅值约束和倾斜角约束（但考虑地表约束），固定初始高度为 1500m，初始位置水平方向从-8000m 到 8000m 内取值，分别选取各种不同的初始速度，可得燃料最优精确着陆轨迹簇如图 2 所示。

图 2 各种不同初始速度对应的火星着陆器动力下降段燃料最优轨迹簇

1)对任意探测器初始位置，特定初始速度对应的燃料最优着陆轨迹在末端必然收敛到一个固定的近似圆锥体内。

2)取决于探测器初始位置和速度的关系，燃料最优轨迹有两种形式：S 型和 C 型，其中 S 型主要对应于期望着陆点位置水平距离较大情况。3)当探测器初始水平速度为零时，圆锥体轴线垂直于火星地表，所有最优轨线关于该轴线中心对称。4)初始速度的大小也直接影响到任务的可靠性，因此需要在超声速进入段和降落伞减速段将着陆器速度下降到合理范围内。

上述结论对上注探测器关键点的选取有着较强的指导意义，比如基于最优轨线的斜率对路径点合并、基于最优轨线簇的对称性对上注轨线进行等效延伸、或者尝试仅将 S 型和 C 型的转折点作为路径点等，这样可以大大降低探测器自主存储与计算需求，进而有效提升任务的可靠性。重力转弯软着陆过程

对于最终着陆点，假设探测器的下降轨迹在一平面内，且月球引力场为垂直于月面XY的均匀引力场，引力加速度g沿-Z，如图1所示，制动推力方向沿探测器的本体轴z。重力转弯软着陆过程中探测器质心动力学方程可表示为

上式中各变量的物理意义如图1中所示，其中m>0为探测器质量；k>0为制动发动机比冲；u表示制动发动机的秒耗量

可通过一定的机构加以调节，故作为软着陆问题的控制变量。假定制动发动机的最大推力与初始质量比大于月面引力加速度，并且制动推进系统能够在一定的初始条件下将探测器停止月面上。

重力转弯过程中，探测器的高度、速度和姿态角度可由雷达高度表、多普勒雷达及惯性仪表测得。令软着陆初始条件探测器到达月面时速度减小到给定的值，故终端条件自由。软着陆燃耗最优问题的描述 对于最终着陆段，可假设

为一小角度。由此可将系统方程（1）化简为

要设计制导律实现软着陆，就是使

着陆时间

对于月球软着陆的燃耗最优控制问题，其性能指标可表示为

对于系统（2）的软着陆过程，燃耗最优问题等价于着陆时间最优问题，性能指标为

在月球重力转弯软着陆过程中，如果存在一个推力控制程序将探测器从初始条件转移到终端条件，并使性能指标（3）或（4）式最大，则称这个推力程序为软着陆燃耗最优或时间最优制导律。根据pontryagin极大值原理，系统的哈密顿函数及其对u的偏导数为

使哈密顿函数（5）式达到极大地控制输入u就是最优控制，科表示为。

如果存在一个有限区间

则最优控制u（t）取值不能由哈密顿函数确定。此时如果最优解存在，则称为奇异解，（8）式称为奇异条件。

最优制导问题的性质：1）对于自治系统（2）的时间最优控制问题，沿最优轨迹其哈密顿函数满足

将其对时间求导并将（2c）和（6c）式代入，得

另外，由于自由，根据横截条件有3）根据（6a）式。又由（9）式可得T（t）=0，4）根据极大值原理，系统的状态变量和共轭变量都是时间的连续可微函数，将切换函数对时间求导，利用（2），（6）式和性质2）得 软着陆最优控制中奇异条件的分析

对于月球重力转弯软着陆问题，最优制导律具有两个很好的性质。

定理一。月球重力转弯软着陆系统（2）的燃耗最优制导或时间最优制导问题不存在奇异条件。证明。用反证法，假设存在奇异条件，则在某个闭区间设，并由（5）式得

。根据反正假将（10）式两边对时间求导，并将（2）和（6）式代入化简得性质2），并考虑到或者情形1.得

下面证明这两种情形均与反证假设矛盾。根据式

及性质2）可知，由性质3）必有

根据

是时间t的斜率非零的线性函数，m和情形2.1）若定，根据横截条件有在区间内为常数。这与反证假设矛盾。

。下面再分三种情况进行分析。

又因为

不与此时由（6b）式有反证假设矛盾。2）若盾。3），与反证假设矛又

因

为

因此有成立，这与

此时（10）式在上根据定理一，重力转弯软着陆的最优制导律是一种开关（Bang-Bang）控制，只须控制发动机开关，不需要调节推力的大小。

定理2.对于月球重力转弯软着陆过程，其开关控制器的最优推力程序（7）最多进行一次切换。

证明。只要证明最多只在一个时间点成立即可。软着陆系统（2）在最优推力控制程序（7）的作用下，按最后轨迹降落。由性质3）知，为常数。根据性质4），若严格单调，因而在上至多有一个零点，即至多进行一次切换；若，则上为常数。由定理1，5 软着陆最优开关制导律

不可能在任何区间上成立，故必有既没有切换点。

对于最优推力控制程序（7），其切换函数中含有共轭变量，它是一个关于状态变量的稳式表达式。为实现实时制导，需求出关于状态变量的切换函数来。

根据定理一和定理二，重力转弯软着陆最优控制程序没有奇异值状态，并且在着陆过程中最多切换一次，其工作方式有4种：1）全开；2）全关；3）先开有关；4）先关后开。对于方式1）软着陆起始点即是开机点；方式2），3）不能实现软着陆；最后一种是通常情况下的最优着陆方式，即探测器先做无制动下降，然后打开发动机软着陆到月面。设开机时刻为到发动机工作时间为

式，在区间

内积分，并考虑

将（11）式中的对数按泰勒展开，忽略

并令

消掉T得到切换函数为

由切换函数（12）式可以看出，速度、位置的误差和制动发动机推动的将直接影响着陆的效果。一种方法是将终端高度从到达月面时实现软着陆设置为离月面还有几米时实现软着陆。另一种方法是考虑制动过程由一个主发动机和一组小推力发动机共同完成，通过调整开启的小发动机的数量，来实现变推力降落。具体地，令切换函数为

式中各符号的含义如图2所示

关机点可取为2m,可取为20m，可取为1m/s。为实现着陆的最优性，减速度

取为

其中T如（12）式中所示，m0为探测器的初始质量。

图三为最优着陆过程与其改进方法按图2降落的次优着陆过程的对比图。由此图中可看出，改进方法提高了着陆的安全性，当探测器的初始质量mo=350kg，发动机着陆过程多消耗燃料2.2kg。

时，改进方法比最优

（a）

（b）

问题三 协方差分析方法的基本原理 对于如下非线性函数关系

yfx1,x2xn（1）

可以使用一阶泰勒级数展开对其进行线性化，有

yyfffx1xnx1xn（2）x1xn其中，x1xn为x1xn的高阶项。从而得到线性化方程

yfxi（3）i1xin或表示为

YPX(4)

这里 P 是偏导数矩阵: Pif（5）xi若自变量x1xn是随机变量，则线性化方程的函数y的协方差矩阵为：

EYYTEPXXTPTPEXXTPT(6)即 CyPCXPT(7)式中Cx是自变量的协方差矩阵；Cy是函数Y的协方差矩阵。

协方差矩阵中对角线元素是方差，非对角线元素为协方差。显然，只要求出传递矩阵 P ,便可确定源误差与欲求量误差之间的关系。若给定各种源误差，如发动机安装误差、敏感器测量误差或发动机推力和点火时间等误差时，便可以分析其对目标轨道误差的影响以及对控制系统精度的影响，进一步对各系统及元部件提出适当的精度要求。计算向月飞行轨道误差的协方差迭代方程

考虑到轨道参数的误差之相对于轨道参数的标称值是小量，因此可以将轨道运动方程进行线性化，从而得到能够反映轨道参数偏差量的传播关系的误差方程。在应用双二体模型且在地球影响球范围内时，对轨道运动产生摄动影响的各项，如月球引力摄动、太阳引力摄动、大气阻力摄动和太阳光压摄动等对误差方程的影响很小，因此在误差方程中将它们忽略掉。反映轨道位置和速度误差的线性化方程如下：

vrg（8）vrrTur，其中u为地球引力常数。式中 gr3rrrx2ry2rz2（9）

写成状态方程形式：

0Irr(10)vG0vg式中 GT

r0Ir令FG0,Xv（11）

则式（9）变为

FX(12)X下面推导矩阵 F 的表达式：

guGTT3rrrruurrT33Trrrruuuur3333I3rrrrryzrxr（13）

式中 r x，r y 和 r z 是探测器在地心惯性坐标系里的轨道位置坐标。则Gu3T(Irr)(14)332rrrx2rxryrxT2rrryrxryrzryrxryrrzrxrzryzrxrzryrz（15）2rz

将式（15）、（14）代入（10），得： 0002-urx(132)Fr3r3urxryr5v3urxrzr5

积分式（11），得到： 0003urxryr520003urxrzr53urzryr5210000ry-u(13)32rr3urzryr5-urz(13)0r3r200100100（16）

0000

XteFtX0

（17）式中

(Ft)2(Ft)3(Ft)4(Ft)neIFt2!3!4!n!

（18）iNtFi.()i!i0Ft取前 6 阶截断，即：

eFttiFi!

（19）i06i

得到计算误差方程的迭代方程：

XtiteFtXti

（20）

eFt相当于式（4）中的 P 阵，由于误差方程是时变方程，因此每一步迭代都需要重新计算 P 阵，计算 P 阵需要利用标称轨道参数数据。

进一步根据式（7），得到协方差矩阵的迭代方程：

T

Ci1PCPiii

（21）向月飞行轨道误差的协方差分析

引起轨道误差的误差源主要是导航误差，包括位 置 误 差 和 速 度 误 差。其 中 ： 位 置 误 差 ：rrx,ry,rz,rx,ry,rz分别为在地心惯性坐标系中 X 轴、Y 轴、Z 轴的分量。速度误差：vvx,vy,vz,vx,vy,vz分别是在地心惯性坐标系 X 轴、Y 轴、Z 轴的分量。向月飞行轨道的初始轨道位置和速度误差由运载火箭的发射入轨精度决定，若探测器在飞行途中进行轨道修正，则经过轨道修正以后的轨道位置误差将由导航误差决定，速度误差将由姿态误差和制导误差决定。

上述误差决定了轨道误差协方差分析的计算初始条件，表 1 给出了在不进行中途轨道修正情况下，在地心惯性坐标系里，初始轨道位置误差和初始速度误差对轨道终点的位置和速度误差的影响。图 1 和图 2 给出了在算例三中探测器从近地轨道入轨点开始至进入月球轨道为止轨道位置的相应的轨道位置和速度总误差（3σ）的时间历程。

表 1 初始轨道位置和速度误差

对轨道终点误差的影响

图 1 轨道位置总误差时间历程（3σ）

图 2 速度总误差时间历程（3σ）基于敏感系数矩阵的制导误差分析

在月球软着陆主制动段，影响制导精度的误差源主要有偏离标准飞行轨迹的初始条件误差和导航与控制传感器误差。初始条件误差由主制动段以前的任务决定，传感器误差则由导航系统和传感器本身决定。此外，影响制导精度的因素还包括月球自转、月球不规则摄动等误差，对它们的研究可单独进行，这里暂不做介绍。2.1 误差模型建立

2.1.1 初始状态误差模型

记着陆器的实际初始状态为Xi，标准初始状态为Xn,则定义初始状态偏差xi为

xiXiXn

(7)对于主制动段这一特定的飞行过程，这些偏差都是确定的；而针对整个月球探测任务，这些偏差就变得具有随机性。在本文中，假定xi 的所有元素均服从零均值高斯分布，相互不独立，其相关性取决于前一阶段任务的特性。2.1.2 传感器误差模型

由于只研究误差对制导律的影响，所以这里假设需要测量的量均可由导航系统直接测得，误差大小

均考虑为典型误差值。由上一目设计的制导律可以看出，需要由导航与控制传感器测量的量主要为着陆器相对于着陆场坐标系的位置、速度和加速度。定义待测量量Q为

QX其估计值记为Q，则传感器误差定义为 YZUVWA

T

qQQ

(8)那么，单个测量量的估计误差模型可用误差向量 q的第j(j =1，2„7)个元素qj 来表示。由参考文献[5]可知，第 j个观测量的总估计误差qj 由以下四部分组成

~~-~qjbsqjnstqtqQtqtQjt

(9)jjbcjnc

j100100~~~~~针对主制动这一特定操作阶段，上述四部分误差具有如下特性：

qjbc—第 j 个观测量的测量误差，恒为常值，其分布服从零均值高斯分布； qjbs—第 j 个观测量的刻度因素误差系数，恒为常值，其分布服从零均值高斯分布； qjnc—第 j 个观测量的随机误差，其为一高斯白噪声；

qjns

—第 j 个观测量的刻度因素随机误差系数，其为一高斯白噪声。

2.2 制导误差分析

由于采用闭环制导，制导控制系统对随机误差具有一定鲁棒性，所以本文将着重对初始偏差和类似于qjbc和qjbs这样的传感器常值误差进行仿真研究，分析它们对制导精度的影响。2.2.1 误差分析系统建立

误差分析系统框图如图 1 所示，下面将对其结构进行分析。~~~~~~

图 1 误差分析系统结构图

图中所示初始状态偏差实际上是加在相应积分器中。

由前面的分析可知，观测量的实际输出值受到初始状态偏差、传感器测量误差以及传感器刻度因素误差的影响，故误差分析系统模拟程序的实际输入应包含以下几部分(以 X通道为例)：

XXxixbc~xbsX

(10)100~~

其中，X为观测量的实际输出值，X 为标准值，xi 为初始状态偏差(只在初始时刻存在)，xbc 为传感器测量偏差，xbs为传感器刻度因素误差系数。由图 1 可以看出，为了更准确地表示传感器误差模型，这里考虑了传感器的动态性能，其传递函数设为一阶惯性环节11Ts，其中，T 为传感器时间常数，因传感器的不同而取不同值。

由误差分析系统结构框图可以看出，其输入量主要包括：标准初始状态向量、初始状态偏差、传感器测量误差、传感器刻度因素误差系数、传感器时间常数、期望终端状态；输出量为加入误差前后的仿真终端状态向量。2.2.2 误差敏感系数矩阵求取

在有形如(7)式误差输入的情况下，首先根据图 1 生成一个模拟整个闭环制导控制系统的数字仿真程序，然后运行该程序，对比程序输出即可得到误差敏感系数矩阵。具体运行过程如下：

第一步：将传感器误差设置为零，初始状态设置为标准值，运行模拟程序。这一步称为标准运行。第二步： 将其中一个传感器误差设置为非零输入或者设置一个非标准初始状态，然后进行一系列运行。

第三步： 将第二步运行的系统输出和标准运行的系统输出进行比较即可确定各误差源的影响。如X 通道标准初始偏差为xi，输入该误差前后，X 通道终端状态分别为X0 和X1，则 X 通道对标准初始偏差xi的敏感性可用(X1X0)/xi来反映。

通过这种方法，可得到一组反映月球软着陆主制动段终端总误差向量pf和两个传感器误差向量~~~qbc、qbs以及初始状态偏差向量pi之间关系的误差敏感系数矩阵。由参考文献[6]可知，其相互关系可表示为

~~pfS1piS2qbcS3qbs（11）

其中，S1、S2和S3分别表示相对于pi、qbc和qbs的误差敏感系数矩阵。

终端误差向量能用这种形式表示的假设条件是动力学的线性化必须在标准轨迹区域内。验证该假设条件的方法有两种： 扩大输入误差仿真法和复合仿真法，这里略去其验证过程。2.2.3 误差分析

假设导航系统采用常规惯性测量单元，表 1 列出了其典型误差值，其中，位置误差能保持在10数量级，速度在10数量级，加速度为 10g 数量级。1-52~~

运用上述方法得到的敏感系数矩阵给出如下：

5.50210-3-4-3.850101.69210-3S1-38.36210-5.86010-4-3-2.57510-2.08010-4-1.05010-31.41810-11.40110-57.30110-5-1.00110-26.41110-53.24010-4-4.40710-2-2.57010-4-1.86210-3-5.58010-11.41010-57.90210-51.31210-55.71010-4-1.15710-38.10010-53.93610-21.73210-2-2.7431017.74610-1-4.02410-2-8.93910-23.21010-34.03010-31.23910-21.83310-2-2-18.742101.41410-1.19610-2-9.90110-3-2-2-2.69010-4.57710-6.81210-1-8.69510-2-5.2031002.11010-14.23510-16.17010-3-3.2811008.20210-2-5.76010-35.63310-1-3.4891022.4431014.401102-9.8331026.86410123.02010-9.85910-1-1.15410-3-40-3.13010-1.00010-1.37910-33.56010-4S2-2-3-5.402101.540101.04510-31.86410-3-34.77010-44.598109.99910-13.408100-7.21010-43.5041005.00010-55.64310-3-1.52710-19.36810-1-6.72110-1-1.30610-1-5.6314100-28.479103.73010-1S30-8.924104.61910-102.03310-5.49410-1-3.53310-1-2.8101001.60010-31.69210-16.75510-18.99610-1-209510-12.47310-21.66410-1-1.0271007.16510-23.344100-1.1121008.61310-17.8521003.246100-1.6181003.54010-14.98210-17.67010-1-1.122100-2.397100-2.38010-1-3.650100-2.5631002.55610-1-4.29110-23.401100-1.88810-1-5.103100-3.23010-13.56610-12.25610-10-1-7.005109.93010A1、A3:12.7592,30.1297j2.1329 A2:11.5522,30.6761j1.8978

由于数值仿真的起始点选为(1，0，-1)，靠近平衡点(1.5，0，-1.05)，仿真实验中混沌系统的基频w0=2.1329，基周期为为T0202.9443S。由前面的数值仿真实验知要使 Chua’s混沌系统保持其类随机性，仿真步长选在(0.0001，0.7)较为合适，用基周期来表达即为129940T015T0 ,15T0内，综观三个连续混沌系统仿真步长的理论计算，我们可以统一选取15000T0这样即可以提高仿真运算速度，又可以使混沌吸引子的形状和类随机性不发生变化，这个选择范围也与通常连续混沌系统数值仿真步长的经验取值相吻合六、模型结果及分析

七、结果分析

八、模型评价与改进方向

九、参考文献

第四篇：全国大学生数学建模感想

数学建模有感

这一暑假，参加学校组织的学习了 “数学建模”。这门课程，数学建模这门课程与其它的数学课程有一个很大的区别，数学建模是把现实世界中的实际问题加以提炼，抽象为数学模型，求出模型的解，验证模型的合理性．并用该数学模型所提供的解答来解释现实问题，它是一门将数学综合应用到实际中解决实际问题的学科。建立模型就是对于现实中的原型．为了某个特定目的，作出一些必要的简化和假设。运后适当的数学工具得到一个数学结构、它或者能解释特定现象的现实状态，或者能预测到对象的未来状况，或者能提供处理对象的最优决策或控制，而现实生活中很多的问题都可以通过建立数学模型来解决，所以学习数学建模，彻底让我体会到了数学在生活中的重要之处。

一次建模，终生受益，这话一点也不假。在没有接触过数学建模这门学科之前，我一直认为数学是一门纯理论的学科，但是数学建模却能把它应用到实际中去，并用它去解决很多来自日常生活及经济、工程、理、化、生、医等学科中的问题。

我觉得数学建模能为学生提供自主学习、自主探索、自主提出问题、自主解决问题的机会，培养学生的数学观念、科学态度和合作精神，激发学生的学习兴趣，培养学生认真求实、崇尚真理、追求完美、讲究效益、联系实际的学习态度和学习习惯。它能提高学生应用所学的数学知识解决实际问题的能力，从过去强调数学知识的“有用、可用”，到使学生所学知识的“想用、能用和会用”，让学生更多自主的实践，把学习知识、应用知识、探索发现、使用计算机、培养良好的科学态度与思维品质更好地结合起来，使学生在问题解决的过程中得到学数学、用数学的实际体验，加深对数学的理解。综合上述可见，开展数学建模活动是非常有必要的。应该在学校大力推广，让更多同学在参与中受益。

通过这次数学建模学到了很多东西，首先，是从现实生活中发现问题，这就需要我们用心观察；然后就是解决问题的方法，由于我们不可能在课堂上学到所有知识，很多东西是要我们自学的，这就培养了我们自学能力，还有自己解决问题的能力：最后，因为模型要建立在真实数据上，就要求我们要有实事求是的态度了。数学建模虽然就告一段落了，学到的方法知识却伴随着我们以后的学习工作。“数学建模”课已使“数学建模”的核心思想深深在我意识中扎根，使我在今后学习中越来越善于发现问题并用数学知识创造性的去分析解决问题。数学建模课对学生思维能力的训练和思维方法上的引导，这就是它的主要魅力所在。一年一度的全国数学建模大赛在在9 月13 日上午8 点拉开战幕，各队 将在3 天72 小时内对一个现实中的实际问题进行模型建立，求解和分析，确定

题目后，我们队三人分头行动，一人去图书馆查阅资料，一人在网上搜索相关信 息，一人建立模型，通过三人的努力，在前两天中建立出两个模型并编程求解，经过艰苦的奋斗，终于在第三天完成了论文的写作，在这三天里我感触很深，现 将心得体会写出，希望与大家交流。1.团队精神：

团队精神是数学建模是否取得好成绩的最重要的因素，一队三个人要相互支 持，相互鼓励。切勿自己只管自己的一部分（数学好的只管建模，计算机好的只 管编程，写作好的只管论文写作），很多时候，一个人的思考是不全面的，只有

大家一起讨论才有可能把问题搞清楚，因此无论做任何板块，三个人要一起齐心 才行，只靠一个人的力量，要在三天之内写出一篇高水平的文章几乎是不可能的。2.有影响力的leader：

在比赛中，leader 是很重要的，他的作用就相当与计算机中的CPU，是全队 的核心，如果一个队的leader 不得力，往往影响一个队的正常发挥，就拿选题来说，有人想做A 题，有人想做B 题，如果争论一天都未确定方案的话，可能就没有足够时间完成一篇论文了，又比如，当队中有人信心动摇时（特别是第三天，人可能已经心力交瘁了），leader 应发挥其作用，让整个队伍重整信心，否则可能导致队伍的前功尽弃。3.合理的时间安排：

做任何事情，合理的时间安排非常重要，建模也是一样，事先要做好一个规 划，建模一共分十个板块（摘要，问题提出，模型假设，问题分析，模型假设，模型建立，模型求解，结果分析，模型的评价与推广，参考文献，附录）。你每 天要做完哪几个板块事先要确定好，这样做才会使自己游刃有余，保证在规定时 间内完成论文，以避免由于时间上的不妥，以致于最后无法完成论文。4.正确的论文格式：

论文属于科学性的文章，它有严格的书写格式规范，因此一篇好的论文一定 要有正确的格式，就拿摘要来说吧，它要包括6 要素（问题,方法,模型,算法,结论,特色），它是一篇论文的概括，摘要的好坏将决定你的论文是否吸引评委的目光，但听阅卷老师 说，这次有些论文的摘要里出现了大量的图表和程序，这都是不符合论文格式的，这种论

文也不会取得好成绩，因此我们写论文时要端正态度，注意书写格式。5.论文的写作：

我个人认为论文的写作是至关重要的，其实大家最后的模型和结果都差不

多，为什么有些队可以送全国，有些队可以拿省奖，而有些队却什么都拿不到，这关键在于论文的写作上面。一篇好的论文首先读上去便使人感到逻辑清晰，有 条例性，能打动评委；其次，论文在语言上的表述也很重要，要注意用词的准确 性；另外，一篇好的论文应有闪光点，有自己的特色，有自己的想法和思考在里 面，总之，论文写作的好坏将直接影响到成绩的优劣。

6.算法的设计：算法的设计的好坏将直接影响运算速度的快慢，建议大家多用数学软件（Mathematice,Matlab,Maple, Mathcad,Lindo,Lingo,SAS 等），这里提供十种数学以上便是我这次参加这次数学建模竞赛的一点心得体会，只当贻笑大方，不过就数学建模本身而言，它是魅力无穷的，它能够锻炼和考查一个人的综合素质，也希望广大同学能够积极参与到这项活动当中来。

第五篇：全国大学生数学建模竞赛

全国大学生数学建模竞赛

1、数模竞赛的起源与历史

数模竞赛是由美国工业与应用数学学会在1985年发起的一项大学生竞赛活动，目的在于激励学生学习数学的积极性，提高学生建立数学模型和运用计算机技术解决实际问题的综合能力，鼓励广大学生踊跃参加课外科技活动，开拓知识面，培养创精神及合作意识，推动大学数学教学体系、教学内容和方法的改革。我国大学生数学建模竞赛是由教育部高教司和中国工业与数学学会主办、面向全国高等院校的、每年一届的通讯竞赛。其宗旨是：创新意 识、团队精神、重在参与、公平竞争。1992载在中国创办，自从创办以来，得到了教育部高教司和中国工业与应用数学协会的得力支持和关心，呈现出迅速的发展发展势头，就2003年来说，报名阶段须然受到“非典”影响，但是全国30个省（市、自治区）及香港的637所院校就有5406队参赛，在职业技术学院增加更快，参赛高校由2002年的1067所上升到了2003年的1410所。可以说：数学建模已经成为全国高校规模最大课外科技活动。

2、什么是数学建模

数学建模（Mathematical Modelling）是一种数学的思考方法，是“对现实的现象通过心智活动构造出能抓住其重要且有

用的特征的表示，常常是形象化的或符号的表示。”从科学，工程，经济，管理等角度看数学建模就是用数学的语言和方法，通过抽象，简化建立能近似刻画并“解决”实际问题的一种强有力的数学工具。顾名思义，modelling一词在英文中有“塑造艺术”的意思，从而可以理解从不同的侧面，角度去考察问题就会有不尽的数学模型，从而数学建模 的创造又带有一定的艺术的特点。而数学建模最重要的特点是要接受实践的检验，多次修改模型渐趋完善的过程。

3、竞赛的内容：

竞赛题目一般来源于工程技术和管理科学等方面经过适当简化加工的实际问题，不要求参赛者预先掌握深入的专门知识，只需要学过普通高校的数学课程。题目有较大的灵活性供参赛者发挥其创造能力。参赛者应根据题目要求，完成一篇包括模型假设、建立和求解、计算方法的设计和计算机实现、结果的分析和检验、模型的改进等方面的论文（即答卷）。竞赛评奖以假设的合理性、建模的创造性、结果的正确性和文字表述的清晰程度为主要标准。

4、竞赛的步骤

建模是一种十分复杂的创造性劳动，现实世界中的事物形形色色，五花八门，不可能用一些条条框 框规定出各种模型如何具体建立，这里只是大致归纳一下建模的一般步骤和原则：

1）模型准备：首先要了解问题的实际背景，明确题目的要求，收集各种必要的信息．

2）模型假设：为了利用数学方法，通常要对问题做必要的、合理的假设，使问题的主要特征凸现出来，忽略问题的次要方面。

3）模型构成：根据所做的假设以及事物之间的联系，构造各种量之间的关系,把它问题化

4）模型求解：利用已知的数学方法来求解上一步所得到的数学问题，此时往往还要作出进一步的简化或假设。注意要尽量采用简单的数学工具。

5）模型分析：对所得到的解答进行分析，特别要注意当数据变化时所得结果是否稳定。

6）模型检验：分析所得结果的实际意义，与实际情况进行比较，看是否符合实际，如果不够理想，应该修改、补充假设，或重新建模，不断完善。

7）模型应用：所建立的模型必须在实际应用中才能产生效益，在应用中不断改进和完善。

5、模型的分类

按模型的应用领域分类： 生物数学模型、医学数学模型、地质数学模型、数量经济学模型、数学社会学模型

按是否考虑随机因素分类 ：确定性模型、随机性模型按是否考虑模型的变化分类 ：静态模型、动态模型按应用离散方法或连续方法 ：离散模型、连续模型

按建立模型的数学方法分类 ：几何模型、微分方程模型、图

论模型、规划论模型、马氏链模型 按人们对事物发展过程的了解程度分类 ：

白箱模型： 指那些内部规律比较清楚的模型。如力学、热学、电学以及相关的工程技术问题。

灰箱模型： 指那些内部规律尚不十分清楚，在建立和改善模型方面都还不同程度地有许多工作要做的问题。如气象学、生态学经济学等领域的模型。

黑箱模型：

指一些其内部规律还很少为人们所知的现象。如生命科学、社会科学等方面的问题。但由于因素众多、关系复杂，也可简化为灰箱模型来研究。

6、数学建模应用

今天，在国民经济和社会活动的以下诸多方面，数学建模都有着非常具体的应用。

1分析与设计： 例如描述药物浓度在人体内的变化规律以分析药物的疗效；建立跨音速空气流和激波的数学模型，用数值模拟设计新的飞机翼型。预报与决策： 生产过程中产品质量指标的预报、气象预报、人口预报、经济增长预报等等，都要有预报模型。使经济效益最大的价格策略、使费用最少的设备维修方案，是决策模型的例子。3 控制与优化： 电力、化工生产过程的最优控制、零件设计中的参数优化，要以数学模型为前提。建立大系统控制与优化的数

学模型，是迫切需要和十分棘手的课题。规划与管理 生产计划、资源配置、运输网络规划、水库优化调度，以及排队策略、物资管理等，都可以用运筹学模型解决 报名时间：从大赛的通知文稿发出后，就可以报名了，报名截止时间一般在开始比赛的前7到10天。

竞赛时间：每年的9月的第三个星期五上午8时至下一个星期一上午8时。

报名方式：如果有分赛区（每个赛区应至少有6所院校的20个队参加），就联系分赛区报名，没有分赛区，则直接向主委会报名。

大学生以队为单位参赛，每队3人（须属于同一所学校），专业不限。竞赛分本科、专科两组进行，本科生参加本科组竞赛，专科生参加专科组竞赛（也可参加本科组竞赛），研究生不得参加。每队可设一名指导教师（或教师组）。