Next Article in Journal
Relative Controllability and Ulam–Hyers Stability of the Second-Order Linear Time-Delay Systems
Previous Article in Journal
The Sense of Cooperation on Interdependent Networks Inspired by Influence-Based Self-Organization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Some Functionals and Approximation Operators Associated with a Family of Discrete Probability Distributions

by
Ana Maria Acu
1,*,†,
Ioan Raşa
2,† and
Hari M. Srivastava
3,4,5,6,†
1
Department of Mathematics and Informatics, Lucian Blaga University of Sibiu, Str. Dr. I. Ratiu, No. 5-7, R-550012 Sibiu, Romania
2
Department of Mathematics, Technical University of Cluj-Napoca, Str. Memorandumului No. 28, R-400114 Cluj-Napoca, Romania
3
Department of Mathematics and Staristics, University of Victoria, Victoria, BC V8W 3R4, Canada
4
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
5
Department of Mathematics and Informatics, Azerbaijan University, 71 Jeyhun Hajibeyli Street, AZ1007 Baku, Azerbaijan
6
Center for Converging Humanities, Kyung Hee University, 26 Kyungheedae-ro, Dongdaemun-gu, Seoul 02447, Republic of Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(4), 805; https://doi.org/10.3390/math11040805
Submission received: 3 December 2022 / Revised: 29 January 2023 / Accepted: 31 January 2023 / Published: 5 February 2023

Abstract

:
A certain discrete probability distribution was considered in [“A discrete probability distribution and some applications”, Mediterr. J. Math., 2023]. Its basic properties were investigated and some applications were presented. We now embed this distribution into a family of discrete distributions depending on two parameters and investigate the properties of the new distributions.

1. Introduction

Let n N , j { 0 , 1 , , n } and α , β > 0 . Define
a n , j , α , β : = n j B ( j + α , n j + β ) B ( α , β ) ,
where B ( λ , μ ) is Euler’s Beta function defined by
B ( λ , μ ) : = 0 1 t λ 1 ( 1 t ) μ 1 d t min { ( λ ) , ( μ ) } > 1 = Γ ( λ ) Γ ( μ ) Γ ( λ + μ )
in terms of the familiar (Euler’s) Gamma function.
The probability distribution ( a n , j , α , α ) j = 0 , 1 , , n was investigated in [1]. In this paper, we extend the results from [1] to the distribution ( a n , j , α , β ) j = 0 , 1 , , n . This enlarges the family of the investigated distributions and the area of applications involving these distributions. In our present investigation, we are motivated also by several related recent developments on approximation operators and probability distributions by (for example) Ong et al. [2].
Let B n : C [ 0 , 1 ] C [ 0 , 1 ] be the classical Bernstein operators, defined as
B n f ( x ) : = j = 0 n n j x j ( 1 x ) n j f k n ( f C [ 0 , 1 ] , x [ 0 , 1 ] ) .
Consider the functional A n , α , β : C [ 0 , 1 ] R ,
A n , α , β ( f ) : = 1 B ( α , β ) 0 1 t α 1 ( 1 t ) β 1 B n f ( t ) d t .
Let e j ( t ) = t j ( j = 0 , 1 , ; t [ 0 , 1 ] ) .
It is well known that B n e 0 = e 0 and B n e 1 = e 1 , so that (2) yields
A n , α , β ( e 0 ) = 1 ,
A n , α , β ( e 1 ) = α α + β .
On the other hand, from (2) it is easy to infer that
A n , α , β ( f ) = j = 0 n a n , j , α , β f j n , f C [ 0 , 1 ] .
Combined with (3) and (4), (5) leads to
j = 0 n a n , j , α , β = 1 ,
j = 0 n j a n , j , α , β = n α α + β .
In particular, for each n N , α , β > 0 , ( a n , j , α , β ) j = 0 , 1 , , n can be considered as a discrete probability distribution, concentrated on a suitable set { x 0 , x 1 , , x n } .
In Section 2, besides the functionals A n , α , β , we consider the functional A α , β from (8). One of the main results is Theorem 1, which shows that for each f C [ 0 , 1 ] the sequence ( A n , α , β ( f ) ) n 1 converges to A α , β ( f ) . The rate of convergence is estimated for f C 2 [ 0 , 1 ] and the convergence for convex functions f is investigated. Equation (15) represents a quadrature formula for the following integral:
0 1 x α 1 ( 1 x ) β 1 f ( x ) d x .
The remainder is estimated for f C 2 [ 0 , 1 ] .
Section 3 is devoted to a sequence of random variables ( Z n , α , β ) n 1 . We describe the sequence of characteristic functions and its limit. Consequently, the sequence ( Z n , α , β ) n 1 converges in law to a Beta-type random variable. This offers a new proof for the convergence to zero of the remainder in the quadrature formula.
In Section 4, we consider two classical sequences of positive linear operators investigated by Lupaş and Lupaş [3] (see also [4,5]). We estimate the difference of these two sequences by using results from Section 2 and from the paper [6] and the references therein.
In Section 5, using the numbers a n , j , α , β , we construct a polynomial logarithmically convex Heun function.
Section 6 is devoted to inequalities between random variables from the preceding sections in the sense of the convex stochastic order. In the particular case α = 1 and β = 1 / 2 formula (33) was proved in [1].
In summary, our distribution ( a n , j , α , β ) j = 0 , 1 , , n has connections with several mathematical objects, including a sequence of positive linear functionals, a quadrature formula and its remainder, a sequence of random variables and their characteristic functions, two sequences of positive linear operators and the differences between them, polynomial logarithmically convex Heun functions, inequalities between random variables in the sense of the stochastic convex order.
In Section 7, we present conclusions and suggestions for further work.

