Next Article in Journal
A Mizoguchi–Takahashi Type Fixed Point Theorem in Complete Extended b-Metric Spaces
Next Article in Special Issue
Non-Parametric Threshold Estimation for the Wiener–Poisson Risk Model
Previous Article in Journal
On Some New Fixed Point Results in Complete Extended b-Metric Spaces
Previous Article in Special Issue
Estimating the Expected Discounted Penalty Function in a Compound Poisson Insurance Risk Model with Mixed Premium Income
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Rate of Convergence for a Characteristic of Multidimensional Birth-Death Process

1
Department of Applied Mathematics, Vologda State University, IPI FRC CSC RAS, VolSC RAS, 160000 Vologda, Russia
2
Department of Mathematics, Vologda State University, 160000 Vologda, Russia
3
Department of Applied Mathematics, Vologda State University, 160000 Vologda, Russia
4
Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, 119991 Moscow, Russia
5
Department of Mathematics, School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(5), 477; https://doi.org/10.3390/math7050477
Submission received: 15 April 2019 / Revised: 18 May 2019 / Accepted: 21 May 2019 / Published: 26 May 2019
(This article belongs to the Special Issue Stochastic Processes: Theory and Applications)

Abstract

:
We consider a multidimensional inhomogeneous birth-death process. In this paper, a general situation is studied in which the intensity of birth and death for each coordinate (“each type of particle”) depends on the state vector of the whole process. A one-dimensional projection of this process on one of the coordinate axes is considered. In this case, a non-Markov process is obtained, in which the transitions to neighboring states are possible in small periods of time. For this one-dimensional process, by modifying the method previously developed by the authors of the note, estimates of the rate of convergence in weakly ergodic and null-ergodic cases are obtained. The simplest example of a two-dimensional process of this type is considered.

1. Introduction and Preliminaries

Multidimensional birth-death processes (BDP) were objects of a number of studies in queueing theory and other applied fields. The authors of these papers studied different special classes of homogeneous multidimensional BDPs under some restrictions and considered fluid approximations [1], simulations [2,3,4,5,6], large deviations [7], stability [8,9], and other features. The problem of the product from solutions for such models was considered, for instance, in [10,11] (also, see the references therein). If the process is inhomogeneous and the transition intensities have a more general form, then the problem of computation of any probabilistic characteristics of the queueing model is much more difficult.
In the general case, it is impossible to obtain explicit solutions and their characteristics, as well as to construct any significant characteristics of the processes, as can be seen from the above list of works. This paper fills this gap and proposes a method of research and evaluation allowing one to estimate the rate of convergence for a one-dimensional projection of the multidimensional birth-death process. The approach also makes it possible to evaluate the main characteristics of the projection, as is demonstrated by the simplest example of an inhomogeneous two-dimensional process.
The background of our approach is the method of investigation of inhomogeneous BDP, see the detailed discussion and some preliminary results in [12,13,14,15]. Estimates for the state probabilities of one-dimensional projections of a multidimensional BDP were studied in [16,17]. However, within that methodology, it was impossible to obtain estimates of the rate of convergence, since the logarithmic norm of the operator cannot be applied to the corresponding nonlinear systems.
Here, we substantially modify that approach so that it can be used for estimation and construction of some explicit bounds on the rate of convergence for one-dimensional projection of a multidimensional BDP. Namely, in Section 2, we develop a simple but efficient method for bounding the rate of convergence for an arbitrary (which may be nonlinear, depending on the number of parameters and so on) differential equation in the space of sequences l 1 , and in Section 3, we apply this method to bounding the rate of convergence for one-dimensional projections of BDP.
Let X ( t ) = ( X 1 ( t ) , , X d ( t ) ) be a d-dimensional BDP such that in the interval ( t , t + h ) , the following transitions are possible with order h: birth of a particle of type j, death of a particle of type j.
Let λ j , m ( t ) be the corresponding birth rate (from the state m = ( m 1 , , m d ) = i = 1 d m i e i to the state m + e j ) and let μ j , m ( t ) be the corresponding death intensity (from the state m = ( m 1 , , m d ) = i = 1 d m i e i to the state m e j ). Denote p m ( t ) = Pr X ( t ) = m .
To consider the existence and uniqueness, we renumber the states (only in this section), transforming the process into a one-dimensional one. Now, let the (finite or countable) state space of the vector process under consideration be arranged in a special order, say 0 , 1 , . Denote by p i ( t ) , the corresponding state probabilities, and by p ( t ) , the corresponding column vector of state probabilities. Applying our standard approach (see details in [12,14,15]), we suppose in addition that all intensities are nonnegative functions locally integrable on [ 0 , ) , and, moreover, in new enumeration,
Pr X ( t + h ) = j / X ( t ) = i = q i j ( t ) h + α i j ( t , h ) , j i , 1 k i q i k ( t ) h + α i ( t , h ) , j = i ,
where q i j ( t ) are the corresponding transition intensities and all α i ( t , h ) are o ( h ) uniformly in i, that is, lim h 0 1 h sup i | α i ( t , h ) | = 0 , for any t 0 .
We suppose that λ j , m ( t ) L < , μ j , m ( t ) M < , for any j, m and almost all t 0 .
The probabilistic dynamics of the process is represented by the forward Kolmogorov system:
d p d t = A ( t ) p ( t ) ,
where A ( t ) is the corresponding infinitesimal (intensity) matrix.
Throughout the paper, we denote the l 1 -norm by · , i.e., x = | x i | , and B = sup j i | b i j | for B = ( b i j ) i , j = 0 .
Let Ω be the set all stochastic vectors, i.e., l 1 -vectors with nonnegative coordinates and unit norm. We have the inequality A ( t ) 2 d L + M < , for almost all t 0 . Hence, the operator function A ( t ) from l 1 into itself is bounded for almost all t 0 and is locally integrable on [ 0 ; ) . Therefore, we can consider (2) as a differential equation in the space l 1 with bounded operator.
It is well known, see [18], that the Cauchy problem for differential Equation (2) has unique solution for an arbitrary initial condition, and p ( s ) Ω implies p ( t ) Ω for t s 0 .
We recall that a Markov chain X ( t ) is called null-ergodic, if all p i ( t ) 0 as t for any initial condition, and it is called weakly ergodic, if p * ( t ) p * * ( t ) 0 as t for any initial condition p * ( 0 ) , p * * ( 0 ) , see for instance [12,14].

