《第12讲 计算机模拟》由会员分享,可在线阅读,更多相关《第12讲 计算机模拟(37页珍藏版)》请在人人文库网上搜索。
1、学习数学建模和数学实验、物流工程学院数学教室、修正机模拟、实验目的、实验内容、修正机模拟的基本过程和方法。 1、模拟的概念。 4、实验工作。 3、计算机模拟实例。 2 .用于生成随机数的校正器命令。连续系统仿真实例3360跟踪问题、离散系统仿真实例3360矩阵问题、用蒙特卡罗法解非线性修正图像问题,来进行修正计算机仿真实例、仿真概念、仿真利用物理、数学的模型表现仿真的基本思想包括研究系统的主要特征通过运行该实验模型,研究系统的必要信息,仿真方法,1,物理仿真:用功能相似的实物系统模拟实际系统及其过程。 例如军事演习、船艇实验、沙盘作业等。 中的组合图层性质变更选项。 物理模拟通常成本昂贵,周期。
2、长,很难在物理模型上更改系统的结构和系数。 另外,很多系统不能进行社会经济系统、生态系统等的物理模拟。 在实际问题中,面对具有一些随机因素的复杂系统,用分析方法建模需要创造许多简化假说,可能与面临的实际问题相差甚远,答案根本无法应用。 在此情况下,校正算机模拟几乎是唯一的选择。 在一定的假设条件下,用数学运算模拟系统的运行,称为数学模拟。 现代数学模拟全部在校正机上进行,称为校正机模拟。 2、数学模拟、校正算机模拟可以重复进行,改变系统的结构和系数比较容易。 蒙特卡洛方法是一种用随机数进行修正机仿真的方法,该方法对研究系统随机观察采样,通过样本值的观察系统修正,求研究系统的一些残奥仪的敌人伪装。
3、其阵地,始终反射我方指挥所对敌人目标的指示50正确,我方火力机关指示正确时,1/3的射击效果损伤敌人一个火炮,1/6的射击效果全部可以消灭敌人,现在想以某种方式表现出我方对敌人实施的20次打击结果。 分析:这是一个概率问题,通过理论修正算法可以得到相应的概率和期望值。 但是,这只是作战行动的最终静态结果,不能显示作战行动的动态过程。 为了显示我方20次射击的过程,现在采用模拟方式。1 .问题分析,2指示正确的话,我方火力单位的射击结果状况,1观察对象的指示是否正确,模拟有2个结果,各自的结果出现的概率是1/2,所以用投1枚硬币的方式确定, 如果硬币出现在正面,指示正确的模拟有三个结果:损伤一个。
4、火炮的可能性是1/3 (即2/6 ),损伤两个火炮的可能性是1/6,不能损伤敌火炮的可能性是1/2 (即3/6 ),此时掷骰子的方法可以确定:如果2 .记号假说,I :模拟的打击次数k1:敌人的火炮击中的射击总数k2:敌人的火炮击中的射击总数k3:敌人的两个火炮击中的射击总数e :有效射击比率E1:20次射击的平均每次伤害敌人的火炮数,3 .模拟6、结果比较、归纳,模拟结果与理论修正算法虽然不完全一致,但它能更真实地表现实际战斗动态过程,用蒙特卡罗方法修正算法模拟的步骤: 2伪随机数现象是借助于具有不同概率分布的伪随机数来伪随机数现象并生成伪随机数的校正器指令,其中在Matlab软件中可直接生。
5、成满足不同分布的随机数,其包括:2mn阶,均匀分布的随机数矩阵: rand (随机数) 一个n )产生均匀分布的b )随机数矩阵: unifrnd (a,b,m,n )产生a,b均匀分布的随机数: unifrnd (a,b ),并且如果知道只有一个随机数在(a,b )内,则例1的校正器模拟研究对象被视为相互独立的多个随机变量的和,在各个变量对总和的影响小的情况下,认为其对象服从正态分布。 机械加工得到的零件尺寸的偏差,射击着落点和目标的偏差,各种测定误差,人的身高,体重等,可以看作几乎服从正态分布。 当连续型随机变量x的概率密度函数把其中的0设为常数时,x被称为遵循残奥仪表的指数分布。 指数分。
6、布的期待值是排队服务系统中的顾客到达率为常数时的到达间隔,故障率为常数时的部件寿命服从指数分布。 指数分布广泛应用于排队论、可靠性分析。 注意:在Matlab中,用于生成残奥仪表的指数分布的命令是exprnd (),例如,根据残奥仪表具有0.1的指数分布,顾客到达某个店铺的间隔时间,指数分布的平均值是1/0.1=10。 两位顾客到达店铺的平均间隔时间为10单位时间,也就是平均10单位时间到达一位顾客。 顾客到达的间隔时间可以用exprnd(10 )模拟。 设离散型随机变量x的所有可能值为0、1、2,取各自值的概率为0,则x服从残奥仪表的个人分布。 个人计算机分布广泛应用于排队系统、产品检测、天。
7、文、物理等领域。 个人分布的期望值是指数分布和个人分布的关系: (1)2个顾客到达店铺的平均间隔时间为10单位时间,即平均10单位时间到达1个顾客;(2)1个单位时间内平均到达0.1个顾客的例子(1) 顾客到达某店铺的间隔时间服从残奥仪表为0.1的指数分布(2)该店铺在单位时间内到达的顾客数服从残奥仪表为0.1的个人分布,返回,例2敌坦克分队对我方阵地实施袭击,其到达规则服从泊松分布用(poissrnd(4)进行模拟。To Matlab(poiss )、(2)坦克到达的间隔时间按照残奥仪表4的负指数分布,用exprnd(1/4)进行模拟。To Matlab(time )、连续系统仿真实例:跟踪。
8、该问题,并且将状态随时连续变化的系统称为连续系统。 对连续系统的校正算法仿真只能进行近似,只要该近似达到一定的精度,就可以满足要求。 例如追踪问题:如图所示,正方形ABCD的4个顶点分别是一个人。在某个时刻,4个人同时出发,以等速v=1m/秒按时间修正的方式追踪下一个人,如果他们总是对准目标,最终以螺旋状曲线在中心点o .求这种情况下的各个人的前进方向y4),2 .将时间间隔设为t,校正各点在各点的坐标,4 .针对各点,求出在各点的位置连结、即运动轨迹。 dt=0.05; x=00十十十; y=0到100; 最后一次打印=133604打印(y(i ),y(i ), ),保持。 威尔(d0.1)。
9、 x (5)=x (1)。 y (5)=y (1) :不等于133604 d=sqrt (x (I1)-x (I ) )2(y (I1)-y (I ) ) 2。 x(i)=x(i) v*dt*(x(i 1)-x(i)/d y(i)=y(i) v*dt*(y(i 1)-y(i)/d 在plot(x(i ),y(i ),),hold on end,校正算法程序3360,To Matlab(chase ),返回,离散系统模拟实例3360队列问题,数组,排队系统中,服务对象排队论通过对各个随机服务现象的统一修正研究,找到反映这些随机现象平均特性的规则,为新服务系统的设置修正和现有服务系统工作的改进提供。
10、依据。 对于排队服务系统,客户关心排队的人是否太多,排队的时间是否太长,服务员关心他的空闲时间是否太短。 因此,人们用排队长度、排队时间、服务利用率等指标来测量系统性能。 (2)行列的长度没有限制(3)到达系统的顾客依次进入服务,即“先到先服务”。 单服务员排队模式:某店店员陆续来,店员陆续接待客人的客人很多时,有些客人就要排队等待。 接待客户离开店铺设定: 1客户的到达间隔时间按照残奥表为0.1的指数分布,按照客户的服务时间,模拟11个工作日内完成的服务个数和客户的平均等待时间t 2,模拟100个工作日,平均每天完成的服务个数和每天的客户平均等待2符号说明w :总等待时间ci :第I位顾客的。
11、到达时刻bi :第I位顾客开始服务的时刻ei :第I位顾客服务结束时刻xi :第i-1位顾客与第I位顾客之间的到达间隔时间yi :第I位顾客的服务时间、c1、b1、c3 c5 1一天模拟To Matlab(simu1),100天模拟To Matlab(simu2),于是用蒙特卡洛法解决非线性修正问题,基本假设,试验点的第j个成分xj遵循aj,bj内的均匀分布,符号假设,求解过程的其他成分原封不动每次制作新的试验点时,新的随机数成分在KMAXK或者PMAXP时停止反复,在Matlab软件包中编程。 需要合并的三个文件:randlp.m、mylp.m、lpconst.m .的%被转换为求最小值的问题,% randlp.m function sol,r1,r2=randlp(a,b,n) %随机模拟a=0; %试验点下界b=10; %试验点上界n=1000; %试验点的个数r1=统一nd (a,b,n,1 ); % n一阶的a,b均匀分布随机数列r2=统一nd (a,b,n,1 ); sol=r1(1) r2(1): z0=inf; I=13360 n x1=r1(I )。 x2=r2(I ):LPC=LP const (x1x2):if LPC=1z=mylp (x1x2):if z。