早高峰电梯系统分析与优化

承 诺 书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 宁 波 工 程 学 院 参赛队员 (打印并签名) ：1. 顾豪 2. 郑启奔 3. 施雪丹 指导教师或指导教师组负责人 (打印并签名)： 数 模 组

日期： 2010 年 9 月 14 日

赛区评阅编号（由赛区组委会评阅前进行编号）：

2010高教社杯全国大学生数学建模竞赛

编 号 专 用 页

评 阅 人 评 分 备 注 赛区评阅编号（由赛区组委会评阅前进行编号）：

赛区评阅记录（可供赛区评阅时使用）：

全国统一编号（由赛区组委会送交全国前编号）：

全国评阅编号（由全国组委会评阅前进行编号）：

早高峰电梯系统分析与优化

摘 要

电梯系统直接关系到乘客与楼房管理者的利益。基于合理假设，本文将乘客到达大厅的时间间隔定为排队问题中服从指数分布的随机变量，并利用计算机产生服从参数为0.1252的指数分布的取值作为相连乘客的间隔时间。针对问题1—4，通过计算机模拟电梯系统，并作统计，对1—4问均作出了准确的回答： 一位乘客等待电梯的平均时间为10.41s，最长时间为36.45s；在电梯中等待的平均时间为48.74s，最长时间为159.95s；从进入大厅直到他到达所要到的楼层的平均时间为49.58s，最长时间为159.95s。4部电梯停留次数分别为226，220，209，214次；实际使用时间的百分比分别为34.586%，34.925%，35.119%，35.878%。

另外，在模拟电梯系统之后，本模型统计出了一些具有参考价值的数据，例如，对于4部电梯，其载客数量经常为1—3人。为电梯管理者改进电梯运行方案提供了十分重要的依据。

在调整电梯系统中，考虑了三个方案：方案一：低楼层的乘客用楼梯来代替电梯；方案二：对电梯关门时间进行；方案三：对电梯的实际运送重新进行安排，分配，建立三个对比模型（1、将楼层进行奇偶分层；2、将楼层按电梯数进行分层；3、将楼层按电梯数的间隔来分层），并对第三个方案进行了具体的求解，最终得出模型三的方案是最优的。

关键词：计算机模拟，排队问题，对比模型

1

一、问题的提出

城市繁华地区有一座12层的写字楼，共有800人在此上班工作，在高峰时间7：50-9：10，人们进入一楼大厅并乘电梯到所在的楼层，有4部电梯为大楼服务，乘客到达大楼的时间间隔在0-30秒内随机变化，到达后每个乘客进入第一部可乘的电梯（1-4号），当某人进入电梯后并选择达到楼层后，电梯在关门前等待15秒，如果另一个人在15秒内到达来，这种等待将重新开始，如果15秒内无人到达，电梯就把全体乘客送上去。假定中途没有其他乘客要上电梯。送完最后一个乘客后，电梯回到大厅，途中也不上客人。

一部电梯的最大容量为12人，当一位乘客来到大厅，没有电梯可乘，就开始在大厅排队等待。

写字楼的管理者希望提高优质服务，但目前有些乘客抱怨在电梯回来之前，他们在大厅等待的时间太长，也有人抱怨他们在电梯里呆的时间太长，还有人说高峰时间大厅太挤，实际情况如何呢？首先对该写字楼电梯系统做理论分析，然后用计算机模拟电梯系统，回答下列问题，：

（1）如果一个人的等待时间是他在排队的时间，即从到达大厅到进入一部可乘电梯的时间，问一个人在队中等待的平均时间和最长时间是多少？

（2）一位乘客实际上呆在电梯中的平均时间和最长时间是多少？

（3）如果运送时间是一位乘客从到达大厅到他到达要去的楼层的时间（包括等电梯的时间）,问平均运送时间和最长的运送时间是多少？

（4）每部电梯停多少次？早高峰时间每部电梯实际上使用时间的百分比是多少？

