Equidistributed sequence

From HandWiki

In mathematics, a sequence (s1, s2, s3, ...) of real numbers is said to be equidistributed, or uniformly distributed, if the proportion of terms falling in a subinterval is proportional to the length of that subinterval. Such sequences are studied in Diophantine approximation theory and have applications to Monte Carlo integration.

Definition

A sequence (s1, s2, s3, ...) of real numbers is said to be equidistributed on a non-degenerate interval [a, b] if for every subinterval [c, d] of [a, b] we have

[math]\displaystyle{ \lim_{n\to\infty}{ \left|\{\,s_1,\dots,s_n \,\} \cap [c,d] \right| \over n}={d-c \over b-a} . }[/math]

(Here, the notation |{s1,...,sn} ∩ [c, d]| denotes the number of elements, out of the first n elements of the sequence, that are between c and d.)

For example, if a sequence is equidistributed in [0, 2], since the interval [0.5, 0.9] occupies 1/5 of the length of the interval [0, 2], as n becomes large, the proportion of the first n members of the sequence which fall between 0.5 and 0.9 must approach 1/5. Loosely speaking, one could say that each member of the sequence is equally likely to fall anywhere in its range. However, this is not to say that (sn) is a sequence of random variables; rather, it is a determinate sequence of real numbers.

Discrepancy

We define the discrepancy DN for a sequence (s1, s2, s3, ...) with respect to the interval [ab] as

[math]\displaystyle{ D_N = \sup_{a \le c \le d \le b} \left\vert \frac{\left|\{\,s_1,\dots,s_N \,\} \cap [c,d] \right|}{N} - \frac{d-c}{b-a} \right\vert . }[/math]

A sequence is thus equidistributed if the discrepancy DN tends to zero as N tends to infinity.

Equidistribution is a rather weak criterion to express the fact that a sequence fills the segment leaving no gaps. For example, the drawings of a random variable uniform over a segment will be equidistributed in the segment, but there will be large gaps compared to a sequence which first enumerates multiples of ε in the segment, for some small ε, in an appropriately chosen way, and then continues to do this for smaller and smaller values of ε. For stronger criteria and for constructions of sequences that are more evenly distributed, see low-discrepancy sequence.

Riemann integral criterion for equidistribution

Recall that if f is a function having a Riemann integral in the interval [a, b], then its integral is the limit of Riemann sums taken by sampling the function f in a set of points chosen from a fine partition of the interval. Therefore, if some sequence is equidistributed in [a, b], it is expected that this sequence can be used to calculate the integral of a Riemann-integrable function. This leads to the following criterion[1] for an equidistributed sequence:

Suppose (s1, s2, s3, ...) is a sequence contained in the interval [a, b]. Then the following conditions are equivalent:

  1. The sequence is equidistributed on [a, b].
  2. For every Riemann-integrable (complex-valued) function f : [a, b] → [math]\displaystyle{ \mathbb{C} }[/math], the following limit holds:
[math]\displaystyle{ \lim_{N \to \infty}\frac{1}{N}\sum_{n=1}^{N} f\left(s_n\right) = \frac{1}{b-a}\int_a^b f(x)\,dx }[/math]

This criterion leads to the idea of Monte-Carlo integration, where integrals are computed by sampling the function over a sequence of random variables equidistributed in the interval.

It is not possible to generalize the integral criterion to a class of functions bigger than just the Riemann-integrable ones. For example, if the Lebesgue integral is considered and f is taken to be in L1, then this criterion fails. As a counterexample, take f to be the indicator function of some equidistributed sequence. Then in the criterion, the left hand side is always 1, whereas the right hand side is zero, because the sequence is countable, so f is zero almost everywhere.

In fact, the de Bruijn–Post Theorem states the converse of the above criterion: If f is a function such that the criterion above holds for any equidistributed sequence in [a, b], then f is Riemann-integrable in [a, b].[2]

Equidistribution modulo 1

A sequence (a1, a2, a3, ...) of real numbers is said to be equidistributed modulo 1 or uniformly distributed modulo 1 if the sequence of the fractional parts of an, denoted by (an) or by an − ⌊an⌋, is equidistributed in the interval [0, 1].

Examples

Illustration of the filling of the unit interval (x-axis) using the first n terms of the Van der Corput sequence, for n from 0 to 999 (y-axis). Gradation in colour is due to aliasing.
0, α, 2α, 3α, 4α, ...
is equidistributed modulo 1.[3]
  • More generally, if p is a polynomial with at least one coefficient other than the constant term irrational then the sequence p(n) is uniformly distributed modulo 1.

