1. Introduction
The polynomial interpolation is a classical subject that has been widely studied under different approaches like Lagrange, Hermite, Birkhoff, Pál-type interpolation and some others. Although it is obvious that the subject is important by itself, its numerous numerical applications like numerical integration or numerical derivation are not less important and indeed the polynomial interpolation continues being a subject of current research.
Lagrange interpolation is a very good tool although it is known that for this interpolation scheme and for good nodal systems such as the classical Chebyshev nodes there exists a continuous function on
for which the Lagrange interpolation polynomial diverges (see [
1]). A similar problem has been posed for arbitrary arrays and it was proved in [
2] that for each nodal array in
, there exists a continuous function such that the Lagrange polynomial interpolation diverges almost everywhere. In any case, recalling the words written by Trefethen in his paper [
3] we can say that there is nothing wrong with Lagrange polynomial interpolation. "Yet the truth is, polynomial interpolants in Chebyshev points always converge if f is a little bit smooth". As a consequence, to obtain better results one needs to assume better properties on the function to be interpolated like bounded variation or a condition on its modulus of continuity. Thus one of the most important questions in relation with the interpolation of functions is the choice of the interpolation arrays or nodal systems for which one can expect to obtain pointwise or uniform convergence to the function to be interpolated and another important issue is to determine the class of functions for which some type of convergence can be guaranteed. The nodal systems strongly normal or normal, introduced by Fejér, play an important role in the interpolation theory, although from a practical point of view, the difficulty of testing the definition makes the applications of these systems quite limited. For these systems Grünwald studies in [
4] the convergence of the polynomials of Lagrange interpolation for functions satisfying a Lipschitz condition.
Most of the research obtain results on the convergence from the distribution properties of the nodal points. Indeed it was Fejér, who was the first to invert the problem, trying to deduce separation properties of the nodal systems from the interpolation results. The importance of this idea is avaled by the fact that the required interpolation properties are easily verified.
In [
5] it is proved that strongly normal distributions on
give quasi uniformly nodal systems on the unit circle, that is the length of the arcs between two consecutive nodes has the order of
. Although the situation more widely studied corresponds to the bounded interval, there are important results in some other situations in the complex plane like the unit circumference, (see [
6]). By taking into account that continuous functions on the unit circle can be approximated by Laurent polynomials, the interpolation polynomials on the unit circumference are constructed in this Laurent space. In this field of research, [
7] deserves to be highlighted. There, the nodal points are the roots of complex numbers with modulus 1 and in this case it is obtained a result about convergence of the interpolants for continuous functions with a suitable modulus of continuity. Moreover, in the same paper the problem of extending the results to general nodal systems. Indeed, since the roots of complex numbers with modulus 1 can be interpreted like the zeros of the para-orthogonal polynomials with respect to the Lebesgue measure, now the natural extension is to consider the zeros of the para-orthogonal polynomials with respect to other measures.
In [
8] we have generalized the results given in [
7] for these new nodal systems. First we work with nodal systems characterized by fulfilling some properties of boundedness, which are suggested by those fulfilled by the roots of unimodular complex numbers, obtaining a result of convergence for continuous functions with a suitable modulus of continuity. Next, by taking into account that the zeros of the para-orthogonal polynomials with respect to measures in the Szego class (see [
9]) with analytic extension up to
satisfy the properties that we need, we obtain a similar result about convergence for these type of nodal systems.
In [
10] we have studied the Lagrange interpolation process for piecewise continuous functions with suitable properties and by using as nodal points the zeros of the para-orthogonal polynomials with respect to analytic weights, which constitutes a novel approach to the Lagrange interpolation theory.
Another extension to more general nodal systems is given in [
11] where it has opened a new trend to interpolation at perturbed roots of unity and the functions to be interpolated belong to the disc algebra.
Now, in the present paper we assume a distribution for the nodes that can be obtained through a perturbation of the uniform distribution and, in particular of the roots of the unity, and which is more general than that given in [
11]. Thus in the present paper we start from a different point of view because we base it on properties satisfied by the nodal systems and we do not need to consider orthogonality nor para-orthogonality with respect to any measure on the unit circle. The interpolation arrays are described by a separation property and the main goal is to obtain the properties that play a role in the Lagrange process, as well as to present some relevant examples.
The organization of the paper is the following. In
Section 2 we introduce the nodal systems that we use throughout all the paper and we prove the main properties that they satisfy in several propositions.
Section 3 is devoted to the Lagrange interpolation problem. First we obtain the Lebesgue constant of the process and then we study the convergence of the Lagrange interpolation polynomials related to continuous functions with appropriate modulus of continuity. Secondly we analyze the rate of convergence when we deal with smooth functions, (see [
12]) and we also deduce analogous behavior for interpolation on the bounded interval. The last section is devoted to give some numerical examples.
2. Some General Nodal Systems on the Unit Circle
The aim of this paper is to study interpolation problems on the unit circle by using nodal systems satisfying some suitable properties.
We denote the nodal polynomials by and their zeros by , that is, we assume that , where for , and for For simplicity, in the sequel we omit the subscript n and we write instead of for We denote the length of the shortest arc between any two points of the unit circle, and , by and we use the Landau’s notation for complex sequences, denoting by if is bounded.
Throughout all the paper we assume that the zeros
of the nodal polynomials
satisfy the following separation property: there exists a positive constant
A such that for
the length of the shortest arc between two consecutive nodes
and
, satisfies:
where
, that is,
We use the same
to denote different sequences. Unless we mention otherwise, the bounds we obtain from (
1) will be uniform.
We also consider other nodal polynomials,
, well connected with
. If we denote
, with
, then
, where
Hence it is clear that the separation property satisfied by the zeros
of
is
In this section we obtain in several propositions the main properties of the nodal polynomials . First we recall the following well known relations between arcs and strings that we are going to use throughout the whole paper and which is based on the convex character of the arcsin function:
If
and
belong to
then
Proposition 1. If and , with , are the nodal points satisfying the separation properties (1) and (2) and we assume they are numbered in the clockwise sense, then Proof. (i) By applying (
1) we have
and for
it holds
Then if we sum, it is straightforward (i).
(ii) Proceeding in the same way we get and for it holds Hence, by computing the sums we have (ii).
(iii) We know that
and we distinguish two possibilities depending on the position of
related to
. If
that is,
and we use that
, which is a consequence of (
2), and we also take into account (i) we get
The second case corresponds to and it can be deduced in the same way.
(iv) We proceed like in (iii) distinguishing the following cases or □
Notice that we can write (iii) and (iv) as follows and
Proposition 2. Let us assume that the zeros of the nodal polynomials satisfy the separation property (1). Then it holds Proof. Since
for every
j, let us take
, such that
z is not a nodal point. In order to obtain the result we renumber the nodes in the clockwise sense in such a way that
is the nodal point nearest to
z. We distinguish two cases depending on whether the node closest to
z is turning in the clockwise sense or in the counter clockwise sense from z. If we assume that the situation is given in
Figure 1, that is,
is turning in the counter clockwise sense from
z, then we have
Now we consider the polynomial
, with
and satisfying (
2), introduced at the beginning of the section. Using property (
1) we have
and then
.
Now, for every
j it holds
and therefore, by using Proposition 1,
Besides, from property (
1) we have
and for
. Hence
Proceeding in the same way and taking into account that for
it holds that
and
we obtain
Therefore we have that
and since
, then we get
Notice that if the node closest to z, , is in the clockwise sense from z, we can proceed in a similar way. Indeed and since then .
The second statement, related to the first and second derivatives of the nodal polynomial, is a consequence of Bernstein’s theorem, (see [
13]). □
Proposition 3. Let us assume that the zeros of the nodal polynomials satisfy the separation property (1). Then there exists a positive constant such that for n large enough and for every , it holds that Proof. For simplicity we take and we try to bound from below .
Thus we consider the polynomial
satisfying (
2), that is,
with
.
Since then and so our aim is to prove that .
Now, by taking into account that
we study the quotients
and
On the one hand,
and by applying (
3) and Proposition 1, we have for
,
and therefore
.
On the other hand,
and since for
we have
then
. Hence
that is,
Thus, given
if
, then for
n large enough it holds that
Notice that for another value of j there is no any significant difference. Indeed to obtain that , we take the auxiliary polynomial with , we renumber the nodes as in the previous proof and we proceed in a similar way. □
Proposition 4. Let us assume that the zeros of the nodal polynomials satisfy the separation property (1). Then there exists a positive constant such that Proof. Following the same steps of the proof of Lemma 2 in [
14] we have
and therefore by using (
4) and its consequences in Proposition 2 we get
□
Remark 1. The nodal systems considered in [15] satisfy condition (1). Indeed they are the para-orthogonal polynomials related to measures in the Baxter class, (see [16]). In that work it is also assumed the additional condition that the sequence is uniformly bounded on , where is the sequence of monic orthogonal polynomials related to the measure and is the sequence of the reciprocal polynomials, (see [9]). In that situation studied in [15], properties (4)–6) also hold. Now, in the present paper we start from a different point of view because we base it on properties satisfied by the nodal systems and we do not need to consider orthogonality nor para-orthogonality with respect to any measure. 3. Lagrange Interpolation. Convergence in Case of Smooth Continuous Functions
To compute the interpolation polynomials, first we recall some well known definitions related to interpolation problems on the unit circle. We work in the space of Laurent polynomials and, in particular, in the subspaces , with p and q integers .
Let us continue denoting by
the zeros of the the nodal polynomial
. If
are arbitrary complex numbers, the Laurent polynomial of Lagrange interpolation
characterized by satisfying
has the following expression
If F is a function and we denote the corresponding Laurent polynomial If n is odd, since , then the interpolation polynomial and when n is even, taking into account that , then the Laurent polynomial of Lagrange interpolation
Without loss of generality, to fix ideas and to simplify the notation we assume that the number of nodes is even,
, in which case the interpolation polynomial
belongs to the space
and it can be written in terms of the fundamental polynomials as follows:
In order to compute the interpolation polynomials for applications and examples it is more convenient to use the barycentric expression, which is given by
with
, (see [
17]).
This last expression has some advantages due to its numerical stability in the sense established in [
18]. In this article author claims literally:
The Lagrange representation of the interpolating polynomial can be rewritten in two more computationally attractive forms: a modified Lagrange form and a barycentric form. We give an error analysis of the evaluation of the interpolating polynomial using these two forms. The modified Lagrange formula is shown to be backward stable. The barycentric formula has a less favourable error analysis, but is forward stable for any set of interpolating points with a small Lebesgue constant. Therefore the barycentric formula can be significantly less accurate than the modified Lagrange formula only for a poor choice of interpolating points.
So with a good Lebesgue constant (see next Theorem 1) we have good accuracy, at least as good as the intensively used Lagrange interpolation on the Chebyshev nodal systems.
Following [
10] we can obtain the Lebesgue constant, (see [
19]), and the convergence of this interpolatory process. Notice that this is a novelty result for our general nodal systems satisfying property (
1), although the techniques that we use are the same as in [
10].
Theorem 1. There exists a positive constant such that for every function F bounded on it holds that for every , where denotes the supremum norm on .
Proof. Let
z be an arbitrary point of
and assume that
z is not a nodal point. Then, if we continue assuming the even case, from (
7) we get
and by our hypothesis about
F and by Proposition 3 we have
If we assume that the nodal points closest to
z are
and
then by applying (
1) we obtain that for
it holds
Thus, by using (
3) we obtain
for some positive constant
E.
Proceeding in the same way we get and therefore
We also obtain
for some positive constant
K.
Indeed
Hence
with
the harmonic number equal to
and
P a positive constant. □
Remark 2. When the values of are affected by any type of error, which we can suppose is bounded by some , then the previous result ensures us, taking into account the linearity of the interpolation process, that the final result is affected by an error bounded by , that is, it is at least so good as the intensively used Lagrange interpolation on the Chebyshev nodal systems.
For applying the interpolation it could be very useful the following results concerning the convergence and the rate of convergence for smooth continuous functions (see [
10,
12]).
Theorem 2. (i) Let be a function defined on . If F is continuous with modulus of continuity , then converges uniformly to F on .
(ii) Let be a function defined on . If with for , with then uniformly converges to F on and the rate of convergence is .
(iii) If is an analytic function in an open annulus containing , then uniformly converges to F on . Besides, the rate of convergence is geometric.
Proof. The results are consequence of the preceding Theorem 1 and they are also based on the properties satisfied by our nodal systems. Thus one can obtain these results following the same steps as in the proof of Theorems 3 and 4 in [
10], where one can see the details. □
The Case of the Bounded Interval
We recall that the Lagrange interpolation polynomial
related to a nodal system
in
and satisfying the conditions
is given by
where
When
for a function
f, we write
In this subsection we consider the nodal polynomial with and numbered in the following way: .
We also assume that the nodes satisfy the following separation property:
There exists a positive constant A such that for
- (i)
If and then with
- (ii)
If and then with and with
- (iii)
If and then with and with
- (iv)
If and then with and with and with
Under the above assumptions we obtain the following results about the convergence and the rate of convergence for the interpolation polynomials with these nodal systems.
Theorem 3. If f is a continuous function on and , then the interpolation polynomial fulfilling the interpolation conditions converges uniformly to f on .
Proof. Let us define a continuous function
F on
by
with
. Then it is clear that its modulus of continuity satisfies
To fix ideas we assume that and , that is, case (iv). By applying Szego transformation where , that is, with . Hence we consider the nodal polynomial .
Now our nodal system is constituted by
and the arguments are
. If we renumber the arguments in such a way that
, then it holds that
with
for
and
.
Let
be the Lagrange interpolation polynomial satisfying the conditions
Since F satisfies the hypothesis of Theorem 2 (i) then converges uniformly to F on .
Analizying the expression of
and by taking into account that
as well as
have real coefficients we get that
has real coefficients. Since it is clear that
satisfies the same interpolation conditions, now we consider the algebraic polynomial
, which has real coefficients and satisfies the interpolation conditions (
9). Since the polynomial satisfying (
9) is uniquely determined, then
and it converges uniformly to
f on
.
When (case (iii)), or (case (ii)), or and (case (i)), one proceeds in a similar way and the auxiliary nodal polynomials are given by or or respectively. Notice that in the four cases the nodal polynomials have real coefficients. □
Notice that some well known results related to Lagrange interpolation with Chebyshev and Chebyshev extended nodes are particular cases of the above theorem (see [
9,
20]).
Notice that from the proof of the above theorem and by applying Theorem 1 we can obtain an analogous bound, the Lebesgue constant, for this interpolatory process, that is, there exists a positive constant
M such that
In order to obtain information concerning the rate of convergence, first we recall the following result about the expansion of an analytic function in a Jacobi series (see [
9,
21]). For a more actual version see [
12].
Theorem 4. Let be analytic on the closed segment . The expansion of f in a Jacobi series, , is convergent in the interior of the greatest ellipse with foci at , in which f is regular. The expansion is divergent in the exterior of this ellipse and the sum R of the semi-axes of the ellipse of convergence is .
Thus, in our conditions we have the following results, which are in concordance with Theorem 2.
Theorem 5. - (i)
If f is a function defined on by , where is the Chebyshev polynomial of degree k and with , with , and , then the Lagrange interpolation polynomial converges to f with rate of convergence
- (ii)
If f is analytic on the closed segment , then the Lagrange interpolation polynomial converges to f with rate of convergence geometric.
Proof. (i) If we decompose
f like
we have that
and
where the last inequality follows from the integral criterion. Thus, by applying the Lebesgue constant obtained before in (
10), we get
Hence
for some
.
(ii) Since f is analytic on , then it can be analytically continued to a neighborhood of in the complex plane. Hence the expansion in Chebyshev series converges to f in the interior of the greatest ellipse with foci at , known as Bernstein ellipse and the sum R of the semi-axes of the ellipse of convergence is . Then it holds that , for some and .
Proceeding in the same way as before we have
and
Thus, by applying the Lebesgue constant we get
and therefore
for some
and
□
4. Numerical Examples
We have carried out different numerical experiments to visualize the main contributions of this article. The first examples correspond to the three cases of Theorem 2 and in all of them we work in the following way:
We construct the nodal systems in a quite random way. We consider four arcs or sections in the unit circumference . The first one begins in and its nodes are constructed in counter clockwise sense separated by an angular length , where the are random errors determined by using the uniform distribution in . The fourth section begins in and its nodes are constructed in clockwise sense with arcs of angular length , where the are random errors determined by using the uniform distribution in . The second section begins after the first one and its nodes are constructed in counter clockwise sense with arcs of angular length , where the are random errors determined by using the uniform distribution in . Finally, in the third section the arcs between the nodes are all equal.
Obviously we use different values of n and we must remark that we obtain always the same results, really we must say similar results because due to our random choice we never have the same nodal system.
We consider a test function , that we detail in each example, and we always plot in black.
We consider the Lagrange interpolation polynomial , which interpolates the test function at the nodal system. We always plot in red.
These examples are devoted to visualize the items (i), (ii) and (iii) respectively of Theorem 2.
Example 1. In this example we work with for , which satisfies the hypotheses of Theorem 2 (i). We take , and we use (8) to obtain . We represent the function which takes real values and, as we have said, the real part of the interpolation polynomial. Notice that due to its variability, F is a quite difficult function to interpolate. Indeed, it is easy to check that is not differentiable at .
We present in Figure 2 two graphics. On the left we have a general panoramic of the interpolation along and we have added the interpolation points in green. We must point up that the interpolatory process is successful where the function has no variability. However, we have an unsuccessful situation where the function has great variability. In the graphic on the right we have a detailed situation between 1.2 and 2, that is near , which can help us to understand the problem. According to the theory presented, we must increase the number of nodes to obtain better results in this region.
Example 2. Now we consider the function defined on by , which satisifies the hypotheses of Theorem 2. In the next Figure 3 we plot on the left and for and . Notice that they are indistinguishable. On the right we plot the errors given by with . We point out that the errors are less or equal than . In the next example we also construct an alternative interpolation polynomial based on the equispaced nodal system on , but using the values of the function on our nodal system. We do this because a natural criticism to our method could be that with errors as we can be so close to the equispaced nodal system to accept this approximation. We denote by this alternative interpolation polynomial.
Example 3. In this example we take , , and we use (8) to obtain the interpolation polynomials and . Taking into account that F is analytic we know that F and must be close. In Figure 4 we plot in black, in red and in green for with . On the left hand side we have a global vision with and we can observe that and are indistinguishable; in fact for this example the maximum error was . Although has a similar shape, see that it can drive us to catastrophic errors. On the right hand side we present a detail of the previous one, which give us an idea of the error. Notice that in general we cannot know the details of the nodal distribution.
We have done the same with the imaginary part and we obtain the same results.
Example 4. Finally we choose defined on as the characteristic function of the superior arc S of , we take , and we use expression (8) to obtain . We know the behavior when the nodal system is related to para-orthogonal polynomials with respect to an analytic positive measure (see [10]), but we do not have a theory about the behavior of in our situation. We plot the results in Figure 5. Notice that the basic ideas of the Gibbs–Wilbraham phenomenon are present in this graphic, that is, the convergence of the interpolator to the function in regions which are far enough from the discontinuities and a heavy oscillation near the discontinuities. A representation of the oscillation and its amplitude, maybe, deserves a detailed study. 5. Discussion
Usually, the nodal systems used for interpolation problems are strongly connected with measures on the bounded interval and on the unit circle and their corresponding orthogonal or paraorthogonal polynomials. We must point out that these choices are very suitable to construct the whole theory but in some numerical applicatons it is possible that the nodal systems do not satisfy this requisite. So, the starting point of the paper is a distribution for the nodes that can be obtained through a perturbation of the uniform distribution and, in particular of the roots of the unity, and which is more general than that related to measures and orthogonality.
The results of this article contribute to elaborate a theory over these type of nodal systems, as well as to the Lagrange interpolation theory based on these interpolatory arrays. Moreover, a theory about the rate of convergence for some types of smooth functions is given. Finally, we translate the results to perturbed Chevyshev nodal systems and to Lagrange interpolation on the bounded interval.
We think that this research could be of interest for some mechanical models that generate these types of nodal systems. As an example we consider the next problem.
Let us suppose that we are studying a equatorial characteristic
of a planet which depends on the angle
and we have a theory which establishes that
is an analytic function. We observe the phenomenon using an observatory in the boundary of a spatial station in an elliptic orbit of period
T which rotates over itself with period
(with
and
n large enough). Moreover, we take our observations when the center of the planet, the observatory and the center of the station are aligned. We can translate the problem thinking that the planet is our Sun, the spatial station is our Earth and the observatory is our city. So the time between observations is the equivalent of a solar day. It is well known that the duration of a solar day is not constant (see
https://en.wikipedia.org/wiki/Equation_of_time for a brief introduction about the so called
Equation of time), in our case have a little oscillation on
and our observations are taken on a nodal system which satisfies (
1). Notice that in this case we do not have a equispaced distribution nor the support of the theory of Orthogonal Polynomials. Therefore, before this paper we did not know how to use our data to reconstruct
and after this paper we can be confident about the use of Lagrange interpolation.
Some future research directions could be the study of other types of interpolation on the unit circle and on the bounded interval by using these general interpolatory arrays; as well as to study the correspondig Gibbs–Wilbraham phenomena.
6. Materials and Methods
The experiments given in the section Numerical examples were obtained by using personal codes elaborated with Mathematica
® 12 (Wolfram Research Europe Ltd, Long Hanborough Oxfordshire, United Kingdom). These programs to obtain the nodal points and to compute the interpolation polynomials as well as the plots of the test functions and their interpolators are available at the public repository
https://www.dropbox.com/sh/0cx9chq3jfzov2w/AAA_SvL2i7HlC7ChMGpuG-Ata?dl=0 There one can find the program related to Example 2. To obtain the other examples some minor changes must be done.