Real reproducing kernel Hilbert spaces

Jordan Bell
October 22, 2015

1 Reproducing kernels

We shall often speak about functions F:X×X→ℝ, where X is a nonempty set. For x∈X, we define Fx:X→ℝ by Fx⁢(y)=F⁢(x,y) and for y∈X we define Fy:X→ℝ by Fy⁢(x)=F⁢(x,y). F is said to be symmetric if F⁢(x,y)=F⁢(y,x) for all x,y∈X and positive-definite if for any x1,…,xn∈X and c1,…,cn∈ℝ it holds that

∑1≤i,j≤nci⁢cj⁢F⁢(xi,xj)≥0.
Lemma 1.

If F:X×X→R is symmetric and positive-definite then

F⁢(x,y)2≤F⁢(x,x)⁢F⁢(y,y),x,y∈X.
Proof.

For α,β∈ℝ define11 1 See Alain Berlinet and Christine Thomas-Agnan, Reproducing Kernel Hilbert Spaces in Probability and Statistics, p. 13, Lemma 3.

C⁢(α,β) =α2⁢F⁢(x,x)+α⁢β⁢F⁢(x,y)+β⁢α⁢F⁢(y,x)+β2⁢F⁢(y,y)
=α2⁢F⁢(x,x)+2⁢α⁢β⁢F⁢(x,y)+β2⁢F⁢(y,y),

which is ≥0. Let

P⁢(α) =C⁢(α,F⁢(x,y))
=α2⁢F⁢(x,x)+2⁢α⁢F⁢(x,y)2+F⁢(x,y)2⁢F⁢(y,y),

which is ≥0. In the case F⁢(x,x)=0, the fact that P≥0 implies that F⁢(x,y)=0. In the case F⁢(x,y)≠0, P⁢(α) is a quadratic polynomial and because P≥0 it follows that the discriminant of P is ≤0:

4⁢F⁢(x,y)4-4⋅F⁢(x,x)⋅F⁢(x,y)2⁢F⁢(y,y)≤0.

That is, F⁢(x,y)4≤F⁢(x,y)2⁢F⁢(x,x)⁢F⁢(y,y), and this implies that F⁢(x,y)2≤F⁢(x,x)⁢F⁢(y,y). ∎

A real reproducing kernel Hilbert space is a Hilbert space H contained in ℝX, where X is a nonempty set, such that for each x∈X the map Λx⁢f=f⁢(x) is continuous H→ℝ. In this note we speak always about real Hilbert spaces.

Let H⊂ℝX be a reproducing kernel Hilbert space. Because H is a Hilbert space, the Riesz representation theorem states that Φ:H→H* defined by

(Φ⁢g)⁢(f)=⟨f,g⟩H,g,f∈H

is an isometric isomorphism. Because H is a reproducing kernel Hilbert space, Λx∈H* for each x∈X and we define Tx=Φ-1⁢Λx∈H, which satisfies

f⁢(x)=Λx⁢(f)=⟨f,Tx⟩H,f∈H.

In particular, because Tx∈H, for y∈X it holds that

Tx⁢(y)=Λy⁢(Tx)=⟨Tx,Ty⟩H.

Define K:X×X→ℝ by

K⁢(x,y)=⟨Tx,Ty⟩H,

called the reproducing kernel of H. For x,y∈X,

Tx⁢(y)=⟨Tx,Ty⟩H=K⁢(x,y)=Kx⁢(y),

which means that Tx=Kx.

A reproducing kernel is symmetric and positive-definite:

K⁢(x,y)=⟨Tx,Ty⟩H=⟨Ty,Tx⟩H=K⁢(y,x)

and

∑1≤i,j≤nci⁢cj⁢K⁢(xi,xj) =∑1≤i,j≤n⟨ci⁢Txi,cj⁢Txj⟩H
=⟨∑1≤i≤nci⁢Txi,∑1≤j≤ncj⁢Txj⟩H
≥0.
Lemma 2.

If E is an orthonormal basis for a reproducing kernel Hilbert space H⊂RX with reproducing kernel K:X×X→R, then

K⁢(x,y)=∑e∈Ee⁢(x)⁢e⁢(y),x,y∈X.
Proof.

Because E is an orthonormal basis for H, Parseval’s identity tell us

⟨Tx,Ty⟩H=∑e∈E⟨Tx,e⟩⁢⟨Ty,e⟩=∑e∈E⟨e,Tx⟩⁢⟨e,Ty⟩=∑e∈Ee⁢(x)⁢e⁢(y).

