Return to Article Details Localization results for the Lagrange max-product interpolation operator based on equidistant knots

Localization Results for the Lagrange Max-Product Interpolation Operator Based on Equidistant Knots

Lucian Coroianu§ Sorin G. Gal

September 15, 2013.

§Department of Mathematics and Computer Science, University of Oradea, Universitatii str. no. 1, 410087 Oradea, Romania, e-mail :

lcoroianu@uoradea.ro

Department of Mathematics and Computer Science, University of Oradea, Universitatii str. no. 1, 410087 Oradea, Romania, e-mail : galso@uoradea.ro

The work of both authors has been supported by a grant of the Romanian National Authority for Scientific Research, CNCS–UEFISCDI, project number PN-II-ID-PCE-2011-3-0861.

In the class of strictly positive functions strong localization results are obtained in approximation by the Lagrange max-product interpolation operators based on equidistant nodes. The results allow to approximate locally bounded strictly positive functions with very good accuracy. Then, it is observed that the results can be extended to bounded functions of variable sign.

MSC. 41A05, 41A40, 41A27, 41A36.

Keywords. Lagrange max-product interpolation operator, localization result, local direct approximation result, Lipschitz function on subintervals.

1 Introduction

Based on the Open Problem 5.5.4, pp. 324-326 in [ 14 ] , in a series of recent papers we have introduced and studied the so-called max-product operators attached to the Bernstein polynomials and to other linear Bernstein-type operators, like those of Favard-Szász-Mirakjan operators (truncated and nontruncated case), see [ 1 ] , [ 3 ] , Meyer-König and Zeller operators, see [ 4 ] , Baskakov operators, see [ 6 ] , [ 7 ] and Bleimann-Butzer-Hahn operators, see [ 5 ] .

For example, in the recent paper [ 2 ] , starting from the linear Bernstein operators Bn(f)(x)=k=0nbn,k(x)f(k/n), where bn,k(x)=(nk)xk(1x)nk, written in the equivalent form

Bn(f)(x)=k=0nbn,k(x)f(k/n)k=0nbn,k(x)

and then replacing the sum operator Σ by the maximum operator , one obtains the nonlinear Bernstein operator of max-product kind

Bn(M)(f)(x)=k=0nbn,k(x)f(kn)k=0nbn,k(x),

where the notation k=0nbn,k(x) means max{bn,k(x);k{0,...,n}} and similarly for the numerator.

For this max-product operator, nice approximation and shape preserving properties were found in the class of positive valued functions, in e.g. [ 2 ] , [ 12 ] .

In other two recent papers [ 9 ] and [ 10 ] , this idea is applied to the Lagrange interpolation based on the Chebyshev nodes of second kind plus the endpoints, and to the Hermite-Fejér interpolation based on the Chebyshev nodes of first kind respectively, obtaining max-product interpolation operators which, in general, (for example, in the class of positive Lipschitz functions) approximates essentially better than the corresponding Lagrange and Hermite-Fejér interpolation polynomials.

Let I=[a,b], a<b and f:[a,b]R. The max-product Lagrange interpolation operator on equidistant knots attached to the function f is given by (see [ 11 ] )

Ln(M)(f)(x)=k=0nln,k(x)f(xn,k)k=0nln,k(x),xI,nN,

where xn,k=a+(ba)k/n for all nN and k{0,1,...,n} and

ln,k(x)=(1)nk(i=0n(xxn,i))1xxn,k

for all xI, nN and k{0,1,...,n}. Note that Ln(M)(f) is a well defined function. Indeed, using the fundamental Lagrange polynomials,

pn,k(x)=(xxn,0)(xxn,1)...(xxn,k1)(xxn,k+1)...(xxn,n)(xn,kxn,0)(xn,kxn,1)...(xn,kxn,k1)(xn,kxn,k+1)...(xn,kxn,n),

we observe that we can rewrite ln,k(x), xI, in the form

ln,k(x)=cn,kpn,k(x),

where

cn,k=(xn,kxn,0)(xn,kxn,1)...(xn,kxn,k1)(xn,k+1xn,k)...(xn,nxn,k).

