Ergodic process

From HandWiki

In physics, statistics, econometrics and signal processing, a stochastic process is said to be in an ergodic regime if an observable's ensemble average equals the time average.[1] In this regime, any collection of random samples from a process must represent the average statistical properties of the entire regime. Conversely, a process that is not in ergodic regime is said to be in non-ergodic regime.[2]

Specific definitions

One can discuss the ergodicity of various statistics of a stochastic process. For example, a wide-sense stationary process [math]\displaystyle{ X(t) }[/math] has constant mean

[math]\displaystyle{ \mu_X= E[X(t)], }[/math]

and autocovariance

[math]\displaystyle{ r_X(\tau) = E[(X(t)-\mu_X) (X(t+\tau)-\mu_X)], }[/math]

that depends only on the lag [math]\displaystyle{ \tau }[/math] and not on time [math]\displaystyle{ t }[/math]. The properties [math]\displaystyle{ \mu_X }[/math] and [math]\displaystyle{ r_X(\tau) }[/math] are ensemble averages (calculated over all possible sample functions [math]\displaystyle{ X }[/math]), not time averages.

The process [math]\displaystyle{ X(t) }[/math] is said to be mean-ergodic[3] or mean-square ergodic in the first moment[4] if the time average estimate

[math]\displaystyle{ \hat{\mu}_X = \frac{1}{T} \int_{0}^{T} X(t) \, dt }[/math]

converges in squared mean to the ensemble average [math]\displaystyle{ \mu_X }[/math] as [math]\displaystyle{ T \rightarrow \infty }[/math].

Likewise, the process is said to be autocovariance-ergodic or d moment[4] if the time average estimate

[math]\displaystyle{ \hat{r}_X(\tau) = \frac{1}{T} \int_0^T [X(t+\tau)-\mu_X] [X(t)-\mu_X] \, dt }[/math]

converges in squared mean to the ensemble average [math]\displaystyle{ r_X(\tau) }[/math], as [math]\displaystyle{ T \rightarrow \infty }[/math]. A process which is ergodic in the mean and autocovariance is sometimes called ergodic in the wide sense.

Discrete-time random processes

The notion of ergodicity also applies to discrete-time random processes [math]\displaystyle{ X[n] }[/math] for integer [math]\displaystyle{ n }[/math].

A discrete-time random process [math]\displaystyle{ X[n] }[/math] is ergodic in mean if

[math]\displaystyle{ \hat{\mu}_X = \frac{1}{N} \sum_{n=1}^{N} X[n] }[/math]

converges in squared mean to the ensemble average [math]\displaystyle{ E[X] }[/math], as [math]\displaystyle{ N \rightarrow \infty }[/math].

Examples

Ergodicity means the ensemble average equals the time average. Following are examples to illustrate this principle.

Call centre

Each operator in a call centre spends time alternately speaking and listening on the telephone, as well as taking breaks between calls. Each break and each call are of different length, as are the durations of each 'burst' of speaking and listening, and indeed so is the rapidity of speech at any given moment, which could each be modelled as a random process.

  • Take N call centre operators (N should be a very large integer) and plot the number of words spoken per minute for each operator over a long period (several shifts). For each operator you will have a series of points, which could be joined with lines to create a 'waveform'.
  • Calculate the average value of those points in the waveform; this gives you the time average.
  • There are N waveforms and N operators. These N waveforms are known as an ensemble.
  • Now take a particular instant of time in all those waveforms and find the average value of the number of words spoken per minute. That gives you the ensemble average for that instant.
  • If ensemble average always equals time average, then the system is ergodic.

Electronics

Each resistor has an associated thermal noise that depends on the temperature. Take N resistors (N should be very large) and plot the voltage across those resistors for a long period. For each resistor you will have a waveform. Calculate the average value of that waveform; this gives you the time average. There are N waveforms as there are N resistors. These N plots are known as an ensemble. Now take a particular instant of time in all those plots and find the average value of the voltage. That gives you the ensemble average for each plot. If ensemble average and time average are the same then it is ergodic.

Examples of non-ergodic random processes

  • An unbiased random walk is non-ergodic. Its expectation value is zero at all times, whereas its time average is a random variable with divergent variance.
  • Suppose that we have two coins: one coin is fair and the other has two heads. We choose (at random) one of the coins first, and then perform a sequence of independent tosses of our selected coin. Let X[n] denote the outcome of the nth toss, with 1 for heads and 0 for tails. Then the ensemble average is ​12  (​12 +  1) = ​34; yet the long-term average is ​12 for the fair coin and 1 for the two-headed coin. So the long term time-average is either 1/2 or 1. Hence, this random process is not ergodic in mean.

See also

Notes

  1. Cherstvy, Andrey; Chechkin, Aleksei V; Metzler, Ralf (2013), "Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes", New J. Phys. 15: 083039, doi:10.1088/1367-2630/15/8/083039 
  2. Originally due to L. Boltzmann. See part 2 of Vorlesungen über Gastheorie. Leipzig: J. A. Barth. 1898. OCLC 01712811. https://archive.org/details/vorlesungenberg02boltgoog.  ('Ergoden' on p. 89 in the 1923 reprint.) It was used to prove equipartition of energy in the kinetic theory of gases
  3. Papoulis, p. 428
  4. 4.0 4.1 Porat, p. 14

References

  • Porat, B. (1994). Digital Processing of Random Signals: Theory & Methods. Prentice Hall. pp. 14. ISBN 0-13-063751-3. 
  • Papoulis, Athanasios (1991). Probability, random variables, and stochastic processes. New York: McGraw-Hill. pp. 427–442. ISBN 0-07-048477-5.