Mathematicaのループ文について

質問

現在Mathematicaを使って交通のシュミレーションをしているのですが、思ったようにループ文が作れずに悩んでいます。

もし条件式が１から４まであったとして、その条件式を毎回ループさせると考えたらどういった式を作ればよいでしょうか？

↓ちなみに条件式です。これを50回ループさせようと考えています。

road = Table[0, {2}, {54}];
road[[1, 1]] = 0;
road[[2, 1]] = 0;
road[[1, 2]] = 0;
road[[2, 2]] = 0;
road[[1, 53]] = 0;
road[[2, 53]] = 0;
road[[1, 54]] = 0;
road[[2, 54]] = 0;
road[[1, 3]] = a;
road[[2, 3]] = b;
If[0.5 < Random[], a = 1, a = 0];
If[0.5 < Random[], b = 1, b = 0];

回答No.2

tinantum

とりあえず，周期境界条件の件は

r = RotateLeft[IdentityMatrix[n]];
r[[n, 1]] = 0;

として，rを次々に掛ければよいのでは？　例えば

keepOnTrucking[p_, t_, n_] :=
Module[{road = Table[0, {n}]},
road[[1]] = a;
If[p < Random[], a = 1, a = 0];
r = RotateLeft[IdentityMatrix[n]];
r[[n, 1]] = 0;
road2 = road;
s = Table[
road2 = (road2.r).r,
{j, 0, t}]; Insert[s, road, {1}];
x = Table[0.5, {n}];
y = RotateRight[Append[s, road]]; y;
z = y;
Table[z = Insert[z, x, {2j}], {j, t + 1}]; z]
Show[Graphics[RasterArray[Reverse[keepOnTrucking[0.2, 30, 50] /. {0 ->
RGBColor[0.8, 0.8, 0.8], 1 -> RGBColor[
Random[], Random[], Random[]], 0.5 -> RGBColor[1, 1,
1]}]], AspectRatio -> Automatic]]

忙しいのと，RasterArrayなど慣れてないものなので，プログラム自体は不細工なものになってます．．

車の台数を複数にしたいのであればもう少し考えて見ます．
そのかわり，車のスピードは同じものでよいのか？
車のスピードは違うものを想定する場合，追い越しはありか？
ということを教えてください．
また，もし既に解決している場合，プログラムを教えていただけると私の勉強にもなりますので，よろしくお願いします．

それと，
drive[1, 0, _, _, _, _, _] = 0;
drive[1, 1, _, _, 0, _, 0] = 0;
drive[1, 1, _, _, 1, _, _] = 1;
drive[0, 1, 1, _, _, _, _] = 1;
drive[0, _, 0, _, 1, 1, _] = 1;
drive[0, _, 0, 1, _, _, _] = 1;
drive[x_, _, _, _, _, _, _] := x;
の使い方がいまいちわからないので教えてくれますか？

回答No.1

tinantum

やりたいことがいまいちわかりません．
書いてあることを素直にループさせると，変わるのはいつも
road[1,3]
road[2,3]
の要素だけですが，それでよいのですか・・？
そうでなければ，５０回のループ中の１，２回くらい書くか，行列roadの要素をどういう繰り返しをさせて変えたいのか書いてくれれば考えてみます．

補足

早速の回答有難うございます。
私が今考えているのは２列の配列を初期値０で設定し、一番端の部分をある一定の確率で１を出現させてセルオートマトン方を用いて交通のシュミレーションを作ろうと考えています。
今作成したプログラムはこんな感じです。

keepOnTrucking[p_, t_, n_] :=
Module[{road=Table[0, {2}, {n}], drive},
road[[1, 1]] = a;
road[[2, 1]] = b;
If[p < Random[], a = 1, a = 0];
If[p < Random[], b = 1, b = 0];
drive[1, 0, _, _, _, _, _] = 0;
drive[1, 1, _, _, 0, _, 0] = 0;
drive[1, 1, _, _, 1, _, _] = 1;
drive[0, 1, 1, _, _, _, _] = 1;
drive[0, _, 0, _, 1, 1, _] = 1;
drive[0, _, 0, 1, _, _, _] = 1;
drive[x_, _, _, _, _, _, _] := x;
NestList[MapThread[drive,
{#, RotateRight[#, {0, -1}], RotateRight[#, {0, 1}],
RotateRight[#, {0, 2}], RotateRight[#, {1, 0}], RotateRight[#, {1, -1}],
RotateRight[#, {1, 1}]}, 2] &, road, t]
]
Map[Show[Graphics[RasterArray[# /.
{0 -> RGBColor[0.8, 0.8, 0.8], 1 -> RGBColor[0.1, 0.1, 0.1]}]],
AspectRatio -> Automatic] &, keepOnTrucking[0.2, 50, 50]]

これで作動したら２つの問題点が出てきてしまうんです。
１つは周期的境界条件が発生して、最後に消えた車が次のステップで出発地点から出てきてしまう。
２つ目は２ステップ目から車が発生しないと言う点です。

もし分かるようでしたらアドバイスをお願いします。