食堂排队问题物流仿真项目计划书.doc

食堂排队问题物流仿真项目计划书
一、 仿真目的
应用仿真技术,对汀香一楼食堂排队问题的进行系统建模,通过仿真进行验证分析。考虑食堂购饭的窗口开设数目是否合适，以达到在高低峰期间能够合理配置资源，减少资源浪费，增加学生就餐满意度的目的。
二、 仿真问题描述
在汀香食堂一楼，经常看见这样的情况：食堂共4个打饭窗口，相当于4个服务窗口，在中午下午下课时间，食堂就餐学生特别多，往往每个窗口都是排着长长的队伍。
食堂的拥挤会造成排队，极大地增加了学生的时间成本，也会影响食堂的服务效率和服务质量。因此解决食堂排队问题，减少排队等待时间，是十分重要的。
然而对于食堂而言，也有更现实的问题，虽然增加窗口数量可减少排队等待时间，但同时也会增加食堂的运营成本，因此如何在两者之间权衡找到最佳的窗口数量，对学生和食堂双方来说是最合适和实用的。
食堂一般实行的是先来先服务原则，且学生可自由在队列间进行转移，并总向最短的队列转移，没有学生会因为队列过长而离去，故可认为排队方式是单一队列等待制。由于周末没课，学生去食堂就餐的时间比较分散，故只考虑周一到周五的情况。据本小组成员的观察，食堂就餐的学生一般都可找到座位就餐，因此食堂的容纳量是足够的，主要解决排队长与服务窗口的问题。
三、 仿真模型与步骤
1.食堂就餐排队系统模型假设
为了更好地研究就餐排队系统模型，本文对系统的组成要素进行假设：
(1) 排队规则：若食堂中有空闲的购饭窗口，则学生到达后可直接开始购饭， 如果有人正在接受服务，学生会选择队伍长度最短的窗口进行等候，直到窗口不再忙碌时再接受业务。
(2) 服务机构：假定食堂开放了c个购饭窗口，每个窗口都可以单独地为学生服务，互不干扰，一起工作，而且在同一时刻同一个窗口下一次只为一位学生服务。
2.食堂购饭排队系统性能指标
为了更好的研究排队系统特性，对得到的数据进行后续分析，需要考虑的系统性能指标有：
(1) 平均排队等待时间 Wqn=1ni=1nWi (2.1)
式中Wi —第 i 个旅客排队等待时间；
(2) 平均队长 Lqn=1Ti=1TmL(t) (2.2)
式中L(t)—t时刻排队等待的学生数目；
Tm—仿真时间上限。
在实际编程中，式(2.2)可采用以下形式：
Lqn=1mi=1mLi
式中m—队长的统计数目；
Li —第 i 次统计的实际队长。
3.食堂购饭排队系统的仿真
由于排队系统内发生的事件主要包括：学生到达事件，学生离去事件，因此采用离散事件对系统进行仿真。
（1）仿真流程分析
分别对排队系统的各个环节及总体模拟流程进行分析，画仿真流程图。
（2）系统内各状态变量的确定
学生到达食堂、服务完成后学生离去。
（3）系统仿真说明，为了清楚模型的参量，便于系统仿真，将该模拟程序中的子程序及主要变量进行定义。
（4）仿真结果分析，按照仿真流程图，利用 C 语言编程，对食堂购饭排队系统进行仿真。
四、拟采用的仿真方法
采用Anylogic软件编程实现仿真模型。多次重复仿真，分析求解方法的有效性，对排队等情况进行分析，求出最符合实际情况且达到系统平衡的方法，并利用软件编程实现不同问题模型的仿真。

目录

第一章 问题描述 1
1.1 食堂排队的影响因素 1
1.2 学生排队行为 2
第二章 数据的收集与整理 3
2.1 学生到达情况 3
2.2 窗口服务时间 4
第三章 模型构建 7
3.1 食堂排队系统的仿真模型 7
3.1.1 学生活动流程。 7
3.1.2 学生到达流程分析。 8
3.1.3 学生离去流程分析。 9
3.1.4 仿真的总体流程。 10
第四章 仿真软件程序 11
4.1 系统内各状态变量的确定 11
4.2 系统仿真说明 11
4.3 仿真条件假定 12
第五章 仿真结果分析 13
5.1 仿真模型图 13
5.2 仿真结果分析 13
5.3 优化调整 14
参考文献 15
成绩评定表 16

