Reproducing kernel Hilbert space explained

In functional analysis (a branch of mathematics), a reproducing kernel Hilbert space (RKHS) is a Hilbert space of functions in which point evaluation is a continuous linear functional. Roughly speaking, this means that if two functions

f and g in the RKHS are close in norm, i.e., \|f - g\| is small, then f and g are also pointwise close, i.e., |f(x) - g(x)| is small for all x . The converse does not need to be true. Informally, this can be shown by looking at the supremum norm: the sequence of functions \sin ^{2n} (x) converges pointwise, but does not converge uniformly i.e. does not converge with respect to the supremum norm. (This is not a counterexample because the supremum norm does not arise from any inner product due to not satisfying the parallelogram law.)

It is not entirely straightforward to construct a Hilbert space of functions which is not an RKHS.^[1] Some examples, however, have been found.^{[2] [3]}

L² spaces are not Hilbert spaces of functions (and hence not RKHSs), but rather Hilbert spaces of equivalence classes of functions (for example, the functions

f and g defined by f(x) =0 and g(x) =1 _{Q} are equivalent in L²). However, there are RKHSs in which the norm is an L²-norm, such as the space of band-limited functions (see the example below).

An RKHS is associated with a kernel that reproduces every function in the space in the sense that for every

x in the set on which the functions are defined, "evaluation at x " can be performed by taking an inner product with a function determined by the kernel. Such a reproducing kernel exists if and only if every evaluation functional is continuous.

The reproducing kernel was first introduced in the 1907 work of Stanisław Zaremba concerning boundary value problems for harmonic and biharmonic functions. James Mercer simultaneously examined functions which satisfy the reproducing property in the theory of integral equations. The idea of the reproducing kernel remained untouched for nearly twenty years until it appeared in the dissertations of Gábor Szegő, Stefan Bergman, and Salomon Bochner. The subject was eventually systematically developed in the early 1950s by Nachman Aronszajn and Stefan Bergman.^[4]

These spaces have wide applications, including complex analysis, harmonic analysis, and quantum mechanics. Reproducing kernel Hilbert spaces are particularly important in the field of statistical learning theory because of the celebrated representer theorem which states that every function in an RKHS that minimises an empirical risk functional can be written as a linear combination of the kernel function evaluated at the training points. This is a practically useful result as it effectively simplifies the empirical risk minimization problem from an infinite dimensional to a finite dimensional optimization problem.

For ease of understanding, we provide the framework for real-valued Hilbert spaces. The theory can be easily extended to spaces of complex-valued functions and hence include the many important examples of reproducing kernel Hilbert spaces that are spaces of analytic functions.^[5]

Definition

Let

X be an arbitrary set and H a Hilbert space of real-valued functions on X , equipped with pointwise addition and pointwise scalar multiplication. The evaluation functional over the Hilbert space of functions H is a linear functional that evaluates each function at a point x ,L _{x} : f \mapsto f(x) \forall f \in H.

We say that H is a reproducing kernel Hilbert space if, for all

x in X , L _{x} is continuous at every f in H or, equivalently, if L _{x} is a bounded operator on H , i.e. there exists some M _{x>0} such that

Although

M _{x<infty} is assumed for all x \in X , it might still be the case that $\backslash sup\_x\; M\_x\; =\; \backslash infty$.

While property is the weakest condition that ensures both the existence of an inner product and the evaluation of every function in

H at every point in the domain, it does not lend itself to easy application in practice. A more intuitive definition of the RKHS can be obtained by observing that this property guarantees that the evaluation functional can be represented by taking the inner product of f with a function K _{x} in H . This function is the so-called reproducing kernel for the Hilbert space H from which the RKHS takes its name. More formally, the Riesz representation theorem implies that for all x in X there exists a unique element K _{x} of H with the reproducing property,

Since

K _{x} is itself a function defined on X with values in the field R (or C in the case of complex Hilbert spaces) and as K _{x} is in H we have thatK _{x(y)} = L _{y(K} _{x)=} \langle K _{x, K} _{y} \rangle _{H,} where K _{y\in} H is the element in H associated to L _{y} .

