• 検索結果がありません。

不確定環境型遺伝的アルゴリズムによる 確率的巡回セールスマン問題の近似解法

N/A
N/A
Protected

Academic year: 2021

シェア "不確定環境型遺伝的アルゴリズムによる 確率的巡回セールスマン問題の近似解法"

Copied!
2
0
0

読み込み中.... (全文を見る)

全文

(1)

2002年日本オペレーションズ・リサーチ学会 秋季研究発表会 2−C−5

不確定環境型遺伝的アルゴリズムによる

確率的巡回セールスマン問題の近似解法

防衛庁

河原一将 KAWAHARAKazumasa

京都府立大学 *吉富康成 YOSHITOMI Yasunari 01704426

団を生成する。

(c)可能な親の対を全て作り、それぞれの親の対

で交叉を行う。ここで行う交叉は、サブツアー

交換交叉【胡であり、二つの親の間で部分的に同

じ遺伝子の集合があれば交換し、遺伝子を並び

替えた新たな個体(子)を作るというのである。

(心全ての子の適応度を計算し、適応度の高い子か

ら順に、所与の交叉確率に従って次世代の個体

として残し、残りの個体は親候補からランダム

に選択する。

b)所与の突然変異確率に従って突然変異を行う。

ここで行う突然変異は子に対し、ランダムに 選んだ二つの遺伝子を入れ替えるというもので

ある。

④全世代での各個体伐牢)の発生頻度を求める。

本研究では確率計画問題において、期待値最大の

解を得ることを目標とした。このため、選択方式と

して、ルーレット戦略をとり、全世代を通じて発生

頻度が最も高い個体を、期待値最大を与える個体(最

適解)の近似最適解とした。

3.実験

3.1 条件

サイズ1.0×1.0の領域にランダムに都市を10ま

たは100配置した。都市間の移動時間の確率分布を、

(め距離の逆数を平均値とした正規分布とする、(0

この平均値を一部の都市間のみ十分大きな備に置き 換える確率を0.5とする、の2つの場合を検討した。 これらの条件では、最適解は確率変動にかかわらず 同じとなる。本法の性能評価のためこれらの条件を

用いる。個体数を10、交叉確率を0.4、突然変異確

率を0.1、世代数を100α10都市)または200α1(刀

都市)とし、正規分布における変動係数α確率変動な

し)、0.1、0.2、0.3の条件で実験を行った。

都市数10の場合3つの都市配置条件で、確率変

動を与えない場合の本法の近似最適解と、ツアーの すべての場合を尽くした解のうち、巡回時間の最/j、

