2015数学建模选修大作业

合集下载

全国大学生数学建模2015年国二a题

全国大学生数学建模2015年国二a题

2015高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。

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

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

如有违反竞赛规则的行为,我们将受到严肃处理。

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

我们参赛选择的题号是(从A/B/C/D中选择一项填写): A我们的参赛报名号为(如果赛区设置报名号的话):所属学校(请填写完整的全名):参赛队员(打印并签名) :1.2.3.指导教师或指导教师组负责人(打印并签名):日期:年月日赛区评阅编号(由赛区组委会评阅前进行编号):编号专用页赛区评阅编号(由赛区组委会评阅前进行编号):全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):太阳影子定位摘要本文研究了太阳影子定位问题,基于天球坐标系相关知识、球面几何理论以及相似度理论,对不同情况下的数据,建立了相应的数学模型并得到了最优的匹配地点与日期。

问题1中,利用球面三角形余弦定理给出了太阳高度角公式,并建立了影子长度变化的数学模型,定性的分析了影子长度关于时角、当地纬度以及赤纬角的变化规律:(1). 时角的绝对值越大,影子长度越大;(2). 在同一经度上(即时角一定),当地纬度与此时的太阳赤纬之差越大,影子长度越大;(3). 在同一纬度不同经度上,当地经度和此时太阳直射点所在的经度之差越大,影子长度越大。

用所建的模型,得到了2015年10月22日北京时间9:00-15:00之间天安门广场3米高的直杆的太阳影子长度的变化曲线。

历年全国数学建模试题及解法归纳(2015)

历年全国数学建模试题及解法归纳(2015)

历年全国数学建模试题及解法归纳赛题解法93A非线性交调的频率设计拟合、规划93B足球队排名图论、层次分析、整数规划94A逢山开路图论、插值、动态规划94B锁具装箱问题图论、组合数学95A飞行管理问题非线性规划、线性规划95B天车与冶炼炉的作业调度动态规划、排队论、图论96A最优捕鱼策略微分方程、优化96B节水洗衣机非线性规划97A零件的参数设计非线性规划97B截断切割的最优排列随机模拟、图论98A一类投资组合问题多目标优化、非线性规划98B灾情巡视的最佳路线图论、组合优化99A自动化车床管理随机优化、计算机模拟99B钻井布局 0-1规划、图论00A DNA序列分类模式识别、Fisher判别、人工神经网络00B钢管订购和运输组合优化、运输问题01A血管三维重建曲线拟合、曲面重建赛题解法01B 公交车调度问题多目标规划02A车灯线光源的优化非线性规划02B彩票问题单目标决策03A SARS的传播微分方程、差分方程03B 露天矿生产的车辆安排整数规划、运输问题04A奥运会临时超市网点设计统计分析、数据处理、优化04B电力市场的输电阻塞管理数据拟合、优化05A长江水质的评价和预测预测评价、数据处理05B DVD在线租赁随机规划、整数规划06A出版社书号问题整数规划、数据处理、优化06B Hiv病毒问题线性规划、回归分析07A 人口问题微分方程、数据处理、优化07B 最佳交通线路查询多目标规划、图论08A 照相机问题非线性方程组、优化08B 大学学费问题数据收集和处理、统计分析、回归分析09A制动器试验台的控制方法分析物理模型,计算机仿真09B 眼科病房的合理安排综合评价,决策与预测10A储油罐的变位识别与罐容标定微积分理论,数值计算10B2010上海世博会影响力的评价综合评价,统计分析11A城市表层重金属污染分析综合评价,统计分析11B交巡警服务平台的设置与调度图论,动态规划12A葡萄酒的评价综合评价,统计分析12B太阳能小屋的设计多目标规划13A车道被占用对城市道路通行能力的影响交通流理论,排队论13B碎纸片的拼接复原算法14A嫦娥三号软着陆轨道设计与控制策略微分方程,最优化问题14B创意平板折叠桌微积分,几何赛题发展的特点:1.对选手的计算机能力提出了更高的要求:赛题的解决依赖计算机,题目的数据较多,手工计算不能完成,需要使用计算机软件。

