2Dランダムウォーク 戻る
確率・拡散

2Dランダムウォーク

多数のウォーカーが原点から広がる様子をリアルタイム描画。⟨r²⟩ vs ステップ数Nの直線関係(拡散の証)を測定し、終点分布がガウス分布へ近づく中心極限定理まで一目で確かめられます。

パラメータ設定

プリセット
歩行タイプ
ウォーカー数 120
ステップ長 a 4
フレーム毎ステップ 2

このツールが見せるもの

原点から出発した多数のウォーカーが各ステップで独立にランダムな向きへ進みます。雲のように広がる軌跡、⟨r²⟩ ∝ N の直線、そして終点分布がガウス分布へ収束する様子(中心極限定理)が同時に観察できます。
ライブ数値
0
ステップ N
0
ウォーカー
0.0
⟨r²⟩
0.00
⟨r²⟩/N(測定)
0.0
RMS変位
リアルタイム・ランダムウォーク(残像つき軌跡)
ウォーカー軌跡 原点 中心の白十字=出発点(原点)。
平均二乗変位 ⟨r²⟩ vs ステップ数 N
実測 ⟨r²⟩ 理論 a²·N(線形拡散則) ⟨r²⟩ ∝ N(拡散):直線の傾き=a²。
終点 x の分布 → ガウス(中心極限定理)
終点xヒストグラム ガウスフィット 終点分布→ガウス(多数の独立ステップの和)。
理論・主要公式

各ステップが独立で平均0・分散一定なら、2Dの平均二乗変位はステップ数に比例して増えます:

$$\langle r^2(N)\rangle = N\,a^2 = 4Dt$$

ここで $a$ はステップ長、$D$ は拡散係数で $D=a^2/(4\,\Delta t)$。ログ-ログプロットでは傾き1の直線になり、これが通常拡散の指紋です。各成分の分散は $\langle x^2\rangle = N a^2/2$ となり、終点分布は多数の独立ステップの和なので中心極限定理によりガウス分布に収束します。原点回帰は2Dでは確率1(Pólyaの定理)。

検証例:a=4、N=200 → ⟨r²⟩≈Na²=3200。傾き⟨r²⟩/N≈a²=16に収束します。

2Dランダムウォークとは

🙋
ランダムウォークって、ただ酔っぱらいがフラフラ歩くだけのモデルですよね?これが何の役に立つんですか?
🎓
大まかに言うと、不規則な動きを理解するための「ものさし」なんだ。例えば、細胞内で栄養分子がどう拡散するか、株価の変動、昆虫の餌探しの動きまで、幅広く応用されているよ。このシミュレーターで「歩行タイプ」を「格子」から「連続」に変えてみると、動きの滑らかさが大きく異なるのがわかる。
🙋
なるほど。でも、グラフにある「平均二乗変位(MSD)」って何を見てるんですか?
🎓
MSDは、粒子が平均的にどれだけ遠くまで広がったかを数値化したものだ。実務では、実験で計測した粒子の動きから拡散の速さ(拡散係数)を求めるのに使うことが多いね。シミュレーター右側のプロットを見て。⟨r²⟩はステップ数Nに比例して直線的に増える——これが拡散の指紋だ。測定した傾き⟨r²⟩/Nがステップ長の2乗(a²)に一致するのを確かめてみよう。
🙋
なるほど!その下の「終点xの分布」のグラフは何を表しているんですか?
🎓
それぞれのウォーカーの最終位置(x座標)のヒストグラムだよ。1歩1歩はランダムでも、多数の独立なステップを足し合わせると、終点の分布はきれいな釣鐘型——ガウス分布に近づく。これが「中心極限定理」だ。ウォーカー数を増やして再生すると、黄色いガウスフィット曲線にぴたりと重なっていくのが見える。標準偏差σは√(N·a²/2)で予測できるよ。

よくある質問

ウォーカー数を増やすとMSDの統計的ばらつきが減り、理論曲線に近づきます。ステップ数を増やすと長時間の拡散挙動を観察でき、ノイズが平滑化されますが、計算負荷が上がるため、目的に応じて調整してください。
格子は固定長の直交方向のみ移動し、通常拡散を示します。ガウスは連続的なランダムな方向と距離でブラウン運動を模倣します。レヴィフライトは稀に長距離ジャンプを含み、超拡散(MSDが時間の2乗以上に比例)を再現します。
初期の過渡領域ではウォーカーの位置が原点付近に集中するため、MSDが理論直線から外れることがあります。また、レヴィフライトでは長距離ジャンプの頻度が低いと、統計的に不安定になり曲線がギザギザになる場合があります。
生体内の分子拡散、生態系での動物の探索行動、金融市場の価格変動など、ランダムな移動パターンの解析に利用できます。特にレヴィフライトは、効率的な探索戦略のモデルとしてロボット経路計画にも応用されています。

