关于电梯系统优化问题的数学模型

1、关于电梯系统优化问题的数学模型摘要在高层商务楼里，电梯承担着将人和货物运送到各个楼层的任务。在当今社会，工作生活节奏愈发加快，因而电梯系统的运行效率对人们的生活的影响不可忽视。目前的高层商务楼等大多数高层建筑中， 一般都使用单井道单轿厢或者单井道双轿厢两种模式的电梯，本文就结合这两种模式，根据实际情况将问题分为两种情况考虑，重点讨论了将电梯运行效率最大化的方法，建立了相关模型，并给出了相应的优化参数。本文将电梯系统的优化分为高峰期和非高峰期两种时期进行讨论。 高峰期时通过对问题的分析， 发现可以设置电梯区间以尽可能减少目标层较高的乘客占用目标层较低的乘客的电梯资源，根据这一思想，我们将其简化为

2、排队问题来考虑，并据此建立了排队模型， 通过实地统计数据以及C 语言的编程， 能够较好地解出模型，得到在高峰期时将一部分电梯区间的顶层设为第 14 层左右的优化方案。非高峰期时通过对这一时期特点的分析，以每台电梯在无乘梯需求时自动停留的楼层为着眼点，采用枚举的方法编程求解，得到在非高峰期将电梯均匀分布在楼层中的优化方案。最后，我们对模型参数进行了灵敏度的分析，发现虽然模型对数据的依赖性较强，但最优方案不随参数的波动而变化，所以这个结果还是可信的。本文提出的方案直观易行，且几乎不需额外的经济投入，可行性很强，具有较好的 参考价值。一问题重述在高层商务楼里，电梯承担着将人和货物运送到各个楼层的任务

3、。目前的高层商务楼等大多数高层建筑中，主要使用单轿厢和双轿厢两种电梯运行系统。单轿厢电梯在向上运行时，只有满足了所有“上行请求”时才会开始满足“下行请求”，反之亦然；而对于双轿厢电梯，乘客在进入轿厢前就通过按钮面板选择了要停靠的楼层，系统迅速整合分析接收到的流量数据，并调度合适的轿箱来应接乘客。现有一座商务楼， 设计地上层数为 28 层， 地下停车楼2 层， 每层的建筑面积为 1500平方米，楼内有6 个用于客梯的电梯井道。电梯按照商务楼建筑面积15至 20平方米每人的标准来设计。第 1 层的楼层高为 4.8 米，其余层均为 3.2 米，设计电梯的平均运行速度 1.6米/ 秒。我们的任务是：1

4、建立一个合适的单轿箱客梯系统的运行方案，使尽可能地提高电梯系统的运行效率；2分别在运行的高峰期与非高峰期，对双轿箱的电梯系统与单轿箱的电梯系统的运行效率等进行对比分析, 评价两种方案的优劣性，估计双轿厢系统运行效率的提高率。二基本假设1 .电梯载客量为13人，且不超载。13人载客量是国内最常见的一种电梯规格，并且 为了乘梯安全，电梯不应超载。2 .电梯在每层停留的时间相等。在假设 1成立的前提下，电梯乘客可以迅速有序地离 开电梯，电梯停留时间受离开人数的影响可以忽略不计。3 .乘客的到达形成泊松流。4 .商务楼工作人员均匀分布在地上 2层到28层的每一层，即电梯乘客在每一层下电梯 的概率相等。