2. A Quadrature Formula

Besides the functionals A n , α , β , consider also the functional A α , β : C [ 0 , 1 ] R ,
A α , β ( f ) : = 1 B ( α , β ) 0 1 t α 1 ( 1 t ) β 1 f ( t ) d t .
Theorem 1. 
(1) If f C [ 0 , 1 ] , then
lim n A n , α , β ( f ) = A α , β ( f ) .
(2) If f C [ 0 , 1 ] is a convex function, then
A n , α , β ( f ) A n + 1 , α , β ( f ) A α , β ( f ) , n 1 .
(3) If f C 2 [ 0 , 1 ] and 2 m f ( x ) 2 M , x [ 0 , 1 ] , then
m α β n ( α + β ) ( α + β + 1 ) A n , α , β ( f ) A α , β ( f ) M α β n ( α + β ) ( α + β + 1 ) .
In particular, if f C 2 [ 0 , 1 ] , then
| A n , α , β ( f ) A α , β ( f ) | α β n ( α + β ) ( α + β + 1 ) f .
Proof. 
1. It is well-known that
lim n B n f ( x ) = f ( x ) , f C [ 0 , 1 ] ,
uniformly with respect to x [ 0 , 1 ] .
Using (2), (8), (13), and the Lebesgue-dominated convergence theorem, we obtain (9).
2. It is also well-known that
B n f B n + 1 f f
for each convex function f C [ 0 , 1 ] .
To prove (10), we combine (2), (8) and (14).
3. If f C 2 [ 0 , 1 ] and 2 m f ( x ) 2 M , x [ 0 , 1 ] , then the functions f m e 2 and M e 2 f are convex. According to Item 2, we have
m A n , α , β ( e 2 ) A α , β ( e 2 ) A n , α , β ( f ) A α , β ( f )
M A n , α , β ( e 2 ) A α , β ( e 2 ) .
Since
A n , α , β ( e 2 ) A α , β ( e 2 ) = α β n ( α + β ) ( α + β + 1 ) ,
we obtain (11) and also (12). □
We now consider the following quadrature formula:
1 B ( α , β ) 0 1 x α 1 ( 1 x ) β 1 f ( x ) d x = j = 0 n a n , j , α , β f j n + R n , α , β ( f ) .
Theorem 2. 
The remainder R n , α , β satisfies
lim n R n , α , β ( f ) = 0 , f C [ 0 , 1 ] .
If f C 2 [ 0 , 1 ] and 2 m f ( x ) 2 M , x [ 0 , 1 ] , then
M α β n ( α + β ) ( α + β + 1 ) R n , α , β ( f ) m α β n ( α + β ) ( α + β + 1 ) .
In particular,
| R n , α , β ( f ) | α β n ( α + β ) ( α + β + 1 ) f .
Proof. 
Using (8) and (5), (15) can be written as
A α , β ( f ) = A n , α , β ( f ) + R n , α , β ( f ) .
Now (16), (17) and (18) are consequences of (9), (11) and (12). □

3. A Random Variable

