Next Article in Journal
A Spatial Agent-Based Model for Studying the Effect of Human Mobility Patterns on Epidemic Outbreaks in Urban Areas
Next Article in Special Issue
Properties of the SURE Estimates When Using Continuous Thresholding Functions for Wavelet Shrinkage
Previous Article in Journal
The Concept of Topological Derivative for Eigenvalue Optimization Problem for Plane Structures
Previous Article in Special Issue
The Implicit Euler Scheme for FSDEs with Stochastic Forcing: Existence and Uniqueness of the Solution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On One Approach to Obtaining Estimates of the Rate of Convergence to the Limiting Regime of Markov Chains

1
Department of Applied Mathematics, Vologda State University, 15 Lenina Str., 160000 Vologda, Russia
2
Federal Research Center “Computer Science and Control”, Russian Academy of Sciences, 44-2 Vavilova Str., 119333 Moscow, Russia
3
Vologda Research Center, Russian Academy of Sciences, 556A Gorky Str., 160014 Vologda, Russia
4
Moscow Center for Fundamental and Applied Mathematics, Moscow State University, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(17), 2763; https://doi.org/10.3390/math12172763
Submission received: 25 July 2024 / Revised: 4 September 2024 / Accepted: 5 September 2024 / Published: 6 September 2024

Abstract

:
We revisit the problem of the computation of the limiting characteristics of (in)homogeneous continuous-time Markov chains with the finite state space. In general, it can be performed only numerically. The common rule of thumb is to interrupt calculations after quite some time, hoping that the values at some distant time interval will represent the sought-after solution. Convergence or ergodicity bounds, when available, can be used to answer such questions more accurately; i.e., they can indicate how to choose the position and the length of that distant time interval. The logarithmic norm method is a general technique that may allow one to obtain such bounds. Although it can handle continuous-time Markov chains with both finite and countable state spaces, its downside is the need to guess the proper similarity transformations, which may not exist. In this paper, we introduce a new technique, which broadens the scope of the logarithmic norm method. This is achieved by firstly splitting the generator of a Markov chain and then merging the convergence bounds of each block into a single bound. The proof of concept is illustrated by simple examples of the queueing theory.

1. Introduction

The number of theoretical and application-oriented results available in the literature concerning Markov chains can (and will) never be counted. The task does not become simpler even if one narrows down the search to the continuous-time Markov chains (CTMCs), which are homogeneous with a finite-state space. Nevertheless, the richness of this mathematical model and its widespread applications give rise to new and old research questions (see, for example, the recent review [1]), which require new, improved, or clarified solutions.
In this paper, we revisit the well-known problem of the computation of the limiting distributions of CTMCs with possibly time-varying transition intensities (for a review, one can refer, for example, to [2,3] and the references therein). Such problems arise in various operation research applications (see [4,5]). But here, for illustration purposes, we will make use of simple examples from the queueing theory, as its models (see, for example, refs. [2,6,7,8,9,10,11,12,13]) fit the mathematical framework of such Markov chains excellently.
The computation of limiting distributions usually has to be performed numerically, without knowing the limiting distribution itself. Therefore, an algorithm must be “told” when to stop the computation, i.e., when the limiting regime is reached. One of the methods which may provide such an answer in the form of an ergodicity (convergence) bound for the considered type of CTMC is the logarithmic norm method. (although ergodic properties of Markov chains have been the subject of extensive research (see, for example, refs. [14,15,16,17,18,19]), practically useful general ergodicity bounds are difficult to obtain and it remains, to a large extent, an open problem).
This implies (see, for example, refs. [20,21,22]) that (under certain conditions on the generator of the CTMC, say, A ( t ) = ( a i j ( t ) ) ) for the vector x ( t ) of state probabilities, which satisfies the system of ordinary differential equations (ODEs), d d t x ( t ) = A ( t ) x ( t ) , t 0 , the upper bound. Henceforth, · is the l 1 -norm; i.e., if x is an ( S + 1 ) -dimensional column vector, then x = k = 0 S | x k | . In the case of a linear operator, the norm is induced by the l 1 -norm on column vectors, for example, A ( t ) = sup 0 j S i = 0 S | a i j ( t ) | . x ( t ) e 0 t α ( u ) d u x ( 0 ) holds, where α ( t ) = sup i a i i ( t ) + j i | a j i ( t ) | is called the logarithmic norm of A ( t ) . The pros and cons of the method can be readily seen. On the one hand, it is quite general and the bounds can be simply computed. On the other hand, it can happen (and in the most interesting cases, it is so) that α ( t ) > 0 and thus the bound is meaningless. Therefore, certain techniques were developed (see, for example, refs. [23,24,25,26]), which allow one to obtain the convergence bounds but in weighted norms. Nevertheless, the success that can be achieved along these lines remains limited due to the fact that the weights need to be guessed for each new CTMC without a “rule of thumb”; they do not have any probabilistic meaning and can be considered as analogues of Lyapunov functions.
In this paper, we present the new technique, which expands the scope of the logarithmic norm method. It implies that in order to obtain the convergence bound for the given CTMC, one firstly partitions its infinitesimal matrix into blocks and the ergodicity properties of the blocks are studied independently. Then, the results are merged into one single ergodicity bound for the original CTMC. Whenever the probability flows between the blocks are not high, the merged bound is good. This technique does not solve the problem of choosing proper weights. But partitioning makes the search for them easier and sometimes unnecessary. Yet it is a matter of pure guessing regarding how to perform the partitioning, and there is no guarantee that the resulting (merged) ergodicity bound will be meaningful (see [27]).
The use of this new technique allows one to study the probabilistic characteristics of quasi-birth–death (QBD) processes, which have become very popular especially in the field of queueing theory. The first works on this class of random processes appeared in the late 1960s, when the first publications of Evans [28] and Wallace [29] appeared and further research was carried out by Neuts [30,31,32,33] and Ost [34]. Ergodicity bounds for QBDs, using easy-to-use techniques like logarithmic norm, usually cannot be found. But, as experiments show, using the technique proposed in this paper, sometimes finite QBD generators can be split into blocks, and the ergodicity condition established for each block can be merged into a single meaningful bound.
In Section 2, we introduce notation and in Section 3 and Section 4, we provide the description of the new technique. Section 5, Section 6 and Section 7 are devoted to special cases, which demonstrate the applicability of the technique and its connections to the known results in the literature.

