乱数の特性について
2007/12/23
これまでのシミュレーションでは、乱数に メルセンヌ・ツイスターを用いてきた。 メルセンヌ・ツイスターは「従来にない長周期, 高次元均等分布を」持った、性質の良い乱数である。 それでは、反対に、周期が短く不均等な、性質の悪い乱数を用いてシミュレーションを行ったらどうなるだろうか。 シミュレーション8、シミュレーション9では、 連成振動子の列を用いてシミュレーションを行ったが、一方向の流れは確認できなかった。 その理由は、おそらく連成振動子の列が”性質の良い乱数”になっていないからではないか。 であれば、”性質の悪い”乱数を用いたならば一方向の流れが生じない、といった予想が成り立つ。 上の予想を確かめるため、シミュレーション12で分子数2個のケースを、性質の悪い乱数に置き換えて実施した。 性質の悪い乱数として、線形合同法による次のような乱数列R[n]を用いた。
R[n+1] = ( R[n] x A + B ) mod MD
乱数列と言うより、むしろ周期的な数列に近い。
A = 9, B = 5, MD = 16 R[0] = 8
結果のグラフ
結果は予想に反して、正(+)の転送が高い頻度で起こっていた。 なぜこのような結果になったのか。 乱数列の周期と、仕掛けの動作周期がたまたま一致した場合には、このように綺麗に直線的な結果が出るのだろうか。 いまのところ良い説明が思い付かない。 1つだけ確かなのは、乱数の善し悪しによって結果が左右される、ということだ。 |