• Views 582
• Citations 5
• ePub 21
• PDF 313
`Journal of Applied MathematicsVolume 2013, Article ID 407128, 8 pageshttp://dx.doi.org/10.1155/2013/407128`
Research Article

## Rate of Convergence of Hermite-Fejér Interpolation on the Unit Circle

1Departamento de Matemática Aplicada I, Facultad de Ciencias, Universidad de Vigo, 32004 Ourense, Spain
2Departamento de Matemática Aplicada I, E. Ingeniería Industrial, Universidad de Vigo, 36310 Vigo, Spain

Received 13 November 2012; Revised 2 March 2013; Accepted 4 March 2013

Academic Editor: Roberto Barrio

Copyright © 2013 E. Berriochoa et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The paper deals with the order of convergence of the Laurent polynomials of Hermite-Fejér interpolation on the unit circle with nodal system, the roots of a complex number with modulus one. The supremum norm of the error of interpolation is obtained for analytic functions as well as the corresponding asymptotic constants.

#### 1. Introduction

The paper is devoted to study the Hermite-Fejér interpolation problem on the unit circle . This topic has attracted the interest of many researchers in recent years, and it has been the subject of several studies. In  Fejér's classical result is extended to the unit circle. It is well known that it ensures uniform convergence of Hermite-Fejér interpolants to continuous functions on taking as nodal system the Chebyshev points (see ). Specifically, in  the authors consider the nodal system of the roots of a complex number with modulus one. Then it is proved that the Laurent polynomials of Hermite-Fejér interpolation for a given continuous function on the unit circle uniformly converge to .

In  second Fejér's theorem concerning the Hermite interpolation with nonvanishing derivatives is extended, to the unit circle. New conditions for the derivatives are obtained in order that the Hermite interpolants uniformly converge to continuous functions on the unit circle.

An algorithm for efficient computing of the coefficients of the Laurent polynomials of Hermite-Fejér and Hermite interpolation with equally spaced nodes on the unit circle was given in . These results were extended to the bounded interval, and the corresponding expressions can be evaluated using the techniques given in . Some results concerning the convergence were obtained in . The convergence of the Laurent polynomials of Hermite-Fejér interpolation has been studied in  for analytic functions defined on open sets containing the unit disk. The results describe the behavior outside and inside the unit disk and are extended to the case of Hermite interpolation, that is, with nonvanishing derivatives.

In the case of bounded interval, the supremum norm of the error of interpolation was studied in several papers (see ). In particular a lower bound for the order of convergence of Hermite-Fejér interpolation was obtained in  for general nodal systems. Now, in the present paper, we study the same problem for the error of interpolation on the unit circle by taking into account the results obtained in .

The organization of the paper is as follows. Section 2 is dedicated to obtain the results for the order of convergence for Laurent polynomials, in other words to the polynomial case. Section 3 contains the extension of the preceding results for analytic functions. The order of convergence and the asymptotic constants are deduced in our main result for analytic functions on an open disk containing the unit circle. As a consequence, the result is generalized to analytic functions outside an open disk with radius less than one,and it is also generalized, in Section 4, to analytic functions on an open annulus containing the unit circle. Finally, the last section is devoted to some numerical experiments to reveal the contributions of our results.

#### 2. The Polynomial Case

Let be the roots of a complex number with modulus . We recall that the Hermite interpolation problem on the unit circle with nodal system consists in obtaining a Laurent polynomial of that satisfies the following interpolation conditions: where and are sets of fixed complex numbers.

The particular case when for all is called the Hermite-Fejér interpolation problem, and the corresponding interpolation polynomial is denoted by . When , (), for a given function defined on , we denote the Hermite-Fejér interpolation polynomial by . To estimate the interpolation error between and we consider their difference that we denote by

It is well known that can be computed in terms of the fundamental polynomial of Hermite interpolation, , as follows: where is given by and it holds that Representation (3) can be seen in , and (4) can be seen in .

We recall that for a continuous function defined on , converges to uniformly on , as it can be seen in .

These results can be improved in case of polynomial functions. Indeed we can obtain nice explicit expressions for and for the polynomial case; that is, in this section we are going to use an algebraic polynomial or a Laurent polynomial in the role of .

Theorem 1. Let be a fixed positive integer number. For the following conditions hold that (i);  (ii); (iii) converges to uniformly on compact subsets of with order of convergence ; (iv); (v); (vi) converges to uniformly on compact subsets of with order of convergence .