This allows us to define the reproducing kernel of

H as a function K: X x X \to R (or C in the complex case) byK(x,y) = \langle K _{x, K} _{y} \rangle _{H.}

From this definition it is easy to see that

K: X x X \to R (or C in the complex case) is both symmetric (resp. conjugate symmetric) and positive definite, i.e.\sum ^{n} _{i,j =1} c _{i} c _{j} K(x _{i,} x _{j)= \sum} ^{n} ^{i =1} c _{i} \left\langle K _{x _{i}} , \sum ^{n} _{j =1} c _{j} K _{x _{j}} \right\rangle _{H} = \left\langle \sum ^{n} _{i =1} c _{i} K _{x _{i}} , \sum ^{n} _{j =1} c _{j} K _{x _{j}} \right\rangle _{H} = \left\|\sum ^{nc} _{iK} ^{x _{i}} \right\| ^{2} _{H} \ge 0

for every

n \in N , x _{1,} ... , x _{n} \in X, and c _{1,} ... , c _{n} \in R . ^[6] The Moore–Aronszajn theorem (see below) is a sort of converse to this: if a function K satisfies these conditions then there is a Hilbert space of functions on X for which it is a reproducing kernel.

Example

The space of bandlimited continuous functions

H is a RKHS, as we now show. Formally, fix some cutoff frequency 0 <a < inf ty and define the Hilbert spaceH = \{ f \in C( R ) \mid \operatorname{supp}(F) \subset [ - a,a] \}

where

C( R ) is the set of continuous square integrable functions, and $F(\backslash omega)\; =\; \backslash int\_^\backslash infty\; f(t)\; e^\; \backslash ,\; dt$ is the Fourier transform of f . As the inner product of this Hilbert space, we use\langle f, g\rangle _{L ^{2}} = \int ^{inf ty} _{-inf ty} f(x) ⋅ \overline{g(x)} dx.

From the Fourier inversion theorem, we have

f(x) = (1)/(2 \pi) \int ^{a} _{- a} F(\omega) e ^{ix} d\omega .

It then follows by the Cauchy–Schwarz inequality and Plancherel's theorem that, for all

x ,|f(x)| \le (1)/(2 \pi) \sqrt{ 2 a\int ^{a} _{- a} |F(\omega)| ^{2} d\omega} = (\sqrt{ 2 a)/(}{ 2 \pi}\sqrt{\int) ^{inf ty} ^{-inf ty} |F(\omega)| ^{2} d\omega} = \sqrt{ (a)/(\pi) } \|f\|_.

This inequality shows that the evaluation functional is bounded, proving that

H is indeed a RKHS.

The kernel function

K _{x} in this case is given byK _{x(y)} = (a)/(\pi) \operatorname{sinc}\left ( (a)/(\pi) (y - x) \right ) = (\sin(a(y - x)))/(\pi(y - x)) .

The Fourier transform of

K _{x(y)} defined above is given by\int ^{inf ty} _{-inf ty} K ^{- i \omega y} _{x(y)e} dy = \begin{cases} e ^{-i} & if \omega \in [ - a, a], \\ 0 & rm{otherwise}, \end{cases}

which is a consequence of the time-shifting property of the Fourier transform. Consequently, using Plancherel's theorem, we have

\langle f, K _{x\rangle} ^{L ^{2}} = \int ^{inf ty} _{-inf ty} f(y) ⋅ \overline{K _{x(y)}} dy = (1)/(2 \pi) \int ^{a} _{- a} F(\omega) ⋅ e ^{i\omega} d\omega = f(x) .

Thus we obtain the reproducing property of the kernel.

K _{x} in this case is the "bandlimited version" of the Dirac delta function, and that K _{x(y)} converges to \delta(y - x) in the weak sense as the cutoff frequency a tends to infinity.

Moore–Aronszajn theorem

