Boid にノイズを加えて相転移させた
はじめに
次の文献の追実験をしてみました.
ちなみに,以下で書いていることの中には,元の論文に書いて無い私の感想や解釈なども混じっているので注意してください.
文献
- タイトル: Novel Type of Phase Transition in a System of Self-Driven Particles
- 著者: Tamas Vicsek, Andras Czirok, Eshel Ben-Jacob, Inon Cohen and Ofer Shochet
- 論文誌: Physical Review Letters, Volume 75, Issue 6, August 7, 1995, pp.1226-1229 (PhRvL Homepage)
- 論文のPDF: http://angel.elte.hu/~vicsek/pdf/vicseketal95prl.pdf
Boid(鳥もどき)は鳥の群れのような動きを,各鳥に簡単な行動ルールを与えるだけで再現するモデルで,複雑系の代表的なモデルの一つです.
やったことは,この記事のタイトル通りで,各鳥の動きにノイズを加えて相転移現象を観察してみました.
この論文では,単純な非平衡自己組織化モデルとして,Boid(自己組織化による鳥の群れ)を用いて,そこに内在している相転移の性質を調べるのが目的のようです.
このモデルそのものは現実の具体的な現象と直接対応するわけではありませんが,私たちの身の回りには非平衡な自己組織化のように見える現象はたくさんあり,非常に重要かつおもしろい現象だと思います.
特に,心臓のポンプのような動き(各細胞が同期して振動),私たちの思考(ニューロンによる動的なパターン形成?)などのように多くの生命現象は非平衡な自己組織化現象だと思います.
モデル
この論文では一般的な Boid をよりシンプルにしたモデルを用いています.シンプルですが Boid の特徴的な振る舞いを十分に表していると思います.
おおまかにどのようなものかというと,
- エージェント(鳥)は一定の速さで移動する.
- 自分の近く(半径 r 以内)にエージェントがいれば,それらの進行方向の平均に自分の進行方向θを変える.
- ただし,進行方向には常に一定のノイズηが加えられる.
というものです.
以下で具体的に記述します.
エージェントはのトーラス空間を移動します.
エージェント i の移動は位置を次式に従って更新することで表現されます.
(1)
ここで,速度のノルムは一定で v とします.
の角度(移動方向)をとし,移動方向の変化は次式に従って更新することで表現します.
(2)
はエージェント i の半径 r 内にいるエージェントの方向の平均です.この項によってエージェントは同じ方向を向いて移動するようになります.
は範囲]の一様乱数です.これはエージェントの移動に関するノイズです.
実験
エージェントの群れ行動の秩序度合いを表す値として以下を考えます.
秩序度は全エージェントが同じ方向を向いていれば 1 に,完全にバラバラならば 0 になります.
よって,この値だけでエージェント群の振る舞いについて大まかに知ることができます.
ということで,ノイズηがエージェント群の秩序度に与える影響について見ていきたいと思います.
その他のパラメータは,v = 0.1,r = 1.0,L = 20,N = 1000,Δt = 1.0 に設定してシミュレーションしました.
十分に時間がたった後のの平均に関するηの影響を下図に示します.
ノイズηが小さいときはは 1 に近く,非常に秩序立った振る舞いをしていることが分かります.
がそれなりの値を持っているのはη = 3.0 ぐらいまでで,それ以上でははほとんど 0 になっています.
この図ではそれほどはっきりしていませんが,Nを非常に大きくすると,がそこそこの値のときと 0 に近いときの境界は明確になり,それは相転移であることが予想されます.
論文中ではこの境界(臨界点)においていろいろ調べていて,相転移の臨界点における特徴的な性質が存在していることを示しており,ほぼ相転移であると言っていいと思います.
次に,η = 0.01のときのスクリーンショットを示します.水色の軌道は適当に選んだエージェントの過去300ステップ(300Δt 時間)の軌道です.
非常に大きな一つの群れとなり,また水色の軌道も方向性を持っていることから,秩序立った振る舞いをしていることがわかります.
以下はη = 2.5のときのスクリーンショットです.
η = 0.01と比べるとややバラバラしてますが,群れがある程度できていること,水色の軌道からそれなりに方向性をもっていることから,ある程度は秩序立った振る舞いをしていると言っていい手思います.
以下はη = 4.0のときのスクリーンショットです.
エージェントの向き・位置関係はバラバラで,また,水色の軌道もランダムに動いているだけのように見え,秩序は無いと言っていいと思います.