计算机模拟课件.pptVIP

2. 符号假设 i：要模拟的打击次数； k1：没击中敌人火炮的射击总数； k2：击中敌人一门火炮的射击总数；k3：击中敌人两门火炮的射击总数． E：有效射击比率； E1：20次射击平均每次毁伤敌人的火炮数． 3. 模拟框图 初始化:i=0, k1=0, k2=0, k3=0 i=i+1 骰子点数? k1=k1+1 k2=k2+1 k3=k3+1 k1=k1+1 i＜20? E= E1= 0× +1 × +2 × 停止 硬币正面? Y N N Y 1，2，3 4，5 6 4. 模拟结果 6. 结果比较 虽然模拟结果与理论计算不完全一致，但它却能更加真实地表达实际战斗动态过程． 用蒙特卡罗方法进行计算机模拟的步骤： [1] 设计一个逻辑框图，即模拟模型．这个框图要正确反映系统各部分运行时的逻辑关系． [2] 模拟随机现象．可通过具有各种概率分布的模拟随机数来模拟随机现象． clear k1=0;k2=0;k3=0; for i=1:20 R1=rand if R1=0.5 R2=rand if R2=3/6 k1=k1+1; else if R25/6 k3=k3+1; else k2=k2+1; end end else k1=k1+1; end end i E=(k2+k3)/20 E1=(0*k1+1*k2+2*k3)/20 离散系统模拟实例: 排队问题 排队论主要研究随机服务系统的工作过程． 在排队系统中，服务对象的到达时间和服务时间都是随机的．排队论通过对每个个别的随机服务现象的统计研究，找出反映这些随机现象平均特性的规律，从而为设计新的服务系统和改进现有服务系统的工作提供依据． 对于排队服务系统, 顾客常常注意排队的人是否太多, 等候的时间是否长, 而服务员则关心他空闲的时间是否太短. 于是人们常用排队的长度、等待的时间及服务利用率等指标来衡量系统的性能. [1] 系统的假设： （1） 顾客源是无穷的； （2） 排队的长度没有限制； （ 3） 到达系统的顾客按先后顺序依次进入服务， 即“先到先服务”． 单服务员的排队模型：在某商店有一个售货员，顾客陆续来到，售货员逐个地接待顾客．当到来的顾客较多时，一部分顾客便须排队等待，被接待后的顾客便离开商店．设： 1．顾客到来间隔时间服从参数为0.1的指数分布． ２．对顾客的服务时间服从［４，１５］上的均匀分布． ３．排队按先到先服务规则，队长无限制． 假定一个工作日为8小时，时间以分钟为单位． [1]模拟一个工作日内完成服务的个数及顾客平均等待时间t． [2]模拟100个工作日，求出平均每日完成服务的个数及每日顾客的平均等待时间． [2] 符号说明 w：总等待时间；ci：第i个顾客的到达时刻； 　　 bi：第i个顾客开始服务时刻； ei：第i个顾客服务结束时刻． xi：第i-1个顾客与第i个顾客到达之间的时间间隔 yi：对第i个顾客的服务时间 c1 b1 c3 c4 c5 c2 e1 b2 e2 b3 e3 b4 e4 b5 ci=ci-1+ xi ei=bi+yi bi=max(ci,ei-1) t [3] 模拟框图 初始化：令i=1，ei-1=0，w=0 产生间隔时间随机数xi服从参数为0.1的指数分布 ci=xi , bi=xi 产生服务时间随机数yi服从[4,15]的均匀分布 ei=bi+yi 累计等待时间：w=w+bi-ci 准备下一次服务：i=i+1 产生间隔时间随机数xi服从参数为0.1的指数分布 ci=ci-1+ xi 确定开始服务时间：bi=max(ci,ei-1) bi＞480? Y N i=i-1,t=w/i 输出结果：完成服务个数：m=i 平均等待时间：t 停止 [1]模拟1日 To MATLAB(simu1) [2]模拟100日 To MATLAB(simu2) clear i=2; w=0; x(i)=exprnd(10); c(i)=x(i); b(i)=x(i); while b(i)=480 y(i)=unifrnd(4,15); e(i)=b(i)+y(i