In this paper, the authors take the first step in the study of constructive methods by using Sobolev polynomials.
To do that, the authors use the connection formulas between Sobolev polynomials and classical Laguerre polynomials, as well as the well-known Fourier coefficients for these latter.
Then, the authors compute explicit formulas for the Fourier coefficients of some families of Laguerre–Sobolev type orthogonal polynomials over a finite interval. The authors also describe an oscillatory region in each case as a reasonable choice for approximation purposes.
In order to take the first step in the study of constructive methods by using Sobolev polynomials, this paper deals with Fourier coefficients for certain families of polynomials orthogonal with respect to the Sobolev type inner product. As far as the authors know, this particular problem has not been addressed in the existing literature.
1. Introduction
Within the framework of spectral approximation, and to recover values of smooth functions with exponential accurate, it is customary to use Fourier series for periodic problems and series of classical orthogonal polynomials for nonperiodic problems. Nevertheless, if it deals with piecewise smooth function, estimates by means of partial sums are unhealthy; oscillations do not decrease near discontinuities with partial sums of higher order; and far of them, convergence order is low. Thus, the global properties from Fourier coefficients are not enough to obtain local information. This lack of uniform convergence is known as Gibbs phenomenon. A priori, this is a serious issue considering the large number of applications modeled through piecewise smooth function. In literature, methods to face the Gibbs phenomenon in reconstruction of piecewise smooth functions from partial sums have been widely studied. For instance, in Refs. [1, 2], the problem to construct piecewise smooth function values with exponential accuracy at all points is solved by means of approximations with Fourier–Gegenbauer coefficients expansions. These are the so-called Gegenbauer reconstruction methods where the expansion of Gegenbauer polynomials in its Fourier series is crucial. In Ref. [3], the Gegenbauer reconstruction methods are revisited and analyzed in order to prove that Gegenbauer reconstruction is also effective for Fourier–Bessel series. To do that, the author obtains coefficients Fourier for Jacobi polynomials and also for classical orthogonal polynomials with unbounded support (Laguerre, Hermite).
On the other hand, consider a vector of Borel positive measures (μ0, μ1, …, μm), on the real line, with finite moments and μ0 with continuous support. Then, we define the Sobolev inner product on the space of polynomials with real coefficients.
A sequence of polynomials , deg Sn = n, is orthogonal with respect to (1.1) if
The sequence is said to be a sequence of Sobolev polynomials orthogonal with respect to (1.1). If μk is discrete, for k = 1, …, m, the above inner product and the sequence are said to be of Sobolev type. Sobolev orthogonal polynomials have been widely studied in the last three decades. The first publication on Sobolev polynomials goes back to 1962 in Ref. [4], which deals with certain extremal problem related to smooth polynomial approximation whose solution is posed by means of Sobolev–Legendre polynomials. Such a problem is formulated previously in Ref. [5], although not in terms of Sobolev orthogonality. It has been documented as the approximations with Sobolev–Fourier series from smooth functions in the corresponding Sobolev space improve approximations made through standard families of orthogonal polynomials (see Ref. [6]). Additional applications include spectral methods in numerical analysis for ordinary differential equations and partial differential equations, and generalization of Gauss quadrature formulas, among others. The nice surveys [7, 8] are highly recommended, as well as the paper [9] and references therein. In order to take the first step in the study of constructive methods by using Sobolev polynomials, this paper deals with Fourier coefficients for certain families of polynomials orthogonal with respect to the Sobolev type inner product (1.1) when μ0 is the classical and absolutely continuous Laguerre measure on [0, ∞). In the next section, we propose the basic background with respect to Laguerre polynomials, and we present the particular Sobolev–Laguerre type families of polynomials to be discussed. In Section 3, we obtain the respective Fourier coefficients by using of similar techniques as the presented in Ref. [10]. Since the orthogonality interval for Laguerre polynomials is unbounded, we will turn special attention to oscillation regions for the Sobolev polynomials.
2. Preliminaries
Let be the space of polynomials with real coefficients
2.1 Classical Laguerre polynomials and generalities
The classical Laguerre polynomials , with α > − 1, are orthogonal with respect to the inner product:
For an arbitrary polynomial p, k(p) will denote the leading coefficient of p. In the sequel, to normalize Laguerre polynomials, we assume that . These polynomials satisfy the three terms recurrence relation (TTRR in short),
for n ≥ 0 with the initial conditions and . For n ≥ 1, the zeros of every are all real, simple and are located in (0, ∞) (see Ref. [11]). In the sequel, will denote the zeros of ordered in increasing order.
Let p be a polynomial with real zeros. An oscillatory region I for p is any bounded interval containing their zeros, in such a way that p is monotone outside I.
With respect to an oscillatory region of classical Laguerre polynomials, we get the next.
We consider, for a nonnegative integer m, the functions defined as (see Ref. [10]),
where em is the m − th partial sum of the Maclaurin series for the exponential function and [a, b] is a bounded interval.
As a consequence of this definition, it is possible to show that if x = ɛξ + δ, and , we get
In this way, the next result for the Fourier series for Laguerre polynomials is presented in Ref. [10].
Let [a, b] be an interval with − ∞ < a < b < ∞ and ξ ∈ [ − 1, 1], , . The Fourier coefficients for , in the local variable ξ, are given by
2.2 Quasi-orthogonality and zeros
Let be a sequence of polynomials orthogonal with respect to a positive Borel measure μ supported on [a, b], with − ∞ ≤ a < b ≤ ∞, i.e.
Let r be a nonnegative integer and Rn a polynomial with degree n ≥ r satisfying for k = 0, 1, 2, …, n − r − 1, and . Then, Rn is said to be quasi-orthogonal of order r on [a, b] and with respect to μ.
Of course, if r = 0, then the orthogonality is recovered. The next result describes a necessary and sufficient condition for quasi-orthogonality.
([12]). Rn is quasi-orthogonal of order r on [a, b] with respect to μ if and only if there exist numbers bn,i, i = 0, 1, …, r, with bn,0bn,r ≠ 0, such that
With respect to zeros of quasi-orthogonal polynomials, the next result is well known.
([12]). If Rn is quasi-orthogonal of order r with respect to μ on [a, b], then Rn has n − r simple zeros on (a, b).
Suppose that Rn is quasi-orthogonal of order r with respect to μ on [a, b] and Rn and Pn are monic. It is well known that the monic orthogonal polynomials can be obtained by means of a TTRR:
and we define and . In the particular case, when r = 2, from (2.5), we get Rn(x) = bn,0Pn(x) + bn,1Pn−1(x) + bn,2Pn−2(x), and we also define and . The next results refer to behavior and localization of zeros of quasi-orthogonal polynomials for r = 2.
([13]). If bn ≤ Cn, the n zeros of Rn are real and simples.
([13]). Suppose that and are the zeros of Pn and Rn, respectively, and ordered in increasing order.
bn < Cn if and only if
(2.7)0 < bn < Cn and if and only if
(2.8)0 < bn < Cn and if and only if
(2.9)
For i = 1, 2, …, n
if and only if
2.3 Laguerre–Sobolev type orthogonal polynomials, nondiagonal case
If and P(x)≔(p(x), p′(x))t, we define the Laguerre–Sobolev type inner product
where , with M0M1 ≥ 0 and λ such that det A ≥ 0. Let be the sequence of polynomials orthogonal with respect to (2.12) such that for n ≥ 0.
2.3.1 Laguerre–Sobolev type polynomials of higher order derivatives
Let be orthogonal with respect to Sobolev inner product
with W > 0, m a nonnegative integer and . Moreover .
With respect to zeros of every , we enunciate the next results.
(See Ref. [16]). For every n, the zeros of are real, simple and at most one of them is outside (0, ∞). If has a zero in (−∞, 0] then n ≥ m + 1. In addition, if for n0, has a negative zero, then has a negative zero for n > n0.
Assume that n ≥ m + 1. If are the zeros of , ordered in increasing order, then υn,i < xn,i for i = 1, …, n.
If ρn is the negative zero of then , where denotes the m − th positive zero of .
2.3.2 Christoffel transformations and Laguerre–Sobolev type inner product with mass outside support
Given ξ ≤ 0, and an integer k ≥ 1, we consider the weight ωα,k(x) = (x − ξ)ke−xxα, on [0, ∞). This is a Christoffel perturbation of the classical Laguerre measure (see Ref. [11]). denotes the respective sequence of orthogonal polynomials, where for every n and . An algebraic connection between polynomials orthogonal with respect to the weight ωα,k(x) is as follows (see Ref. [11]),
Assume are the zeros of in increasing order, with .
Now we consider the Sobolev–Laguerre type inner product:
with M, N ≥ 0, ξ ≤ 0. Let be the respective sequence of orthogonal polynomials such that for n ≥ 0.
([17]). There exist constants Dn,0, Dn,1 and Dn,2 such that
If M, N > 0, then , and .
If M = 0 and N > 0, then , and .
If N = 0 and M > 0, then , and Dn,2 = 0.
Let be the zeros of in increasing order. To describe results on zeros of every , we present the next results.
([17]). The zeros of are real, simple and at most one of them is outside [ξ, ∞).
3. Fourier coefficients for Laguerre–Sobolev type polynomials
In this section, we describe the Fourier coefficients associated to Laguerre–Sobolev type polynomials presented in the above section, computed on any finite interval [a, b]. For approximation purposes, we will find an oscillatory region for every family of Sobolev–Laguerre polynomials, in order to exhibit a reasonable choose for the interval [a, b].
3.1 Nondiagonal case
Let be the sequence of Sobolev polynomials orthogonal with respect to the inner product (2.12). From (2.13), this sequence is quasi-orthogonal of order 2 with respect to the classical Laguerre polynomials with parameter α + 2. Then, we consider (2.1) for α + 2, and from (2.13) we define
and
Then, in the language of Theorem 3, we get the next result.
If
Finally, since the zeros of and are interlaced, we obtain
The above inequalities imply that
From the Proposition 1, we get the result. □
On the other hand, we suppose that x ∈ [a, b], and we make the transformation x = ɛξ + δ, where ξ ∈ [ − 1, 1], and Then
and by using of (2.13) we get
and by using of (2.4) we obtain
We summarize in the next.
Let [a, b] be a bounded interval. Assume x in [a, b], and x = ɛξ + δ, where ξ ∈ [ − 1, 1], and The coefficients of Fourier for , in the local variable ξ, are giving by
3.2 Higher order derivatives
represents the sequence of polynomials orthogonal with respect to (2.14). As before, we propose a bounded interval that containing the n zeros of for n large enough.
For n ≥ m + 1, the zeros of are located in [ − mxn,m+1, ζn,α].
As before, we assume x ∈ [a, b] and x = ɛξ + δ, where ξ ∈ [ − 1, 1], and
From (2.15) we have
and from (2.4) we arrive to the next.
Let [a, b] a bounded interval and n ≥ m + 1. Consider the transformation x = ɛξ + δ, where x ∈ [a, b], ξ ∈ [ − 1, 1], and . The Fourier series for , with the local variable ξ, is given by
3.3 Mass outside support
Assume that
with
or equivalently
then we get the equations
and
If
If ξ = 0 then
and forj = 1, …, n
According to (2.17), we can deduce that
and
and as consequence we get the next.
The zeros of are located in [xn,1, βn + k,α].
Since the Fourier series for in the local variable η is determined by the coefficients
by using of Lemma 1, (2.2) and (2.3) we have the next.
Fourier coefficients for on a finite interval [a, b], in the local variable η, with x = ɛη + δ, η ∈ [ − 1, 1], and , are defined by
Let be the sequence of polynomials orthogonal with respect to (2.18).
For n ≥ 2, the zeros of are into [2ξ − xn+1,3, xn+1,n+1].
In the same way, and , thus
□
On the one hand, and on a finite interval [a, b], we compute the Fourier coefficients for every and in terms of the local variable η. Indeed, if
by using of (2.19)
Then, from (3.5), for k = 2 and k = 4, we obtain
and
respectively.
As a consequence, (2.19) can be written as
where
Then, Fourier coefficients are given by
and if we use (2.4) and (3.6) we get
We summarize in the next.
The authors thank the referees for the careful revision of the manuscript. Their suggestions have contributed to improving the presentation. This work has been partially supported by Dirección de Investigaciones of Universidad Pedagógica y Tecnológica de Colombia.