Then, since for any xI we have i=0npn,i(x)=1 it follows the existence of i(x){0,1,...,n} such that pn,i(x)(x)>0 and noting that cn,i(x)>0 it easily results that ln,i(x)(x)>0 and this implies that k=0nln,k(x)>0 for all xI, which means that indeed Ln(M)(f) is a well defined function on [a,b].

The max-product operator Ln(M)(f)(x) is continuous on [a,b] and has the interpolation properties Ln(M)(f)(xn,j)=f(xn,j) for all j{0,1,...,}.

Also, according to Corollary 3.2, (i), in [ 11 ] , for positive valued functions, i.e. for f:[a,b]R+, it satisfies the Jackson-type estimate

|Ln(M)(f)(x)f(x)|2ω1(f;ban)[a,b], for all x[a,b],nN,

where ω1(f;ban)[a,b] denotes the modulus of continuity of f on [a,b]. This estimate for the Lagrange max-product operator essentially improves for positive valued functions the order of approximation by the classical Lagrange interpolation polynomials on equidistant nodes, when as it is well-know, we can also have a very pronounced divergence phenomenon in [a,b] (see e.g. Chapter 4 in the book [ 17 ] , see also [ 16 ] , [ 8 ] ).

It is worth noting that saturation and local inverse results for Ln(M)(f)(x) were obtained in [ 13 ] .

The plan of the paper goes as follows. In Section 2 an interesting strong localization result for the Lagrange max-product operator Ln(M) is obtained. At the end of the section and as consequences of this localization result, a local direct result and an interesting local shape preserving property are proved. Section 3 contains comparisons with some linear interpolation operators of rational type.

It is worth noting in Section 2 the strong localization result expressed by Theorem 2.1, that shows that if the continuous strictly positive functions f and g coincide on a subinterval [α,β][0,1], then for sufficiently large values of n, Ln(M)(f) and Ln(M)(g) coincide on subintervals sufficiently close to [α,β]. Clearly, Corollary 2.4 shows that Ln(M)(f) is very suitable to approximate continuous functions which are constant on some subintervals. Namely, if f is a continuous strictly positive function which is constant on some subintervals [αi,βi], i=1,...,p, of [a,b], then for sufficiently large n, Ln(M)(f) takes the same constant values on subintervals sufficiently close to each [αi,βi], i=1,...,p.

2 Localization Results

Let  ln,k denote the fundamental Lagrange polynomials attached to the knots xn,k=k/n, k{0,1,...,n}, nN.

The main result of this section is the following localization result.

Theorem 2.1

Let f,g:[0,1][0,) be both bounded on [0,1] with strictly positive lower bounds and suppose that there exist a,b[0,1], 0<a<b<1 such that f(x)=g(x) for all x[a,b]. Then for all c,d[a,b] satisfying a<c<d<b there exists n~N which depends only on f,g,a,b,c,d such that Ln(M)(f)(x)=Ln(M)(g)(x) for all x[c,d] and nN, nn~.

Proof â–¼
Let us choose arbitrary x[c,d] and for each nN let jx{0,1,...,n} ( jx depends on n too, but there is no need at all to complicate on the notations) be such that x[jx/n,(jx+1)/n]. Then we know that
Ln(M)(f)(x)=k=0nln,k(x)f(kn)k=0nln,k(x)=kJn(x)ln,k(x)f(kn)kJn(x)ln,k(x)

where Jn(x)={k{0,1,...,n}:ln,k(x)>0} and ln,k, k{0,1,...,n} are the Lagrange fundamental polynomials attached to the knots xn,k=k/n, k{0,1,...,n}. Since x[c,d][jx/n,(jx+1)/n] and since a<c<d<b it is immediate that for nn0 where n0 is chosen such that 1/n0<min{ca,db}, we obtain a<jx/n<b which gives na<jx<nb for all nn0 (indeed, if we would suppose that there exists n>n0 which does not satisfy the previous double inequalities, then we would easily get a contradiction).

It is important to notice here that n0 does not depend on x. From the inequalities na<jx<b it follows that if nn0 then for any x[c,d] there exists αx[a,b] such that jx=nαx.

