Next Article in Journal
Position and Attitude Tracking Finite-Time Adaptive Control for a VTOL Aircraft Using Global Fast Terminal Sliding Mode Control
Next Article in Special Issue
Differential Properties of Jacobi-Sobolev Polynomials and Electrostatic Interpretation
Previous Article in Journal
Pricing Variance Swaps under MRG Model with Regime-Switching: Discrete Observations Case
Previous Article in Special Issue
Sequentially Ordered Sobolev Inner Product and Laguerre–Sobolev Polynomials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Look at Generalized Degenerate Bernoulli and Euler Matrices

by
Juan Hernández
1,*,
Dionisio Peralta
1 and
Yamilet Quintana
2,3
1
Escuela de Matemáticas, Facultad de Ciencias, Universidad Autónoma de Santo Domingo, Santo Domingo 10105, Dominican Republic
2
Departamento de Matemáticas, Universidad Carlos III de Madrid, Avenida de la Universidad 30, Leganés, 28911 Madrid, Spain
3
Instituto de Ciencias Matemáticas (ICMAT), Campus de Cantoblanco UAM, 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2731; https://doi.org/10.3390/math11122731
Submission received: 24 May 2023 / Revised: 14 June 2023 / Accepted: 14 June 2023 / Published: 16 June 2023

Abstract

:
In this paper, we consider the generalized degenerate Bernoulli/Euler polynomial matrices and study some algebraic properties for them. In particular, we focus our attention on some matrix-inversion formulae involving these matrices. Furthermore, we provide analytic properties for the so-called generalized degenerate Pascal matrix of the first kind, and some factorizations for the generalized degenerate Euler polynomial matrix.

1. Introduction

Matrices play an important role in all branches of science, engineering, social science, and management. In many settings (see, e.g., [1,2,3,4] and the references therein), a number of interesting and useful identities involving binomial (q-binomial or λ -binomial) coefficients can be obtained from a matrix representation of a particular counting sequence. Such a matrix representation provides a powerful computational tool for deriving identities and an explicit formula related to the sequence.
There are many special types of matrices such as Pascal, Vandermonde, Stirling, Riordan arrays, and others. These matrices are of specific importance in many scientific and engineering applications. For instance, Pascal matrices appear in combinatorics, image processing, signal processing, numerical analysis, probability, and surface reconstruction.
In the case of generalized Pascal matrices of the first kind, extensive research has been devoted to them (cf., e.g., [3,4,5,6,7,8,9,10] and the references therein). Situations with a matrix representation—including analogs of generalized Pascal matrices of the first kind and degenerate versions of special classes of polynomials (e.g., Bernstein, Bernoulli, and Euler polynomials, etc.)—are of particular interest.
Motivated by recent articles [1,2,3,4,11,12,13,14] that consider degenerate Bernstein polynomials, degenerate Euler polynomials, generalized degenerate Euler–Genocchi polynomials of order α , and algebraic properties of the generalized Euler and generalized Apostol-type polynomial matrices, in the present article, we consider the generalized degenerate Bernoulli/Euler polynomial matrix. In particular, we focus our attention on some inversion-type formulae from a matrix framework. Furthermore, we show some analytic properties for the so-called generalized degenerate Pascal matrix of the first kind. Furthermore, some factorizations for the generalized degenerate Euler polynomial matrix in terms of such a matrix are given.
The paper is organized as follows. Section 2 is a preliminary section containing the definitions, notations, and terminology needed. Section 3 contains the main results of this paper. First, we provide the corresponding inversion-type formulae for the degenerate Bernoulli and Euler polynomials, respectively (Theorems 1 and 2). Second, we show that the generalized degenerate Pascal matrix of the first kind is a matrix exponential (Theorem 4), and, as a consequence, we obtain an Appell-type property for this matrix (Corollary 5). In addition, factorizations for the generalized degenerate Pascal matrix of the first kind in terms of the degenerate Bernoulli/Euler matrices are given (Theorems 6 and 7, respectively). The remainder of this section is devoted to establishing the corresponding product formulae for generalized degenerate Euler polynomial matrices and their factorizations in terms of generalized degenerate Pascal matrices of the first kind (Theorems 8 and 9).

2. Background and Previous Results