（5）为了减少乘客的排队等待时间及运送时间，提高电梯的使用效率，应如何来安排、调整电梯运行系统？

二、问题的分析

对于前四个小题的求解，由于电梯系统中乘客的到达和乘客的运送都是复杂的随机过程，因此，电梯系统是一个极为典型的提供成批服务的随即服务系统，对于此类问题的可以用蒙特卡洛方法进行计算机模拟来解决。

若假设乘客进入一楼大厅的时间在高峰时间7：50-9：10是均匀分布的，则可以用计算机随机产生乘客进入一楼大厅的时间，并对其进行排序，由此求出相邻两个时刻之间的间隔t，并求出其均值，通过理论分析可知t服从参数为1的指数分布。计算机按求得的指数分布，随机产生乘客到达一楼大厅的时间间隔，经过累加得到乘客到达一楼大厅的时刻。

结合题中给出的信息，假设800名工作人员的上班地点均匀分布在2-12楼，由电脑按均匀分布产生乘客的目标楼层，再由电脑模拟电梯系统，统计数据，即可解出前四个小题。

对于第五小题，考虑到高峰期电梯运送的不平衡，对电梯运送模式重新进行定义：如奇偶分层，按电梯数进行分类等等。如果可以的话对电梯的软件控制进行修改（电梯关门的时间）。楼层低的乘客还可以通过走楼梯来缓解这个压力。

2

三、基本假设

1.电梯限载12人，但对总载重无； 2.乘客都是单独进入大厅的；

3.乘客的工作单位在2-12楼层均匀分布；

4.在目标层已选的情况下，乘客进入电梯后不再重复选择；

5.若第十二个乘客进入电梯并选择新楼层，则电梯需再等15秒； 6.若有多部可乘电梯，乘客作出任意选择。

四、定义符号说明

1.h：楼层的高度；

2.v：电梯匀速运行时的速度； 3.Ti：i方案的总耗时；

4.xm：电梯实际上经过m层楼的次数； 5.mij：第i种方案中第j部电梯的目标层； 6.Tp：乘客出电梯的时刻；

7.Tm：从第i种方案中第j部电梯的目标层启动的时间；

ij8.ti：加速到第m层楼的时间。

五、模型的分析、建立

5.1指数分布的参数确定[1]

乘客进入一楼大厅的时间在高峰时间7：50-9：10是均匀分布的，首先用计算机随机产生乘客进入一楼大厅的时间，并对其进行排序，由此求出相邻两个时刻之间的间隔，得到时间间隔的均值。

1乘客到达大厅的时间间隔服从参数为的指数分布（其中），设乘客

到达大厅的时间间隔为X，则随机变量X的密度函数为

ex,fx0,x0;x0,

根据随机变量X的密度函数，用计算机模拟出乘客到达一楼大厅的时间间隔，经过累加得到乘客到达一楼大厅的时刻。

5.2 计算机模拟流程图

3

开始 i1;j1 j进入大厅 第位乘客选择电梯 jj1 N iii1 4 Y 有无其他排队乘客 1 2 3 电梯1 电梯2 电梯3 电梯4 Y 目标层是否已有人选 N 重新计时 电梯是否已满12人 N Y Y 载送 15s内有无新乘客 N Y 有无新乘客进入大厅 N 结束 图1 计算机模拟流程图

4

在计算电梯运行时间时，本文将电梯运行过程分为：匀加速运动，匀速运动，匀减速运动三部分，查找相关资料得到具体参数如下：

电梯匀速时的速度[2]v2.5ms； 加速度[3]a1.2ms2； 楼层高度[4] h3.85m。

根据流程图，用MATLAB进行编程，通过计算机模拟电梯系统，并作统计，即可得出1-4题的答案。

5.3 电梯运行方案优化

