Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
NCTM or The Math Forum.



Re: Challenge of the week: random numbers
Posted:
Jun 30, 2014 8:08 PM


On 06/30/2014 06:37 PM, Herman Rubin wrote: > On 20140629, Vincent Granville <datashaping@gmail.com> wrote: > >> Most random number generators use an algorithm a(k+1) = f(a(k)) > to produce a sequence of integers a(1), a(2), etc. that behaves like > random numbers. The function f is integervalued and bounded; because > of these two conditions, the sequence a(k) eventually becomes periodic > for k large enough. This is an undesirable property, and many public > random number generators (those built in Excel, Python, and other > languages) are poor and not suitable for cryptographic applications, > Markov Chains MonteCarlo associated with hierarchical Bayesian models, > or largescale MonteCarlo simulations to detect extreme events (example: > fraud detection, big data context). > >> This challenge involves identifying natural transcendental numbers with > expansions (series, continued fractions) converging very fast to easily > produce the decimals of the numbers in question. It starts with Log 10  > 2 Log 3 an a series that produces one new digit for each new term. > >> In this challenge, you will learn how to build nonperiodic, > highquality random number generators, test the strength of your random > number generator, and design generators that are easy to implement, fit > for big data and extreme events simulation, and run fast. Dr. Vincent > Granville has promised to participate, so this should add more fun to > this challenge. > > The above condition, combined with the standard mathematical definition > of a random number, is impossible. If one can compute the number ti a > specified accuracy in a finite number of steps, one can easily give > a stopping rule so that the digits at the stopped position are not > asymptotically uniform. But what does asymptotic tell one about the > first 10^10 anyhow? > > I mistrust all computer genrated sequences as random; von Neumann > stated that to even think of using systematic numbers as random > puts one in a state of sin. For SOME purposes of computing integrals > over the unit kcell, there are what are called quasirandom numbers, > which will, at least asymptotically, do better than random numbers > for sufficiently smooth functions. > > What I have suggested is that one take a "good" pseudorandom sequence, > possibly using some of Mr. Granville's suggestions, as a binary sequence > and exclusive OR it with a reasonable random sequence. The pseudorandom > part will guarantee the uniformity and more, and the physical random numbers > will break up any periodicities of the pseudorandom ones. > > >> Here's the URL: > http://www.analyticbridge.com/forum/topics/challengeoftheweekrandomnumbers > >
I take your idea seriously:
mix: (a) "strong" pseudorandom bit generator (deterministic) with: (b) a "highentropy" source of (what else?) physically generated "random bits" . (nondeterministic).
via XOR.
 http://meditationatae.wordpress.com/