2. Bounds on the Rate of Convergence for a Differential Equation

Consider a general (linear or nonlinear) differential equation
d y d t = H y ( t ) ,
in the space of sequences l 1 under the assumption of existence and uniqueness of a solution for any initial condition y ( 0 ) .
Let H = ( h i j ) , where all h i j depend on some parameters (for instance, on y, t , ).
We have
d y i d t = h i i y i + j i h i j y j .
Now, if y i > 0 , then
d | y i | d t = d y i d t = h i i | y i | + j i h i j y j h i i | y i | + j i | h i j | | y j | ,
and if y i < 0 , then we also have
d | y i | d t = d y i d t = h i i y i j i h i j y j h i i | y i | + j i | h i j | | y j | .
Finally, using the continuity of all coordinates of the solution and the absolute convergence of all series, we obtain the estimate
d y d t = i d | y i | d t i h i i | y i | + j i | h i j | | y j | β * y ,
where
β * = sup i h i i + j i | h j i | .
Remark 1.
One can see that inequality (4) implies the bound
y ( t ) e 0 t β * d u y ( 0 ) .
Moreover, if H is bounded for any t linear operator function from l 1 to itself, then β * ( t ) = γ ( H ( t ) ) is the corresponding logarithmic norm of H ( t ) , see [12,13,14,15].
On the other hand, in a nonlinear situation, β * ( t ) yields a generalization of this notion.

3. Bounds on the Rate of Convergence for a Projection of Multidimensional BDP

