Transcript eil51 - 同志社大学

温度並列シミュレーテッドアニーリング
における重要温度
三木 光範 （同志社大学）
廣安 知之 （同志社大学）
○ 窪田 耕明 （同志社大学大学院）
吉田 武史 （同志社大学）
研究背景
最適化問題
ヒューリスティック探索手法
シミュレーテッドアニーリング（SA） ，GA，etc.
高い計算負荷で並列化は必須
並列アルゴリズム
並列SA ，分散GA，etc.
温度並列SA（TPSA）
SAの基本アルゴリズム
SA：焼きなましのアナロジーを模倣
・生成処理
・受理判定
減少
エネルギーが
・温度のクーリング
小
の確率で受理
SAの欠点
大
Energy
増加
必ず受理
－(Ｅnext－Ｅcurrent)
exp
Ｔemperature
・計算時間がかかる
・温度設定が困難
局所解
最適解
温度並列SA（TPSA）
T1
Temparature
TPSA
SA
・生成処理
・受理判定
・解交換
T2
T3
T4
T5
T6
Time
各プロセスに温度を割り当てる
TPSAの利点
・並列処理との親和性が高い
・温度スケジュールの自動化
温度パラメータの設定が重要
Temparature
TPSA
T1
T2
T3
T4
T5
T6
Time
重要温度
Temparature
一定温度でのアニーリングが有効
T1
T2
T3
T4
T5
T6
研究目的
無駄な探索
効率の良い探索
（重要温度）
局所解に収束
Time
・巡回セールスマン問題において重要温度の存在を確認する．
・重要温度の規則性を検証する．
TPSAの適切な温度パラメータ設定が可能！
巡回セールスマン問題（TSP）
eil51（５１都市，最適解：４２６）
初期解
最適解
80
80
70
70
60
60
50
50
40
40
30
30
20
20
10
10
0
0
0
10
20
30
40
50
60
70
80
0
10
20
30
40
50
60
70
80
重要温度を確認する実験
解交換を行なわないTPSA
使用するパラメータ
Temparature
T1
T2
最高温度：1.0E＋6
T3
最低温度：1.0E－2
T4
温度数：32
T5
T6
計算回数：160000回
Time
実験結果
経路長
eil51　(TSPLIB)
750
700
650
600
550
500
450
400
1.0E-02 1.0E-01 1.0E+00 1.0E+01 1.0E+02 1.0E+03 1.0E+04 1.0E+05 1.0E+06
温度
重要温度を確認できた
各問題での重要温度
問題
都市数
eil51
最適解
重要温度
426
2
kroA100
100
21282
50
pｒ152
152
73682
130
bier127
127
118282
180
76
108159
300
pr76
１．最適解に比例？
２．都市数に反比例？
最適解／都市数
＝
51
重要温度は
平均経路長
重要温度は平均経路長と関係がある？
問題
eil51
平均経路長と重要温度
都市数
最適解
重要温度
51
426
2
kroA100
100
21282
50
pr152
152
73682
130
bier127
127
118282
180
76
108159
300
pr76
（pr76）
（bier127）
（pr152）
（kroA100）
（eil51）
1
傾き 4 の直線
平均経路長と重要温度
・一般的なTSPでは
最適解と都市数（平均経路長）から
重要温度を導くことができる
1
重要温度 ＝ ×（最適解／都市数）
4
特殊な問題ではどうか？
テスト問題による検証
eｉｌ51を使用してテスト問題を作成
都市数を2倍にしたｅｉｌ51
・各都市に近接して都市を付加 （ｅｉｌ51-A）
・各都市の中間点に都市を付加 （ｅｉｌ51-B）
解交換を行なわないTPSAを適用
eil51-A
各都市に近接して都市を付加(102都市)
eil51-B 最適解：426
80
70
60
50
40
30
20
10
0
0
10
20
30
40
50
60
70
80
eil51-B
各都市の中間点に都市を付加(102都市)
eil51 最適解：426
80
80
70
70
60
60
50
50
40
40
30
30
20
20
10
10
0
0
0
10
20
30
40
50
60
70
eil51-C 最適解：426
80
0
10
20
30
40
50
60
70
80
実験結果
実験結果
・eil51-A → 重要温度は変化しない．
重要温度は都市数に反比例していない？
・eil51-B → 重要温度の幅が広くなる．
重要温度は高くてもいい？
考察
1
1
最適解：８
都市数：８
平均経路長 ＝ 1
実際の各経路長
＝
最適解：８
都市数：８
平均経路長 ＝ 1
実際の各経路長
≠
1
1
2.1
1.8
0.2
計算で求めた
平均経路長
計算で求めた
平均経路長
eil51の経路長分布
10
平均経路長＝426÷51＝8.4
頻度 (回 )
8
6
4
2
0
0.1
2.1
4.1
6.1
8.1
10.1
12.1
各経路長
14.1
16.1
18.1
eil51-A(近接)の経路長分布
35
平均経路長＝426÷102＝4.2
30
頻度 (回 )
25
20
15
10
5
0
0.1
2.1
4.1
6.1
8.1
10.1 12.1 14.1 16.1 18.1
各経路長
考察（eil51-A）
35
平均経路長＝426÷102＝4.2
30
頻度 (回 )
25
20
温度に関係なく結ばれる経路長
15
10
5
0
0.1
2.1
4.1
6.1
8.1
10.1
12.1
14.1
16.1
18.1
各経路長
経路長にばらつきがあり
平均経路長と異なる
平均経路長から重要温度は
計算できない！
eil51-B（中間点）の経路長分布
10
平均経路長＝426÷102＝4.2
頻度(回)
8
6
経路長にばらつきがなく
平均経路長と等しい
4
2
0
0.1
2.1
4.1
6.1
8.1
10.1
12.1
各経路長
14.1
16.1
18.1
考察（eil51-B）
各都市の中間点に都市を付加
80
70
60
50
問題が簡単になった
（迷うような経路がない）
40
30
20
10
多少の温度差に影響を受けない
0
0
10
20
30
40
50
60
問題が簡単になる
70
80
局所探索で解が求まる
（重要温度はなくなる）
考察
eil51-A 小さい経路長は考慮されない
特殊な問題
重要温度が変わらない
eil51-B 問題が簡単になる
重要温度の幅が広くなる
平均経路長（最適解／都市数）が使えない
実際にどのような大きさの組み換えがあるのか？
経路長の組み換えによって生じる
エネルギー(ΔE)が重要
結論
・TSPにおいて重要温度がある
・特殊なTSPに関しては，組み換えで生じる
エネルギー(ΔE)が重要
・一般的なTSPにおいて，重要温度は
平均経路長（最適解／都市数）に依存している
1
重要温度 ＝ ×（最適解／都市数）
4
TPSAの適切な温度パラメータを決める指針を示した．