数学建模选修大作业

数学建模选修大作业

中华女子学院成绩2014 — 2015学年第二学期期末考试(论文类)论文题目数学建模算法之蒙特卡罗算法课程代码 01课程名称数学建模学号 9姓名陈可心院系计算机系专业计算机科学与技术考试时间 2015年5月27日一、数学建模十大算法1、蒙特卡罗算法该算法又称随机性模拟算法,是通过计算机仿真来解决问题的算法,同时可以通过模拟可以来检验自己模型的正确性,是比赛时必用的方法。

接下来本文将着重介绍这一算法。

2、数据拟合、参数估计、插值等数据处理算法比赛中通常会遇到大量的数据需要处理,而处理数据的关键就在于这些算法,通常使用Matlab作为工具。

3、线性规划、整数规划、多元规划、二次规划等规划类问题建模竞赛大多数问题属于最优化问题,很多时候这些问题可以用数学规划算法来描述,通常使用Lindo、Lingo软件实现。

这个也是我们数学建模选修课时主要介绍的问题,所以对这方面比较熟悉,也了解了Lindo、Lingo软件的基本用法。

4、图论算法这类算法可以分为很多种,包括最短路、网络流、二分图等算法,涉及到图论的问题可以用这些方法解决,上学期数据结构课程以及离散数学课程中都有介绍。

它提供了对很多问题都很有效的一种简单而系统的建模方式。

5、动态规划、回溯搜索、分治算法、分支定界等计算机算法这些算法是算法设计中比较常用的方法,很多场合可以用到竞赛中6、最优化理论的三大非经典算法:模拟退火法、神经网络、遗传算法这些问题是用来解决一些较困难的最优化问题的算法,对于有些问题非常有帮助,但是算法的实现比较困难,需慎重使用。

7、网格算法和穷举法网格算法和穷举法都是暴力搜索最优点的算法,在很多竞赛题中有应用,当重点讨论模型本身而轻视算法的时候,可以使用这种暴力方案,最好使用一些高级语言作为编程工具。

8、一些连续离散化方法很多问题都是实际来的,数据可以是连续的,而计算机只认的是离散的数据,因此将其离散化后进行差分代替微分、求和代替积分等思想是非常重要的。

数学建模选修课作业

数学建模选修课作业
⑹Editor/Debugger M文件编辑调试器可以用来处理复杂的问题。
二﹑MATLAB解决线性代数问题
⑴①矩阵A的逆矩阵-0.1250 -0.1250 0.4000 -0.0500
0.1250 0.1250 0.2000 -0.1500
1.7500 0.7500 -1.0000 0.5000
1.5000 0.5000 -1.2000 0.4000
3.利用三维曲面绘图命令绘制3个三维曲面的图形
(1)绘制平面z=x+y的图形
程序:x=0:1:2;
y=0:1:2;
[x,y]=meshgrid(x,y)%生成xy平面上网格坐标矩阵
z=x+y;
axis([0,2,0,2,0,4])
mesh(x,y,z)%做三维曲面命令
(2)绘制螺旋抛物面z=x2+y2的图形
ans secx*tanx
(4)求导exyz,
程序:syms u x y z;
u=exp(x*y*z);
diff(diff(u,(x),y),z)
ans xexyz+x2zyexyz
(5)求导e2xy,
程序:syms u x y
u=exp(2*x*y)
diff(diff(u,(2x),y),z)
ans 2xexyz+x2zyexyz
-0.1616 - 1.4968i
1.0452
-1.0554
⑷利用MATLAB矩阵访问命令对下面的矩阵做初等行变换使其成为一个上三角矩阵
1 -2 2 -1 1
0 -4 8 0 2
0 0 -2 3 3
0 0 0 -6 4
三﹑MATLAB绘图绘制下面的曲线
1.在区间0≤x≤2中绘制cosx的图形

2015年全国大学生数学建模竞赛A题全国二等奖优秀论文设计

2015年全国大学生数学建模竞赛A题全国二等奖优秀论文设计