Proof. In order to obtain (i), take into account that when we evaluate the proposed Laurent polynomial at we have that is, the interpolation conditions for the function are fulfilled. In the same way, when we evaluate the corresponding derivative at we obtain Thus the existence and uniqueness of the Hermite interpolation polynomial ensures (i). (ii) It is an immediate consequence of (i) and the definition of . (iii) Take into account that where the last expression is uniformly bounded if and is large enough.
(iv), (v), and (vi) can be proved proceeding in the same way.

Remark 2. The resulting expressions for and , given in the preceding theorem, can be rewritten as follows:

Corollary 3. The following hold.  (i) If is a Laurent polynomial with nonnegative powers of , that is, is an algebraic polynomial, then  (a).  (b) If is a compact subset of, then where is the supremum norm on . (c) If is a compact subset of with no isolated points, then  (ii) If is a Laurent polynomial with only negative powers of , then  (a).  (b) If is a compact subset of , then  (c) If is a compact subset of with no isolated points, then

Proof. (i) (a) It is a straightforward consequence of the previous remark.
(i) (b) First of all take into account that for each it holds that Therefore, if is attained at , then for each the following relation holds:
On the other hand, since for we have , then we obtain the result.
(i) (c) If is the point where is attained, then for each it holds that which implies that
Due to the continuity of , for each there exists a neighborhood of , , such that for it is . On the other hand, for large enough any arc of contains points with , and therefore for some , with , we have Then we obtain To obtain (ii) (a), (b), and (c) proceed in the same way.

Theorem 4. Let be a Laurent polynomial with positive and negative powers of , and , respectively. It holds that (i);  (ii) if is a compact subset of with no isolated points, then

Proof. It is clear that exists, it is positive, and it is attained at a point . We denote this maximum by . Besides, since can be represented as with , then
Due to the continuity of , for each there exists a neighborhood of , , such that for it is which implies, in particular, . Moreover, taking into account that for large enough any arc , with , contains points with , then we have

Remark 5. The previous result can be rewritten as follows: where means that the sequences are equivalent.

#### 3. Rate of Convergence for Analytic Functions on a Disk

In this section we extend the previous results to analytic functions. Indeed we study the supremum norm of for analytic functions on an open disk containing along a circumference of radius .

Theorem 6. Let be a nonconstant analytic function defined on an open disk , with , and let be its Hermite-Fejér interpolation polynomial corresponding to the roots of . If is the circumference with radius , , and is the supremum norm on , then there exist satisfying

Proof. We exclude the constant case because we have in this situation. So let be a nonconstant analytic function. Taking into account that the evaluation of at is , we have . Then by using Theorem 1 we can write
Let be a point belonging to with and , such that . Then there exists a positive constant , such that for large enough we have
Now we consider (2) and (3) as follows:
So we obtain that there exists a positive constant , such that for large enough
In order to obtain the lower bound we use the following inequality:
For an arbitrary and large enough we have
On the other hand, as the zeros of cannot have accumulation points on , then for there exist an arc and a positive constant , such that Next we study two different cases and .
(i) If it holds that and, as before, for large enough there exist positive constants and , such that for we have Thus for large enough and we have Then there exists satisfying  (ii) If , we consider an arc with for and a positive constant . Then
Proceeding in the same way as in the previous case we have . Thus there exists , such that
Taking into account (35) and (37) we have for large enough, and
Then it is straightforward that there exists , such that for large enough.

Remark 7. Notice that  (i) the constants and are closely related to the supremum norm ; (ii) clearly we can obtain an analogous result for nonconstant analytic function defined on with .

#### 4. Rate of Convergence for Analytic Functions on an Annulus

Next we deal with the case of analytic functions on an open annulus containing . We obtain explicit expressions for and the asymptotic behavior of its supremum norm; that is, we obtain the order of convergence and the asymptotic constant.

Throughout this section we consider a function with Laurent expansion at given by which converges on an annulus containing . Then there exist and , , such that for every .

For each we denote by and by . In the same way we denote by and by (if).

Furthermore we denote by and by .

Then we have the following decompositions of for each : By using this notation for the decompositions of we obtain the following results.

Lemma 8. In our conditions it holds that

Proof. If we have and . On the other hand, by using (3) and (4) we have for Then we can write and the result is proved.

Lemma 9. In our conditions the following holds. (i) For and , such that  (ii) For and , such that