5、5 .在上班高峰期无人下电梯，在下班高峰期无人上电梯。6 .使用每层地下停车楼的人数相等。三符号及名词说明输入层：有需要乘电梯的人流入的楼层。目标层：乘客想要到达的楼层。服务：在上班高峰期电梯由输入层出发到载完13个人回到输入层称为一次服务。%：=3岗),：第k个电梯或电梯井道的运行区间，即被限制只能从p层运行至U q层。A =(%，%，%，%Qj：高峰期电梯系统运行的一种安排方案。年： 第k个电梯在无乘梯需求是停留的楼层。田31，m个电梯在非高峰期的一种运行方案， m=6或12。f(A)：安排方案A下乘客等待时间的期望。f( :安排方案B下乘客等待时间的期望。W(k): 乘坐第k个电梯的乘客

6、等待时间的期望。入，A:乘客形成的泊松流的强度。t(p,q):电梯从p层运行到q层所用的时间为:电梯在每层停留的时间。t(M :在高峰期第k个电梯完成一次服务所用的时间m 1:使用地下停车楼的人数比例。但2：不使用地下停车楼的人数比例。nA):第k个电梯一次服务中所能运行到的最高层。P(n):在上班高峰期电梯在一次服务中停留n次的概率。四问题分析本题是对电梯系统的优化问题，优化的标准就是找到一种方案 A使所有乘客等待时 问的期望f(A)最小。这里为了叙述方便，将地下 1层、2层分别记为-1层、-2层，地 上1层、2层、28层分别记为0层、1层、27层。我们发现，不管是单轿厢电梯系统，还是双轿厢

7、电梯系统，在上班高峰期，0层、-1层和-2层为输入层，1层至27层为目标层，在下班高峰期，1层至27层为输入层，0 层、-1层和-2层为目标层，也就是说，在高峰期，输入层和目标层分别有所集中；而在 非高峰期，输入层和目标层都是随机分散的。所以，为了合理优化电梯系统的效率，应 把这两种时期分开考虑。4.1 高峰期的分析4.1.1 上班高峰期的分析上班高峰期的输入层为 0， -1， -2 层，则电梯的初始位置只能集中分布在这三层。目标层越大，电梯需要上升的高度就越高，一次服务的时间就会越多。由于乘客想要到达的目标层是随机的，因而一次服务中只要有人的目标层较大，相应电梯的等待人群需要等待的时间就越多

8、，而一些目标层较低的乘客同样需要等待这样的时间，可以理解为高目标层乘客占用了低目标层乘客的“资源” 。这就造成了等待时间的增加。所以我们提出一种电梯区间的思想，即在上班高峰期将每个电梯所能运行的范围加以限制，同时令目标层不同的乘客乘坐不同区间的电梯，这样目标层较低的乘客乘坐区间较小的电梯，等待的时间就会有所降低，而目标层较高的乘客乘坐区间较大的电梯，等待时间影响不大。在这种情况下，单轿厢电梯系统和双轿厢电梯系统的模型一致，考虑到这一过程符合排队过程的特点，可以将其简化为排队模型，并编程求得最优解。4.1.2 下班高峰期的分析下班高峰期的输入层为 1 层至 27 层，目标层为0， -1， -2

9、层，电梯的初始位置无法集中。输入层越高，电梯需要运行到很低的目标层再回到输入层，经过的楼层数越多，所用的时间也就越多。因而只要高输入层的乘客有乘梯需求，那么低输入层的乘客就会大大增加，可以理解为高输入层乘客占用了低输入层乘客的“资源” 。所以沿用 4.1.1 中的思想，利用电梯区间将下班高峰期电梯的运行范围加以限制，同时令输入层不同的乘客乘坐不同区间的电梯，这样输入层较低的乘客乘坐区间较小的电梯，等待时间就会有所降低，而输入层较高的乘客乘坐区间较大的电梯，等待时间影响不大。在这种情况下，单轿厢电梯系统每个输入层都符合排队过程的特点，可将其简化为排队模型;4.2 非高峰期的分析非高峰期的输入层和

10、目标层都是随机分散的，且人流量小，因而不同于高峰期的分析。对于每个单轿厢电梯和双轿厢电梯，其初始位置应在-2层至27层之间，在某一时刻，有人需乘电梯，则他在1层至27层的概率相等，只需简化为安排 6个单轿厢电梯 或者12个双轿厢电梯的初始位置，使乘客等待电梯的时间期望尽可能小即可。这一模 型可以通过编程完成。五模型的建立与求解5.1 单轿厢电梯系统的求解5.1.1 上班高峰期单轿厢电梯系统的求解对于上班高峰期，每个输入层都要有一个区间从本层到27层的电梯以保证乘客能I1到达任何目标层，贝U叽(0,27)丁，叼=(-1,27)/，% = -2,27)同时令出独力)，rr%=( - 2%),I -

11、 Io那么对于每个电梯及其乘客，都可以简化为如图模型111顾客随机到达顾客 排队服务机构(服务时间随机)顾客离开其中电梯为“服务机构”，且服务时间随机，乘客被送往目标层后可视为“顾客离开”，则这一模型与排队模型类似，但排队模型中服务机构是从等待的顾客中随机取其一进行服务121 o为了使模型与排队模型相符，这里把13个乘客看作一个“乘客集合”，则“乘客集合”输入的泊松流强度为耳此时模型符合排队模型，且符合 M/G/1排队，可用 排队论公式求解。对于输入层为0层的吗，t(吗)为电梯停留所用时间与电梯运行所用时间之和，电梯运行所用时间为2(2N(z)+1)=4N(0)+2,电梯停留所用时间为n%P(

12、n),其中Q 3M X4；nC 1,min13,N( , P(n)=%*, Q(13,n)为把 13个人分为 n 组的可能数。则Q6声)x4：t q。 a o 。13t( )=4N( ) +2+ n由排队论公式，乘第2个电梯的乘客等待时间的期望W-2尸 改),(尸，方(3)且 W(5)=W(%)(1=27)。对于输入层为0层，当=0,乘坐2号电梯的概率为0,当W=27,乘坐2号电梯的概率为1/2,假设次概率服从线性关系，则乘坐 2号电梯的概率为篇，那么乘坐1、2号电梯的乘客等待时间的期望为W(aVa2it)=W(电)+(1西W(%)+狈()+(1-)2(1叫A,故同时，记A为所有乘客到达的泊松

13、强度，则乘 1、2号电梯乘客的泊松强度为 1、2号电梯“乘客集合”的泊松强度分别为为了解出模型，我们需要(JA和叫三组参数=(1国对于,我们实地做了实验，统计记录下了一组电梯停留时间的数据，如图所示：19我们发现，数据大致都集中在一条平行于 x轴的直线上，对数据求均值得=6 =6.7s。对于3，我们找到了一家与问题中商务楼规模类似的公司，调查得到开车上班的人所占比例为42.3%,这里认为%=42.3%, =57.7%对于A,我们同样是在这家公司大厅实地做了统计，得到30分钟内到达329人，这 里认为 A= 0.183。取力=1 , 227,得到W(优也)与灯的关系如图00从图中可以看出，当力=

14、14时，W(4%)最小，即-1：/)=27 141时为最优方案同样，对于输入层为-1层,+(1-)t7(13用 x 血；(12 %A外32A且 t(%)=4N(%) +4+ n 或 ,乙=(1-耘)而，4f 区, 得到W(%，比4)与明的关系如图从图中可以看出，当华=14时,W(%)最小,即网R尸2714时为最优方案。对于输入层为-2层，有W(尸口(1 A)又从： 叮CZ (X% 杼 13且 t( 6)=4N( 6) +6+ n+(1-)得到W(%)与的的关系如图2 00320i SO- 2 - 2从图中可以看出，当的=14时，W（%）最小，即（%）=2714 1时为最优方案。001 1 2.

15、 2于是我们得到，当A=1427 142714 时，f（A）最小，为如Z侬2f（A）=%W（%七|） + W W（%声4） +2 W55，/）= 33.34。5.1.2下班高峰期单轿厢电梯系统的求解对于下班高峰期，每个目标层都要有一个区间从本层到27层的电梯以保证任何输入层的乘客都能到达目标层，则=（0,27）丁，%=（- 127）7，% = - 227），同时令.丁，F|p%=（。闯1）%=（ - 1 闯2）%1- 2%）C对于每个输入层的乘客，都有刚好没乘上电梯的乘客需要等待电梯一次服务之后才可以接受服务，和5.1.1类似，同样符合排队模型的特点。将乘坐同一电梯的各输入层的乘客合在一起看作

16、同一个排队，并且将 13个乘客视为一个“乘客集合”，则该模型可 简化为排队模型，并且和 5.1.1的模型完全相同。参数方面，、和但】应当保持不变，而 A则会发生变化，于是我们在同一家公司于下班高峰期做了统计， 得到30分钟离开391人，这里认为A = 0.217。故我们得到Wd，优2）与囚、w（q，%）与牝、WC%）与心的关系分别如图X BOO IL &-OCJ 工4。口 1=0*0 :Ld=cra 。口口 &-O-O 4 30 工口0O。之。1,0 0由图可知，当1=13时，W(r%)最小，即0 %)=27 131时为最优方案。_ 1_ li由图可知，当取=14时，W(3%)最小，即(%)=

17、 27I”时为最优方案。-2 - 2由图可知，当心=14时，W(%6)最小，即(5%)= 27141时为最优方案。00-1-1- 2- 2于是我们得到，当A=271327 142714 时，解)最小，为% 力f(A)=1W(5)+2W( J 4) +2 W( 5，6) = 45.06。5.1.3 非高峰期单轿厢电梯系统的求解非高峰期的输入层和目标层都是随机分散的，且人流量小，因此不应分析电梯的区 间安排，而应从电梯在无乘梯需求时自动停留的位置入手分析 如4.2所说，记 后电电也，与如织)，设某乘客所在楼层为n,则他所要等待的时 间为mint(l,n)(i=1,2,3,4,5,6)。并且我们认为

18、此乘客在-2层到27层的概率相等，故等待时间的期望271f( b)=E/l =- 2而11s 他再) ,(i=1,2,3,4,5,6)通过编程枚举，可以得出，当 后(12*2712.17.22)时，f(B)最小，为f(护“=2而皿血：3四：=2.47。5.1.4 模型结论至此，我们得出了单轿厢电梯系统运行效率最优化的运行方案，即在高峰期采取方001 1 2 2|案A=271427 142714,上班时乘客等待时间的期望为33.34s,下班时等待时间的期望为45.06s ;非高峰期采取方案“-2,2712,17,22)，等待时间期望为2.47so5.2双轿厢电梯系统的求解5.2.1 上班高峰期双

19、轿厢电梯系统的求解对于上班高峰期，每个输入层都要有一个区间从本层到27层的电梯井道以保证乘客能到达任何目标层，和 5.1.1类似，令同一井道内两个电梯的区间相同，这样可以避免控制台的混乱，则 叫=(0,27), %=(-1,27/，工(-2,27)，同时令心=。/),7T%_( 一 1/)优(- 2%)=，-o此时，同一井道内两个电梯一次服务一共可以运载26个人，这里把26个乘客视为A一个“乘客集合”，相应的泊松流强度为赤，则此模型可以简化为排队模型。同5.1.1,我们得到,?4加/(%) +见限电) 比 :铲。(叩)+仇(丐)W(尸+(1-)心kOjA也叫八1=(1-而恢，氏=示故我们得到W

20、(町丐)与乞的关系如图302 52000由图可知，当1=12时，W(r%)最小，即-1 %尸27 121时为最优方案。同理可得W(5%与宝、W(%，%与的的关系分别如图AQ3 0由图可知，当牝=14时，W(34)最小，即(5区*)=-1 -114 时为最优方案。0 0 -1-1 - 2 - 2于是我们得到，当A=时,f（A）最小,为27 12 27 14 27 14f(A)=】W(叼知 +W（吗/ 4）+ 灸 W5，%）= 12.24。5.2.2下班高峰期双轿厢电梯系统的求解对于下班高峰期，每个目标层都要有一个区间从本层到27层的电梯井道以保证任何输入层的乘客都能到达目标层，则 =（。，27）

21、丁，%4-1,27）, %=-2,27）、同时令 y%=（0闻1）值4=（ - 1 闯2）%-2%） .同5.1.2,将26个乘客视为一个“乘客集合”，则此模型可简化为排队模型，参数中的泊松流强度沿用5.1.2中的A。故我们得到W（外电）与支、W（%与叼、W（%/）与的的关系分别如图曰口3口时为最优方案。由图可知，当1=12时，W(r%)最小，即00 .17 12)Q27141由图可知，当牝=14时，W(34)最小，即 (*)=时为最优方案。-2 - 2由图可知，当心=14时，W(%)最小，即-5,%)=27141时为最优方案。001 1 2 2于是我们得到，当A=271227 142714

22、时，f(A)最小，为始)=叫亚(&)+W(%)+Tw(u5(ug)= 15.24。5.2.3 非高峰期双轿厢电梯系统的求解非高峰期的输入层和目标层都是随机分散的，且人流量小，因此同一井道中的电梯在无乘梯需求时自动停留的位置可以不同，则同 5.1.3,记 片血%J,设某乘 客所在楼层为n,则他所要等待的时间为mint(a,n)(i=1,2,3,，校)。并且我们认为此乘 客在-2层到27层的概率相等，故等待时间的期望股嬴n附，叫 g,2,3,.,i2)通过编程枚举可以得出，当 芹(-2,025710,12,15/7,202326)时，f(份最小，为E乙焉回叫叽I.5.2.4 模型结论至此，我们得出

23、了双轿厢电梯系统运行效率最优化的运行方案，即在高峰期采取方001 1 2 2案A=271427 142714,上班时乘客等待时间的期望为12.24s,下班时等待时间的期望为15.24s；非高峰期采取方案 芹(-2,025710,1215,1720,23,26)、等待 时间期望为1.33s。六模型的比较6.1 高峰期电梯系统效率的比较上班高峰期，双轿厢电梯系统平均等待时间为12.24s,单轿厢电梯系统平均等待时间为33.34s,双轿厢电梯系统比单轿厢系统效率提高了 -1224X100%=172.4%;下班高峰期，双轿厢电梯系统平均等待时间为15.24s,单轿厢电梯系统平均等待时间为45,06-1

24、5.2445.06s,双轿厢电梯系统比单轿厢系统效率提高了 一国一X 100%=195.7%。1.33s,单轿厢电梯系统平均等待时间为2.17- L33 X100%=85.7% 。6.2 非高峰期电梯系统效率的比较非高峰期，双轿厢电梯系统平均等待时间为2.47s,双轿厢电梯系统比单轿厢系统效率提高了七 模型的灵敏度分析因为本文的模型所需参数几乎都是通过小范围的统计得到，因此还需考虑参数波动对模型结果的影响。先考虑A的波动对结果的影响。这里将A的值作 0.05的波动，得到等待时间期望 随楼层的变化，结果如图toW=O.577A=O 2 3SW=a.577A=0.1B3UJ = n E77 A=Q

25、 13 3J41 I I I I I I I I I I I I I I I I I I I I I Illi I I125 4 5678 9 10 1112131-4 15-1617 18 19 2021 22 23 24 2 5 26 27我们发现虽然期望值均随A的波动而变化，但整体增减趋势没有改变再考虑叫的波动对结果的影响。这里将叫的值作 0.05的波动，得到等待时间期望 随楼层的变化，结果如图我们发现虽然期望值均随 的波动而变化，但整体增减趋势没有改变所以我们认为模型结果是可信的八 模型的优缺点8.1 模型优点本模型最显著的优点就是简单直观， 能很好地借助现有模型对问题进行分析和求解，

26、便于编程计算。同时将问题根据实际情况作不同考虑，建立不同的模型，使结果更具实际参考意义，而且提出的解决方案简单易行，在经济上几乎不会造成额外的支出，可行性很强。8.2 模型缺点与本模型最显著优点简单相伴而来的缺点就是参数过多，对数据的依赖性强，需要统计大量的真实数据才能更加准确地求解模型，而由于时间有限，我们这里统计的数据量还不够，参数的波动虽然对方案整体设计基本上没有影响，但对相关的数据结果可能会造成一些影响，还需要进一步加大数据统计量，以对模型作进一步完善。参考文献【 1】 张莹，运筹学基础，北京：清华大学出版社，2010。【 2】 宋荣兴，孙海涛，运筹学，北京：经济科学出版社， 2011

27、。【 3】 孟玉柯，排队论基础及应用，上海：同济大学出版社，1989。附录1 .求解期望值的C 语言程序#include #include #include #include #define T1 2#define T2 6.7#define W 0.577#define A 0.183/* 分别取 0.183 和 0.217 进行计算 */intmcn(int n, int m);intzuhe(int m, int n);intjiecheng(int m);int min(int a, int b);int main(intargc, char *argv)int i, j, p, q;

28、p = 0; /* p 取-2 -1 0 分别计算 */for (q = 1; q 28; q+)int TT = 0;int TN = 0;int data3030;/* data ml【n】代表电梯最高向上走m层返回，其中n代表运行过程中电梯停留的层数所对应的次数，若上升一层的时间是T1 ，停留一层的时间是T2 ，则其对应的时间是 2*m*T1+n*T2 */memset(data, 0, sizeof(data);for (i = 1; i q - p + 1; i+)for (j = 1; j min(i + 1, 14); j+) dataij =mcn(13, j) * zuhe(

29、i - 1, j - 1);/* 下一步应该是确定对应的时间中每个时间所对应的次数来计算平均数和方差*/for (i = 1; i 30; i+)for (j = 1; j 30; j+) if (dataij != 0)TT +=(2 * i * T1 +j * T2) * dataij;TN += dataij;floatave = 0.0;ave = TT / TN;floatvar = 0.0;float temp = 0.0;for (i = 1; i 30; i+)for (j = 1; j 30; j+)if (dataij != 0)temp +=dataij * pow(2

30、* i * T1 +j * T2 - ave, 2);var = temp / TN;float E = 0.0;E = ave;float D = 0.0;D = var;float lmt1 = 0.0;lmt1 = (1 - (float) q / 27) * W * A / 26;float lmt2 = 0.0;lmt2 = (float) q *W * A / 702;float key = 0.0;key =(float) q / 54 * lmt2 * (E * E + D) / (1 - lmt2 * E) + (1 - (float) q / 27) * lmt1 * 5287 / (1 - 101 * lmt1);printf(q=%d 平均数 =%f 方差 =%f %f %f 结果 =%fn, q, E, D, lmt1, lmt2, key);system(PAUSE);return 0;intmcn(int n, int m)int f2020, i;for (i = 1; i = 15; i+)for (m = 1; m = i; m+)if (m = 1 | i = m | i = 2)fim = 1;elsefim =fi - 1m - 1 + m * fi - 1m;return fnm; intzuhe(int m, int n)re