Throughout this paper, let N , N 0 , Z , R , and C denote, respectively, the set of all natural numbers, the set of all non-negative integers, the set of all integers, the set of all real numbers, and the set of all complex numbers. As usual, we will always use the principal branch for complex powers, in particular, 1 α = 1 for α C . Furthermore, the convention 0 0 = 1 will be adopted.
For w C and k Z , we use the notations w ( k ) and ( w ) k for the rising and falling factorials, respectively, i.e.,
w ( k ) = 1 , if   k = 0 , i = 1 k ( w + i 1 ) , if   k 1 , 0 , if   k < 0 ,
and
( w ) k = 1 , if   k = 0 , i = 1 k ( w i + 1 ) , if   k 1 , 0 , if   k < 0 .
Any matrix is assumed an element of M n + 1 ( R ) , the set of all ( n + 1 ) -square matrices over the real field R . Moreover, for i , j , any nonnegative integers, and any matrix A M n + 1 ( R ) we adopt, respectively, the following conventions
i j = 0 , w h e n e v e r   j > i , a n d A 0 = I n + 1 = diag ( 1 , 1 , , 1 ) ,
where I n + 1 denotes the identity matrix of order n + 1 .
For λ , x R and z C , the degenerate exponentials are defined as follows (cf., [15]):
e λ x ( z ) = ( 1 + λ z ) x λ , if   λ R { 0 } , e x z , if   λ = 0 .
As usual, for x = 1 , we use the notation e λ ( z ) = e λ x ( z ) .
It follows immediately from (1) that
e λ x ( z ) = n = 0 ( x ) n , λ z n n ! , | λ z | < 1 , if   λ R { 0 } , n = 0 x n z n n ! , if   λ = 0 .
where the generalized falling factorials ( x ) n , λ , are given by (cf., [1,2,12,13,14,15]):
( x ) n , λ = 1 , if   n = 0 , i = 1 n ( x ( i 1 ) λ ) , if   n 1 , 0 , if   n < 0 ,
where x , λ R and n Z .
It is clear that lim λ 0 e λ x ( z ) = e 0 x ( z ) = e x z , and for n N 0 , the polynomial in two variables Q n ( x , λ ) , given by
Q n ( x , λ ) = 1 , if   n = 0 , i = 1 n ( x ( i 1 ) λ ) , if   n 1 ,
is a continuous function on R 2 , and consequently, ( x ) n , 0 = x n .
In [16,17], Carlitz introduced the degenerate Bernoulli (Euler) and the generalized degenerate Bernoulli (Euler) polynomials of order α C , respectively, by means of the generating functions and series expansions:
z e λ ( z ) 1 e λ x ( z ) = n = 0 B n , λ ( x ) z n n ! ,
2 e λ ( z ) + 1 e λ x ( z ) = n = 0 E n , λ ( x ) z n n ! ,
z e λ ( z ) 1 α e λ x ( z ) = n = 0 B n , λ ( α ) ( x ) z n n ! ,
2 e λ ( z ) + 1 α e λ x ( z ) = n = 0 E n , λ ( α ) ( x ) z n n ! .
These are valid in a suitable neighborhood of z = 0 and represent degenerate versions of the classical Bernoulli and Euler polynomials, respectively. In [8], the notation β n ( λ , x ) is used for the degenerate Bernoulli (3).
Since the degenerate exponentials (1) satisfy the same exponent product law as the exponentials functions, i.e.,
e λ x + y ( z ) = e λ x ( z ) e λ y ( z ) ,
we can use the generating relations (2), (5) and (6) to deduce the following addition formulas:
( x + y ) n , λ = k = 0 n n k ( x ) k , λ ( y ) n k , λ , n 0 ,
B n , λ ( α + β ) ( x + y ) = k = 0 n n k B k , λ ( α ) ( x ) B n k , λ ( β ) ( y ) , n 0 ,
E n , λ ( α + β ) ( x + y ) = k = 0 n n k E k , λ ( α ) ( x ) E n k , λ ( β ) ( y ) , n 0 .
For a treatment of diverse aspects of some summation formulas and their applications, the interested reader is referred to the relatively recent works [18,19,20].
For r N 0 , λ R , and α C , definitions of generalized degenerate Euler–Genocchi and generalized degenerate Euler–Genocchi polynomials of order α , respectively, have recently been introduced in [14] (Section 2):
2 z r e λ ( z ) + 1 e λ x ( z ) = n = 0 A n , λ ( r ) ( x ) z n n ! ,
z r 2 e λ ( z ) + 1 α e λ x ( z ) = n = 0 A n , λ ( r , α ) ( x ) z n n ! .
Remark 1.
Notice that:
 (i) 
If r N , then it follows immediately from (2), (4) and (10), that
A 0 , λ ( r ) ( x ) = A 1 , λ ( r ) ( x ) = = A r 1 , λ ( r ) ( x ) = 0 , and
A n , λ ( r ) ( x ) = n ! ( n r ) ! ( x ) n , λ = n ( r ) E n r , λ ( 0 ) ( x ) , n r .
Furthermore, A n , λ ( 0 ) ( x ) = E n , λ ( x ) , n 0 .
The first above identities guarantee that, up to multiplicative constants, it suffices to take generalized degenerate Euler polynomials of order 0 instead of the so-called generalized degenerate Euler–Genocchi polynomials as the main family to study. Similarly, the last identity tells us that the generalized degenerate Euler polynomials coincides with the generalized degenerate Euler–Genocchi polynomials of order 0.
 (ii) 
In [14], Theorem 4 proves the following reduction formula:
A n , λ ( r , α ) ( x ) = n ( r ) E n r , λ ( α ) ( x ) , n r , n , r N 0 .
In particular, we obtain that up to multiplicative constants, the generalized degenerate Euler–Genocchi polynomials of order α = 1 can be reduced to the generalized degenerate Euler polynomials (4).
Hence, in order to avoid essentially redundant definitions (cf., [21]), the families of polynomials eqrefeul-gen1 and (11) will not be considered in this paper.

3. The Generalized Degenerate Bernoulli and Euler Matrices and Their Properties