物流仿真软件及应用课程设计

第一章 问题描述

排队系统已经越来越广泛地被应用在生活中，随着考虑因素的复杂化，传统的数学方法已经很难解决排队问题，而计算机模拟可以有效、快速地分析系统特征与性能，模拟实际排队过程，可以很好地解决排队问题。因此，对排队模型进行仔细分析，并建立具有普遍性且可行的排队模拟，为理论上复杂又难解的排队问题提供了解决依据，进而可以拓展排队论的应用范围，对现实生活具有重大意义。
1.1 食堂排队的影响因素
对实际食堂排队系统整体运行过程与环节的观察与了解，可以发现，在食堂大厅内，设置有多个打饭窗口，每个运行的窗口前可以排一列队，学生因为就餐到达食堂，基本上是一个接着一个到来的，人数多时，需要排队，选择在哪个窗口打饭也是一个随机量，一般情况下会选择服务效率高或最短队伍进行排队，先来的学生先接受服务，当学生开始打饭时，服务人员一次只为一个学生提供服务，作业时间是不确定的。因此，从这整个排队的过程来看，学生随机到食堂就餐形成的排队系统是输入过程、服务机构的指标均为随机变化的排队系统。总结影响因素如下。
（1）学生到达情况分析
考虑食堂购饭窗口排队系统的特点，可以发现，进入食堂并且开始排队的学生是窗口服务的对象。学生到达数目较少、间隔较长，低于窗口开设数目及服务时间，则不需排队可直接接受服务，此时并无排队现象，学生到达数目较多、间隔较短，则需要进入队列排队，若达到的学生数目远远超过窗口的开设数目，会出现很严重的排队现象，此时会严重影响窗口的服务效率和学生的满意程度，造成学生的不耐烦情绪，因此，学生的到达情况，即排队系统的输入过程，是研究排队系统需要考虑的重要因素之一。 研究输入过程即是研究输入数据的概率分布特性，不难发现，在食堂，学生的到来是源源不断的，因此可以认为学生源是无限的。大多数学生达到食堂都是随机的，即前后两个学生是独立的，互不干扰的。
（2）服务效率分析
食堂一般开放多个窗口，同时提供服务，每个窗口都可以单独地为学生服务，互不干扰，而且在同一时刻同一个窗口下一次只为一位学生办理业务。学生按照来到食堂的先后顺序购饭，并且学生如果不是队首需要排队等待接受服务。在实际观察中，可以发现工作人员对每个学生的服务时间是不同的。
1.2 学生排队行为
学生在食堂就餐的过程中，或多或少的会受到来自外界环境因素、其自身等因素的影响。外界环境因素主要包括了食堂的结构设计、窗口在食堂内的位置及数量多少以及在窗口前现有顾客的排队情况；自身影响因素主要有：学生的性格、年龄、口味等。学生在食堂购饭时产生的行为主要有：对窗口的选择，是否选择排队，选择排队后接受来自工作人员为其提供的服务以及包括在排队过程中可能由于受到学生自身条件因素的影响或是受到其他排队队伍服务效率等的影响会选择离开本次排队队列等，下面就学生选择在食堂购饭时产生的行为做具体的分析。
(1) 窗口选择。学生到食堂后会按照现在食堂开放窗口的数量、每个窗口前排队等候的学生数量选择在哪个窗口排队等待。学生通常会进入窗口前没有人或者人数少、服务效率高的窗口购饭或者排队等待，这样学生就可以减少等待时间，节约出行时间成本。在本文中学生到达食堂时即根据队列长度选择较短队列进行排队。
(2) 跟随。是指学生进入食堂并且选择好队列后，不再变动，一直跟着前方排队人员，不会因为外界环境，如其他窗口是否行进速度快、服务效率高、队伍短而选择离开当前队伍，会一直排队直到办理业务的行为。
综上所述，根据食堂窗口的特征以及排队论的知识，发现食堂排队系统属于一个先到先服务，等待制的多服务台排队系统。

第二章 数据的收集与整理

