Matlab-离散事件系统仿真实验

分析一个单服务台排队系统

1、排队系统的一般理论

一般的排队系统都有三个基本组成部分：

（1） 到达模式：指动态实体（顾客）按怎样的规律到达，描写实体到达的统计特性。通常假定顾客总体是无限的。

（2） 服务机构：指同一时刻有多少服务设备可以接纳动态实体，它们的服务需要多少时间。它也具有一定的分布特性。通常，假定系统的容量（包括正在服务的人数加上在等待线等待的人数）是无限的。（3） 排队规则：指对下一个实体服务的选择原则。通用的排队规则包括先进先出（FIFO），后进先出（LIFO），随机服务（SIRO）等。

2、离散系统常用的仿真策略

（1）事件调度法（Event Scheduling）：

基本思想：离散事件系统中最基本的概念是事件，事件发生引起系统状态的变化，用事件的观点来分析真实系统。通过定义事件或每个事件发生系统状态的变化，按时间顺序确定并执行每个事件发生时有关逻辑关系。

（2）活动扫描法：

基本思想：系统有成分组成，而成分又包含活动。活动的发生必须满足某些条件，且每一个主动成分均有一个相应的活动例程。仿真过程中，活动的发生时间也作为条件之一，而且较之其他条件具有更高的优先权。

（3）进程交互法：

基本思想：将模型中的主动成分历经系统所发生的事件及活动，按时间发生的顺序进行组合，从而形成进程表。系统仿真钟的推进采用两张进程表，一是当前事件表，二是将来事件表。

3、本实验采用单服务台模型

（1） 到达模式：顾客源是无限的，顾客单个到达，相互独立，一定时间的到达数服从指数分布。

（2） 排队规则：单队，且对队列长度没有限制，先到先服务的 FIFO规则。

（3） 服务机构：单服务台，各顾客的服务时间相互独立，服从相同的指数分布。

（4） 到达时间间隔和服务时间是相互独立的。

4、仿真运行方式

仿真运行方式可分为两大类：

（1）终止型仿真：仿真的运行长度是事先确定的由于仿真运行时间长度有限，系统的性能与运行长度有关，系统的初始状态对系统性能的影响是不能忽略的。为了消除由于初始状态对系统性能估计造成的影响，需要多次独立运行仿真模型。

（2）稳态型仿真：这类仿真研究仅运行一次，但运行长度却是足够长，仿真的目的是估计系统的稳态性能。

5、系统设计

采用事件调度法，事件调度法共有三种事件：即“顾客到达”、“服务开始”、“服务结束”。其中“服务开始”为条件事件，其条件是“顾客队列长度不为零且服务员空闲”。在用事件调度法时，不单独考虑条件事件，而将其并入非条件事件中。因此，需要考虑的事件例程有“顾客到达时间例程”和 “服务结束时间例程”。“顾客到达时间例程”和 “服务结束时间例程”如图

事件调度法的程序结构设计如下：

（1）初始化。

 设置仿真的开始时间 t0和结束时间 tf ；

 设置实体的初始化状态；

 设置初始事件及其发生时间 ts。

（2）仿真时钟 TIME = ts。

（3）确定在当前时钟 TIME 下发生的事件类型 E（i = 1,2,3，…，n），并按解结规则排序。

（4）如果 TIME <= tf执行。

{case E1：执行 E1 的事件例程；产生后续事件类型及发生时间；……case En：执行 En 的事件例程；产生后续事件类型及发生时间；

}

否则，转(6).

（5）将仿真时钟 TIME 推进到下一最早事件发生时刻；转(3)。（6）结束仿真。

6、思路分析

采用事件调度法来研究单服务台排队系统。顾客逐个到达服务台，且相邻两个顾客到达服务台的时间间隔服从参数为 3 min 的指数分布。到达服务台后，若这时服务员空闲，则为其提供服务，若此时服务员正在为其他顾客服务，则刚到的顾客排队等待。服务员为每位顾客服务的时间长度服从参数为 4 min 的指数分布。使用 Matlab 软件进行建模仿真，用 exprnd 函数生成符合指数分布的随机数。用三个空白数组分别存储第 i 个顾客引起的三种事件先后发生的时刻，对获得的参数按照时间顺序进行整理和分析，可以得出平均队长、平均等待时间等重要参数，流程样例如下

单服务台排队系统的流程框架图如下：

clc
close all
clear
rng default
T = 1000; % T-仿真长度（min）
mu1=3; %顾客到达时间间隔（指数分布）的均值
mu2=4; %服务时间（指数分布）均值

arriveGap = []; % 到达时间间隔
serveGap = [];  % 服务时间
Arrive=[]; %顾客到达时间
Serve=[]; %服务开始时间
Leave=[]; %服务结束时间
%===========初始化============%
i=0;  % 第几个顾客抵达
t0=0;  % 仿真开始时间
TIME=0; % 系统时间
arriveGap = [arriveGap, exprnd(mu1)];
Tarrv=t0+arriveGap(i+1); % 第i个抵达时间
serveGap = [serveGap, exprnd(mu2)];
Tleave=Tarrv+serveGap(i+1);  % 第i个服务结束时间
Arrive=[Arrive,Tarrv];
Serve=[Serve,Tarrv]; 
Leave=[Leave,Tleave];
% 开始仿真
while TIME < T
    i=i+1;
    arriveGap = [arriveGap, exprnd(mu1)];
    Tarrv=Tarrv+arriveGap(i+1); %确定下一顾客到达时刻
    % 留个空白思考 
    if Leave(i)<=Arrive(i+1) %服务员空闲，无需排队
        % 留个空白思考 
    else %服务员忙碌，需要排队
        % 留个空白思考 
    end
    % 留个空白思考 
    Leave=[Leave,Tleave];
    
end
%%

Twait=[]; %每个顾客的等待时间
Nwait=[]; %每个顾客接受服务时的队长
WaitNum=0;  % 平均等待人数
WaitTime=0;  % 平均等待时间
%平均等待时间
for j=1:i
    Twait(j)=Serve(j)-Arrive(j);
    WaitTime=WaitTime+Twait(j);
end
WaitTime=WaitTime/i;
%平均队长
for m=1:i
    k=0;
    for n=m+1:i
        if Leave(m)>Arrive(n)
            k=k+1;
        else
            break;
        end
    end
    Nwait(m)=k;
    WaitNum=WaitNum+Nwait(m); 
end 
WaitNum=WaitNum/i;
%绘图
figure
plot(arriveGap)
xlabel 顾客序号
ylabel 到达间隔时间/min
figure
plot(serveGap)
xlabel 顾客序号
ylabel 服务时间/min
figure
plot(Twait);
xlabel 顾客序号
ylabel 等待时间/min
figure
plot(Nwait);
xlabel 顾客序号
ylabel 队伍长度/人

disp(['平均到达间隔时间：', num2str(mean(arriveGap)),' min'])
disp(['平均服务时间：', num2str(mean(serveGap)),' min'])
disp(['平均等待时间：', num2str(WaitTime),' min'])
disp(['平均队长指标：', num2str(WaitNum),' 人'])