In this section, we present some novel properties for the generalized degenerate Bernoulli and Euler matrices. Before that, we show the corresponding inversion-type formulae for the generalized degenerate Bernoulli and Euler polynomials, respectively.
Theorem 1.
For every n 0 and λ R , the degenerate Bernoulli polynomials satisfy the following inversion-type formula:
( x ) n , λ = 1 n + 1 k = 0 n n + 1 k + 1 ( 1 ) k + 1 , λ B n k , λ ( x )
= 1 n + 1 k = 0 n n + 1 k + 1 ( 1 λ ) k , λ B n , λ ( x ) .
Proof. 
Let λ R . In view of (2) and (3), and the identity
z n = 0 ( x ) n , λ z n n ! = n = 0 ( n + 1 ) ( x ) n , λ z n + 1 ( n + 1 ) ! ,
we have
n = 0 ( n + 1 ) ( x ) n , λ z n + 1 ( n + 1 ) ! = n = 0 ( 1 ) n , λ z n n ! 1 n = 0 B n , λ ( x ) z n n ! = n = 0 ( 1 ) n + 1 , λ z n + 1 ( n + 1 ) ! n = 0 B n , λ ( x ) z n n ! .
From the use of the Cauchy product rule on the right-hand side of (14), it follows that
n = 0 ( n + 1 ) ( x ) n , λ z n + 1 ( n + 1 ) ! = n = 0 k = 0 n n + 1 k + 1 ( 1 ) k + 1 , λ B n k , λ ( x ) z n + 1 ( n + 1 ) ! .
Hence, comparing the coefficients of z n + 1 on both sides of (15), we obtain (12).
Finally, (13) is a simple consequence of the identity ( 1 ) k + 1 , λ = ( 1 λ ) k , λ , for all k N 0 . □
Remark 2.
Notice that the substitution of λ = 0 into (12) recovers the inversion formula for the classical Bernoulli polynomials (cf., [22] (Equation (9))).
From a matrix framework, Theorem 1 has the following consequence.
Corollary 1.
For n N 0 and λ R , the matrix T λ ( x ) = 1 ( x ) 1 , λ ( x ) n , λ T can be expressed as follows:
T λ ( x ) = M λ B λ ( x ) = 1 1 ( 1 ) 1 , λ 0 0 0 1 2 2 2 ( 1 ) 2 , λ 1 2 2 1 ( 1 ) 1 , λ 0 0 1 3 3 3 ( 1 ) 3 , λ 1 3 3 2 ( 1 ) 2 , λ 1 3 3 1 ( 1 ) 1 , λ 0 1 n + 1 n + 1 n + 1 ( 1 ) n + 1 , λ 1 n + 1 n + 1 n ( 1 ) n , λ 1 n + 1 n + 1 n 1 ( 1 ) n 1 , λ 1 n + 1 n + 1 1 ( 1 ) 1 , λ B λ ( x ) = 1 0 0 0 1 2 ( 1 ) 2 , λ 1 0 0 1 3 ( 1 ) 3 , λ ( 1 ) 2 , λ 1 0 1 n + 1 ( 1 ) n + 1 , λ ( 1 ) n , λ 1 2 ( 1 ) n 1 , λ 1 B λ ( x ) ,
where B λ ( x ) = B 0 , λ ( x ) B 1 , λ ( x ) B n , λ ( x ) T .
Theorem 2.
For every n 0 and λ R . The degenerate Euler polynomials satisfy the following inversion-type formula:
( x ) n , λ = 1 2 k = 0 n n k ( 1 + a k ( λ ) ) ( 1 ) k , λ E n k , λ ( x )
where
a k ( λ ) = 1 , i f   k = 0 , 0 , i f   1 k n ,
Proof. 
From (2) and (4) we have
2 n = 0 ( x ) n , λ z n n ! = n = 0 ( 1 ) n , λ z n n ! + 1 n = 0 E n , λ ( x ) z n n ! = n = 0 ( 1 + a k ( λ ) ) ( 1 ) n , λ z n n ! n = 0 E n , λ ( x ) z n n ! = n = 0 k = 0 n ( 1 + a k ( λ ) ) n k ( 1 ) k , λ E n k , λ ( x ) z n n ! ,
where
a k ( λ ) = 1 , if   k = 0 , 0 , if   1 k n .
Therefore, by comparing the coefficients of z n on both sides, we obtain the identity. □
Remark 3.
Notice that if λ = 0 in (17), then we recover the inversion formula for the classical Euler polynomials (cf., [22] (Equation (27))).
Theorem 2 has the following consequence.
Corollary 2.
For n N 0 and λ R , the matrix T λ ( x ) = 1 ( x ) 1 , λ ( x ) n , λ T can be expressed as follows:
T λ ( x ) = 1 2 N λ E λ ( x ) = 1 2 0 0 ( 1 + a 0 ( λ ) ) ( 1 ) 0 , λ 0 0 1 1 ( 1 + a 1 ( λ ) ) ( 1 ) 1 , λ 1 0 ( 1 + a 0 ( λ ) ) ( 1 ) 0 , λ 0 2 2 ( 1 + a 2 ( λ ) ) ( 1 ) 2 , λ 2 1 ( 1 + a 1 ( λ ) ) ( 1 ) 1 , λ 0 n n ( 1 + a n ( λ ) ) ( 1 ) n , λ n n 1 ( 1 + a n 1 ( λ ) ) ( 1 ) n 1 , λ n 0 ( 1 + a 0 ( λ ) ) ( 1 ) 0 , λ E λ ( x ) = 1 2 2 0 0 0 0 ( 1 ) 1 , λ 2 0 0 0 ( 1 ) 2 , λ 2 ( 1 ) 1 , λ 2 0 0 ( 1 ) n , λ n ( 1 ) n 1 , λ ( n ) 2 2 ! ( 1 ) n 2 , λ ( n ) 3 3 ! ( 1 ) n 3 , λ 2 E λ ( x ) ,
where E λ ( x ) = E 0 , λ ( x ) E 1 , λ ( x ) E n , λ ( x ) T and a k ( λ ) = 1 , i f   k = 0 , 0 , i f   1 k n .
Clearly, when λ R , the matrix N λ is an invertible matrix.
Corollary 3.
For n N 0 and λ R , we have
E λ ( x ) = 2 N λ 1 M λ B λ ( x ) .
The degenerate Pascal matrices corresponding to the generalized falling factorials can be defined as follows:
Definition 1.
Let x be any nonzero real number. For λ R , the generalized degenerate Pascal matrix of the first kind P λ [ x ] , is an ( n + 1 ) × ( n + 1 ) matrix whose entries are given by
p i , j , λ ( x ) : = i j ( x ) i j , λ , i j , 0 , o t h e r w i s e .
Remark 4.
 (i) 
