Next Article in Journal
Modified Mann-Type Subgradient Extragradient Rules for Variational Inequalities and Common Fixed Points Implicating Countably Many Nonexpansive Operators
Previous Article in Journal
Zeroing Neural Network Approaches Based on Direct and Indirect Methods for Solving the Yang–Baxter-like Matrix Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Polynomials Orthogonal with Respect to an Inner Product Involving Higher-Order Differences: The Meixner Case

by
Roberto S. Costas-Santos
1,
Anier Soria-Lorente
1,* and
Jean-Marie Vilaire
2
1
Facultad de Ciencias Técnicas, Universidad de Granma, Km. 17.5 de la carretera de Bayano-Manzanillo, Bayamo 85100, Cuba
2
Institut des Sciences, des Technologies et des Études Avancées d’Haïti, #10, Rue Mercier-Laham, Delmas 60, Musseau, Port-au-Prince 15953, Haiti
*
Author to whom correspondence should be addressed.
Submission received: 15 April 2022 / Revised: 16 May 2022 / Accepted: 3 June 2022 / Published: 6 June 2022
(This article belongs to the Section Difference and Differential Equations)

Abstract

:
In this contribution we consider sequences of monic polynomials orthogonal with respect to the Sobolev-type inner product f , g = u M , f g + λ T j f ( α ) T j g ( α ) , where u M is the Meixner linear operator, λ R + , j N , α 0 , and T is the forward difference operator Δ or the backward difference operator ∇. Moreover, we derive an explicit representation for these polynomials. The ladder operators associated with these polynomials are obtained, and the linear difference equation of the second order is also given. In addition, for these polynomials, we derive a ( 2 j + 3 ) -term recurrence relation. Finally, we find the Mehler–Heine type formula for the particular case α = 0 .

1. Introduction

Meixner orthogonal polynomials, usually denoted in the literature as M n ( x ; β , c ) , constitute a family of classical orthogonal polynomials introduced by J. Meixner in 1934 in his seminal paper [1]. When β > 0 and 0 < c < 1 , they are orthogonal with respect to the well-known negative binomial distribution of the probability theory, i.e., in such a case:
x = 0 x + β 1 x c x M n ( x ; β , c ) M m ( x ; β , c ) = 0 , m n , m , n = 0 , 1 , ,
so the Meixner linear functional is
u M = x = 0 x + β 1 x c x δ x .
Therefore, they are orthogonal on the uniform lattice in the interval [ 0 , + ) , and they satisfy a hypergeometric-type difference equation on the aforesaid uniform lattice. Due to its classical character, their finite differences also constitute an orthogonal polynomial family, their corresponding orthogonality weights satisfy a Pearson-type difference equation, and they even satisfy two different kinds of structure relations. A distinctive and interesting characteristic of this family is that they have certain dual character, that is, every formula one can derive for M n ( x ; β , c ) has a dual formula with x and n interchanged:
c m n n ! ( 1 + β ) m 1 M n ( m ; β , c ) = c n m m ! ( 1 + β ) n 1 M m ( n ; β , c ) .
(See, for example [2,3,4] and the references therein).
On the other hand, since the first paper [5] on Sobolev’s orthogonal polynomials was published by Althammer, the results connected to these polynomials have attracted the attention of several mathematicians. The name of the Sobolev orthogonal polynomials is given to those families of polynomials orthogonal with respect to inner products involving positive Borel measures supported on infinite subsets of the real lines and also involving regular derivatives. Moreover, in the case that the derivatives appear only on function evaluations on a finite discrete set, the corresponding families are called Sobolev-type or discrete Sobolev orthogonal polynomial sequences. For a recent and comprehensive survey on the subject, see [6] and the references therein. In the last decade of the past century, H. Bavinck introduced the study of orthogonal polynomials with respect to the inner product involving differences instead of derivatives:
f , g λ = R f ( x ) g ( x ) d ψ ( x ) + λ ( Δ f ) ( d ) ( Δ g ) ( d ) ,
where λ R + , d R , and ψ is a distribution function with infinite spectrum, see [7,8]. Moreover, in these works Bavinck obtained algebraic properties and some results connected to the location of the zeros of the orthogonal polynomials with respect to the inner product (1). On the other hand, in [8], he proved that the orthogonal polynomials with respect to inner product defined in Equation (1) satisfy a five-term recurrence relation. Furthermore, in [9], the authors considered the following inner product:
f , g = ( 1 c ) β x = 0 f ( x ) g ( x ) c x Γ ( β + x ) Γ ( β ) Γ ( x + 1 ) + λ f ( 0 ) g ( 0 ) ,
where β > 0 , 0 < c < 1 , λ > 0 , and P denotes the linear space of all polynomials with real coefficients. Here, they obtained a second-order difference equation with polynomial coefficients, which the orthogonal polynomials with respect to (2) satisfy. Then, in [10], the authors showed that the Sobolev-type Meixner polynomials orthogonal with respect to the inner product are eigenfunctions of a difference operator:
f , g = ( 1 c ) β x = 0 f ( x ) g ( x ) c x Γ ( β + x ) Γ ( β ) Γ ( x + 1 ) + M f ( 0 ) g ( 0 ) + N ( Δ f ) ( 0 ) ( Δ g ) ( 0 ) ,
where β > 0 , 0 < c < 1 , and M , N 0 . Other slightly more recent results, connected with the Sobolev–Meixner polynomials, can be found in [11,12]. Note that our work is in someways connected to the paper [13], where the authors consider the Charlier case.
The structure of the paper is the following: In Section 2, we introduce some preliminary results about Meixner polynomials, which will be very useful in the analysis presented. In Section 3, we obtain the connection formula between the Meixner polynomials and the polynomials orthogonal with respect to the Sobolev-type inner product:
f , g = u M , f g + λ T j f ( α ) T j g ( α ) ,
where u M is the Meixner linear operator, λ R + , j N , α 0 , and T is the forward or the backward difference operator. In addition, we deduce the hypergeometric representation of such polynomials. In Section 4, we find the ladder (creation and annihilation) operators for the sequence of Sobolev-type orthogonal polynomials. As a consequence, the second-order linear difference equation associated with them is deduced. On the other hand, in Section 5, we determine the ( 2 j + 3 ) -term recurrence relation that these polynomials satisfy. Finally, in Section 6, we determine the Mehler–Heine-type formula for the special case of α = 0 . Indeed, the techniques used in Section 3, Section 4 and Section 5 are based on those used in [3,14,15], respectively.