太阳影子定位摘要如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面,太阳影子定位技术就是通过分析视频中物体的太阳影子变化,确定视频拍摄的地点和日期的一种方法,该技术的日益成熟将有利于对视频中的场景进行大致定位和推算出拍摄时间。

可能会对部分案件的破解等事件产生极大的帮助。

为了更精确的计算视频中的拍摄地点和摄影时间,本文主要基于 MATLAB 与Excel处理软件,运用了遗传优化算法与模拟退火算法等,采用了视频数据化法、图片灰度化等处理手法,使计算更简便精确,使模型更完整可靠。

针对问题一,根据权威文献给出的太阳高度角算法建立模型一,先计算出太阳时角和太阳赤纬角后得到太阳高度角,再经过三角函数转换得到直杆的影长。

随后我们还考虑到因地球的大气状态并非真空状态会使到达地球的阳光折射,于是对太阳高度角进行了修正,使结果更加精确。

针对问题二,可以把这个问题当做是第一问的逆过程。

直杆影子的理论值与实际值的最小误差所对应的经纬度即为最优解。

在模型一的基础上,建立模型二并利用遗传算法计算此优化模型。

利用所给的21组坐标数据得到最优的直杆地点若干。

针对问题三,相较于问题二多了一个未知参数,在问题二的模型中加入这个未知参数即可得到模型三,得到最优的直杆地点与日期若干。

针对问题四,第一问中,利用 MATLAB 将视频每隔1min截取一张图片,把图片灰度化,测出影子、直杆底端与顶端的坐标,算得图中影长。

再根据已知图中影长、直杆实际长度与图中直杆长度的比例算出影长,运用模型二并进行优化后得出结果。

第二问中,运用模型三得到最优的视频的拍摄地点与日期若干,再进行优化得到最后结果关键词:遗传算法太阳高度角模拟退火算法最小二乘拟合问题粒子群算法1一、问题重述如何确定视频的拍摄地点和拍摄日期是视频数据分析技术的一个重要方面。

太阳影子定位技术就是通过分析视频中物体的太阳影子变化,确定视频拍摄的地点和日期的一种方法。

现需通过数学建模解决以下四个问题。

2015年全国大学生数学建模A题全国二等奖

2015年全国大学生数学建模A题全国二等奖