2. Preliminaries

Let { X ( t ) , t 0 } be a (possibly inhomogeneous) CTMC with a finite state space X = { 0 , 1 , , S } . The transition probabilities will be denoted by p i j ( s , t ) , i.e., p i j ( s , t ) = Pr X ( t ) = j X ( s ) = i , i , j 0 , 0 s t . Let the state probability also be denoted by p i ( t ) = Pr X ( t ) = i , and the vector of state probabilities by p ( t ) = p 0 ( t ) , p 1 ( t ) , , p S ( t ) T . In what follows, it is assumed that
Pr X t + h = j | X t = i = q i j t h + α i j t , h , if   j i 1 k i q i k t h + α i t , h , if   j = i ,
where all α i ( t , h ) are o ( h ) uniformly in i, that is, sup i | α i ( t , h ) | = o ( h ) . Moreover, if X ( t ) is inhomogeneous, then all its infinitesimal characteristics (intensity functions) q i j t are integrable in t on any interval [ a , b ] , 0 a b .
Denote a i j ( t ) = q j i ( t ) as j i and a i i ( t ) = j i a j i ( t ) = j i q i j ( t ) .
Further, in order to be able to obtain tighter estimates, it is assumed that
| a i i ( t ) | L <
for almost all t 0 . The state probabilities then satisfy the forward Kolmogorov system:
d d t p ( t ) = A ( t ) p ( t ) , t 0 ,
where the initial distribution p ( 0 ) is given, A ( t ) = Q T ( t ) , and Q ( t ) is the infinitesimal matrix of the CTMC.
Let · be the usual l 1 -norm, i.e., x = | x i | , and B = sup j i | b i j | for B = ( b i j ) i , j = 0 . Denote Ω = x : x l 1 + & x = 1 .
Therefore, A ( t ) = 2 sup k a k k ( t ) 2 L for almost all t 0 , and we can apply the results of [35] to the Equation (3) in the space l 1 . Namely, in [35] it was shown that the Cauchy problem for the Equation (3) has a unique solution for an arbitrary initial condition. Moreover, if p ( s ) Ω , then p ( t ) Ω , for any 0 s t and any initial condition p ( s ) .
In what follows, we will make use of the fact that the subtraction of a given function from any row of the matrix A ( t ) does not change the CTMC’s convergence rate (see, for example, ref. [26]). Indeed, let 1 denote the vector of ones and ⊗ be the Kronecker product. Then, for the vector r ( t ) = ( r 0 ( t ) , r 1 ( t ) , , r S ( t ) T of arbitrary functions r i ( t ) , the system (3) can be rewritten in the following form:
d d t p ( t ) = A ( t ) 1 r ( t ) p ( t ) + 1 r ( t ) p ( t ) .
Since i = 0 S p i ( t ) = 1 for any t 0 , one has that 1 r ( t ) p ( t ) = r ( t ) . Therefore, for any two solutions p * ( t ) and p * * ( t ) of the system (3), one has that
d d t p * ( t ) = ( A ( t ) 1 r ( t ) ) p * ( t ) + r ( t ) ,
d d t p * * ( t ) = ( A ( t ) 1 r ( t ) ) p * * ( t ) + r ( t ) ,
and hence,
d d t x ( t ) = A ˜ ( t ) x ( t ) ,
where x ( t ) = p * ( t ) p * * ( t ) and A ˜ ( t ) = A ( t ) 1 r ( t ) .
Henceforth, it will be convenient to use the block representation of A ˜ ( t ) for the matrix (which is always feasible).
A ˜ ( t ) = C ˜ ( t ) D ˜ ( t ) E ˜ ( t ) F ˜ ( t ) ,
in which the diagonal blocks are square matrices of size m > 0 and N > 0 , respectively, subject to the constraint m + N = S + 1 .

3. Upper Bounds for the Rate of Convergence

In the previous section, it was shown that the rate of convergence of the system d d t p ( t ) = A ( t ) p ( t ) coincides with the rate of convergence of the system d d t x ( t ) = A ˜ ( t ) x ( t ) , where A ˜ ( t ) is given by (5). By denoting
x ( t ) = ( x 0 ( t ) , x 1 ( t ) , , x m 1 ( t ) , x m ( t ) , x m + 1 ( t ) , , x S ( t ) ) T , x 0 ( t ) = ( x 0 ( t ) , x 1 ( t ) , , x m 1 ( t ) ) T , x 1 ( t ) = ( x m ( t ) , x m + 1 ( t ) , , x S ( t ) ) T ,
this new system can be split into the following two coupled systems of equations:
d d t x 0 ( t ) = C ˜ ( t ) x 0 ( t ) + D ˜ ( t ) x 1 ( t ) ,
d d t x 1 ( t ) = E ˜ ( t ) x 0 ( t ) + F ˜ ( t ) x 1 ( t ) .
The solutions of these systems can be written in the form
x 0 ( t ) = U 0 ( t , 0 ) x 0 ( 0 ) + 0 t U 0 ( t , s ) D ˜ ( s ) x 1 ( s ) d s ,
x 1 ( t ) = U 1 ( t , 0 ) x 1 ( 0 ) + 0 t U 1 ( t , s ) E ˜ ( s ) x 0 ( s ) d s .
Assume without the loss of generality that U 0 ( t , 0 ) U 1 ( t , 0 ) and K = sup t D ˜ ( t ) ,   L = sup t E ˜ ( t ) . Let x 1 ( 0 ) = 0 , where 0 denotes the vector of zeros. By sequentially substituting the expressions for x 0 ( t ) and x 1 ( t ) , given by (8) and (9), into the integrals on the right-hand side of (8) and (9), and then taking norms of both sides, one obtains the following upper bound for x 0 ( t ) :
x 0 ( t ) U 0 ( t , 0 )   x 0 ( 0 ) + 0 t U 0 ( t , s 1 )   D ˜ ( s 1 )   U 1 ( s 1 , 0 )   x 1 ( 0 ) d s 1 +                           U 1 ( t , 0 )   x 0 ( 0 ) + 0 t U 1 ( t , s 1 )   D ˜ ( s 1 )   U 1 ( s 1 , 0 )   x 1 ( 0 ) d s 1 + U 1 ( t , 0 ) 1 + K L t 2 2 ! + K L 2 t 4 4 ! + x 0 ( 0 ) .
Since, due to (2), both K and L are finite, from the last inequality, it follows that
x 0 ( t ) U 1 ( t , 0 ) cosh t K L x 0 ( 0 ) .
Accordingly, by putting x 0 ( 0 ) = 0 , one has that
x 1 ( t ) U 1 ( t , 0 ) cosh t K L x 1 ( 0 ) .
Now, since the vector x ( t ) can be represented as x ( t ) = 1 0 x 0 ( t ) + 0 1 x 1 ( t ) and
x 0 ( t ) = 1 0 x 0 ( t ) , x 1 ( t ) = 0 1 x 1 ( t ) ,
we obtain from (11) and (12)
x ( t ) 2 U 1 ( t , 0 ) cosh t K L x ( 0 ) .
The derivations made above can be summarized in the following theorem.
Theorem 1. 
Let there exist a function α ( t ) and a positive M such that the rate of convergence for the matrices C ˜ ( t ) and F ˜ ( t ) satisfies the inequalities U 0 ( t , 0 ) U 1 ( t , 0 ) e 0 t α ( s ) d s . Then, if 0 ( α ( s ) K L ) d s = + for some K sup t D ˜ ( t ) and L sup t E ˜ ( t ) , the CTMC { X ( t ) , t 0 } is exponentially ergodic and the bound
p * ( t ) p * * ( t ) 2 M e α t cosh t K L p * ( 0 ) p * * ( 0 )
holds for any two initial conditions p * ( 0 ) and p * * ( 0 ) .
A number of corollaries can be drawn from the theorem. For example, the bound (14) is applicable even if the initial distributions are unknown, since p * ( 0 ) p * * ( 0 ) 2 . Next, making use of the assumption that the state space of the CTMC is finite, one immediately obtains the bound 4 M S n e α t cosh t K L for the rate of convergence of any positive moment n of { X ( t ) , t 0 } . This enables one to use concentration inequalities to obtain some insight into the the behaviour of the probabilities of events related to { X ( t ) , t 0 } .

4. Relationship between Cauchy Operators

Denote the Cauchy operator of the ODE system d d t x ( t ) = A ( t ) x ( t ) , t 0 , by U x ( t , s ) . Thus, x ( t ) = U x ( t , 0 ) x ( 0 ) . For y ( t ) = T x ( t ) , where T is any linear transform, we have that d d t y ( t ) = T A ( t ) T 1 y ( t ) and y ( t ) = U y ( t , 0 ) y ( 0 ) , where U y ( t , s ) is the corresponding Cauchy operator. Assuming that T has a unique inverse, we have that
x ( t ) = T 1 y ( t ) = T 1 U y ( t , 0 ) T T 1 y ( 0 )
and therefore,
x ( t ) = T 1 U y ( t , 0 ) T x ( 0 ) .
Due to the uniqueness of the solution of the initial ODE system, we have that U x ( t , 0 ) = T 1 U y ( t , 0 ) T . Therefore, the Cauchy operators of the two systems satisfy the following inequalities:
U x ( t , 0 ) T 1   T   U y ( t , 0 ) ,
U y ( t , 0 ) T 1   T   U x ( t , 0 ) .
Since U y ( t , 0 )   exp γ ( T A ( t ) T 1 ) , the following theorem holds.
Theorem 2. 
If the matrix T has a unique inverse, the Cauchy operator of the ODE system d d t x ( t ) = A ( t ) x ( t ) satisfies the following inequality:
U ( t , 0 )     T 1   T exp γ ( T A ( t ) T 1 ) .

5. Homogeneous BDPs

Let us consider a homogeneous BDP process with birth and death rates λ and μ , respectively, and S + 1 states. From its infinitesimal matrix,
A = λ μ 0 0 0 λ ( λ + μ ) μ 0 0 0 λ ( λ + μ ) μ 0 0 0 λ ( λ + μ ) 0 μ 0 0 0 0 μ ,
one can obtain the matrix A ˜ by (5), for example, by subtracting μ from the first row, and λ from the second row. Therefore, the blocks of the matrix A ˜ will have the following form:
C ˜ = λ μ , D ˜ = ( 0 , μ , μ , , μ ) , E ˜ = ( 0 , 0 , , 0 ) T ,
F ˜ = ( 2 λ + μ ) μ λ λ λ ( λ + μ ) 0 μ 0 0 μ .
Consider the matrix T of the form
T = 1 1 1 1 0 μ λ μ λ μ λ 0 0 μ λ 2 μ λ 2 0 0 0 μ λ S 1 .
Then, the similarity transformation of the matrix F ˜ yields
T F ˜ T 1 = λ μ μ λ 0 0 μ λ ( 2 λ + μ ) μ λ 0 0 μ λ ( λ + μ ) μ λ 0 0 0 λ μ .
Since K = max i E ˜ = 0 and L = max i D ˜ = μ , the logarithmic norm of the matrix F ˜ is equal to ( μ λ ) 2 . Therefore, from (14), we obtain the bound
p * ( t ) p * * ( t ) 2 μ λ S 1 max 1 , 2 λ μ e ( μ λ ) 2 t p * ( 0 ) p * * ( 0 ) ,
for any two initial conditions p * ( 0 ) and p * * ( 0 ) . We note that the exponential rate of convergence ( μ λ ) 2 coincides with the well-known results in the literature [36].

6. One Example of the Quasi-Birth–Death Process

Consider the model of a production line (see Model 3 in [37]) as a single server with the finite capacity FIFO-queue, single Poisson flow of customers of intensity λ ( t ) , and Erlang-2 service times with the intensity μ ( t ) at each stage. It gives rise to the QBD on the state space E = { ( i , j ) , 0 i M 1 j n } , with the following infinitesimal generator (the argument t is omitted for brevity):
Q = B 0 R 0 0 0 0 A B R 0 0 0 0 A B R 0 0 0 0 0 A B R 0 0 0 0 A B M ,
where
A = 0 0 μ 0 , B 0 = ( λ + μ ) μ μ ( λ + μ ) , B = B 0 A ,
R = λ 0 0 λ , B M = B + R .
Let us split the infinitesimal generator A ( t ) = Q T ( t ) into four blocks of the following form:
C ( t ) = λ μ μ μ λ μ ,
D ( t ) = 0 μ 0 0 0 0 0 0 0 0 0 0 ,
E ( t ) = λ 0 0 λ 0 0 0 0 ,
F ( t ) = λ μ 0 0 μ 0 0 0 0 μ λ μ 0 0 0 0 0 0 λ 0 λ μ 0 0 μ 0 0 0 λ μ λ μ 0 0 0 0 0 0 λ 0 λ μ 0 0 0 0 0 0 λ μ λ μ 0 0 0 0 0 0 0 0 μ 0 0 0 0 0 0 0 μ μ .
By subtracting μ ( t ) 2 from both rows of the matrix C ( t ) , we obtain
C ˜ ( t ) = λ 3 2 μ 1 2 μ 1 2 μ λ 3 2 μ ,
D ˜ ( t ) = μ 2 μ 2 μ 2 μ 2 μ 2 μ 2 μ 2 μ 2 μ 2 μ 2 μ 2 μ 2 .
E ˜ ( t ) = E ( t ) , F ˜ ( t ) = F ( t ) .
Let us find the rate of convergence to the limiting regime for the block F ˜ ( t ) . Fix h > 0 and k > 0 , and consider the square matrix T of size N of the following form:
T = 0 1 0 0 0 0 0 0 k 0 0 0 0 0 0 0 0 0 0 h 0 0 0 0 0 0 h k 0 0 0 0 0 0 0 0 0 0 h 2 0 0 0 0 0 0 h 2 k 0 0 0 0 0 0 0 0 0 0 h 2 0 0 0 0 0 0 h 2 k 0 .
Then, the similarity transformation T F ˜ T 1 yields the following matrix:
T F ˜ T 1 = λ μ μ k 0 0 0 0 0 0 μ λ μ k μ h 0 0 0 0 0 h λ 0 λ μ μ k 0 0 0 0 0 h λ μ λ μ k μ h 0 0 0 0 0 h λ 0 λ μ μ k 0 0 0 0 0 h λ μ λ μ 0 0 0 0 0 0 0 0 μ μ k 0 0 0 0 0 0 μ μ ,
the logarithmic norm of which is equal to
γ ( T F ˜ T 1 ) = max λ μ + μ k + h λ , λ μ + k μ h + h λ .
Now, if one puts h = k 2 , k > 1 , then the rate of convergence to the limiting regime is equal to λ μ + μ k + λ k 2 . Therefore, for the norm U 1 ( t , 0 ) of the Cauchy operator, one has the upper bounded (note that the size of the block N coincides with the size of the block F ˜ ) by
U 1 ( t , 0 ) k N 1 exp 0 t λ ( s ) μ ( s ) + μ ( s ) k + λ ( s ) k 2 d s .
The convergence rate for the block C ˜ ( t ) is equal to λ μ and therefore, the Cauchy operator U 0 ( t , 0 ) is bounded by
U 0 ( t , 0 ) exp 0 t λ ( s ) μ ( s ) d s .
Clearly, K = sup t D ˜ = sup t μ ( t ) and L = sup t E ˜ = sup t λ ( t ) . Since U 1 ( t , 0 )     U 0 ( t , 0 ) , the Theorem 1 yields the following bound:
p * ( t ) p * * ( t ) 2 k N 1 e 0 t λ ( s ) μ ( s ) + μ ( s ) k + λ ( s ) k 2 d s + K L t p * ( 0 ) p * * ( 0 ) .

7. Illustrative Examples

Let us take the model from the previous section as the basis for the two illustrative numerical examples: one with time-independent intensities and the other with time-dependent intensities.
Let λ = 1 , μ = 16 . It is straightforward to check that in such a case, K = 16 , L = 1 and
min k λ μ + μ k + λ k 2 + K L = min k 13 + 16 k + k 2 = 1 ,
which is attained when k = 2 . Then, (33) provides us with the following bound:
p * ( t ) p * * ( t ) 2 N e t p * ( 0 ) p * * ( 0 ) .
Now assume that the Poisson arrival process has the periodic intensity λ ( t ) = 0.5 + 0.5 sin ( π t ) , μ ( t ) = 15 + sin ( π t ) . Then, K = sup t μ ( t ) = 16 and L = sup t λ ( t ) = 1 . Hence, the minimum is again attained at k = 2 , i.e.,
min k > 1 sup t > 0 ( 0.5 + 0.5 sin ( π t ) ) ) ( 15 + sin ( π t ) ) + 15 + sin ( π t ) k + ( 0.5 + 0.5 sin ( π t ) ) k 2 + K L = min k > 1 sup t > 0 2 + sin ( π t ) 1 ,
and from (33), we obtain the bound
p * ( t ) p * * ( t ) 2 N e t p * ( 0 ) p * * ( 0 ) .
By plugging N = 100 into the last inequality, one has that p * ( t ) p * * ( t )   10 8 for t 90 ; i.e., starting from the instant t = 90 , the distribution of the considered CTMC can be regarded as limiting (with the error not greater than 10 8 ). The behaviour of the (two largest) probabilities p 0 ( t ) and p 1 ( t ) under two different initial conditions can be seen in Figure 1 and Figure 2. Since the upper bounds in Theorem 1 are not tight, the value t = 90 is clearly a pessimistic estimate. The system “forgets” its initial state much earlier and already, beyond t > 10 , the influence of the initial conditions vanishes (see Figure 3 and Figure 4). The limiting distribution of X ( t ) is periodic with the period equal to 2 (being the smallest common multiple of the periods of λ ( t ) and μ ( t ) ). The behaviour of the probabilities p 0 ( t ) and p 1 ( t ) on a single period, which covers the transition to the limiting regime, is depicted in Figure 5 and Figure 6.