2. Preliminaries

We adopt the following set of notations: N 0 : = { 0 } N = { 0 , 1 , 2 , 3 , } , and we use the sets Z , R , and C , which represent the integers, real numbers, and complex numbers, respectively. P denotes the vector space of univariate, complex-valued polynomials, and let P denote its algebraic dual space.
We also adopt the following notation and conventions. We denote by u , p the duality bracket for u P and p P .
Definition 1.
For u P , π P , and c C , let π u , ( x c ) 1 u , u , and Δ u be the linear functionals defined by:
π u , p : = u , π p , p P , ( x c ) 1 u , p : = u , p ( x ) p ( c ) x c , p P , u , p : = u , Δ p ,
and thus: Δ u , p : = u , p , whereand Δ are the backward and forward difference operators, defined as:
Δ f ( x ) : = f ( x + 1 ) f ( x ) , f ( x ) : = f ( x ) f ( x 1 ) .
The Dirac delta functional, δ c , is the functional defined by δ c , p : = p ( c ) , p P .
In order to obtain our derived identities, we rely on properties of the Pochhammer symbol (shifted factorial). For any n N 0 , a C , the Pochhammer symbol is defined as:
( a ) n : = ( a ) ( a + 1 ) ( a + n 1 ) , n N 0 .
Furthermore, define for all a , b C :
( a ) b : = Γ ( a + b ) Γ ( a ) ,
where a + b N 0 , and we will also use the common notational product convention
( a 1 , , a k ) b : = ( a 1 ) b ( a k ) b .
The hypergeometric series, which we will use often, is defined for z C such that | z | < 1 , s , r N 0 , with s r 1 , as [16]:
r F s a 1 , , a r b 1 , , b s ; z : = k = 0 ( a 1 , , a r ) k ( b 1 , , b s ) k z k k ! ,
wherethe parameters must be chosen such that the denominator factors in the terms of the series are never zero.

The Meixner Polynomials