(1)
因为题目中所给的不是刚好的正午时刻,所以需要推导一般时间太阳高度角的计算 公式。假设某时刻,太阳没有处于当地的子午面上,此时太阳高度角 h 的值取决于此地 的纬度 、此时刻的太阳赤纬角 和太阳时角 t ' (如图 2) 。
图 2 文献中太阳高度角求解示意图
其中,S 为太阳移动的位置,MFD 为赤道截面, 太阳与赤道截面的夹角为太阳赤纬 角 ,即 FS , SN 90 。因为赤纬值日变化很小,因此一年中任意一天的赤纬 角 可由以下公式求得:
(6)
对于线性无关的方程组,当方程组个数多于未知数个数的时候,方程组就为超定方 程组,超定方程组理论上没有解满足所有方程,仅有最小二乘解。 对于本问题,需要找到三个参数使得代入数据之后其残差平方和最小,则有目标函 数:
min (
i
sin Li
2 L2 i h
1 sin
2

h Li
sin 1 sin
图 5 太阳方位角示意图
如果太阳方位角的第一种方法进行定义,求解较为繁琐,因此采用第二种定义方法 (如图 5) ,则方位角近似满足下公式:
sin sinh s sin cos As 0 cosh s cos
(5)
其中 cos As 中 As 为测量时的太阳方位角。因为题中所给的坐标系并未指明 x, y 正半 轴方向,所以我们假设测量坐标系的 y 正半轴和当地位置的正北方向夹角为 ,影子在 测量坐标系中与 y 正半轴的夹角为 As ' ,如下图(图 6)所示:
H
L
h



t
t'
N
3
五、模型的建立与求解

数学建模课后习题作业

数学建模课后习题作业

选修课——数学建模部分习题详细解答【陈文滨】1、在稳定的椅子问题中,如设椅子的四脚连线呈长方形,结论如何?【模型假设】(1)椅子四条腿一样长,椅脚与地面接触处视为一点,四脚的连线呈长方形.(2)地面高度是连续变化的,沿任何方向都不会出现间断 (没有像台阶那样的情况),即从数学的角度看,地面是连续曲面.这个假设相当于给出了椅子能放稳的必要条件.(3)椅子在任何位置至少有三只脚同时着地.为保证这一点,要求对于椅脚的间距和椅腿的长度而言,地面是相对平坦的.因为在地面上与椅脚间距和椅腿长度的尺寸大小相当的范围内,如果出现深沟或凸峰(即使是连续变化的),此时三只脚是无法同时着地的。

【模型建立】在上述假设下,解决问题的关键在于选择合适的变量,把椅子四只脚同时着地表示出来.首先,引入合适的变量来表示椅子位置的挪动.生活经验告诉我们,要把椅子通过挪动放稳,通常有拖动或转动椅子两种办法,也就是数学上所说的平移与旋转变换.然而,平移椅子后问题的条件没有发生本质变化,所以用平移的办法是不能解决问题的.于是可尝试将椅子就地旋转,并试图在旋转过程中找到一种椅子能放稳的情形.注意到椅脚连线呈长方形,长方形是中心对称图形,绕它的对称中心旋转180度后,椅子仍在原地.把长方形绕它的对称中心O旋转,这可以表示椅子位置的改变。

于是,旋转角度θ这一变量就表示了椅子的位置.为此,在平面上建立直角坐标系来解决问题.如下图所示,设椅脚连线为长方形ABCD,以对角线AC所在的直线为x轴,对称中心O为原点,建立平面直角坐标系.椅子绕O点沿逆时针方向旋转角度θ后,长方形ABCD转至A1B1C1D1 的位置,这样就可以用旋转角θ(0≤θ≤π)表示出椅子绕点O旋转θ后的位置.其次,把椅脚是否着地用数学形式表示出来.我们知道,当椅脚与地面的竖直距离为零时,椅脚就着地了,而当这个距离大于零时,椅脚不着地.由于椅子在不同的位置是θ的函数,因此,椅脚与地面的竖直距离也是θ的函数.由于椅子有四只脚,因而椅脚与地面的竖直距离有四个,它们都是θ的函数.而由假设(3)可知,椅子在任何位置至少有三只脚同时着地,即这四个函数对于任意的θ,其函数值至少有三个同时为0.因此,只需引入两个距离函数即可.考虑到长方形ABCD是中心对称图形,绕其对称中心 O沿逆时针方向旋转180°后,长方形位置不变,但A,C和B,D对换了.因此,记A、B两脚与地面竖直距离之和为f(θ),C、D两脚与地面竖直距离之和为g(θ),其中θ∈[0,π],从而将原问题数学化。

数学建模大作业习题答案

数学建模大作业习题答案

数学建模大作业习题答案数学建模大作业习题答案作为一门应用数学课程,数学建模在现代科学研究和工程技术中具有重要的地位和作用。

通过数学建模,我们可以将实际问题转化为数学模型,从而利用数学方法进行分析和求解。

在数学建模的学习过程中,我们经常会遇到一些习题,下面我将为大家提供一些数学建模大作业题目的答案,希望能对大家的学习有所帮助。

1. 题目:某城市的交通拥堵问题解答:针对这个问题,我们可以采用图论的方法进行建模和求解。

首先,我们将城市的道路网络抽象为一个图,图的节点表示交叉口,边表示道路。

然后,我们可以给每条边赋予一个权重,表示道路的通行能力。

接着,我们可以使用最短路径算法,比如Dijkstra算法,来计算从一个交叉口到另一个交叉口的最短路径,从而找到最优的交通路线。

此外,我们还可以使用最小生成树算法,比如Prim算法,来构建一个最小的道路网络,以减少交通拥堵。

2. 题目:某工厂的生产调度问题解答:对于这个问题,我们可以采用线性规划的方法进行建模和求解。

首先,我们可以将工厂的生产任务抽象为一个线性规划模型,其中目标函数表示最大化生产效益,约束条件表示生产能力、物料供应和市场需求等方面的限制。

然后,我们可以使用线性规划求解器,比如Simplex算法或内点法,来求解这个线性规划模型,得到最优的生产调度方案。

此外,我们还可以引入一些启发式算法,比如遗传算法或模拟退火算法,来寻找更好的解决方案。

3. 题目:某股票的价格预测问题解答:对于这个问题,我们可以采用时间序列分析的方法进行建模和求解。

首先,我们可以将股票的价格序列抽象为一个时间序列模型,比如ARIMA模型。

然后,我们可以使用历史数据来拟合这个时间序列模型,并进行参数估计。

接着,我们可以利用这个时间序列模型来预测未来的股票价格。

此外,我们还可以引入其他的预测方法,比如神经网络或支持向量机,来提高预测的准确性。

通过以上的例子,我们可以看到,在数学建模的过程中,我们需要将实际问题抽象为数学模型,然后利用数学方法进行分析和求解。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

中华女子学院成绩2014 — 2015学年第二学期期末考试(论文类)论文题目数学建模算法之蒙特卡罗算法课程代码1077080001课程名称数学建模学号*********姓名陈可心院系计算机系专业计算机科学与技术考试时间2015年5月27日一、数学建模十大算法1、蒙特卡罗算法该算法又称随机性模拟算法,是通过计算机仿真来解决问题的算法,同时可以通过模拟可以来检验自己模型的正确性,是比赛时必用的方法。

接下来本文将着重介绍这一算法。

2、数据拟合、参数估计、插值等数据处理算法比赛中通常会遇到大量的数据需要处理,而处理数据的关键就在于这些算法,通常使用Matlab作为工具。

3、线性规划、整数规划、多元规划、二次规划等规划类问题建模竞赛大多数问题属于最优化问题,很多时候这些问题可以用数学规划算法来描述,通常使用Lindo、Lingo软件实现。

这个也是我们数学建模选修课时主要介绍的问题,所以对这方面比较熟悉,也了解了Lindo、Lingo软件的基本用法。

4、图论算法这类算法可以分为很多种,包括最短路、网络流、二分图等算法,涉及到图论的问题可以用这些方法解决,上学期数据结构课程以及离散数学课程中都有介绍。

它提供了对很多问题都很有效的一种简单而系统的建模方式。

5、动态规划、回溯搜索、分治算法、分支定界等计算机算法这些算法是算法设计中比较常用的方法,很多场合可以用到竞赛中6、最优化理论的三大非经典算法:模拟退火法、神经网络、遗传算法这些问题是用来解决一些较困难的最优化问题的算法,对于有些问题非常有帮助,但是算法的实现比较困难,需慎重使用。

7、网格算法和穷举法网格算法和穷举法都是暴力搜索最优点的算法,在很多竞赛题中有应用,当重点讨论模型本身而轻视算法的时候,可以使用这种暴力方案,最好使用一些高级语言作为编程工具。

8、一些连续离散化方法很多问题都是实际来的,数据可以是连续的,而计算机只认的是离散的数据,因此将其离散化后进行差分代替微分、求和代替积分等思想是非常重要的。

9、数值分析算法如果在比赛中采用高级语言进行编程的话,那一些数值分析中常用的算法比如方程组求解、矩阵运算、函数积分等算法就需要额外编写库函数进行调用。

10、图象处理算法赛题中有一类问题与图形有关,即使与图形无关,论文中也应该要不乏图片的,这些图形如何展示以及如何处理就是需要解决的问题,通常使用Matlab进行处理。

二、蒙特卡罗方法2.1算法简介蒙特·卡罗方法(Monte Carlo method),也称统计模拟方法,1946年,美国拉斯阿莫斯国家实验室的三位科学家John von Neumann,Stan Ulam 和 NickMetropolis共同发明了,蒙特卡罗方法。

此算法被评为20世纪最伟大的十大算法之一。

是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明,而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。

是指使用随机数来解决很多计算问题的方法。

由于传统的经验方法由于不能逼近真实的物理过程,很难得到满意的结果,而蒙特卡罗方法由于能够真实地模拟实际物理过程,故解决问题与实际非常符合,可以得到很圆满的结果。

与它对应的是确定性算法。

蒙特·卡罗方法在金融工程学,宏观经济学,计算物理学(如粒子输运计算、量子热力学计算、空气动力学计算)等领域应用广泛。

2.2蒙特卡罗方法的特点蒙特卡罗方法通过抓住事物运动的几何数量和几何特征,利用数学方法来加以模拟,即进行一种数字模拟实验。

它是以一个概率模型为基础,按照这个模型所描绘的过程,通过模拟实验的结果,作为问题的近似解。

蒙特卡罗方法与一般计算方法有很大区别,一般计算方法对于解决多维或因素复杂的问题非常困难,而蒙特卡罗方法对于解决这方面的问题却比较简单。

其特点如下:1、直接追踪粒子,物理思路清晰,易于理解。

2、采用随机抽样的方法,较真切的模拟粒子输运的过程,反映了统计涨落的规律。

3、不受系统多维、多因素等复杂性的限制,是解决复杂系统粒子输运问题的好方法。

蒙特卡罗方法的基本原理及思想如下:当所要求解的问题是某种事件出现的概率,或者是某个随机变量的期望值时,它们可以通过某种“试验”的方法,得到这种事件出现的频率,或者这个随机变数的平均值,并用它们作为问题的解。

这就是蒙特卡罗方法的基本思想。

蒙特卡罗方法通过抓住事物运动的几何数量和几何特征,利用数学方法来加以模拟,即进行一种数字模拟实验。

它是以一个概率模型为基础,按照这个模型所描绘的过程,通过模拟实验的结果,作为问题的近似解。

2.3适用模型假设我们要计算一个不规则图形的面积,那么图形的不规则程度和分析性计算(比如,积分)的复杂程度是成正比的。

蒙特卡洛方法是怎么计算的呢?,我们可以想象把图形画在一块方形布上,然后找来一袋豆子,然后将所有豆子洒在布上,落在图形内豆子的重量比上那块布上所有豆子的重量再乘以布的面积就是他所要求的图形的面积。

这确实是一个求面积的好方法,将整个坐标轴看成一个固定的面积,然后均匀的这个分成N(N的大小取决于划分的步长)个点,然后找出N个点中有多少个点是属于阴影部分中,假设这个值为k,则阴影部分的面积就求出来了。

此方法是利用蒙特卡罗方法计算阴影部分面积,是把豆子均匀分布在布上;就计算结果的精度而言,取决点的分割是否够密,即N是否够大;在数值积分法中,利用求单位圆的1/4的面积来求得Pi/4从而得到Pi。

单位圆的1/4面积是一个扇形,它是边长为1单位正方形的一部分。

只要能求出扇形面积S1在正方形面积S中占的比例K=S1/S就立即能得到S1,从而得到Pi的值。

怎样求出扇形面积在正方形面积中占的比例K呢?一个办法是在正方形中随机投入很多点,使所投的点落在正方形中每一个位置的机会相等看其中有多少个点落在扇形内。

将落在扇形内的点数m与所投点的总数n的比m/n作为k的近似值。

P落在扇形内的充要条件是221 x y+≤。

已知:K=1ss,K≈mn,s=1,s1=4Pi,求Pi。

num++; /*统计落在四分之一圆之内的点数*/}printf("Pi值等于:%f\n",num*4.0/COUNT);printf("RAND_MAX=%d\n",RAND_MAX);}结果:测试6次的结果显示:(可以看出: 随着点数的增加,求得的Pi值渐渐接近真实值。

)此外,蒙特·卡罗方法在金融工程学,宏观经济学,计算物理学(如粒子输运计算、量子热力学计算、空气动力学计算)等领域应用广泛。