Again, consider the forward Kolmogorov system (2) in the original vector form. Then, we have
d p m d t = l λ l , m e l ( t ) p m e l + l μ l , m + e l ( t ) p m + e l l λ l , m + μ l , m ( t ) p m ,
for any m .
In this section, we consider the one-dimensional process X j ( t ) for a fixed j. Denote x k ( t ) = Pr X j ( t ) = k . Then, x k ( t ) = m , m j = k p m ( t ) . The process X j ( t ) has nonzero jump rates only for unit jumps ( ± 1 ), namely, if X j ( t ) = k , then for small positive h only the jumps X j ( t + h ) = k ± 1 are possible with positive intensities, say λ ˜ k and μ ˜ k , respectively. Moreover, (7) implies the equalities
λ ˜ k x k ( t ) = m , m j = k λ j , m ( t ) p m ( t ) ,
μ ˜ k x k ( t ) = m , m j = k μ j , m ( t ) p m ( t ) ,
and hence
λ ˜ k = m , m j = k λ j , m ( t ) p m ( t ) m , m j = k p m ( t ) ,
and
μ ˜ k = m , m j = k μ j , m ( t ) p m ( t ) m , m j = k p m ( t ) .
Then, X j ( t ) is a (in general, non-Markovian) birth and death process with birth and death intensities λ ˜ k and μ ˜ k , respectively, (that is, it is a process with possible infinitesimal jumps ± 1 , the intensities of which depend on t and on the initial condition for the original multidimensional process X ( t ) .)
For any fixed initial distribution p ( 0 ) and any t > 0 , the probability distribution p ( t ) is unique. Hence, λ ˜ k = λ k p ( 0 ) , t and μ ˜ k = μ k p ( 0 ) , t uniquely define the system
d x d t = A ˜ x ( t ) ,
for the vector x ( t ) of state probabilities of the projection X j ( t ) under the given initial condition. Obviously, different initial conditions specify different systems.
Here, A ˜ is the corresponding three-diagonal “birth-death” transposed intensity matrix such that all off-diagonal elements are nonnegative and all column-wise sums are equal to zero.
Let for all m and any t 0
l j λ j , m ( t ) L j , m j μ j , m ( t ) M j .
Then, from (10) and (11), we obtain the two-sided bounds
l j λ ˜ k L j , m j μ ˜ k M j ,
for any k, any t, and any initial conditions.
1. Let the state space of X j ( t ) be countable and
M j < l j .
Put σ = M j / l j < 1 , δ n = σ n , n 0 , x ˜ n = δ n x n , and x ˜ = x ˜ 0 , x ˜ 1 , . Let Λ be a diagonal matrix, Λ = d i a g δ 0 , δ 1 , .
Note that in this situation, x ˜ ( t ) = i = 0 δ k x k ( t ) , and x ˜ ( t ) 0 as t implying null ergodicity of X j ( t ) , that is p k ( t ) = Pr X j ( t ) = k 0 as t for any k.
Then,
d x ˜ d t = Λ A ˜ Λ 1 x ˜ ( t ) .
Then, we have
λ ˜ k + μ ˜ k δ k + 1 δ k λ ˜ k δ k 1 δ k μ ˜ k λ ˜ k 1 σ μ ˜ k 1 / σ 1 l j 1 σ M j 1 / σ 1 = l j M j 2 = α * ,
implying the estimate
d x ˜ d t sup k δ k + 1 δ k λ ˜ k + δ k 1 δ k μ ˜ k λ ˜ k μ ˜ k x ˜ = inf k λ ˜ k + μ ˜ k δ k + 1 δ k λ ˜ k δ k 1 δ k μ ˜ k x ˜ α * x ˜ ,
and the following statement.
Theorem 1.
Let (15) hold for some j. Then, X j ( t ) is null-ergodic and the following bounds hold:
x ˜ ( t ) e α * t x ˜ ( 0 ) ,
and
Pr X j ( t ) n / X j ( 0 ) = k σ k n · e α * t .
Hence,
Pr X j ( t ) > n / X j ( 0 ) = k > 1 σ k n · e α * t ,
and Pr X j ( t ) > n / X j ( 0 ) = k 1 as t , for any n , k .
Remark 2.
It should be noted that the above requirements are imposed only on this one coordinate.
2. Let
L j < m j , α * = l j + m j 2 L j M j > 0 .
We have x ( t ) Ω for any t 0 . Set x 0 ( t ) = 1 i 1 x i ( t ) . Then, from (12), we obtain the system
d z d t = B ˜ z + f ˜ ,
where z = ( x 1 , x 2 , ) , f ˜ = λ ˜ 0 , 0 , 0 , , and the corresponding matrix B ˜ = b ˜ i j i , j = 1 , where b ˜ i j = a ˜ i j a ˜ i 0 for the corresponding elements of the matrix A ˜ .
For the solutions of system (23), the rate of convergence is determined by the system
d w d t = B ˜ w ,
where all elements of B ˜ depend on t and the initial condition of the original process.
Now, let β = M j L j > 1 in accordance with (22). Let d k + 1 = β k , k 0 . Denote by D, the upper triangular matrix
D = d 1 d 1 d 1 0 d 2 d 2 0 0 d 3 .
Let w ˜ = D w . Then, the following bound holds:
d w ˜ d t sup i 0 d i + 1 d i λ ˜ i + 1 + d i 1 d i μ ˜ i λ ˜ i + μ ˜ i + 1 w ˜ = inf i 0 λ ˜ i + μ ˜ i + 1 β λ ˜ i + 1 μ ˜ i / β w ˜ α * w ˜ .
Note that w ˜ = D w 1 2 w , see detailed discussion in [15], therefore, if w ˜ ( t ) 0 as t , then X j ( t ) is weakly ergodic.
Thus, we obtain the following statement.
Theorem 2.
Let (22) hold for some j. Then, X j ( t ) is weakly ergodic and the following bound holds:
D w ( t ) e α * t D w ( 0 ) ,
for any t 0 and any corresponding initial conditions.

