srand（）でシードを修正すると、rand（）と異なる結果が出る

Question

私のJulia 0.5スクリプトでは、スクリプトを再実行するたびにrand（）と同じ結果を得るためにsrand（1234）を使用します。しかし、私は別の結果を得る。何が悪いですか？

Answer 1

@Dan Getzがこのコメントで述べたように、これはあなたが知らないうちにランダム関数を呼び出すコードがあるためです。

あなたは同じ種セットと同じrand()関数を呼び出す場合、期待どおり、あなたは同じ結果を得る：

julia> for i in 1:3 
      srand(1) 
      println(rand()) 
     end 
0.23603334566204692 
0.23603334566204692 
0.23603334566204692 

しかし、あなたはしてもしなくてもよいかもしれランドにスクリプト内の別のコールを持っている場合調査したrand()コールに行くと、乱数ジェネレータはさまざまな段階にあります。ループの2回目の反復で異なる値の乱数発生器と結果をオフセットし、余分なrand()コールがありますか

julia> for i in 1:3 
      srand(1) 
      if i == 2 
       rand() 
      end 
      println(rand()) 
     end 
0.23603334566204692 
0.34651701419196046 
0.23603334566204692 

注意：ここではこれを説明するための例です。

Answer 2

@ niczky12によって与えられた答えに加えて、あなた自身のジェネレータを定義し、より良い再現性のためにそのジェネレータを使用することをお勧めします。あなたのコントロール）は、あなたが得る乱数には影響しません。例えば

、シード1234のMersenneTwister作成：次に、あなたは、単にあなたのrand呼び出しに、この発電機を渡す

rng = MersenneTwister(1234) 

を：

julia> rng = MersenneTwister(1234); 

julia> rand(rng) 
0.5908446386657102 

julia> rand(rng, 2, 3) 
2×3 Array{Float64,2}: 
0.766797 0.460085 0.854147 
0.566237 0.794026 0.200586