2.4算法应用实例例:在我方某前沿防守地域,敌人以一个炮排(含两门火炮)为单位对我方进行干扰和破坏.为躲避我方打击,敌方对其阵地进行了伪装并经常变换射击地点.经过长期观察发现,我方指挥所对敌方目标的指示有50%是准确的,而我方火力单位,在指示正确时,有1/3的概率能毁伤敌人一门火炮,有1/6的概率能全部消灭敌人.现在希望能用某种方式把我方将要对敌人实施的1次打击结果显现出来,利用频率稳定性,确定有效射击(毁伤一门炮或全部消灭)的概率.分析:这是一个复杂概率问题,可以通过理论计算得到相应的概率. 为了直观地显示我方射击的过程,现采用模拟的方式。

1. 问题分析需要模拟出以下两件事:[1]观察所对目标的指示正确与否模拟试验有两种结果,每一种结果出现的概率都是1/2。

因此,可用投掷一枚硬币的方式予以确定,当硬币出现正面时为指示正确,反之为不正确。

[2] 当指示正确时,我方火力单位的射击结果情况。

模拟试验有三种结果:毁伤一门火炮的可能性为1/3,毁伤两门的可能性为1/6,没能毁伤敌火炮的可能性为1/2。

这时可用投掷骰子的方法来确定:如果出现的是1、2、3三个点:则认为没能击中敌人;如果出现的是4、5点:则认为毁伤敌人一门火炮;若出现的是6点:则认为毁伤敌人两门火炮。