2.1 学生到达情况
通过统计调查，得到汀香一楼食堂11：:30-12:00高峰期 2016年 5 月24日-26 日 3 天的学生到达人数的数据，对其进行整理和分析，得到各时段的样本均值，如表2-1所示。
表 2-1 学生到达人数汇总表
时间段
每1分钟学生到达数
11：30-11：40
2
8
9
7
10
9
7
9
8
13
11:41-11：50
14
12
15
15
14
16
15
20
20
27
11：51-12：00
34
28
33
29
17
14
10
9
10
7
利用x2检验法对每个时段内的学生到达分布是否为泊松分布进行检验，具体过程如下：根据每个时间段内的学生到达数据，对单位时间内即每1分钟内到达的学生数是否服从泊松分布进行假设检验。随机调查该时间段的30个单位时间内到达的旅客情况整理如图2-1所示：

图2-1 单位时间学生到达数直方图
提出假设H0：该时段内的学生到达数服从参数为 l 的泊松分布。 计算 l 的最大似然估计值为：
l=1ni=1nXi=X
由表2-1中的数据可知该时段的单位时间内的旅客到达率X=43.25，所以l =43.25。根据相关理论知识，将 X 的取值分为五组，记为A1=0，5，A2=6，10 ，A3=11，15，A4=16，20， A5=[21，∞) 。计算 X = i 的概率Pi ， Pi的估计是:
Pi=i=aj-1ajlii!e-l
计算过程如下：
表 2-2 学生到达分布检验过程
Ai
fi
Pi
nPi
　fi-nPi2nPi
A2
12
0.3524
10.572
0.193
A3
8
0.3976
11.928
1.294
A1
1
0.0135
15.688
0.405
0.874
A4
10
5
0.1198
11.98
0.550
A5
4
0.1101
3.303
0.147

30

3.058
由表 2-2 的计算结果可知x2=3.058，将理论频数小于5的组合并，此时组数为 k = 3，r =1 ，故x2的自由度为 k - r -1= 3 -1-1= 1，取显著性水平 a = 0.05 ，查表可知x0.0521=3.841，则有x2<x0.0521，故可以认为该时段的单位时间内的学生到达数服从 l =14.7的泊松分布。因此，本文以每个时间段为单位对开设窗口数进行优化。
2.2 窗口服务时间
（1）随机调查的60位学生在窗口4的购饭时间即食堂窗口服务时间，如表2-3：

表2-3 服务时间 (单位：秒)
64
72
80
106
67
136
130
102
77
63
76
150
138
76
63
72
93
155
114
89
125
61
104
147
99
82
79
128
132
104
144
83
94
113
143
117
91
86
115
85
169
150
122
134
84
61
74
69
83
73
62
159
72
85
103
116
68
97
108
64
其中，
将上表中数据分组，得频数表如下：
表2-4 窗口服务时间频数
对学生服务时间及其频数
服务时间（s）
频数
0-59.5
0
59.5-69.5
10
69.5-79.5
9
79.5-89.5
9
89.5-99.5
5
99.5-109.5
6
109.5-119.5
5
119.5-129.5
3
129.5-139．5
5
139.5-149.5
3
149.5以上
5
合计
60
假设对学生的服务时间符合负指数分布。
负指数分布公式模型如下：
pt=λ2e-λ2t-μ2，t≥μ20 ,t<μ2;λ2>0,μ2≠0
由于负指数分布有两个未知参数，因此我们利用样本对其进行极大似然估计。μ2=t1=61(s),1λ2=39.13,λ2=0.0256,将T分为5组，记A1=0，59.5，A2=59.5，89.5 ，A3=89.5，119.5,A4=119.5，149.5，A5=[149.5，∞)。
表 2-5 服务时间分布检验过程
Ai
fi
Pi
nPi
　fi-nPi2nPi
A2
0
0
0
0
A3
28
0.557
33.42
0.879
A1
16
0.2584
15.504
0.016
A4
11
0.1199
7.194
2.014
A5
5
0.0647
3.882
0.322

60
1