In what follows, it will serve to our purpose to use the sequence (an)n1, an=n. For this sequence there exists n1N such that naan>0 for all nn1.

Our intention is to prove as an intermediate result, that there exists an absolute constant N0N which does not depend of x[c,d] such that for any nN0 and x[c,d] we have k=0nln,k(x)f(kn)=kIn,xln,k(x)f(kn) where In,x={kJn(x):jxankjx+an}. In order to obtain this conclusion, for nmax{n0,n1} let us choose kJn(x)In,x. We have two cases: i) k+an<jx and ii) jx+an<k.

Case i) Firstly, note that jxJn(x), because sign(ln,j(x))=(1)nj(1)nj=1. Noting that k/n<(jxan)/n and nxjx, we get   

ln,jx(x)f(jxn)ln,k(x)f(kn)=|xk/n||xjx/n|f(jxn)f(kn)=xk/nxjx/nf(jxn)f(kn)x(jxan)/n1/nf(jxn)f(kn)=(nxjx+an)f(jxn)f(kn)nf(jxn)f(kn).

Then, denoting the infimum and the supremum of  f on [a,b] with mf and Mf respectively (according to the hypotheses these values are strictly positive), we get that

ln,jx(x)f(jxn)ln,k(x)f(kn)nmfMf

Since limnnmfMf=, it follows that there exists n2N, n2max{n0,n1} such that ln,jx(x)f(jxn)ln,k(x)f(kn)>1 for all x[c,d], nn2 and k{0,1,...,n}, k<jxan (as kIn,x). In addition, it is important to notice that n2 does not depend on x[c,d], but of course it depends on f.

Case ii) The proof is identical with the proof of the above Case i) and therefore we conclude that there exists an absolute constant n3N which depends only on a,b,c,d,f such that

ln,jx(x)f(jx/n)ln,k(x)f(k/n)>1

for all x[c,d], nn3 and k{0,1,...,n}, k>jx+an.

Analyzing the results obtained in cases i)-ii), it results that for all x[c,d], nN0, N0=max{n2,n3} and k{0,1,...,n}, with k<jxan or k>jx+an, we have

ln,jx(x)f(jx/n)ln,k(x)f(k/n)>1.

Since from the Case i) we know that jxJn(x) and since this easily implies that actually jxIn,x, we obtain our preliminary result, that is

k=0nln,k(x)f(kn)=kIn,xln,k(x)f(kn),

where In,x={kJn(x):jxankjx+an}.

Next, let us choose arbitrary x[c,d] and nN so that nN0. If there exists k In,x such that k/n[c,d] then we distinguish two cases. Either k/n<c or k/n>d. In the first case we observe that

0<cknxkn(jx+1)n+1kn(jx+1)nknan+1n.

Since limnan+1n=0, it results that for sufficiently large n we necessarily have an+1n<ca which clearly implies that k/n[a,c]. In the same manner, when k/n>d, for sufficiently large n we necessarily have k/n[d,b].

Summarizing, there exists N~1N independent of any x[c,d], such that

k=0nln,k(x)f(kn)=kIn,xln,k(x)f(kn), nN~1

and in addition for any x[c,d], nN~1 and kIn,x, we have k/n[a,b]. Also, it is easy to check that N~1 depends only on a,b,c,d,f. We thus obtain that

Ln(M)(f)(x)=kIn,xln,k(x)f(kn)k=0nln,k(x),nN~1,x[c,d]

and in addition for any x[c,d], nN~1 and kIn,x, we have k/n[a,b].

Reasoning for the function g exactly as in the case of the function f, it follows that there exists N~2N which depends only on a,b,c,d,g such that

Ln(M)(g)(x)=kIn,xln,k(x)g(kn)k=0nln,k(x),nN~2,x[c,d]

and in addition for any x[c,d], nN~2 and kIn,x, we have k/n[a,b]. Taking n~=max{N~1,N~2}, we easily obtain the desired conclusion.

We can easily extend the above result to arbitrary intervals, as follows.

Theorem 2.2

