交通仿真讲义第三章离散仿真讲义

1、第三章 离散系统仿真概述离散系统的状态只是在离散时间点上发生变化，而且这些离散时间点一般是随机的。离散系统的数学模型通常用流程图或网络图来描述。离散仿真的目的是实体的活动以找出(分析)系统的(潜在)行为。离散仿真事件(event)时间因变量随(事件)时间离散变化。因变量时间基本概念实体 构成系统的各种成分称为实体，系统的研究对象。 临时实体、永久实体。属性 反映实体的某些性质。状态 在某一确定时间点，系统的状态是系统中所有实体属性的集合。基本概念事件 引起系统状态发生变化的行为，它是在某一时间点上的瞬间行为。活动 实体在某一状态的持续过程。进程 进程由和实体相关的事件及若干活动组成，一个进程描

2、述了它所包括的事件及活动间的相互逻辑关系和时序关系。车辆到达事件服务开始事件服务结束事件服务活动进程排队活动事件、活动、进程三者之间的关系基本概念仿真时钟 仿真时钟用于表示仿真时间的变化，仿真时钟推进的时间间隔称为时间步长。 时间步长法、事件步长法统计计数器 离散系统的状态随事件的不断发生呈现动态变化过程，这种动态变化过程在统计意义下才有参考价值。 统计计数器用于记录仿真规程中系统性能的统计信息。离散仿真模型建立步骤定义系统的参变量集合，构造系统映象；定义事件类型及其发生时点；定义每一事件时间发生状态变化的点；描述系统中实体的活动；构造状态转移函数或算法；通过系统流图描述整个过程。构造初始映象

3、找一个尽可能简单的系统状态作为初始状态；从一个远离平衡状态但容易构造的状态开始进行模拟，当运行一段时间之后，系统的状态会接近或处于稳定状态（平衡状态）。离散仿真中的关键问题事件取舍： 确定改变系统状态的事件集，并将它们用逻辑关系联系起来。 系统仿真就是靠按顺序执行联系这些事件的逻辑来实现的。 活动扫描： 要描述系统中实体的活动，设定引起活动开始及结束的条件。启动或结束活动的事件不由建模者设定，而随仿真推进而自动进行。为保证活动得到记录，每一次推进均需要扫描活动中的实体集。 由于需要在每一次推进中扫描每一活动，这种方法的效率不太高。过程定位： 提供整个仿真过程中实体流的一种用仿真语言描述的方法。

4、系统仿真的推进时间步长法 以固定的时间间隔进行驱动；事件步长法 按下一类最早发生事件的发生时间推进。t时间步长法在进行系统仿真的同时，把整个仿真过程分为许多相等的时间间隔，程序按此步长前进的时钟就是仿真时钟。在每个时间间隔做如下处理：该步内若无事件发生，则仿真时钟再推进一个单位时间；若在该步内有若干个事件发生，则认为这些事件均发生在这一步的结束时刻，同时必须规定当出现这种情况时各类事件处理的优先顺序。扫描与处理方法对每一类事件或每一个主导实体设置一个模拟时钟，以此记录和控制实体活动的延续时间。对系统实体进行扫描；对系统事件进行扫描；对事件和实体结合起来进行扫描。改变状态，预测下一事件。时间步长

5、法流程图事件步长法仿真时钟不断地从一个事件发生时间推进到下一个最早发生事件的发生时间。以事件发生的时间点相互间隔作为步长，按照时间的进展，一步一步地对系统的行为进行仿真，直到预定的仿真时间为止。事件表按照事件发生时间先后顺序安排事件，将仿真过程看作一个事件点序列。事件控制部件始终从事件表中选择最早发生时间的事件记录，然后将仿真时钟该事件发生的时刻。时钟推进举例模拟运行150个时间单位，顾客到达事件、顾客服务完毕离去事件，Ti到达时间间隔，Si第i个顾客服务时间，Di第i个顾客等待时间，Ci= Ti+Si+Di第i个顾客离开系统的时间, qi第i个顾客排队的对长，bi模拟时钟推进到第i次的时间，