It is clear that the matrix P λ [ x ] tends to the generalized Pascal matrix of the first kind P [ x ] as λ 0 .
 (ii) 
For n N 0 , x R { 0 } , λ R , it is clear that P λ [ x ] = P n , λ [ x ] , where P n , λ [ x ] is the Pascal functional matrix introduced in [5]. Hence, all results corresponding to P λ [ x ] given in [5] hold in this setting.
 (iii) 
It is worth mentioning that the matrix entries (19) coincide with the entries of the variation of Pascal functional matrix P n [ x , λ ] introduced by Can and Cihat-Dağli in [8]. Hence, all results corresponding to factorizing the matrix P n [ x , λ ] by the summation matrices also hold for P λ [ x ] , taking into account the suitable shift on the respective order for these matrices (cf., [8] (Lemma 1 and Theorem 2)).
 (iv) 
If for x R { 0 } , λ R we consider the truncated exponential generating function for the binomial-type polynomial sequence { ( x ) n , λ } n 0 (cf., [9]):
f ( t ; x ) = k = 0 n ( x ) k , λ t k k ! ,
then, it is easy to see that
P λ [ x ] = P n [ f ( x , t ) ] t = 0 = P n k = 0 n ( x ) k , λ t k k ! t = 0 ,
where P n [ f ( t ; x ) ] denotes the generalized Pascal functional matrix introduced by Yang and Micek in [9].
From now on, we denote P λ = P λ [ 1 ] . The following theorem summarizes some properties of P λ [ x ] .
Theorem 3.
Let P λ [ x ] M n + 1 ( R ) be the generalized degenerate Pascal matrix of the first kind. Then, the following statements hold.
 (a) 
Special value. If the convention ( 0 ) 0 , λ = 1 is adopted, then it is possible to define
P λ [ 0 ] : = I n + 1 .
 (b) 
For x , y R , we have
P λ [ x + y ] = P λ [ x ] P λ [ y ] .
 (c) 
P λ [ x ] is an invertible matrix and its inverse is given by
P λ 1 [ x ] : = P λ [ x ] 1 = P λ [ x ] .
Proof. 
Since part (a) is a straightforward consequence of the extension of Definition 1 for the case x = 0 , we shall omit its proof. Thus, we focus our efforts on the proof of parts (b) and (c).
Let A i , j , λ ( x , y ) be the ( i , j ) -th entry of the matrix product P λ [ x ] P λ [ y ] . Then, by (7), we have
A i , j , λ ( x , y ) = k = 0 n i k ( x ) i k , λ k j ( y ) k j , λ = k = j i i k ( x ) i k , λ k j ( y ) k j , λ = k = j i i j i j i k ( x ) i k , λ ( y ) k j , λ = i j k = 0 i j i j k ( x ) i j k , λ ( y ) k , λ = i j ( x + y ) i j , λ ,
which implies (20).
The substitution y = x into (20) yields
P λ [ 0 ] = P λ [ x ] P λ [ x ] = P λ [ x ] P λ [ x ] .
By part (a), we have P λ [ 0 ] = I n + 1 , thus
P λ [ x ] P λ [ x ] = I n + 1 = P λ [ x ] P λ [ x ] ,
and (21) follows. □
Corollary 4.
For any λ R , r Z and s Z { 0 } we have
 (a) 
P λ r = P λ [ r ] .
 (b) 