Let f,g:[a,b][0,) (a<b) be both bounded on [a,b] with strictly positive lower bounds and suppose that there exist a,b[a,b], a<a<b<b such that f(x)=g(x) for all x[a,b]. Then for all c,d[a,b] satisfying a<c<d<b , there exists n~N which depends only on f,g,a,b,a,bc,d , such that Ln(M)(f)(x)=Ln(M)(g)(x) for all x[c,d] and nN, nn~.

Proof â–¼
We obtain the desired conclusion as a direct consequence of the previous theorem. Indeed, firstly to make a distinction, we denote with Ln(M) the Lagrange max-product operator attached to functions defined on the interval [0,1]. In addition, in what follows, for all all nN and k{0,1,...,n} we denote with ln,k1 the fundamental Lagrange polynomials defined on the interval [0,1].

Suppose now that for the two functions f,gC([a,b]) we have f(x)=g(x), for all x[a,b]. Let us define the function h:[0,1][a,b], h(y)=a+(ba)y. It is immediate that for any x[a,b] there exists an unique y(x)=h1(x)[0,1] such that f(x)=(fh)(y(x)) and g(x)=(gh)(y(x)).

Then we observe that for any x[a,b] we have

ln,k(x)=(ba)nln,k1(y(x)), nNk{0,1,...,n}.

The above equalities imply

Ln(M)(f)(x)=k=0nln,k(x)f(xn,k)k=0nln,k(x)=(ba)nk=0nln,k1(y(x))(fh)(kn)(ba)nk=0nln,k1(y(x))=Ln(M)(fh)(y(x))

and analogously Ln(M)(g)(x)=Ln(M)(gh)(y(x)), for all x[a,b].

Then, our result is immediate by applying Theorem 2.1 to Ln(M)(gh)(y(x)) and Ln(M)(gh)(y(x)), where recall that fh,gh:[0,1][0,+), y(x)=h1(x) and h:[0,1][a,b], h(x)=a+(ba)x.

The next direct approximation result is now an immediate consequence of the localization result in Theorem 2.2, as follows.

Corollary 2.3

Let f:[a,b][0,) (a<b) be bounded on [a,b] with strictly positive lower bound and suppose that there exist a,b[a,b], a<a<b<b and the constant C0 which depends only on a,b,a,b, such that

|f(x)f(y)|C0|xy|,for all x[a,b],
2.1

that is f|[a,b]Lip[a,b]. Then, for any c,d[a,b] satisfying a<c<d<b, we have

|Ln(M)(f)(x)f(x)|Cn,for all nN and x[c,d],

where the constant C depends only on f and the values a,b,a,b,c,d.

