PROBABILITY-STATISTICAL ANALYSIS OF QUEUING NETWORKS WITH INFINITE NUMBER OF SERVERS
* G.Tsitsiashvili, ** M. Osipova •
IAM, FEB RAS, Vladivostok, Russia e-mails: * guram@ iam.dvo.ru. ** mao1975@list.ru
ABSTRACT
In this paper queuing networks with infinite number of servers and exponentially distributed service times are considered. Limit distributions for numbers of customers in aggregated nodes are calculated. Using these calculations statistical estimates of limit distributions parameters are constructed. Mean times of customer sojourn time in opened network and its analog for closed network are estimated also. Obtained results are spread onto estimates of populations and life cycle parameters.
Keywords: ruin probability, transportation problem, asymptotic formula, enumeration problem
INTRODUCTION
In (Tsitsiashvili & Osipova, 2009) a problem of estimates of product limit distributions parameters is formulated and is solved for opened and closed queuing networks with single server nodes. In this paper opened and closed queuing networks with infinite number of servers in each node are considered. A problem of an aggregation of nodes in these networks is analyzed and limit product distributions are represented via some analogs of nodes load coefficients. Such approach allows to obtain formulas for efficiency indexes of networks and to construct algorithms of their estimates by observations of numbers of customers in different nodes. We speak about a customer mean sojourn time in opened network and about mean time between customer successive appearances in some node for closed network. Obtained results allow analyzing a dynamics of a population with individuals which have many states (Hun & Gulevich, 2008) and life cycle of technical systems with many states (Yung-Wen Liu & Kailash C. Kapur, 2008).
1. OPENED EXPONENTIAL NETWORKS
Consider opened queuing network with the states set S = {0,1,...,m} and irresolvable matrix of
transition probabilities © = =0 and Poisson input flow with intensity X> 0. In each node i e I,
I = {1,..., m}, there is infinite number of servers with service time which has exponential distribution with parameter fj,f > 0 . The node 0 is dummy, it is a source of arriving customers and is a runoff of customers which depart from network. For fixed X > 0 the system of linear algebraic equations (Basharin & Tolmachev, 1983)
(M,-,X ) = (M,..,X, )© (1)
has single solution Xi > 0, i e I. Markov process (n1 (t),...,nm (t)) which describes numbers of
customers in nodes of opened network is ergodic (Foss, 1991) and its limit distribution has the following form (Serfozo, 1999, Example 1.29) (see (Basharin & Tolmachev, 1983) also)
m pni X
P(nl,...,nm ) = n exp(~Pi, Pi = —, i e1 , (2)
i=1 ni! ¡i
Theorem 1. Almost surely we have
Km nj = p, nj , ; e I. (3)
T J
Theorem 1 statement is obtained from the law of large numbers for discrete Markov processes
(Serfozo, 1999, Theorem 1.2).
Construct a model of aggregated network uniting nodes from not intersected subsets I1,..., Ir
of set I and conserve initial order of customers service and routing. Denote (N1 (t),...,Nr (t))
random process (non Markov) describing numbers of customers in united nodes, Nk (t) = £ ni (t),
ieIk
k = 1,...r.
Theorem 2. The process (N1 (t),...,Nr (t)) has limit distribution
r RNk
P(N1,...,Nr)=nexp(-R)Nrr, Rk = £ Pr (4)
k=1 Nk ! ieIk
and almost surely
lim NJ = R,, NJ =JjNj (t)dt, j = 1,...,r . (5)
j j j j j v '
Proof. Formula (4) is a corollary of formula (2) and well known fact that (Feller, 1984, Chapter XI) a sum of independent random variables with Poisson distributions and parameters a, b has Poisson distribution with parameter a + b . Formula (5) is a corollary of formulas (3), (4).
Suppose that f is mean sojourn time of customer in opened network which equals with mean time of customer service in network nodes. The quantity f may be interpreted as a mean life time of individual represented by a customer of input flow. Then we have (Borovkov, 1972, § 31, Theorem 6)
f=X, R=£ Pi. (6)
X ieI
So if we estimate the input intensity X then theorem 2 allows to estimate the parameter R and consequently the parameter f .
2. DYNAMICS OF POPULATION OF INDIVIDUALS WITH MANY STATEDS
In (Hun & Gulevich, 2008) a model of a dynamics of a population with n individuals is considered. A state of k -th individual (a stage of illness) is described by Markov chain xk (t),
t = 0,1,..., with the states set I = {1,..., m} and with the transition probabilities
0iii, 0i1 > 0, i = 1,...,m , 6i i+1 > 0, i = 1,...,m -1.
A transition of an individual into state 1 means that it dies and immediately new individual burns with same number. If the process xk (t) is interpreted as a displacement of particle then it is in i -th
state random time with geometric distribution and with parameter 1 -Qi i. For ergodic process
xk (t) it is possible easily to calculate stationary distribution n(i), i e I.
In (Hun & Gulevich, 2008) if Markov chains xk (t), k = 1,...,n, are independent then numbers
of individuals which are in state i -th satisfy the law of large numbers (Serfozo, 1999, Theorem 1.2) and almost surely
T n'
lim — = K;, i e I. T n
These formulas allowed to analyze optimal strategies of a prevention of non infectious diseases.
In this paper a model of a population dynamics in continuous time is constructed. In constructed model non realistic assumption of constant number of individuals in a population is cancelled.
We consider population model as opened queuing network with infinite number of servers in each node. Then a vector of numbers of individuals in different stages of diseases is described by Markov process (n1 (t),...,nm (t)). In this case sojourn time of individual in i-th state has exponential distribution with the parameter ¡if. A presence of infinite number of servers in each node guarantees independent individuals displacements along network nodes.
Analogously to (Hun & Gulevich, 2008) assume that the route matrix © of opened queuing network contains following non zero elements:
#i,i +1, #i,0 = 1- #i,i+1, i = 1,•••, m - 1 #m,0 = #0,1 = 1 •
Then it is easy to obtain from (1) that
X+1 = Xiдi,i+l, i = m, And consequently from (6) we have that mean individual life time is
f = 2^ e( k, k+1).
i=1k=0 ¡¡i
Remark 1. Suppose that customer-individual passing through i -th node may receive additional service-prevention in parallel node i' with intensity ¡ii < a. That increases mean sojourn time of customer in this node. If customer arrives in node i then with probability ai it is served on server of this node and with probability bt = 1 - a^ it is served in node i'. As result equality p^ = /¡if is replaced by equality
aiXi bX
Pi = —, Pr = —.
¡i a
Using these equalities it is possible to estimate influence of prevention on main indexes of this network.
3. CLOSED EXPONENTIAL NETWORKS
Consider closed queuing network with the states set S, the route matrix ©, n customers circulating along this network and n servers in each node with service intensity a on a server of i -th node, i e S. Suppose that for fixed Xq > 0 the system of linear algebraic equations (1) has single solution X1 > 0,...,Xm > 0. Then Markov process (n0 (t),...,nm (t)) characterizing numbers of
customers in nodes of closed network is ergodic (Ivchenko & Kashtanov & Kovalenko, 1982, Theorem 2.4) and has limit multinomial distribution (Serfozo, 1999, Example 1.29)
P(no,...,nm) = n!n^-r, di = Pi , Po = X0, ieS, 2 n = n. (7)
i=0 ni! Po + - + Pm ¡0 i=0
Theorem 3. Almost surely
nT
lim = dt, i e S. (8)
tn
Proof. From binomial theorem we have that the process
(no (t), «1 (t) m- 2 (t), nm -1 + nm (t))
has limit distribution
(d + d \ m-\ m m-2 d i
P("0,n1,...nm-2,nm-1 + nm ) = n!( m> )! n TTT. (9)
(nm-1 + nm )! i=0 ni!
Using mathematical induction it is easy to prove that the process I n0 (t), £ nk (t) I has limit
V keI J
distribution
P(n n n ) = n!d0n° (d1 +... + dm Tn°
p (n0, nl,..., nm )= n0!( n - n0)! • Consequently the process n0 (t) has binomial limit distribution with mean nd0 . Analogously it is possible to prove that random process ni (t) has binomial limit distribution with mean ndi, i = 1,...,m. Then from the law of large numbers for ergodic Markov processes (Serfozo, 1999, Theorem 1.2) we obtain the formula (8).
Construct a model of aggregated network uniting nodes in not intersected subsets I1,...,Ir of the set I , conserving initial characteristics of customers service and routing. Denote (N1 (t),...,Nr (t)) random process (non Markov) which describes numbers of customers in united
nodes, Nk (t)= £ ni (t), k = 1,...,r .
ieIk
Theorem 4. The process (N1 (t),...,Nr (t)) has multinomial limit distribution
rDNk
! n Nr, Dk = £
k =1 Nk ! ieI,
DNk
P(N1,...Nr) = n! n ^T, Dk = £ di, k = 1,...,r, (10)
and almost surely
lim Nt = Dk, k = 1,...,r. (11)
Tn
Proof. The formula (10) may be obtained from the formula (9) by mathematical induction. The formula (11) is a corollary of the formula (10) and the theorem 3.
Denote fi mean time between successive appearances of a customer in i -th node of closed network. To calculate f0 consider opened network with route matrix 0 and with input intensity X0. It is obvious that
f = 1 + f = 1 + R = p0 +... + pm = p0 = 1 (12)
¡0 ¡0 X0 X0 d0X0 ¡0d0 Consequently if we know statistical estimate of service intensity u0 then using the theorem 3 it is
possible to estimate d0 and consequently f0. Analogously it is possible to obtain the equality
f = -d, (13)
Uidi
and using it to estimate mean fi by estimate of service intensity u , i e I.
CONCLUSION
Theorems 2, 4 allow to aggregate nodes of initial network and consequently to simplify a procedure of statistical estimates of its parameters. It means that we may delete a solution of the system (1) and to simplify a procedure of the network observation. Besides of product distribution parameters suggested approach allows to estimate mean sojourn times of customers in the network. Obtained results may be generalized onto queuing networks with non exponential service time distributions using (Harrison & Lemoine, 1981) and may be applied in technical gerontology and demography to estimate mean time of endowment for persons of fixed age.
REFERENCES
1. Tsitsiashvili G. Sh., Osipova M.A. 2009. Parameter Estimation for Product-Form Distributions of Queueing Networks. Problems of Information Transmission 45 (4). 400-405.
2. Yung-Wen-Liu, Kailash C. Kapur. 2008. New Models and Measures for Reliability of Multi-state Systems. In Handbook of Performability Engineering. Krishna B. Misra (ed). New York: SpringerVerlag. Pp. 431-445.
3. Basharin G.P., Tolmachev A.L. 1983. The theory of queuing networks and its applications to the analysis of information and computation systems. The totals of science and engineering, the probability theory. M. VINITI. Pp. 3-119. (In Russian).
4. Foss S.G. 1991. Ergodicity of queueing networks. Siberian Mathematical Journal 32 (4): 184203. (In Russian).
5. Serfozo R. 1999. Introduction to Stochastic Networks. New York: Springer Verlag.
6. Feller W. 1984. Introduction to probability theory and its applications. M.: World, Part 1 (in Russian).
7. Borovkov A.A. 1972. Probability processes in queueing theory. M.: Science (In Russian).
8. Hun I.S., Gulevich A.P. 2008. Optimization of prevention strategies for not infection diseases. Bulletin of Far Eastern State Economical University (2): 10-22. (In Russian).
9. Ivchenko G.I., Kashtanov V.A., Kovalenko I.N. 1982. Queueing theory. M.:High school (in Russian).
10. Harrison J.M., Lemoine A.J. 1981. A note on networks of infinite-server queues. J. Appl. Probab. (2): 561-567.