Asymptotic equipartition property

In information theory the asymptotic equipartition property (AEP) is a general property of the output samples of a stochastic source. It is fundamental to the concept of typical set used in theories of compression.

Roughly speaking, the theorem states that although there are many series of results that may be produced by a random process, the one actually produced is most probably from a loosely-defined set of outcomes that all have approximately the same chance of being the one actually realized. (This is a consequence of the law of large numbers and ergodic theory.) Although there are individual outcomes which have a higher probability than any outcome in this set, the vast number of outcomes in the set almost guarantees that the outcome will come from the set.

In the field of Pseudorandom number generation, a candidate generator of undetermined quality whose output sequence lies too far outside the typical set by some statistical criteria is rejected as insufficiently random. Thus, although the typical set is loosely defined, practical notions arise concerning sufficient typicality.

Definition

Given a discrete-time stationary ergodic stochastic process <math>X</math> on the probability space <math>(\Omega,B,p)</math>, AEP is an assertion that

where <math>X_1^n</math> denotes the process limited to duration <math>\{1,\dots,n\}</math>, and <math>H(X)</math> or simply <math>H</math> denotes the entropy rate of <math>X</math>, which must exist for all discrete-time stationary processes including the ergodic ones. AEP is proved for finite-valued (i.e. <math>|\Omega|<\infty</math>) stationary ergodic stochastic processes in the Shannon-McMillan-Breiman theorem using the ergodic theory and for any i.i.d. sources directly using the law of large numbers in both the discrete-valued case (where <math>H</math> is simply the entropy of a symbol) and the continuous-valued case (where <math>H</math> is the differential entropy instead). The definition of AEP can also be extended for certain classes of continuous-time stochastic processes for which a typical set exists for long enough observation time. The convergence is proven almost sure in all cases.

AEP for discrete-time i.i.d. sources

Given <math>X</math> is an i.i.d. source, its time series X₁, ..., X_n is i.i.d. with entropy H(X) in the discrete-valued case and differential entropy in the continuous-valued case. The weak law of large numbers gives the AEP with convergence in probability,

<math>

\lim_{n\to\infty}\Pr\left[\left|-\frac{1}{n} \log p(X_1, X_2, ..., X_n) - H(X)\right|> \epsilon\right]=0 \qquad \forall \epsilon>0. </math> since the entropy is equal to the expectation of <math>\frac{1}{n} \log p(X_1, X_2, ..., X_n)</math>. The strong law of large number asserts the stronger almost sure convergence,

<math>

\Pr\left[\lim_{n\to\infty} - \frac{1}{n} \log p(X_1, X_2, ..., X_n) = H(X)\right]=1 </math> which implies the result from the weak law of large numbers.

AEP for discrete-time finite-valued stationary ergodic sources

Consider a finite-valued sample space <math>\Omega</math>, i.e. <math>|\Omega|<\infty</math>, for the discrete-time stationary ergodic process <math>X:=\{X_n\}</math> defined on the probability space <math>(\Omega,B,p)</math>. The AEP for such stochastic source, known as the Shannon-McMillan-Breiman theorem, can be shown using the sandwich proof by Algoet and Cover outlined as follows:

Let <math>x</math> denote some measurable set <math>x=X(A)</math> for some <math>A\in B</math>
Parameterize the joint probability by <math>n</math> and x as <math>j(n,x):=p(x_0^{n-1})</math>
Parameterize the conditional probability by <math>i,k</math> and <math>x</math> as <math>c(i,k,x) := p(x_i|x_{i-k}^{i-1})</math>.
Take the limit of the conditional probability as <math>k\to\infty</math> and denote it as <math>c(i,x):=p(x_i|x_{-\infty}^{i-1})</math>
Argue the two notions of entropy rate <math>\lim_{n\to\infty} E[-\log j(n,X)]</math> and <math>\lim_{n\to\infty} E[-\log c(n,n,X)]</math> exist and are equal for any stationary process including the stationary ergodic process <math>X</math>. Denote it as <math>H</math>.
Argue that both <math>c(i,k,X):=\{p(X_i|X_{i-k}^{i-1})\}</math> and <math>c(i,X):=\{p(X_i|X_{-\infty}^{i-1})\}</math>, where <math>i</math> is the time index, are stationary ergodic processes, whose sample means converge almost surely to some values denoted by <math>H^k</math> and <math>H^{\infty}</math> respectively.
Define the <math>k</math>-th order Markov approximation to the probability <math>a(n,k,x)</math> as

