第６回 アリの採餌モデル（１１月１５日）
＊出席確認
＜出席確認メール＞ 添付ファイルなし 受付：本日１４：３０〜１６：１５
本日のキーワードは授業中に連絡します。
宛先：miwamoto[at]riko.shimane-u.ac.jp
件名：MMM2 20161115
本文： 学生番号 名前 「本日のキーワード」
ランダムウォークを利用した蟻の行動
ランダムウォークを利用して蟻の行動をシミュレーションしよう。
餌を探しに出た蟻はどのようにして巣に帰るのか？
次のような観察結果がある。
- 蟻は巣から出て、餌を探す時にはランダムウォークする。
- 巣に帰るときには何故か巣の方向に向かってほぼ一直線に戻る。
さらに観察すると、次のような事実がわかる。
- 蟻は餌を求めるためにランダムウォークで動き回るが、動いたところにある種の化学物質（フェロモン）をつけて、足跡を残す。
- （餌を見つけるか、あるいはあきらめて）巣に帰るときには、付けた化学物質（足跡）の多い方へ向かって進む。（このとき、餌を見つけた蟻はそのことを知らせるために、別の異なる化学物質を足跡としてつける。）
（注意）蟻と言えば蟻の行列を思い浮かべる。上記の話は、行列を作る前の話である。
蟻は餌を探す時には上記の様なランダムウォークを行う。そのとき足跡フェロモンを分泌する。餌を見つけると、上記のようなルールで巣に戻るのであるが、そのときには道しるべフェロモンという別の物質を分泌し、巣に帰った後にその道しるべフェロモンをたどることで餌の場所に行列をつくりつつ、さらに道しるべフェロモンを強化すると言われている。
練習１ ２次元配列を使って、１匹の蟻がランダムウォークした軌跡を描画するシミュレーションを作成せよ。
ヒント１） ２次元配列の宣言
int x[][] = new int[particlenum][step];
ヒント２） 時間ステップ数
frameCount
ヒント３） 用意した配列より大きい数が入るとエラーになってしまうので、例えば1000ステップで止めると決めて、繰り返しループを止める。
if(frameCount >= 1000) noLoop();
補足：２次元配列について
例えば、
int a[2][3];
と宣言すると、
a[0][0], a[1][0],
a[0][1], a[1][1],
a[0][2], a[1][2]
の計６つ(2×3) の整数型の変数が用意される。
６つの変数に２重for文で値を代入することができる。
for(int i = 0; i < 2; i++)
{
for(int j = 0; j < 3; j++)
{
a[i][j] = i*j*sqrt(2.0);
}
}
これによって、
a[0][0] = 0
a[1][0] = 0
a[0][1] = 0
a[1][1] =
a[0][2] = 0
a[1][2] =
が代入される。
足跡を数える
蟻は移動するたびに足跡フェロモンを残すということが観察結果からわかっているので、足跡を数える仕組みをランダムウォークシミュレーションのプログラムに組み込もう。
考えている格子点全てについて、足跡を保持している変数を用意してやればよい。この場合、２次元配列変数を使うと簡単である。具体的には、
int count[][] = new int[x方向の格子の数][y方向の格子の数];
として宣言した配列変数 count に足跡の個数を記録する。
練習２ 蟻の足跡フェロモンの量を表示するシミュレーションを作成せよ。
ヒント１） 初期値として、count[i][j]の全てに0を代入せよ。
ヒント２） 蟻がその場所に来たらcount[蟻のx格子][蟻のy格子] += 1;
ヒント３） 塗りつぶしにcount[i][j]を使って、フェロモンの量を色で表現せよ。
for(int i = 0; i < 100; i++)
{
for(int j = 0; j < 100; j++)
{
fill(255-min(255,2*count[i][j]));
rect(i*5, j*5, 5,5);
}
最小値を求める
min(a,b)
最大値を求める
max(a,b)
最大値最小値を求めるときは、maxやminを使用することが可能である。a,bに入るのは、int型かfloat型（同じである必要がある）で、最大３つまでの数を比較できる。
練習３ M匹の蟻が一斉に巣から出てくる状況を考える。このとき、各格子の足跡の合計をカウントし、画面に表示するプログラムを作成せよ。蟻の衝突は考えない。
巣に帰る蟻
足跡の分布は巣から離れるにつれて減少しているように見え、この傾向は蟻の数Ｍが大きくなるにつれ強くなった。
蟻の帰巣行動について、蟻は次の様なルールに従い、巣に帰るとする。
- 蟻はある回数動くと、餌を見つけられたか否かによらず、（疲れたので）巣に帰る。
- 巣に帰るときには化学物質量（足跡数）の多い方に向かって移動する。
蟻の帰巣行動をシミュレーションするための手順を考えよう。
最終的に到達した地点から帰巣行動をはじめるため、Ｍ匹の蟻それぞれがランダムウォークの結果到達した最終到達点の座標を記録しておく必要がある。
ヒント１）最終到達点の座標を記録しておく配列変数を用意する。
ヒント２）delayで指定した時間計算を停めることができる。括弧の中には数字を入れる。単位はミリ秒なので、次のようにかくと100ステップ目で３秒間計算がとまり、その後また始まる。
if(frameCount == 100) delay(3000);
本日の課題１ 巣からランダムウォークして出て行った蟻が5秒間止まったのち、巣に帰るシミュレーションを作成せよ。蟻は４方向をみてフェロモンの量が多い格子へ進むとする。
以下の点を工夫せよ。
① ４方向のうち、２方向以上のフェロモンの量が同じであった場合どうするか。
② ４方向に比べて自分のいるところの方がフェロモンの量が多い場合はどうするか。
本日の課題２（発展） ランダムウォークしたM匹の蟻のうち、ある時間内に何匹が帰巣に成功するか調べよ。蟻の数を増やしたら帰巣に成功する蟻の数はどうなるか調べよ。