3.231
检验统计量：x2=3.231，据x2分布表，由于pt中含有两个个未知参数λ、μ，故r=2。将A2，A5合并，则k=4.取检验水平，故x2的自由度为 k - r -1=1.查x2表得临界值x0.0521=3.841，x2<x0.0521，故应肯定原假设H0，即认为对顾客的服务时间服从负指数分布。
（2）随机调查的60位学生在窗口1-3的购饭时间即食堂窗口服务时间，如表2-6：
表2-6 服务时间 (单位：秒)
12.00
9.60
12.80
14.70
9.40
18.30
14.70
9.60
9.80
14.70
8.70
9.90
9.00
10.70
13.00
9.34
14.40
14.52
9.10
14.20
13.40
10.70
14.60
14.30
9.00
11.10
14.00
15.30
16.00
9.52
10.00
14.40
20.30
14.00
9.50
11.40
19.60
11.70
14.10
20.30
17.20
12.70
9.40
12.00
14.10
15.40
14.10
20.70
9.67
9.50
14.30
16.40
11.10
21.70
10.40
14.10
14.70
15.40
18.10
22.00

图 2-1 窗口服务时间曲线图
根据数据整理，窗口1-3的服务时间相近连续。我们可以把它看成是低限为 a=9、众数为 c=14、上限为 b =20的三角形分布。

第三章 模型构建

3.1 食堂排队系统的仿真模型
由于排队系统内发生的事件主要包括：学生到达事件，学生离去事件，窗口选择效率事件，因此采用离散事件仿真策略中的事件调度法对系统进行仿真。
仿真流程分析仿真系统中，需要设置的模块有：系统初始化，包括仿真时钟，队列，窗口，事件表等数据的初始化，学生到达模块，学生排队模块，队伍变化状态模块，窗口状态模块，学生就餐完毕离开模块等。
3.1.1 学生活动流程。
分析一名学生为了就餐而来到食堂，选择窗口前的队列排队，开始购饭，直到就餐完成后离开食堂的一系列过程。学生因为就餐来到食堂，此时系统将会按照当前所有开放窗口的排队 情况将这名学生进行分派，放至排队长度最小的窗口，若此刻该窗口前没有人排队，也没有人正在购饭，即窗口空闲，那么学生直接接受服务，反之，需要排队等候。在排队过程中，队伍第一个学生完成服务后离开，窗口空闲，队伍前移，下一位学生开始服务，此时，窗口繁忙，依次循环，直到该学生成为队伍第一个，而窗口空闲，才开始接受服务。过程见图3-1。

排队队列空？
学生到达
选择队列排队

N

Y

Y
N
窗口空

排队等待
接受服务，窗口忙

有学生离去
学生离去，窗口闲

位置前移
结束

图3-1 学生活动流程图

3.1.2 学生到达流程分析。
当需要就餐的学生来到食堂，仿真程序会使用相应的统计函数自动计算出这时食堂各个窗口前队伍的长度，并且会自动存储学生到达后开始排队的时间，接着使用相应函数对这名学生进行分派，使学生进入最短队列排队，同时判断是否队伍为0且窗口空闲，若是，则立刻开始接受服务，并记录下时刻，反之，学生在队伍中排队等候服务，此时，队列长度增加一人，当队伍的第一人离开系统时，学生所处位置前移，直到成为队伍第一个，当窗口空闲时，记录此时学生接受服务的时间，并计算系统中该学生的等待时间。

排队队列空？
学生到达
记录数据
排队等待

N

Y

窗口空？
N
窗口空？

Y

结束
接受服务

图3-2 学生到达的程序操作流程

3.1.3 学生离去流程分析。
当学生接受完服务后便离开系统，此时队长减少一人，记录相关数据，即离开时间和此刻系统中队伍长度，窗口状态置为空闲。
学生离去
记录数据，统计队长，队长减一一

置窗口为空闲

服务结束

图3-3 学生离去的程序操作流程

3.1.4 仿真的总体流程。
在仿真程序中，包括系统的初始化模块、学生到达模块、学生排队模块、学生离去模块等。学生到达及服务时间的分布情况是已知的，根据各自的概率分布编写相应程序模拟学生到达及服务时间产生的情形，直到达到仿真时间结束，结束仿真，输出系统各项指标数据，为后续的优化做准备，具体流程如图 3-5所示。
Y

