Mercer's theorem

In mathematics, specifically functional analysis, Mercer's theorem is a representation of a symmetric positive-definite function on a square as a sum of a convergent sequence of product functions. This theorem, presented in (Mercer 1909), is one of the most notable results of the work of James Mercer. It is an important theoretical tool in the theory of integral equations; it is used in the Hilbert space theory of stochastic processes, for example the Karhunen–Loève theorem; and it is also used to characterize a symmetric positive semi-definite kernel.[1]

Introduction

To explain Mercer's theorem, we first consider an important special case; see below for a more general formulation. A kernel, in this context, is a symmetric continuous function that maps

 K: [a,b] \times [a,b] \rightarrow \mathbb{R}

where symmetric means that K(x, s) = K(s, x).

K is said to be non-negative definite (or positive semidefinite) if and only if

 \sum_{i=1}^n\sum_{j=1}^n \langle K(x_i, x_j) c_i, c_j \rangle \geq 0

for all finite sequences of points x1, ..., xn of [a, b] and all choices of real numbers c1, ..., cn (cf. positive definite kernel).

Associated to K is a linear operator on functions defined by the integral

 [T_K \varphi](x) =\int_a^b  K(x,s) \varphi(s)\, ds.

For technical considerations we assume φ can range through the space L2[a, b] (see Lp space) of square-integrable real-valued functions. Since T is a linear operator, we can talk about eigenvalues and eigenfunctions of T.

Theorem. Suppose K is a continuous symmetric non-negative definite kernel. Then there is an orthonormal basis {ei}i of L2[a, b] consisting of eigenfunctions of TK such that the corresponding sequence of eigenvalues {λi}i is nonnegative. The eigenfunctions corresponding to non-zero eigenvalues are continuous on [a, b] and K has the representation

 K(s,t) = \sum_{j=1}^\infty \lambda_j \, e_j(s) \, e_j(t)

where the convergence is absolute and uniform.

Details

We now explain in greater detail the structure of the proof of Mercer's theorem, particularly how it relates to spectral theory of compact operators.

To show compactness, show that the image of the unit ball of L2[a,b] under TK equicontinuous and apply Ascoli's theorem, to show that the image of the unit ball is relatively compact in C([a,b]) with the uniform norm and a fortiori in L2[a,b].

Now apply the spectral theorem for compact operators on Hilbert spaces to TK to show the existence of the orthonormal basis {ei}i of L2[a,b]

 \lambda_i e_i(t)= [T_K e_i](t)  = \int_a^b K(t,s) e_i(s)\, ds.

If λi 0, the eigenvector ei is seen to be continuous on [a,b]. Now

 \sum_{i=1}^\infty \lambda_i |e_i(t) e_i(s)| \leq \sup_{x \in [a,b]} |K(x,x)|^2,

which shows that the sequence

 \sum_{i=1}^\infty \lambda_i e_i(t) e_i(s)

converges absolutely and uniformly to a kernel K0 which is easily seen to define the same operator as the kernel K. Hence K=K0 from which Mercer's theorem follows.

Trace

The following is immediate:

Theorem. Suppose K is a continuous symmetric non-negative definite kernel; TK has a sequence of nonnegative eigenvalues {λi}i. Then

 \int_a^b K(t,t)\, dt = \sum_i \lambda_i.

This shows that the operator TK is a trace class operator and

 \operatorname{trace}(T_K) = \int_a^b K(t,t)\, dt.

Generalizations

Mercer's theorem itself is a generalization of the result that any positive semidefinite matrix is the Gramian matrix of a set of vectors.

The first generalization replaces the interval [a, b] with any compact Hausdorff space and Lebesgue measure on [a, b] is replaced by a finite countably additive measure μ on the Borel algebra of X whose support is X. This means that μ(U) > 0 for any nonempty open subset U of X.

A recent generalization replaces this conditions by that follows: the set X is a first-countable topological space endowed with a Borel (complete) measure μ. X is the support of μ and, for all x in X, there is an open set U containing x and having finite measure. Then essentially the same result holds:

Theorem. Suppose K is a continuous symmetric positive definite kernel on X. If the function κ is L1μ(X), where κ(x)=K(x,x), for all x in X, then there is an orthonormal set {ei}i of L2μ(X) consisting of eigenfunctions of TK such that corresponding sequence of eigenvalues {λi}i is nonnegative. The eigenfunctions corresponding to non-zero eigenvalues are continuous on X and K has the representation

 K(s,t) = \sum_{j=1}^\infty \lambda_j \, e_j(s) \, e_j(t)

where the convergence is absolute and uniform on compact subsets of X.

The next generalization deals with representations of measurable kernels.

Let (X, M, μ) be a σ-finite measure space. An L2 (or square integrable) kernel on X is a function

 K \in L^2_{\mu \otimes \mu}(X \times X).

L2 kernels define a bounded operator TK by the formula

 \langle T_K \varphi, \psi \rangle = \int_{X \times X} K(y,x) \varphi(y) \psi(x) \,d[\mu \otimes \mu](y,x).

TK is a compact operator (actually it is even a Hilbert–Schmidt operator). If the kernel K is symmetric, by the spectral theorem, TK has an orthonormal basis of eigenvectors. Those eigenvectors that correspond to non-zero eigenvalues can be arranged in a sequence {ei}i (regardless of separability).

Theorem. If K is a symmetric positive definite kernel on(X, M, μ), then

 K(y,x) = \sum_{i \in \mathbb{N}} \lambda_i e_i(y) e_i(x)

where the convergence in the L2 norm. Note that when continuity of the kernel is not assumed, the expansion no longer converges uniformly.

See also

References

Notes