Consider the random variable Z n , α , β defined by P Z n , α , β = j n = a n , j , α , β , j = 0 , 1 , , n .
According to (7), E Z n , α , β = α α + β , where E stands for mathematical expectation.
Let B α , β be the Beta-type random variable with density
t α 1 ( 1 t ) β 1 B ( α , β ) ( 0 < t < 1 ) 0 ( otherwise ) .
Theorem 3. 
Each of the following assertions holds true:
1. The characteristic function of Z n , α , β is given by
g n , α , β ( s ) = 1 B ( α , β ) 0 1 t α 1 ( 1 t ) β 1 1 + t ( e i s / n 1 ) n d t ( s R ) .
2. It is asserted that
lim n g n , α , β ( s ) = 0 1 t α 1 ( 1 t ) β 1 B ( α , β ) e i s t d t ( s R ) .
3. ( Z n , α , β ) n N converges in law to B α , β , as n .
Proof. 
The characteristic function Z n , α , β is by definition
g n , α , β ( s ) : = E e i s Z n , α , β = j = 0 n a n , j , α , β e i s j / n .
We have
j = 0 n a n , j , α , β z j = 1 B ( α , β ) j = 0 n z j n j 0 1 t j + α 1 ( 1 t ) n j + β 1 d t = 1 B ( α , β ) 0 1 t α 1 ( 1 t ) β 1 ( 1 + t ( z 1 ) ) n d t .
From (21) and (22) with z : = e i s / n , we get (19).
Now (20) is a consequence of (19) and Statement 3 follows from (20). □
Corollary 1. 
For f C [ 0 , 1 ] , it is asserted that
lim n E ( Z n , α , β ) = E ( B α , β ) .
Proof. 
The relation (23) is a consequence of 3) from Theorem 3. □
Remark 1. 
The relation (23) can be written as
lim n j = 0 n a n , j , α , β f j n = 1 B ( α , β ) 0 1 t α 1 ( 1 t ) β 1 f ( t ) d t ,
i.e.,
lim n A n , α , β ( f ) = A α , β ( f ) .
So, we have another proof of (9).

4. Two Sequences of Operators

Let B ¯ n : C [ 0 , 1 ] C [ 0 , 1 ] ,
B ¯ n f ( x ) : = f ( 0 ) , x = 0 , 1 B ( n x , n ( 1 x ) ) 0 1 t n x 1 ( 1 t ) n ( 1 x ) 1 f ( t ) d t ( 0 < x < 1 ) , f ( 1 ) ( x = 1 ) .
This operator was introduced by Mühlbach [7,8] and Lupaş [3,9].
We use the notation ( a ) k : = a ( a + 1 ) ( a + k 1 ) , k 1 , ( a ) 0 = 1 .
The operators L n : C [ 0 , 1 ] C [ 0 , 1 ] ,
L n f ( x ) : = j = 0 n n j ( n x ) j ( n n x ) n j ( n ) n f j n ,
were investigated in [3] (see also [4,5]). Clearly, L n f ( 0 ) = f ( 0 ) , L n f ( 1 ) = f ( 1 ) .
Theorem 4. 
If f C 2 [ 0 , 1 ] and 2 m f ( x ) 2 M , x [ 0 , 1 ] , then
m x ( 1 x ) n + 1 L n f ( x ) B ¯ n f ( x ) M x ( 1 x ) n + 1 , x [ 0 , 1 ] .
Proof. 
Clearly, (24) is satisfied for x = 0 and x = 1 . Let 0 < x < 1 . Then, according to (8),
B ¯ n f ( x ) = A n x , n ( 1 x ) ( f ) .
On the other hand, using (2) we obtain
A n , n x , n ( 1 x ) ( f ) = 1 B ( n x , n ( 1 x ) ) j = 0 n n j B ( j + n x , n j + n ( 1 x ) ) f j n = j = 0 n n j Γ ( j + n x ) Γ ( 2 n j n x ) Γ ( 2 n ) · Γ ( n ) Γ ( n x ) Γ ( n ( 1 x ) ) f j n = j = 0 n n j ( n x ) j ( n n x ) n j ( n ) n f j n ,
and so
L n f ( x ) = A n , n x , n ( 1 x ) ( f ) .
Now, (24) is a consequence of (25), (26) and (11). □

5. A Heun Function

