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 xX, we define Fx:X by Fx(y)=F(x,y) and for yX 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,yX and positive-definite if for any x1,,xnX and c1,,cn it holds that

1i,jncicjF(xi,xj)0.
Lemma 1.

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

F(x,y)2F(x,x)F(y,y),x,yX.
Proof.

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

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

which is 0. Let

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

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

4F(x,y)4-4F(x,x)F(x,y)2F(y,y)0.

That is, F(x,y)4F(x,y)2F(x,x)F(y,y), and this implies that F(x,y)2F(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 xX the map Λxf=f(x) is continuous H. In this note we speak always about real Hilbert spaces.

Let HX be a reproducing kernel Hilbert space. Because H is a Hilbert space, the Riesz representation theorem states that Φ:HH* defined by

(Φg)(f)=f,gH,g,fH

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

f(x)=Λx(f)=f,TxH,fH.

In particular, because TxH, for yX it holds that

Tx(y)=Λy(Tx)=Tx,TyH.

Define K:X×X by

K(x,y)=Tx,TyH,

called the reproducing kernel of H. For x,yX,

Tx(y)=Tx,TyH=K(x,y)=Kx(y),

which means that Tx=Kx.

A reproducing kernel is symmetric and positive-definite:

K(x,y)=Tx,TyH=Ty,TxH=K(y,x)

and

1i,jncicjK(xi,xj) =1i,jnciTxi,cjTxjH
=1inciTxi,1jncjTxjH
0.
Lemma 2.

If E is an orthonormal basis for a reproducing kernel Hilbert space HRX with reproducing kernel K:X×XR, then

K(x,y)=eEe(x)e(y),x,yX.
Proof.

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

Tx,TyH=eETx,eTy,e=eEe,Txe,Ty=eEe(x)e(y).

If HX 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 HRX be a reproducing kernel Hilbert space with reproducing kernel K:X×XR, let V be a closed linear subspace of H with reproducing kernel G:X×XR, and let PV:HV be the projection onto V. Then

Gx=PVKx,xX.
Proof.

H=VV, thus for fH there are unique gV,hV such that f=g+h, and PVf=g.33 3 http://individual.utoronto.ca/jordanbell/notes/pvm.pdf Then f-PVfV. Therefore for yX, as GyV it holds that

f,GyH=f-PVf+PVf,GyH=PVf,GyH=(PVf)(y).

In particular, for x,yX and f=Kx,

(PVKx)(y)=Kx,GyH=Gy,TxH=Gy(x)=G(y,x)=G(x,y)=Gx(y),

which means that PVKx=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 HX 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], fL1, and

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

    There is some gL1 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 fL2. For f,gH define

f,gH=0Tf(s)g(s)𝑑s.

If fH=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 gL2. 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 fH. Then fnf 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|2T0Tf(s)2𝑑s=TfH2,

i.e. |Ltf|T1/2fH, which shows that LtH*. 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 fH,

f,gaH=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,KbH=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,KbH=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 fg when f=g almost everywhere, with

f,gL2=fg𝑑λ.

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

f,gH1=f,gL2+f,gL2.

Define K:× by

K(x,y)=12exp(-|x-y|),x,y.

Let x. For y<x, Kx(y)=Kx(y) and for y>x, Kx(y)=-Kx(y), which shows that KxH1(). For fH1(), doing integration by parts,

f,KxH1 =-fKx𝑑λ+-xf(y)Kx(y)𝑑λ(y)-xf(y)Kx(y)𝑑λ(y)
=-fKx𝑑λ+f(x)Kx(x)--xf(y)Kx(y)𝑑λ(y)
+f(x)Kx(x)+xf(y)Kx(y)𝑑λ(y)
=-fKx𝑑λ+12f(x)--xf(y)Kx(y)𝑑λ(y)
+12f(x)-xf(y)Kx(y)𝑑λ(y)
=f(x)
=Txf.

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

Tx,TyH1 =-xKxKy𝑑λ+xyKxKy𝑑λ+yKxKy𝑑λ
+-xKxKy𝑑λ-xyKxKy𝑑λ+yKxKy𝑑λ
=418exp(x-y)
=K(x,y).

This shows that K(x,y)=12exp(-|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.