We have seen how a reproducing kernel Hilbert space defines a reproducing kernel function that is both symmetric and positive definite. The Moore–Aronszajn theorem goes in the other direction; it states that every symmetric, positive definite kernel defines a unique reproducing kernel Hilbert space. The theorem first appeared in Aronszajn's Theory of Reproducing Kernels, although he attributes it to E. H. Moore.

Theorem. Suppose K is a symmetric, positive definite kernel on a set X. Then there is a unique Hilbert space of functions on X for which K is a reproducing kernel.

Proof. For all x in X, define K_x = K(x, ⋅). Let H₀ be the linear span of . Define an inner product on H₀ by

\left\langle \sum ^{n} _{j =1} b _{j} K _{y _{j}} , \sum ^{m} _{i =1} a _{i} K _{x _{i}} \right \rangle _{H _{0}} = \sum ^{m} _{i =1} \sum ^{n} _{j =1} {a _{i}} b _{j} K(y _{j,} x _{i),}

which implies

K(x,y) = \left\langle K _{x} , K _{y} \right\rangle _{H _{0}} .The symmetry of this inner product follows from the symmetry of K and the non-degeneracy follows from the fact that K is positive definite.

Let H be the completion of H₀ with respect to this inner product. Then H consists of functions of the form

f(x) = \sum ^{inf ty} _{i =1} a _{i} K _{x _{i}} (x) where \lim _{n} \sup _{p\geq0} \left\|\sum ^{n + p} _{i = n} a _{i} K _{x _{i}} \right\| _{H _{0}} = 0 .

Now we can check the reproducing property :

\langle f, K _{x} \rangle _{H} = \sum ^{inf ty} _{i =1} a _{i\left} \langle K _{x _{i}} , K _{x} \right \rangle _{H _{0}} = \sum ^{inf ty} _{i =1} a _{i} K (x _{i,} x) = f(x).

To prove uniqueness, let G be another Hilbert space of functions for which K is a reproducing kernel. For every x and y in X, implies that

\langle K _{x,} K _{y} \rangle _{H} = K(x, y) = \langle K _{x,} K _{y} \rangle _{G.}

By linearity,

\langle ⋅ , ⋅ \rangle _{H} = \langle ⋅ , ⋅ \rangle _{G} on the span of \{K _{x} : x \in X\} . Then H \subset G because G is complete and contains H₀ and hence contains its completion.

Now we need to prove that every element of G is in H. Let

f be an element of G. Since H is a closed subspace of G, we can write f = f _{H} + f _{H ^{\bot}} where f _{H} \in H and f _{H ^{\bot}} \in H ^{\bot} . Now if x \in X then, since K is a reproducing kernel of G and H:f(x) = \langle K _{x} , f \rangle _{G} = \langle K _{x,} f _{H} \rangle _{G} + \langle K _{x,} f _{H ^{\bot}} \rangle _{G} = \langle K _{x} , f _{H} \rangle _{G} = \langle K _{x} , f _{H} \rangle _{H} = f _{H(x),}

where we have used the fact that

K _{x} belongs to H so that its inner product with f _{H ^{\bot}} in G is zero.This shows that f = f _{H} in G and concludes the proof.

Integral operators and Mercer's theorem

We may characterize a symmetric positive definite kernel

K via the integral operator using Mercer's theorem and obtain an additional view of the RKHS. Let X be a compact space equipped with a strictly positive finite Borel measure \mu and K: X x X \to \R a continuous, symmetric, and positive definite function. Define the integral operator T _{K:} L _{2(X)} \to L _{2(X)} as[T _{K} f]( ⋅ ) = \int _{X} K( ⋅ ,t) f(t) d\mu(t)

where

L _{2(X)} is the space of square integrable functions with respect to \mu .

Mercer's theorem states that the spectral decomposition of the integral operator

T _{K} of K yields a series representation of K in terms of the eigenvalues and eigenfunctions of T _{K} . This then implies that K is a reproducing kernel so that the corresponding RKHS can be defined in terms of these eigenvalues and eigenfunctions. We provide the details below.

Under these assumptions