4. Example

Consider a simple two-dimensional BDP with finite state space { i , j } , 0 i 10 , 0 j 10 and the following transition intensities:
(i)
λ 1 , i , 0 ( t ) = i + 1 11 λ 1 ( t ) from ( i , 0 ) to ( i + 1 , 0 ) ;
(ii)
λ 1 , i , j ( t ) = λ 1 ( t ) from ( i , j ) to ( i + 1 , j ) if j 0 ;
(iii)
λ 2 , i , j ( t ) = λ 2 ( t ) from ( i , j ) to ( i + 1 , j ) ;
(iv)
μ 1 , i , j ( t ) = μ 1 ( t ) from ( i , j ) to ( i 1 , j ) ;
(v)
μ 2 , i , j ( t ) = μ 2 ( t ) from ( i , j ) to ( i , j 1 ) ;
where λ 1 ( t ) = 1 + cos ( 2 π t ) , λ 2 ( t ) = 5 + sin ( 2 π t ) , μ 1 ( t ) = 11 + sin ( 2 π t ) , μ 2 = 3 .
Then, β = M 1 L 1 = 6 , and Theorem 2 gives bound (27) with α * = 10 4 6 .
We computed some important characteristics for the original process and its projection X 1 ( t ) , namely:
Figure 1, Figure 2 and Figure 3 show the behaviour of the state probabilities for X 1 ( t ) , namely Pr ( X 1 ( t ) = 0 ) , Pr ( X 1 ( t ) = 1 ) , and Pr ( X 1 ( t ) = 2 ) under two initial conditions for the original BDP:
(i)
p i , j ( 0 ) = 1 121 , for any i , j (blue); and
(ii)
p 0 , 0 ( 0 ) = 109 121 , p i , j ( 0 ) = 1 1210 , for any i , j such that i + j > 0 (green).
Note that the corresponding initial conditions for the projection are x ( 0 ) = 1 11 , , 1 11 T , and x ( 0 ) = 10 / 11 , 1 / 110 , , 1 / 110 T .
These Figures illustrate the rate of convergence in a weak ergodic situation.
Figure 4 and Figure 5 show the ’birth intensities’ λ ˜ 0 and λ ˜ 1 for X 1 ( t ) under the same initial conditions.
Note that all the quantities are found by numerically solving the Cauchy problem for the forward Kolmogorov system (2) and the corresponding system (12) for its projection on the corresponding interval.
As can be seen from Theorem 2 and the figures below, to construct all the characteristics of interest with good accuracy, it suffices to carry out the numerical solution on the interval [ 0 , 5 ] .
As was already noted, the projection of the original process is not a Markov process, and all probabilistic characteristics depend on the initial conditions of the original process.
It can be seen that these characteristics present comprehensive information concerning the behavior of the projection of the original process.