Proof â–¼
Let us define the function F:[a,b]R,
F(x)={f(x),ifx[a,b],f(a),ifx[a,a],f(b),ifx[b,b].

The hypothesis imply that F is continuous and strictly positive on [a,b] and according to Corollary 3.2 in [ 11 ] it results that

|Ln(M)(F)(x)F(x)|2ω1(F,ban)[a,b], x[a,b],nN,

Since by the definition of F we have ω1(F,ban)[a,b]ω1(f,ban)[a,b] and since by the relation (2.1) it easily follows ω1(f,ban)[a,b]C0(ba)/n, we get

|Ln(M)(F)(x)F(x)|2C0(ba)/n, x[a,b],nN.

Now, let us choose arbitrary c,d[a,b] such that a<c<d<b. Then, by Theorem 2.2 (applicable to f and F) it results the existence of n~N which depends only on a,b,a,bc,d,f,F such that Ln(M)(F)(x)=Ln(M)(f)(x) for all x[c,d]. But since actually the function F depends on the function f, it is clear that in fact n~ depends only on a,b,a,bc,d and f.

Therefore, for arbitrary x[c,d] and nN with nn~ we obtain

|Ln(M)(f)(x)f(x)|=|Ln(M)(F)(x)F(x)|2C0(ba)/n,

where C0 and n~ depend only on a,b,a,bc,d and f.

Now, denoting

C1=max1n<n~{nLn(M)(f)f[c,d]},

we finally obtain

|Ln(M)(f)(x)f(x)|Cn, for all nN,x[c,d],

with C=max{2C0(ba),C1} depending only on a,b,c,d and f. This proves the corollary. At the end of this section, as a consequence of the localization result in Theorem 2.2 we present a locally constant preserving property.

Corollary 2.4

Let f:[a,b][0,) be bounded on [a,b] with strictly positive lower bound and suppose that there exists a,b[a,b], a<a<b<b such that f is constant on [a,b] with the constant value α>0. Then for any c,d[a,b] with a<c<d<b, there exists n~N which depends only on a,b,a,b,c,d and f such that Ln(M)(f)(x)=α for all x[c,d] and nN, nn~.

Proof â–¼
Let g:[a,b]R+ be given by g(x)=α>0 for all x[a,b]. Since f(x)=g(x) for all x[a,b] and since obviously Ln(M)(g)(x)=α for all x[a,b], by Theorem 2.2 we easily obtain the desired conclusion.

3 Final Remarks

Let us note that in Hermann-Vértesi [ 15 ] , starting from a Lagrange interpolatory process (convergent or not)

Pn(f)(x)=k=0npn,k(x)f(xn,k),

with

pn,k(x)=(xxn,0)...(xxn,k1)(xxn,k+1)...(xxn,n)(xn,kxn,0)...(xn,kxn,k1)(xn,kxn,k+1)...(xn,kxn,n),

new linear interpolatory (rational) operators of the form

Rn(f)(x)=k=0nf(xn,k)|pn,k(x)|rk=0nf(xn,k)|pn,k(x)|r,

are constructed, for which in the case when r>2 and the knots xn,k satisfy some special requirements (e.g. some Jacobi knots), the Jackson-type order of approximation

Rn(f)fCω1(f;1/n),

is obtained (see e.g. Theorem 3.2 in Hermann-Vértesi [ 15 ] ).

In other words, for the linear rational construction Rn(f)(x), we get the same order of approximation as for the interpolatory max-product operator (which is piecewise rational)

Ln(M)(f)(x)=k=0nln,k(x)f(xn,k)k=0nln,k(x).

But clearly that with respect to Rn(f)(x), the max-product rational operator Ln(M)(f)(x) presents several advantages, pointed out by the next remarks.

Remark 3.1

For positive continuous functions, it provides an estimate in terms of ω1(f;1/n) for the simplest systems of knots (that is for the equidistant nodes). But, in fact, as it was mentioned in the last Remark in the paper [ 7 ] , the estimate holds for any kind of interpolatory systems of points with the property that the distance between two consecutive nodes converges to zero as n. It is worth noting that the operator Rn(f)(x) provides the same Jackson-type estimate, but for systems of interpolatory points satisfying additional requirements (e.g. Theorem 3.2 in [ 15 ] for the Jacobi knots).â–¡

Remark 3.2

In our best knowledge, the strong localization results in Theorem 2.1 and Corollary 2.4, have not equivalence for Rn(f)(x).â–¡

Remark 3.3

Although the expression of Ln(M)(f)(x) theoretically looks more complicated than that of Rn(f)(x), however from practical/computational point of view, there not exists any difference between the usage of computer softwares (like Matlab or Mathematica) to trace the graphs of Ln(M)(f)(x) and Rn(f)(x), for any concrete choices of f. In fact, in Computer Science, the sum () operation and the maximum () operator have similar levels of computability.â–¡

Remark 3.4

The results in Theorem 2.2 and Corollary 2.4 show the nice property of the max-product interpolation operator Ln(M) to reproduce locally with great accuracy the graph of a strictly positive non-smooth continuous function f. For example, Corollary 2.4 shows that Ln(M)(f) is very suitable to approximate continuous functions which are strictly positive constants on some subintervals. Namely, if f is a continuous strictly positive function which is constant on some subintervals [αi,βi], i=1,...,p, of [a,b], then for sufficiently large n, Ln(M)(f) takes the same constant values on subintervals sufficiently close to each [αi,βi], i=1,...,p.â–¡

Remark 3.5

It is easy to see that the results expressed by Theorem 2.2 and Corollaries 2.3-2.4 can be extended to bounded functions of variable sign, for the new max-product operators of the form L¯n(M)(f)(x)=Ln(M)(f+c¯)(x)c¯, where c¯ is a constant such that f(x)+c¯>0, for all x[a,b]. Note that, for example in the case of Theorem 2.2, for f and g bounded and of variable sign, evidently that we may choose a constant c¯ such that f(x)+c¯>0 and g(x)+c¯>0, for all x[a,b].â–¡

Bibliography

1

B. Bede and S.G. Gal, Approximation by nonlinear Bernstein and Favard-Szász-Mirakjan operators of max-product kind, J. Concr. Applicable Math., 8 (2010), no. 2, pp.  193–207.

2

B. Bede, L. Coroianu and S.G. Gal, Approximation and shape preserving properties of the Bernstein operator of max-product kind, Intern. J. Math. Math. Sci., 2009 (2009), Article ID 590589, 26 pages. doi:10.1155/2009/590589.

3

B. Bede, L. Coroianu and S.G. Gal, Approximation by truncated Favard-Szász-Mirakjan operator of max-product kind, Demonstratio Mathematica, XLIV (2011), no. 1, pp.  105–122.

4

B. Bede, L. Coroianu and S.G. Gal, Approximation and shape preserving properties of the nonlinear Bleimann-Butzer-Hahn operators of max-product kind, Comment. Math. Univ. Carol., 51 (2010), no. 3, pp.  397–415.

5

B. Bede, L. Coroianu and S.G. Gal, Approximation and shape preserving properties of the nonlinear Meyer-Konig and Zeller operator of max-product kind, Numer. Funct. Anal. Optim., 31 (2010), no. 3, pp.  232–25.

6

B. Bede, L. Coroianu and S.G. Gal, Approximation and shape preserving properties of the truncated Baskakov operator of max-product kind, Revista Union Mat. Argentina, 52 (2011), no. 1, pp.  89–107.

7

B. Bede, L. Coroianu and S.G. Gal, Approximation and shape preserving properties of the nonlinear Baskakov operator of max-product kind, Studia Univ. Babes-Bolyai, Ser. Math., LV (2010), pp.  193–218.

8

S. Cobzaş and I. Muntean, Condensation of singularities and divergence results in approximation theory, J. Approx. Theory, 31 (1980), no. 2, pp.  138–153.

9

L. Coroianu and S.G. Gal, Approximation by nonlinear Lagrange interpolation operators of max-product kind on Chebyshev knots of second kind, J. Comp. Anal. Appl., 13 (2010), no. 2, pp.  211–224.

10

L. Coroianu and S.G. Gal, Approximation by nonlinear Hermite-Fejér interpolation operators of max-product kind on Chebyshev nodes, Rev. Anal. Numér. Théor. Approx. (Cluj), 39 (2010), no. 1, pp.  29–39. \includegraphics[scale=0.1]{ext-link.png}

11

L. Coroianu and S.G. Gal, Approximation by max-product Lagrange interpolation operators, Studia Univ. “Babes-Bolyai", Ser. Math., LVI (2011), no. 2, pp.  1–11.

12

L. Coroianu and S.G. Gal, Classes of functions with improved estimates in approximation by the max-product Bernstein operator, Anal. Appl. (Singapore), 9 (2011), no. 3, pp.  249–274.

13

L. Coroianu and S.G. Gal, Saturation results for the Lagrange max-product interpolation operator based on equidistant knots, Rev. Anal. Numér. Théor. Approx. (Cluj), 41 (2012), no. 1, pp.  27–41. \includegraphics[scale=0.1]{ext-link.png}

14

S.G. Gal, Shape-Preserving Approximation by Real and Complex Polynomials, Birkhäuser, Boston-Basel-Berlin, 2008.

15

T. Hermann and P. Vértesi, On the method of Somorjai, Acta Math. Hung., 54 (1989) no. (3-4), pp.  253–262.

16

I. Muntean, The Lagrange interpolation operators are densely divergent, Studia Univ. “Babes-Bolyai" (Cluj), Ser. Math., 21 (1976), pp.  28–30.

17

J. Szabados and P. Vértesi, Interpolation of Functions, World Scientific, Singapore, New Jersey, London, Hong Kong, 1990.