Let β and c be two complex numbers, such that c 0 , 1 and β is not a negative integer. We write { M n ( x ; β , c ) } n 0 for the sequence of monic Meixner polynomials defined by [2]:
M n ( x ; β , c ) = c n n ! ( 1 c ) n j = 0 n c j x j x β n j .
These polynomials are orthogonal with respect to the linear functional u M P , which is a classical functional since it fulfills the Pearson difference equation:
Δ x u M = x ( c 1 ) + β c u M ,
which is equivalent to the Pearson difference equation:
c ( x + β ) u M = x ( c 1 ) + β c u M .
Remark 1.
Observe that when β > 0 and 0 < c < 1 , then:
u M , f : = x = 0 x + β 1 x c x f ( x ) ,
which is a positive definite linear functional. Note that this definition can be extended to | c | < 1 , and β C is not a negative integer. Moreover, since
M n ( x ; β , c ) = ( 1 ) n M n ( x β ; β , c 1 ) ,
then one can extend the Meixner functional to | c | > 1 . In [17], the authors obtained an integral representation for this operator for β , c C , with c [ 0 , + ) and β N :
u M , f = C f ( z ) Γ ( z ) Γ ( β + z ) ( c ) z d z ,
where C is a complex contour from i to i separating the increasing poles { 0 , 1 , 2 , } from the decreasing poles { β , β 1 , β 2 , } .
When c 0 , 1 , they satisfy the following three-term recurrence formula:
x M n ( x ; β , c ) = M n + 1 ( x ; β , c ) + ( c + 1 ) n + β c 1 c M n ( x ; β , c ) + γ n M n 1 ( x ; β , c ) ,
where
γ n = c n ( n + β 1 ) ( 1 c ) 2 ,
which can be given explicitly in terms of hypergeometric series as [16]:
M n ( x ; β , c ) = c n ( β ) n ( c 1 ) n 2 F 1 n , x β ; 1 1 c .
Next, we summarize some basic properties of Meixner orthogonal polynomials to be used in the sequel.
Proposition 1
(cf. Proposition 2.1 in [18]). Let β , c C such that β N and | c | < 1 . Then, the following identities hold true for the Meixner polynomials:
1. 
Second-order difference equation:
c ( x + β ) y n ( x + 1 ) ( x ( c + 1 ) + β c ) y n ( x ) + x y n ( x 1 ) = n ( c 1 ) y n ( x ) ;
2. 
Structure relations. For every n N :
( x + β ) Δ M n ( x ; β , c ) = n M n ( x ; β , c ) + n ( n + β 1 ) 1 c M n 1 ( x ; β , c ) ,
x M n ( x ; β , c ) = n M n ( x ; β , c ) + n c ( n + β 1 ) 1 c M n 1 ( x ; β , c ) ;
3. 
Squared norm. For every n N :
d n 2 = M n ( x ; β , c ) 2 = u M , M n 2 ( x ; β , c ) = ( β ) n c n n ! ( 1 c ) β + 2 n ;
4. 
The value at the left endpoint initial extreme:
M n ( 0 ; β , c ) = ( β ) n c n ( c 1 ) n ;
5. 
Forward and backward difference operators. For every n , k N :
Δ k M n ( x ; β , c ) = ( m k + 1 ) k M n k ( x ; β + k , c ) ,
where Δ k f ( x ) = Δ k 1 Δ f ( x ) for k = 1 , 2 , , and Δ 0 f ( x ) = f ( x ) ;
6. 
Mehler–Heine type formula (Equation (35) in [19]):
lim n ( 1 c ) n + β + z M n ( z ; β , c ) ( z n + 1 ) n = 1 , z C N .
To complete this section, we present some useful results we will need throughout the paper.
Proposition 2
(Christoffel-Darboux formula). Let p n n N 0 be a sequence of monic polynomials orthogonal with respect to the linear functional u . If we denote the n-th reproducing kernel by
K n ( x , y ) : = k = 0 n 1 p k ( x ) p k ( y ) u , p k 2 ,
then, for all n N and x y :
K n ( x , y ) = 1 u , p n 1 2 p n ( x ) p n 1 ( y ) p n ( y ) p n 1 ( x ) x y .
Taking into account the inner product we have considered, then it is natural to consider the partial derivatives of K n ( x , y ) . We will use the following notation:
K n , 1 ( i , j ) ( x , y ) : = k = 0 n 1 i p k ( x ) j p k ( y ) u , p k 2 ,
and
K n , 2 ( i , j ) ( x , y ) : = k = 0 n 1 Δ i p k ( x ) Δ j p k ( y ) u , p k 2 .
Remark 2.
Since this work is related to the Meixner polynomials, we will have to substitute p k ( x ) (resp. p k ( y ) ) in (16) by M k ( x ; β , c ) (resp. M k 1 ( x ; β , c ) ) and in all expressions that depend on it, such as (18) and (19).
Proposition 3.
The following identities hold:
K n , 1 ( 0 , j ) ( x , y ) = j ! d n 1 2 k = 0 j M n ( x ; β , c ) k M n 1 ( y ; β , c ) M n 1 ( x ; β , c ) k M n ( y ; β , c ) k ! ( x y + k ) j + 1 k ,
K n , 2 ( 0 , j ) ( x , y ) = j ! d n 1 2 k = 0 j M n ( x ; β , c ) Δ k M n 1 ( y ; β , c ) M n 1 ( x ; β , c ) Δ k M n ( y ; β , c ) k ! ( x y j ) j + 1 k ,
where d n 1 2 = u M , M n 1 2 ( x ; β , c ) .
Proof. 
We are going to prove the first identity. After applying the difference operator j to (17) with respect to y, we obtain:
K n , 1 ( 0 , j ) ( x , y ) = 1 d n 1 2 M n ( x ; β , c ) y j M n 1 ( x ; β , c ) x y M n 1 ( x ; β , c ) y j M n ( y ; β , c ) x y .
Using an analogue of the Leibnitz’s rule
n f ( x ) g ( x ) = k = 0 n n k k f ( x ) n k g ( x k ) ,
and since for any positive integer k we have
y k 1 x y = k ! ( x y ) k + 1 ,
the result follows after a straightforward calculation. The proof of the second identity is analogous, and it will be omitted. Hence, the result follows. □
Proposition 4.
The following identity holds for the kernel associated with the Meixner polynomials:
K n , 2 ( j , j ) ( 0 , 0 ) = j ! ( 1 c ) β + 2 j c j ( β ) j k = 0 n j 1 ( j + 1 ) k ( β + j ) k ( 1 ) k c k k ! .
Proof. 
Taking (13)–(16) into account, and, by definition, (18), we have:
K n , 2 ( j , j ) ( 0 , 0 ) = k = 0 n 1 ( k j + 1 ) j 2 ( β + j ) j 2 c 2 k 2 j ( 1 c ) β + 2 k ( c 1 ) 2 k 2 j ( β ) j c k k ! = ( 1 c ) 2 j + β c 2 j ( β ) j 2 k = 0 n 1 ( k j + 1 ) j 2 ( β ) k c k k ! .
By using some identities of the Pochhammer symbols, we obtain:
K n , 2 ( j , j ) ( 0 , 0 ) = ( 1 c ) 2 j + β c 2 j ( β ) j 2 k = 0 n 1 ( k j + 1 ) j 2 ( β ) k c k ( 1 ) k j 2 = ( 1 c ) 2 j + β c 2 j ( β ) j 2 k = 0 n 1 j ( 1 ) k + j ( β ) k + j c k + j ( 1 ) k 2 ,
from this expression, it is a direct calculation to obtain the desired expression. □
Corollary 1.
The following limit for the kernels associated to the Meixner polynomials holds:
lim n K n , 2 ( j , j ) ( 0 , 0 ) = j ! ( 1 c ) β + 2 j c j ( β ) j 2 F 1 1 + j , β + j 1 ; c .

3. The Sobolev-Type Meixner Polynomials

We start this section by introducing the Sobolev-type inner product (3):
f , g λ , j , = u M , f g + λ T j f ( α ) T j g ( α ) ,
where u M is the Meixner linear operator, j N , α 0 , and T is the operator ∇ when = 1 , and it is the operator Δ when = 2 .
We denote by M n j , ( x ; β , c ; λ ) n N 0 the sequence of monic polynomials, orthogonal with respect to the inner product (24). These polynomials are said to be Sobolev-type Meixner polynomials.

Connection Formula and Hypergeometric Representation

We first express the Sobolev-type Meixner polynomials in terms of the monic Meixner polynomials and the Kernel polynomials associated with the Meixner polynomials.
Taking into account the Fourier expansion and using the orthogonality conditions of ( M n ( x ) ) and ( M n j , ( x ; β , c ; λ ) ) , we obtain (see cf. Equation (2.8) in [20]):
M n j , ( x ; β , c ; λ ) = M n ( x ; β , c ) λ T j M n ( α ; β , c ) 1 + λ K n , ( j , j ) ( α , α ) K n , ( 0 , j ) ( x , α ) , = 1 , 2 .
We can express the Sobolev-type Meixner polynomials in terms of the Meixner and their associated Kernel polynomials. Moreover, starting from (25) and by using the recurrence relation of the Meixner polynomials (7), we have:
M n j , ( x ; β , c ; λ ) = A 1 , n j , ( x ) M n ( x ; β , c ) + B n , 1 j , ( x ) M n 1 ( x ; β , c ) ,
where:
A 1 , n j , 1 ( x ) = 1 λ T j M n ( α ; β , c ) 1 + λ K n , ( j , j ) ( α , α ) j ! d n 1 2 k = 0 j k M n 1 ( α ; β , c ) k ! ( x α + k ) j + 1 k ,
A 1 , n j , 2 ( x ) = 1 λ T j M n ( α ; β , c ) 1 + λ K n , ( j , j ) ( α , α ) j ! d n 1 2 k = 0 j Δ k M n 1 ( α ; β , c ) k ! ( x α j ) j + 1 k ,
B 1 , n j , 1 ( x ) = λ T j M n ( α ; β , c ) 1 + λ K n , ( j , j ) ( α , α ) j ! d n 1 2 k = 0 j k M n ( α ; β , c ) k ! ( x α + k ) j + 1 k ,
B 1 , n j , 2 ( x ) = λ T j M n ( α ; β , c ) 1 + λ K n , ( j , j ) ( α , α ) j ! d n 1 2 k = 0 j Δ k M n ( α ; β , c ) k ! ( x α j ) j + 1 k ,
where, in such a case, T f ( α ) means
T f ( α ) = T f ( x ) | x = α .
From these identities, we can express the Sobolev-type Meixner polynomials in terms of hypergeometric series.
Theorem 1.
The monic Sobolev-type Meixner polynomial M n j , ( x ; β , c ; λ ) has the following hypergeometric representation for = 1 , 2 :
M n j , ( x ; β , c ; λ ) = ( β ) n 1 c n 1 ( c 1 ) n 1 h n ( x ) 3 F 2 n , x , f n ( x ) β , f n ( x ) 1 ; 1 1 c ,
where f n ( x ) is given in (33) and
h n ( x ) = c ( β + n 1 ) 1 c A 1 , n j , ( x ) B 1 , n j , ( x ) .
Proof. 
Taking into account ( x ) k = 0 if x < k , as well as (8) and (26), we deduce:
M n j , ( x ; β , c ; λ ) = c n ( β ) n ( c 1 ) n A 1 , n j , ( x ) k = 0 n ( n ) k ( x ) k ( β ) k k ! 1 1 c k + c n 1 ( β ) n 1 ( c 1 ) n 1 B 1 , n j , ( x ) k = 0 n 1 ( 1 n ) k ( x ) k ( β ) k k ! 1 1 c k .
By using the identity:
( a + k ) ( a ) k = a ( a + 1 ) k ,
we obtain:
M n j , ( x ; β , c ; λ ) = c n ( β ) n ( c 1 ) n A 1 , n j , ( x ) k = 0 n ( n ) k ( x ) k ( β ) k k ! 1 1 c k + c n 1 ( β ) n 1 ( c 1 ) n 1 B 1 , n j , ( x ) n k = 0 n ( n k ) ( n ) k ( x ) k ( β ) k k ! 1 1 c k .
Thus, we have:
M n j , ( x ; β , c ; λ ) = ( β ) n 1 c n 1 ( c 1 ) n 1 B 1 , n j , ( x ) n k = 0 n ( f n ( x ) k + 1 ) ( n ) k ( x ) k ( β ) k k ! 1 1 c k ,
where
f n ( x ) = n 1 n c ( β + n 1 ) A 1 , n j , ( x ) ( 1 c ) B 1 , n j , ( x ) .
In addition, after a straightforward calculation and by using (32) with a f n ( x ) 1 , the identity (31) follows. This completes the proof. □

4. Second-Order Linear Difference Equation

In this section, we obtain a second-order linear difference equation that the sequence { M n j , ( x ; β , c ; λ ) } n 0 satisfies. In order to achieve this, we will find the ladder (creation and annihilation) operators, using the connection Formula (26), the three-term recurrence relation (7), and the structure relations (10) and (11) satisfied by them.
From (26) and the recurrence relation (7), we deduce the following result:
M n 1 j , ( x ; β , c ; λ ) = A 2 , n j , ( x ) M n ( x ; β , c ) + B 2 , n j , ( x ) M n 1 ( x ; β , c ) ,
where
A 2 , n j , ( x ) = ( c 1 ) B 1 , n 1 j , ( x ) ( c + 1 ) ( n 1 ) + β c , and B 2 , n j , ( x ) = A 1 , n 1 j , ( x ) + A 2 , n j , ( x ) ( 1 x ) .
Applying the operator T to (26) and by using (23), we have:
T M n j , ( x ; β , c ; λ ) = M n ( x ; β , c ) T A 1 , n j , ( x ) + A 1 , n j , ( x + ( 1 ) ) T M n ( x ; β , c ) + M n 1 ( x ; β , c ) T B 1 , n j , ( x ) + A 1 , n j , ( x + ( 1 ) ) T M n 1 ( x ; β , c ) .
Then, multiplying the previous expression by x and using the structure relation (11) if = 1 and x + β , using the structure relation (10) if = 2 , as well as using the recurrence relation (7), we deduce the following expressions:
x M n j , 1 ( x ; β , c ; λ ) = C 1 , n 1 ( x ) M n ( x ; β , c ) + D 1 , n 1 ( x ) M n 1 ( x ; β , c ) ,
( x + β ) Δ M n j , 2 ( x ; β , c ; λ ) = C 1 , n 2 ( x ) M n ( x ; β , c ) + D 1 , n 2 ( x ) M n 1 ( x ; β , c ) ,
x M n 1 j , 1 ( x ; β , c ; λ ) = C 2 , n 1 ( x ) M n ( x ; β , c ) + D 2 , n 1 ( x ) M n 1 ( x ; β , c ) ,
and
( x + β ) Δ M n 1 j , 2 ( x ; β , c ; λ ) = C 2 , n 2 ( x ) M n ( x ; β , c ) + D 2 , n 2 ( x ) M n 1 ( x ; β , c ) ,
respectively, where all the coefficients can be computed explicitly. Moreover, from (26)–(34), for = 1 , 2 , we have:
Θ n ( x ; ) M n ( x ; β , c ) = B 2 , n ( x ) M n j , ( x ; β , c ; λ ) B 1 , n ( x ) M n 1 j , ( x ; β , c ; λ ) ,
and
Θ n ( x ; ) M n 1 ( x ; β , c ) = A 1 , n ( x ) M n 1 j , ( x ; β , c ; λ ) A 2 , n ( x ) M n j , ( x ; β , c ; λ ) ,
where
Θ n ( x ; ) = det A 1 , n ( x ) B 1 , n ( x ) A 2 , n ( x ) B 2 , n ( x ) , = 1 , 2 .
After replacing the above in (35)–(38), we conclude:
Θ ˜ n ( x ; ) T + Λ 2 , n ( 1 ) ( x ; ) M n j , ( x ; β , c ; λ ) = Λ 1 , n ( 1 ) ( x ; ) M n 1 j , ( x ; β , c ; λ ) .
and
Θ ˜ n ( x ; ) T + Λ 1 , n ( 2 ) ( x ; ) M n 1 j , ( x ; β , c ; λ ) = Λ 2 , n ( 2 ) ( x ; ) M n j , ( x ; β , c ; λ ) ,
respectively, where
Θ ˜ n ( x ; ) = x Θ n ( x ; ) , = 1 , ( x + β ) Θ n ( x ; ) , = 2 ,
and
Λ j , n ( k ) ( x ; ) = ( 1 ) k det C k , n ( x ) A ν , n j , ( x ) D k , n ( x ) B ν , n j , ( x ) , ν = 1 , 2 , k = 1 , 2 , = 1 , 2 .
Proposition 5.
Let M n j , ( x ; β , c ; λ ) n N 0 be the sequence of monic Sobolev-type Meixner polynomials defined by (31) and let I be the identity operator. Then, the ladder (destruction and creation) operators a , a are defined by:
a = Θ ˜ n ( x ; ) T + Λ 2 , n ( 1 ) ( x ; ) ,
a = Θ ˜ n ( x ; ) T + Λ 1 , n ( 2 ) ( x ; ) ,
which verify
a ( M n j , ( x ; β , c ; λ ) ) = Λ 1 , n ( 1 ) ( x ; ) M n 1 j , ( x ; β , c ; λ ) ,
a ( M n 1 j , ( x ; β , c ; λ ) ) = Λ 2 , n ( 2 ) ( x ; ) M n j , ( x ; β , c ; λ ) ,
where Θ ˜ n ( x ; ) and Λ j , n ( k ) ( x ; ) with j , k , = 1 , 2 are given in (39) and (40).
Theorem 2.
The monic Sobolev-type Meixner polynomial sequence, which is orthogonal with respect to the inner product (24), fulfills the second-order difference equation:
F n ( x ; ) T 2 y ( x ) + G n ( x ; ) T y ( x ) + H n ( x ; ) y ( x ) = 0 ,
where
F n ( x ; ) = Θ ˜ n ( x ; ) Θ ˜ n ( x + ( 1 ) ; ) Λ 1 , n ( 1 ) ( x + ( 1 ) ; ) ,
G n ( x ; ) = Θ ˜ n ( x ; ) Λ 1 , n ( 1 ) ( x + ( 1 ) ; ) T Θ ˜ n ( x ; ) Θ ˜ n ( x ; ) T Λ 1 , n ( 1 ) ( x ; ) Λ 1 , n ( 1 ) ( x ; ) + Λ 2 , n ( 1 ) ( x + ( 1 ) ; ) ) + Θ ˜ n ( x ; ) Λ 1 , n ( 2 ) ( x ; ) Λ 1 , n ( 1 ) ( x ; ) ,
and
H n ( x ; ) = Θ ˜ n ( x ; ) T Λ 2 , n ( 1 ) ( x ; ) Λ 1 , n ( 1 ) ( x + ( 1 ) ; ) Θ ˜ n ( x ; ) Λ 2 , n ( 1 ) ( x ; ) T Λ 1 , n ( 1 ) ( x ; ) Λ 1 , n ( 1 ) ( x ; ) Λ 1 , n ( 1 ) ( x + ( 1 ) ; i ) + Λ 1 , n ( 2 ) ( x ; ) Λ 2 , n ( 1 ) ( x ; ) Λ 1 , n ( 1 ) ( x ; ) Λ 2 , n ( 2 ) ( x ; ) ,
where Θ ˜ n ( x ; ) and Λ j , n ( k ) ( x ; ) with j , k , = 1 , 2 are given in (39) and (40).
Proof. 
From (41), we have:
1 Λ 1 , n ( 1 ) ( x ; ) a ( M n j , ( x ; β , c ; λ ) ) = M n 1 j , ( x ; β , c ; λ ) .
Next, applying the operator a to both members of the previous expression, we obtain:
a 1 Λ 1 , n ( 1 ) ( x ; ) a ( M n j , ( x ; β , c ; λ ) ) = Λ 2 , n ( 2 ) ( x ; ) M n j , ( x ; β , c ; λ ) .
Thus, by using the definitions of the operators a and a , taking into account the identity:
T f ( x ) g ( x ) = g ( x ) T f ( x ) f ( x ) T g ( x ) g ( x ) g ( x + ( 1 ) ) ,
and after tedious calculations we obtain (42). Hence, the result follows. □

5. The ( 2 j + 3 ) -Term Recurrence Relation

In this section, we find the ( 2 j + 3 ) -term recurrence relation that the sequence of monic Sobolev-type Meixner polynomials (31) fulfill. For this purpose, we use the fact, which is a straightforward consequence of (3), that the multiplication operator by x α j + 1 is a symmetric operator with respect to such a discrete Sobolev inner product. Indeed, for any p , q P and for = 1 we have the following:
( x α ) j + 1 p ( x ) , q ( x ) λ , j , = ( x α ) j + 1 u M , p ( x ) q ( x ) = p ( x ) , ( x α ) j + 1 q ( x ) λ , j , ,
and for = 2 :
( x α j ) j + 1 p ( x ) , q ( x ) λ , j , = ( x α j ) j + 1 u M , p ( x ) q ( x ) = p ( x ) , ( x α j ) j + 1 q ( x ) λ , j , .
Taking these identities into account and by using the three-term recurrence relation (7), we can state the following result.
Lemma 1.
The following identities related to the monic Sobolev-type Meixner polynomials hold:
( x α ) j + 1 M n j , 1 ( x ; β , c ; λ ) = A n 1 ( x ) M n ( x ; β , c ) + B n 1 ( x ) M n 1 ( x ; β , c ) ,
( x α j ) j + 1 M n j , 2 ( x ; β , c ; λ ) = A n 2 ( x ) M n ( x ; β , c ) + B n 2 ( x ) M n 1 ( x ; β , c ) ,
where A n ( x ) , B n ( x ) are polynomials which can be computed explicitly.
Theorem 3.
Let λ R + and j N , and let M n j , ( x ; β , c ; λ ) n N 0 be the sequence of monic Sobolev-type Meixner polynomials defined by (31).
Then, the norm of these polynomials fulfills the following identity:
M n j , ( x ; β , c ; λ ) λ , j , 2 = M n ( x ; β , c ) 2 + b n j , M n 1 ( x ; β , c ) 2 ,
where
b n j , = λ M n 1 ( x ; β , c ) 2 T j M n ( α ; β , c ) T j M n ( α ; β , c ) 1 + λ K n , ( j , j ) ( α , α ) 0 .
Proof. 
We will consider the = 1 case. The = 2 case is analogous.
By the property of orthogonality of Sobolev-type Meixner polynomials, we have:
M n j , ( x ; β , c ; λ ) λ , j , 2 = M n j , ( x ; β , c ; λ ) , ( x α ) j + 1 π n j 1 ( x ) λ , j , ,
for any monic polynomial π of degree n j 1 . From (43), we have:
M n j , ( x ; β , c ; λ ) , ( x α ) j + 1 π ( x ) λ , j , = ( x α ) j + 1 M n j , ( x ; β , c ; λ ) , π ( x ) λ , j , .
By using the connection Formula (45) and taking into account that A n ( x ) is a monic polynomial of degree j + 1 and B n ( x ) is a polynomial of degree j with the leading coefficient b n j , , we deduce:
M n j , ( x ; β , c ; λ ) λ , j , 2 = ( x α ) j + 1 M n j , ( x ; β , c ; λ ) , π ( x ) λ , j , = u M , A n ( x ) M n ( x ; β , c ) π ( x ) + u M , B n ( x ) M n 1 ( x ; β , c ) π ( x ) = u M , M n ( x ; β , c ) x n + b n j , u M , M n 1 ( x ; β , c ) x n 1 ,
which coincides with (47). □
Remark 3.
Observe that a direct consequence is the following:
M n j , ( x ; β , c ; λ ) λ , j , 2 M n ( x ; β , c ) 2 = 1 + λ K n + 1 , ( j , j ) ( α , α ) 1 + λ K n , ( j , j ) ( α , α ) .
Theorem 4
( ( 2 j + 3 ) -term recurrence relation). Let λ R + , j N 0 . Then, the monic Sobolev-type Meixner orthogonal polynomials sequence with respect to the inner product (24) satisfies the following ( 2 j + 3 ) -term recurrence relation:
( x α ) j + 1 M n j , 1 ( x ; β , c ; λ ) = M n + j + 1 j , 1 ( x ; β , c ; λ ) + k = n j 1 n + j c n , k j , 1 M k j , 1 ( x ; β , c ; λ ) ,
and
( x α j ) j + 1 M n j , 2 ( x ; β , c ; λ ) = M n + j + 1 j , 2 ( x ; β , c ; λ ) + k = n j 1 n + j c n , k j , 2 M k j , 2 ( x ; β , c ; λ ) ,
where the constant coefficients c n , k j , can be explicitly computed for = 1 , 2 .
Proof. 
In such a case, we will consider the = 2 case.
Since the Sobolev-type Meixner polynomials form a basis in L 2 ( · , · λ , j , ) , if we consider the Fourier expansion of ( x α j ) j + 1 M n j , ( x ; β , c ; λ ) in terms of the Sobolev-type Meixner polynomials, then:
( x α j ) j + 1 M n j , ( x ; β , c ; λ ) = M n + j + 1 j , ( x ; β , c ; λ ) + k = 0 n + j c n , k j , M k j , ( x ; β , c ; λ ) .
Thus, by using the property of orthogonality of the sequence { M n j , ( x ; β , c ; λ ) } , we obtain:
c n , k j , = ( x α j ) j + 1 M n j , ( x ; β , c ; λ ) , M k j , ( x ; β , c ; λ ) λ , j , M k j , ( x ; β , c ; λ ) λ , j , 2 , k = 0 , , n + j .
Using (44) and the property of orthogonality of { M n j , ( x ; β , c ; λ ) } , we deduce that c n , k j , = 0 for k = 0 , , n j 2 . Observe that the rest of the coefficients can be computed by using again the same orthogonality conditions. The proof of the another identity is similar, and it will be omitted. □

6. Mehler–Heine-Type Formula