実世界での応用

生物物理学・生態学:細胞内でのタンパク質やmRNAの拡散、微生物や動物の採餌行動の分析に使われます。レヴィフライトは効率的な探索戦略として注目されています。

材料科学:多孔質材料内でのガスや液体の浸透・拡散過程のモデル化。不規則な構造の中での粒子の動きを予測するのにランダムウォークが応用されます。

金融工学:株価や為替レートの変動をモデル化する基礎として使われます。特に、大きな価格変動(ジャンプ)を含むモデルではレヴィフライトが参考にされます。

高分子科学:ポリマー鎖の構造や動的特性を、鎖を構成するモノマーのランダムウォークの集合体として統計的に理解するアプローチがあります。

よくある誤解と注意点

まず、「拡散係数Dは粒子の種類だけで決まる定数」と思いがちですが、それは誤解です。 例えば、同じ分子でも、水の中と油の中では拡散の速さは全く異なります。このシミュレーターで「歩行タイプ」を「連続角度」にし、「ステップ長」のスライダーを動かしてみてください。粒子の動きの「活発さ」が変わりますよね?本ツールでは拡散係数は D=a²/(4Δt)(aはステップ長)に対応し、実際の拡散係数が「環境(粘性など)」と「温度」に強く依存することを意味しています。実務では、シミュレーションで使うDの値を実験値や文献値から正しく設定することが第一歩です。

次に、MSDのグラフを1回の実行で判断しないこと。 特にウォーカー数が少ない場合、MSDのプロットはかなりばらつきます。「格子」歩行でウォーカー数1で実行すると、MSDが綺麗な直線にならないことが確認できるはずです。これは統計的なゆらぎです。信頼性のある結果を得るには、十分な数のウォーカー(例えば100以上)でシミュレーションを走らせ、平均を取る必要があります。

最後に、終点分布がガウスに見えるには十分なステップ数が必要です。 中心極限定理は「多数の独立な確率変数の和」が正規分布に近づくことを保証しますが、ステップ数Nが小さいうちは、特に格子歩行で離散的な凹凸が残ります。Nを増やしていくと、終点xのヒストグラムが黄色いガウスフィット曲線に滑らかに重なっていくのが観察できます。実現象への適用時は、観測時間が十分長く「拡散領域」に入っているかを確認するプロセスが不可欠です。

使い方ガイド

  1. ウォーカー数(valN)を1~400の範囲で設定。統計的ばらつきを抑え終点分布をガウスに近づけるには100以上が有効です
  2. ステップ長(valStep)で1歩あたりの移動距離(1~12)を決定。総ステップ数Nは「再生」で時間とともに増えていきます
  3. 歩行タイプで格子(4方向/8方向)と連続角度を切り替え、軌跡の滑らかさと⟨r²⟩の傾き(a²)の関係を比較
  4. フレーム毎ステップで時間発展の速さを調整し、⟨r²⟩ vs N プロットが直線(傾き1)に乗るのを確認します
  5. プリセット(単一/雲/格子/連続)で典型シナリオを即座に再現。終点xヒストグラムがガウスフィットへ重なる様子を観察します

具体的な計算例

ウォーカー数N=120、ステップ長a=4、格子(4方向)タイプの場合、t=200ステップ後のMSD期待値はおよそ⟨r²⟩≈N·a²=200×16=3200(正規ランダムウォーク:⟨r²⟩∝N)。ドリフトなしならログ-ログプロットの傾きは約1で、⟨r²⟩=4Dt から拡散係数D≈a²/4=4が読み取れます。測定した⟨r²⟩/Nはステップ長の2乗a²=16に収束します。終点x座標の標準偏差は√(N·a²/2)=√1600=40となり、黄色いガウスフィット曲線と一致します。

実務での注意点

  1. 拡散係数測定では十分なウォーカー数を確保。N<10では統計ノイズが支配的になり、MSDの信頼区間が±30%超に拡大
  2. 長時間シミュレーション(ステップ>10000)では浮動小数点誤差蓄積に注意。周期的な位置ラップアラウンドが必要な場合は格子模型を選択
  3. 生体分子の異常拡散(核内DNA運動)やコロイド粒子のブラウン運動解析では、観測時間が拡散領域に入っているかを確認し、過渡領域を除外して傾きを推定
  4. ドリフト存在下では見かけの拡散係数が過大評価。MSD曲線から線形成分を差引いた残差で真の拡散のみを抽出

🎬 動画で見る

ブラウン運動|でたらめな揺れが時間に比例して広がる #Shorts
ブラウン運動|でたらめな揺れが時間に比例して広がる #Shorts
拡散律速凝集|ランダムな粒子が樹枝の結晶を編む #Shorts
拡散律速凝集|ランダムな粒子が樹枝の結晶を編む #Shorts