T _{K} is a compact, continuous, self-adjoint, and positive operator. The spectral theorem for self-adjoint operators implies that there is an at most countable decreasing sequence (\sigma _{i)} _{i} \geq 0 such that $\backslash lim\_\backslash sigma\_i\; =\; 0$ andT _{K\varphi} _{i(x)} = \sigma _{i\varphi} _{i(x)} , where the \{\varphi _{i\}} form an orthonormal basis of L _{2(X)} . By the positivity of T _{K,} \sigma _{i} > 0 for all i. One can also show that T _{K} maps continuously into the space of continuous functions C(X) and therefore we may choose continuous functions as the eigenvectors, that is, \varphi _{i} \in C(X) for all i. Then by Mercer's theorem K may be written in terms of the eigenvalues and continuous eigenfunctions asK(x,y) = \sum ^{inf ty} _{j =1} \sigma _{j} \varphi _{j(x)} \varphi _{j(y)}

for all

x, y \in X such that\lim _{n} \sup _{u,v} \left |K(u,v) - \sum ^{n} _{j =1} \sigma _{j} \varphi _{j(u)} \varphi _{j(v)} \right | = 0 .

This above series representation is referred to as a Mercer kernel or Mercer representation of

K .

Furthermore, it can be shown that the RKHS

H of K is given byH = \left \{ f \in L _{2(X)} \ vert \sum ^{inf ty} _{i =1} (\left\langle f,\varphi _{i} \right \rangle ^{2} _{L _{2}})/(\sigma _{i}) < inf ty \right\}

where the inner product of

H given by\left\langle f,g \right\rangle _{H} = \sum ^{inf ty} _{i =1} (\left\langle f,\varphi _{i} \right\rangle _{L _{2}} \left\langle g,\varphi _{i} \right\rangle _{L _{2}})/(\sigma _{i}) .

This representation of the RKHS has application in probability and statistics, for example to the Karhunen-Loève representation for stochastic processes and kernel PCA.

Feature maps

A feature map is a map

\varphi\colon X → F , where F is a Hilbert space which we will call the feature space. The first sections presented the connection between bounded/continuous evaluation functions, positive definite functions, and integral operators and in this section we provide another representation of the RKHS in terms of feature maps.

Every feature map defines a kernel via

Clearly

K is symmetric and positive definiteness follows from the properties of inner product in F . Conversely, every positive definite function and corresponding reproducing kernel Hilbert space has infinitely many associated feature maps such that holds.

For example, we can trivially take

F = H and \varphi(x) = K _{x} for all x \in X . Then is satisfied by the reproducing property. Another classical example of a feature map relates to the previous section regarding integral operators by taking F = \ell ^{2} and \varphi(x) = (\sqrt{\sigma _{i}} \varphi _{i(x))} _{i} .

This connection between kernels and feature maps provides us with a new way to understand positive definite functions and hence reproducing kernels as inner products in

H . Moreover, every feature map can naturally define a RKHS by means of the definition of a positive definite function.

Lastly, feature maps allow us to construct function spaces that reveal another perspective on the RKHS. Consider the linear space

H _{\varphi} = \{ f: X \to R \mid \exists w \in F, f(x) = \langle w, \varphi(x) \rangle _{F} , \forall x \in X \} .

We can define a norm on

H _{\varphi} by\|f\| _{\varphi} = inf \{\|w\| _{F} : w \in F, f(x) = \langle w, \varphi(x)\rangle _{F,} \forall x \in X \} .

It can be shown that

H _{\varphi} is a RKHS with kernel defined by K(x,y) = \langle\varphi(x), \varphi(y)\rangle _{F} . This representation implies that the elements of the RKHS are inner products of elements in the feature space and can accordingly be seen as hyperplanes. This view of the RKHS is related to the kernel trick in machine learning.^[7]

Properties

Useful properties of RKHSs:

• Let

(X _{i)} ^{p} ^{i =1} be a sequence of sets and (K _{i)} ^{p} ^{i =1} be a collection of corresponding positive definite functions on (X _{i)} ^{p.} ^{i =1} It then follows thatK((x _{1,\ldots} ,x _{p),(y} _{1,\ldots,y} _{p))} = K _{1(x} _{1,y} _{1) …} K _{p(x} _{p,y} _{p)}

