Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Equidistributed sequence

From Wikipedia, the free encyclopedia
Type of number sequence

Inmathematics, asequence (s1,s2,s3, ...) ofreal numbers is said to beequidistributed, oruniformly distributed, if the proportion of terms falling in a subinterval is proportional to the length of that subinterval. Such sequences are studied inDiophantine approximation theory and have applications toMonte Carlo integration.

Definition

[edit]

A sequence (s1,s2,s3, ...) ofreal numbers is said to beequidistributed on a non-degenerateinterval [a,b] if for every subinterval [c,d] of [a,b] we have

limn|{s1,,sn}[c,d]|n=dcba.{\displaystyle \lim _{n\to \infty }{\left|\{\,s_{1},\dots ,s_{n}\,\}\cap [c,d]\right| \over n}={d-c \over b-a}.}

(Here, the notation |{s1,...,sn} ∩ [c,d]| denotes the number of elements, out of the firstn elements of the sequence, that are betweenc andd.)

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], asn becomes large, the proportion of the firstn 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 ofrandom variables; rather, it is a determinate sequence of real numbers.

Discrepancy

[edit]

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

DN=supacdb||{s1,,sN}[c,d]|Ndcba|.{\displaystyle D_{N}=\sup _{a\leq c\leq d\leq b}\left\vert {\frac {\left|\{\,s_{1},\dots ,s_{N}\,\}\cap [c,d]\right|}{N}}-{\frac {d-c}{b-a}}\right\vert .}

A sequence is thus equidistributed if the discrepancyDN tends to zero asN 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, seelow-discrepancy sequence.

Riemann integral criterion for equidistribution

[edit]

Recall that iff is afunction having aRiemann integral in the interval [a,b], then its integral is the limit ofRiemann sums taken by sampling the functionf in aset 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) functionf : [a,b] →C{\displaystyle \mathbb {C} }, the following limit holds:
limN1Nn=1Nf(sn)=1baabf(x)dx{\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}
Proof
First note that the definition of an equidistributed sequence is equivalent to the integral criterion wheneverf is theindicator function of an interval: Iff =1[c,d], then the left hand side is the proportion of points of the sequence falling in the interval [c,d], and the right hand side is exactlydcba.{\displaystyle \textstyle {\frac {d-c}{b-a}}.}

This means 2 ⇒ 1 (since indicator functions are Riemann-integrable), and 1 ⇒ 2 forf being an indicator function of an interval. It remains to assume that the integral criterion holds for indicator functions and prove that it holds for general Riemann-integrable functions as well.

Note that both sides of the integral criterion equation arelinear inf, and therefore the criterion holds forlinear combinations of interval indicators, that is,step functions.

To show it holds forf being a general Riemann-integrable function, first assumef is real-valued. Then by usingDarboux's definition of the integral, we have for everyε > 0 two step functionsf1 andf2 such thatf1 ≤ f ≤ f2 andab(f2(x)f1(x))dxε(ba).{\displaystyle \textstyle \int _{a}^{b}(f_{2}(x)-f_{1}(x))\,dx\leq \varepsilon (b-a).} Notice that:

1baabf1(x)dx=limN1Nn=1Nf1(sn)lim infN1Nn=1Nf(sn){\displaystyle {\frac {1}{b-a}}\int _{a}^{b}f_{1}(x)\,dx=\lim _{N\to \infty }{\frac {1}{N}}\sum _{n=1}^{N}f_{1}(s_{n})\leq \liminf _{N\to \infty }{\frac {1}{N}}\sum _{n=1}^{N}f(s_{n})}
1baabf2(x)dx=limN1Nn=1Nf2(sn)lim supN1Nn=1Nf(sn){\displaystyle {\frac {1}{b-a}}\int _{a}^{b}f_{2}(x)\,dx=\lim _{N\to \infty }{\frac {1}{N}}\sum _{n=1}^{N}f_{2}(s_{n})\geq \limsup _{N\to \infty }{\frac {1}{N}}\sum _{n=1}^{N}f(s_{n})}