系统初始化
推动仿真时钟化
学生到达处理化
所有队列为空？化
服务台为空？
接受服务，置服务台为忙
等待服务，队长加一空？化
学生离去，队长减一空？化
仿真时间到？
输出各项系统值空？化
结束
队首空？化
N
N
Y
N

图3-4 仿真流程图

第四章 仿真软件程序

4.1 系统内各状态变量的确定
（1）定义事件，引起主动成分状态变化的事件有学生到达食堂、服务完成后学生离去，因此本模型应定义如下三类事件：
表 4-1 事件表及其定义
成分
事件定义
事件类型
学生
学生到达食堂
1
队列
学生排队，中途可能换队
2
学生
服务完成后离开
3
（2）建立模拟记录表，便于收集、统计数据。表中 n 个队列对应有 n 个记录相关属性，如到达时间，所属窗口等。该模型中设置如表4-2所示：
表 4-2 队列中学生的相关属性
表的编号及其定义
记录数
属性1
属性2
属性3
……
1
队列
n
到达时间1
到达时间2
到达时间3
……
2
窗口
n
忙闲状态
开始空闲的时间
忙闲状态
……
3
事件表
不定
发生时间
事件类型
窗口编号
……
4.2 系统仿真说明
(1) exit1程序段（排队程序）：
int N=0,m,i;
int a [];
a=new int [4];
a [0] =queue.size ();a [1] =queue1.size ();
a [2] =queue2.size ();a [3] =queue3.size ();
m=min(min(a [0],a [1]),min(a [2],a [3]));for(i=0;i<4;i++) if(a [i] ==m)N=i;
switch(N)
{
case 0:enter.take(entity);break;
case 1:enter1.take(entity);break;
case 2:enter2.take(entity);break;
case 3:enter3.take(entity);break;
}
程序说明 ：将0-3队（对于窗口1-4）的当前队长给数据组a[]，找出最短队，将下标赋值给N，并让当前实体进入enterN。该程序在有实体离开exit1时触发。
（2）仿真模型的参量说明。
表4-3 参量说明
模型变量名
定义
模型变量名
定义
a[]
数组
queue.size()
队列排队长度函数
TotalServiceTime
总服务时间
ServiceRate
窗口利用率
exponential()
负指数分布函数
double EnterSystemTime
双精度变量,用于记录实体进入系统的时间
double[]
BeginServiceTime
双精度数组,用于记录实体在各窗口的开始服务时间
TimeInSystem.add(time()
-entity.EnterSystemTime)
表示把当前时间减去实体进入系统的时间送人TimeInSystem统计表
4.3 仿真条件假定
假定某个食堂窗口数为 4，仿真时间为 60min。设学生到达时间间隔是服从参数l = 15的泊松分布的随机变量，每个窗口前排有一个队列，队列数与服务台数相等，每队最大队长均为 45，新到达的学生加入到最短队列中，排队之后不换队。

第五章 仿真结果分析

5.1 仿真模型图
按照仿真流程图，利用 C 语言编程，对食堂排队系统进行仿真。利用Anylogic建立的仿真模型如图5-1。

图5-1 仿真模型图
5.2 仿真结果分析
（1）首先分析学生在系统的时间。其最短时间为0.16min，最长为15.18min，平均3.61min，对学生来说是属于可接受的范围，但是其中窗口4 因为菜品不同，因此服务时间较长。
（2）仿真结束时，各窗口排队长分别为26、26、27、28,队长均值分别为12、15、15、14。由于学生选择最短队排队，所以各队长较相近。但是数值却不小，有排长队的现象。因为学生窗口的服务不是特别满意。
（3）窗口使用率分析，各窗口使用率分别是97.4%、98.7%、99.4%、99.9%。4个窗口的使用率都非常高，没有造成资源浪费。
综合以上分析，我们为仿真进行优化，增加窗口。
5.3 优化调整
因为仿真各窗口都存在排队现象，因此增加了一个与窗口4相同的窗口（系统限制只能增加一个进行优化）。