∎

If H⊂ℝX is a reproducing kernel Hilbert space with reproducing kernel K:X×X→ℝ and V is a closed linear subspace of H, then V is itself a reproducing kernel Hilbert space, with some reproducing kernel G:X×X→ℝ. The following theorem expresses G in terms of K.22 2 Ward Cheney and Will Light, A Course in Approximation Theory, p. 234, Chapter 31, Theorem 4.

Theorem 3.

Let H⊂RX be a reproducing kernel Hilbert space with reproducing kernel K:X×X→R, let V be a closed linear subspace of H with reproducing kernel G:X×X→R, and let PV:H→V be the projection onto V. Then

Gx=PV⁢Kx,x∈X.
Proof.

H=V⊕V⟂, thus for f∈H there are unique g∈V,h∈V⟂ such that f=g+h, and PV⁢f=g.33 3 http://individual.utoronto.ca/jordanbell/notes/pvm.pdf Then f-PV⁢f∈V⟂. Therefore for y∈X, as Gy∈V it holds that

⟨f,Gy⟩H=⟨f-PV⁢f+PV⁢f,Gy⟩H=⟨PV⁢f,Gy⟩H=(PV⁢f)⁢(y).

In particular, for x,y∈X and f=Kx,

(PV⁢Kx)⁢(y)=⟨Kx,Gy⟩H=⟨Gy,Tx⟩H=Gy⁢(x)=G⁢(y,x)=G⁢(x,y)=Gx⁢(y),

which means that PV⁢Kx=Gx, proving the claim. ∎

The Moore-Aronszajn theorem states that if X is a nonempty set and K:X×X→ℝ is a symmetric and positive-definite function, then there is a unique reproducing kernel Hilbert space H⊂ℝX for which K is the reproducing kernel.

We now prove that given a symmetric positive-definite kernel there is a unique reproducing Hilbert space for which it is the reproducing kernel.44 4 Alain Berlinet and Christine Thomas-Agnan, Reproducing Kernel Hilbert Spaces in Probability and Statistics, p. 19, Theorem 3.

2 Sobolev spaces on [0,T]

Let f∈ℝ[0,T]. The following are equivalent:55 5 Elias M. Stein and Rami Shakarchi, Real Analysis, p. 130, Theorem 3.11.

  1. 1.

    f is absolutely continuous.

  2. 2.

    f is differentiable at almost all t∈[0,T], f′∈L1, and

    f⁢(t)=f⁢(0)+∫0tf′⁢(s)⁢𝑑s,t∈[0,T].
  3. 3.

    There is some g∈L1 such that

    f⁢(t)=f⁢(0)+∫0tg⁢(s)⁢𝑑s,t∈[0,T].

In particular, if f is absolutely continuous and f′=0 almost everywhere then ∫0tf′⁢(s)⁢𝑑s=0 and so f⁢(t)=f⁢(0) for all t∈[0,T]. That is, if f is absolutely continuous and f′=0 almost everywhere then f is constant.

Let H be the set of those absolutely continuous functions f∈ℝ[0,T] such that f⁢(0)=0 and f′∈L2. For f,g∈H define

⟨f,g⟩H=∫0Tf′⁢(s)⁢g′⁢(s)⁢𝑑s.

If ∥f∥H=0 then ∫0Tf′⁢(s)2⁢𝑑s=0, which implies that f′=0 almost everywhere and hence that f is constant, and therefore f=0. Thus ⟨⋅,⋅⟩H is indeed an inner product on H.

If fn is a Cauchy sequence in H then fn′ is a Cauchy sequence in L2 and hence converges to some g∈L2. Then the function f∈ℝ[0,T] defined by

f⁢(t)=∫0tg⁢(s)⁢𝑑s,t∈[0,T],

is absolutely continuous, f⁢(0)=0, and satisfies f′=g almost everywhere, which shows that f∈H. Then fn→f in H, which proves that H is a Hilbert space. For t∈[0,T], by the Cauchy-Schwarz inequality,

|f⁢(t)|2=|∫0tf′⁢(s)⁢𝑑s|2≤|∫0Tf′⁢(s)⁢𝑑s|2≤T⁢∫0Tf′⁢(s)2⁢𝑑s=T⁢∥f∥H2,