By subtracting, we see that thelimit superior and limit inferior of1Nn=1Nf(sn){\displaystyle \textstyle {\frac {1}{N}}\sum _{n=1}^{N}f(s_{n})} differ by at most ε. Since ε is arbitrary, we have the existence of the limit, and by Darboux's definition of the integral, it is the correct limit.

Finally, for complex-valued Riemann-integrable functions, the result follows again from linearity, and from the fact that every such function can be written asf = u +vi, whereu,v are real-valued and Riemann-integrable. 

This criterion leads to the idea ofMonte-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 theLebesgue integral is considered andf is taken to be inL1, then this criterion fails. As acounterexample, takef to be theindicator 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 iscountable, sof is zeroalmost everywhere.

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

Equidistribution modulo 1

[edit]

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

Examples

[edit]
Illustration of the filling of the unit interval (x-axis) using the firstn terms of the Van der Corput sequence, forn from 0 to 999 (y-axis). Gradation in colour is due to aliasing.
0,α, 2α, 3α, 4α, ...
is equidistributed modulo 1.[3]
  • More generally, ifp is apolynomial with at least one coefficient other than theconstant term irrational then the sequencep(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) isnot uniformly distributed modulo 1.[3] This fact is related toBenford's law.
  • The sequence of all multiples of an irrationalα by successiveprime numbers,
2α, 3α, 5α, 7α, 11α, ...
is equidistributed modulo 1. This is a famous theorem ofanalytic number theory, published byI. M. Vinogradov in 1948.[5]

Weyl's criterion

[edit]

Weyl's criterion states that the sequencean is equidistributed modulo 1if and only if for all non-zerointegers ℓ,

limn1nj=1ne2πiaj=0.{\displaystyle \lim _{n\to \infty }{\frac {1}{n}}\sum _{j=1}^{n}e^{2\pi i\ell a_{j}}=0.}

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

Sketch of proof
If the sequence is equidistributed modulo 1, then we can apply the Riemann integral criterion (described above) on the functionf(x)=e2πix,{\displaystyle \textstyle f(x)=e^{2\pi i\ell x},} which has integral zero on the interval [0, 1]. This gives Weyl's criterion immediately.

Conversely, suppose Weyl's criterion holds. Then the Riemann integral criterion holds for functionsf as above, and by linearity of the criterion, it holds forf being anytrigonometric polynomial. By theStone–Weierstrass theorem and an approximation argument, this extends to anycontinuous functionf.

Finally, letf be the indicator function of an interval. It is possible to boundf from above and below by two continuous functions on the interval, whose integrals differ by an arbitrary ε. By an argument similar to the proof of the Riemann integral criterion, it is possible to extend the result to anyinterval indicator functionf, thereby proving equidistribution modulo 1 of the given sequence. 

Generalizations

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

The sequencevn of vectors inRk is equidistributed modulo 1 if and only if for any non-zero vector ℓ ∈ Zk,

limn1nj=0n1e2πivj=0.{\displaystyle \lim _{n\to \infty }{\frac {1}{n}}\sum _{j=0}^{n-1}e^{2\pi i\ell \cdot v_{j}}=0.}

Example of usage

[edit]

Weyl's criterion can be used to easily prove theequidistribution 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 byaj =  (wherej starts from 0, to simplify the formula later). Let ≠ 0 be an integer. Sinceα is irrational,ℓα can never be an integer, soe2πiα{\textstyle e^{2\pi i\ell \alpha }} can never be 1. Using the formula for the sum of a finitegeometric series,

|j=0n1e2πijα|=|j=0n1(e2πiα)j|=|1e2πinα1e2πiα|2|1e2πiα|,{\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|\leq {\frac {2}{\left|1-e^{2\pi i\ell \alpha }\right|}},}

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

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

Complete uniform distribution

[edit]

