How to set custom seed for pseudo-random number generator

I need to perform few tests where I use randn pseudo random number generator. How can I set the seed on my own, so every time I run this test I will get the same results? (yeah, I know it's a little bit weird, but that's the problem).

I've found the RANDSTREAM object that has the seed property, but it's read only. Is there any way to use it for seeding the generator?


Solution 1:

The old way of doing it:

randn('seed',0)

The new way:

s = RandStream('mcg16807','Seed',0)
RandStream.setDefaultStream(s)

Note that if you use the new way, rand and randn share the same stream so if you are calling both, you may find different numbers being generated compared to the old method (which has separate generators). The old method is still supported for this reason (and legacy code).

See http://www.mathworks.com/help/techdoc/math/bsn94u0-1.html for more info.

Solution 2:

You can just call rng(mySeed) to set the seed for the global stream (tested in Matlab R2011b). This affects the rand, randn, and randi functions.

The same page that James linked to lists this as the recommended alternative to various old methods (see the middle cell of the right column of the table).

Here's some example code:

format long;             % Display numbers with full precision
format compact;          % Get rid of blank lines between output
mySeed = 10;
rng(mySeed);             % Set the seed
disp(rand([1,3]));
disp(randi(10,[1,10]));
disp(randn([1,3]));
disp(' ');
rng(mySeed);             % Set the seed again to duplicate the results
disp(rand([1,3]));
disp(randi(10,[1,10]));
disp(randn([1,3]));

Its output is:

   0.771320643266746   0.020751949359402   0.633648234926275
     8     5     3     2     8     2     1     7    10     1
   0.060379730526407   0.622213879877005   0.109700311365407

   0.771320643266746   0.020751949359402   0.633648234926275
     8     5     3     2     8     2     1     7    10     1
   0.060379730526407   0.622213879877005   0.109700311365407

Solution 3:

mySeed=57; % an integer number
rng(mySeed,'twister') %You can replace 'twister' with other generators