5. Conclusions

In the paper, some estimates of the rate of convergence were discussed for one-dimensional projections of multidimensional inhomogeneous birth and death processes. Some specific queueing models were considered. The applied approach allows one to use an analogue of the logarithmic norm of an operator function for a nonlinear system of differential equations, as was shown in Section 2. In addition, similar estimates can be obtained for other one-dimensional processes related to the original one. For example, the total number of “particles” of all types can be studied. Moreover, it is possible to study multidimensional processes with possible transformations of particles from one type to another. Such processes play a very important role in stochastic models of epidemics, see, for example, References [19,20] and the references therein.

Author Contributions

Conceptualization, A.Z., Y.S.; Methodology, A.Z., Y.S.; Software, Y.S.; Validation, A.Z., Y.S., K.K., V.K.; Investigation, A.Z., V.K.; Writing—Original Draft Preparation, A.Z., K.K., Y.S., V.K.; Writing—Review and Editing, A.Z., K.K., V.K.; Supervision, A.Z.; Project Administration, K.K.

Funding

The work of A.Z. and V.K. was supported by the Russian Science Foundation under grant 18-11-00155.

Acknowledgments

The authors thank the referees for useful comments that improved the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Michaelides, M.; Hillston, J.; Sanguinetti, G. Geometric fluid approximation for general continuous-time Markov chains. arXiv 2019, arXiv:1901.11417. [Google Scholar]
  2. Goswami, B.B.; Khouider, B.; Phani, R.; Mukhopadhyay, P.; Majda, A. Improving synoptic and intraseasonal variability in CFSv2 via stochastic representation of organized convection. Geophys. Res. Lett. 2017, 44, 1104–1113. [Google Scholar] [CrossRef]
  3. Goswami, B.B.; Khouider, B.; Phani, R.; Mukhopadhyay, P.; Majda, A.J. The Stochastic Multi-cloud Model (SMCM) Convective Parameterization in the CFSv2: Scopes and Opportunities. In Current Trends in the Representation of Physical Processes in Weather and Climate Models; Springer: Singapore, 2019; pp. 157–181. [Google Scholar]
  4. Deng, Q.; Khouider, B.; Majda, A.J. The MJO in a coarse-resolution GCM with a stochastic multicloud parameterization. J. Atmos. Sci. 2015, 72, 55–74. [Google Scholar] [CrossRef]
  5. Marsan, M.A.; Marano, S.; Mastroianni, C.; Meo, M. Performance analysis of cellular mobile communication networks supporting multimedia services. Mob. Netw. Appl. 2000, 5, 167–177. [Google Scholar] [CrossRef]
  6. Stevens-Navarro, E.; Mohsenian-Rad, A.H.; Wong, V.W. Connection admission control for multiservice integrated cellular/WLAN system. IEEE Trans. Veh. Technol. 2008, 57, 3789–3800. [Google Scholar] [CrossRef]
  7. Jonckheere, M.; Lypez, S. Large deviations for the stationary measure of networks under proportional fair allocations. Math. Oper. Res. 2013, 39, 418–431. [Google Scholar] [CrossRef]
  8. Jonckheere, M.; Shneer, S. Stability of multi-dimensional birth-and-death processes with state-dependent 0-homogeneous jumps. Adv. Appl. Probab. 2014, 46, 59–75. [Google Scholar] [CrossRef]
  9. Lee, C. On moment stability properties for a class of state-dependent stochastic networks. J. Korean Stat. Soc. 2011, 40, 325–336. [Google Scholar] [CrossRef]
  10. Mather, W.H.; Hasty, J.; Tsimring, L.S.; Williams, R.J. Factorized time-dependent distributions for certain multiclass queueing networks and an application to enzymatic processing networks. Queueing Syst. 2011, 69, 313–328. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Tsitsiashvili, G.S.; Osipova, M.A.; Koliev, N.V.; Baum, D. A product theorem for Markov chains with application to PF-queueing networks. Ann. Oper. Res. 2002, 113, 141–154. [Google Scholar] [CrossRef]
  12. Granovsky, B.; Zeifman, A. Nonstationary queues: Estimation of the rate of convergence. Queueing Syst. 2004, 46, 363–388. [Google Scholar] [CrossRef]
  13. Zeifman, A.I. On the estimation of probabilities for birth and death processes. J. Appl. Probab. 1995, 32, 623–634. [Google Scholar] [CrossRef] [PubMed]
  14. Zeifman, A.I. Upper and lower bounds on the rate of convergence for nonhomogeneous birth and death processes. Stoch. Proc. Appl. 1995, 59, 157–173. [Google Scholar] [CrossRef]
  15. Zeifman, A.; Leorato, S.; Orsingher, E.; Satin, Y.; Shilova, G. Some universal limits for nonhomogeneous birth and death processes. Queueing Syst. 2006, 52, 139–151. [Google Scholar] [CrossRef] [Green Version]
  16. Zeifman, A.; Sipin, A.; Korolev, V.; Bening, V. Estimates of some characteristics of multidimensional birth-and-death processes. Dokl. Math. 2015, 92, 695–697. [Google Scholar] [CrossRef]
  17. Zeifman, A.I.; Sipin, A.S.; Korotysheva, A.V.; Panfilova, T.L.; Satin, Y.A.; Shilova, G.N.; Korolev, V.Y. Estimation of Probabilities for Multidimensional Birth-Death Processes. J. Math. Sci. 2016, 218, 238–244. [Google Scholar] [CrossRef]
  18. Daleckij, J.L.; Krein, M.G. Stability of solutions of differential equations in Banach space. Am. Math. Soc. Transl. 2002, 43, 1–386. [Google Scholar]
  19. Britton, T. Stochastic epidemic models: A survey. Math. Biosci. 2010, 225, 24–35. [Google Scholar] [CrossRef] [PubMed]
  20. Zeifman, A.I. Processes of birth and death and simple stochastic epidemic models. Autom. Remote Control 1985, 6, 128–135. [Google Scholar]