The main result of this section will be to establish a Mehler–Heine-type formula for the polynomial M n j , ( x ; β , c ; λ ) for the α 0 case. The following results will allow us to prove such an expression.
Lemma 2.
Let β , c C , with | c | < 1 and β N , and let m be a positive integer. Then, the following limit holds:
lim n ( β ) n n m c n ( n 1 ) ! = 0 .
Proof. 
If we use the identity (p. 23) [21]
Γ ( z ) = lim n ( n 1 ) ! n z ( z ) n ,
we deduce:
lim n ( β ) n n m c n ( n 1 ) ! = 1 Γ ( β ) lim n n β + m c n .
Therefore, if Re ( β + m ) > 0 , then, applying to (51) the L’Hôpital’s rule several times, we obtain the desired result; otherwise, the limit is zero. Hence, the result holds. □
Lemma 3.
Let β , c R , with | c | < 1 , and β is not a negative integer, and let k , j be integers, with 0 k j . If we set α = 0 in (25), then the following limits hold:
lim n A 1 , n j , ( x ) = 1 a n d lim n B 1 , n j , ( x ) = 0 , = 1 , 2 .
Proof. 
By starting with (25) and using (21), we obtain:
A i , n j , ( x ) = 1 λ T j M n ( α ; β , c ) 1 + λ K n , ( j , j ) ( α , α ) j ! d n 1 2 k = 0 j T k M n 1 ( α ; β , c ) k ! ( x α + k ) j + 1 k ,
and
B i , n j , ( x ) = λ T j M n ( α ; β , c ) 1 + λ K n , ( j , j ) ( α , α ) j ! d n 1 2 k = 0 j T k M n ( α ; β , c ) k ! ( x α + k ) j + 1 k .
where i = 1 , 2 . Then, to prove this result, it is enough to check:
lim n ( 1 c ) n a k , n j , = lim n ( 1 c ) n b k , n j , = 0 , k = 0 , 1 , , n , i = 1 , 2 , = 1 , 2 ,
where
a k , n j , = λ j ! T j M n ( α ; β , c ) T k M n 1 ( α ; β , c ) k ! ( 1 + λ K n , ( j , j ) ( α , α ) ) d n 1 2 ,
and
b k , n j , = λ j ! T j M n ( α ; β , c ) T k M n ( α ; β , c ) k ! ( 1 + λ K n , ( j , j ) ( α , α ) ) d n 1 2 .
After a straightforward calculation, by using (15), we have that for any 0 k j
T k M n ( α ; β , c ) T j M n ( α ; β , c ) ,
for n large, and since c 1 < c , then by using Lemma 2, it is clear that both limits related to such coefficients tend to zero. Hence, we deduce (52). □
Theorem 5.
Let β , c R , with | c | < 1 and β is not a negative integer, and let m be a positive integer. Then, we have:
lim n ( 1 c ) n + β + z M n j , ( z ; β , c ; λ ) ( z n + 1 ) n = 1 , z C N ,
uniformly on compact subsets of the complex plane.
Proof. 
If we multiply (26) by ( 1 c ) n + β + z / ( z n + 1 ) n , we obtain:
( 1 c ) n + β + z M n j , ( z ; β , c ; λ ) ( z n + 1 ) n = A i , n ( j , ) ( z ) ( 1 c ) n + β + z M n ( z ; β , c ) ( z n + 1 ) n + B 1 , n ( j , ) ( z ) ( 1 c ) n + β + z M n 1 ( z ; β , c ) ( z n + 1 ) n .
Then, applying the previous Lemma as well as the (15), we arrive at the desired result. □
Remark 4.
Observe that we can extend some of the previous results even for c C , so that | c 1 | < | c | < 1 , or even into a wider region of the complex plane taking into account (6).
Finally, we show some graphical experiments of the limit function in (53) for several values of n using Mathematica software at the masspoint α = 0 . Since
( z n + 1 ) n = ( 1 ) n Γ ( z + n ) Γ ( z ) ,
then (53) can be rewritten as:
lim n ( c 1 ) n M n j , ( z ; β , c ; λ ) Γ ( n z ) = 1 ( 1 c ) β + z Γ ( z ) , z C N .
The Figure 1, Figure 2 and Figure 3 show how such limit converges to the right value.

Author Contributions

R.S.C.-S., A.S.-L. and J.-M.V. conceived the mathematics; R.S.C.-S., A.S.-L. and J.-M.V. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

The research of RSCS was funded by Agencia Estatal de Investigación of Spain, grant number PGC-2018-096504-B-C33. The work of ASL is supported by Dirección General de Investigación e Innovación, Consejería de Educación e Investigación of the Comunidad de Madrid (Spain) and Universidad de Alcalá under grant CM/JIN/2021-014, Proyectos de I+D para Jóvenes Investigadores de la Universidad de Alcalá 2021.

Acknowledgments