方案一：低楼层的乘客用楼梯来代替电梯，这样能减少一部分的乘客的等待时间，使电梯的运送效率增加。从电梯角度来说，低楼层电梯的速度比较慢，时间较长，而如果从一定楼层开始运送后，能使电梯速度增加，乘客等待的时间可以减少。

方案二：电梯的关门需要15秒，发现耗时太长，如果可以的话对电梯的软件控制进行修改（电梯关门的时间）。通过摄像头等操作，当发现乘客饱和时即可关门。

方案三：对电梯的实际运送重新进行安排，分配，建立三个模型，并和电梯前几问中的运送模型进行对比。考虑电梯运送的时间和出电梯的时间，以及电梯下楼的时间。得出一个总时间Ti（i=1,2,3,4）。并通过这个总时间的对比，得出电梯用哪种方案比较合适。由于电梯速度的加快，可以间接地减少等待时间。

模型一：将楼层进行奇偶分层，然后求出此模型的总时间：

T1x1t1xmtmTm11Tm12Tm13Tm14xmTpm2m11111

模型二：将楼层按电梯数进行分层（例如：1—3层用电梯1），然后求出此

模型的总时间：

T2xmtmx1t1Tm21Tm22Tm23Tm24xmTpm2m11111

模型三：将楼层按电梯数的间隔来分层（例如：1—5—9用电梯1，2—6—10

用电梯2），然后求出此模型的总时间：

T3xmtmx1t1Tm31Tm32Tm33Tm34xmTpm2m11111

原问题中的模型：

T0x1t1xmtmTm01Tm02Tm03Tm04xmTp

m2m11111六、模型的求解

1、计算机模拟求解

通过计算机模拟，得到各乘客进入电梯时刻。其中第101至第200位乘客具体情况如下表所示：

5

乘客 进入大 编号 厅时刻 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136

655.41 656.16 683.26 688.49 693.11 711.66 725.5 743.73 755.72 758.68 765.08 767.77 775.43 780.02 788.04 798.32 800.19 805.73 808. 823.28 826.03 828.68 845.23 846.02 848.04 858.72 860.56 870.6 883.01 0.94 908.79 918.25 918.81 922.66 939.57 945.11 表1 部分乘客乘电梯的详细情况 是否 进入电 目标 乘客 进入大 进入电 需要 梯时刻 楼层 编号 厅时刻 梯时刻 排队 655.41 6 否 151 1076.3 1076.3 656.16 7 否 152 1085.3 1085.3 683.26 2 否 153 1087.3 1087.3 688.49 12 否 154 1095 1095 693.11 8 否 155 1102.2 1102.2 713.21 12 是 156 1108.3 1108.3 725.5 10 否 157 1109.9 1109.9 743.73 10 否 158 1114.9 1114.9 755.72 9 否 159 1122.3 1122.3 758.68 9 否 160 1123.1 1123.1 765.08 5 否 161 1123.4 1123.4 767.77 2 否 162 1123.4 1123.4 775.43 9 否 163 1147 1147 780.02 7 否 1 1162 1162 788.04 3 否 165 1163.2 1163.2 798.32 9 否 166 1172.6 1172.6 800.19 5 否 167 1182.5 1182.5 805.73 10 否 168 1185 1185 808. 7 否 169 1188.4 1188.4 827.74 6 是 170 1193.4 1193.4 827.74 2 是 171 1199.8 1199.8 828.68 3 否 172 1203.9 1203.9 845.23 10 否 173 1209 1209 846.02 2 否 174 1209.5 1209.5 848.04 8 否 175 1214.9 1214.9 858.72 10 否 176 1214.9 1214.9 860.56 5 否 177 1217.1 1217.1 870.6 6 否 178 1231.3 1231.3 883.01 12 否 179 1248.2 1248.2 0.94 8 否 180 1255.2 1255.2 908.79 7 否 181 1261.7 1261.7 918.25 11 否 182 1268.6 1268.6 918.81 5 否 183 1271.3 1271.3 922.66 9 否 184 1283.2 1283.2 940.17 5 是 185 12.1 12.1 945.11 6 否 186 1290.5 1290.5 6