Consider the function
h n , α , β ( x ) : = j = 0 n a n , j , α , β ( 2 x 1 ) 2 j , x R .
Theorem 5. 
If 0 < α β , then h n , α , β is a solution to the Heun differential equation
x ( x 1 ) x 1 2 u ( x ) + α + β x + α + β x 1 + 1 2 n 2 β x 1 / 2 u ( x ) 4 n α x 2 n α x ( x 1 ) ( x 1 / 2 ) u ( x ) = 0 .
Moreover, h n , α , β is a logarithmically convex function.
Proof. 
It was proved in [1] that if 0 < 2 θ γ , then the function
H n , γ , θ ( x ) : = j = 0 n n j ( θ ) j ( γ θ ) n j ( γ ) n ( 2 x 1 ) 2 j
is a logarithmically convex solution to the Heun differential equation
x ( x 1 ) x 1 2 u ( x ) + γ x + γ x 1 + 2 θ + 1 2 n 2 γ x 1 / 2 u ( x ) 4 n θ x 2 n θ x ( x 1 ) ( x 1 2 ) u ( x ) = 0 .
Setting θ = α , γ = α + β , we obtain
n j ( θ ) j ( γ θ ) n j ( γ ) n = a n , j , α , β ,
So, (29) becomes (28), and H n , γ , θ ( x ) becomes h n , α , β ( x ) . This concludes the proof. □
Remark 2. 
From (27) and (6), we see that h n , α , β ( 0 ) = 1 . Therefore h n , α , β is a polynomial, logarithmically convex, Heun function.
Remark 3. 
The function h n , α , β ( x ) = H n , α + β , α ( x ) can be expressed also in terms of Appell polynomials. For further details, (see [1], Section 6).

6. Stochastic Convex Orderings

Let X and Y be random variables on the same probability space. We say that X is dominated by Y (and write X c x Y ) in the sense of the convex stochastic order if
E f ( X ) E f ( Y )
for all convex functions f : R R such that the expectations exist (see [10,11]).
Theorem 6. 
Let 0 < β α and n N . Then, with respect to the convex stochastic order, we have
B α , α c x B β , β ,
and
Z n , α , α c x Z n , β , β .
Proof. 
It was proved that (see [1] (Theorem 10.1))
1 B ( α , α ) 0 1 x α 1 ( 1 x ) α 1 f ( x ) d x 1 B ( β , β ) 0 1 x β 1 ( 1 x ) β 1 f ( x ) d x ,
for each convex function f [ 0 , 1 ] , provided that 0 < β α .
This proves (30). Now, let f C [ 0 , 1 ] be convex. Then, B n f is also convex, and (32) shows that
1 B ( α , α ) 0 1 x α 1 ( 1 x ) α 1 B n f ( x ) d x 1 B ( β , β ) 0 1 x β 1 ( 1 x ) β 1 B n f ( x ) d x .
Using (2) and (5), we obtain
A n , α , α ( f ) A n , β , β ( f ) ,
i.e.,
j = 0 n a n , j , α , α f j n j = 0 n a n , j , β , β f j n .
This means that E f ( Z n , α , α ) E f ( Z n , β , β ) , therefore, (31) is proved. □

7. Conclusions and Directions for Further Work

The probability distribution ( a n , j , α ) j = 0 , 1 , , n was investigated from several points of view in [1]. In this paper, we generalize the corresponding results by considering the distribution ( a n , j , α , β ) j = 0 , 1 , , n such that ( a n , j , α , α ) j = 0 , 1 , , n is ( a n , j , α ) j = 0 , 1 , , n from [1]. A sequence of positive linear functionals is constructed in terms of the probability distribution. This sequence is convergent to another functional and this gives rise to a quadrature formula. The remainder of this formula is estimated for functions in C 2 [ 0 , 1 ] in terms of the uniform norm of f . We intend to extend this result to functions in C [ 0 , 1 ] by considering suitable moduli of continuity or K functionals. We also estimate the difference between two classical operators acting on functions in C 2 [ 0 , 1 ] and we study the same problem for functions in C [ 0 , 1 ] . A sequence of random variables is constructed using the probability distribution and is investigated from the point of view of the characteristic functions and their convergence. The probability distribution is useful for constructing a polynomial, logarithmically convex, Heun function. An inequality in the sense of the convex stochastic order is also established.
The numbers a n , j , α , β satisfy the following recurrence relations:
a 0 , 0 , α , β : = 1 , a n + 1 , 0 , α , β = n + β n + α + β a n , 0 , α , β , n 0 , a n + 1 , j + 1 , α , β = n + 1 n + α + β · j + α j + 1 a n , j , α , β , 0 j n .
For certain values of α , β , n , j , we present, in Table 1, the numerical values of the numbers a n , j , α , β .
We will study the possibility of extending the definitions of the numbers a n , j , α , β , when α and β tend individually or simultaneously to 0 or to . This would increase the family of results and examples related to the probability distribution which we have considered herein.

Author Contributions

These authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a Hasso Plattner Excellence Research Grant (LBUS-HPI-ERG-2020-07), financed by the Knowledge Transfer Center of the Lucian Blaga University of Sibiu.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Acu, A.M.; Raşa, I. A discrete probability distribution and some applications. Mediterr. J. Math. 2023, 20, 34. [Google Scholar] [CrossRef]
  2. Ong, S.H.; Ng, C.M.; Yap, H.K.; Srivastava, H.M. Some probabilistic generalizations of the Cheney-Sharma and Bernstein approximation operators. Axioms 2022, 10, 537. [Google Scholar] [CrossRef]
  3. Lupaş, A.; Lupaş, L. Polynomials of binomial type and approximation operators. Stud. Univ. Babeş-Bolyai Math. 1987, 32, 61–69. [Google Scholar]
  4. Lupaş, A. The Approximation by Means of Some Linear Positive Operators. In Approximation Theory, Proc. IDoMAT 95; Müller, M.W., Felten, M., Mache, D.H., Eds.; Mathematical Research; Academic Verlag: Berlin, Germany, 1995; Volume 86, pp. 201–229. [Google Scholar]
  5. Stancu, D.D. Approximation of functions by a new class of linear polynomial operators. Rev. Roum. Math. Pures Appl. 1968, 13, 1173–1194. [Google Scholar]
  6. Acu, A.M.; Raşa, I. Estimates for the differences of positive linear operators and their derivatives. Numer. Algor. 2020, 85, 191–208. [Google Scholar] [CrossRef]
  7. Mühlbach, G. Verallgemeinerungen der Bernstein- und der Lagrangepolynome, Bemerkungen zu einer Klasse linearer Polynomoperatoren von D.D. Stancu. Rev. Roum. Math. Pure Appl. 1970, 15, 1235–1252. [Google Scholar]
  8. Mühlbach, G. Rekursionsformeln für die zentralen Momente der Polya- und der Beta-Verteilung. Metrika 1972, 19, 171–177. [Google Scholar] [CrossRef]
  9. Lupaş, A. Die Folge der Beta-Operatoren. Ph.D. Thesis, Universität Stuttgart, Stuttgart, Germany, 1972. [Google Scholar]
  10. Rajba, T. On Some Recent Applications of Stochastic Convex Ordering Theorems to Some Functional Inequalities for Convex Functions: A Survey. In Developments in Functional Equations and Related Topics; Brzdek, J., Cieplinski, K., Rassias, T.M., Eds.; Springer Optimization and Its Applications; Springer: Cham, Switzerland, 2017; Chapter 11; Volume 124, pp. 231–274. [Google Scholar]
  11. Shaked, M.; Shanthikumar, J.G. Stochastic Orders; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
Table 1. Values of a n , j , α , β , α = 1 / 2 , β = 2 and n , j = 0 , 1 , 9 .
Table 1. Values of a n , j , α , β , α = 1 / 2 , β = 2 and n , j = 0 , 1 , 9 .
n j 0123456789
01
1 4 5 1 5
2 24 35 8 35 3 35
3 64 105 8 35 4 35 1 21
4 128 231 256 155 48 385 16 231 1 33
5 512 1001 640 3003 128 1001 80 1001 20 429 3 143
6 1024 2145 1024 5005 128 1001 256 3003 8 143 24 715 1 65
7 16384 36465 7168 36465 1536 12155 640 7293 448 7293 504 12155 28 1105 1 85
8 98304 230945 131072 692835 28672 230945 4096 46189 8960 138567 10752 230945 672 20995 32 1615 3 323
9 131072 323323 294912 1616615 196608 1616615 4096 46189 3072 46189 2304 46189 768 20995 288 11305 36 2261 1 133
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

Acu, A.M.; Raşa, I.; Srivastava, H.M. Some Functionals and Approximation Operators Associated with a Family of Discrete Probability Distributions. Mathematics 2023, 11, 805. https://doi.org/10.3390/math11040805

AMA Style

Acu AM, Raşa I, Srivastava HM. Some Functionals and Approximation Operators Associated with a Family of Discrete Probability Distributions. Mathematics. 2023; 11(4):805. https://doi.org/10.3390/math11040805

Chicago/Turabian Style

Acu, Ana Maria, Ioan Raşa, and Hari M. Srivastava. 2023. "Some Functionals and Approximation Operators Associated with a Family of Discrete Probability Distributions" Mathematics 11, no. 4: 805. https://doi.org/10.3390/math11040805

APA Style

Acu, A. M., Raşa, I., & Srivastava, H. M. (2023). Some Functionals and Approximation Operators Associated with a Family of Discrete Probability Distributions. Mathematics, 11(4), 805. https://doi.org/10.3390/math11040805

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