This was proven by Weyl and is an application of van der Corput's difference theorem.[4]

  • The sequence log(n) is not uniformly distributed modulo 1.[3] This fact is related to Benford's law.
  • The sequence of all multiples of an irrational α by successive prime numbers,
2α, 3α, 5α, 7α, 11α, ...
is equidistributed modulo 1. This is a famous theorem of analytic number theory, published by I. M. Vinogradov in 1948.[5]

Weyl's criterion

Weyl's criterion states that the sequence an is equidistributed modulo 1 if and only if for all non-zero integers ℓ,

[math]\displaystyle{ \lim_{n\to\infty}\frac{1}{n}\sum_{j=1}^{n}e^{2\pi i \ell a_j}=0. }[/math]

The criterion is named after, and was first formulated by, Hermann Weyl.[7] It allows equidistribution questions to be reduced to bounds on exponential sums, a fundamental and general method.

Generalizations

  • A quantitative form of Weyl's criterion is given by the Erdős–Turán inequality.
  • Weyl's criterion extends naturally to higher dimensions, assuming the natural generalization of the definition of equidistribution modulo 1:

The sequence vn of vectors in Rk is equidistributed modulo 1 if and only if for any non-zero vector ℓ ∈ Zk,

[math]\displaystyle{ \lim_{n\to\infty}\frac{1}{n}\sum_{j=0}^{n-1}e^{2\pi i \ell \cdot v_j}=0. }[/math]

Example of usage

Weyl's criterion can be used to easily prove the equidistribution theorem, stating that the sequence of multiples 0, α, 2α, 3α, ... of some real number α is equidistributed modulo 1 if and only if α is irrational.[3]

Suppose α is irrational and denote our sequence by aj =  (where j starts from 0, to simplify the formula later). Let  ≠ 0 be an integer. Since α is irrational, ℓα can never be an integer, so [math]\displaystyle{ e^{2\pi i \ell \alpha} }[/math] can never be 1. Using the formula for the sum of a finite geometric series,

[math]\displaystyle{ \left|\sum_{j=0}^{n-1}e^{2\pi i \ell j \alpha}\right| = \left|\sum_{j=0}^{n-1}\left(e^{2\pi i \ell \alpha}\right)^j\right| = \left| \frac{1 - e^{2\pi i \ell n \alpha}} {1 - e^{2\pi i \ell \alpha}}\right| \le \frac 2 { \left|1 - e^{2\pi i \ell \alpha}\right|}, }[/math]

a finite bound that does not depend on n. Therefore, after dividing by n and letting n tend to infinity, the left hand side tends to zero, and Weyl's criterion is satisfied.

Conversely, notice that if α is rational then this sequence is not equidistributed modulo 1, because there are only a finite number of options for the fractional part of aj = .

Complete uniform distribution

