拡散のシミュレーション このときはランダムウォークから拡散方程式まで このときは解析的に解いたあとでグラフを書いたもの 今回は拡散方程式を解析的にではなく数値的に計算していく いろいろな項を考えているがまずは拡散項のみ 独立変数： 従属変数： の…

こちらで解を探すときに振動してしまった ノイズをいれると振動しなくなった この振動の原因は何か 振動しているとき、特定の状態を行ったり来たりしている 一方ランダムな要素を入れたら行ったり来たりしなくなる なので推移確率行列を使って、特定の場所を…

これは拡散の様子のシミュレーション 密度に書き換えた z軸の上限を0.02にしてプロットしてみた ここまではランダムウォークのシミュレーション その場にとどまる確率： 各方向距離動く確率： その他： 拡散方程式にする 仮定など 距離関数： 初期分布： 境…

昨日のつづき（こちら） 昨日は独立の運動 今回は次ステップが周囲の状況に依存する いくつか入れてみたいパラメータも書いてみた 局所ルールが全体にどれだけ影響するか 組織的な形、動きとなるかどうか N<-40 M<-40 T<-30 C<-B<-A<-matrix(0,N,M) for (i i…

ImageJというソフト（ImageJ）（Wiki） アニメーションを作ってみた こちらを参考にRでファイルを書き出す 中心に集中していた点がランダムに動いて散らばっていく様子 #NかけM行列 時間T N<-30 M<-30 T<-50 B<-A<-matrix(0,N,M) A[N/2,M/2]<-200 #最初の分…