A sequence(a1,a2,){\displaystyle (a_{1},a_{2},\dots )} of real numbers is said to bek-uniformly distributed mod 1 if not only the sequence of fractional partsan:=an[an]{\displaystyle a_{n}':=a_{n}-[a_{n}]} is uniformly distributed in[0,1]{\displaystyle [0,1]} but also the sequence(b1,b2,){\displaystyle (b_{1},b_{2},\dots )}, wherebn{\displaystyle b_{n}} is defined asbn:=(an+1,,an+k)[0,1]k{\displaystyle b_{n}:=(a'_{n+1},\dots ,a'_{n+k})\in [0,1]^{k}}, is uniformly distributed in[0,1]k{\displaystyle [0,1]^{k}}.

A sequence(a1,a2,){\displaystyle (a_{1},a_{2},\dots )} of real numbers is said to becompletely uniformly distributed mod 1 it isk{\displaystyle k}-uniformly distributed for each natural numberk1{\displaystyle k\geq 1}.

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

van der Corput's difference theorem

[edit]

A theorem ofJohannes van der Corput[8] states that if for eachh the sequencesn+hsn is uniformly distributed modulo 1, then so issn.[9][10][11]

Avan der Corput set is a setH of integers such that if for eachh inH the sequencesn+hsn is uniformly distributed modulo 1, then so is sn.[10][11]

Metric theorems

[edit]

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

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

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

Well-distributed sequence

[edit]

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

limn|{sk+1,,sk+n}[c,d]|n=dcba{\displaystyle \lim _{n\to \infty }{\left|\{\,s_{k+1},\dots ,s_{k+n}\,\}\cap [c,d]\right| \over n}={d-c \over b-a}}

uniformly ink. 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

[edit]

For an arbitraryprobability measure space(X,μ){\displaystyle (X,\mu )}, a sequence of points(xn){\displaystyle (x_{n})} is said to be equidistributed with respect toμ{\displaystyle \mu } if the mean ofpoint measuresconverges weakly toμ{\displaystyle \mu }:[13]

k=1nδxknμ .{\displaystyle {\frac {\sum _{k=1}^{n}\delta _{x_{k}}}{n}}\Rightarrow \mu \ .}

In anyBorelprobability measure on aseparable,metrizable space, there exists an equidistributed sequence with respect to the measure; indeed, this follows immediately from the fact that such a space isstandard.

The general phenomenon of equidistribution comes up a lot for dynamical systems associated withLie groups, for example in Margulis' solution to theOppenheim conjecture.

See also

[edit]

Notes

[edit]
  1. ^Kuipers & Niederreiter (2006) pp. 2–3
  2. ^http://math.uga.edu/~pete/udnotes.pdf, Theorem 8
  3. ^abcKuipers & Niederreiter (2006) p. 8
  4. ^Kuipers & Niederreiter (2006) p. 27
  5. ^Kuipers & Niederreiter (2006) p. 129
  6. ^Kuipers & Niederreiter (2006) p. 127
  7. ^abWeyl, H. (September 1916)."Über die Gleichverteilung von Zahlen mod. Eins" [On the distribution of numbers modulo one](PDF).Math. Ann. (in German).77 (3):313–352.doi:10.1007/BF01475864.S2CID 123470919.
  8. ^van der Corput, J. (1931), "Diophantische Ungleichungen. I. Zur Gleichverteilung Modulo Eins",Acta Mathematica,56, Springer Netherlands:373–456,doi:10.1007/BF02545780,ISSN 0001-5962,JFM 57.0230.05,Zbl 0001.20102
  9. ^Kuipers & Niederreiter (2006) p. 26
  10. ^abMontgomery (1994) p. 18
  11. ^abMontgomery, Hugh L. (2001)."Harmonic analysis as found in analytic number theory"(PDF). In Byrnes, James S. (ed.).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. Vol. 33. Dordrecht: Kluwer Academic Publishers. pp. 271–293.doi:10.1007/978-94-010-0662-0_13.ISBN 978-0-7923-7169-4.Zbl 1001.11001.
  12. ^Koksma, J. F. (1935),"Ein mengentheoretischer Satz über die Gleichverteilung modulo Eins",Compositio Mathematica,2:250–258,JFM 61.0205.01,Zbl 0012.01401
  13. ^Kuipers & Niederreiter (2006) p. 171

References

[edit]

Further reading

[edit]

External links

[edit]
Retrieved from "https://en.wikipedia.org/w/index.php?title=Equidistributed_sequence&oldid=1310219003"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp