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

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

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

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

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

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

6、:乘坐第k个电梯的乘客等待时间的期望。入，A:乘客形成的泊松流的强度。t(p,q):电梯从p层运行到q层所用的时间?(?电梯在每层停留的时间。t(?9 :在高峰期第k个电梯完成一次服务所用的时间。?:使用地下停车楼的人数比例。?2 :不使用地下停车楼的人数比例。N(?) :第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层，

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

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

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

12、(?)+2,电梯停留所用时间为 n?0?P(n),其中 nC1,min13,N( ?), P(n尸?，母：；嚼,Q(13,n)为把13个人分为n组的可能数。则?13,?) X 需t(?)=4N(?) +2+ n?7-?由排队论公式，乘第2个电梯的乘客等待时间的期望 且 w(?)=w(?)(?=27)oW(?)=?吊+?72?(耳?2)2?(1- p)(尸?(t(?)对于输入层为0层，当？?=0,乘坐2号电梯的概率为0,当？?=27,乘坐2号电梯的概率为1/2,假设次概率服从线性关系，则乘坐 2号电梯的概率为那么乘坐1、2号电梯的乘客等待时间的期望为 ?W(?,?) =51w(?)+(1-51)

13、W(?) 5454? ?!(?2(t(?!)+?(t(?2)? ? (?(t(?)+?(t(?1)542(1-? 2?。(?勿)( -54)2(1-? 1?(t(?1)同时，记A为所有乘客到达的泊松强度，则乘 1、2号电梯乘客的泊松强度为？A 1、2号电梯“乘客集合”的泊松强度分别为？ ？1A?阳14号? ?1 A?=画出。为了解出模型，我们需要猊A和？三组参数。对于？?，我们实地做了实验，统计记录下了一组电梯停留时间的数据，如图所示:我们发现，数据大致都集中在一条平行于x轴的直线上，对数据求均值得？?=6.7s 。对于？，我们找到了一家与问题中商务楼规模类似的公司，调查得到开车上班的人 所占

14、比例为42.3%,这里认为？=42.3% , ?2=57.7%对于A,我们同样是在这家公司大厅实地做了统计，得到 30分钟内到达329人, 这里认为 A= 0.183。取？?=1 , 2 27,得至ij W(?,?)与？?的关系如图从图中可以看出，当？?=14时，W(?,?)最小，即(?？,?)=27 104时为最优方案。 ？0?(?2(t(?0)+?(t(?3)+(1-54)2(1-? 3?吟)同样，对于输入层为-1层，有八8 ._? ?2(?2。(？)+?。(?？4)(, ?, ?)=542(1-?4?(t(?%)且 t(?)=4N(?) +4+ n?b13J；温,?=(瑶).，？?=54

15、第 得到W(?,?)与?的关系如图ICO工口口从图中可以看出，当？?=14时，w(?,?)最小，即(?？,?尸；7；4时为最优方案。对于输入层为-2层，有WP? ?)= ? ?&?2(t(?)+?(t(?6)(? ?, ? ?)=542(1-?6?(t(?6)?® ?(?2(t(?)+?(t(?5) 共1百)2(1-? 5?(t(?0)得到W(?,?)与?的关系如图124时为最优方案。?13'且 t(?)=4N(?) +6+ n?从图中可以看出，当？?=14时,?=(仔爰，？符等，2w(?,?)最小,即(?,?尸227于是我们得到，当a=2710 2714 27 -1

16、24时，始)最小，为?2?2f(A尸？W(?,?) +W(?,?) + W(?,?) = 33.34。5.1.2下班高峰期单轿厢电梯系统的求解对于下班高峰期，每个目标层都要有一个区间从本层到27层的电梯以保证任何输入层的乘客都能到达目标层，则 ？=(0,27)?, ?=(-1,27) ?, ?=(-2,27) ?,同时令 ?=(0, ?)? ?=(-1, ?)?, ?=(-2, ?)?0对于每个输入层的乘客，都有刚好没乘上电梯的乘客需要等待电梯一次服务之后才由图可知,当？?=14时，w(?,?)最小,即(?,?尸27 -24时为最优方案。口,壬门7目不" 、.0I, 00-1 -1-

17、2-2 1 nPf/ 甲-彳寸至,A=A= 271327 142714盯，f(A)取，为f(A尸？W(?, ?) +-22W(?,?) + 甘W(?,?) = 45.06。5.1.3 非高峰期单轿厢电梯系统的求解非高峰期的输入层和目标层都是随机分散的，且人流量小，因此不应分析电梯的区 间安排，而应从电梯在无乘梯需求时自动停留的位置入手分析。如4.2所说，记3=(?,?,?,?,?,?)?,设某乘客所在楼层为n,则他所要等待的 时间为mint(?n)(i=1,2,3,4,5,6)。并且我们认为此乘客在-2层到27层的概率相等，故 等待时间的期望一,-0-71_f( 0=11?=-2 mint(?

18、n) ,(i=1,2,3,4,5,6)30通过编程枚举，可以得出，当月-2,2,7,12,17,22) ?W, f(份最小，为f( S=I|?=-2 Amint(?n) = 2.47。305.1.4 模型结论至此，我们得出了单轿厢电梯系统运行效率最优化的运行方案，即在高峰期采取方 案人二* 104 27、27 124,上班时乘客等待时间的期望为33.34s,下班时等待时间的期望为45.06s ;非高峰期采取方案3=(-2,2,7,12,17,22) ?，等待时间期望为2.4745.2双轿厢电梯系统的求解5.2.1 上班高峰期双轿厢电梯系统的求解对于上班高峰期，每个输入层都要有一个区间从本层到2

19、7层的电梯井道以保证乘客能到达任何目标层，和 5.1.1类似，令同一井道内两个电梯的区间相同，这样可以避 免控制台的混乱，贝U ?=(0,27)?, ?=(-1,27) ?, ?=(-2,27) ?同时令?=(0, ?)?, ?=(-1, ?)?, ?=(-2, 2球。此时，同一井道内两个电梯一次服务一共可以运载26个人，这里把26个乘客视为一个“乘客集合”，相应的泊松流强度为三，则此模型可以简化为排队模型。26同5.1.1,我们得到,W(? ?)= ? ?(?.(t(?)+?(t(?2) +? ?(?2(t(?)+?(t(?1)(? ?,? ?)=542(1-? 2?(t(?2)(1-54)

20、2(1-? 1 ?(t(?)?=(1 ?三?=?1A'? ( 54) 26 5 2 54 26故我们得到 W(?,?)与?的关系如图由图可知，当？?=12时，w(?,?)最小，即(?,？?尸127 1c21时为最优方案。同理可得 W(?,?)与？?、W(?,?)与？?的关系分别如图由图可知，当？?=14时，W(?,?)最小，即(?,？?)=27 一，时为最优方案。-24时为最优方案。由图可知,当？?=14时，W(?,?)最小,即(?,?尸27于是我们得到，当A=27 1c2 217-114 -27一 一 一?2一f(A)= ?W(?, ?) +22W(?,?) +力时，f(A)最小，为

21、?2_ _22W(?,?) = 12.24。5.2.2下班高峰期双轿厢电梯系统的求解对于下班高峰期，每个目标层都要有一个区间从本层到27层的电梯井道以保证任何输入层的乘客都能到达目标层，则？=(0,27)?,?=(0, ?)?, ?=(-1, ?)?, ?=(-2, ?茶?。同5.1.2,将26个乘客视为一个“乘客集合”， 的泊松流强度沿用5.1.2中的A'。?=(-1,27) ?, ?=(-2,27) ?同日寸令则此模型可简化为排队模型，参数W(?, ?)与?的关系分别如图SiO3口工。由图可知，当？?=12时，w(?,?)最小，即(?？,？?尸27 1c2时为最优方案。故我们得到

22、W(?,?)与?、W(?,?)与？?、由图可知，当？?=14时，W(?,?)最小，即(?,？?)=4 ；4时为最优方案。 - "* .30 至-一 20.XCo H111n11O5IOX520NG3-0由图可知，当？?=14时，w(?,?)最小,即(?,?尸成-124时为最优方案。T 口,壬门4目不" 、.0I, 00-1 -1-2-2 1 nP/八、甲.J、于我至'* 一271227 1427存打，f(A)取，为?f(A尸？W(?, ?) +-?2W(?,?) +即(？,？?)= 15.24。5.2.3非高峰期双轿厢电梯系统的求解非高峰期的输入层和目标层都是随机分

23、散的，且人流量小，因此同一井道中的电梯在无乘梯需求时自动停留的位置可以不同，则同 5.1.3,记 田(?,?-,?,？?2)-,设某 乘客所在楼层为n,则他所要等待的时间为 mint(?n)(i=1,2,312)。并且我们认为此 乘客在-2层到27层的概率相等，故等待时间的期望f( T7=-2 卷 mint(?n),(i=1,2,3,12)30通过编程枚举可以得出，当3=(-2,0,2,5,7,10,12,15,17,20,23,26) ?时，f(最小，为f( 9=2.?=-2 mint(?n) = 1.330305.2.4模型结论至此，我们得出了双轿厢电梯系统运行效率最优化的运行方案，即在高

24、峰期采取方案 A= 00-1 -1-2-2 卜班时乘客等彳寺时间的期望为12 24s 下班时等待l271427 142714' 口力、" ir口it目j/力/j乙.4rs, 口i、-时间的期望为15.24s；非高峰期采取方案 片(-2,0,2,5,7,10,12,15,17,20,23,26) ?,等待时 间期望为1.33s六模型的比较6.1高峰期电梯系统效率的比较上班高峰期，双轿厢电梯系统平均等待时间为12.24s,单轿厢电梯系统平均等待间为33.34s,双轿厢电梯系统比单轿厢系统效率提高了333X100%=172.4%班高峰期，双轿厢电梯系统平均等待时间为15.24s,单

25、轿厢电梯系统平均等待时间为45.06s,双轿厢电梯系统比单轿厢系统效率提高了5FX100%=195.7%。6.2非高峰期电梯系统效率的比较非高峰期，双轿厢电梯系统平均等待时间为2.47s,双轿厢电梯系统比单轿厢系统效率提高了1.33s,单轿厢电梯系统平均等待时间为2.47-1.33X 100%=85.7% 。1.33七 模型的灵敏度分析因为本文的模型所需参数几乎都是通过小范围的统计得到，因此还需考虑参数波动 对模型结果的影响。先考虑A的波动对结果的影响。这里将 A的值作± 0.05的波动，得到等待时间 望随楼层的变化，结果如图60W=0.577A=0.2 3W=a.577A=0 18

26、 3577D.l 331 2 5 4 5 6 7 S D 101112131 15-1617 IB25 26 27我们发现虽然期望值均随 A的波动而变化，但整体增减趋势没有改变。再考虑？的波动对结果的影响。这里将？的值作± 0.05的波动，得到等待时间期望 随楼层的变化，结果如图我们发现虽然期望值均随？的波动而变化，但整体增减趋势没有改变 所以我们认为模型结果是可信的。八模型的优缺点8.1 模型优点本模型最显著的优点就是简单直观，能很好地借助现有模型对问题进行分析和求 解，便于编程计算。同时将问题根据实际情况作不同考虑，建立不同的模型，使结果更 具实际参考意

27、义，而且提出的解决方案简单易行，在经济上几乎不会造成额外的支出， 可行性很强。8.2 模型缺点与本模型最显著优点一一简单相伴而来的缺点就是参数过多，对数据的依赖性强， 需要统计大量的真实数据才能更加准确地求解模型， 而由于时间有限，我们这里统计的 数据量还不够，参数的波动虽然对方案整体设计基本上没有影响，但对相关的数据结果可能会造成一些影响，还需要进一步加大数据统计量，以对模型作进一步完善。参考文献11 张莹，运筹学基础，北京：清华大学出版社， 201002 宋荣兴，孙海涛，运筹学，北京：经济科学出版社，2011。31 孟玉树，排队论基础及应用，上海：同济大学出版社，1989。附录1.求解期望

28、值的C语言程序#include <stdio.h>#include <stdlib.h>#include <math.h>#include <string.h># 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

29、*argv口)int i, j, p, q;p = 0; /* p取-2 -1 0分别计算*/for (q = 1; q < 28; q+)int TT = 0;int TN = 0;int data3030; /* data mJ 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+

30、)dataij= mcn(13, j) * zuhe(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 &l

31、t; 30; j+)if (dataij != 0)temp +=dataij * pow(2 * i * T1 + j * T2 - ave, 2);var = temp / TN;float E = 0.0;E = ave;float D = 0.0;D = var;float lmtl = 0.0;lmtl = (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) /

32、(1 - lmt2E) +(1 - (float) q / 27) * lmtl * 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；

33、elsefim=fi - 1m - 1 + m * fi - 1m;return fnm;intzuhe(int m, int n) jiecheng(mreturnjiecheng(m) / (jiecheng(n) * n)； intjiecheng(int m)int result = 1;int i;for (i = 1; i < m + 1; i+) result = result * i;return result;int min(int a, int b)if (a > b)return b;elsereturn a;2.求解枚举的C语言程序#include <stdio.h>#include <stdlib.h>void building(intg,inth,inti,int