Queueing Networks G. Rubino INRIA / IRISA, Rennes, France
|
|
- Avis Reed
- 6 years ago
- Views:
Transcription
1 Queueing Networks G. Rubino INRIA / IRISA, Rennes, France February 2006
2 Index 1. Open nets: Basic Jackson result 2 2. Open nets: Internet performance evaluation Closed nets: Basic Gordon-Newell result Closed nets: Convolution algorithm Closed nets: Mean Value Analysis (MVA) Closed nets: Applications Loss nets: The model Loss nets: Product-form Loss nets: Erlang fixed point approach 72 Open QN G. Rubino, INRIA / IRISA slide 1
3 1: Basic Jackson result The basic model! N nodes; each one is a./m/1 queue! in node i, service rate: µ i ; external arrival flow: Poisson with rate λ i 0; arrival flows are independent! routing probabilities (transfer time = 0): r ij = Pr(a customer leaving i goes to j) (i, j = 1, 2,..., N) r i0 = Pr(a customer leaving i, leaves the network), (i = 1, 2,..., N) for i = 1, 2,..., N, N j=0 r ij = 1! total arrival rate: λ = N i=1 λ i > 0 Open QN G. Rubino, INRIA / IRISA slide 2
4 (1: Basic Jackson result) Example λ 1 1 λ Here, we assume λ 1 > 0, λ 2 > 0, λ 3 = 0. For instance, we have r 1,2 = 1/3, r 1,3 = 2/3, r 2,3 = 1/4, r 2,0 = 3/4, r 3,0 = 1. Open QN G. Rubino, INRIA / IRISA slide 3
5 (1: Basic Jackson result) The assumptions and the goal! arrival processes, service processes and switching processes are independent of each other! at each server, service times are independent of each other! at each switching point, switching results are also independent of each other! the goal is to evaluate performance metrics in equilibrium, such as mean # of customers at each node, or in the whole network, mean delay at each node, or at the whole network, utilization rate of some specific server, etc. Open QN G. Rubino, INRIA / IRISA slide 4
6 (1: Basic Jackson result) Computing the mean throughput (in equilibrium) λ 1 1 λ Define T i = mean throughput at i, i = 1, 2, 3, in customers per sec. We have T 1 = λ 1. Then, T 2 = λ 2 + T 1 r 1,2 = λ 2 + λ 1 /3, and T 3 = T 1 r 1,3 + T 2 r 2,3 = (3λ 1 + λ 2 )/4. Open QN G. Rubino, INRIA / IRISA slide 5
7 (1: Basic Jackson result) Extended network A node 0, representing the environment of the net, with service rate λ, is added as the figure shows. def The routing probability from 0 to i is r 0i = λ i /λ. 0 λr 0,1 λr 0, To see that the two networks are equivalent, we use the decomposition properties of Poisson processes. Open QN G. Rubino, INRIA / IRISA slide 6
8 (1: Basic Jackson result) Eliminating pathological and/or trivial cases! Seeing the extended network as a graph, we assume it is strongly connected. This way, we eliminate the pathological situations in which nodes will eventually empty, nodes will eventually explode. Open QN G. Rubino, INRIA / IRISA slide 7
9 (1: Basic Jackson result) A pathological example Assume an arrow entering a node means that the external arriving flow has a non-null rate. Node 2 will explode (at t =, customers inside). 2 1 Here, eliminate node 2. 3 Open QN G. Rubino, INRIA / IRISA slide 8
10 (1: Basic Jackson result) Another pathological example In this example, node 1 receives no customer. So, at t =, it is empty. 2 1 Here, eliminate node 1. 3 Open QN G. Rubino, INRIA / IRISA slide 9
11 (1: Basic Jackson result) Jackson s theorem! Resuming, from T i = λ i + N T j r ji, the strongly connected property of the extended graph implies that there is an unique (and strictly positive) solution, = we get the T i s. j=1 This can be seen, from instance, from Markov chain theory, seeing the extended graph as the graph of a Markov chain.! Let ϱ i = T i. µ i Open QN G. Rubino, INRIA / IRISA slide 10
12 (1: Basic Jackson result)! If for all i we have ϱ i < 1, then, in equilibrium (Jackson result) Pr(n 1 in 1,..., n N in N) = N (1 ϱ i )ϱ n i i=1 i. This is called a product form result.! The simplest proof of this result is straightforward: argue that the multidimensional stochastic process (L (1) t,, L (N) t ) is Markov (continuous time, infinite denumerable state space, homogeneous), that it is irreducible (so, when ergodic, only one stationary distribution), and then check that the given expression (product form) satisfies the equilibrium equations (Chapman-Kolmogorov equations). Open QN G. Rubino, INRIA / IRISA slide 11
13 (1: Basic Jackson result) Proof of the result concerning the T i s! Consider a Markov chain in continuous time, with state space {0, 1,, N} and transition probabilities (r i,j ). The departure rate of state i is some d i > 0.! The strong connectivity condition means that the chain is irreducible, and then, we know that the system of N + 1 unknowns (x 0,, x N ) and equations x i = N x j r ji, i = 0, 1,, N j=0 has a unique non-trivial solution up to a factor, and that it is composed of strictly positive numbers. Open QN G. Rubino, INRIA / IRISA slide 12
14 (1: Basic Jackson result)! Then, choose the solution such that x 0 = λ. It follows that the remaining x i s verify the throughtput equations. Reciprocally, given a solution (T 1,, T N ) to the throughput equations, then (λ, T 1,, T N ) is a solution to the auxiliary homogeneous system, so, the solution with x 0 = λ, proving the unicity of the solution in the throughputs.! Observe that this also gives for any node i, T i > 0.! Given the fact that if π x is the stationary probability of state x in an ergodic Markov chain, 1/π x is the relative average frequency of visits to that state, we deduce that T i /λ is the expected # of visits a customer will make to node i before leaving.! For the same reason, the number i T i/λ is the mean number visits to nodes a customer does before leaving. Open QN G. Rubino, INRIA / IRISA slide 13
15 (1: Basic Jackson result) The marginal distributions Pr(n 1 in 1) = n 2 0,,n N 0 = = (1 ϱ 1 )ϱ n 1 1 Pr(n 1 in 1,..., n N in N) So, Pr(n 1 in 1,..., n N in N) = N Pr(n i in i). i=1 We say that the occupation processes at each node (the processes L (1) t,, L (N) t ) are asymptotically independent. Open QN G. Rubino, INRIA / IRISA slide 14
16 (1: Basic Jackson result) An equivalent view! Everything happens as if the network was composed of N queues M/M/1 independent, in parallel, where the ith has arrival rate T i and service rate µ i (same joint stationary distributions).! So, mean # of packets at i: L i = ϱ i 1 ϱ i = T i µ i T i! Since in the two models, each node i has the same L i and the same T i, then by Little the mean delay to pass through it is also the same: 1 D i = µ i T i! mean delay to pass through the whole network (Little): D = 1 λ N i=1 T i µ i T i Open QN G. Rubino, INRIA / IRISA slide 15
17 (1: Basic Jackson result) Questions! Since the factors in the product form are geometric, we could conjecture that arrival flows at each node are Poisson: if the network is feedforward, this can be shown to be true: it is based in Burke s Theorem: in equilibrium, the output of a M/M/1 queue is Poisson (true also for M/M/c or M/M/ queues) if not, the conjecture is false! We still don t understand completely what is happening here. Open QN G. Rubino, INRIA / IRISA slide 16
18 (1: Basic Jackson result) Extensions! Jackson s result also holds if the nodes are./m/c or./m/ queues, if the global network has finite storage capacity, if, more generally, the arrival rates depend on the total # of units in the network, and if the service rate at i is a function µ i (n i ) of the backlog at i.! There are also different multiclass extensions. Open QN G. Rubino, INRIA / IRISA slide 17
19 2: Internet performance evaluation Analyzing a network of routers! A packet switching network has N switches, working in datagram mode.! Denote by c ij the speed (or capacity) of the line from i to j (0 means that there is no connexion from i to j).! Let C i denote the set of switches that can directly receive packets from i, that is, C i = {j : c ij > 0}.! Inside switch i there are C i = M i output queues, each corresponding to an output link. Denote by Q ij the queue associated with the link connected to j C i ; it has service rate in pps µ ij = c ij /B. Open QN G. Rubino, INRIA / IRISA slide 18
20 2: Internet performance evaluation! A packet leaving i to j, goes to queue Q jk with probability r jk and leaves the network (to go to some host connected to switch j) with probability r j0.! This leads to a Jackson model having M = i M i nodes (one for each output line leaving a switch).! Problem: service times are certainly not indep.! Kleinrock approximation: assume they are!!! Observation: if there are many packets traveling and the speeds are not too low, the approximation is very good. Open QN G. Rubino, INRIA / IRISA slide 19
21 (2: Internet performance evaluation) Other metrics! The utilization rate U i of node i is the probability that, in equilibrium, server i is busy. We have U i = ϱ i.! Assume a customer arrives at a node i (even for a node where λ i = 0). Let E i = mean delay to travel through the whole network for the customer.! In the example, we have E 3 = D 3. Then, E 2 = D E 3 = D D 3, and E 1 = D E E 3 = D D D 3. Open QN G. Rubino, INRIA / IRISA slide 20
22 (2: Internet performance evaluation)! In general, the system of equations leading to the E i s is E i = D i + N r ij E j, j=1 which solution is unique (observe that the system matrix is the transpose of the one leading to the T i s). Open QN G. Rubino, INRIA / IRISA slide 21
23 (2: Internet performance evaluation) Exercise Consider the network represented here After leaving node 0, packets can follow either of the two paths (1) or (2,3). Working on the routing algorithm, we can control the proportion of packets using each path. The output link speed at router i is c i bps; the packets arrive at a rate of λ pps; the mean packet size is B bits. We denote by p the fraction of packets following path (1). Problem: find p so that the global delay is minimal, when λ = 10 4, c 0 = c 1 = , c 2 = , c 3 = and B = Open QN G. Rubino, INRIA / IRISA slide 22
24 (2: Internet performance evaluation) Example (cont.)! With the usual assumptions, denoting µ i = c i /B, we have the Jackson network p 1 λ µ p! The mean delay at each node is D 0 = 1/(µ 0 λ), µ 0 2 µ 2 3 µ 3 D 1 = 1 µ 1 pλ, D 2 = 1 µ 2 (1 p)λ, D 3 = 1 µ 3 (1 p)λ.! Simplifying and using as unit of time 100 µs, D 0 = 1, D 1 = 1 2 p, D 2 = p, D 3 = p. Open QN G. Rubino, INRIA / IRISA slide 23
25 (2: Internet performance evaluation)! Global delay: D = 1 + pd 1 + (1 p)(d 2 + D 3 ). Solving dd/dp = 0 we obtain (numerically) the minimal value for p , and the min delay is ! For the same prize, let us compare D 1 and D 2 + D 3 : D 1? D 2 + D p? p p 3p 2 + 6p 4? 0.! So, with p 0 = , p < p 0 (resp. p > p 0 ) (1) shorter (resp. longer)! Remark: when p = p 0, D 1 = D 2 + D 3 and the total delay is Open QN G. Rubino, INRIA / IRISA slide 24
26 3: Basic Gordon-Newell result The basic model! N nodes; each one is a./m/1 queue! in node i, mean output rate: µ i ;! no external arrivals, no departures to outside! routing probabilities (transfer time = 0): r ij = Pr(a customer leaving i goes to j) (i, j = 1, 2,..., N) N for i = 1, 2,..., N, j=1 r ij = 1! the graph defined by the routing probabilities is assumed to be strongly connected (again, everything reduces to this case)! H customers circulating in the network Open QN G. Rubino, INRIA / IRISA slide 25
27 (3: Basic Gordon-Newell result) Example H customers Open QN G. Rubino, INRIA / IRISA slide 26
28 (3: Basic Gordon-Newell result) The assumptions and the goal! service processes and switching processes are independent of each other! service times and switching results are also independent of each other and of the previous mentioned processes! the goal is to evaluate performance metrics in equilibrium, such as mean # of customers at each node, or in the whole network, mean delay at each node, or at the whole network, utilization rate of some specific server, etc. Open QN G. Rubino, INRIA / IRISA slide 27
29 (3: Basic Gordon-Newell result) Computing the mean throughput (in equilibrium) H customers Define T i = mean throughput at i, i = 1, 2, 3, in customers per sec. We have T 1 = T 3 and T 2 = T 1 /3. Up to this point, the throughputs are known up to a multiplicative constant. Open QN G. Rubino, INRIA / IRISA slide 28
30 (3: Basic Gordon-Newell result) Gordon-Newell s theorem! From T i = N T j r ji, j=1 the strongly connected property of the extended graph unique (and strictly positive) solution up to a factor. For instance, let (V i ) be the only solution such that N i=1 V i = 1.! V i is the visit ratio of node i. Recall that V i > 0. If ν i (t) is the number of visits made by a customer to node i in [0, t], then ν i (t) lim t ν j (t) = V i. V j Open QN G. Rubino, INRIA / IRISA slide 29
31 (3: Basic Gordon-Newell result)! Let X i = V i µ i.! Since we are assuming that all service times are indep., then (Gordon-Newell s result) Pr(n 1 in 1,..., n N in N) = This is called a product form result. N i=1 Xn i i G.! G is the normalizing constant. It is given by G = n 1 + +n N =H N i=1 X n i i.! Immediate to verify that we can express the product-form using X i = KV i /µ i for any K > 0. Open QN G. Rubino, INRIA / IRISA slide 30
32 (3: Basic Gordon-Newell result) Example In the example, H customers from V 1 = V 3 and V 2 = V 1 /3 we have V 1 = 3 7, V 2 = 1 7, V 3 = 3 7. Open QN G. Rubino, INRIA / IRISA slide 31
33 (3: Basic Gordon-Newell result) Example (cont.) Let us look for the probability that the server in node 3 is busy, for H = 2 customers and in the case µ i = µ, i = 1, 2, 3.! Let us denote a = V 1 /µ and b = V 2 /µ.! G = 3a 2 + 2ab + b 2! Pr(node 3 busy) = 1 Pr(node 3 empty) Pr(node 3 empty) = a2 + ab + b 2 G! So, Pr(node 3 busy) = 2a 2 + ab 3a 2 + 2ab + b 2 = Open QN G. Rubino, INRIA / IRISA slide 32
34 (3: Basic Gordon-Newell result) Example (end)! Since we know the probability that node 3 is not empty, we have T 3 = µ ! From this, we obtain all the throughputs: T 1 = T 3 = 21µ 25, T 2 = 7µ 25. Open QN G. Rubino, INRIA / IRISA slide 33
35 (3: Basic Gordon-Newell result) Proof of Gordon-Newell s theorem! The simplest proof of Gordon-Newell s result is straightforward: argue that the multidimensional stochastic process (L (1) t,, L (N) t ) is Markov (continuous time, infinite denumerable state space, homogeneous), where L (i) t is the occupation process (# of customers) of node i at time t; the chain is irreducible and the state space is finite; the chain is then ergodic; there is only one stationary distribution; then, check that the given expression (product form) satisfies the equilibrium equations (Chapman-Kolmogorov equations). Remark: see that, trivially, the variables L (1) t,, L (N) t are never independent (not even for t = ), since H i=1 L(i) t = H. Open QN G. Rubino, INRIA / IRISA slide 34
36 (3: Basic Gordon-Newell result) The state space! Vector (L (1) t! There are ( ) N+H 1 H states.,, L (N) t ) lives in the state space { (n 1,..., n N ) IN N : } N n i = H. i=1! For instance, N = 10, H = states N = 20, H = states Open QN G. Rubino, INRIA / IRISA slide 35
37 (3: Basic Gordon-Newell result) The convolution operator! Denote G = G(H).! Define f i (n) = X n i.! Clearly, G() is the convolution of f 1 (),, f N () (recall that convolution is commutative and associative). Open QN G. Rubino, INRIA / IRISA slide 36
38 (3: Basic Gordon-Newell result) The marginal distribution Denote by π i (n) = Pr(n in i) (always, in equilibrium). We have π i (n) = n 1 + +n N =H, n i =n = Xn i G(i) (H n) G(H) Pr(n 1 in 1,..., n N in N) where G (i) () is the convolution of f 1 (),, f i 1 (), f i+1 (),, f N () Open QN G. Rubino, INRIA / IRISA slide 37
39 (3: Basic Gordon-Newell result) P r(at least k at i) = Pr(L (i) k) = H Xi ng(i) (H n) G(H) = n=k = Xk i X k i G(H) H k n=0 G(H k). G(H) X n i G (i) (H k n) Open QN G. Rubino, INRIA / IRISA slide 38
40 (3: Basic Gordon-Newell result) Throughputs! We know that T i = T V i for some T. Since N i=1 T i = T, T is called mean global throughput.! We have Using the previous slide, T i = (1 π i (0))µ i = Pr(L (i) 1)µ i. T i = X ig(h 1) G(H)! The mean global throughput is, then, G(H 1) µ i = V i. G(H) T = G(H 1) G(H). Open QN G. Rubino, INRIA / IRISA slide 39
41 (3: Basic Gordon-Newell result)! In the example, G(1) = 2a + b T = G(1) G(2) = 2a + b 3a 2 + 2ab + b 2 = µ. Open QN G. Rubino, INRIA / IRISA slide 40
42 4: Convolution algorithm Main idea! Recall that π n1,,n N = N i=1 Xn i i G(H).! The problem is G(). Using its definition G(H) = n 1 + +n N =H N i=1 X n i i. is useless except for very small models.! Key idea: G() is the convolution of f 1 (),, f N () with f i (n) = Xi n, and the convolution is commutative and associative. Open QN G. Rubino, INRIA / IRISA slide 41
43 (4: Convolution algorithm) Convolution algorithm! Denote by G i () the convolution of f 1 (),, f i ().! We have that G 1 () = f 1 (), G() = G N () and that G i () is the convolution of G i 1 () and f i () for i 2.! Working ( this way, the computation of G() moves from a complexity (N+H 1 ) ) Θ to θ(nh). H Open QN G. Rubino, INRIA / IRISA slide 42
44 (4: Convolution algorithm) Using convolutions Recall that everything can be computed from these convolutions:! Thoughputs: T = G(H 1) G(H) and T i = T V i! Marginals: π i (n) = Xn i G(i) (H n) G(H) (or H π i (k) = Xk i n=k G(H k) G(H) Open QN G. Rubino, INRIA / IRISA slide 43 )
45 (4: Convolution algorithm)! Mean number of customers: L i = H k=1 X k i G(H k) G(H)! Mean response time: R i = L i /T i (Little) Open QN G. Rubino, INRIA / IRISA slide 44
46 (4: Convolution algorithm) Implementing! We can exhibit another recurrence while computing G i (): G i (n) = G i 1 (n) + X i G i (n 1).! This means that to compute G() we can use the scheme for (i=1; i<=n; i++) for (h=1; h=h; h++) G[h] += X[i]*G[h-1] Open QN G. Rubino, INRIA / IRISA slide 45
47 (4: Convolution algorithm) Possible problems! Main difficulty: overflows or underflows while computing convolutions.! For instance, consider a ring composed of N identical./m/1 nodes, with H customers. Let µ i = 1 for all i. Then, By symmetry, L i = H/N, but G(n) = 1 N n ( n + H 1 Consider N = 20 and H = 100: L i = 5 but G() has very small values, most of them too small to be stored as double.! Some solutions exist (dynamic scaling, for instance). H ). Open QN G. Rubino, INRIA / IRISA slide 46
48 5: Mean Value Analysis (MVA) Main idea! Instead of a recurrence on the nodes, as in the convolution approach, we look for a recurrence on the customers.! Make explicit the dependency on the number of customers: L i (H), R i (H). Write S i = 1/µ i the mean service time.! Key result: R i (H) = S i [1 + L i (H 1)].! From L i (H) = T i (H)R i (H), T i (H) = T (H)V i and N i=1 L i(h) = H we obtain T (H) = H N i=1 V ir i (H). Open QN G. Rubino, INRIA / IRISA slide 47
49 (5: Mean Value Analysis (MVA)) MVA scheme! Start with L i (0) = 0! Call R(h) the network with h customers.! The parameters for R(h) are obtained from the following scheme: Compute the response times: R i (h) = S i [1 + L i (h 1)] for i = 1,, N Then, compute the global throughput: T (h) = h N i=1 V ir i (h). Finally, compute the mean queue lengths: L i (h) = T (h)v i R i (h) for i = 1,, N Open QN G. Rubino, INRIA / IRISA slide 48
50 (5: Mean Value Analysis (MVA)) Properties! MVA is robust (no overflow or underflow risk)! Complexity in θ(nh) (as for convolution)! The probabilistic interpretation of the key relationships leads to many other approximation algorithms, bounding techniques, etc. Open QN G. Rubino, INRIA / IRISA slide 49
51 (5: Mean Value Analysis (MVA)) An example of approximation technique based on MVA! In some problems, we have to evaluate many times a network with a large number H of customers. This can be costly.! Schweitzer s method: for H large, assume L i (H 1) H 1 H L i(h).! Write L i (H) as a function of L 1 (H 1),, L N (H 1) using the MVA recurrences: L i (H) = HX i [1 + L i (H 1)] N j=1 X j [1 + L j (H 1)] Open QN G. Rubino, INRIA / IRISA slide 50
52 (5: Mean Value Analysis (MVA))! Use Schweitzer s idea, solve the fixed-point equation [ HX i 1 + H 1 ] H y i y i = [ N j=1 X j 1 + H 1 ] H y j! It can be proved (using Brower s fixed point theorems) results about existence and unicity of a fixed point for the previous equation.! Many improvements of this scheme have been published. Open QN G. Rubino, INRIA / IRISA slide 51
53 (5: Mean Value Analysis (MVA)) Extensions! Gordon-Newell s result also holds if the nodes are./m/c or./m/ queues, for different multiclass extensions.! For instance, if node i is of the./m/ type, we have that for all population size H, R i (H) = 1 µ i.! If node i is./m/c i, then 1 if c i H, R i (H) = µ i 1 (1 L i (H 1)) if c i < H. µ i Open QN G. Rubino, INRIA / IRISA slide 52
54 6: Applications Client-server models A standard model of a client-server architecture is depicted in the following figure: C clients requests server services (answers) Open QN G. Rubino, INRIA / IRISA slide 53
55 (6: Applications) Client-server models (cont.)! This leads to use a closed queueing network as a representative model.! For instance, the goal can be to evaluate the maximum of the number of customers that can be served in order to keep the response time low enough.! The model can be made more precise, if relevant, as for instance in the next picture. Open QN G. Rubino, INRIA / IRISA slide 54
56 (6: Applications) Client-server models (cont.) the clients (a./m/k node) p a database (say) p the server Here, a fraction p of the request sent to the server are followed by a request to another resource (a database), and then the answer arrives at the client. Open QN G. Rubino, INRIA / IRISA slide 55
57 (6: Applications) Flow control using windows mechanisms! Consider a network connection where some window mechanism is used to work with a maximum of H packets.! If the connection uses N routers, a simple model allowing to perform some global performance analysis is a closed queuing network with N nodes and H customers organized as a ring.! As an exercise, assume that we use such a model with identical routers, represented as./m/1 nodes with same mean service time S. Compute the mean RTT (Round Trip Time). Open QN G. Rubino, INRIA / IRISA slide 56
58 (6: Applications) Flow control using windows mechanisms (cont.) 1 2 N H customers! We have already observed that, by symmetry, the mean # of customers at node i is L i = H/N.! The main MVA recurrence gives us the response time: [ R i = S[1 + L i (H 1)] = S 1 + H 1 ]. N! The mean throughput in i is given by T i = L i R i = H S(N + H 1).! RTT = N i=1 R i = S(N + H 1) Open QN G. Rubino, INRIA / IRISA slide 57
59 7: Loss nets Telephone switches! A switch in a telephone system can be seen as a G/G/c/c queue. Integer c is the # of channels or circuits.! The blocking probability or loss probability is the probability that an arriving call finds all channels busy.! Computing this blocking probability is the basic task in telephony dimensioning. Open QN G. Rubino, INRIA / IRISA slide 58
60 (7: Loss nets) Blocking probabilities! Assuming that calls arrive according to a Poisson process with rate λ, and that the duration of a call is exponentially distributed with mean S, the blocking probability (M/M/c/c model) is E(ϱ, c) = ϱ c /c! c i=0 ϱi /i!, where ϱ = λs.! Insensitivity: recall that if we keep the same mean call duration but we change the distribution (that is, if we consider a M/G/c/c model with same ϱ), the result is the same. Open QN G. Rubino, INRIA / IRISA slide 59
61 (7: Loss nets) Networks of switches! In telephone systems, switches are interconnected by links. link l (it has c l circuits) switch i! Each link has several channels or circuits.! Calls arrive at a switch and follow a path along the network before leaving. Open QN G. Rubino, INRIA / IRISA slide 60
62 (7: Loss nets) Networks of switches (cont.)! We associate a call with a route along the graph, that can be seen as a sequence of consecutive switches or links, without repetitions (acyclic).! Each call occupies an available circuit on each link of its route. If there is no available link, that is, if at least at one link of the route all circuits are busy, the call is blocked or lost.! As in the single-node case, the most important engineering task in telephony is the dimensioning of the system, which needs to evaluate the blocking probability. Open QN G. Rubino, INRIA / IRISA slide 61
63 8: Product-form Notation and assumptions! # of switches: N! # of links: L! link l has c l circuits! we assume that calls duration has mean 1 with exponential distribution Open QN G. Rubino, INRIA / IRISA slide 62
64 (8: Product-form) Notation and assumptions (cont.)! there are R routes! the flow of calls arriving to route r is Poisson, with rate λ r! route r is characterized by the sequence of consecutive links that are followed by the associated calls! we define matrix A = (A rl ), which is R L; A rl = 1 {link l belongs to route r}. Open QN G. Rubino, INRIA / IRISA slide 63
65 (8: Product-form) State space! We analyze the stochastic process L t = (L (1) t L (r) t is the # of calls in progress following route r.! We must respect the constraint,, L (R) t ), where for every link l, R r=1 L (r) t A rl c l.! L t is Markov (homogeneous) with finite state space { } R S = (n 1,, n R ) IN R : for every link l, n r A rl c l. r=1! Denoting n = (n 1,, n R ) and c = (c 1,..., c L ), S = { (n 1,, n R ) IN R : na c }. Open QN G. Rubino, INRIA / IRISA slide 64
66 (8: Product-form) Stationary (product-form) solution! Assume that for any l we have c l =. In this case, the r.v. L (1) t,, L (R) t are independent, the distribution of L (r) is Pr(L (r) = n) = e λ r λ n r n! = p(λ r, n), and the joint stationary distribution is ( ) Pr (L (1),, L (R) ) = (n 1,, n R ) = p(λ 1, n 1 ) p(λ R, n R ). Open QN G. Rubino, INRIA / IRISA slide 65
67 (8: Product-form)! In our model, routes interfere, but, remarkably enough, the joint distribution is still proportional to the previous one: Pr ( ) (L (1),, L (R) ) = (n 1,, n R ) = p(λ 1, n 1 ) p(λ R, n R ) G.! G is a normalizing constant: G = (n 1,,n R ) S p(λ 1, n 1 ) p(λ R, n R ).! The proof can be done in a straightforward manner as for Jackson or Gordon-Newell models. Open QN G. Rubino, INRIA / IRISA slide 66
68 (8: Product-form) Blocking probabilities! Let us denote by S r the following space: S r = {(n 1,, n R ) S : (n 1,, n r + 1,, n R ) S}.! A call in route r is lost if it arrives when the network state is in S r.! Using PASTA, where blocking probability in route r = B r = G r G, G r =! Writing S( c) and S r ( c), (n 1,,n R ) S r p(λ 1, n 1 ) p(λ R, n R ). S r ( c) = S( c 1 r A), where 1 r is the vector (δ r1,, δ rr ). Open QN G. Rubino, INRIA / IRISA slide 67
69 (8: Product-form) Computational aspects! The irregular shape of the state space makes difficult the analysis of this kind of models if its size is not very small.! There are recurrences similar to those of Gordon-Newell networks, but they are more complex due to the model s structure.! Other possibilities exist, for instance, using Monte Carlo to approximate the numerator and the denominator of the expression in B r = G r /G. Open QN G. Rubino, INRIA / IRISA slide 68
70 (8: Product-form) Back to a single node! Consider the case of a single node, but with K flows (for instance, K different types or classes or calls) arriving to a switch with c channels.! Class-k calls arrive with mean intensity λ k and their calls last a random time with exponential distribution and mean 1/µ k.! This is a multiclass M/M/c/c queue with K classes {! State space: S = n = (n 1,, n K ) IN K K : k=1 n k c}. Open QN G. Rubino, INRIA / IRISA slide 69
71 (8: Product-form)! The analysis of this model gives π n = 1 G ϱ n 1 1 ϱn K K n 1! n K!. The proof consists of verifyng that the evolution of the state with time is Markov irreducible, of writing the balance equations and then verifying that the given expression satisfies them.! See that the loss probability is the same for any class, and that it is given by p L = π n. n: n 1 + +n K =c! Let us define N = total number of units in the queue. Open QN G. Rubino, INRIA / IRISA slide 70
72 (8: Product-form) We have Pr(N = m) = = = 1 G = 1 G n: n 1 + +n K =m n: n 1 + +n K =m 1 m! where ϱ = ϱ ϱ K. 1 m! ϱm, π n 1 G n: n 1 + +n K =m ϱ n 1 1 ϱn K K n 1! n K! m! n 1! n K! ϱn 1 1 ϱn K K! So, N behaves as the # of ongoing calls in the classic Erlang B model, and thus the loss probability wee are looking for is simply p L = E(ϱ, c) = ϱ c c! c j=0 Open QN G. Rubino, INRIA / IRISA slide 71 ϱ j j!.
73 9: Erlang fixed point approach Asymptotic behavior! It can be shown that if c l and λ r such that c l /λ r is kept constant, B r 1 (1 B l ) l route r where B l is the locking probability of link l.! In the limit, the links behave independently. Open QN G. Rubino, INRIA / IRISA slide 72
74 (9: Erlang fixed point approach) The Erlang fixed point idea! Assume that the links behave independently. The flow of calls in route r arrives at link l with mean rate (1 B j ). λ r j route r, j l! The global rate arriving at node l is then ν l = R λ r r=1 j route r, j l (1 B j ).! Now, use B l = E(ν l, c l ), to obtain the fixed point equations for all link l, ν l = R λ r r=1 j route r, j l (1 E(ν j, c j )). Open QN G. Rubino, INRIA / IRISA slide 73
75 (9: Erlang fixed point approach)! Equivalently, we have for all link l, B l = E R λ r r=1 j route r, j l (1 B j ), c l.! It has been shown that these fixed point equations have a solution and that it is unique.! The accuracy of this approximation is, in general, very good if the network is large. Open QN G. Rubino, INRIA / IRISA slide 74
Performance 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 informationQueuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem. Wade Trappe
Queuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem Wade Trappe Lecture Overview Network of Queues Introduction Queues in Tandem roduct Form Solutions Burke s Theorem What
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 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 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 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 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 informationTCOM 501: Networking Theory & Fundamentals. Lecture 6 February 19, 2003 Prof. Yannis A. Korilis
TCOM 50: Networking Theory & Fundamentals Lecture 6 February 9, 003 Prof. Yannis A. Korilis 6- Topics Time-Reversal of Markov Chains Reversibility Truncating a Reversible Markov Chain Burke s Theorem Queues
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 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 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 Chain Model for ALOHA protocol
Markov Chain Model for ALOHA protocol Laila Daniel and Krishnan Narayanan April 22, 2012 Outline of the talk A Markov chain (MC) model for Slotted ALOHA Basic properties of Discrete-time Markov Chain Stability
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 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 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 informationExercises Stochastic Performance Modelling. Hamilton Institute, Summer 2010
Exercises Stochastic Performance Modelling Hamilton Institute, Summer Instruction Exercise Let X be a non-negative random variable with E[X ]
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 informationBuzen s algorithm. Cyclic network Extension of Jackson networks
Outline Buzen s algorithm Mean value analysis for Jackson networks Cyclic network Extension of Jackson networks BCMP network 1 Marginal Distributions based on Buzen s algorithm With Buzen s algorithm,
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 informationNon 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 informationSolutions to Homework Discrete Stochastic Processes MIT, Spring 2011
Exercise 6.5: Solutions to Homework 0 6.262 Discrete Stochastic Processes MIT, Spring 20 Consider the Markov process illustrated below. The transitions are labelled by the rate q ij at which those transitions
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 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 informationTime Reversibility and Burke s Theorem
Queuing Analysis: Time Reversibility and Burke s Theorem Hongwei Zhang http://www.cs.wayne.edu/~hzhang Acknowledgement: this lecture is partially based on the slides of Dr. Yannis A. Korilis. Outline Time-Reversal
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 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 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 informationCOMP9334 Capacity Planning for Computer Systems and Networks
COMP9334 Capacity Planning for Computer Systems and Networks Week 2: Operational Analysis and Workload Characterisation COMP9334 1 Last lecture Modelling of computer systems using Queueing Networks Open
More informationNICTA Short Course. Network Analysis. Vijay Sivaraman. Day 1 Queueing Systems and Markov Chains. Network Analysis, 2008s2 1-1
NICTA Short Course Network Analysis Vijay Sivaraman Day 1 Queueing Systems and Markov Chains Network Analysis, 2008s2 1-1 Outline Why a short course on mathematical analysis? Limited current course offering
More information8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains
8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains 8.1 Review 8.2 Statistical Equilibrium 8.3 Two-State Markov Chain 8.4 Existence of P ( ) 8.5 Classification of States
More informationreversed chain is ergodic and has the same equilibrium probabilities (check that π j =
Lecture 10 Networks of queues In this lecture we shall finally get around to consider what happens when queues are part of networks (which, after all, is the topic of the course). Firstly we shall need
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 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 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 queuing theory
Introduction to queuing theory Claude Rigault ENST claude.rigault@enst.fr Introduction to Queuing theory 1 Outline The problem The number of clients in a system The client process Delay processes Loss
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 information1 IEOR 4701: Continuous-Time Markov Chains
Copyright c 2006 by Karl Sigman 1 IEOR 4701: Continuous-Time Markov Chains A Markov chain in discrete time, {X n : n 0}, remains in any state for exactly one unit of time before making a transition (change
More informationEfficient Nonlinear Optimizations of Queuing Systems
Efficient Nonlinear Optimizations of Queuing Systems Mung Chiang, Arak Sutivong, and Stephen Boyd Electrical Engineering Department, Stanford University, CA 9435 Abstract We present a systematic treatment
More informationCOMP9334: Capacity Planning of Computer Systems and Networks
COMP9334: Capacity Planning of Computer Systems and Networks Week 2: Operational analysis Lecturer: Prof. Sanjay Jha NETWORKS RESEARCH GROUP, CSE, UNSW Operational analysis Operational: Collect performance
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 informationMAS275 Probability Modelling Exercises
MAS75 Probability Modelling Exercises Note: these questions are intended to be of variable difficulty. In particular: Questions or part questions labelled (*) are intended to be a bit more challenging.
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 informationOn Tandem Blocking Queues with a Common Retrial Queue
On Tandem Blocking Queues with a Common Retrial Queue K. Avrachenkov U. Yechiali Abstract We consider systems of tandem blocking queues having a common retrial queue. The model represents dynamics of short
More informationCHAPTER 14. Quantifying the Quality of Audio and Video Transmissions over the Internet: The PSQA Approach. Gerardo Rubino
CHAPTER 14 Quantifying the Quality of Audio and Video Transmissions over the Internet: The PSQA Approach Gerardo Rubino INRIA / IRISA, Campus de Beaulieu 35042 Rennes Cedex, France E-mail: rubino@irisa.fr
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 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 informationLink Models for Packet Switching
Link Models for Packet Switching To begin our study of the performance of communications networks, we will study a model of a single link in a message switched network. The important feature of this model
More informationChapter 16 focused on decision making in the face of uncertainty about one future
9 C H A P T E R Markov Chains Chapter 6 focused on decision making in the face of uncertainty about one future event (learning the true state of nature). However, some decisions need to take into account
More informationLecture 11: Introduction to Markov Chains. Copyright G. Caire (Sample Lectures) 321
Lecture 11: Introduction to Markov Chains Copyright G. Caire (Sample Lectures) 321 Discrete-time random processes A sequence of RVs indexed by a variable n 2 {0, 1, 2,...} forms a discretetime random process
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 informationConvolution Algorithm
Convolution Algorithm Raj Jain Washington University in Saint Louis Saint Louis, MO 63130 Jain@cse.wustl.edu Audio/Video recordings of this lecture are available at: http://www.cse.wustl.edu/~jain/cse567-08/
More informationElementary queueing system
Elementary queueing system Kendall notation Little s Law PASTA theorem Basics of M/M/1 queue M/M/1 with preemptive-resume priority M/M/1 with non-preemptive priority 1 History of queueing theory An old
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 informationPerformance analysis of queueing systems with resequencing
UNIVERSITÀ DEGLI STUDI DI SALERNO Dipartimento di Matematica Dottorato di Ricerca in Matematica XIV ciclo - Nuova serie Performance analysis of queueing systems with resequencing Candidato: Caraccio Ilaria
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 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 informationλ λ λ In-class problems
In-class problems 1. Customers arrive at a single-service facility at a Poisson rate of 40 per hour. When two or fewer customers are present, a single attendant operates the facility, and the service time
More informationQueueing systems. Renato Lo Cigno. Simulation and Performance Evaluation Queueing systems - Renato Lo Cigno 1
Queueing systems Renato Lo Cigno Simulation and Performance Evaluation 2014-15 Queueing systems - Renato Lo Cigno 1 Queues A Birth-Death process is well modeled by a queue Indeed queues can be used to
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 informationReadings: Finish Section 5.2
LECTURE 19 Readings: Finish Section 5.2 Lecture outline Markov Processes I Checkout counter example. Markov process: definition. -step transition probabilities. Classification of states. Example: Checkout
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 informationSTOCHASTIC PROCESSES Basic notions
J. Virtamo 38.3143 Queueing Theory / Stochastic processes 1 STOCHASTIC PROCESSES Basic notions Often the systems we consider evolve in time and we are interested in their dynamic behaviour, usually involving
More informationThe Erlang Model with non-poisson Call Arrivals
The Erlang Model with non-poisson Call Arrivals Thomas Bonald To cite this version: Thomas Bonald. The Erlang Model with non-poisson Call Arrivals. Sigmetrics / Performance 2006 - Joint International Conference
More informationComputer Networks More general queuing systems
Computer Networks More general queuing systems Saad Mneimneh Computer Science Hunter College of CUNY New York M/G/ Introduction We now consider a queuing system where the customer service times have a
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 information2. Transience and Recurrence
Virtual Laboratories > 15. Markov Chains > 1 2 3 4 5 6 7 8 9 10 11 12 2. Transience and Recurrence The study of Markov chains, particularly the limiting behavior, depends critically on the random times
More informationMarkov Processes and Queues
MIT 2.853/2.854 Introduction to Manufacturing Systems Markov Processes and Queues Stanley B. Gershwin Laboratory for Manufacturing and Productivity Massachusetts Institute of Technology Markov Processes
More informationModel reversibility of a two dimensional reflecting random walk and its application to queueing network
arxiv:1312.2746v2 [math.pr] 11 Dec 2013 Model reversibility of a two dimensional reflecting random walk and its application to queueing network Masahiro Kobayashi, Masakiyo Miyazawa and Hiroshi Shimizu
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 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 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 informationG-networks with synchronized partial ushing. PRi SM, Universite de Versailles, 45 av. des Etats Unis, Versailles Cedex,France
G-networks with synchronized partial ushing Jean-Michel FOURNEAU ;a, Dominique VERCH ERE a;b a PRi SM, Universite de Versailles, 45 av. des Etats Unis, 78 05 Versailles Cedex,France b CERMSEM, Universite
More information6 Markov Chain Monte Carlo (MCMC)
6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution
More informationDynamic resource sharing
J. Virtamo 38.34 Teletraffic Theory / Dynamic resource sharing and balanced fairness Dynamic resource sharing In previous lectures we have studied different notions of fair resource sharing. Our focus
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 informationInput-queued switches: Scheduling algorithms for a crossbar switch. EE 384X Packet Switch Architectures 1
Input-queued switches: Scheduling algorithms for a crossbar switch EE 84X Packet Switch Architectures Overview Today s lecture - the input-buffered switch architecture - the head-of-line blocking phenomenon
More informationBirth-Death Processes
Birth-Death Processes Birth-Death Processes: Transient Solution Poisson Process: State Distribution Poisson Process: Inter-arrival Times Dr Conor McArdle EE414 - Birth-Death Processes 1/17 Birth-Death
More informationOn Tandem Blocking Queues with a Common Retrial Queue
On Tandem Blocking Queues with a Common Retrial Queue K. Avrachenkov U. Yechiali Abstract We consider systems of tandem blocking queues having a common retrial queue, for which explicit analytic results
More informationMarkov Processes Hamid R. Rabiee
Markov Processes Hamid R. Rabiee Overview Markov Property Markov Chains Definition Stationary Property Paths in Markov Chains Classification of States Steady States in MCs. 2 Markov Property A discrete
More informationT. Liggett Mathematics 171 Final Exam June 8, 2011
T. Liggett Mathematics 171 Final Exam June 8, 2011 1. The continuous time renewal chain X t has state space S = {0, 1, 2,...} and transition rates (i.e., Q matrix) given by q(n, n 1) = δ n and q(0, n)
More informationECE 3511: Communications Networks Theory and Analysis. Fall Quarter Instructor: Prof. A. Bruce McDonald. Lecture Topic
ECE 3511: Communications Networks Theory and Analysis Fall Quarter 2002 Instructor: Prof. A. Bruce McDonald Lecture Topic Introductory Analysis of M/G/1 Queueing Systems Module Number One Steady-State
More informationA New Technique for Link Utilization Estimation
A New Technique for Link Utilization Estimation in Packet Data Networks using SNMP Variables S. Amarnath and Anurag Kumar* Dept. of Electrical Communication Engineering Indian Institute of Science, Bangalore
More informationλ, µ, ρ, A n, W n, L(t), L, L Q, w, w Q etc. These
Queuing theory models systems with servers and clients (presumably waiting in queues). Notation: there are many standard symbols like λ, µ, ρ, A n, W n, L(t), L, L Q, w, w Q etc. These represent the actual
More informationSession-Based Queueing Systems
Session-Based Queueing Systems Modelling, Simulation, and Approximation Jeroen Horters Supervisor VU: Sandjai Bhulai Executive Summary Companies often offer services that require multiple steps on the
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 informationMulti Stage Queuing Model in Level Dependent Quasi Birth Death Process
International Journal of Statistics and Systems ISSN 973-2675 Volume 12, Number 2 (217, pp. 293-31 Research India Publications http://www.ripublication.com Multi Stage Queuing Model in Level Dependent
More information57:022 Principles of Design II Final Exam Solutions - Spring 1997
57:022 Principles of Design II Final Exam Solutions - Spring 1997 Part: I II III IV V VI Total Possible Pts: 52 10 12 16 13 12 115 PART ONE Indicate "+" if True and "o" if False: + a. If a component's
More informationCDA5530: Performance Models of Computers and Networks. Chapter 4: Elementary Queuing Theory
CDA5530: Performance Models of Computers and Networks Chapter 4: Elementary Queuing Theory Definition Queuing system: a buffer (waiting room), service facility (one or more servers) a scheduling policy
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 informationCPSC 531: System Modeling and Simulation. Carey Williamson Department of Computer Science University of Calgary Fall 2017
CPSC 531: System Modeling and Simulation Carey Williamson Department of Computer Science University of Calgary Fall 2017 Motivating Quote for Queueing Models Good things come to those who wait - poet/writer
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 informationIrreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1
Irreducibility Irreducible every state can be reached from every other state For any i,j, exist an m 0, such that i,j are communicate, if the above condition is valid Irreducible: all states are communicate
More informationFlow Equivalence and Stochastic Equivalence in G-Networks
Flow Equivalence and Stochastic Equivalence in G-Networks Jean-Michel Fourneau Laboratoire PRISM Univ. de Versailles Saint-Quentin 45 Avenue des Etats-Unis 78000 Versailles, France jmf@prism.uvsq.fr Erol
More informationExample: physical systems. If the state space. Example: speech recognition. Context can be. Example: epidemics. Suppose each infected
4. Markov Chains A discrete time process {X n,n = 0,1,2,...} with discrete state space X n {0,1,2,...} is a Markov chain if it has the Markov property: P[X n+1 =j X n =i,x n 1 =i n 1,...,X 0 =i 0 ] = P[X
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 informationQueueing Theory and Simulation. Introduction
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, Japan
More informationName of the Student:
SUBJECT NAME : Probability & Queueing Theory SUBJECT CODE : MA 6453 MATERIAL NAME : Part A questions REGULATION : R2013 UPDATED ON : November 2017 (Upto N/D 2017 QP) (Scan the above QR code for the direct
More informationA Study on Performance Analysis of Queuing System with Multiple Heterogeneous Servers
UNIVERSITY OF OKLAHOMA GENERAL EXAM REPORT A Study on Performance Analysis of Queuing System with Multiple Heterogeneous Servers Prepared by HUSNU SANER NARMAN husnu@ou.edu based on the papers 1) F. S.
More informationLecture 20 : Markov Chains
CSCI 3560 Probability and Computing Instructor: Bogdan Chlebus Lecture 0 : Markov Chains We consider stochastic processes. A process represents a system that evolves through incremental changes called
More information