is a kernel on

X = X _{1} x ... x X _{p.}

• Let

X _{0} \subset X, then the restriction of K to X _{0} x X _{0} is also a reproducing kernel.

• Consider a normalized kernel

K such that K(x, x) = 1 for all x \in X . Define a pseudo-metric on X asd _{K(x,y)} = \|K _{x} - K _{y\|} ^{2} ^{H} = 2 ( 1- K(x,y)) \forall x \in X .

By the Cauchy–Schwarz inequality,

K(x,y) ^{2} \le K(x, x)K(y, y) =1 \forall x,y \in X.

This inequality allows us to view

K as a measure of similarity between inputs. If x, y \in X are similar then K(x,y) will be closer to 1 while if x,y \in X are dissimilar then K(x,y) will be closer to 0.

• The closure of the span of

\{ K _{x} \mid x \in X \} coincides with H .^[8]

Common examples

Bilinear kernels

K(x,y) = \langle x,y\rangle The RKHS H corresponding to this kernel is the dual space, consisting of functions f(x) = \langle x,\beta\rangle satisfying \|f\| ^{2= \|\beta\|} _{H} ^{2} .

Polynomial kernels

K(x,y) = (\alpha\langle x,y \rangle + 1 ) ^{d,} \alpha \in \R, d \in \N

Radial basis function kernels

These are another common class of kernels which satisfy

K(x,y) = K(\|x - y\|) . Some examples include:

• Gaussian or squared exponential kernel:

K(x,y) = e ^{- (\|x - y\| ^{2})/(2 \sigma ^{2})} , \sigma > 0

• Laplacian kernel:

K(x,y) = e ^{- (\|x - y\|)/(\sigma)} , \sigma > 0

The squared norm of a function

f in the RKHS H with this kernel is:^{[9] [10]} \|f\| ^{2= \int} _{R} ( (1 {\sigma})/(f(x)) ^{2} + \sigma f'(x) ^{2)} d x.

Bergman kernels

We also provide examples of Bergman kernels. Let X be finite and let H consist of all complex-valued functions on X. Then an element of H can be represented as an array of complex numbers. If the usual inner product is used, then K_x is the function whose value is 1 at x and 0 everywhere else, and

K(x,y) can be thought of as an identity matrix sinceK(x,y) = \begin{cases} 1 & x = y \ 0 & x ≠ y \end{cases}

In this case, H is isomorphic to

\Complex ^{n} .

The case of

X = D (where D denotes the unit disc) is more sophisticated. Here the Bergman space H ^{2(D)} is the space of square-integrable holomorphic functions on D . It can be shown that the reproducing kernel for H ^{2(D)} isK(x,y) = (1)/(\pi) (1)/(( 1- x\overline{y) ) ^{2}.}

Lastly, the space of band limited functions in

L ^{2(\R)} with bandwidth 2 a is a RKHS with reproducing kernelK(x,y) = (\sin a (x - y))/(\pi (x - y)) .

Extension to vector-valued functions

In this section we extend the definition of the RKHS to spaces of vector-valued functions as this extension is particularly important in multi-task learning and manifold regularization. The main difference is that the reproducing kernel

\Gamma is a symmetric function that is now a positive semi-definite matrix for every x,y in X . More formally, we define a vector-valued RKHS (vvRKHS) as a Hilbert space of functions f: X \to R ^{T} such that for all c \in R ^{T} and x \in X \Gamma _{xc(y)} = \Gamma(x, y)c \in H for y \in X

and

\langle f, \Gamma _{x} c \rangle _{H} = f(x) ^{\intercal} c.

This second property parallels the reproducing property for the scalar-valued case. This definition can also be connected to integral operators, bounded evaluation functions, and feature maps as we saw for the scalar-valued RKHS. We can equivalently define the vvRKHS as a vector-valued Hilbert space with a bounded evaluation functional and show that this implies the existence of a unique reproducing kernel by the Riesz Representation theorem. Mercer's theorem can also be extended to address the vector-valued setting and we can therefore obtain a feature map view of the vvRKHS. Lastly, it can also be shown that the closure of the span of

\{ \Gamma _{xc} : x \in X, c \in R ^{T} \} coincides with H , another property similar to the scalar-valued case.

We can gain intuition for the vvRKHS by taking a component-wise perspective on these spaces. In particular, we find that every vvRKHS is isometrically isomorphic to a scalar-valued RKHS on a particular input space. Let

Λ = \{ 1 , ... , T \} . Consider the space X x Λ and the corresponding reproducing kernel

As noted above, the RKHS associated to this reproducing kernel is given by the closure of the span of

\{ \gamma _{(x,t)} : x \in X, t \in Λ \} where \gamma _{(x,t)} (y,s) = \gamma( (x,t), (y,s)) for every set of pairs

The connection to the scalar-valued RKHS can then be made by the fact that every matrix-valued kernel can be identified with a kernel of the form of via

\Gamma(x,y) _{(t,s)} = \gamma((x,t), (y,s)).

Moreover, every kernel with the form of defines a matrix-valued kernel with the above expression. Now letting the map

D: H _{\Gamma} \to H _{\gamma} be defined as(Df)(x,t) = \langle f(x), e _{t} \rangle _{R ^{T}}

where

e _{t} is the t ^{th} component of the canonical basis for R ^{T} , one can show that D is bijective and an isometry between H _{\Gamma} and H _{\gamma} .

While this view of the vvRKHS can be useful in multi-task learning, this isometry does not reduce the study of the vector-valued case to that of the scalar-valued case. In fact, this isometry procedure can make both the scalar-valued kernel and the input space too difficult to work with in practice as properties of the original kernels are often lost.^{[11] [12] [13]}

An important class of matrix-valued reproducing kernels are separable kernels which can factorized as the product of a scalar valued kernel and a

T -dimensional symmetric positive semi-definite matrix. In light of our previous discussion these kernels are of the form\gamma((x,t),(y,s)) = K(x,y) K _{T(t,s)}

for all

x,y in X and t,s in T . As the scalar-valued kernel encodes dependencies between the inputs, we can observe that the matrix-valued kernel encodes dependencies among both the inputs and the outputs.

We lastly remark that the above theory can be further extended to spaces of functions with values in function spaces but obtaining kernels for these spaces is a more difficult task.^[14]

Connection between RKHSs and the ReLU function

The ReLU function is commonly defined as

f(x) =max \{ 0 , x\} and is a mainstay in the architecture of neural networks where it is used as an activation function. One can construct a ReLU-like nonlinear function using the theory of reproducing kernel Hilbert spaces. Below, we derive this construction and show how it implies the representation power of neural networks with ReLU activations.

We will work with the Hilbert space

l{H}=L ^{1} _{2 ( 0 )[ 0 ,} inf ty) of absolutely continuous functions with f( 0 ) = 0 and square integrable (i.e. L _{2} ) derivative. It has the inner product \langle f,g \rangle _{l{H} } = \int_0^\infty f'(x)g'(x) \, dx .

To construct the reproducing kernel it suffices to consider a dense subspace, so let

f\in C ^{1[0,} inf ty) and f( 0 ) =0 . The Fundamental Theorem of Calculus then gives f(y) = \int ^{y} _{0} f'(x) dx = \int ^{inf ty} _{0} G(x,y) f'(x) dx = \langle K _{y,f} \rangle

where

G(x,y) = \begin{cases} 1 , & x < y\\ 0 , & otherwise \end{cases}

and

K _{y'(x)=} G(x,y), K _{y(0)} = 0 i.e.K(x, y) = K _{y(x)=\int} ^{x} ^{0} G(z, y) dz = \begin{cases} x, & 0 \leq x<y \\ y, & otherwise. \end{cases} =min (x, y)

This implies

K _{y=K( ⋅ ,} y) reproduces f .

Moreover the minimum function on

X x X = [ 0 , inf ty) x [ 0 , inf ty) has the following representations with the ReLu function: min (x,y) = x - \operatorname{ReLU}(x - y) = y - \operatorname{ReLU}(y - x).

Using this formulation, we can apply the representer theorem to the RKHS, letting one prove the optimality of using ReLU activations in neural network settings.

See also

• Positive definite kernel
• Mercer's theorem
• Kernel trick
• Kernel embedding of distributions
• Representer theorem

References

• Alvarez, Mauricio, Rosasco, Lorenzo and Lawrence, Neil, “Kernels for Vector-Valued Functions: a Review,” https://arxiv.org/abs/1106.6251, June 2011.
• Nachman . Aronszajn . Nachman Aronszajn . Theory of Reproducing Kernels . . 68 . 3 . 337–404 . 1950 . 1990404 . 51437 . 10.1090/S0002-9947-1950-0051437-7. free.
• Berlinet, Alain and Thomas, Christine. Reproducing kernel Hilbert spaces in Probability and Statistics, Kluwer Academic Publishers, 2004.
• Felipe . Cucker . Steve . Smale . On the Mathematical Foundations of Learning . . 39 . 1 . 1–49 . 2002 . 10.1090/S0273-0979-01-00923-5 . 1864085. free.
• De Vito, Ernest, Umanita, Veronica, and Villa, Silvia. "An extension of Mercer theorem to vector-valued measurable kernels,", June 2013.
• Durrett, Greg. 9.520 Course Notes, Massachusetts Institute of Technology, https://www.mit.edu/~9.520/scribe-notes/class03_gdurett.pdf, February 2010.
• George . Kimeldorf . Grace . Wahba . Grace Wahba . Some results on Tchebycheffian Spline Functions . Journal of Mathematical Analysis and Applications . 33 . 1 . 1971 . 82–95 . 10.1016/0022-247X(71)90184-3 . 290013. free.
• Okutmustur, Baver. “Reproducing Kernel Hilbert Spaces,” M.S. dissertation, Bilkent University, http://www.thesis.bilkent.edu.tr/0002953.pdf, August 2005.
• Paulsen, Vern. “An introduction to the theory of reproducing kernel Hilbert spaces,” https://citeseerx.ist.psu.edu/document?repid=rep1&type=pdf&doi=440218056738e05b5ab43679f932a9f33fccee87.
• Ingo . Steinwart . Clint . Scovel . Mercer's theorem on general domains: On the interaction between measures, kernels, and RKHSs . Constr. Approx. . 35 . 3. 2012 . 363–417 . 2914365 . 10.1007/s00365-012-9153-3. 253885172 .
• Rosasco, Lorenzo and Poggio, Thomas. "A Regularization Tour of Machine Learning – MIT 9.520 Lecture Notes" Manuscript, Dec. 2014.
• Wahba, Grace, Spline Models for Observational Data, SIAM, 1990.
• Zhang . Haizhang . Xu . Yuesheng . Zhang . Qinghui . 2012 . Refinement of Operator-valued Reproducing Kernels . Journal of Machine Learning Research . 13 . 91–136 .

Notes and References

1. Alpay, D., and T. M. Mills. "A family of Hilbert spaces which are not reproducing kernel Hilbert spaces." J. Anal. Appl. 1.2 (2003): 107–111.
2. Z. Pasternak-Winiarski, "On weights which admit reproducing kernel of Bergman type", International Journal of Mathematics and Mathematical Sciences, vol. 15, Issue 1, 1992.
3. T. Ł. Żynda, "On weights which admit reproducing kernel of Szegő type", Journal of Contemporary Mathematical Analysis (Armenian Academy of Sciences), 55, 2020.
4. Okutmustur
5. Paulson
6. Durrett
7. Rosasco
8. Rosasco
9. Berlinet, Alain and Thomas, Christine. Reproducing kernel Hilbert spaces in Probability and Statistics, Kluwer Academic Publishers, 2004
10. Thomas-Agnan C . Computing a family of reproducing kernels for statistical applications. Numerical Algorithms, 13, pp. 21-32 (1996)
11. De Vito
12. Zhang
13. Alvarez
14. Rosasco