i.e. |Lt⁢f|≤T1/2⁢∥f∥H, which shows that Lt∈H*. Therefore H is a reproducing kernel Hilbert space.

For a∈[0,T] define ha:[0,T]→ℝ by ha⁢(s)=1[0,a]⁢(s), which belongs to L2, and define ga:[0,T]→ℝ by

ga⁢(t)=∫0tha⁢(s)⁢𝑑s=min⁡(t,a),

which belongs to H. For f∈H,

⟨f,ga⟩H=∫0Tf′⁢(s)⁢ga′⁢(s)⁢𝑑s=∫0Tf′⁢(s)⁢1[0,a]⁢(s)⁢𝑑s=∫0af′⁢(s)⁢𝑑s=f⁢(a).

This means that Ka=ga. For a,b∈[0,T],

⟨Ka,Kb⟩H=∫0Tga′⁢(s)⁢gb′⁢(s)⁢𝑑s=∫0T1[0,a]⁢(s)⁢1[0,b]⁢(s)⁢𝑑s=∫0T1[0,min⁡(a,b)]⁢(s)⁢𝑑s.

That is, the reproducing kernel of H is K:[0,T]×[0,T]→ℝ,

K⁢(a,b)=⟨Ka,Kb⟩H=min⁡(a,b).

3 Sobolev spaces on ℝ

Let λ be Lebesgue measure on ℝ. Let ℒ2⁢(λ) be the collection of Borel measurable functions f:ℝ→ℝ such that |f|2 is integrable, and let L2⁢(λ) be the Hilbert space of equivalence classes of elements of ℒ2⁢(λ) where f∼g when f=g almost everywhere, with

⟨f,g⟩L2=∫ℝf⁢g⁢𝑑λ.

Let H1⁢(ℝ) be the set of locally absolutely continuous functions f:ℝ→ℝ such that f,f′∈L2⁢(λ). This is a Hilbert space with the inner product66 6 http://individual.utoronto.ca/jordanbell/notes/sobolev1d.pdf

⟨f,g⟩H1=⟨f,g⟩L2+⟨f′,g′⟩L2.

Define K:ℝ×ℝ→ℝ by

K⁢(x,y)=12⁢exp⁡(-|x-y|),x,y∈ℝ.

Let x∈ℝ. For y<x, Kx′⁢(y)=Kx⁢(y) and for y>x, Kx′⁢(y)=-Kx⁢(y), which shows that Kx∈H1⁢(ℝ). For f∈H1⁢(ℝ), doing integration by parts,

⟨f,Kx⟩H1 =∫-∞∞f⁢Kx⁢𝑑λ+∫-∞xf′⁢(y)⁢Kx⁢(y)⁢𝑑λ⁢(y)-∫x∞f′⁢(y)⁢Kx⁢(y)⁢𝑑λ⁢(y)
=∫-∞∞f⁢Kx⁢𝑑λ+f⁢(x)⁢Kx⁢(x)-∫-∞xf⁢(y)⁢Kx′⁢(y)⁢𝑑λ⁢(y)
+f⁢(x)⁢Kx⁢(x)+∫x∞f⁢(y)⁢Kx′⁢(y)⁢𝑑λ⁢(y)
=∫-∞∞f⁢Kx⁢𝑑λ+12⁢f⁢(x)-∫-∞xf⁢(y)⁢Kx⁢(y)⁢𝑑λ⁢(y)
+12⁢f⁢(x)-∫x∞f⁢(y)⁢Kx⁢(y)⁢𝑑λ⁢(y)
=f⁢(x)
=Tx⁢f.

This shows that H1⁢(ℝ) is a reproducing kernel Hilbert space. We calculate, for x<y,

⟨Tx,Ty⟩H1 =∫-∞xKx⁢Ky⁢𝑑λ+∫xyKx⁢Ky⁢𝑑λ+∫y∞Kx⁢Ky⁢𝑑λ
+∫-∞xKx⁢Ky⁢𝑑λ-∫xyKx⁢Ky⁢𝑑λ+∫y∞Kx⁢Ky⁢𝑑λ
=4⋅18⁢exp⁡(x-y)
=K⁢(x,y).

This shows that K⁢(x,y)=12⁢exp⁡(-|x-y|) is the reproducing kernel of H1⁢(ℝ).77 7 cf. Alain Berlinet and Christine Thomas-Agnan, Reproducing Kernel Hilbert Spaces in Probability and Statistics, pp. 8–9, Example 5.