8. Conclusions

Although the convergence bounds, if obtained using the proposed technique, are, in general, not tight, they can cover quite general CTMCs, which are beyond the scope of other methods. For CTMCs with countable state spaces, the logarithmic norm method proved to be useful in establishing the truncation bounds for numerical solutions of systems of ODEs with an infinite number of equations. Therefore, the generalization of the proposed technique to the case of the countable state space is a promising direction for further research.

Author Contributions

Conceptualization, supervision, A.Z.; methodology, Y.S. and R.R.; software, validation, visualization, I.U.; investigation, writing, I.U., Y.S., R.R. and A.Z. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Ministry of Science and Higher Education of the Russian Federation, project No. 075-15-2024-544.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mitrophanov, A.Y. The Arsenal of Perturbation Bounds for Finite Continuous-Time Markov Chains: A Perspective. Mathematics 2024, 12, 1608. [Google Scholar] [CrossRef]
  2. Giorno, V.; Nobile, A. On a class of birth-death processes with time-varying intensity functions. Appl. Math. Comput. 2020, 379, 125255. [Google Scholar] [CrossRef]
  3. Viswanath, N. Transient study of Markov models with time-dependent transition rates. Oper. Res. 2020, 1, 35. [Google Scholar] [CrossRef]
  4. Schwarz, J.A.; Selinka, G.; Stolletz, R. Performance analysis of time-dependent queueing systems: Survey and classification. Omega 2016, 63, 170–189. [Google Scholar] [CrossRef]
  5. Kwon, S.; Gautam, N. Guaranteeing performance based on time-stability for energy-efficient data centers. IIE Trans. 2016, 48, 812–825. [Google Scholar] [CrossRef]
  6. Ivanova, D.V.; Markova, E.V.; Shorgin, S.Y.; Gaidamaka, Y.V. Priority-based eMBB and URLLC traffic coexistence models in 5G NR industrial deployments. Inform. Appl. 2023, 17, 64–70. [Google Scholar] [CrossRef]
  7. Samoylov, A.K.; Platonova, A.A.; Shorgin, V.S.; Gaidamaka, Y.V. On modeling the effects of multicast traffic servicing in 5G NR networks. Inform. Appl. 2023, 17, 71–77. [Google Scholar] [CrossRef]
  8. Daraseliya, A.V.; Sopin, E.S.; Moltchanov, D.A.; Samouylov, K.E. Analysis of 5G NR base stations offloading by means of NR-U technology. Inform. Appl. 2021, 15, 98–111. [Google Scholar] [CrossRef]
  9. Sopin, E.S.; Maslov, A.R.; Shorgin, V.S.; Begishev, V.O. Modeling insistent user behavior in 5G New Radio networks with rate adaptation and blockage. Inform. Appl. 2023, 17, 25–32. [Google Scholar] [CrossRef]
  10. Rasmi, K.; Jacob, M.J.; Alexander Dudin, A. Krishnamoorthy. Queueing Inventory System with Multiple Service Nodes and Addressed Retrials from a Common Orbit. Methodol. Comput. Appl. Probab. 2024, 26, 1–15. [Google Scholar] [CrossRef]
  11. D’Apice, C.; D’arienzo, M.P.; Dudin, A.; Manzo, R. Optimal Hysteresis Control via a Queuing System with Two Heterogeneous Energy-Consuming Servers. Mathematics 2023, 11, 4515. [Google Scholar] [CrossRef]
  12. Dudin, A.N.; Dudin, S.A.; Klimenok, V.I.; Dudina, O.S. Stability of Queueing Systems with Impatience, Balking and Non-Persistence of Customers. Mathematics 2024, 12, 2214. [Google Scholar] [CrossRef]
  13. Barabanova, E.A.; Vishnevsky, V.M.; Vytovtov, K.A.; Semenova, O.V. Methods of analysis of information-measuring system performance under fault conditions. Phys. Bases Instrum. 2022, 11, 49–59. Available online: https://elibrary.ru/item.asp?id=54046598 (accessed on 20 March 2024). [CrossRef]
  14. Motyer, A.J.; Taylor, P.G. Decay rates for quasi-birth-and-death processes with countably many phases and tridiagonal block generators. Adv. Appl. Probab. 2006, 38, 522. [Google Scholar] [CrossRef]
  15. Wang, W.; Zhu, Q. Stability analysis of semi-Markov switched stochastic systems. Automatica 2018, 94, 72–80. [Google Scholar] [CrossRef]
  16. Earnshaw, B.A.; Keener, J.P. Global asymptotic stability of solutions of nonautonomous master equations. SIAM J. Appl. Dyn. Syst. 2010, 9, 220–237. [Google Scholar] [CrossRef]
  17. Down, D.; Meyn, S.P.; Tweedie, R.L. Exponential and uniform ergodicity of Markov processes. Ann. Probab. 1995, 23, 1671–1691. [Google Scholar] [CrossRef]
  18. Meyn, S.P.; Tweedie, R.L. Markov Chains and Stochastic Stability; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  19. Kartashov, N.V. Strongly stable Markov chains. J. Soviet Math. 1986, 34, 1493–1498. [Google Scholar] [CrossRef]
  20. 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]
  21. Zeifman, A.; Razumchik, R.; Satin, Y.; Kiseleva, K.; Korotysheva, A.; Korolev, V. Bounds on the Rate of Convergence for One Class of Inhomogeneous Markovian Queueing Models with Possible Batch Arrivals and Services. Int. J. Appl. Math. Comput. Sci. 2018, 28, 141–154. [Google Scholar] [CrossRef]
  22. Zeifman, A.; Satin, Y.; Kiseleva, K.; Korolev, V.; Panfilova, T. On limiting characteristics for a non-stationary two-processor heterogeneous system. Appl. Math. Comput. 2019, 351, 48–65. [Google Scholar] [CrossRef]
  23. Usov, I.; Satin, Y.; Zeifman, A.; Korolev, V. Ergodicity Bounds and Limiting Characteristics for a Modified Prendiville Model. Mathematics 2022, 1, 4401. [Google Scholar] [CrossRef]
  24. Kovalev, I.A.; Satin, Y.A.; Sinitcina, A.V.; Zeifman, A.I. On an approach for estimating the rate of convergence for nonstationary Markov models of queueing systems. Inform. Appl. 2022, 16, 75–82. [Google Scholar] [CrossRef]
  25. Zeifman, A.I.; Satin, Y.A.; Kovalev, I.A. On one nonstationary service model with catastrophes and heavy tails. Inform. Appl. 2021, 15, 20–25. [Google Scholar] [CrossRef]
  26. Satin, Y.A. On the bounds of the rate of convergence for Mt/Mt/1 model with two different requests. Syst. Means Inform. 2021, 31, 17–27. [Google Scholar] [CrossRef]
  27. Usov, I.; Satin, Y.; Zeifman, A. On the rate of convergence and limiting characteristics for one quasi-birth-death process. Inform. Appl. 2023, 7, 49–57. [Google Scholar] [CrossRef]
  28. Evans, R.V. Geometric distribution in some two-dimensional queuing systems. Oper. Res. 1967, 15, 830–846. [Google Scholar] [CrossRef]
  29. Wallace, V. The Solution of Quasi Birth and Death Processes Arising from Multiple Access Computer Systems. Ph.D. Thesis, Systems Engineering Laboratory, University of Michigan, Ann Arbor, MI, USA, 1969. [Google Scholar]
  30. Neuts, M.F. Matrix-Geometric Solutions in Stochastic Models: An Algorithmic Approach; Johns Hopkins Press: Baltimore, MD, USA, 1981. [Google Scholar]
  31. Bladt, M.; Neuts, M.F. Matrix-Exponential Distributions: Calculus and Interpretations via Flows. Stoch. Model. 2003, 19, 113–124. [Google Scholar] [CrossRef]
  32. Neuts, M.F. Structured Stochastic Matrices of M/G/1 Type and Their Applications; Marcel Dekker: New York, NY, USA, 1989. [Google Scholar]
  33. He, Q.; Neuts, M.F. Two M/M/1 queues with transfers of customers. Queueing Syst. 2002, 42, 377–400. [Google Scholar] [CrossRef]
  34. Ost, A. Quasi-Birth-and-Death Processes. In Performance of Communication Systems; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar] [CrossRef]
  35. Daleckii, J.L.; Krein, M.G. Stability of solutions of differential equations in Banach space. Am. Math. Soc. 2002, 43, 1024–1102. [Google Scholar]
  36. Van Doorn, E.A. Conditions for exponential ergodicity and bounds for the decay parameter of a birth-death process. Adv. Appl. Probab. 1985, 17, 514–530. [Google Scholar] [CrossRef]
  37. Fadiloglu, M.M.; Yeralan, S. Models of production lines as quasi-birth-death processes. Math. Comput. Model. 2002, 35, 913–930. [Google Scholar] [CrossRef]