値を与える解(最適角覇とを比較した。そして、確率

変動を与えた場合での本法の近似最適解を求めた。 1.緒言

著者らは、確率計画問題の近似解法として、遺伝

的アルゴリズム(G刃の環境(目的関数、制約条件)に

確率変動を導入した手法(不確定環境型G刃を提案

した【1】。本法では、世代ごとに、目的関数、制約条

件で定義される適応度関数を所与の確率分布に応じ

て変化させ、全世代を通じての個体の集合とその出

現頻度を算出する。そして、まずこれにより、期待

値最大の解が得られるかどうかの検討を行なった。

その結果、選択方式として、ルーレット戦略の下で、

出現頻度が最も高い個体伐劉を選べば、それが期待

値最大を与える個体となることを実証し、本法を確

率的画像圧縮問題、確率的ジョブショップ問題など

へ適用し、その有効性を示した【1,2,3]。

本報では、確率的巡回セールスマン問題に不確定

環境型GAを適用し、その有効性を検討した。本研

究では、都市間の移動時間を不確定とし、全ての都

市の巡回時間の期待値を最小にする問題を対象にし た。

2.GAの環境への確率変動の導入と計算手順

確率計画問題において、確率変数の変動にともな い解の目的関数値が変動することを、GAにおいて は、同じ個体の適応度が確率的変動を含んでいると

考えることとする。本研究では、セールスマンが、

決められた数の都市を巡回する。そして、都市の順

序(ツアー)を与えられた個体の、巡回時間の逆数を

適応度とし、世代ごとに都市間の移動時間を確率的

に変動させることにより、適応度の確率的変動を実

現する。すなわち、GAの各世代の環境が不確定で

あるとして取り扱う。そして、全世代を通じての個

体の集合とその出現頻度を算出した。計算手順を以

下に示す。個体の各遺伝子は、都市と対応している。

①都市を配置する。

②初期集団の生成を行う。

③所与の世代まで以下の処理を行う。

扉確率分布に従う乱数を用いて、都市間の移動時

間を確定し、適応度を計算する。

他)適応度の高い個体脚が確率的に選ばれや すいルーレット戦略によって、新しい個体の集 ー202− © 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

(2)

次に、都市数100の場合に、確率変動を与えない

場合の、初期集団の最良ツアーと、10(カ世代目の最

良ツアーを比較した。また、各確率変動条件で、最 大適応度の収束を調べた。

3.2 結果と考察

都市数10の場合、最適解と本法の近似最適解は

11/16の割合で一致した。不一致の場合でも、最適

解の頻度は上位4番以内であった。そして、10(X沿仁

個体数×世代数)個の総個体のうち、494∼6414は

最頻度角牢であった。この最頻度解の個数は、変動係

数が大きい程少なかった。また、この最頻度解の個 数が少ない場合に、本法の近似最適解と最適解が一 致しにくい傾向があった。不確定環境型GAで高頻

度解を所定量選び、その解の巡回時間の期待値の近

似値をモンテカルロ法で求めることにより、それら の解から最適解を選択することができた。 4 2 1 00 6 4 2 0 ﹂ ﹂ 〇..〇 ∩岬 0 0 0 ▲U O O O O 塑︽咄Q世擾憎 0 500 1000 1500 2000 世代 図2 適応度の変化 都市を100とした場合、平面上の都市とツアーは 図1のようになり、初期集団と比べると10(カ世代 目では経路の交わる部分がなくなり、より良いツア ーが得られた(変動係数=功。また、図2に示すよう に、世代が進むにつれ、適応度の最大値が収束した (変動係数=0.2)。 4.結論 確率的巡回セールスマン問題に不確定環境型GA を適用した。巡回時間の期待値を最小にする問題を 対象にし、本法の有効性を実証した。今後の課題と しては、所定の時間以上に巡回時間がかかる確率を 最小にする解なども決定できるように選択等GAの オペレータを検討することなどがあげられる。 参考文献 【1】Y.Yoshitomi,H.Ikenoue,T.Takebaand S.Tbmita, ‘‘GeneticjugorithminUncertain Envirt)nmentSfor SolvingStochasticProgrammingProblem”,日本オペ レーションズ・リサーチ学会論文誌,43(2000),266−29α 【2】吉富康成,“不確定環境型遺伝的アルゴリズムによる確 率的ジョブショップ問題の解法”,スケジューリング・ シンポジウム’99講演論文集(1999),119・124, 【3】吉富康成,山口理絵,“不確定環境型GAとヒューリステ ィック法による確率的ジョブショップ問題の近似解 法”,日本オペレーションズ・リサーチ学会秋季研究発 表会アブストラクト集,(200¢),%・97. k仙村雅章,小野貴久小林重信,“形質の遺伝を重視した遺 伝的アルゴリズムに基づく巡回セールスマン問題の解法,, 7(1992),1049・1059. 図1 最良ツアー −203− © 日本オペレーションズ・リサーチ学会. 無断複写・複製・転載を禁ず.

参照

関連したドキュメント

2 E-LOCA を仮定した場合でも,ECCS 系による注水流量では足りないほどの原子炉冷却材の流出が考

・逆解析は,GA(遺伝的アルゴリズム)を用い,パラメータは,個体数 20,世 代数 100,交叉確率 0.75,突然変異率は

・ 津波高さが 4.8m 以上~ 6.5m 未満 ( 津波シナリオ区分 3) において,原

炉心損傷 事故シーケンスPCV破損時期RPV圧力炉心損傷時期電源確保プラント損傷状態 後期 TW 炉心損傷前 早期 後期 長期TB 高圧電源確保 TQUX 早期 TBU

表4.1.1.f-1代表炉心損傷シーケンスの事故進展解析結果 PDS 炉心溶融 RPV下部プレナム リロケーションRPV破損 PCV破損 TQUV (TBP) TQUX (TBU、TBD) TQUX (RPV破損なし)

【A2】 ROV 北回りル ートから ペデスタ

【A2】 ROV 北回りル ートから ペデスタ

地震 L1 について、状態 A+α と状態 E の評価結果を比較すると、全 CDF は状態 A+α の 1.2×10 -5 /炉年から状態 E では 8.2×10 -6 /炉年まで低下し