首先要知道一点计算机生成的随机数都是伪随机数,当然malab也不例外
而这些伪随机数是怎样生成的呢?大多是用同余来生成比如x_n+1=(25173*x_n+13849) (mod 65536),而这些伪随机数都需要一个初值x_0,这样才能迭代这个x_0就是此处的’state’参数,所以rand(‘state’,sum(clock))的作用是定义一个随时间变化的初值x_0。clock返回一个6个元素的向量分辨是年月日时分秒,sum(clock),sum加起来就作为随机数的状态,因为你每次运行程序的时间不同,所以得到的随机数序列就不同
1.设置初始状态
rand('seed', S)
rand('state', S)
rand('twister', S)
seed表示采用v4版本的随机数产生器,state是v5版本的随机数产生器,最后的twister用的则是Mersenne Twister随机数产生器
2.更新
2012及以后的版本默认采用Mersenne Twister随机数产生器rng(S)函数表示设定初始状态,rng(‘shuffle’) 表示随机分配一个初始状态
>> rng(1);
>> rand(3,1)
ans =
0.4170
0.7203
0.0001
>> rand(3,1)
ans =
0.3023
0.1468
0.0923
>> rng(1)
>> rand(3,1)
ans =
0.4170
0.7203
0.0001