Proof. (i) Take into account (i) (a) in Corollary 3 in order to prove the first equality. To obtain the second equality take into account that for , , and is bounded.
(ii) Proceed in the same way.

Theorem 10. In our conditions the following holds. (i) For and such that  (ii) If is a compact subset of   with no isolated points, then

Proof. (i) It is a straightforward from previous lemmas.
To prove (ii) use the same technique as in Theorem 4.

#### 5. Numerical Tests

In this section we present some numerical experiments concerning the main results in Sections 3 and 4.

Theorem 10 ensures that under appropriate assumptions, Moreover from the proofs of Lemmas 8 and 9 we can predict where this limit can be observed. In fact near , when is the point where the maximum of is attained, we would observe the convergence.

A second interesting point is that when the maximum is attained at a unique point, then for a compact set with : We have developed some numerical examples to see these phenomena about Theorem 10.

Example 11. Let be with . It is easy to see that the corresponding maximum with is attained at . Furthermore the maximum value is , and it is unique. For , we obtain the corresponding Hermite- Fejér approximants (based on the roots of ) the corresponding , and we evaluate the quotient in 5000 random points of the arc . As we have said the maximum of the quotients must converge to a value less than . As a second part of the example we evaluate the quotients in 1000 random points of the arc . This second sequence must converge to ; notice that the great number of evaluations gives an estimate of the supremum norm.
Table 1 shows the results observed for .
Figure 1 shows the graphic of for , , and . In this graphic we can see given in thw following we can see that the maximum of the error is attained near to or near and its corresponding points in are near to .

Table 2 shows the results observed for .

Next we are going to apply Theorem 6. This result and the details of its proof claim that for an analytic function , a compact arc of radius and under the corresponding assumptions, and the convergence can be increasing or decreasing; really it depends on the sign of .

We must point out that outside the unit disc the algorithms for Hermite-Fejér interpolation can be unstable, so we deal with a compact set near .

Example 12. Let be , and let , , and be the arcs , , and , respectively. It is easy to see that is attained at . So we can observe tending to and tending to a number . For , we obtain the corresponding Hermite-Fejér approximants, and the corresponding , and we obtain evaluations for the quotient in random points of the arc . As we have said the maximum of quotients must converge to a value less than . As a second part of the example we obtain evaluations for the quotients in random points of the arc , and this second sequence must converge to . Notice that the great number of evaluations gives an estimation of the supremum norm. Table 3 shows the observed results.

#### Acknowledgments

The authors want to give thanks to the referee for the valuable suggestions to improve the paper. The research was supported by Ministerio de Educación y Ciencia of Spain under Grant no. MTM2011-22713.

#### References

1. L. Daruis and P. González-Vera, “A note on Hermite-Fejér interpolation for the unit circle,” Applied Mathematics Letters, vol. 14, no. 8, pp. 997–1003, 2001.
2. P. J. Davis, Interpolation and Approximation, Dover Publications, New York, NY, USA, 1975.
3. J. C. Mason and D. C. Handscomb, Chebyshev Polynomials, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2003.
4. T. J. Rivlin, The Chebyshev Polynomials, John Wiley & Sons, New York, NY, USA, 1974, Pure and Applied Mathematics.
5. E. Berriochoa, A. Cachafeiro, and J. M. García-Amor, “An extension of Fejér's condition for Hermite interpolation,” Complex Analysis and Operator Theory, vol. 6, no. 3, pp. 651–664, 2012.
6. E. Berriochoa and A. Cachafeiro, “Algorithms for solving Hermite interpolation problems using the fast Fourier transform,” Journal of Computational and Applied Mathematics, vol. 235, no. 4, pp. 882–894, 2010.
7. C. W. Clenshaw, “A note on the summation of Chebyshev series,” Mathematical Tables and Other Aids to Computation, vol. 9, pp. 118–120, 1955.
8. E. Berriochoa, A. Cachafeiro, and E. Martínez Brey, “Some improvements to the Hermite-Fejér interpolation on the circle and bounded interval,” Computers & Mathematics with Applications, vol. 61, no. 4, pp. 1228–1240, 2011.
9. J. Szabados and P. Vértesi, Interpolation of Functions, World Scientific Publishing, Teaneck, NJ, Singapore, 1990.
10. J. Szabados, “Optimal order of convergence of Hermite-Fejér interpolation for general systems of nodes,” Acta Scientiarum Mathematicarum, vol. 57, no. 1–4, pp. 463–470, 1993.