A Simple Solution for the M/D/c Waiting Time Distribution
|
|
- Aubrey Stanley
- 6 years ago
- Views:
Transcription
1 A Simple Solution for the M/D/c Waiting Time Distribution G.J.Franx, Universiteit van Amsterdam November 6, 998 Abstract A surprisingly simple and explicit expression for the waiting time distribution of the M/D/c queueing system is derived by a full probabilistic analysis, requiring neither generating functions nor Laplace transforms. Unlike the solutions known so far, this expression presents no numerical complications, not even for high traffic intensities. Finally, the result is proved explicitly to satisfy Erlang s integral equation for the M/D/c queue, which has been somewhat problematic for the expressions known so far. M/D/c waiting time distribution, numerically convenient explicit formula, full probabilistic analysis, Erlang s M/D/c integral equation. Introduction The M/D/c system is one of the classical queueing models that have already been studied since the beginning of this century. In this model we assume the customers to arrive according to a Poisson process with rate λ. There are c identical servers, serving each customer on a first come first serve basis during a constant time D. We assume the traffic intensity ρ λd/c <. Erlang [909] derived an explicit expression for the waiting time distribution F (y) of the M/D/ queue, by solving for c the following integral equation: F (y) 0 F (x + y D) λc x c (c )! e λx dx, y 0. () This equation is obtained by comparing the waiting time of an arbitrary customer A to the waiting time of the c th customer after A, who will arrive between x and x + dx time units after A with probability λc x c (c )! e λx dx. If this customer has to wait, he will be served exactly D time units after A. Since the M/D/c queue is equivalent to the E c /D/ queue with respect to the waiting time distribution, Erlang s equation can be interpreted as a special case of Lindley s integral
2 equation for the GI/G/ queueing system. Although his equation proved quite successful for the case c, Erlang [920] realized that for c > it would hardly lead to an explicit mathematical solution. By an ingenious probabilistic argument and the use of generating functions Crommelin [932] derived a general expression for the waiting time distribution of the M/D/c queue for all c N, which for c corresponds to Erlang s result. If P n denotes the stationary probability of the system containing no more than n customers, Crommelin s result reads: c P {W x} n0 P n m { λ(x md)} (k+)c n e λ(x md), md x < (m + )D (2) {(k + )c n}! Prabhu [962] proved that for c Erlang s integral equation yields a solution of type (2), where P n is replaced by some alternatively defined constant α n. However, it is unresolved how to interpret α n as the cumulative state probability. Apart from this, Crommelin s result is not really practical for numerical purposes, due to alternating terms which are in general much larger than their sum. As a way to get around the problem of round off errors, a recursion scheme based on Crommelin s argument is described in Tijms [994]. However, for increasing c and ρ this recursion scheme (too) will ultimately be hampered by round off errors, for which case an asymptotic expansion is recommended. This paper presents an alternative probabilistic approach, leading to a simple formula for the waiting time distribution, which is numerically stable for all ρ <. This formula (0 or ) is the main result of this paper, and is derived without the use of generating functions or Laplace transforms. Before doing so, some preliminary results about the number of customers in the system are introduced, and an essential lemma about the queue length at service initiations is proved. Obviously, the derived expression (0) must satisfy Erlang s integral equation (. In order to complete the circle, an explicit proof of this is given in section 5. 2 The number of customers present in the system Let p i (t) denote the probability of the system holding i customers at time t. Observe that all customers in service at time t will have left the system at time t + D. Consequently, customers present at time t + D either arrived during the time interval (t, t + D] or were already waiting for service at time t. Hence, by conditioning on the number of customers present at time t we find: p i (t + D) c p j (t) (λd)i e λd + i+c jc+ p j (t) (λd)i+c j (i + c j)! e λd, t R, i N 0. (3) 2
3 The stationary distribution p i lim t p i (t) is found by letting t in (3): p i c (λd) i p j e λd + i+c jc+ p j (λd) i+c j (i + c j)! e λd, for all i N 0. (4) Analogously, by conditioning on the queue length, we can derive a similar set of equations for q i, the stationary probability of the queue containing i customers: q 0 c c j q j m0 (λd) m i+c e λd, and q i m! q j (λd) i+c j (i + c j)! e λd for i > 0. (5) We can also derive (5) from (4) by using q 0 c p j, and q i p i+c for i > 0. Together with the normalization equation and the capacity utilization equation, (4) constitutes an infinite system of linear equations that can be solved in several ways. In Tijms [994] a Fast Fourier Transform method and a geometric tail approach to solving (4) are described, both of which are quite efficient. 3 The queue length at service initiations For our analysis we will define the following random variables: A n : the arrival time of the n th customer after t 0 S n W n : the starting time of the service of the n th customer : S n A n : the waiting time of the n th customer L + q (t) : the queue length at time t + (immediately after time t) q + i (t) : P {L+ q (t)i} : probability of finding a queue of length i at time t + ω i q i : lim q + i (S n n) lim P {L + q (S n )i}. n : lim P {L q (t)i} : stationary probability of finding a queue of length i t Since the probabilistic argument will concentrate on the service initiations S n, we are interested in the stationary probabilities ω i of having i waiting customers in the queue immediately after some arbitrary service start. Therefore we will prove the following lemma: Lemma. ω i lim n P {L + q (S n )i} lim t P {L q (t)i} q i, for all i N 0. Proof. With probability q + i (S n) there are i customers left in the queue immediately after the service start S n of the n th customer. During the service interval (S n, S n + D] there will 3
4 be j new arrivals with probability (λd)j e λd. We distinguish two cases. If i + j c, the server will start serving the (n + c) th customer immediately after finishing the service of the n th customer. So for this case we find L + q (S n+c ) L + q (S n + D) i + j c. On the other hand, if i + j < c, the (n + c) th customer has not yet arrived at time S n + D. So immediately after finishing the service of the n th customer, the server in question will be idle. Consequently, on arrival the (n + c) th customer will find at least one server idle, implying L + q (S n+c ) 0. Combining both cases, we conclude that L + q (S n+c ) max(0, L + q (S n )+j c), irrespective of the residual service times at t S n. Therefore, the probability q + i (S n+c) of having i waiting customers in the queue immediately after the service start of the (n+c) th customer can be determined by conditioning on L + q (S n ) : q + 0(S n+c ) c j c q + j (S (λd) m n) e λd, and q + i m! (S i+c n+c) q + j (S n) (λd)i+c j (i+c j)! e λd for i>0 (6) m0 The stationary distribution ω i lim j q + i (S j) is found by letting n in (6): ω 0 c c j ω j m0 (λd) m i+c e λd, and ω i m! ω j (λd) i+c j (i + c j)! e λd for i > 0. (7) This set of equations is identical to the set of equations (5) for the stationary probabilities q i. Since the corresponding embedded Markov chain is ergodic for ρ <, there is a unique solution to (7). Therefore ω i equals q i, for all i 0. Remark Observe that (6) is derived by conditioning on merely L + q (S n ). This is possible because of the crucial property that the residual service times of the busy servers at t S n have no influence on L + q (S n+c ). This is not a trivial property. For instance, it does not hold for L, the total number of customers in the system. Therefore lemma holds only for the queue length L q, not for L. 4 The Waiting Time Distribution Without loss of generality, we assume that all customers are assigned to the servers in cyclic order. This will not violate the FCFS discipline because of the deterministic service times. Our first objective will be to determine W n, the waiting time of the n th customer to arrive after t 0. This customer will be called the marked customer and will be served by the so called marked server. 4
5 For (k )D x < kd the analysis of P {W n < x} will concentrate on the time instant S n kc, when the marked server will start to serve the (n kc) th customer. Observe that the marked customer will be the k th customer to be served by the marked server from S n kc onwards. Let u be some positive time lapse D. If the marked customer arrives no earlier than t S n kc + u, the marked server has been serving the (n kc) th customer for at least u time units on the arrival instant A n. Therefore the marked customer will find the marked server with an amount of unfinished work kd u, implying W n kd u. (Here it is essential that u D, guaranteeing that the marked server is busy during all of the time interval (S n kc, S n kc + u]. ) On the other hand, if the marked customer arrives before t S n kc + u, he will find the marked server with an amount of unfinished work > kd u, implying W n > kd u. Combining both cases, we conclude: ( k N) ( u (0, D]) A n S n kc + u W n kd u The corresponding probabilities can be found by conditioning on L + q (S n kc ), the queue length immediately after S n kc. If at this instant the queue contains kc or more customers, the marked customer has already entered the queue and is waiting there in (kc) th position, implying A n < S n kc. On the other hand, if the queue contains i < kc customers at time S n kc, the marked customer will be the (kc i) th customer to arrive from S n kc onwards. With probability kc i Thus, conditioning on L + q (S n kc ) yields: (λu) j e λu this arrival takes place no earlier than t S n kc + u. P {W n kd u} P {A n S n kc + u} kc kc i q + i (S n kc) (λu) j e λu. (8) By letting n, using lemma, we find the stationary waiting time distribution P {W kd u} kc kc i q i (λu) j kc e λu e λu (λu) j kc j q i. (9) Defining the cumulative probability Q m : m 0 q i, we can summarize the result as: kc ( k N) ( u (0, D]) P {W kd u} e λu Q kc j (λu) j Substitution of x kd u gives us the waiting time distribution in terms of x :. (0) kc P {W x} e λ(kd x) Q kc j λ j (kd x) j, for (k )D x < kd. () As this expression contains only a finite number of positive terms, it does not present any numerical complications, regardless of the traffic intensity ρ. 5
6 5 Erlang s integral equation for the M/D/c queue Since it represents the waiting time distribution of the M/D/c queue, expression (0) must satisfy Erlang s integral equation. In order to complete the circle, this section presents an explicit proof, which will make use of the following two lemmas: Lemma 2. n+c (λd) i Q n Q n+c i e λd, for all n N 0. Proof. We define Q n (t) as the probability of the queue containing no more than n customers at time t. Note that Q n lim t Q n (t). Let (t, t + D] be an arbitrary time interval of length D, during which i new arrivals will take place with probability (λd)i e λd. By conditioning on i we find: n+c Q n (t + D) Q n+c i (t) (λd)i e λd, for all n N 0. (2) By letting t in (2) we find lemma 2. Lemma 3. k (j + i + )(k i)! (j + k + )!, for all j, k N 0. Proof. For j 0 lemma 3 holds for all k N 0, which can be verified as follows: k (i + )!(k i)! k+ ( ) l (k + l)! (k + )! k+ l (k + )! k+ ( ) l (k + )! k + ( ) l (k + l)! ( ) l k+ l (k + )!. For general j the validity of the lemma is proved by induction. Suppose the lemma holds for some j and all k N 0, then (j + k + )! k (k i)!(j + i + ) k i (k i)!(j + i + ) + k!(j+)!. (3) 6
7 Since the lemma has already been proved for j 0, we can rewrite the last term for any positive k as (j + )! k! (j + )! k { k n k (i + )!(k i)! k (i+)!(k i)!(j+i+2) + (i + )!(k i)!(j + ) } (j+i+2)(k i)!(j+) ( ) n n!(k n)!(j + n + ) + k (j+)! (j + i + 2)(k i)!. If this expression for k!(j+)! is substituted in the right hand side of (3) we find: k (j + )! (j + i + 2)(k i)! (j + k + )!. Thus we proved the lemma to be valid for j +, k, if it is valid for j, k N. Since there were no restrictions on the value of k in our starting supposition, this finishes our proof. See also exercise 5, p 29, Riordan [968]. In the formulation () of Erlang s integral equation it is required that F (t) 0 for all t < 0. We can get rid of this requirement by reformulating the equation as: F (y) max(0,d y) F (x + y D) λc x c (c )! e λx dx, for all y R +. (4) Because of the special structure of our result (), we will substitute: y md u, with m y D +, and 0 < u D. Abbreviating the integrand as H(x, y), we can rewrite the right hand side of (4): u max(0,u (m )D) H(x, md u)dx + u+(k+)d u+kd This format enables us to substitute expression () in the integrand: + u max(0,u (m )D) u+(k+)d (k+m)c u+kd (m )c λ j+c (u x) j x c Q (m )c j e λu dx (c )! H(x, md u)dx. (5) λ j+c {(k+)d+u x} j x c Q (k+m)c j e λ{(k+)d+u} dx. (6) (c )! 7
8 By carrying out the integration, and realizing that max(0, u (m )D) 0 except for m, in which case the first integral vanishes, we can write (6) as: + (m )c e λu e λ{(k+)d+u} (k+m)c Q (m )c j (λu) j+c λ j+c Q (k+m)c j c c (j + i + )(c i)! {(k+)d+u} c i D j+i+ (j + i + )(c i)!. (7) By expanding {(k + )D + u} c i and subsequently rearranging the summations, using c c i n0 A i,n c n0 A i,n, the second term of (7) can be rewritten as: with β n (k+m)c c e λu n0 (λd) j+c n Q (k+m)c j (λu) n β n, (8) n! e λd (k + ) c i n (j+i+)(c i n)!. (9) Using lemma 3, the last summation in (9) can be transformed as follows: (k + ) c i n (j + i + )(c i n)! i (j + i + ) k l (c n i l)! k l l (j + i + )(c n i l)! k l (j + c n l)!. Substitution of this result in (9) in combination with lemma 2, yields: β n (k+m)c Q (k+m)c j (k+m+)c n l (λd) j+c n l (j + c n l)! e λd (λd) i Q (k+m+)c n l i e λd ic n l { } c n l (λd) i Q (k+m)c n l Q (k+m+)c n l i e λd 8
9 e λ(k+)d Q (k+m)c n l jl Q (k+m+)c n j (λd) j l (j l)!. We can simplify the last term by observing that: jl Q (k+m+)c n j Q (k+m+)c n j (λd) j l (j l)! j (λd) j l (j l)! {λ(k + )D} j Q (k+m+)c n j, which enables us to rewrite β n as: β n 0 e λ(k+)d Q m. jo Q (k+m)c n l {λ(k + )D} j Q (k+m)c n l Q (k+m+)c n j Putting the parts together, using lemma 3, we can reconstruct (7) as follows: (m )c e λu { mc e λu ic mc e λu (λu) i Q mc i (λu) j+c Q (m )c j (j + c)! + c (λu) n e λu n! n0 } (λu) i c (λu) n Q mc i + Q m n! n0 This implies that the integral equation (4) reduces to. β n mc F (y) e λu Q mc i (λu) i, 9
10 which, for y md u, corresponds exactly to expression (0). Acknowledgements The author wishes to thank Onno Boxma, Nico van Dijk, Flora Spieksma, and Henk Tijms for their fruitful remarks and suggestions. References [] Erlang, A.K. (909, 920), The theory of probabilities and telephone conversations, and Telephone waiting times, first published in Nyt Tidsskrift for Matematik, B Vol. 20, p. 33, and Matematisk Tidsskrift, B Vol. 3, p 25. English translations of both articles in Brockmeyer, E. et. al. (948), The Life and Works of A.K.Erlang, The Copenhagen Telephone Company, Copenhagen. [2] Crommelin, C.D. (932), Delay probability formulas when the holding times are constant, P.O. Electr. Engr. J., 25, [3] Prabhu, N.U. (962), Elementary methods for some waiting time problems, Opns. Res. 0, [4] Prabhu, N.U. (965), Queues and Inventories, Wiley, New York. [5] Riordan, J. (968), Combinatorial Identities, Wiley, NewYork. [6] Tijms, H.C. (994), Stochastic Models, an Algorithmic Approach, Wiley, New York. 0
Non Markovian Queues (contd.)
MODULE 7: RENEWAL PROCESSES 29 Lecture 5 Non Markovian Queues (contd) For the case where the service time is constant, V ar(b) = 0, then the P-K formula for M/D/ queue reduces to L s = ρ + ρ 2 2( ρ) where
More informationCHAPTER 4. Networks of queues. 1. Open networks Suppose that we have a network of queues as given in Figure 4.1. Arrivals
CHAPTER 4 Networks of queues. Open networks Suppose that we have a network of queues as given in Figure 4.. Arrivals Figure 4.. An open network can occur from outside of the network to any subset of nodes.
More informationPerformance Evaluation of Queuing Systems
Performance Evaluation of Queuing Systems Introduction to Queuing Systems System Performance Measures & Little s Law Equilibrium Solution of Birth-Death Processes Analysis of Single-Station Queuing Systems
More informationWaiting time characteristics in cyclic queues
Waiting time characteristics in cyclic queues Sanne R. Smits, Ivo Adan and Ton G. de Kok April 16, 2003 Abstract In this paper we study a single-server queue with FIFO service and cyclic interarrival and
More informationOther properties of M M 1
Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected
More informationQueues and Queueing Networks
Queues and Queueing Networks Sanjay K. Bose Dept. of EEE, IITG Copyright 2015, Sanjay K. Bose 1 Introduction to Queueing Models and Queueing Analysis Copyright 2015, Sanjay K. Bose 2 Model of a Queue Arrivals
More informationStatistics 150: Spring 2007
Statistics 150: Spring 2007 April 23, 2008 0-1 1 Limiting Probabilities If the discrete-time Markov chain with transition probabilities p ij is irreducible and positive recurrent; then the limiting probabilities
More informationON THE LAW OF THE i TH WAITING TIME INABUSYPERIODOFG/M/c QUEUES
Probability in the Engineering and Informational Sciences, 22, 2008, 75 80. Printed in the U.S.A. DOI: 10.1017/S0269964808000053 ON THE LAW OF THE i TH WAITING TIME INABUSYPERIODOFG/M/c QUEUES OPHER BARON
More informationThe Transition Probability Function P ij (t)
The Transition Probability Function P ij (t) Consider a continuous time Markov chain {X(t), t 0}. We are interested in the probability that in t time units the process will be in state j, given that it
More informationLIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974
LIMITS FOR QUEUES AS THE WAITING ROOM GROWS by Daniel P. Heyman Ward Whitt Bell Communications Research AT&T Bell Laboratories Red Bank, NJ 07701 Murray Hill, NJ 07974 May 11, 1988 ABSTRACT We study the
More informationIntroduction to Markov Chains, Queuing Theory, and Network Performance
Introduction to Markov Chains, Queuing Theory, and Network Performance Marceau Coupechoux Telecom ParisTech, departement Informatique et Réseaux marceau.coupechoux@telecom-paristech.fr IT.2403 Modélisation
More informationIntroduction to Queuing Networks Solutions to Problem Sheet 3
Introduction to Queuing Networks Solutions to Problem Sheet 3 1. (a) The state space is the whole numbers {, 1, 2,...}. The transition rates are q i,i+1 λ for all i and q i, for all i 1 since, when a bus
More informationContents Preface The Exponential Distribution and the Poisson Process Introduction to Renewal Theory
Contents Preface... v 1 The Exponential Distribution and the Poisson Process... 1 1.1 Introduction... 1 1.2 The Density, the Distribution, the Tail, and the Hazard Functions... 2 1.2.1 The Hazard Function
More information2905 Queueing Theory and Simulation PART III: HIGHER DIMENSIONAL AND NON-MARKOVIAN QUEUES
295 Queueing Theory and Simulation PART III: HIGHER DIMENSIONAL AND NON-MARKOVIAN QUEUES 16 Queueing Systems with Two Types of Customers In this section, we discuss queueing systems with two types of customers.
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationLink Models for Circuit Switching
Link Models for Circuit Switching The basis of traffic engineering for telecommunication networks is the Erlang loss function. It basically allows us to determine the amount of telephone traffic that can
More informationPart I Stochastic variables and Markov chains
Part I Stochastic variables and Markov chains Random variables describe the behaviour of a phenomenon independent of any specific sample space Distribution function (cdf, cumulative distribution function)
More informationClass 11 Non-Parametric Models of a Service System; GI/GI/1, GI/GI/n: Exact & Approximate Analysis.
Service Engineering Class 11 Non-Parametric Models of a Service System; GI/GI/1, GI/GI/n: Exact & Approximate Analysis. G/G/1 Queue: Virtual Waiting Time (Unfinished Work). GI/GI/1: Lindley s Equations
More informationGI/M/1 and GI/M/m queuing systems
GI/M/1 and GI/M/m queuing systems Dmitri A. Moltchanov moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/tlt-2716/ OUTLINE: GI/M/1 queuing system; Methods of analysis; Imbedded Markov chain approach; Waiting
More informationStationary Probabilities of Markov Chains with Upper Hessenberg Transition Matrices
Stationary Probabilities of Marov Chains with Upper Hessenberg Transition Matrices Y. Quennel ZHAO Department of Mathematics and Statistics University of Winnipeg Winnipeg, Manitoba CANADA R3B 2E9 Susan
More informationData analysis and stochastic modeling
Data analysis and stochastic modeling Lecture 7 An introduction to queueing theory Guillaume Gravier guillaume.gravier@irisa.fr with a lot of help from Paul Jensen s course http://www.me.utexas.edu/ jensen/ormm/instruction/powerpoint/or_models_09/14_queuing.ppt
More informationA Queueing System with Queue Length Dependent Service Times, with Applications to Cell Discarding in ATM Networks
A Queueing System with Queue Length Dependent Service Times, with Applications to Cell Discarding in ATM Networks by Doo Il Choi, Charles Knessl and Charles Tier University of Illinois at Chicago 85 South
More informationProbability Models in Electrical and Computer Engineering Mathematical models as tools in analysis and design Deterministic models Probability models
Probability Models in Electrical and Computer Engineering Mathematical models as tools in analysis and design Deterministic models Probability models Statistical regularity Properties of relative frequency
More informationQueueing Theory I Summary! Little s Law! Queueing System Notation! Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K "
Queueing Theory I Summary Little s Law Queueing System Notation Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K " Little s Law a(t): the process that counts the number of arrivals
More informationRecap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks
Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution
More informationChapter 1. Introduction. 1.1 Stochastic process
Chapter 1 Introduction Process is a phenomenon that takes place in time. In many practical situations, the result of a process at any time may not be certain. Such a process is called a stochastic process.
More informationM/G/1 and M/G/1/K systems
M/G/1 and M/G/1/K systems Dmitri A. Moltchanov dmitri.moltchanov@tut.fi http://www.cs.tut.fi/kurssit/elt-53606/ OUTLINE: Description of M/G/1 system; Methods of analysis; Residual life approach; Imbedded
More informationThe discrete-time Geom/G/1 queue with multiple adaptive vacations and. setup/closedown times
ISSN 1750-9653, England, UK International Journal of Management Science and Engineering Management Vol. 2 (2007) No. 4, pp. 289-296 The discrete-time Geom/G/1 queue with multiple adaptive vacations and
More informationChapter 5. Continuous-Time Markov Chains. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan
Chapter 5. Continuous-Time Markov Chains Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan Continuous-Time Markov Chains Consider a continuous-time stochastic process
More informationHITTING TIME IN AN ERLANG LOSS SYSTEM
Probability in the Engineering and Informational Sciences, 16, 2002, 167 184+ Printed in the U+S+A+ HITTING TIME IN AN ERLANG LOSS SYSTEM SHELDON M. ROSS Department of Industrial Engineering and Operations
More informationRelating Polling Models with Zero and Nonzero Switchover Times
Relating Polling Models with Zero and Nonzero Switchover Times Mandyam M. Srinivasan Management Science Program College of Business Administration The University of Tennessee Knoxville, TN 37996-0562 Shun-Chen
More information1 Basic concepts from probability theory
Basic concepts from probability theory This chapter is devoted to some basic concepts from probability theory.. Random variable Random variables are denoted by capitals, X, Y, etc. The expected value or
More informationAsymptotics for Polling Models with Limited Service Policies
Asymptotics for Polling Models with Limited Service Policies Woojin Chang School of Industrial and Systems Engineering Georgia Institute of Technology Atlanta, GA 30332-0205 USA Douglas G. Down Department
More informationDerivation of Formulas by Queueing Theory
Appendices Spectrum Requirement Planning in Wireless Communications: Model and Methodology for IMT-Advanced E dite d by H. Takagi and B. H. Walke 2008 J ohn Wiley & Sons, L td. ISBN: 978-0-470-98647-9
More informationA tandem queueing model with coupled processors
A tandem queueing model with coupled processors Jacques Resing Department of Mathematics and Computer Science Eindhoven University of Technology P.O. Box 513 5600 MB Eindhoven The Netherlands Lerzan Örmeci
More informationLecture 20: Reversible Processes and Queues
Lecture 20: Reversible Processes and Queues 1 Examples of reversible processes 11 Birth-death processes We define two non-negative sequences birth and death rates denoted by {λ n : n N 0 } and {µ n : n
More information15 Closed production networks
5 Closed production networks In the previous chapter we developed and analyzed stochastic models for production networks with a free inflow of jobs. In this chapter we will study production networks for
More informationContinuous Time Processes
page 102 Chapter 7 Continuous Time Processes 7.1 Introduction In a continuous time stochastic process (with discrete state space), a change of state can occur at any time instant. The associated point
More informationAnswers to selected exercises
Answers to selected exercises A First Course in Stochastic Models, Henk C. Tijms 1.1 ( ) 1.2 (a) Let waiting time if passengers already arrived,. Then,, (b) { (c) Long-run fraction for is (d) Let waiting
More informationBIRTH DEATH PROCESSES AND QUEUEING SYSTEMS
BIRTH DEATH PROCESSES AND QUEUEING SYSTEMS Andrea Bobbio Anno Accademico 999-2000 Queueing Systems 2 Notation for Queueing Systems /λ mean time between arrivals S = /µ ρ = λ/µ N mean service time traffic
More informationJ. MEDHI STOCHASTIC MODELS IN QUEUEING THEORY
J. MEDHI STOCHASTIC MODELS IN QUEUEING THEORY SECOND EDITION ACADEMIC PRESS An imprint of Elsevier Science Amsterdam Boston London New York Oxford Paris San Diego San Francisco Singapore Sydney Tokyo Contents
More informationEconomy of Scale in Multiserver Service Systems: A Retrospective. Ward Whitt. IEOR Department. Columbia University
Economy of Scale in Multiserver Service Systems: A Retrospective Ward Whitt IEOR Department Columbia University Ancient Relics A. K. Erlang (1924) On the rational determination of the number of circuits.
More information11 The M/G/1 system with priorities
11 The M/G/1 system with priorities In this chapter we analyse queueing models with different types of customers, where one or more types of customers have priority over other types. More precisely we
More informationOn the static assignment to parallel servers
On the static assignment to parallel servers Ger Koole Vrije Universiteit Faculty of Mathematics and Computer Science De Boelelaan 1081a, 1081 HV Amsterdam The Netherlands Email: koole@cs.vu.nl, Url: www.cs.vu.nl/
More information15 Closed production networks
5 Closed production networks In the previous chapter we developed and analyzed stochastic models for production networks with a free inflow of jobs. In this chapter we will study production networks for
More informationMarkov Chains. X(t) is a Markov Process if, for arbitrary times t 1 < t 2 <... < t k < t k+1. If X(t) is discrete-valued. If X(t) is continuous-valued
Markov Chains X(t) is a Markov Process if, for arbitrary times t 1 < t 2
More informationKendall notation. PASTA theorem Basics of M/M/1 queue
Elementary queueing system Kendall notation Little s Law PASTA theorem Basics of M/M/1 queue 1 History of queueing theory An old research area Started in 1909, by Agner Erlang (to model the Copenhagen
More informationQueueing Theory II. Summary. ! M/M/1 Output process. ! Networks of Queue! Method of Stages. ! General Distributions
Queueing Theory II Summary! M/M/1 Output process! Networks of Queue! Method of Stages " Erlang Distribution " Hyperexponential Distribution! General Distributions " Embedded Markov Chains M/M/1 Output
More informationQueueing Review. Christos Alexopoulos and Dave Goldsman 10/6/16. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA
1 / 24 Queueing Review (mostly from BCNN) Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/6/16 2 / 24 Outline 1 Introduction 2 Queueing Notation 3 Transient
More informationAnalysis of an M/G/1 queue with customer impatience and an adaptive arrival process
Analysis of an M/G/1 queue with customer impatience and an adaptive arrival process O.J. Boxma 1, O. Kella 2, D. Perry 3, and B.J. Prabhu 1,4 1 EURANDOM and Department of Mathematics & Computer Science,
More informationIEOR 6711, HMWK 5, Professor Sigman
IEOR 6711, HMWK 5, Professor Sigman 1. Semi-Markov processes: Consider an irreducible positive recurrent discrete-time Markov chain {X n } with transition matrix P (P i,j ), i, j S, and finite state space.
More informationOne important issue in the study of queueing systems is to characterize departure processes. Study on departure processes was rst initiated by Burke (
The Departure Process of the GI/G/ Queue and Its MacLaurin Series Jian-Qiang Hu Department of Manufacturing Engineering Boston University 5 St. Mary's Street Brookline, MA 2446 Email: hqiang@bu.edu June
More informationNetworks of Queues Models with Several. Classes of Customers and Exponential. Service Times
Applied Mathematical Sciences, Vol. 9, 2015, no. 76, 3789-3796 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.53287 Networks of Queues Models with Several Classes of Customers and Exponential
More informationAn engineering approximation for the mean waiting time in the M/H 2 b /s queue
An engineering approximation for the mean waiting time in the M/H b /s queue Francisco Barceló Universidad Politécnica de Catalunya c/ Jordi Girona, -3, Barcelona 08034 Email : barcelo@entel.upc.es Abstract
More informationP (L d k = n). P (L(t) = n),
4 M/G/1 queue In the M/G/1 queue customers arrive according to a Poisson process with rate λ and they are treated in order of arrival The service times are independent and identically distributed with
More informationTHE INTERCHANGEABILITY OF./M/1 QUEUES IN SERIES. 1. Introduction
THE INTERCHANGEABILITY OF./M/1 QUEUES IN SERIES J. Appl. Prob. 16, 690-695 (1979) Printed in Israel? Applied Probability Trust 1979 RICHARD R. WEBER,* University of Cambridge Abstract A series of queues
More informationPart II: continuous time Markov chain (CTMC)
Part II: continuous time Markov chain (CTMC) Continuous time discrete state Markov process Definition (Markovian property) X(t) is a CTMC, if for any n and any sequence t 1
More informationIntroduction to Queueing Theory
Introduction to Queueing Theory Raj Jain Washington University in Saint Louis Jain@eecs.berkeley.edu or Jain@wustl.edu A Mini-Course offered at UC Berkeley, Sept-Oct 2012 These slides and audio/video recordings
More informationIEOR 6711: Stochastic Models I, Fall 2003, Professor Whitt. Solutions to Final Exam: Thursday, December 18.
IEOR 6711: Stochastic Models I, Fall 23, Professor Whitt Solutions to Final Exam: Thursday, December 18. Below are six questions with several parts. Do as much as you can. Show your work. 1. Two-Pump Gas
More informationMultiserver Queueing Model subject to Single Exponential Vacation
Journal of Physics: Conference Series PAPER OPEN ACCESS Multiserver Queueing Model subject to Single Exponential Vacation To cite this article: K V Vijayashree B Janani 2018 J. Phys.: Conf. Ser. 1000 012129
More informationMarkov processes and queueing networks
Inria September 22, 2015 Outline Poisson processes Markov jump processes Some queueing networks The Poisson distribution (Siméon-Denis Poisson, 1781-1840) { } e λ λ n n! As prevalent as Gaussian distribution
More informationClassical Queueing Models.
Sergey Zeltyn January 2005 STAT 99. Service Engineering. The Wharton School. University of Pennsylvania. Based on: Classical Queueing Models. Mandelbaum A. Service Engineering course, Technion. http://iew3.technion.ac.il/serveng2005w
More informationA Direct Approach to Transient Queue-Size Distribution in a Finite-Buffer Queue with AQM
Appl. Math. Inf. Sci. 7, No. 3, 99-915 (213) 99 Applied Mathematics & Information Sciences An International Journal A Direct Approach to Transient Queue-Size Distribution in a Finite-Buffer Queue with
More informationAdvanced Computer Networks Lecture 3. Models of Queuing
Advanced Computer Networks Lecture 3. Models of Queuing Husheng Li Min Kao Department of Electrical Engineering and Computer Science University of Tennessee, Knoxville Spring, 2016 1/13 Terminology of
More informationLecture 10: Semi-Markov Type Processes
Lecture 1: Semi-Markov Type Processes 1. Semi-Markov processes (SMP) 1.1 Definition of SMP 1.2 Transition probabilities for SMP 1.3 Hitting times and semi-markov renewal equations 2. Processes with semi-markov
More informationQueueing Networks and Insensitivity
Lukáš Adam 29. 10. 2012 1 / 40 Table of contents 1 Jackson networks 2 Insensitivity in Erlang s Loss System 3 Quasi-Reversibility and Single-Node Symmetric Queues 4 Quasi-Reversibility in Networks 5 The
More informationQueueing Review. Christos Alexopoulos and Dave Goldsman 10/25/17. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA
1 / 26 Queueing Review (mostly from BCNN) Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/25/17 2 / 26 Outline 1 Introduction 2 Queueing Notation 3 Transient
More informationQueueing. Chapter Continuous Time Markov Chains 2 CHAPTER 5. QUEUEING
2 CHAPTER 5. QUEUEING Chapter 5 Queueing Systems are often modeled by automata, and discrete events are transitions from one state to another. In this chapter we want to analyze such discrete events systems.
More informationInventory Ordering Control for a Retrial Service Facility System Semi- MDP
International Journal of Engineering Science Invention (IJESI) ISS (Online): 239 6734, ISS (Print): 239 6726 Volume 7 Issue 6 Ver I June 208 PP 4-20 Inventory Ordering Control for a Retrial Service Facility
More informationBudapest University of Tecnology and Economics. AndrásVetier Q U E U I N G. January 25, Supported by. Pro Renovanda Cultura Hunariae Alapítvány
Budapest University of Tecnology and Economics AndrásVetier Q U E U I N G January 25, 2000 Supported by Pro Renovanda Cultura Hunariae Alapítvány Klebelsberg Kunó Emlékére Szakalapitvány 2000 Table of
More informationChapter 8 Queuing Theory Roanna Gee. W = average number of time a customer spends in the system.
8. Preliminaries L, L Q, W, W Q L = average number of customers in the system. L Q = average number of customers waiting in queue. W = average number of time a customer spends in the system. W Q = average
More informationM/M/3/3 AND M/M/4/4 RETRIAL QUEUES. Tuan Phung-Duc, Hiroyuki Masuyama, Shoji Kasahara and Yutaka Takahashi
JOURNAL OF INDUSTRIAL AND doi:10.3934/imo.2009.5.431 MANAGEMENT OPTIMIZATION Volume 5, Number 3, August 2009 pp. 431 451 M/M/3/3 AND M/M/4/4 RETRIAL QUEUES Tuan Phung-Duc, Hiroyuki Masuyama, Shoi Kasahara
More informationStochastic Models in Computer Science A Tutorial
Stochastic Models in Computer Science A Tutorial Dr. Snehanshu Saha Department of Computer Science PESIT BSC, Bengaluru WCI 2015 - August 10 to August 13 1 Introduction 2 Random Variable 3 Introduction
More informationSYMBOLS AND ABBREVIATIONS
APPENDIX A SYMBOLS AND ABBREVIATIONS This appendix contains definitions of common symbols and abbreviations used frequently and consistently throughout the text. Symbols that are used only occasionally
More informationarxiv: v1 [math.pr] 11 May 2018
FCFS Parallel Service Systems and Matching Models Ivo Adan a, Igor Kleiner b,, Rhonda Righter c, Gideon Weiss b,, a Eindhoven University of Technology b Department of Statistics, The University of Haifa,
More informationTHE ON NETWORK FLOW EQUATIONS AND SPLITTG FORMULAS TRODUCTION FOR SOJOURN TIMES IN QUEUEING NETWORKS 1 NO FLOW EQUATIONS
Applied Mathematics and Stochastic Analysis 4, Number 2, Summer 1991, III-I16 ON NETWORK FLOW EQUATIONS AND SPLITTG FORMULAS FOR SOJOURN TIMES IN QUEUEING NETWORKS 1 HANS DADUNA Institut flit Mathematische
More informationChapter 2 Queueing Theory and Simulation
Chapter 2 Queueing Theory and Simulation Based on the slides of Dr. Dharma P. Agrawal, University of Cincinnati and Dr. Hiroyuki Ohsaki Graduate School of Information Science & Technology, Osaka University,
More informationIntroduction to Queueing Theory with Applications to Air Transportation Systems
Introduction to Queueing Theory with Applications to Air Transportation Systems John Shortle George Mason University February 28, 2018 Outline Why stochastic models matter M/M/1 queue Little s law Priority
More informationLIMITS AND APPROXIMATIONS FOR THE M/G/1 LIFO WAITING-TIME DISTRIBUTION
LIMITS AND APPROXIMATIONS FOR THE M/G/1 LIFO WAITING-TIME DISTRIBUTION by Joseph Abate 1 and Ward Whitt 2 April 15, 1996 Revision: January 2, 1997 Operations Research Letters 20 (1997) 199 206 1 900 Hammond
More informationSection 1.2: A Single Server Queue
Section 12: A Single Server Queue Discrete-Event Simulation: A First Course c 2006 Pearson Ed, Inc 0-13-142917-5 Discrete-Event Simulation: A First Course Section 12: A Single Server Queue 1/ 30 Section
More information6 Solving Queueing Models
6 Solving Queueing Models 6.1 Introduction In this note we look at the solution of systems of queues, starting with simple isolated queues. The benefits of using predefined, easily classified queues will
More informationQueueing Systems: Lecture 3. Amedeo R. Odoni October 18, Announcements
Queueing Systems: Lecture 3 Amedeo R. Odoni October 18, 006 Announcements PS #3 due tomorrow by 3 PM Office hours Odoni: Wed, 10/18, :30-4:30; next week: Tue, 10/4 Quiz #1: October 5, open book, in class;
More informationSynchronized Queues with Deterministic Arrivals
Synchronized Queues with Deterministic Arrivals Dimitra Pinotsi and Michael A. Zazanis Department of Statistics Athens University of Economics and Business 76 Patission str., Athens 14 34, Greece Abstract
More informationAn M/G/1 Queue with Adaptable Service Speed R. Bekker a ; O. J. Boxma b a
This article was downloaded by: [Vrije Universiteit, Library] On: 13 June 211 Access details: Access Details: [subscription number 9721892] Publisher Taylor & Francis Informa Ltd Registered in England
More informationIntro Refresher Reversibility Open networks Closed networks Multiclass networks Other networks. Queuing Networks. Florence Perronnin
Queuing Networks Florence Perronnin Polytech Grenoble - UGA March 23, 27 F. Perronnin (UGA) Queuing Networks March 23, 27 / 46 Outline Introduction to Queuing Networks 2 Refresher: M/M/ queue 3 Reversibility
More informationTHE HEAVY-TRAFFIC BOTTLENECK PHENOMENON IN OPEN QUEUEING NETWORKS. S. Suresh and W. Whitt AT&T Bell Laboratories Murray Hill, New Jersey 07974
THE HEAVY-TRAFFIC BOTTLENECK PHENOMENON IN OPEN QUEUEING NETWORKS by S. Suresh and W. Whitt AT&T Bell Laboratories Murray Hill, New Jersey 07974 ABSTRACT This note describes a simulation experiment involving
More informationPBW 654 Applied Statistics - I Urban Operations Research
PBW 654 Applied Statistics - I Urban Operations Research Lecture 2.I Queuing Systems An Introduction Operations Research Models Deterministic Models Linear Programming Integer Programming Network Optimization
More informationProbability and Statistics Concepts
University of Central Florida Computer Science Division COT 5611 - Operating Systems. Spring 014 - dcm Probability and Statistics Concepts Random Variable: a rule that assigns a numerical value to each
More informationNEW FRONTIERS IN APPLIED PROBABILITY
J. Appl. Prob. Spec. Vol. 48A, 209 213 (2011) Applied Probability Trust 2011 NEW FRONTIERS IN APPLIED PROBABILITY A Festschrift for SØREN ASMUSSEN Edited by P. GLYNN, T. MIKOSCH and T. ROLSKI Part 4. Simulation
More informationNUMERICAL INVERSION OF PROBABILITY GENERATING FUNCTIONS
NUMERICAL INVERSION OF PROBABILITY GENERATING FUNCTIONS by Joseph Abate Ward Whitt 900 Hammond Road AT&T Bell Laboratories Ridgewood, NJ 07450-2908 Room 2C-178 Murray Hill, NJ 07974-0636 July 26, 1991
More informationThe M/G/1 FIFO queue with several customer classes
The M/G/1 FIFO queue with several customer classes Boxma, O.J.; Takine, T. Published: 01/01/2003 Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume
More informationContinuous-Time Markov Chain
Continuous-Time Markov Chain Consider the process {X(t),t 0} with state space {0, 1, 2,...}. The process {X(t),t 0} is a continuous-time Markov chain if for all s, t 0 and nonnegative integers i, j, x(u),
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationThe Status Update Problem: Optimal Strategies for Pseudo-Deterministic Systems
The Status Update Problem: Optimal Strategies for Pseudo-Deterministic Systems Jared Tramontano under the direction of Dr. Shan-Yuan Ho Mr. Siddharth Venkatesh Department of Mathematics Massachusetts Institute
More informationNANYANG TECHNOLOGICAL UNIVERSITY SEMESTER I EXAMINATION MH4702/MAS446/MTH437 Probabilistic Methods in OR
NANYANG TECHNOLOGICAL UNIVERSITY SEMESTER I EXAMINATION 2013-201 MH702/MAS6/MTH37 Probabilistic Methods in OR December 2013 TIME ALLOWED: 2 HOURS INSTRUCTIONS TO CANDIDATES 1. This examination paper contains
More informationLecture 7: Simulation of Markov Processes. Pasi Lassila Department of Communications and Networking
Lecture 7: Simulation of Markov Processes Pasi Lassila Department of Communications and Networking Contents Markov processes theory recap Elementary queuing models for data networks Simulation of Markov
More informationCS 798: Homework Assignment 3 (Queueing Theory)
1.0 Little s law Assigned: October 6, 009 Patients arriving to the emergency room at the Grand River Hospital have a mean waiting time of three hours. It has been found that, averaged over the period of
More informationQueueing Theory. VK Room: M Last updated: October 17, 2013.
Queueing Theory VK Room: M1.30 knightva@cf.ac.uk www.vincent-knight.com Last updated: October 17, 2013. 1 / 63 Overview Description of Queueing Processes The Single Server Markovian Queue Multi Server
More informationA Regression Analysis Approach to Queueing System Modelling: a Case of Banks
Journal of Applied Sciences Research, 7(3): 200-212, 2011 ISSN 1819-544X This is a refereed journal and all articles are professionally screened and reviewed 200 ORIGINAL ARTICLES A Regression Analysis
More informationIntroduction to Queueing Theory
Introduction to Queueing Theory Raj Jain Washington University in Saint Louis Saint Louis, MO 63130 Jain@cse.wustl.edu Audio/Video recordings of this lecture are available at: 30-1 Overview Queueing Notation
More information