【導讀】的提高對建筑物的有效利用和性能發(fā)揮將產生極為重要的影響。為提高服務質量和。行合理的調度和管理，即電梯群控。所謂電梯群控系統(tǒng)EGCS(ElevatorGroup. 物中的實際交通狀況，產生最優(yōu)派梯決策的控制系統(tǒng)[1]。而我國在電梯群控方面的起步比較晚，現(xiàn)階段對于電梯群控的關鍵技術尚未能。比還有較大的差距。因此，很有必要在電梯群控方面展開研究。歷史上第一臺真正的電梯出現(xiàn)在1889年12月，由美國Otis電梯公司研制，它是由電力驅動，齒輪直接傳動的。此時的電梯必須由司機操作運行，既浪費人力。為了改善這一問題，人們逐漸發(fā)展出了以下。下兩個方向的呼梯按鈕以供選擇，電梯能夠同時記住轎內選層和層站呼梯信號。廂應答啟動運行，在順向運動中，依次應答順向的呼梯，在呼梯層站?？?。的情況，對各臺電梯進行協(xié)調控制，采用最優(yōu)的輸送方式。交通模式由上行高峰模式、下行高峰模式及非高峰期模式等組成，采用。1975年至1982年是現(xiàn)代電梯群控系統(tǒng)的第一代，這一時期的電梯群控系統(tǒng)，

　　

【正文】 。 %維數 SwarmSize=400。 %粒子群規(guī)模 MaxIter=35。 %最大迭代次數 Vmax=3。 Vmin=3。 Ub=4*ones(1,Dim)。 %粒子群位置上限 清華 大學 2020 屆畢業(yè)論文 第 27 頁 共 43 頁 Lb=ones(1,Dim)。 %粒子群位置下限 wmax=。 %慣性權重 wmin=。 %粒子群初始化 Swarm=randint(SwarmSize,Dim,[1,NUMBER1])。 %初始化粒子群 VStep=randint(SwarmSize,Dim,[3,3])。 %初始化速度 fSwarm=zeros(SwarmSize,1)。 %初始化粒子群適應度值 for i=1:SwarmSize [averageWT,averageRT,averageE]=diantiqunkong_PSO_fitness2(Swarm(i,:),FO,DO,FC,DC,NUMBER2,TimeRun,TimeStop)。 AWT(i)=averageWT。 %第 i個粒子群的平均候梯時間 ART(i)=averageRT。 %第 i個粒子群的平均乘梯時間 AE(i)=averageE。 %第 i個粒子群的平均能量損耗 end %計算粒子的適應度值 for i=1:SwarmSize fSwarm(i,:)=weight1*(AWT(i)/max(AWT))+weight2*(ART(i)/max(ART))+weight3*(AE(i)/max(AE))。 end %個體極值和群體極值 [bestf,bestindex]=min(fSwarm)。 zbest=Swarm(bestindex,:)。 gbest=Swarm。 fgbest=fSwarm。 fzbest=bestf。 清華 大學 2020 屆畢業(yè)論文 第 28 頁 共 43 頁 [averageWT,averageRT,averageE]=diantiqunkong_PSO_fitness2(zbest,FO,DO,FC,DC,NUMBER2,TimeRun,TimeStop)。 y_fitness=fzbest。 AVERAGEwaittime=[averageWT]。 AVERAGEruntime=[averageRT]。 AVERAGEenergy=[averageE]。 %迭代尋優(yōu) iter=0。 while(iterMaxIter) for j=1:SwarmSize %速度更新 w=wmax(wmaxwmin)/MaxIter*iter。 VStep(j,:)=w*VStep(j,:)+c1*rand*(gbest(j,:)Swarm(j,:))+c2*rand*(zbestSwarm(j,:))。 if VStep(j,:)Vmax,VStep(j,:)=Vmax。end if VStep(j,:)Vmin,VStep(j,:)=Vmin。end %位置更新 Swarm(j,:)=Swarm(j,:)+VStep(j,:)。 for k=1:Dim if Swarm(j,k)Ub(k),Swarm(j,k)=Ub(k)。end if Swarm(j,k)Lb(k),Swarm(j,k)=Lb(k)。end end Swarm=round(Swarm)。 end for j=1:SwarmSize [averageWT,averageRT,averageE]=diantiqunkong_PSO_fitness2(Swarm(j,:),FO,DO,FC,DC,NUMBER2,TimeRun,TimeStop)。 AWT(j)=averageWT。 %第 i個粒子群的平均候梯時間 清華 大學 2020 屆畢業(yè)論文 第 29 頁 共 43 頁 ART(j)=averageRT。 %第 i個粒子群的平均乘梯時間 AE(j)=averageE。 %第 i個粒子群的平均能量損耗 end %計算粒子的適應度值 for j=1:SwarmSize fSwarm(j,:)=weight1*(AWT(j)/max(AWT))+weight2*(ART(j)/max(ART))+weight3*(AE(j)/max(AE))。 end %個體最優(yōu)更新 for j=1:SwarmSize if fSwarm(j)fgbest(j) gbest(j,:)=Swarm(j,:)。 fgbest(j)=fSwarm(j)。 end end %群體最優(yōu)更新 [bestf,bestindex]=min(fSwarm)。 if bestffzbest zbest=Swarm(bestindex,:)。 fzbest=bestf。 end [averageWT,averageRT,averageE]=diantiqunkong_PSO_fitness2(zbest,FO,DO,FC,DC,NUMBER2,TimeRun,TimeStop)。 iter=iter+1。 y_fitness=[y_fitness fzbest]。 AVERAGEwaittime=[AVERAGEwaittime,averageWT]。 AVERAGEruntime=[AVERAGEruntime,averageRT]。 AVERAGEenergy=[AVERAGEenergy,averageE]。 清華 大學 2020 屆畢業(yè)論文 第 30 頁 共 43 頁 end disp(39。最優(yōu)派梯方案 39。) zbest disp(39。最優(yōu)適應度值 39。) fzbest %繪圖輸出 figure(1) plot(y_fitness,39。LineWidth39。,2) title(39。最優(yōu)個體適應值 39。,39。fontsize39。,12)。 xlabel(39。迭代次數 39。,39。fontsize39。,12)。ylabel(39。適應值39。,39。fontsize39。,12)。 set(gca,39。fontsize39。,12)。 figure(2) plot(AVERAGEwaittime,39。*b39。) hold on。 plot(AVERAGEruntime,39。*g39。) plot(AVERAGEenergy,39。*r39。) title(39。AVERAGEwaittime161。162。AVERAGEruntime161。162。AVERAGEenergy優(yōu)化曲線 39。,39。fontsize39。,8)。 xlabel(39。迭代次數 39。,39。fontsize39。,8)。ylabel(39。參數值39。,39。fontsize39。,8)。 set(gca,39。fontsize39。,8)。 legend(39。AVERAGEwaittime39。,39。AVERAGEruntime39。,39。AVERAGEenergy39。,1)。 計算適應度程序 %計算電梯調度方案評價函數程序 function [averageWT,averageRT,averageE]=diantiqunkong_PSO_fitness2(X,FO,DO,FC,DC,NUMBER2,TimeRun,TimeStop) 清華 大學 2020 屆畢業(yè)論文 第 31 頁 共 43 頁 for i=1:NUMBER2 n=X(i)。 %響應第 i個呼梯信號的電梯序號 fo=FO(n)。 %n號電梯當前所在層 do=DO(n)。 %n號電梯當前運行方向 %第 n號電梯已登記內呼響應 if n==1 load DJNH1 DJNH=DJNH1。 elseif n==2 load DJNH2 DJNH=DJNH2。 elseif n==3 load DJNH3 DJNH=DJNH3。 else load DJNH4 DJNH=DJNH4。 end fc=FC(i)。%第 i個呼梯信號所在樓層 dc=DC(i)。%第 i個呼梯信號去往方向 %計算候梯時間 %計算候梯時間時的類屬情況 if dc==do if do==1 STATE=1。 else STATE=4。 end elseif do==1 STATE=2。 else 清華 大學 2020 屆畢業(yè)論文 第 32 頁 共 43 頁 STATE=3。 end %根據不同情況計算候梯時間 switch STATE case {1} %同向上行 if fc==fo WT(i)=0。 elseif fcfo if fcfo==1 m11=0。 else edge1=fo+1。 edge2=fc1。 Nstop=histc(DJNH,[edge1,edge2])。 m11=Nstop(1)+Nstop(2)。 %梯內指定停靠次數 end m12=randint(1,1,[1,fcfo1])。 %上升過程中隨機停靠次數 m12=round(m12/2)。 mw=m11+m12。 WT(i)=abs(fcfo)*TimeRun+mw*TimeStop。 else FMAX=16。 %上行可達最遠端 FMIN=1。 %下行可達最遠端 edge1=fo+1。 edge2=FMAX。 Nstop=histc(DJNH,[edge1,edge2])。 m11=Nstop(1)+Nstop(2)。 %梯內指定?？看螖? m12=randint(1,1,[1,(FMAXfo)m11+(15)+(fc2)])。%上升中隨機?？看螖? m12=round(m12/2)。 mw=m11+m12。 清華 大學 2020 屆畢業(yè)論文 第 33 頁 共 43 頁 WT(i)=(abs(FMAXfo)+abs(FMAXFMIN)+abs(fcFMIN))*TimeRun+mw*TimeStop。 end case {4} %同向下行 if fc==fo WT(i)=0。 elseif fcfo