Argue that <math>a(n,k,X(\Omega))</math> is finite from the finite-value assumption.
Express <math>-\frac1n\log a(n,k,X)</math> in terms of the sample mean of <math>c(i,k,X)</math> and show that it converges almost surely to <math>H^k</math>
Define <math>a(n,x):=p(x_0^{n-1}|x_{-\infty}^{-1})</math>, which is a probability measure.
Express <math>-\frac1n\log a(n,X)</math> in terms of the sample mean of <math>c(i,X)</math> and show that it converges almost surely to <math>H^{\infty}</math>
Argue that <math>H^k\searrow H</math> as <math>k\to\infty</math> using the stationarity of the process.
Argue that <math>H=H^{\infty}</math> using the Lévy's martingale convergence theorem and the finite-value assumption.
Show that <math>E\left[\frac{a(n,k,X)}{j(n,X)}\right]=a(n, k,X(\Omega))</math> which is finite as argued before.
Show that <math>E\left[\frac{j(n,X)}{a(n,X)}\right]=1</math> by conditioning on the infinite past <math>X_{-\infty}^{-1}</math> and iterating the expectation.
Show that <math>(\forall \alpha\in\mathbb{R})\left(\Pr\left[\frac{a(n,k,X)}{j(n,X)}\geq \alpha \right]\leq \frac{a(n, k,X(\Omega))}{\alpha}\right)</math> using the Markov's inequality and the expectation derived previously.
Similarly, show that <math>(\forall \alpha\in\mathbb{R})\left(\Pr\left[\frac{j(n,X)}{a(n,X)}\geq \alpha \right]\leq \frac{1}{\alpha}\right)</math>, which is equivalent to <math>(\forall \alpha\in\mathbb{R})\left(\Pr\left[\frac1n\log\frac{j(n,X)}{a(n,X)}\geq \frac1n\log\alpha \right]\leq \frac{1}{\alpha}\right)</math>.
Show that <math>\limsup_{n\to\infty}</math> of both <math>\frac1n \log \frac{a(n,k,X)}{j(n,X)}</math> and <math>\frac1n \log\frac{j(n,X)}{a(n,X)}</math> are non-positive almost surely by setting <math>\alpha=n^{\beta}</math> for any <math>\beta>1</math> and applying the Borel-Cantelli lemma.
Show that <math>\liminf</math> and <math>\limsup</math> of <math>-\frac1n \log j(n,X)</math> are lower and upper bounded almost surely by <math>H^{\infty}</math> and <math>H^k</math> respectively by breaking up the logarithms in the previous result.
Complete the proof by pointing out that the upper and lower bounds are shown previously to approach <math>H</math> as <math>k\to\infty</math>.

AEP for non-stationary discrete-time source producing independent symbols

The assumptions of stationarity/ergodicity/idential distribution of random variables is not essential for the AEP to hold. Indeed, as is quite clear intuitively, the AEP requires only some form of the law of large numbers to hold, which is fairly general. However, the expression needs to be suitably generalized, and the conditions need to be formulated precisely.

We assume that the source is producing independent symbols, with possibly different output statistics at each instant. We assume that the statistics of the process are known completely, that is, the marginal distribution of the process seen at each time instant is known. The joint distribution is just the product of marginals. Then, under the condition (which can be relaxed) that <math>Var[Logp(X i)<M</math> for all i, for some M>0, the following holds (AEP):

<math> \lim_{n\to\infty}\Pr\left[\left|-\frac{1}{n} \log p(X_1, X_2, ..., X_n) - \bar{H_n}(X)\right|< \epsilon\right]=1\qquad \forall \epsilon>0 </math> where, <math>\bar{H_n}(X)=\frac{H(X_1,X_2,\ldots,X_n)}{n}</math>

Proof

The proof follows from a simple application of Markov's inequality (applied to second moment of <math>log(p(X_i))</math>.

<math> \Pr\left[\left|-\frac{1}{n} \log p(X_1, X_2, ..., X_n) - \bar{H}(X)\right|> \epsilon\right]\leq \frac{E[\sum_{i=1}^n [log(p(X_i)]^2]}{n^2\times \epsilon^2}\leq \frac{M}{n\times \epsilon^2} \rightarrow 0 \;\;\;\mbox{as}\;\;\; n\rightarrow\infty </math>

It is obvious that the proof holds if any moment <math>E[|logp(X i)|^r]</math> is uniformly bounded for <math>r>1</math> (again by Markov's inequality applied to rth moment). <math>\Box{} </math>

Even this condition is not necessary, but given a non-stationary random process, it should not be difficult to test whether the AEP holds using the above method.

Applications for AEP for non-stationary source producing independent symbols

The AEP for non-stationary discrete-time independent process leads us to (among other results) source coding theorem for non-stationary source (with independent output symbols) and channel coding theorem for non-stationary memoryless channels.

Source Coding Theorem

The source coding theorem for discrete time non-stationary independent sources can be found here: source coding theorem

Channel Coding Theorem

Channel coding theorem for discrete time non-stationary memoryless channels can be found here: noisy channel coding theorem

AEP for certain continuous-time stationary ergodic sources

Discrete-time functions can be interpolated to continuous-time functions. If such interpolation <math>f</math> is measurable, we may define the continuous-time stationary process accordingly as <math>\tilde{X}:=f\circ X</math>. If AEP holds for the discrete-time process, as in the i.i.d. or finite-valued stationary ergodic cases shown above, it automatically holds for the continuous-time stationary process derived from it by some measurable interpolation. i.e. <math> -\frac{1}{n} \log p(\tilde{X}_0^\tau) \to H(X) </math> where <math>n</math> corresponds to the degree of freedom in time <math>\tau</math>. <math>nH(X)/\tau</math> and <math>H(X)</math> are the entropy per unit time and per degree of freedom respectively, defined by Shannon.

An important class of such continuous-time stationary process is the bandlimited stationary ergodic process with the sample space being a subset of the continuous <math>\mathcal{L}_2</math> functions. AEP holds if the process is white, in which case the time samples are i.i.d., or there exists <math>T>1/2W</math>, where <math>W</math> is the nominal bandwidth, such that the <math>T</math>-spaced time samples take values in a finite set, in which case we have the discrete-time finite-valued stationary ergodic process.

Any time-invariant operations also preserves AEP, stationarity and ergodicity and we may easily turn a stationary process to non-stationary without losing AEP by nulling out a finite number of time samples in the process.

References

The Classic Paper

Claude E. Shannon. A Mathematical Theory of Communication. Bell System Technical Journal, July/October 1948.

Textbooks on Information Theory

Thomas M. Cover, Joy A. Thomas. Elements of Information Theory New York: Wiley, 1991. ISBN 0-471-06259-6

Template:WS