P λ r s s = P λ r .
Proof. 
Making the corresponding modifications, we apply the same reasoning as in the proof of [7] (Corollary 3). Since P λ = P λ [ 1 ] , P λ [ 0 ] , and P λ 0 coincide with the identity matrix, it follows from Theorem 3, by induction on r, that P λ [ r ] = P λ r , for all r N 0 . Again, by Theorem 3, we have that P λ [ 1 ] = P λ 1 , and a similar induction on | r | shows P λ [ r ] = P λ r , for all r < 0 .
Next, by Theorem 3 and part (a), we obtain P λ r s s = P λ [ r ] = P λ r . □
Remark 5.
Part (b) of Corollary 4 shows that for a fixed λ R and any rational number x, P λ [ x ] is the x-th power of P λ . Indeed, this property could be expected in the sense that it is satisfied for the generalized Pascal matrix of the first kind P [ x ] (cf., [7]).
From the addition Formula (20), we proceed according to [7] and conclude that the degenerate Pascal matrix P λ [ x ] has an exponential form as follows: Assume that for λ R , there is a matrix L λ , such that P λ [ x ] = e x L λ . Then,
d d x P λ [ x ] = L λ e x L λ = L λ P λ [ x ] ,
and
d d x P λ [ x ] x = 0 = L λ P λ [ 0 ] = L λ I n + 1 = L λ .
Thus, there is at most one matrix L λ such that P λ [ x ] = e x L λ . For instance, in the case n = 3 , we can find the only possible value as follows:
d d x P λ [ x ] = 0 0 0 0 1 0 0 0 λ + 2 x 2 0 0 x ( 2 λ + x ) + x ( λ + x ) + ( 2 λ + x ) ( λ + x ) 3 ( λ + 2 x ) 3 0 ,
and
L λ = d d x P λ [ x ] x = 0 = 0 0 0 0 1 0 0 0 λ 2 0 0 2 λ 2 3 λ 3 0 .
While, in the case n = 7 , we have
L λ = d d x P λ [ x ] x = 0 = 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 λ 2 0 0 0 0 0 0 2 λ 2 3 λ 3 0 0 0 0 0 6 λ 3 8 λ 2 6 λ 4 0 0 0 0 24 λ 4 30 λ 3 20 λ 2 10 λ 5 0 0 0 120 λ 5 144 λ 4 90 λ 3 40 λ 2 15 λ 6 0 0 720 λ 6 840 λ 5 504 λ 4 210 λ 3 70 λ 2 21 λ 2 7 0 .
This suggests a general way of choosing L λ . More precisely, the entries of L λ are given by
L λ i , j = s λ ( i j , 1 ) i j , if   i j + 1 , 0 , o t h e r w i s e ,
where s λ ( n , k ) denotes the degenerate Stirling number of the first kind, defined as follows (cf., [17,23] or [24] (Ch. 5)):
k = 0 n s λ ( n , k ) x k = ( x ) n , λ .
Furthermore, the entries of the matrix L λ k , for 1 k n and n N can be explicitly represented as follows.
Lemma 1.
For every n N and 1 k n , the entries of L λ k are given by the formula
L λ k i , j = k ! s λ ( i j , k ) i j , i f   i j + k , 0 , o t h e r w i s e ,
where s λ ( n , k ) is the degenerate Stirling number of the first kind (22).
Proof. 
It suffices to proceed by induction on k, taking into account that for k > n , we have L λ k = 0 . □
Theorem 4.
For every real numbers x , λ R , P λ [ x ] = e x L λ .
Proof. 
By part (a) of Theorem 3, if x = 0 , then e x L λ = I n + 1 = P λ [ x ] . Now, assume that x 0 since L λ k = 0 for k > n , the infinite series for e x L λ reduces to the finite sum
e x L λ = I n + 1 + x L λ + x 2 2 L λ 2 + + x n n ! L λ n .
Applying Lemma 1, we can now read off the entries in e x L λ . Clearly, it is a lower triangular matrix, and the diagonal entries are all 1. Now suppose i > j , and let 0 k i j . Then, using (22), we have that the ( i , j ) -th entry in the sum (23) is
e x L λ i , j = k = 0 i j x k k ! L λ k i , j = i j k = 0 i j s λ ( i j , k ) x k = i j ( x ) i j , λ = p i , j , λ ( x ) .
This completes the proof. □
As a consequence of Lemma 1 and Theorem 4, we obtain the following Appell-type property.
Corollary 5.
The generalized degenerate Pascal matrix of the first kind P λ [ x ] satisfies the following differential equations:
D x k P λ [ x ] = L λ k P λ [ x ] , 1 k n ,
where D x k P λ [ x ] is the matrix resulting from the k-th derivative with respect to x of each entry of P λ [ x ] .
Definition 2.
The generalized degenerate ( n + 1 ) × ( n + 1 ) Bernoulli matrix B λ ( α ) ( x ) of (real or complex) order α is defined by the entries
B i , j , λ ( α ) ( x ) = i j B i j , λ ( α ) ( x ) , i j , 0 , o t h e r w i s e .
Remark 6.
 (i) 
It is worth mentioning that the entries (2) of B λ ( α ) ( x ) coincide with the entries of the generalized degenerate Bernoulli matrix B m ( α ) [ λ , x ] introduced in [8], when these matrices are the same order.
 (ii) 
