0
$\begingroup$

Gentlemen! I have to reproduce MATLAB random normal generator (randn) with some fixed seed.

Also, I have to reproduce it in an external script, without MATLAB itself. Has anyone done this? What random normal generator implementation is better to use?

P.S. Preferred languages are Python and C/C++, but other languages may be useful too.

Improve this question
$\endgroup$
6
  • $\begingroup$ You can save/set the state of Matlab's random number generator with rng. Using it, you can configure Matlab to use the Mersenne Twister, which has implementations in many languages. However, the MT generates a uniform distribution; I don't know the details of how Matlab converts it to normal. $\endgroup$
    – MBaz
    Commented Mar 24, 2015 at 16:13
  • $\begingroup$ Well, I know about rng(), I want to emulate MATLAB randn generator outside MATLAB. I need it for some of my tests. I know, that I can just save pre-calculated noise data, but it looks like a perversion. $\endgroup$
    – Felix
    Commented Mar 24, 2015 at 16:53
  • 2
    $\begingroup$ Already answered here. $\endgroup$
    – jojeck ♦
    Commented Mar 24, 2015 at 18:01
  • 1
    $\begingroup$ @jojek How is that the same question? $\endgroup$
    – endolith
    Commented Mar 25, 2015 at 0:12
  • $\begingroup$ do you need to reproduce randn() precisely? to do the very same thing MATLAB does with the very same seed values? $\endgroup$
    – robert bristow-johnson
    Commented May 24, 2015 at 6:28

2 Answers 2

Reset to default
1
$\begingroup$

Its hard to say without knowing exactly how Matlab's randn is processed from the RNG you're using and how Matlab uses the seed.

The most reasonable way is to save pre-calculated noise data or just deal with a different random number generator (Boost among others has a normal random number generator, and you can pull it into Matlab to use via MEX if necessary).

Improve this answer
$\endgroup$
1
  • $\begingroup$ +1 even if we know that MATLAB uses the Mersenne Twister (i think that's only one mode), we don't know exactly what they do with the seed or seeds. and we don't know how they make a Gaussian random variable out if it. they could pass the uniform RV through a non-linear mapping or they could generate a dozen uniform RVs, add them up, and then subtract 6 to get something that is close to Gaussian. there are all sorts of ways of making a uniform RV into Gaussian. $\endgroup$
    – robert bristow-johnson
    Commented May 24, 2015 at 6:27
0
$\begingroup$

python: using the numpy plugin (numpy, scipy, matplotlib will make python incredibly matlab like) has the function already defined

numpy.random.normal

I don't know if this is 'cheating' but unless this is a school project where you must explicitly write out an algorithm I would assume it's ok. I say why reinvent the wheel if you don't need to?

edit

as endolith pointed out numpy provides another predefined that is even better suited for your situation

numpy.random.randn

Improve this answer
$\endgroup$
1
  • $\begingroup$ thanks for that info. I suppose both are equivalent when you supply numpy.random.normal with mean=0 and std=1 but your suggestion would be more correct, thanks $\endgroup$
    – andrew
    Commented Mar 25, 2015 at 2:34

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.