Figure 1. The behaviour of the probability p 0 ( t ) under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Figure 1. The behaviour of the probability p 0 ( t ) under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Mathematics 12 02763 g001
Figure 2. The behaviour of the probability p 1 ( t ) under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Figure 2. The behaviour of the probability p 1 ( t ) under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Mathematics 12 02763 g002
Figure 3. The behaviour of the probability p 0 ( t ) on the interval [20, 22] under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Figure 3. The behaviour of the probability p 0 ( t ) on the interval [20, 22] under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Mathematics 12 02763 g003
Figure 4. The behaviour of the probability p 1 ( t ) on the interval [20, 22] under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Figure 4. The behaviour of the probability p 1 ( t ) on the interval [20, 22] under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Mathematics 12 02763 g004
Figure 5. The single period behaviour of the probability p 0 ( t ) under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Figure 5. The single period behaviour of the probability p 0 ( t ) under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Mathematics 12 02763 g005
Figure 6. The single period behaviour of the probability p 1 ( t ) under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Figure 6. The single period behaviour of the probability p 1 ( t ) under two different initial conditions X ( 0 ) = 0 and X ( 0 ) = 99 .
Mathematics 12 02763 g006
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

Satin, Y.; Razumchik, R.; Zeifman, A.; Usov, I. On One Approach to Obtaining Estimates of the Rate of Convergence to the Limiting Regime of Markov Chains. Mathematics 2024, 12, 2763. https://doi.org/10.3390/math12172763

AMA Style

Satin Y, Razumchik R, Zeifman A, Usov I. On One Approach to Obtaining Estimates of the Rate of Convergence to the Limiting Regime of Markov Chains. Mathematics. 2024; 12(17):2763. https://doi.org/10.3390/math12172763

Chicago/Turabian Style

Satin, Yacov, Rostislav Razumchik, Alexander Zeifman, and Ilya Usov. 2024. "On One Approach to Obtaining Estimates of the Rate of Convergence to the Limiting Regime of Markov Chains" Mathematics 12, no. 17: 2763. https://doi.org/10.3390/math12172763

APA Style

Satin, Y., Razumchik, R., Zeifman, A., & Usov, I. (2024). On One Approach to Obtaining Estimates of the Rate of Convergence to the Limiting Regime of Markov Chains. Mathematics, 12(17), 2763. https://doi.org/10.3390/math12172763

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