We denote by B λ ( x ) the degenerate Bernoulli matrix B λ ( 1 ) ( x ) .
The following result was established in [8] (Theorem 4).
Theorem 5.
The generalized degenerate Bernoulli matrices B λ ( α ) ( x ) satisfy the following product formulas.
B λ ( α + β ) ( x + y ) = B λ ( α ) ( x ) B λ ( β ) ( y ) = B λ ( β ) ( x ) B λ ( α ) ( y ) = B λ ( α ) ( y ) B λ ( β ) ( x ) .
Definition 2 and the inversion-type Formula (12) lead to the following result:
Theorem 6.
The generalized degenerate Pascal matrix of the first kind P λ [ x ] can be factorized in terms of B λ ( x ) as follows:
P λ [ x ] = B λ ( x ) H λ ,
where H λ is an ( n + 1 ) × ( n + 1 ) invertible matrix with entries
H i , j , λ = i i j ( 1 ) i j + 1 , λ i j + 1 , i j , 0 , o t h e r w i s e .
Proof. 
Let us consider n N 0 and 0 i , j n such that i j . From Definition 2 and the inversion-type Formula (12), we have
p i , j , λ ( x ) = i j ( x ) i j , λ = i j i j + 1 k = 0 i j i j + 1 k + 1 ( 1 ) k + 1 , λ B i j k , λ ( x ) = k = 0 i j i j k B i j k , λ ( x ) i i j ( 1 ) k + 1 , λ k + 1 .
Since the right hand member of (27) is the ( i , j ) -th entry of matrix product B λ ( x ) H λ , we conclude that (26) holds. □
The following example shows the validity of Theorem 6.
Example 1.
Let us consider n = 2 . It follows from Definition 1, (26), and a simple computation that
P λ [ x ] = 0 0 ( x ) 0 , λ 0 0 1 0 ( x ) 1 , λ 1 1 ( x ) 0 , λ 0 2 0 ( x ) 2 , λ 2 1 ( x ) 1 , λ 2 2 ( x ) 0 , λ = 0 0 B 0 , λ ( x ) 0 0 1 0 B 1 , λ ( x ) 1 1 B 0 , λ ( x ) 0 2 0 B 2 , λ ( x ) 2 1 B 1 , λ ( x ) 2 2 B 0 , λ ( x ) B λ ( x ) 0 0 ( 1 ) 1 , λ 0 0 1 1 ( 1 ) 2 , λ 2 1 0 ( 1 ) 1 , λ 0 2 2 ( 1 ) 3 , λ 3 2 1 ( 1 ) 2 , λ 2 2 0 ( 1 ) 1 , λ H λ
Definition 3.
The generalized degenerate ( n + 1 ) × ( n + 1 ) Euler matrix E λ ( α ) ( x ) is defined by the entries
E i , j , λ ( α ) ( x ) = i j E i j , λ ( α ) ( x ) , i j , 0 , o t h e r w i s e .
We denote by E λ ( x ) the degenerate Euler matrix E λ ( 1 ) ( x ) .
Definition 3 and the inversion-type Formula (17) lead to the following result:
Theorem 7.
The generalized degenerate Pascal matrix of the first kind P λ [ x ] can be factorized in terms of E λ ( x ) as follows:
P λ [ x ] = E λ ( x ) T λ ,
where T λ is an ( n + 1 ) × ( n + 1 ) invertible matrix with entries
T i , j , λ = i i j 1 + a i j ( λ ) ( 1 ) i j , λ 2 , i j , 0 , o t h e r w i s e .
Proof. 
Let us consider n N 0 and 0 i , j n such that i j . From Definition 3 and the inversion-type Formula (17), we have
p i , j , λ ( x ) = i j ( x ) i j , λ = 1 2 i j k = 0 i j i j k 1 + a k ( λ ) ( 1 ) k , λ E i j k , λ ( x ) = k = 0 i j i j k E i j k , λ ( x ) i j 1 + a k ( λ ) ( 1 ) k , λ 2 .
Since the right-hand member of (29) is the ( i , j ) -th entry of matrix product E λ ( x ) T λ , we conclude that (28) holds. □
Combining Theorems 6 and 7 gives the following connection formula.
Corollary 6.
For any λ , x R , we have
E λ ( x ) = B λ ( x ) H λ T λ 1 .
The next result is an immediate consequence of Definition 3 and the addition Formula (9).
Theorem 8.
The generalized degenerate Euler matrices E λ ( α ) ( x ) satisfy the following product formulas.
E λ ( α + β ) ( x + y ) = E λ ( α ) ( x ) E λ ( β ) ( y ) = E λ ( β ) ( x ) E λ ( α ) ( y ) = E λ ( α ) ( y ) E λ ( β ) ( x ) .
Proof. 
Let C i , j , λ ( α , β ) ( x , y ) be the ( i , j ) -th entry of the matrix product E λ ( α ) ( x ) E λ ( β ) ( y ) , then, by the addition Formula (9), we have
C i , j , λ ( α , β ) ( x , y ) = k = 0 n i k E i k , λ ( α ) ( x ) k j E k j , λ ( β ) ( y ) , n 0 = k = j i i j i j i k E i k , λ ( α ) ( x ) E k j , λ ( β ) ( y ) = i j k = 0 i j i j k E i j k , λ ( α ) ( x ) E k , λ ( β ) ( y ) , = i j E i j , λ ( α + β ) ( x + y ) , for i j ,
which implies the first equality of (30). The second and third equalities of (30) can be derived in a similar way. □
Corollary 7.
Let ( x 1 , , x k ) R k . For α j real or complex parameters, the generalized degenerate Euler matrices E λ ( α ) ( x ) satisfy the following product formulas, j = 1 , , k .
E λ ( α 1 + α 2 + + α k ) ( x 1 + x 2 + + x k ) = j = 1 k E λ ( α j ) ( x j ) .
Proof. 
The application of induction on k gives the desired result. □
Taking x = x 1 = x 2 = = x k and α = α 1 = α 2 = = α k , we obtain the following simple formula for the powers of the generalized degenerate Euler matrices E λ ( α ) ( x ) .
Corollary 8.
The generalized degenerate Euler matrices E λ ( α ) ( x ) satisfy the following identity.
E λ ( α ) ( x ) k = E λ ( α ) ( k x ) , k N .
Remark 7.
Analogously, the above corollaries hold, mutatis mutandis, for the generalized degenerate Bernoulli matrices. More precisely, from Theorem 5, and using the same assumptions as Corollaries 7 and 8, we obtain
B λ ( α 1 + α 2 + + α k ) ( x 1 + x 2 + + x k ) = j = 1 k B λ ( α j ) ( x j ) , B λ ( α ) ( x ) k = B λ ( α ) ( k x ) .
Theorem 9.
The generalized degenerate Euler matrices E λ ( α ) ( x ) satisfy the following relations.
E λ ( α ) ( x + y ) = E λ ( α ) ( x ) P λ [ y ] = P λ [ x ] E λ ( α ) ( y ) = E λ ( α ) ( y ) P λ [ x ] .
Proof. 
The substitution β = 0 into (30) yields
E λ ( α ) ( x + y ) = E λ ( α ) ( x ) E λ ( 0 ) ( y ) = E λ ( 0 ) ( x ) E λ ( α ) ( y ) = E λ ( α ) ( y ) E λ ( 0 ) ( x ) .
Since E λ ( 0 ) ( x ) = P λ [ x ] , we obtain
E λ ( α ) ( x + y ) = P λ [ x ] E λ ( α ) ( y ) .
A similar argument allows us to show that E ( α ) ( x + y ) = E ( α ) ( x ) P λ [ y ] and E λ ( α ) ( x + y ) = E λ ( α ) ( y ) P λ [ x ] . This completes the proof of (31). □