A sequence [math]\displaystyle{ (a_1,a_2,\dots) }[/math] of real numbers is said to be k-uniformly distributed mod 1 if not only the sequence of fractional parts [math]\displaystyle{ a_n':=a_n-[a_n] }[/math] is uniformly distributed in [math]\displaystyle{ [0,1] }[/math] but also the sequence [math]\displaystyle{ (b_1,b_2,\dots) }[/math], where [math]\displaystyle{ b_n }[/math] is defined as [math]\displaystyle{ b_n:= (a'_{n+1}, \dots, a'_{n+k})\in [0,1]^k }[/math], is uniformly distributed in [math]\displaystyle{ [0,1]^k }[/math].

A sequence [math]\displaystyle{ (a_1,a_2,\dots) }[/math] of real numbers is said to be completely uniformly distributed mod 1 it is [math]\displaystyle{ k }[/math]-uniformly distributed for each natural number [math]\displaystyle{ k\ge 1 }[/math].

For example, the sequence [math]\displaystyle{ (\alpha, 2\alpha,\dots) }[/math] is uniformly distributed mod 1 (or 1-uniformly distributed) for any irrational number [math]\displaystyle{ \alpha }[/math], but is never even 2-uniformly distributed. In contrast, the sequence [math]\displaystyle{ (\alpha, \alpha^2,\alpha^3,\dots) }[/math] is completely uniformly distributed for almost all [math]\displaystyle{ \alpha\gt 1 }[/math] (i.e., for all [math]\displaystyle{ \alpha }[/math] except for a set of measure 0).

van der Corput's difference theorem

A theorem of Johannes van der Corput[8] states that if for each h the sequence sn+hsn is uniformly distributed modulo 1, then so is sn.[9][10][11]

A van der Corput set is a set H of integers such that if for each h in H the sequence sn+hsn is uniformly distributed modulo 1, then so is sn.[10][11]

Metric theorems

Metric theorems describe the behaviour of a parametrised sequence for almost all values of some parameter α: that is, for values of α not lying in some exceptional set of Lebesgue measure zero.

  • For any sequence of distinct integers bn, the sequence (bnα) is equidistributed mod 1 for almost all values of α.[12]
  • The sequence (αn) is equidistributed mod 1 for almost all values of α > 1.[13]

It is not known whether the sequences (en) or (πn) are equidistributed mod 1. However it is known that the sequence (αn) is not equidistributed mod 1 if α is a PV number.

Well-distributed sequence

A sequence (s1, s2, s3, ...) of real numbers is said to be well-distributed on [a, b] if for any subinterval [c, d] of [a, b] we have

[math]\displaystyle{ \lim_{n\to\infty}{ \left|\{\,s_{k+1},\dots,s_{k+n} \,\} \cap [c,d] \right| \over n}={d-c \over b-a} }[/math]

uniformly in k. Clearly every well-distributed sequence is uniformly distributed, but the converse does not hold. The definition of well-distributed modulo 1 is analogous.

Sequences equidistributed with respect to an arbitrary measure

For an arbitrary probability measure space [math]\displaystyle{ (X,\mu) }[/math], a sequence of points [math]\displaystyle{ (x_n) }[/math] is said to be equidistributed with respect to [math]\displaystyle{ \mu }[/math] if the mean of point measures converges weakly to [math]\displaystyle{ \mu }[/math]:[14]

[math]\displaystyle{ \frac{\sum_{k=1}^n \delta_{x_k}}{n}\Rightarrow \mu \ . }[/math]

In any Borel probability measure on a separable, metrizable space, there exists an equidistributed sequence with respect to the measure; indeed, this follows immediately from the fact that such a space is standard.

The general phenomenon of equidistribution comes up a lot for dynamical systems associated with Lie groups, for example in Margulis' solution to the Oppenheim conjecture.

See also

Notes

  1. Kuipers & Niederreiter (2006) pp. 2–3
  2. http://math.uga.edu/~pete/udnotes.pdf, Theorem 8
  3. 3.0 3.1 3.2 Kuipers & Niederreiter (2006) p. 8
  4. Kuipers & Niederreiter (2006) p. 27
  5. Kuipers & Niederreiter (2006) p. 129
  6. Kuipers & Niederreiter (2006) p. 127
  7. Weyl, H. (September 1916). "Über die Gleichverteilung von Zahlen mod. Eins" (in de). Math. Ann. 77 (3): 313–352. doi:10.1007/BF01475864. http://www.fuchs-braun.com/media/3ed54b58b68a224cffff80dffffffff1.pdf. 
  8. van der Corput, J. (1931), "Diophantische Ungleichungen. I. Zur Gleichverteilung Modulo Eins", Acta Mathematica (Springer Netherlands) 56: 373–456, doi:10.1007/BF02545780, ISSN 0001-5962 
  9. Kuipers & Niederreiter (2006) p. 26
  10. 10.0 10.1 Montgomery (1994) p. 18
  11. 11.0 11.1 Montgomery, Hugh L. (2001). "Harmonic analysis as found in analytic number theory". in Byrnes, James S.. Twentieth century harmonic analysis–a celebration. Proceedings of the NATO Advanced Study Institute, Il Ciocco, Italy, July 2–15, 2000. NATO Sci. Ser. II, Math. Phys. Chem.. 33. Dordrecht: Kluwer Academic Publishers. pp. 271–293. doi:10.1007/978-94-010-0662-0_13. ISBN 978-0-7923-7169-4. http://www.nato-us.org/analysis2000/papers/montgomery.pdf. 
  12. See Bernstein, Felix (1911), "Über eine Anwendung der Mengenlehre auf ein aus der Theorie der säkularen Störungen herrührendes Problem", Mathematische Annalen 71 (3): 417–439, doi:10.1007/BF01456856, https://zenodo.org/record/1701508 .
  13. Koksma, J. F. (1935), "Ein mengentheoretischer Satz über die Gleichverteilung modulo Eins", Compositio Mathematica 2: 250–258, http://www.numdam.org/item?id=CM_1935__2__250_0 
  14. Kuipers & Niederreiter (2006) p. 171

References

  • Kuipers, L.; Niederreiter, H. (2006). Uniform Distribution of Sequences. Dover Publications. ISBN 0-486-45019-8. 
  • Kuipers, L.; Niederreiter, H. (1974). Uniform Distribution of Sequences. John Wiley & Sons Inc.. ISBN 0-471-51045-9. 
  • Montgomery, Hugh L. (1994). Ten lectures on the interface between analytic number theory and harmonic analysis. Regional Conference Series in Mathematics. 84. Providence, RI: American Mathematical Society. ISBN 0-8218-0737-4. 

Further reading

External links