我想使用mersenne twister在10到50之间生成'N'个随机数。我希望能够一遍又一遍地生成相同的序列。
我写了以下代码:(种子= 50,a = 10,b = 50,N =所需随机数的数量)
s = rng(seed, 'twister');
r = a + (b-a)*rand(N,1);
rng(s);
r1 = a + (b-a)*rand(N,1);
现在我打印
r1 - r
我没有得到零。我期望得到零,因为我已经将随机数生成器重置为我的代码第三行中的初始状态。
我的问题是我哪里出错了?