4. Conclusions

The aim of our research was to determine novel properties of generalized degenerate Bernoulli and Euler matrices. First, we focused our attention on some matrix-inversion formulae involving these matrices. Secondly, we showed some analytic properties for the generalized degenerate Pascal matrix of the first kind and provided some factorizations for the generalized degenerate Euler polynomial matrix in terms of the generalized degenerate Pascal matrix of the first kind.
Finally, it is worth mentioning that the use of the Cauchy product of the power series is the technique behind some of our formulations. This approach is not a novelty; however, it has been useful for generating new families of special polynomials (satisfying or not Appell-type conditions), even very recently. In this regard, we refer the interested reader to [25,26] and the references therein for a detailed exposition about very recent trends in this broad field.

Author Contributions

Conceptualization, J.H. and Y.Q.; methodology, J.H., D.P. and Y.Q.; formal analysis, J.H., D.P. and Y.Q.; investigation, J.H., D.P. and Y.Q.; writing—original draft preparation, Y.Q.; writing—review and editing, J.H., D.P. and Y.Q.; supervision, Y.Q.; project administration, Y.Q.; funding acquisition, J.H. and Y.Q. All authors have read and agreed to the published version of the manuscript.

Funding

The research of J. Hernández has been partially supported by the Fondo Nacional de Innovación y Desarrollo Científico y Tecnológico (FONDOCYT), Dominican Republic, under grant 2020-2021-1D1-135. The research of Y. Quintana has been partially supported by the grant CEX2019-000904-S funded by MCIN/AEI/10.13039/501100011033, and by the Madrid Government (Comunidad de Madrid-Spain) under the Multiannual Agreement with UC3M in the line of Excellence of University Professors (EPUC3M23), in the context of the Fifth Regional Programme of Research and Technological Innovation (PRICIT).