We are grateful for the exhaustive work of the referees. Their comments and suggestions have improved the presentation of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Meixner, J. Orthogonale polynomsysteme mit einem besonderen Gestalt der erzeugenden funktion. J. Lond. Math. Soc. 1934, 9, 6–13. [Google Scholar] [CrossRef]
  2. Chihara, T.S. An Introduction to Orthogonal Polynomials; Mathematics and its Applications Series; Gordon and Breach: New York, NY, USA, 1978. [Google Scholar]
  3. Ismail, M.E.H. Classical and Quantum Orthogonal Polynomials in One Variable; Encyclopedia of Mathematics and Its Applications; Cambridge University Press: Cambridge, CA, USA, 2005; Volume 98. [Google Scholar]
  4. Nikiforov, A.F.; Uvarov, V.B.; Suslov, S.K. Classical Orthogonal Polynomials of a Discrete Variable; Springer Series in Computational Physics; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  5. Althammer, P. Eine Erweiterung des Orthogonalitätsbegriffes bei Polynomen und deren Anwendung auf die beste Approximation. J. Reine Angew. Math. 1962, 211, 192–204. [Google Scholar] [CrossRef]
  6. Marcellán, F.; Xu, Y. On Sobolev orthogonal polynomials. Expo. Math. 2015, 33, 308–352. [Google Scholar] [CrossRef] [Green Version]
  7. Bavinck, H. On polynomials orthogonal with respect to an inner product involving differences. J. Comp. Appl. Math. 1995, 57, 17–27. [Google Scholar] [CrossRef] [Green Version]
  8. Bavinck, H. On polynomials orthogonal with respect to an inner product involving differences (the general case). Appl. Anal. 1995, 59, 233–240. [Google Scholar] [CrossRef]
  9. Bavinck, H.; Haeringen, H.V. Difference equations for generalized Meixner polynomials. J. Math. Anal. Appl. 1994, 184, 453–463. [Google Scholar] [CrossRef] [Green Version]
  10. Bavinck, H.; Koekoek, R. Difference operators with sobolev type meixner polynomials as eigenfunctions. Comput. Math. Appl. 1998, 36, 163–177. [Google Scholar] [CrossRef] [Green Version]
  11. Khwaja, S.F.; Olde-Daalhui, A.B. Uniform asymptotic approximations for the Meixner–Sobolev polynomials. Anal. Appl. 2012, 10, 345–361. [Google Scholar] [CrossRef]
  12. Moreno-Balcázar, J. δ-Meixner-Sobolev orthogonal polynomials: Mehler–heine type formula and zeros. J. Comput. Appl. Math. 2015, 284, 228–234. [Google Scholar] [CrossRef]
  13. Huertas, E.J.; Soria, A. On Second Order q-Difference Equations Satisfied by Al-Salam–Carlitz I-Sobolev Type Polynomials of Higher Order. Mathematics 2020, 8, 1300. [Google Scholar]
  14. Nodarse, R.x.; Marcellán, F. Difference equation for modifications of Meixner polynomials. J. Math. Anal. Appl. 1995, 194, 250–258. [Google Scholar]
  15. Huertas, E.J.; Marcellán, F.; Pérez-Valero, M.F.; Quintana, Y. Asymptotics for Laguerre–Sobolev type orthogonal polynomials modified within their oscillatory regime. Appl. Math. Comput. 2014, 236, 260–272. [Google Scholar] [CrossRef] [Green Version]
  16. Koekoek, R.; Lesky, P.A.; Swarttouw, R.F. Hypergeometric Orthogonal Polynomials and Their q-Analogues; Springer Monographs in Mathematics; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  17. Costas, R.S.; Sanchez-Lara, J.F. Extensions of discrete classical orthogonal polynomials beyond the orthogonality. J. Comput. Appl. Math. 2009, 225, 440–451. [Google Scholar] [CrossRef] [Green Version]
  18. Álvarez de Morales, M.; Pérez, T.E.; Piñar, M.A.; Ronveaux, A. Non-standard orthogonality for Meixner polynomials. Orthogonal polynomials: Numerical and symbolic algorithms (Leganés, 1998). Electron. Trans. Numer. Anal. 1999, 9, 1–25. [Google Scholar]
  19. Dominici, D. Mehler–Heine type formulas for Charlier and Meixner polynomials. Ramanujan J. 2016, 39, 271–289. [Google Scholar] [CrossRef] [Green Version]
  20. Marcellán, F.; Ronveaux, A. On a class of polynomials orthogonal with respect to a discrete Sobolev inner product. Indag. Math. 1990, 4, 451–464. [Google Scholar] [CrossRef] [Green Version]
  21. Rainville, D. Special Functions; The Macmillan Company: New York, NY, USA, 1960. [Google Scholar]
Figure 1. The gray-color curve represents the function in the left hand side in (54) for n = 50 , and the black curve is the limiting function. Data: β = 7 , c = 1 / 5 , λ = 10 21 and j = 2 .
Figure 1. The gray-color curve represents the function in the left hand side in (54) for n = 50 , and the black curve is the limiting function. Data: β = 7 , c = 1 / 5 , λ = 10 21 and j = 2 .
Mathematics 10 01952 g001
Figure 2. The gray-color curve represents the function in the left hand side in (54) for n = 70 , and the black curve is the limiting function. Data: β = 7 , c = 1 / 5 , λ = 10 21 and j = 2 .
Figure 2. The gray-color curve represents the function in the left hand side in (54) for n = 70 , and the black curve is the limiting function. Data: β = 7 , c = 1 / 5 , λ = 10 21 and j = 2 .
Mathematics 10 01952 g002
Figure 3. The gray-color curve represents the function in the left hand side in (54) for n = 100 , and the black curve is the limiting function. Data: β = 7 , c = 1 / 5 , λ = 10 21 and j = 2 .
Figure 3. The gray-color curve represents the function in the left hand side in (54) for n = 100 , and the black curve is the limiting function. Data: β = 7 , c = 1 / 5 , λ = 10 21 and j = 2 .
Mathematics 10 01952 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Costas-Santos, R.S.; Soria-Lorente, A.; Vilaire, J.-M. On Polynomials Orthogonal with Respect to an Inner Product Involving Higher-Order Differences: The Meixner Case. Mathematics 2022, 10, 1952. https://0-doi-org.brum.beds.ac.uk/10.3390/math10111952

AMA Style

Costas-Santos RS, Soria-Lorente A, Vilaire J-M. On Polynomials Orthogonal with Respect to an Inner Product Involving Higher-Order Differences: The Meixner Case. Mathematics. 2022; 10(11):1952. https://0-doi-org.brum.beds.ac.uk/10.3390/math10111952

Chicago/Turabian Style

Costas-Santos, Roberto S., Anier Soria-Lorente, and Jean-Marie Vilaire. 2022. "On Polynomials Orthogonal with Respect to an Inner Product Involving Higher-Order Differences: The Meixner Case" Mathematics 10, no. 11: 1952. https://0-doi-org.brum.beds.ac.uk/10.3390/math10111952

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop