1. Introduction
Frequently, positive kernels reproducing Hilbert spaces of continuous functions appear in some applications, and they are presented as radial basis functions (RBFs),
where
denotes the Euclidean norm in
, and
is a given smooth univariate function.
The Wendland functions [
1] yield compactly supported and differentiable functions in
that reproduce kernels of Hilbert spaces isomorphic to the Sobolev space
. Thus, when the dimension
n is even, the order of this Sobolev space is not an integer.
Robert Schaback [
2] extends the classical Wendland functions to the missing Wendland functions that reproduce kernels of Hilbert spaces isomorphic to the Sobolev spaces of integer order in even dimensions. Moreover, they have compact support. In this context, in [
3], Schaback and Wendland used compactly supported radial basis functions in order to solve some partial differential equations.
In [
4], Argáez, Hafstein, and Giesl provided a numerical code in C++ in order to calculate explicitly the Wendland function with any given parameters. Previously, in [
2], Schaback and Zhu, in [
5], used instead a code written in MAPLE. In [
6], Chen and other authors proposed a study of a surrogate model assisted by an evolutionary algorithm for high-dimensional expensive optimization problems also using this type of radial basis functions. Saberi et al. in [
7] provided the required formulas in one dimension for the Riemann Liouville fractional derivative of five kinds of RBFs, including the Powers, Gaussian, Multiquadric, Matérn, and Thin-plate splines. After, they also considered the discretization of the fractional diffusion equation with the RBF collocation method.
Radial basis function (RBF) approximations have been used for some time to interpolate data on a sphere. In this context, Fornberg and Piret, in [
8], extended the earlier works for computations in three aspects: firstly, tested with a large number of different types of radial functions; secondly, calculated in a stable way for e-values all the way down to the parameter equal to zero; thirdly, results presented at both short and long times, in order to contrast time scales appropriate for weather and for climate modeling, respectively.
In [
9], Rosenfeld and Dixon developed a pseudo spectral method for the estimation of the fractional Laplacian function using the approach by RBF interpolation.
Buhmann and Jager, in [
10], presented the connections of the monotonicity properties and the strict positive definiteness of vectorial functions. They studied a technique to construct positive definite functions from multiple monotone functions.
Chernih et al., in [
11], also demonstrated that with an appropriate rescaling of the variables, both the original and the missing Wendland functions converge uniformly to Gaussian, as the smoothness parameter tends to infinity.
To better understand the objective of this work, we believe that we should cite a brief history of the theory of the approximation problem using variational spline functions. The theory of the approach using variational splines was introduced by Attéia [
12], based on the
-splines functions, after Duchon [
13] developed the idea, using the technique of the minimization of quadratic functionals. We enriched this generic idea by minimizing various types of quadratic functionals, first in Hilbert spaces and secondly in a finite element space, such as in [
14] by Kouibia et al. We studied some interpolation and smoothing methods for constructing free-form curves and surfaces from a given Lagrangian and/or Hermite data set. These methods consist of the minimization of a certain quadratic functional in a Sobolev space.
In [
15], Kouibia et al. presented an approximation method from a given scattered data set, by minimizing a quadratic functional in a parametric finite element space. In [
16], Kouibia and collaborators considered the same problem from a given noisy data set; meanwhile, in [
17], they studied these problems in a bicubic spline functional space, and the optimal solution was obtained by a suitable optimization of some parameters that appear in the minimization functional.
In recent years, some of the authors of this article started to work on some problems of approximation using the Wendland radial basis functions. Recent publications include, for example, [
18], where González et al. proposed an approximation method for solving second-kind Volterra integral equation systems by radial basis functions. Recently, in [
19], Noorizadegan and Schaback introduced the evaluation condition number by a novel assessment of conditioning in radial basis function methods.
In this work, we deal with the smoothing problem in a finite-dimensional generalized Wendland functions space; formulating the problem of smoothing variational splines by generalized Wendland functions, we show how to compute, in practice, the solution of such a problem, and the method is justified by proving the corresponding convergence result. In order to illustrate the method, some graphical and numerical examples are presented in , and a comparison with another work is analyzed.
The remainder of this manuscript is organized as follows. In
Section 2, we present some notations and preliminaries that are necessary to formulate the problem.
Section 3 is devoted to studying the generalized Wendland compactly supported radial basis functions, while
Section 4 is dedicated to developing the problem of the smoothing variational splines by generalized Wendland functions. In the last section, we finish this article by illustrating some numerical and graphical examples and presenting a comparison with another work.
3. Generalized Wendland Compactly Supported Radial Basis Functions
Definition 1. Let there be as a continuous function, a set , and a finite set of points of Ω;
the linear space generated by the functions setis called the radial basis functions space relative to the function ψ and the centers set , where is the Euclidean inner product in . Definition 2. Consider a function and the radial basis function given bywhere are determined by the interpolating conditions Then, , if it exists, is called the interpolation RBF of u in (relative to ψ and ).
Remark 1. The interpolation RBF exists, and it is unique if and only if Robert Schaback in [
2] considered the integral operator
for all
,
.
Consider the truncated power functions for all
.
Since the operators preserve compact supports and are applicable to for all , we can define
Definition 3. We call generalized Wendland functions to given bywhich are well defined and supported in . Remark 2. Taking into account the above definition, we have In [
2], the author deduces an algorithm for constructing the generalized Wendland functions for even dimensions
in the following way (
Table 1):
for any integers
, with
and
as two associated polynomials of degree
and
, respectively.
Theorem 1. Let there be , as a centers set, and . Let be the interpolation RBF of relative to from , with .
Letbe the fill distance of in Ω, where denotes the Euclidean norm in . Then,where C is independent of f. Proof. Applying ([
20], Proposition 3.2) for
,
, and
, it is verified that
; thus, there exists a real constant
, independent of
f, such that
From Madych-Nelson ([
21], Theorem 6), it is verified that
and
where
denotes the inner product in the dual space of
.
Then, , for all , and we have that is orthogonal to .
Thus, for any
, it is verified that
, and we obtain that
Hence, we have
and taking
, we conclude that
From (
5), (
6), and Jiayin ([
22], Lemma 3.3.3), we can affirm that there exists
, independent of
f, such that
Then, there exists
, independent of
f, such that
and (
4) holds. □
4. Smoothing Variational Splines by Generalized Wendland Functions
Given a function
with
and a finite set of points
, we consider the functional
given by
and for any
, let
be the functional defined on
by
Remark 3. The first term of indicates how well v approaches f in a least discrete square sense. The second term represents a classical smoothness measure weighted by the parameter ε.
Let
be the radial basis functions space relative to the function
and the centers set
, and consider the following minimization problem: find
such that
Suppose that
A is a
-unisolvent set; that is,
and suppose that
Theorem 2. Problem (7) has a unique solution, called the smoothing variational spline in associated with A, , and ε, which is the unique solution of the following variational problem: find , such that Proof. From (
8), we have that the bilinear application
, given by
is continuous and
-elliptic. Applying the Lax–Milgram Lemma ([
23], Theorem 3.8.2) for
and the continuous linear application
given by
, there exists
, such that
and (
10) holds. Moreover,
minimizes the functional
; thus,
is the solution to Problem (
7). □
To compute the solution function
, for
, let
be the function
then,
. Applying Theorem 2, we obtain that
is the solution to the linear system
where its coefficients are given as follows:
and
Now, we prove that the smoothing variational spline converges to the function f under suitable hypotheses.
Theorem 3. Suppose the hypotheses (8) and (9) hold and thatand Proof. Let
be the interpolation RBF of
f relative to
from
; then,
, and one has
From (
4), there exists
, such that
and
Thus, from (
13)–(
15), we have that
and from (
12), we conclude that there exists
and
, such that
Moreover, from (
13)–(
15), we have that
and from (
11) and (
12), there exists
and
, such that
From (
16) and (
17), we can deduce that there exists a real constant
and
, such that
which means that the family
is bounded in
. It follows that there exists a subsequence
with
and an element
, such that
Finally, reasoning as in the points (3), (4), and (5) of the proof of ([
24], Theorem VI-3.2), we obtain the result. □
5. Numerical and Graphical Examples
To show the effectiveness of the method, we computed two relative error estimations given by
with
as five thousand distinct random points, which are some approximations of the relative error of
and
, respectively, with respect to
f in
.
From Theorems 1 and 3, these relative error estimations and tend to 0 as n tends to , under adequate conditions.
Consider the Franke function (see [
25]), given by
for any
.
Moreover, the discrete space that we use to calculate the approximated solution
is the RBFs space constructed from the generalized Wendland function
and the centers set
with
.
Table 2 shows the relative error estimation
with
(
) and
for different values of
; this specific parameter is introduced to avoid any oscillation. In this case,
. We observe that there exists an optimum value of
that could be estimated minimizing
.
Table 3 shows the relative error estimation
with
(
) and
for different values of
n. In this case,
. We observe that
decreases when
n increases, and it seems that it tends to stabilize.
Table 4 shows the relative error estimations
and
with
and
for different values of
r. We observe that
and
decrease when
r increases.
Figure 1 shows the graphs of the function
f, and
Figure 2 shows the interpolation RBF
and the smoothing variational spline
for
,
, and
, from left to right. We obtained that
and
.
6. Conclusions
While the method developed in this work is known, the use of the generalized Wendland compactly supported RBFs in this context is totally new. In fact, the question that one can ask is why use these functions? The answer is that the time cost of programming these functions is quite reduced, if we compare it, for example, to the variational splines mentioned in the references [
14,
15,
16,
17]. Moreover, the order of the degree of approximation, represented with the calculation of the estimate of the interpolation error
and the smoothing error
, with 500–1000 approximation points are of an order between
and
in most cases, as shown in
Table 2,
Table 3 and
Table 4, while in Table 2 subsection 5.2.2 of [
14], the degree of approximation with 900 points of approximation is
. All this shows the improvement and the effectiveness of the approximation method studied in this manuscript.
As a subject for another manuscript in the future or as an open topic, we think it is possible to extend the study to higher dimensions.