6、Zi第i个事件发生时服务员的状态。 Ti=15,32,24,40,22,； Si=43,36,34,28, 初始状态： q0=0, Z0=0求： bi 、Ci、离散系统仿真程序的结构状态变量；时钟变量；事件表（按时间顺序记录仿真过程中将要发生的事件）统计计数器初始化子程序时钟推进子程序 （由事件表确定下一事件，然后将将仿真时钟推进到该事件发生的时间）调度子程序（将仿真过程中产生的未来事件插入事件表）离散系统仿真程序的结构事件子程序 每一类事件对应一个事件子程序，相应的事件发生时就转入该事件子程序进行处理，更新系统状态，产生新的事件。统计报告子程序随机数发生器主程序 调用时钟推进子程序，控制转移

7、到相应的事件子程序，完成仿真程序的总体控制。离散系统仿真程序流程图离散系统仿真策略建立描述系统行为的仿真模型 由于系统采用伪随机数，可以得到确定的状态转移函数，模型采用流程图或网络图的形式。仿真策略 仿真策略决定仿真模型的结构模型描述中采用的主要术语成分 相当于系统中的实体，用于构造模型中的各个部分。主动成分 （可以主动产生活动的成分） 被动成分 （本身不激发活动，只有在主动成分作用下才能产生状态变化）描述变量 成分状态、属性的描述。成分间的相互关系 描述成分之间相互影响的规律。模型描述中采用的符号C=a1,a2,an为成分集合； CA=a1,a2,am为主动成分子集合；CP=a1,a2,ah

8、为被动成分子集合；Sa为成分a的状态变量；P=p1,p2,pr为参数（属性）集合；ta为成分a的状态下一发生变化的时刻;Da(S)为成分a在状态变量值S时的条件变量;TIME为模拟时钟的值。典型仿真策略事件调度法活动描述法进程交互法事件调度法通过定义事件及每个事件发生对系统状态的变化，按时间顺序确定并执行每个事件发生时有关的逻辑关系。所有事件均放在事件表中，模型中设有一个时间控制机构，该机构从事件表中选取最早发生时刻的事件。以事件种类为控制依据，不同种类事件的处理进入相应的事件处理模块，并在时间处理完毕返回时间控制机构。事件调度法模型的基本结构事件调度算法初始时间t= t0、事件表初始化、置系

