2017数学建模计算机仿真解决方案.pptVIP

计算机仿真 一个问题 我们做一个实验：把一个硬币掷一万次，统计两个面出现的次数。这样做很简单但却需要大量时间，有没有一种较快的办法把这个实验完成呢？ 利用计算机可以实现这一想法 生成一个在 [ 0, 1] 中的随机数a， 如果a0.5，则认为是掷硬币出现了正面，给计数变量k1增加1； 如果，则认为是掷硬币出现了反面，给计数变量k2增加1。 将该过程循环一万次即可。 计算机仿真的定义 计算机仿真就是根据已知的信息和知识，利用计算机模拟现实情况或系统演变过程，发现新的知识和规律，从而解决问题的一种方法。 计算机仿真被称为独立于理论研究和实验研究的第三种方法。 计算机仿真的特点 代价小，时间短，可重复，参数设置灵活 是一种独特的“数”学模型。 是一种求解许多实际问题和数学模型的简单方法，由于它不需要太多的数学知识，非常适合各类工程技术人员。 计算机仿真仿的是“象”、是“数”，要忽略许多具体的事物特征。 如何把计算机仿真的过程作为一个“数学模型”表述出来呢？ 描述计算机仿真模型要包括两个内容，一是对系统关键数据计算方法的清晰表述，二是对仿真的程序流程的描述，可以用算法步骤的形式，也可以用算法流程图。 计算机仿真要靠一个计算机程序来实现，然而程序代码是不能作为模型，而且由于选择的系统语言不同，表述上也会有较大差异。 计算机仿真的分类 物理系统仿真——电系统、机械系统等的仿真，大坝承受力仿真，原子弹爆炸威力…… 系统演变仿真——河堤垮塌后洪水蔓延程度仿真，海啸蔓延程度仿真，种群生长仿真，战争推演仿真…… 蒙特卡洛方法——不规则图形的面积、体积，圆周率的计算，电脑围棋……核心：生成随机数，……被列为20世纪最伟大的10大算法之首。 离散事件仿真——企业经营策略…… 有些仿真需要一些设备工具甚至人的参与，这里不涉及此类，只考虑与完全可用数学推演描述的问题。 计算机仿真的核心思想方法 过程明确，机理清晰 连续问题离散化 蒙特卡洛方法 遍历 计算机仿真案例2 例2 （赶火车过程仿真）一列火车从A站经过B站开往C站，某人每天赶往B站乘这趟火车。已知火车从A站到B站的运行时间是均值为30min、标准差为2min的正态随机变量。火车大约在下午1点离开Ａ站。火车离开时刻的频率分布和这个人到达Ｂ站时刻的频率分布如下表所示。问他能赶上火车的概率有多大？ 计算机仿真案例4 某自行车商店的仓库管理人员采取一种简单的订货策略，当库存量降低到P辆自行车时就向厂家订货，每次订货Q辆，如果某一天的需求量越过了库存量，商店就有销售损失和信誉损失，但如果库存量过多，会导致资金积压和保管费增加。该问题的已知条件是： （1）从发出订货到收到货物需隔三天； （2）每辆自行车保管费为0.75元/天，每辆自行车的缺货损失为1.8元/天，每次的订货费为75元； （3）每天自行车的需求量服从0到99之间的均匀分布； （4）原始库存为115辆，并假设第一天没有发出订货。 若现在已有如下表所示的五种库存策略，请选择一种总费用最少的策略。 评价函数的设置方法 第i天销售量为 辆 缺货量为 辆 销售完后的库存量为 辆 0-1变量 表示第i天销售完后是否有库存，若库存量大于等于0，则 ＝1，否则为0； 用0-1变量 表示第i天销售是否会缺货，如果缺货量大于0，则 ＝1； 用0-1变量 表示第i天是否要订货。 * * 上面就是一个计算机仿真最简单的例子！ 可以说计算机仿真的适用于几乎所有的社会生活领域！ 产生模拟随机数的计算机命令 在MATLAB软件中，可以直接产生满足各种分布的随机数，命令如下： 2．产生m×n阶［０，１］均匀分布的随机数矩阵： rand (m, n) 产生一个［０，１］均匀分布的随机数：rand 1．产生m×n阶［a，b］上均匀分布U（a，b）的随机数矩阵： unifrnd (a,b,m, n) 产生一个［a，b］均匀分布的随机数：unifrnd (a,b) 当只知道一个随机变量取值在（a，b）内，但不知道（也没理由假设）它在何处取值的概率大，在何处取值的概率小，就只好用U（a，b）来模拟它． 当研究对象视为大量相互独立的随机变量之和，且其中每一种变量对总和的影响都很小时，可以认为该对象服从正态分布． 机械加工得到的零件尺寸的偏差、射击命中点与目标的偏差、各种测量误差、人的身高、体重等，都可近似看成服从正态分布． 若连续型随机变量X的概率密度函数为 其中 0为常数，则称X服从参数为 的指数分布． 指数分布的期望值为 排队服务系统中顾客到达率为常数时的到达间隔、故障率为常数时零件的寿命都服从指数分布． 指数分