2. 符号假设i:要模拟的打击次数;k1:没击中敌人火炮的射击总数;k2:击中敌人一门火炮的射击总数;k3:击中敌人两门火炮的射击总数;E:有效射击(毁伤一门炮或两门炮)的概率;3. 在Matlab中编辑:function liti6(p,mm)efreq=zeros(1,mm);randnum1 = binornd(1,p,1,mm);randnum2 = unidrnd(6,1,mm);k1=0;k2=0;k3=0;for i=1:mmif randnum1(i)==0k1=k1+1;elseif randnum2(i)<=3k1=k1+1;else if randnum2(i)==6k3=k3+1;elsek2=k2+1;endendefreq(i)=(k2+k3)/i;endnum=1:mm;plot(num,efreq)4.在Matlab命令行中输入以下命令:liti6(0.5,2000)liti6(0.5,20000)5.理论计算设:⎩⎨⎧=观察所对目标指示正确确观察所对目标指示不正10j A 0:射中敌方火炮的事件;A 1:射中敌方一门火炮的事件;A 2:射中敌方两门火炮的事件.则由全概率公式:P(A 0) = P(j=0)P(A 0∣j=0) + P(j=1)P(A 0∣j=1) = 25.02121021=⨯+⨯ P(A 1) = P(j=0)P(A 1∣j=0) + P(j=1)P(A 1∣j=1) =613121021=⨯+⨯ P(A 2) = P(j=0)P(A 2∣j=0) + P(j=1)P(A 2∣j=1) = 1216121021=⨯+⨯ E= 25.012161=+6. 结果比较模拟结果与理论计算近似一致,能更加真实地表达实际战斗动态过程.三、思考和体会它所教给我们的不单是一些数学方面的知识,更多的其实是综合能力的培养、锻炼与提高。

它培养了我们全面、多角度考虑问题的能力,使我们的逻辑推理能力和量化分析能力得到很好的锻炼和提高。

它还让我了解了多种数学软件,以及运用数学软件对模型进行求解。

其实,数学建模对我们来说并不陌生,在我们的日常生活和工作中,经常会用到有关建模的概念。

例如,我们平时出远门,会考虑一下出行的路线,以达到既快速又经济的目的;一些工厂为了获得更大的利润,往往会策划出一个合理安排生产和销售的最优方案……这些问题和建模都有着很大的联系。

数学建模所要解决的问题决不是单一学科问题,它除了要求我们有扎实的数学知识外,还需要我们不停地去学习和查阅资料,除了我们要学习许多数学分支问题外,还要了解工厂生产、经济投资、社会生活等方面的知识,这些知识决不是任何专业中都能涉猎得到的。

它能极大地拓宽和丰富我们的内涵,让我们感到了知识的重要性,这些知识必将为我们将来的学习工作打下坚实的基础。

相关文档
最新文档