图 5-2 优化后仿真模型图
（1）首先学生在系统的时间，最短时间为0.16min，最长为6.72min，平均1.11min，对学生来说是各窗口服务时间都减少，窗口服务效率也明显提高。
（2）仿真结束时，各窗口排队长分别为0、0、1、1、1，队长均值1、1、2、2、2。数值都较小，没有了排长队的现象。因此增加窗口后会让学生增加对窗口服务的满意度。
（3）窗口使用率分析，各窗口使用率分别是80.3%、92.5%、98%、99.9%。5个窗口的使用率都非常高，没有造成资源浪费。
综合以上分析，我们进行的仿真优化结果是增加1个窗口，对食堂和学生都是最适合的和实用。

参考文献

[1] 李维铮,李梅生,胡运权等.运筹学[M].北京:清华大学出版社,2000.
[2] 宋卫斌,苏秦.虚拟顾客服务系统排队模型[J].管理科学学报,2001,4(3):53-64.
[3] 郑欢,古福文.大型超市顾客交费排队系统优化分析[J].管理学报.2005,2(2):171-173.
[4]纪莹,徐行方.基于排队论的售票厅售票组织优化[J].运输经济.2008(12),134-136.
[5]薛春荣.超市排队系统不同时段收银台数目优化研究[J].河南科学.2013,31(5):608-610.
[6] Levy Y,Yechiali U,Utilization of idle time in an M/G/1 queuing system[J], Management Science, 1996(2),202-211.
[7] 张坚.医院门诊排队的模拟和优化建设[J].系统工程理论与实践.1998(3).140-144.
[8] 李深.基于计算机仿真的排队系统优化问题研究[D].潘阳工业大学.2007.
[9] 张博,王亚超.基于 Witness 的超市收款服务系统的仿真与优化[J].商场现代化.2007(122), 39-40.
[10] 杨树国,杨晓研.基于排队论的优化银行服务方面的研究应用[J].科技视界,2015(2):89-90.

目 录
第一章 总论 1
1.1项目名称与承办单位 1
1.2研究工作的依据、内容及范围 1
1.3编制原则 3
1.4项目概况 3
1.5技术经济指标 5
1.6结论 6
第二章 项目背景及建设必要性 8
2.1项目背景 8
2.2建设的必要性 9
第三章 建设条件 11
3.1项目区概况 11
3.2建设地点选择 错误！未定义书签。
3.3项目建设条件优劣势分析 错误！未定义书签。
第四章 市场分析与销售方案 13
4.1市场分析 13
4.2营销策略、方案、模式 14
第五章 建设方案 15
5.1建设规模和产品方案 15
5.2建设规划和布局 15
5.3运输 18
5.4建设标准 18
5.5公用工程 20
5.6工艺技术方案 21
5.7设备方案 21
5.8节能减排措施 24
第六章 环境影响评价 25
6.1环境影响 25
6.2环境保护与治理措施 26
6.3评价与审批 28
第七章 项目组织与管理 29
7.1组织机构与职能划分 29
7.2劳动定员 29
7.3经营管理措施 30
7.4技术培训 30
第八章 劳动、安全、卫生与消防 31
8.1编制依据及采用的标准 31
8.2安全卫生防护原则 31
8.3自然灾害危害因素分析及防范措施 32
8.4生产过程中产生的危害因素分析及防范措施 32
8.5消防编制依据及采用的标准 34
8.6消防设计原则 35
8.7火灾隐患分析 35
8.8总平面消防设计 35
8.9消防给水设计 36
8.10建筑防火 36
8.11火灾检测报警系统 37
8.12预期效果 37
第九章 项目实施进度 38
9.1实施进度计划 38
9.2项目实施建议 38
第十章 项目招投标方案 40
10.1招标原则 40
10.2项目招标范围 40
10.3投标、开标、评标和中标程序 40
10.4评标委员会的人员组成和资格要求 42
第十一章 投资估算和资金筹措 43
11.1投资估算 43
11.2资金筹措及使用计划 45
第十二章 财务评价 47
12.1费用与效益估算 47
12.2财务分析 48
12.3不确定性分析 49
12.5财务评价结论 50
第十三章 建设合理性分析 51
13.1产业政策符合性分析 51
13.2清洁生产符合性分析 51
13.3规划符合性分析 51
13.4项目建设环保政策符合性分析 51
13.5环境承载性分析 51
13.6结论 52
第十四章 结论与建议 53

第19页共16页