Data Availability Statement

Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank the reviewers for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, T.; Kim, D.S. Degenerate Bernstein polynomials. Rev. R. Acad. Cienc. Exactas Fís. Nat. Ser. A Math. RACSAM 2019, 113, 2913–2920. [Google Scholar] [CrossRef] [Green Version]
  2. Kim, T.; Kim, D.S.; Hyeon, S.-H.; Park, J.-W. Some new formulas of complete and incomplete degenerate Bell polynomials. Adv. Differ. Equ. 2021, 2021, 326. [Google Scholar] [CrossRef]
  3. Quintana, Y.; Ramírez, W.; Urieles, A. Generalized Apostol-type polynomial matrix and its algebraic properties. Math. Rep. 2019, 21, 249–264. [Google Scholar]
  4. Quintana, Y.; Ramírez, W.; Urieles, A. Euler matrices and their algebraic properties revisited. Appl. Math. Inf. Sci. 2020, 14, 583–596. [Google Scholar] [CrossRef]
  5. Bayat, M.; Teimoori, H. The linear algebra of the generalized Pascal functional matrix. Linear Algebra Appl. 1999, 295, 81–89. [Google Scholar] [CrossRef]
  6. Brawer, R.; Pirovino, M. The linear algebra of the Pascal matrix. Linear Algebra Appl. 1992, 174, 13–23. [Google Scholar] [CrossRef]
  7. Call, G.S.; Velleman, D.J. Pascal’s Matrices. Am. Math. Mon. 1993, 100, 372–376. [Google Scholar] [CrossRef]
  8. Can, M.; Cihat-Dağli, M. Extended Bernoulli and Stirling matrices and related combinatorial identities. Linear Algebra Appl. 2014, 444, 114–131. [Google Scholar] [CrossRef] [Green Version]
  9. Yang, Y.; Micek, C. Generalized Pascal functional matrix and its applications. Linear Algebra Appl. 2007, 423, 230–245. [Google Scholar] [CrossRef] [Green Version]
  10. Zhang, Z. The linear algebra of the generalized Pascal matrix. Linear Algebra Appl. 1997, 250, 51–60. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, T.; Kim, D.S. Identities involving degenerate Euler numbers and polynomials arising from non-linear differential equations. J. Nonlinear Sci. Appl. 2016, 9, 2086–2098. [Google Scholar] [CrossRef]
  12. Kim, T.; Kim, D.S. Identities of symmetry for degenerate Euler polynomials and alternating generalized falling factorial sums. Iran. J. Sci. Technol. Trans. Sci. 2017, 41, 939–949. [Google Scholar] [CrossRef]
  13. Kim, T.; Kim, D.S.; Jang, L.-C.; Lee, H.; Kim, H. Representations of degenerate Hermite polynomials. Adv. Appl. Math. 2022, 139, 102359. [Google Scholar] [CrossRef]
  14. Kim, T.; Kim, D.S.; Kim, H.K. On generalized degenerate Euler-Genocchi polynomials. Appl. Math. Sci. Eng. 2022, 31, 2159958. [Google Scholar] [CrossRef]
  15. Kim, T.; Kim, D.S. On some degenerate differential and degenerate difference operators. Russ. J. Math. Phys. 2022, 29, 37–46. [Google Scholar] [CrossRef]
  16. Carlitz, L. A degenerate Staudt-Clausen theorem. Arch. Math. 1956, 7, 28–33. [Google Scholar] [CrossRef]
  17. Carlitz, L. Degenerate Stirling, Bernoulli and Eulerian numbers. Util. Math. 1979, 15, 51–88. [Google Scholar]
  18. Chandragiri, S.; Shishkina, O.A. Generalized Bernoulli numbers and polynomials in the context of the Clifford analysis. J. Sib. Fed. Univ.-Math. Phys. 2018, 11, 127–136. [Google Scholar]
  19. Grigoriev, A.A.; Leinartas, E.K.; Lyapin, A.P. Summation of functions and polynomial solutions to a multidimensional difference equation. J. Sib. Fed. Univ.-Math. Phys. 2023, 16, 153–161. [Google Scholar]
  20. Leinartas, E.K.; Shishkina, O.A. The discrete analog of the Newton-Leibniz formula in the problem of summation over simplex lattice points. J. Sib. Fed. Univ.-Math. Phys. 2019, 12, 503–508. [Google Scholar] [CrossRef] [Green Version]
  21. Navas, L.; Ruiz, F.J.; Varona, J.L. Existence and reduction of generalized Apostol-Bernoulli, Apostol-Euler and Apostol-Genocchi polynomials. Arch. Math. 2019, 55, 157–165. [Google Scholar] [CrossRef]
  22. Nørlund, N.E. Vorlesungen über Differenzenrechnung; Springer: Berlin, Germany, 1924. [Google Scholar]
  23. Howard, F.T. Degenerate weighted Stirling numbers. Discrete Math. 1985, 57, 45–58. [Google Scholar] [CrossRef] [Green Version]
  24. Sándor, J.; Crstici, B. Handbook of Number Theory II; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2004. [Google Scholar]
  25. Albosaily, S.; Quintana, Y.; Iqbal, A.; Khan, W. Lagrange-based hypergeometric Bernoulli polynomials. Symmetry 2022, 14, 1125. [Google Scholar] [CrossRef]
  26. Quintana, Y. Generalized mixed type Bernoulli-Gegenbauer polynomial. Kragujev. J. Math. 2023, 47, 245–257. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hernández, J.; Peralta, D.; Quintana, Y. A Look at Generalized Degenerate Bernoulli and Euler Matrices. Mathematics 2023, 11, 2731. https://doi.org/10.3390/math11122731

AMA Style

Hernández J, Peralta D, Quintana Y. A Look at Generalized Degenerate Bernoulli and Euler Matrices. Mathematics. 2023; 11(12):2731. https://doi.org/10.3390/math11122731

Chicago/Turabian Style

Hernández, Juan, Dionisio Peralta, and Yamilet Quintana. 2023. "A Look at Generalized Degenerate Bernoulli and Euler Matrices" Mathematics 11, no. 12: 2731. https://doi.org/10.3390/math11122731

APA Style

Hernández, J., Peralta, D., & Quintana, Y. (2023). A Look at Generalized Degenerate Bernoulli and Euler Matrices. Mathematics, 11(12), 2731. https://doi.org/10.3390/math11122731

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