是否 目标 需要 楼层 排队 3 否 12 否 7 否 3 否 9 否 10 否 3 否 12 否 11 否 2 否 5 否 11 否 7 否 6 否 10 否 12 否 12 否 6 否 10 否 4 否 6 否 6 否 2 否 2 否 4 否 6 否 4 否 12 否 4 否 9 否 6 否 10 否 12 否 3 否 10 否 3 否 137 967.23 976.71 5 是 187 1297.1 1297.1 5 否 138 9.4 9.4 2 否 188 1313.4 1313.4 8 否 139 1013.4 1013.4 5 否 1 1317.1 1317.1 6 否 140 1013.9 1013.9 8 否 190 1318.2 1318.2 11 否 141 1021.2 1021.2 4 否 191 1324.3 1324.3 8 否 142 1039.8 1039.8 9 否 192 1354 1354 6 否 143 1043.6 1043.6 3 否 193 1357.2 1357.2 11 否 144 1043.8 1043.8 6 否 194 1366.9 1366.9 12 否 145 1050.9 1050.9 12 否 195 1369.6 1369.6 11 否 146 1053.5 1053.5 7 否 196 1376.1 1376.1 10 否 147 1056.3 1056.3 9 否 197 1377.3 1377.3 2 否 148 1058.4 1058.4 6 否 198 1386.3 1386.3 11 否 149 1065.8 1065.8 3 否 199 1387.6 1387.6 2 否 150 1066.8 1066.8 3 否 200 1395.9 1395.9 10 否 可以看出，该写字楼电梯系统运行状况良好，以上100名乘客中只有5名乘客需等待电梯。

4部电梯该日高峰期总载客次数统计如下：

图2 四部电梯载客次数

第一部电梯载客数量频数 第二部电梯载客数量频数

7

第三部电梯载客数量频数 第四部电梯载客数量频数

图3 四部电梯载客数量频数分布图

下表列出了电梯1高峰时期共49次载客的具体停留楼层与停留次数：

表2 电梯1停留楼层详情 载客 234567101112停留 次数 楼 楼 楼 楼 楼 楼 楼 楼 楼 楼 楼 次数 1 1 0 0 0 0 0 0 0 0 0 0 1 2 0 1 1 0 0 0 0 0 1 0 0 3 3 0 0 1 1 0 0 0 0 0 0 1 3 4 0 0 1 0 0 1 0 0 0 1 0 3 5 0 1 1 0 1 1 0 0 1 1 1 7 6 0 0 1 0 0 0 0 0 1 0 1 3 7 0 0 0 0 0 1 0 0 0 0 0 1 8 1 0 0 0 0 0 1 0 0 0 1 3 9 0 0 0 0 0 0 0 1 0 0 0 1 10 1 1 0 0 1 0 0 0 0 0 0 3 11 0 0 0 0 0 0 0 0 0 1 0 1 12 1 0 0 1 0 0 0 0 0 0 0 2 13 0 1 0 0 0 1 0 1 0 0 0 3 14 0 0 0 0 1 1 0 0 0 0 0 2 15 0 0 1 0 1 0 0 1 1 0 1 5 16 1 0 0 0 0 0 0 0 0 1 0 2 17 0 0 0 0 1 0 0 0 1 1 1 4 18 0 0 1 1 0 0 0 0 0 0 0 2 19 1 1 0 0 1 1 1 0 1 0 0 6 20 0 0 0 0 0 0 0 0 0 1 0 1 21 0 0 0 1 0 1 0 1 0 0 0 3 22 0 0 1 0 0 0 0 0 0 0 0 1 23 0 1 0 0 0 1 1 0 1 1 0 5 24 0 0 0 0 0 0 0 1 0 0 0 1 25 1 0 0 1 0 0 1 1 0 0 0 4 26 0 0 1 0 0 0 1 1 0 0 0 3

8

27 0 0 1 0 0 0 0 0 0 0 0 28 0 0 1 0 0 0 1 0 0 1 0 29 1 0 0 0 0 0 0 1 0 0 0 30 0 1 0 0 0 0 0 0 1 0 0 31 0 0 0 0 0 1 0 0 1 0 0 32 1 0 0 0 0 0 1 0 1 0 0 33 0 0 0 1 1 0 1 0 0 1 0 34 0 1 0 0 0 1 0 0 1 1 0 35 0 0 0 0 0 0 1 0 0 0 0 36 0 0 0 1 0 0 1 0 0 0 0 37 1 0 0 0 1 1 1 0 0 1 1 38 0 0 0 0 1 0 0 1 1 0 0 39 1 0 0 0 0 0 0 0 0 0 0 40 0 0 0 0 0 0 0 0 1 0 0 41 0 0 0 0 0 0 1 1 0 1 1 42 0 0 0 0 0 0 0 0 1 0 0 43 1 0 1 0 0 0 1 0 1 0 1 44 0 1 1 0 0 1 0 1 1 0 0 45 1 0 0 0 0 0 0 0 1 0 0 46 0 0 1 0 1 0 0 1 0 0 0 47 0 0 0 0 1 1 0 1 1 1 0 48 1 1 1 1 0 1 0 1 1 0 1 49 1 1 0 1 1 1 1 1 1 1 1 其中，1代表相应楼层要停留，0代表不停留。

另外，四部电梯相连两次停留楼层差分布情况如下图所示：

1 3 2 2 2 3 4 4 1 2 6 3 1 1 4 1 5 5 2 3 5 8 10

图4 四部电梯相连两次停留楼层差分布图

可以看出，该写字楼电梯载客时常出现只上升1层的状况。

对于问题1，模拟某一天电梯系统结果如下：

总等待人数：54人；

9

总等待时间：562.15s； 平均等待时间：10.41s；

平均等待时间（含不需要等待的人）：0.84s； 最长等待时间：36.45s；

模拟100天后各自取平均值结果如下：

平均等待时间：10.51s 最长等待时间：33.13s

对于问题2、3、4，模拟得出结果如下：

乘客实际上位于电梯内的平均时间为48.74s，最长时间为159.95s； 运送时间平均值为49.58s，最长时间为159.95s；

第1-4部电梯停留次数分别为226，220，209，214次；使用时间分别为1660.1，1676.4，1685.7，1722.1s；使用时间百分比分别为34.586%，34.925%，35.119%，35.878%.

2、三种方案求解

考虑到电梯在运行的过程中为先加速再匀速最后再减速，对电梯运行时间进行了求解，并做如下规定，只有一层时，只进行加速，和减速，没有匀速过程，而两层或两层以上时，就开始有匀速加速的过程了。

根据这个对电梯运行的时间进行了求解： 只有一层时：

1212hatat112222a1t1a2t22a2h2a1ht12 解得：ttt112a1a2a1a1a2a22两层或两层以上时：

；

1212mhatvtat1342522atvmhvv13 解得：tmv2a12a2a2t5vt4t5tmt3由上述分析可得如下不等式：

；

t1t2t3t4t5t6t7t8t9t10t11； （1） 12345671011对方案三中的式子进行联立得到：

10

1111Tm02Tm03Tm04xmTpT0x1t1xmtmTm01m2m11111Tm12Tm13Tm14xmTpT1x1t1xmtmTm11m2m11111TxtxtTTTTxT （2）

mm11m21m22m23m24mp2m2m11111TxtxtTTTTxT43mm11m31m32m33m3mpm2m1因为各种方案得到的总的楼层数S是一样的，所以对每个乘客对应的楼层数ym有：

y 方案 电梯 电梯1 电梯2 电梯3 电梯4

mS （3）

由于方案的不同，电梯在实际运行中xm的次数是不同的，联立（1）（2）（3）得：T3T1T2T0，所以得出的新的运送方案为：

新运送方案 1——5——9 2——6——10 3——7——11 4——8——12 七、结果分析

计算机模拟并没有理论计算精确，但其模拟结果更符合实际情况。实际上，对于电梯运行方案在理论上进行优劣评判是十分困难的，因此，能否通过计算机模拟电梯运行过程，得出实际而且可靠的结果显得尤为重要。

本模型对模拟电梯系统运转设计了合理而又正确的流程图，不仅对问题1-4作出了的回答，同时保留了模拟过程中的各个具体数据，通过对数据的统计，可以得出许多具有参考价值的结果，为电梯管理员优化电梯运行方案提供了重要的参考依据。

而问题5求解的答案和运送方案在实际中也可以实现的，可以提高电梯的效率。

11

八、模型推广、评价与改进

此问题还可以推广到m个电梯，n层楼的情况，也可以通过建立类似于问题5中的模型来进行分类。从而提高在高峰期电梯的使用率，使乘客排队时间和运送时间缩短。

而在本题中，假设了乘客在每层楼都是均匀分布的，这简化了这道题的求解，但和实际情况相比，有些不太相符。当考虑到实际情况时，每层楼的人数不一定是统一的，而且只考虑上楼人数，中途不上楼，但实际中更多的可能确实中途上楼。这时可以借助计算机模拟来进行具体的模拟，求解。而且根据实际人数的不同，最终得到的方案可能也会不同。

也可以把模型推广到对新建写字楼的电梯交通设计中，用计算机计算出不同数量和不同额定容量的电梯系统在不同的服务强度下的平均队长、平均等待和服务台的平均占有率，并根据结果画出对应的曲线，得到较优的设计。

现实中，电梯运行高峰期不止只有这么一段时间，还有下班的高峰期和平常期，可以根据实际情况，不同时刻，对电梯的运行采用不同的方案。

参考文献：

[1]茆诗松，《概率论与数理统计教程》，高等教育出版社，2004.7； [2]中国菱电股份有限公司，《多层次微电脑省电智慧型电梯》，

http://www.tmec.com.tw/images/p3/Download/HGPS4-5.pdf，2010.8.16； [3]佚名，《交流双速电梯的加速,启动,满速,减速过程!》，

http://zhidao.baidu.com/question/76609373.html，2010.8.16； [4]远洋·光华国际，《楼盘信息》，

http://www.wolai.com/newhouse/15388/material/，2010.8.16. [5]姜启源，《数学模型》，高等教育出版社，2003.8.

12

附件

附件一：确定指数分布均值的MATLAB程序代码

for j=1:1000

x=4800*rand(1,600); x=sort(x); for i=1:599

derta(i)=x(i+1)-x(i); end

miu(j)=mean(derta); end

miu0=mean(miu)

附件二：模拟电梯运行的MATLAB程序代码

clear;clc miu0=7.9869; i=1; time=0;

while time<4800

x(i)=exprnd(miu0); while x(i)>30

x(i)=exprnd(miu0); end

time=time+x(i); i=i+1; end

x=x(1:length(x)-1); %hist(x,40)

mean(x) %由于乘客到达的时间间隔不允许超过30s，所以实际期望值小于miu0

d=length(x) %实际乘客总数大致600-700人之间

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

x=cumsum(x)'; for i=1:d

x(i,2)=floor(2+11*rand); %确定第i位乘客目标楼层 end

13

flag=0;

y=zeros(4,3);

i=1;g=ones(1,4);g1=ones(1,4); louceng=zeros(4,12); for j=1:d for m=i:j

for k=randperm(4) if y(k,1)>x(j,1) continue; end

if x(m,1)x(m,3)=max(x(m,1),y(k,1)); %进入电梯时刻 go(1:3,g1(k),g(k),k)=[m;x(m,2);x(m,3)]; g1(k)=g1(k)+1; i=i+1; flag=1;

y(k,3)=y(k,3)+1;

if louceng(k,x(m,2))==0 louceng(k,x(m,2))=1;

y(k,2)=x(m,3); %重新计时 end end

if y(k,3)==12&sum(louceng(k,:))>0 go(4,1,g(k),k)=y(k,2)+15; g1(k)=1;g(k)=g(k)+1; nex=next(louceng(k,:));

y(k,1)=y(k,2)+nex+15; %下一次电梯返回底层时刻

y(k,2)=y(k,1); louceng(k,:)=0; y(k,3)=0; end else

go(4,1,g(k),k)=y(k,2)+15; g1(k)=1;g(k)=g(k)+1; nex=next(louceng(k,:)); y(k,1)=y(k,2)+nex+15; y(k,2)=y(k,1); louceng(k,:)=0; y(k,3)=0; end

if flag flag=0;

14

break; end end end end n=d;

while x(n,3)==0 n=n-1; end

wr=d-sum(x(1:n,1)==x(1:n,3)) wt=sum(x(1:n,3)-x(1:n,1)) wait=wt/wr

waitmax=max(x(1:n,3)-x(1:n,1))

附件三：问题（2）、（3）及问题（4）前一问的MATLAB程序代码

load data for k=1:4 num=1;

for l=1:g(k)-1

ceng=zeros(1,12); n=1;

while n<=12&go(2,n,l,k)~=0 ceng(go(2,n,l,k))=1; n=n+1; end

q=next0(ceng);

ceng42(num,1:12,k)=ceng; num=num+1; for i=1:n-1

x(go(1,i,l,k),4)=go(4,1,l,k)+q(go(2,i,l,k)); end end end n=d;

while x(n,4)==0 n=n-1; end

mean(x(1:n,4)-x(1:n,3)) %乘客实际上位于电梯内的平均时间 max(x(1:n,4)-x(1:n,3)) %最长时间

mean(x(1:n,4)-x(1:n,1)) %平均运送时间 max(x(1:n,4)-x(1:n,1)) %最长运送时间

15

stp=zeros(1,4); for k=1:4

for l=1:g(k)-1 n=1;

while n<=12&go(1,n,l,k)>0 n=n+1; end

stp(k)=stp(k)+n; end end stp

附件四：计算电梯实际上使用百分比的MATLAB程序代码

dianti0=[49 52 54 51]; for dianti=1:4 n=1;

for i=1:dianti0(dianti) k=1;

for j=2:12

if ceng42(i,j,dianti) hld(dianti,n)=j-k; n=n+1; k=j; end end

hld(dianti,n)=k-1; n=n+1; end end

t=[3.204 5.1638 6.7038 8.2438 9.7838 11.324 12.8 14.404 15.944 17.484 19.024]; use=zeros(1,4); for dianti=1:4 m=1;

while m<=size(hld,2)&hld(dianti,m)>0

use(dianti)=use(dianti)+t(hld(dianti,m)); m=m+1; end end use

usep=use/4800

16

附件五：MATLAB程序调用函数1

function y=next(ceng) j=1;

for i=2:12 if ceng(i) x(j)=i; j=j+1; end end

t=[3.204 5.1638 6.7038 8.2438 9.7838 11.324 12.8 14.404 15.944 17.484 19.024]; x=[1 x]; x

for k=2:length(x)

dx(k-1)=x(k)-x(k-1); end

y=sum(t(dx))+8.5*length(dx)+t(x(end)-1);

附件六：MATLAB程序调用函数2

function q=next0(ceng)

t=[3.204 5.1638 6.7038 8.2438 9.7838 11.324 12.8 14.404 15.944 17.484 19.024]; i=1;q(1)=0;flag=1; for j=2:12

if ceng(j)==1 if flag==1 flag=0;

q(j)=q(i)+t(j-i); else

q(j)=q(i)+t(j-i)+8.5; end i=j; end end

17