Figure 1. Pr ( X 1 ( t ) = 0 ) under initial conditions (i) and (ii).
Figure 1. Pr ( X 1 ( t ) = 0 ) under initial conditions (i) and (ii).
Mathematics 07 00477 g001
Figure 2. Pr ( X 1 ( t ) = 1 ) under initial conditions (i) and (ii).
Figure 2. Pr ( X 1 ( t ) = 1 ) under initial conditions (i) and (ii).
Mathematics 07 00477 g002
Figure 3. Pr ( X 1 ( t ) = 2 ) under initial conditions (i) and (ii).
Figure 3. Pr ( X 1 ( t ) = 2 ) under initial conditions (i) and (ii).
Mathematics 07 00477 g003
Figure 4. λ ˜ 0 under initial conditions (i) and (ii).
Figure 4. λ ˜ 0 under initial conditions (i) and (ii).
Mathematics 07 00477 g004
Figure 5. λ ˜ 1 under initial conditions (i) and (ii).
Figure 5. λ ˜ 1 under initial conditions (i) and (ii).
Mathematics 07 00477 g005

Share and Cite

MDPI and ACS Style

Zeifman, A.; Satin, Y.; Kiseleva, K.; Korolev, V. On the Rate of Convergence for a Characteristic of Multidimensional Birth-Death Process. Mathematics 2019, 7, 477. https://doi.org/10.3390/math7050477

AMA Style

Zeifman A, Satin Y, Kiseleva K, Korolev V. On the Rate of Convergence for a Characteristic of Multidimensional Birth-Death Process. Mathematics. 2019; 7(5):477. https://doi.org/10.3390/math7050477

Chicago/Turabian Style

Zeifman, Alexander, Yacov Satin, Ksenia Kiseleva, and Victor Korolev. 2019. "On the Rate of Convergence for a Characteristic of Multidimensional Birth-Death Process" Mathematics 7, no. 5: 477. https://doi.org/10.3390/math7050477

APA Style

Zeifman, A., Satin, Y., Kiseleva, K., & Korolev, V. (2019). On the Rate of Convergence for a Characteristic of Multidimensional Birth-Death Process. Mathematics, 7(5), 477. https://doi.org/10.3390/math7050477

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