9、统初始事件;成分表初始化S =(Sa1,ta1), (Sam,tam), Sam+1, San);操作事件表，取出t =min ta|aCA ，修改事件表;推进时钟 TIME= t(s);While (TIMETIMEFRESENT(S)=a|ta=TIMEPAST(S)=a|taTIMETIME=min(ta|a FUTURE(S)活动扫描法的算法初始时间t= t0、设置主动成分的模拟时钟ta(i);成分状态初始化S =(Sa1,ta1), (Sam,tam), Sam+1, San);设置系统时钟 TIME= t0;While (TIME= t)，执行扫描 for j=最高优先数到最低优先

10、数 将优先数为j的成分置成i if （ tai(i)= TIME 且 Dai(S)= true） 执行活动子例程 endif endfor TIME=min(ta|aFUTURE(S) endwhile按活动扫描法建立的排对系统模型进程交互法进程由事件的时间序列及若干活动组成具有上述两种方法的特点，接近实际系统，编程实现非常复杂采用进程描述系统，将模型的主动成分所发生的事件及活动按照时间顺序进行组合形成进程表，一个成分一旦进入进程，它将完成进程的全部活动。进程交互法采用两张事件表，当前事件表、将来事件表；当仿真时钟推进，满足条件的所有事件记录从将来事件表移到当前事件表，取出每个事件记录，判断所

11、属进程与位置，当发生条件真，发生包含该事件的活动，并让该进程尽可能地推进，直至结束。时间控制以主导实体进入该进程的的时间序列及其经历该进程的各项活动的时间顺序，走向控制主要以断点为依据。以进程为基础的排队系统模型几种仿真策略的比较系统描述事件调度法中，只有主动成分才能施加作用；事件调度法中，系统的动态特性表现为主动成分不断产生事件；活动扫描法中表现为主动成分产生活动；进程交互法中则是通过成分在其进程中一步一步地推进描述。 几种仿真策略的比较建模要点事件调度法中，用户要对所定义的全部事件进行建 模，条件测试只能在事件处理子程序中进行； 活动扫描法设置了一个条件子例程用于条件测试，还设置了一个活动

12、扫描模块，该模块对所定义的活动进行建模； 进程交互法将一个进程分成若干步，每一步包括条件测试及执行活动两部分。 几种仿真策略的比较时钟推进事件调度法中，控制机构从事件表中取出最早发生时间的事件记录，将时钟推进到该时刻，执行该事件处理子程序； 活动扫描法除系统时钟外，每一个主动成分还有成分模拟子时钟，控制机构选取那些大于当前系统时钟且所有成分模拟时钟最小的那个成分模拟时钟，将系统时钟推进到该时钟； 进程交互法中，一旦某个进程被执行，要求尽可能走下去，但并不改变系统时钟。如果该进程未完成，记录中断时间及事件类型放入将来事件表。几种仿真策略的比较执行控制事件调度法按下一最早发生时间选择事件记录； 活

13、动扫描法对全部活动扫描，只有Dai(S)=true且taiTIME的活动才能被执行； 进程交互法对当前事件表中所有的记录扫描，根据该事件在其进程中的指针进行条件判断。当Dai(S)=true 执行该进程，并一直执行下去，否则记下断点。 几种仿真策略的比较事件调度法建模灵活，建模工作量大；活动扫描法对于各成分相关性很强的系统来说模型效率较高，但执行程序结构复杂；进程交互法建模最直观，模型表示接近实际系统特别适用于可以预测、顺序比较确定的系统，但流程控制复杂，建模灵活性不好。适用性系统中的各个成分相关性较少，宜采用事件调度法，反之宜采用活动扫描法；系统成分的活动比较规则，宜采用进程交互法。排队系统

14、仿真某个时刻要求服务的的数量超过服务机构的容量。到达模式：顾客按怎样的规律到达；服务机构：同一时刻有多少服务台可接纳顾客；排队规则：服务台完成当前的服务后，从对列选择下一个实体服务的原则。到达模式平均到达间隔时间Ta=T/n；平均到达率 =1/Ta ；到达间隔时间分布函数；1-F(t)到达时间变化系数Sa/Ta ；服务机构一个或多个服务员，没有服务员；多个服务台，并列、串列、混合；单独服务，成批服务；服务时间，确定、随机；排队规则先到先服务FIFO；后到先服务LIFO；随机服务SIRO；优先权服务PR；最短处理时间先服务SPT。系统容量有限；无限。排队模型分类A/B/C/D/E A:到达模式；

15、B:服务模式；C:并行服务员的数目； D:系统容量；E:排队规则。常见的到达和服务间隔时间分布：M(指数分布)、 G(确定型分布)、 M(一般随机分布)；M/M/1/FIFO排队系统的特征量服务员利用率P=平均服务时间/平均到达间隔时间；系统中平均顾客数P/1- P；系统内排队等待的顾客数P2/1- P ；顾客在系统内的停留时间；平均等待时间；系统出现大于n个顾客的概率。研讨多队多服务台排队系统模拟 以多出纳台银行系统为例（换对原则）系统中的成分成分状态变量系统中的初始状态顾客到达间隔时间服从的分布 服务时间服从的分布 事件顾客到达银行顾客完成服务后离开顾客换对银行关闭系统指标平均等待时间平均对长最大等待时间最大对长多级串联封闭式排队系统模拟基于主导实体活动扫描法的模拟模型从一种称为主导实