JPH0790995B2 - 群管理エレベータの最適割当手法 - Google Patents
群管理エレベータの最適割当手法Info
- Publication number
- JPH0790995B2 JPH0790995B2 JP3124959A JP12495991A JPH0790995B2 JP H0790995 B2 JPH0790995 B2 JP H0790995B2 JP 3124959 A JP3124959 A JP 3124959A JP 12495991 A JP12495991 A JP 12495991A JP H0790995 B2 JPH0790995 B2 JP H0790995B2
- Authority
- JP
- Japan
- Prior art keywords
- solution
- call
- hall
- calls
- car
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Elevator Control (AREA)
Description
【0001】
【産業上の利用分野】この発明は群管理エレベータの呼
び割当てに係り、特に既知の呼び発生に対する割当ての
最適解を探索する手法に関するものである。
び割当てに係り、特に既知の呼び発生に対する割当ての
最適解を探索する手法に関するものである。
【0002】
【従来の技術】現在の群管理エレベータは、乗場呼びを
選択された特定のかごにだけ与える呼び割当て制御方式
が主流であるが、この呼び割当て制御方式として、従来
から評価関数による呼び割当て制御に加え、最近では、
例えば特開平1−197287号「エレベータの群管理
制御方法」のように、ファジー理論を用いた呼び割当て
制御や、或いは特開平1−275381号「エレベータ
の群管理装置」のように、ニューラルネットを用いた呼
び割当て制御など種々の方式が提案されている。
選択された特定のかごにだけ与える呼び割当て制御方式
が主流であるが、この呼び割当て制御方式として、従来
から評価関数による呼び割当て制御に加え、最近では、
例えば特開平1−197287号「エレベータの群管理
制御方法」のように、ファジー理論を用いた呼び割当て
制御や、或いは特開平1−275381号「エレベータ
の群管理装置」のように、ニューラルネットを用いた呼
び割当て制御など種々の方式が提案されている。
【0003】一方、既知の呼び発生に対する割当ての最
適解を求め、すなわち乗客の呼びの状況(すべての乗場
呼びの発生時刻,発生階床,行先階床等)が既知である
場合に、複数台のかごに対する呼び割当ての最適解(各
かごの最適運行)を求め、これを上記のようなファジー
理論やニューラルネットを用いた呼び割当て制御に活用
することが行われている。
適解を求め、すなわち乗客の呼びの状況(すべての乗場
呼びの発生時刻,発生階床,行先階床等)が既知である
場合に、複数台のかごに対する呼び割当ての最適解(各
かごの最適運行)を求め、これを上記のようなファジー
理論やニューラルネットを用いた呼び割当て制御に活用
することが行われている。
【0004】勿論、実際の群管理エレベータにおいて
は、呼びを割当てる際には将来の呼びの発生が未確実な
ため、上記のような既知の呼び発生に対する最適解をそ
のまま呼び割当ての制御アルゴリズムに適用することは
できないが、例えばファジー理論を用いた呼び割当ての
ルール抽出に活用したり、或いはニューラルネットを用
いた呼び割当てにおいて、群管理に有効な制御法則を学
習する際の教師信号として使用することができ、その他
にも制御アルゴリズムを考える上で大いに参考になるデ
ータを得ることができるので種々の活用が考えられる。
は、呼びを割当てる際には将来の呼びの発生が未確実な
ため、上記のような既知の呼び発生に対する最適解をそ
のまま呼び割当ての制御アルゴリズムに適用することは
できないが、例えばファジー理論を用いた呼び割当ての
ルール抽出に活用したり、或いはニューラルネットを用
いた呼び割当てにおいて、群管理に有効な制御法則を学
習する際の教師信号として使用することができ、その他
にも制御アルゴリズムを考える上で大いに参考になるデ
ータを得ることができるので種々の活用が考えられる。
【0005】ところで、すべての呼びの発生状況が既知
であり、それらの呼びに複数台のかごを割当てる場合、
すべての割当ての組合せの中から最適解を見つけ出すと
いう問題は、巡回セールスマン問題などでよく知られて
いる組合せ最適化問題としてとらえることができる。こ
の組合せ最適化問題の解法には種々あるが、近似解法の
1つとしてよく知られているものにシミュレーテッド・
アニーリング法(Simulated Annealing 法,以下SA法
と略す)があり、このSA法を用いて群管理エレベータ
における呼び割当ての最適解(最適運行)を探索する方
法が、例えば特開昭62−205972号「エレベータ
群最適運行解析システム」に示されている。
であり、それらの呼びに複数台のかごを割当てる場合、
すべての割当ての組合せの中から最適解を見つけ出すと
いう問題は、巡回セールスマン問題などでよく知られて
いる組合せ最適化問題としてとらえることができる。こ
の組合せ最適化問題の解法には種々あるが、近似解法の
1つとしてよく知られているものにシミュレーテッド・
アニーリング法(Simulated Annealing 法,以下SA法
と略す)があり、このSA法を用いて群管理エレベータ
における呼び割当ての最適解(最適運行)を探索する方
法が、例えば特開昭62−205972号「エレベータ
群最適運行解析システム」に示されている。
【0006】ここでSA法とは、例えば1986年1月
発行の「オペレーションズ・リサーチ第31巻第1号4
3頁〜48頁」に記載されているように周知の手法であ
るが、簡単に説明するとその手順は次に示すようにな
る。 初期解Xを生成する。 初期温度Tを設定する。 解の変換(摂動とも呼ばれる)を行う。 解の変換による解の値の変化量Δを計算する。 Δ>0すなわち解の値が改善された場合は解の変換を
受理する。 Δ<0すなわち解の値が改善されなかった場合は P(Δ)=exp(−Δ/T)の確率で解の変換を受理
する。 解の値の分布が定常状態になるまで上記〜を繰り
返す。 温度Tを少し低い温度に再設定する。 上記〜を繰り返し、解の値の変動が十分に小さく
なると停止する。 このようにして、温度Tをうまく制御していくと最適解
または最適解に近い最良解を得ることができるというも
のである。
発行の「オペレーションズ・リサーチ第31巻第1号4
3頁〜48頁」に記載されているように周知の手法であ
るが、簡単に説明するとその手順は次に示すようにな
る。 初期解Xを生成する。 初期温度Tを設定する。 解の変換(摂動とも呼ばれる)を行う。 解の変換による解の値の変化量Δを計算する。 Δ>0すなわち解の値が改善された場合は解の変換を
受理する。 Δ<0すなわち解の値が改善されなかった場合は P(Δ)=exp(−Δ/T)の確率で解の変換を受理
する。 解の値の分布が定常状態になるまで上記〜を繰り
返す。 温度Tを少し低い温度に再設定する。 上記〜を繰り返し、解の値の変動が十分に小さく
なると停止する。 このようにして、温度Tをうまく制御していくと最適解
または最適解に近い最良解を得ることができるというも
のである。
【0007】
【発明が解決しようとする課題】ところで、このSA法
を実行するに当っては上記手順すなわちAnnealing のス
ケジュールとして、初期温度、温度の系列、各温度にお
ける定常状態の判定条件、停止の判定条件等を決めなけ
ればならないが、このスケジュールの決め方が非常に困
難であるという問題があり、また常に確率の計算を伴う
ので計算に非常に時間がかかるといった問題点があっ
た。
を実行するに当っては上記手順すなわちAnnealing のス
ケジュールとして、初期温度、温度の系列、各温度にお
ける定常状態の判定条件、停止の判定条件等を決めなけ
ればならないが、このスケジュールの決め方が非常に困
難であるという問題があり、また常に確率の計算を伴う
ので計算に非常に時間がかかるといった問題点があっ
た。
【0008】このため、最近SA法を改良するものとし
て、例えば1990年発行の「JOURNAL OF COMPUTATION
AL PHYSICS 90,P161-175」にThreshold Accepting 法
(以下TA法と略す)なるものが提案されており、その
手順は次のようになる。 初期解Xを生成する。 しきい値Tを設定する。 解の変換を行う。 解の変換による値の変化量Δを計算する。 Δ>−Tのとき解の変換を受理し、そうでなければ受
理しない。すなわち解の値が改善されたか、若しくは改
悪の場合でもその変化量がしきい値Tより小さい場合は
解の変換を受理する。 解の値の分布が定常状態になるまで上記〜を繰り
返す。 しきい値Tを少し小さい値に再設定する。 上記〜を繰り返し、解の値の変動が十分に小さく
なると停止する。
て、例えば1990年発行の「JOURNAL OF COMPUTATION
AL PHYSICS 90,P161-175」にThreshold Accepting 法
(以下TA法と略す)なるものが提案されており、その
手順は次のようになる。 初期解Xを生成する。 しきい値Tを設定する。 解の変換を行う。 解の変換による値の変化量Δを計算する。 Δ>−Tのとき解の変換を受理し、そうでなければ受
理しない。すなわち解の値が改善されたか、若しくは改
悪の場合でもその変化量がしきい値Tより小さい場合は
解の変換を受理する。 解の値の分布が定常状態になるまで上記〜を繰り
返す。 しきい値Tを少し小さい値に再設定する。 上記〜を繰り返し、解の値の変動が十分に小さく
なると停止する。
【0009】このようにTA法は、解の値が改悪された
場合にSA法ではある確率で解の変換を受理するのに対
し、しきい値より小さい範囲内であればすべて受理する
ようにした点でSA法と異なる。この結果、TA法はス
ケジュールの決定において、また確率の計算が不要とな
る点でSA法に比べてかなり単純化され、しかもSA法
より良い結果の得られることが確かめられている。
場合にSA法ではある確率で解の変換を受理するのに対
し、しきい値より小さい範囲内であればすべて受理する
ようにした点でSA法と異なる。この結果、TA法はス
ケジュールの決定において、また確率の計算が不要とな
る点でSA法に比べてかなり単純化され、しかもSA法
より良い結果の得られることが確かめられている。
【0010】しかしながらエレベータの場合、かごがm
台、呼びの発生数をn個とすると、割当ての組合せすな
わち解の数はmn となるので、例えば或るビルの1日の
呼びを対象にその1日の最適運行を求めようとすると大
規模な組合せ最適化問題となり、例えTA法を用いたと
してもあまりにも時間がかかりすぎるという問題点があ
った。
台、呼びの発生数をn個とすると、割当ての組合せすな
わち解の数はmn となるので、例えば或るビルの1日の
呼びを対象にその1日の最適運行を求めようとすると大
規模な組合せ最適化問題となり、例えTA法を用いたと
してもあまりにも時間がかかりすぎるという問題点があ
った。
【0011】本発明は上記の点に鑑みなされたもので、
群管理エレベータにおける呼び割当てを組合せ最適化問
題として捉え、TA法を適用してより質の高い最適解を
得るとともに、エレベータに特有の性質を利用すること
により演算回数を減らし、収束の早い探索方法を提供す
るものである。
群管理エレベータにおける呼び割当てを組合せ最適化問
題として捉え、TA法を適用してより質の高い最適解を
得るとともに、エレベータに特有の性質を利用すること
により演算回数を減らし、収束の早い探索方法を提供す
るものである。
【0012】
【課題を解決するための手段】本発明は上記の問題点を
解決するため、複数の既知の乗場呼びに複数台のかごを
割当てる場合の最適解を探索するに際し、各かごに数が
ほぼ均等となるように前記乗場呼びを任意に割当ててこ
れを初期解とし、解の変換は前記乗場呼びのうち異なる
かごに割当てられた任意の2つの乗場呼びを入れ替える
ことにより行い、これを探索の基本ステップとしてThre
shold Accepting 法を適用することにより、最適解を探
索するようにしたことを特徴とするものである。
解決するため、複数の既知の乗場呼びに複数台のかごを
割当てる場合の最適解を探索するに際し、各かごに数が
ほぼ均等となるように前記乗場呼びを任意に割当ててこ
れを初期解とし、解の変換は前記乗場呼びのうち異なる
かごに割当てられた任意の2つの乗場呼びを入れ替える
ことにより行い、これを探索の基本ステップとしてThre
shold Accepting 法を適用することにより、最適解を探
索するようにしたことを特徴とするものである。
【0013】
【作用】かごの台数をm,乗場呼びの数をnとすると、
その割当ての解の数はmn 通りとなるが、本発明によれ
ば、まず各かごに数が均等となるように乗場呼びを割当
て初期解を生成し、その後の変換は前記乗場呼びのうち
異なるかごに割当てられた2つの乗場呼びを入れ替える
ことにより行うようにしたので、探索の対象となる解の
数はmn から約n!/((n/m)!)m 通りへと大幅
に減少し、収束が非常に速くなる。しかも各かごに割当
てられた呼びの数が均等となる解の中には、エレベータ
の性質上、解の値が極端に悪いものは存在せず、しかも
最適解が含まれている可能性が高いというだけでなく、
解の変換はこの呼びを入れ替えるだけなので解の変換に
よって解の値が極端に悪化したり大きく変動したりする
ことがなくなり、探索空間はTA法に適した非常に滑ら
かなものとなる。
その割当ての解の数はmn 通りとなるが、本発明によれ
ば、まず各かごに数が均等となるように乗場呼びを割当
て初期解を生成し、その後の変換は前記乗場呼びのうち
異なるかごに割当てられた2つの乗場呼びを入れ替える
ことにより行うようにしたので、探索の対象となる解の
数はmn から約n!/((n/m)!)m 通りへと大幅
に減少し、収束が非常に速くなる。しかも各かごに割当
てられた呼びの数が均等となる解の中には、エレベータ
の性質上、解の値が極端に悪いものは存在せず、しかも
最適解が含まれている可能性が高いというだけでなく、
解の変換はこの呼びを入れ替えるだけなので解の変換に
よって解の値が極端に悪化したり大きく変動したりする
ことがなくなり、探索空間はTA法に適した非常に滑ら
かなものとなる。
【0014】
【実施例】以下、本発明の一実施例を図面に基いて説明
する。図1は、一例として7階床のビルにおける乗場呼
びの発生状況(乗客の到着状況)を時間を横軸にとって
示したもので、これらの既知の呼びに対する割当てを組
合せ最適化問題として考える。
する。図1は、一例として7階床のビルにおける乗場呼
びの発生状況(乗客の到着状況)を時間を横軸にとって
示したもので、これらの既知の呼びに対する割当てを組
合せ最適化問題として考える。
【0015】図1において、H1 は1階の乗場に到着し
た乗客によって発生した上昇方向の乗場呼び、H2 は同
じく4階の乗場に到着した乗客によって発生した下降方
向の乗場呼び、H3 〜H9 は同様に図示した階床位置で
発生した乗場呼びを示している。また、B1 〜B3 は乗
場呼びH1 が発生後に1階の乗場に遅れて到着した乗客
をそれぞれ示し、他の〇印も同様に各乗場に図示の時点
で到着した乗客を表わしている。
た乗客によって発生した上昇方向の乗場呼び、H2 は同
じく4階の乗場に到着した乗客によって発生した下降方
向の乗場呼び、H3 〜H9 は同様に図示した階床位置で
発生した乗場呼びを示している。また、B1 〜B3 は乗
場呼びH1 が発生後に1階の乗場に遅れて到着した乗客
をそれぞれ示し、他の〇印も同様に各乗場に図示の時点
で到着した乗客を表わしている。
【0016】図1の例では、例えば1階では乗場呼びH
1 に割当てられた或るかご(図示せず)に乗客B1 〜B
3 も同時に乗り込んで出発し、その後に到着した乗客に
よって乗場呼びH5 が発生する様子を示しているが、も
し各かごの初期位置が異なったり乗場呼びH1 に他のか
ごが割当てられた場合にはその応答時点が異なるため、
例えば乗客B1 とB2 を乗せて出発し、後から到着した
乗客B3 によって1階の乗場呼びが発生するというよう
に、各乗客の発生時点が同じでも割当てや各かごの初期
位置によって乗場呼びの発生状況も図1とは異なる状況
となる可能性がある。しかし、各乗客の発生時点と行先
階床及び各かごの初期位置はすべて既知であり不変とし
て問題を取り扱うため、各かごの呼び割当てを発生順に
順次決定していけば、各かごの運行が自動的に定まるこ
とになり、後述のシミュレーションを容易に実行するこ
とができる。
1 に割当てられた或るかご(図示せず)に乗客B1 〜B
3 も同時に乗り込んで出発し、その後に到着した乗客に
よって乗場呼びH5 が発生する様子を示しているが、も
し各かごの初期位置が異なったり乗場呼びH1 に他のか
ごが割当てられた場合にはその応答時点が異なるため、
例えば乗客B1 とB2 を乗せて出発し、後から到着した
乗客B3 によって1階の乗場呼びが発生するというよう
に、各乗客の発生時点が同じでも割当てや各かごの初期
位置によって乗場呼びの発生状況も図1とは異なる状況
となる可能性がある。しかし、各乗客の発生時点と行先
階床及び各かごの初期位置はすべて既知であり不変とし
て問題を取り扱うため、各かごの呼び割当てを発生順に
順次決定していけば、各かごの運行が自動的に定まるこ
とになり、後述のシミュレーションを容易に実行するこ
とができる。
【0017】図2は、本発明による最適解の探索手法を
説明するためのフローチャートであり、各ステップの内
容を図1の例に適用した場合について説明する。まず、
ステップS1において、次のように割当ての初期化を行
う。 Ai =((i−1)mod m)+1 (i=1・・・ n) ここでAi は乗場呼びHi に割当てられるかごの号機番
号を表わし、nは呼びの個数、mはかご台数である。な
お、((i−1)mod m)は(i−1)をmで割った余
りを示す。つまりこの式はn個の乗場呼びをm台のかご
に順番に割当てていくことを示している。
説明するためのフローチャートであり、各ステップの内
容を図1の例に適用した場合について説明する。まず、
ステップS1において、次のように割当ての初期化を行
う。 Ai =((i−1)mod m)+1 (i=1・・・ n) ここでAi は乗場呼びHi に割当てられるかごの号機番
号を表わし、nは呼びの個数、mはかご台数である。な
お、((i−1)mod m)は(i−1)をmで割った余
りを示す。つまりこの式はn個の乗場呼びをm台のかご
に順番に割当てていくことを示している。
【0018】図1の例ではn=9であり、またかごの台
数m=3とすると、上記の式により、各乗場呼びHi に
割当てられるかごの初期解は、図3(a) に示すように、
各乗場呼びHi に対して1号機〜3号機が順番に繰り返
し割当てられることになる。この結果、1号機〜3号機
に割当てられた呼びをC1 〜C3 としてそれぞれ号機別
に示すと、図3(b) のように各号機はそれぞれ3つずつ
乗場呼びを割当てられたことになり、これを初期解とす
る。
数m=3とすると、上記の式により、各乗場呼びHi に
割当てられるかごの初期解は、図3(a) に示すように、
各乗場呼びHi に対して1号機〜3号機が順番に繰り返
し割当てられることになる。この結果、1号機〜3号機
に割当てられた呼びをC1 〜C3 としてそれぞれ号機別
に示すと、図3(b) のように各号機はそれぞれ3つずつ
乗場呼びを割当てられたことになり、これを初期解とす
る。
【0019】なお、前述のように初期のかご位置や割当
てられたかごが異なる場合には呼びの発生状況が異なる
可能性があるため、図1のH1 〜H9 と図3(b) のH1
〜H9 とは必ずしも一致するものではなく、ここではH
1 〜H9 はあくまでも乗場呼びの発生した順序を表わ
し、図3(b) はC1 すなわち1号機に割当てられる呼び
は1番目と4番目と7番目に発生する乗場呼びであるこ
とを表わしている。
てられたかごが異なる場合には呼びの発生状況が異なる
可能性があるため、図1のH1 〜H9 と図3(b) のH1
〜H9 とは必ずしも一致するものではなく、ここではH
1 〜H9 はあくまでも乗場呼びの発生した順序を表わ
し、図3(b) はC1 すなわち1号機に割当てられる呼び
は1番目と4番目と7番目に発生する乗場呼びであるこ
とを表わしている。
【0020】次に、ステップS2ではTA法で重要な役
割を果たすしきい値Tの初期値をTo として初期化を行
う。
割を果たすしきい値Tの初期値をTo として初期化を行
う。
【0021】次に、ステップS3では上記の初期解に対
してシミュレーションを実行し、最適化の目標関数の値
を計算する。ここで最適化の目標関数としては種々考え
られるが、ここでは一般的な各乗場呼びに対する待時間
の平均値ta を最適化の目標関数として考える。
してシミュレーションを実行し、最適化の目標関数の値
を計算する。ここで最適化の目標関数としては種々考え
られるが、ここでは一般的な各乗場呼びに対する待時間
の平均値ta を最適化の目標関数として考える。
【0022】シミュレーションによる各乗場呼びの待時
間の計算は次のようにして行う。まず図3(b) に示すよ
うに、1番最初に発生した乗場呼びH1 は1号機に割当
てられるが、図1では1階の上昇方向の乗場呼びH1 が
これに相当し、その発生時点と1号機のかごの初期位置
とから乗場呼びH1 に対する待時間が容易に求まる。ま
た、これにより1号機の応答時点が確定するので、図1
の1階の各乗客の発生時点とから、1号機の応答時点ま
でに到着している乗客はそのかごに乗り、その直後に到
着した乗客によって次に発生する1階の乗場呼びの発生
時点が確定する。
間の計算は次のようにして行う。まず図3(b) に示すよ
うに、1番最初に発生した乗場呼びH1 は1号機に割当
てられるが、図1では1階の上昇方向の乗場呼びH1 が
これに相当し、その発生時点と1号機のかごの初期位置
とから乗場呼びH1 に対する待時間が容易に求まる。ま
た、これにより1号機の応答時点が確定するので、図1
の1階の各乗客の発生時点とから、1号機の応答時点ま
でに到着している乗客はそのかごに乗り、その直後に到
着した乗客によって次に発生する1階の乗場呼びの発生
時点が確定する。
【0023】この時点と図1の4階下降乗場呼びH2 と
の発生時点を比較し、もし4階の乗場呼びの方が早けれ
ばこれが図3(b) に示すH2 となるので2号機に割当て
られ、この乗場呼びH2 と2号機のかご位置とから乗場
H2 に対する待時間が求まり、2号機の応答時点も確定
する。
の発生時点を比較し、もし4階の乗場呼びの方が早けれ
ばこれが図3(b) に示すH2 となるので2号機に割当て
られ、この乗場呼びH2 と2号機のかご位置とから乗場
H2 に対する待時間が求まり、2号機の応答時点も確定
する。
【0024】このようにして図3(b) に従い、発生した
乗場呼びを順番に1〜3号機に割当てていくと、各乗客
の発生時点は既知なので各かごの運行と呼びの発生が確
定していき、各乗場呼びに対する待時間をそれぞれ求め
ることができる、なお、図1において、割当ての結果、
乗場呼びの発生数が10個以上となる場合もあるが同様
にして10番目以降を各号機に順番に割当てていけばよ
い。
乗場呼びを順番に1〜3号機に割当てていくと、各乗客
の発生時点は既知なので各かごの運行と呼びの発生が確
定していき、各乗場呼びに対する待時間をそれぞれ求め
ることができる、なお、図1において、割当ての結果、
乗場呼びの発生数が10個以上となる場合もあるが同様
にして10番目以降を各号機に順番に割当てていけばよ
い。
【0025】こうして初期解における各乗場呼びに対す
る待時間が求まると、その平均値ta を計算し、ステッ
プS4においてこれを現時点における最適値t* として
保存する。これで初期化が完了し、ステップS5以下で
最適解の探索を開始する。
る待時間が求まると、その平均値ta を計算し、ステッ
プS4においてこれを現時点における最適値t* として
保存する。これで初期化が完了し、ステップS5以下で
最適解の探索を開始する。
【0026】ステップS5では解の変換を行う。ここで
解の変換はすべての解を対象とするのではなく、直前の
解に対して異なるかごに割当てられた任意の2つの乗場
呼びだけを入れ替えることにより行う。すなわち次の条
件により入れ替えるべき乗場呼びを選択し、これを入れ
替えて新たな解とする。 P(j=i)=1/n (i=1・・・ n) P(k=i)=1/n (i=1・・・ n) ただし j≠k , Aj ≠Ak
解の変換はすべての解を対象とするのではなく、直前の
解に対して異なるかごに割当てられた任意の2つの乗場
呼びだけを入れ替えることにより行う。すなわち次の条
件により入れ替えるべき乗場呼びを選択し、これを入れ
替えて新たな解とする。 P(j=i)=1/n (i=1・・・ n) P(k=i)=1/n (i=1・・・ n) ただし j≠k , Aj ≠Ak
【0027】ここで、P(j=i)或いはP(k=i)
は、乗場呼びHi のうち、入れ替えの対象としてHj 或
いはHk (j,k=1・・・ n)が選択される確率を表わ
す。いま、入れ替え前のHj とHk に割当てられたかご
をそれぞれAj ,Ak とし、入れ替え後のHj とHk に
割当てられたかごをそれぞれA′j ,A′k で表わすと
すると、Hj とHk の割当てを入れ替えるとA′j =A
k ,A′k =Aj となる。
は、乗場呼びHi のうち、入れ替えの対象としてHj 或
いはHk (j,k=1・・・ n)が選択される確率を表わ
す。いま、入れ替え前のHj とHk に割当てられたかご
をそれぞれAj ,Ak とし、入れ替え後のHj とHk に
割当てられたかごをそれぞれA′j ,A′k で表わすと
すると、Hj とHk の割当てを入れ替えるとA′j =A
k ,A′k =Aj となる。
【0028】例えば、いまj,kとしてそれぞれj=
8,j=9が選択され、すなわち図3(b) の割当てに対
し乗場呼びH8 とH9 の割当てが入れ替えられたとする
と、A′9 =A8 ,A′8 =A9 となり、H9 は2号機
にH8 は3号機に割当てられることになるので、変換後
の新たな解は図4(a),(b) に示すようになる。
8,j=9が選択され、すなわち図3(b) の割当てに対
し乗場呼びH8 とH9 の割当てが入れ替えられたとする
と、A′9 =A8 ,A′8 =A9 となり、H9 は2号機
にH8 は3号機に割当てられることになるので、変換後
の新たな解は図4(a),(b) に示すようになる。
【0029】ステップS6ではこの新たな解に対して前
述と同様にシミュレーションを実行し、各乗場呼びに対
する待時間の平均値ta を算出する。そしてステップS
7では新たな解における待時間の平均値ta と現時点に
おける最適値t* との差をt* で除して0〜1の範囲に
正規化し、これがしきい値Tより小さいか否かを判断す
る。すなわちステップS7では、解を変換したことによ
り、解の値が改善されたか(待時間の平均値ta が現時
点の最適値t* より小さくなったか)或いは改悪された
としてもその改悪量が所定値より小さいか否かを判断す
る。
述と同様にシミュレーションを実行し、各乗場呼びに対
する待時間の平均値ta を算出する。そしてステップS
7では新たな解における待時間の平均値ta と現時点に
おける最適値t* との差をt* で除して0〜1の範囲に
正規化し、これがしきい値Tより小さいか否かを判断す
る。すなわちステップS7では、解を変換したことによ
り、解の値が改善されたか(待時間の平均値ta が現時
点の最適値t* より小さくなったか)或いは改悪された
としてもその改悪量が所定値より小さいか否かを判断す
る。
【0030】この結果、解の値が改善されたか或いは所
定値以内の改悪の場合にはステップS8へと進んでその
解における待時間の平均値ta を最適値t* とし、一
方、解の値が所定値以上の改悪となった場合にはステッ
プS9へと進んで解の変換を元に戻す。
定値以内の改悪の場合にはステップS8へと進んでその
解における待時間の平均値ta を最適値t* とし、一
方、解の値が所定値以上の改悪となった場合にはステッ
プS9へと進んで解の変換を元に戻す。
【0031】ステップS10では探索回数が所定値に達
したか否かを判断し、達していなければステップS11
でしきい値Tの値を一定量δTで僅かに小さくして再設
定し、その後解の変換を再び繰り返す。こうして探索を
繰り返し、所定回数に達すると探索を終了し、その時点
におけるt* の値が待時間の平均値の最適値であり、そ
の最適値t* に対する解が呼び割当ての最適解となる。
したか否かを判断し、達していなければステップS11
でしきい値Tの値を一定量δTで僅かに小さくして再設
定し、その後解の変換を再び繰り返す。こうして探索を
繰り返し、所定回数に達すると探索を終了し、その時点
におけるt* の値が待時間の平均値の最適値であり、そ
の最適値t* に対する解が呼び割当ての最適解となる。
【0032】なお、以上の説明において階床数やかごの
台数或いは呼び(乗客)の発生数は上記実施例に限定さ
れないことは言うまでもないが、その他、最適化の目標
関数としても、上記の待時間の平均値に限らず、例えば
最大待時間であるとか、消費エネルギーや最大輸送量、
或いはそれらの複合(論理和、重み付き合計)等、種々
の関数を用いることができる。
台数或いは呼び(乗客)の発生数は上記実施例に限定さ
れないことは言うまでもないが、その他、最適化の目標
関数としても、上記の待時間の平均値に限らず、例えば
最大待時間であるとか、消費エネルギーや最大輸送量、
或いはそれらの複合(論理和、重み付き合計)等、種々
の関数を用いることができる。
【0033】また、上記のステップS5の説明におい
て、乗場及びHj 或いはHk を選択する確率Pは1/n
と均等であったが、例えば探索の初期段階ではP(j=
1)>P(j=n)すなわち発生順序の早い乗場呼びの
選択確率を大きくし、後半の段階ではP(j=1)<P
(j=n)すなわち発生順序の遅い乗場呼びの選択確率
が大きくなるように不均等にすることもできる。発生順
序の早い乗場呼びほど入れ替えにより後の呼びの発生に
与える影響が大きいので、このように探索の初期の段階
ほど発生順序の早い乗場呼びが選択されるようにする
と、収束性の向上を図ることができる。
て、乗場及びHj 或いはHk を選択する確率Pは1/n
と均等であったが、例えば探索の初期段階ではP(j=
1)>P(j=n)すなわち発生順序の早い乗場呼びの
選択確率を大きくし、後半の段階ではP(j=1)<P
(j=n)すなわち発生順序の遅い乗場呼びの選択確率
が大きくなるように不均等にすることもできる。発生順
序の早い乗場呼びほど入れ替えにより後の呼びの発生に
与える影響が大きいので、このように探索の初期の段階
ほど発生順序の早い乗場呼びが選択されるようにする
と、収束性の向上を図ることができる。
【0034】また上記の実施例では、乗場呼びを各かご
に均等に割当て、解の変換も任意の2つの呼びを入れ替
えるだけなので、常に各かごに割当てられた呼びの数が
均等になるようにしている。これは充分長い時間のエレ
ベータ動作及び乗客発生は、それぞれ平均到着率μ、平
均発生率λのポアソンプロセスで近似することができ、
その時の各かごの利用率ρi は、各かごの分担する乗客
発生率をλi とすると ρi =λi /μ (i=1・・・ m,Σλi =λ) となる。そして全体の平均待時間tt は tt =Σ(1/mμ)・(ρi /(1−ρi )) であり、その最小値はρ1 =ρ2 =・・・・=ρm すなわ
ち、各かごの呼び割当てが均等となる場合である(この
ことは容易に証明できる)。
に均等に割当て、解の変換も任意の2つの呼びを入れ替
えるだけなので、常に各かごに割当てられた呼びの数が
均等になるようにしている。これは充分長い時間のエレ
ベータ動作及び乗客発生は、それぞれ平均到着率μ、平
均発生率λのポアソンプロセスで近似することができ、
その時の各かごの利用率ρi は、各かごの分担する乗客
発生率をλi とすると ρi =λi /μ (i=1・・・ m,Σλi =λ) となる。そして全体の平均待時間tt は tt =Σ(1/mμ)・(ρi /(1−ρi )) であり、その最小値はρ1 =ρ2 =・・・・=ρm すなわ
ち、各かごの呼び割当てが均等となる場合である(この
ことは容易に証明できる)。
【0035】また上記の式によらずとも、一般的に考え
て充分長い時間を考えればその間の呼びに対する最適割
当ての結果は、各かごへの割当てが均等となることは容
易に予測できる。このため呼びの数が常に均等となる場
合だけを対象とし探索を行うようにしている。しかし、
最適解が常に呼びを完全に均等に割り当てた場合に存在
するとは限らず、また割当てが多少不均等であっても最
適解と同等の解が含まれる可能性も高いので、割当てに
多少の不均等を含むようにしてもよい。
て充分長い時間を考えればその間の呼びに対する最適割
当ての結果は、各かごへの割当てが均等となることは容
易に予測できる。このため呼びの数が常に均等となる場
合だけを対象とし探索を行うようにしている。しかし、
最適解が常に呼びを完全に均等に割り当てた場合に存在
するとは限らず、また割当てが多少不均等であっても最
適解と同等の解が含まれる可能性も高いので、割当てに
多少の不均等を含むようにしてもよい。
【0036】そのためには、例えばある確率でA′j =
A k ,A′ k=A k となるように、すなわち、乗場
呼びHj の1つだけを他の号機に割当て変更するような
解の変換を許すことにより達成することができる。
A k ,A′ k=A k となるように、すなわち、乗場
呼びHj の1つだけを他の号機に割当て変更するような
解の変換を許すことにより達成することができる。
【0037】また、上記実施例の手順により得られた最
適解を1つの解とし、更に同一の手順をそれぞれ異なっ
た乱数を用いてそれぞれにつき最適解を得、これらの複
数の最適解の中から最良のものを最終的な最適解として
選択するようにすることもできる。
適解を1つの解とし、更に同一の手順をそれぞれ異なっ
た乱数を用いてそれぞれにつき最適解を得、これらの複
数の最適解の中から最良のものを最終的な最適解として
選択するようにすることもできる。
【00038】
【発明の効果】本発明によれば、群管理エレベータに特
有の性質を利用し、均等な呼び割当て又はその近傍だけ
を探索するようにしたので、探索空間は全ての解を対象
とした探索空間に比べて非常に小さくなり、極めて効率
よく最適な呼び割当てを探索することができる。
有の性質を利用し、均等な呼び割当て又はその近傍だけ
を探索するようにしたので、探索空間は全ての解を対象
とした探索空間に比べて非常に小さくなり、極めて効率
よく最適な呼び割当てを探索することができる。
【0039】また探索空間の解の質は比較的に良好であ
り、しかも、解の変換は任意の2つの乗場呼びを入れ替
えるだけか或いは1つの乗場呼びの割当変更だけにより
行うようにしているので、解の変換によって解の値が極
端に悪化したり大きく変動する可能性が小さく、滑らか
な探索空間となるため、局所解を乗り越える可能性が高
くなり、より最適解に到達する確率が高くなる。
り、しかも、解の変換は任意の2つの乗場呼びを入れ替
えるだけか或いは1つの乗場呼びの割当変更だけにより
行うようにしているので、解の変換によって解の値が極
端に悪化したり大きく変動する可能性が小さく、滑らか
な探索空間となるため、局所解を乗り越える可能性が高
くなり、より最適解に到達する確率が高くなる。
【0040】すなわち、TA法はしきい値を用いて解の
変換の受理、不受理を決定するため、解の変換による値
の変化量がしきい値より大きくなる局所解を乗り越える
ことができず、探索空間が滑らかでない最適化問題に適
用することは困難であったが、本発明によればTA法の
適用が可能となり、より最適な解を効率よく求めること
ができる。
変換の受理、不受理を決定するため、解の変換による値
の変化量がしきい値より大きくなる局所解を乗り越える
ことができず、探索空間が滑らかでない最適化問題に適
用することは困難であったが、本発明によればTA法の
適用が可能となり、より最適な解を効率よく求めること
ができる。
【0041】このように本発明によれば、効率的にしか
も質の高い最適割当解をえることができるので、ファジ
ー理論を用いた呼び割当てのルール抽出に或いはニュー
ラルネットを用いた呼び割当てにおいて、群管理に有効
な制御法則を学習する際の教師信号としてより有効に活
用することができる。
も質の高い最適割当解をえることができるので、ファジ
ー理論を用いた呼び割当てのルール抽出に或いはニュー
ラルネットを用いた呼び割当てにおいて、群管理に有効
な制御法則を学習する際の教師信号としてより有効に活
用することができる。
【図1】組合せ最適化問題としての乗場呼びの発生状況
の一例を示す図である。
の一例を示す図である。
【図2】本発明による最適解の探索手法を説明するため
のフローチャートである。
のフローチャートである。
【図3】図1の乗場呼びに対して本発明を適用した場合
の初期解の一例を示す図である。
の初期解の一例を示す図である。
【図4】図3の初期解に対する解の変換の一例を示す図
である。
である。
H1 〜H9 乗場呼び B1 〜B3 乗客 Ai 乗場呼びHi に割当てられるかごの号機番号 T しきい値 ta 各乗場呼びに対する待時間の平均値 t* 平均値ta の最適値
Claims (2)
- 【請求項1】 複数の既知の乗場呼びに対し、複数台の
かごを割当てる場合の最適解を探索し抽出する群管理エ
レベータの最適割当手法において、各かごに数がほぼ均
等となるように前記乗場呼びを任意に割当ててこれを初
期解とし、解の変換は前記乗場呼びのうち異なるかごに
割当てられた任意の2つの乗場呼びを入れ替えることに
より行い、これを探索の基本ステップとしてThreshold
Accepting 法を適用するこにより、最適解を探索するこ
とを特徴とする群管理エレベータの最適割当手法。 - 【請求項2】 任意の乗場呼びの1つを所定の確率で他
のかごに割当て変更することを解の変換として加える請
求項1記載の群管理エレベータの最適割当手法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3124959A JPH0790995B2 (ja) | 1991-04-26 | 1991-04-26 | 群管理エレベータの最適割当手法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3124959A JPH0790995B2 (ja) | 1991-04-26 | 1991-04-26 | 群管理エレベータの最適割当手法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH05319706A JPH05319706A (ja) | 1993-12-03 |
JPH0790995B2 true JPH0790995B2 (ja) | 1995-10-04 |
Family
ID=14898457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP3124959A Expired - Fee Related JPH0790995B2 (ja) | 1991-04-26 | 1991-04-26 | 群管理エレベータの最適割当手法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH0790995B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4690703B2 (ja) * | 2004-11-17 | 2011-06-01 | 株式会社東芝 | エレベータ群管理方法及び装置 |
US7546905B2 (en) * | 2006-03-27 | 2009-06-16 | Mitsubishi Electric Research Laboratories, Inc. | System and method for scheduling elevator cars using pairwise delay minimization |
JP2011144025A (ja) * | 2010-01-15 | 2011-07-28 | Toshiba Elevator Co Ltd | エレベータの群管理装置 |
-
1991
- 1991-04-26 JP JP3124959A patent/JPH0790995B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH05319706A (ja) | 1993-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110171753B (zh) | 一种电梯调度策略处理方法、装置、设备和存储介质 | |
JP4870863B2 (ja) | エレベータ群最適管理方法、及び最適管理システム | |
US6000504A (en) | Group management control method for elevator | |
CN102264620B (zh) | 电梯的群管理控制方法及装置 | |
JP4372681B2 (ja) | エレベータシステムを制御する方法および装置 | |
JP4494696B2 (ja) | エレベーター群管理装置 | |
EP1560778B1 (en) | Method for controlling an elevator system and controller for an elevator system | |
US5672853A (en) | Elevator control neural network | |
JP2008538737A (ja) | 予想交通情報を用いる群エレベータのスケジューリング | |
JP2002509848A (ja) | 仮想乗客トラヒックを生成するエレベータ群制御方法 | |
CN102139824A (zh) | 电梯群管理装置 | |
JPH02265876A (ja) | エレベータの群管理装置 | |
CN117252318B (zh) | 一种智能网联汽车群机协同拼车调度方法及*** | |
JP4621620B2 (ja) | エレベータ群管理システム、方法およびプログラム | |
Cortes et al. | Double deck elevator group control systems using evolutionary algorithms: Interfloor and lunchpeak traffic analysis | |
JPH06305649A (ja) | エレベータの運行制御方法及び制御装置 | |
JPH0790995B2 (ja) | 群管理エレベータの最適割当手法 | |
CN101356106A (zh) | 电梯控制方法及其装置 | |
WO2001072622A1 (fr) | Dispositif de commande de gestion d'un groupe d'ascenseurs | |
JPH07117941A (ja) | エレベータの群管理制御装置 | |
JPH11209010A (ja) | エレベータ群管理装置及び方法 | |
Sun et al. | Estimation of optimal elevator scheduling performance | |
JP3714343B2 (ja) | エレベータ群管理簡易シミュレータならびにエレベータ群管理装置 | |
Sorsa | A real-time genetic algorithm for the bilevel double-deck elevator dispatching problem | |
JPWO2019016937A1 (ja) | エレベーター装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071004 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101004 Year of fee payment: 15 |
|
LAPS | Cancellation because of no payment of annual fees |