第六章 离散系统仿真3(事件步长法)

第六章事件步长法一．事件步长法的一般方法事件步长法是以时间为增量，按照时间的进展，一步一步地对系统的行为进行仿真，直到预定的时间结束为止事件步长法与时间步长法的主要区别在于：1。在时间步长法中仿真时钟等步长的前进，而在事件步长法中，仿真时钟的步长由事件之间的时间间隔决定。2时间步长法在一个步长内认为系统所有的状态是相同的，因此所选的步长的大小将会影响到仿真的精度，而在事件步长法中，每一个事件发生都有确切的时间，步长的大小对他的影响很小。3。时间步长法每前进一步就要扫描系统的所有的状态，事件步长法只在事件的发生点上扫描系统事件步长法流程图二．事件表法的基本思路事件表法的主要思路是将系统的仿真过程看成是一个事件点序列。根据事件出现的时序，用一个称为时间表的表格来调度事件执行的顺序。将当前要处理的事件列入表里，从表里取出最接近的事件处理。处理完后自动退出事件表，并产生一个后继的事件。问题：投资可行性问题某个港口有一个万吨级的泊位，根据长期观察记录，依次到港的两艘船只的时间间隔有如下图的规律：到港船只时间间隔分布船只到港时间间隔/h151015202530频率0.150.100.120.140.170.260.06港口现有一台装卸机。根据其他港口的经验，若用两台装卸机可以节约装卸时间，两种情况下的装卸情况如下：船只装卸时间每条船只装卸时间/h一台装卸机两台装卸机频率14100.0516120.5018140.2020150.2022170.05船只装卸时，按照先到先装卸的原则进行。船只到港口，若泊位有空则立即停泊卸货；若泊位不空则排队等候。按照规定，到港的船只必须在15～30h内装卸完毕，其中包括等候和装卸的时间，若超过30小时，港口每小时支付200元的赔偿费；若能少于15小时,提前一小时奖励250元。港口在没有船只装卸时每小时的经济损失为400元，而每一艘船在港口停泊1小时损失200元。已知一台装卸机的购置和安装费用为60万元，折旧期为10年，每台装卸机每月的维修及油费的开支为3000元。请用计算机仿真的方法分析该港口添置第二台装卸机在经济上是否合算？购置设备的经济可行性以投资回收期来衡量，若其短于标准投资回收期，则添加设备是可行的；否则是不可行的。标准投资回收期为：T=k/c其中k＝60万元为添加设备的投资，c是一台装卸机和两台装卸机两种情况下的经营费用之差，即经营费用的节约值。经营费用包括:船只等待与卸货时间之和小于15小时的奖励费c1;船只等待与卸货时间之和多于30小时的赔偿费c2;船只停港损失费c3;港口空闲损失费c4;装卸机折旧费c5;维修与油料费c6;其中c5和c6两项是确定性费用,c1、c2、c3、c4这四项费用和船只到港间隔时间及卸货时间有关,因而是随机性的,可以由仿真来确定.船只的到港间隔时间的模拟抽样规则，即到港时间与均匀分布随机数的对应规则如下表所示：船只到港时间间隔/h频率随机数区间10.15(0.00,0.15）50.10[0.15，0.25）100.12[0.25，0.37）150.14[0.37，0.51）200.17[0.51，0.68）300.26[0.68，0.94）400.06[0.94，1.00）一台及两台装卸机装卸时间的模拟抽样的规则，即装卸时间于随机数的对应关系如下图：每条船只装卸时间/h一台装卸机两台装卸机频率随机数区间14100.05（0.00，0.05）16120.50[0.05,0.55）18140.20[0.55,0.75）20150.20[0.75,0.95）22170.05[0.95,1.00）我们可以产生两组随机数,一组用来模拟相邻两船到港的时间间隔,另一组随机数用来模拟各船的卸货时间,并累计船只停泊超过30小时和不足15小时的时间,累计船只在港时间和港口空闲时间.从而得到经营费用。仿真中先根据以上的两个表的随机数与到港时间和装卸时间的对应关系产生出如下表的到港时间间隔，到港时刻和装船的时间表到港时间和装卸时间装卸时间/h船号到港间隔到港时刻一台装卸机两台装卸机115151612230451612320651612457020155171161261081161274012120158301511612930181161210302111612113024116121210251161213302811612141529614101540336161216133716121715352221718403921612到港时间和装卸时间装卸时间/h船号到港间隔到港时刻一台装卸机两台装卸机19304221814203045216122115467161222547214102320492161224205121612251552716122640567221727155822015282060214102915617161230106271612312064720153220667181433156821612342070218143520722161236207421612根据上图的数据我们可以产生初始数据，利用事件步长法进行仿真。记A为船只到达事件，B为装卸结束事件对一台装卸机的情况前几步的仿真情况如下：1.产生初始事件表：初始事件表序号事件类型发生时刻1A152B452.处理1号事件A由上表知最早发生在第15小时的1号事件，置仿真时钟的时间为t=15将装卸机的状态由闲置为忙，产生结局事件B和下一个船只的到港时间A.计算港口的空闲损失费为15×400＝6000元。删除1号事件A，刷新事件表序号事件类型发生时刻2A403B304A65刷新后事件表3.处理3号事件B.又下表可以看出，最早的是3号事件B,为装卸的结束事件,它发生在第31小时，所以置t=31。判断装卸结束的事件是否小于15小时,或大于30小时,累加装卸机的工作时间。4.重复2.3每一步都找出最早的事件处理，判断是否到达结束仿真的时间，最后可以得到一台装卸机的各种费用。同理可以仿真的到两台装卸机的费用，其情况为：经营费用比较费用一台装卸机两台装卸机1c250190002c1200022003c1608001082004c552001164005c494698916c30116021总的经营费用235707223712可以看出，用两台装卸机时，每月可以节约的经营费用为11950元，所以投资回收期为T＝600000/(11950×12)＝4.17年所以，如果规定的投资回收期大于4.17年时，港口添加一台装卸机在经济上是可行的3.1排队系统的组成1.到达模式，顾客群体到达系统的情况。2.服务机构，系统对顾客群体中每一个成员服务时花去的时间长短和服务台个数。3.排队规则，系统提供服务的先后次序。三、排队论基础到达模式排队服务机构到达接受服务离去3.2排队系统的可解决的问题1.系统的最优设计（静态最优化）2.系统的运行控制（动态最优化）3.3排队系统的主要参数指标1.队长与等待队长。2.顾客在系统中的滞留时间与等待时间。3.系统的忙、闲比例。4.离去过程。5.顾客损失率。3.4到达模式表示平均到达时间间隔Ta＝T/n;平均到达速率λ＝1/Ta＝n/T;到达时间间隔分布函数A0(t)＝1-F(t)。到达时间变换系数：描述数据围绕平均值的分散程度。到达模式的形式：规则到达；完全随机到达；一般相同而独立的到达；成批到达；非平稳性到达；依态到达；连续到达。3.5服务机构平均服务时间：Ta;平均服务速率：μ;服务台形式：无服务台、单服务台、多服务台，无限服务台（IS），纵列系统，网络队列系统。3.6排队规则先到先服务（FIFO）；后到先服务（LCFS）；优先服务（PR）；共同服务（PS）；循环服务；占而不用。3.7排队系统模型的描述A/B/n/S/ZA:到达间隔时间分布；B:服务时间分布；n:服务台个数；S:排队室大小，排队队列的最大容量；Z:服务规则。3.8常用统计参数平均排队时间（※）；实体平均滞留时间；平均排队队长（※）；系统平均实体个数；设备利用率西安工业大学－计算机科学与工程学院三、单服务台排队系统的仿真举例个体序号到达时刻服务时间离开时刻等候时间12350265110376174493208513222761822447233271827－－－计算：平均排队时间；平均队长。状态变化点排队队长7192111132171182201220231240西安工业大学－计算机科学与工程学院例子：交换机系统的排队仿真N口交换机的N×N交换结构。西安工业大学－计算机科学与工程学院输入端口输出端口等待时间Hol阻塞时间交换时间服务时间逗留时间输入队列排头