Queueing Theory II. Summary. ! M/M/1 Output process. ! Networks of Queue! Method of Stages. ! General Distributions

Similar documents
Queueing Theory II. Summary. M/M/1 Output process Networks of Queue Method of Stages. General Distributions

Contents Preface The Exponential Distribution and the Poisson Process Introduction to Renewal Theory

Queuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem. Wade Trappe

P (L d k = n). P (L(t) = n),

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 "

Link Models for Circuit Switching

Non Markovian Queues (contd.)

GI/M/1 and GI/M/m queuing systems

J. MEDHI STOCHASTIC MODELS IN QUEUEING THEORY

6 Solving Queueing Models

Computer Networks More general queuing systems

Performance Evaluation of Queuing Systems

Queueing Review. Christos Alexopoulos and Dave Goldsman 10/6/16. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA

Queues and Queueing Networks

M/G/1 and M/G/1/K systems

Introduction to Markov Chains, Queuing Theory, and Network Performance

Queueing Review. Christos Alexopoulos and Dave Goldsman 10/25/17. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA

Name of the Student:

CPSC 531: System Modeling and Simulation. Carey Williamson Department of Computer Science University of Calgary Fall 2017

Class 11 Non-Parametric Models of a Service System; GI/GI/1, GI/GI/n: Exact & Approximate Analysis.

Introduction to Queueing Theory

Intro to Queueing Theory

Lecture 7: Simulation of Markov Processes. Pasi Lassila Department of Communications and Networking

INDEX. production, see Applications, manufacturing

Introduction to Queueing Theory

Kendall notation. PASTA theorem Basics of M/M/1 queue

Queueing Theory and Simulation. Introduction

Elementary queueing system

Solutions to Homework Discrete Stochastic Processes MIT, Spring 2011

HITTING TIME IN AN ERLANG LOSS SYSTEM

NEW FRONTIERS IN APPLIED PROBABILITY

Chapter 5: Special Types of Queuing Models

Advanced Queueing Theory

Chapter 2 Queueing Theory and Simulation

ECE 3511: Communications Networks Theory and Analysis. Fall Quarter Instructor: Prof. A. Bruce McDonald. Lecture Topic

Probability Models in Electrical and Computer Engineering Mathematical models as tools in analysis and design Deterministic models Probability models

Introduction to Queueing Theory

QUEUING SYSTEM. Yetunde Folajimi, PhD

Intro Refresher Reversibility Open networks Closed networks Multiclass networks Other networks. Queuing Networks. Florence Perronnin

Contents LIST OF TABLES... LIST OF FIGURES... xvii. LIST OF LISTINGS... xxi PREFACE. ...xxiii

THE KINLEITH WEIGHBRIDGE

Data analysis and stochastic modeling

Markov 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

Part II: continuous time Markov chain (CTMC)

Statistics 150: Spring 2007

Queueing systems. Renato Lo Cigno. Simulation and Performance Evaluation Queueing systems - Renato Lo Cigno 1

CDA5530: Performance Models of Computers and Networks. Chapter 4: Elementary Queuing Theory

A Queueing System with Queue Length Dependent Service Times, with Applications to Cell Discarding in ATM Networks

Time Reversibility and Burke s Theorem

15 Closed production networks

UNIVERSITY OF YORK. MSc Examinations 2004 MATHEMATICS Networks. Time Allowed: 3 hours.

Part I Stochastic variables and Markov chains

1 Basic concepts from probability theory

2905 Queueing Theory and Simulation PART III: HIGHER DIMENSIONAL AND NON-MARKOVIAN QUEUES

Stationary remaining service time conditional on queue length

15 Closed production networks

SYMBOLS AND ABBREVIATIONS

Queueing Networks G. Rubino INRIA / IRISA, Rennes, France

Page 0 of 5 Final Examination Name. Closed book. 120 minutes. Cover page plus five pages of exam.

Exercises Stochastic Performance Modelling. Hamilton Institute, Summer 2010

Queueing Systems: Lecture 3. Amedeo R. Odoni October 18, Announcements

EP2200 Queueing theory and teletraffic systems

LIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974

Introduction to Queueing Theory with Applications to Air Transportation Systems

Stochastic Models of Manufacturing Systems

Introduction to queuing theory

Bulk input queue M [X] /M/1 Bulk service queue M/M [Y] /1 Erlangian queue M/E k /1

Lecture 20: Reversible Processes and Queues

1 IEOR 4701: Continuous-Time Markov Chains

Classical Queueing Models.

Solutions to COMP9334 Week 8 Sample Problems

Queueing Networks and Insensitivity

Operations Research II, IEOR161 University of California, Berkeley Spring 2007 Final Exam. Name: Student ID:

Glossary availability cellular manufacturing closed queueing network coefficient of variation (CV) conditional probability CONWIP

Chapter 8 Queuing Theory Roanna Gee. W = average number of time a customer spends in the system.

Stochastic process. X, a series of random variables indexed by t

Simple queueing models

Introduction to Queuing Theory. Mathematical Modelling

TCOM 501: Networking Theory & Fundamentals. Lecture 6 February 19, 2003 Prof. Yannis A. Korilis

Introduction to Queuing Networks Solutions to Problem Sheet 3

Exact Simulation of the Stationary Distribution of M/G/c Queues

Figure 10.1: Recording when the event E occurs

M/G/1 and Priority Queueing

Queuing Theory. Richard Lockhart. Simon Fraser University. STAT 870 Summer 2011

Lecturer: Olga Galinina

CS 798: Homework Assignment 3 (Queueing Theory)

reversed chain is ergodic and has the same equilibrium probabilities (check that π j =

NICTA Short Course. Network Analysis. Vijay Sivaraman. Day 1 Queueing Systems and Markov Chains. Network Analysis, 2008s2 1-1

IEOR 6711, HMWK 5, Professor Sigman

M/G/1 queues and Busy Cycle Analysis

Probability and Statistics Concepts

Buzen s algorithm. Cyclic network Extension of Jackson networks

Jitter Analysis of an MMPP 2 Tagged Stream in the presence of an MMPP 2 Background Stream

Chapter 10. Queuing Systems. D (Queuing Theory) Queuing theory is the branch of operations research concerned with waiting lines.

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science

ECE-517: Reinforcement Learning in Artificial Intelligence. Lecture 4: Discrete-Time Markov Chains

Continuous-Time Markov Chain

Control of Fork-Join Networks in Heavy-Traffic

IEOR 4106: Spring Solutions to Homework Assignment 7: Due on Tuesday, March 22.

Stochastic Models in Computer Science A Tutorial

Transcription:

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 Process! Burke s Theorem:! The Departure process of a stable M/M/1 queueing system with arrival rate! is also a Poisson process with rate!. A k D k-1 I k Z k Z k D k D D k-1 k Y k Y k Interdeparture Time Burke s Theorem PASTA PASTA

Two Queues in Series! µ 1 µ 2! Let the state of this system be (X 1,X 2 ) where X i is the number of customers in queue i.! 02 12 µ 2 µ 2 µ 1 µ 1!! 01 11 21 µ 2 µ 2 µ 2 µ 1 µ 1 µ 1!!! 00 10 20 30 Two Queues in Series! Balance Equations 02 12 µ 2 µ 2! µ 1! µ 1 01 11 21 µ 2 µ 2 µ 2 µ 1 µ 1 µ 1!!! 00 10 20 30!

Product Form Solution! Let " 1 =!/µ 1, " 2 =!/µ 2! Recall the M/M/1 Results! Therefore the two queues can be decoupled and studied in isolation. Jackson Networks! The product form decomposition holds for all open queueing networks with Poisson input processes that do not include feedback! Even though customer feedback causes the total input process (external Poisson and feedback) become non- Poisson, the product form solution still holds!! These types on networks are referred to as Jackson Networks! The total input rate to each note is given by Aggregate Rate External Arrival Rate Routing Prob

Closed Networks! The product form decomposition holds for also for closed networks Aggregate Rate Routing Prob Non-Poisson Processes! Assume that the service time distribution can be approximated by the sum of m iid exponential random variables with rate mµ. That is! You can show that the density of Z is given by! Z is an Erlang random variable with parameters (m, µ).

Erlang Distribution! You can show that the distribution function of Z! The expected value of Z is given by! The variance of Z is given by! Note that the variance of an Erlang is always less than or equal to the variance of the exponential random variable M/Er m /1 Queueing System! Meaning: Poisson Arrivals, Erlang distributed service times (order m), single server and infinite capacity buffer. mµ mµ mµ Stage 1 Stage 2 Stage m Erlang Server! Only a single customer is allowed in the server at any given time.! The customer has to go through all m stages before it is released.

M/Er m /1 Queueing System! The state of the system needs to take into account the stage that the customer is in, so one could use (x,s) where x=0,1,2 is the number of customers and s=1,,m is the stage that the customer in service is currently in.! Alternatively, one can use a single variable y to be the total number of stages that need to be completed before the system empties.! Let m=2!!!!!!! 0 1 2 3 4 n-1 n 2µ 2µ 2µ 2µ 2µ 2µ 2µ 2µ Er m /!/1 Queueing System! Similar to the M/Er m /1 System we can let the state be the number of arrival stages in the system, so for example, for m=2 2! 2! 2! 2! 2! 2! 2! 0 1 2 3 4 n-1 n µ µ µ µ µ µ 2! µ

Hyperexponential Distribution! Recall that the variance of the Erlang distribution is always less than or equal than the variance of the exponential distribution.! What if we need service times with higher variance? p 1 p 2 µ 1 µ 2 p m µ m Hyperexponential Server Hyperexponential Distribution! Expected values if the Hyperexponential distribution! The variance of Z is given by

M/H m /1 Queueing System! In this case, the state of the system should include both, the number of customers in the system and the stage that the customer in service is in. For example, for m=2!!!!!p 0 11 µ 1 21 31 n1 pµ 1 pµ 1 pµ 1 pµ 1 (1-p)µ 1 (1-p)µ 1 (1-p)µ 1 µ 2 pµ 2 pµ 2 pµ 2!(1-p) 12 (1-p)µ 2 (1-p)µ 2 (1-p)µ 2 22 32 n2!!!! (1-p)µ 2 M/G/1 Queueing System! When the arrival and service processes do not possess the memoryless property, we are back to the GSMP framework.! In general, we will need to keep track of the age or residual lifetime of each event.! Embedded Markov Chains " Some times it may be possible to identify specific points such that the Markov property holds " For example, for the M/G/1 system, suppose that we choose to observe the state of the system exactly after each customer departure. " For this problem, it doesn t matter how long ago the previous arrival occurred since arrivals are generated from Poisson processes. Furthermore, at the point of a departure, we know that the age of the event is always 0.

Embedded Markov Chains! So, right after each departure we observe the following chain (only the transitions from state n are drawn)! Let a j be the probability that j arrivals will occur during the interval Y defined by two consecutive departures a 1 a k-n+1 n-2 n-1 n n+1 n+2 a 2 a 3 k a 0! Let N be a random variable that indicates the number of arrivals during Y, then a j = Pr{N = j}. Embedded Markov Chains! Suppose we are given the density of Y, f Y (y) then! Since we have a Poisson process! Therefore

State Iteration! Let X k be the state of the system just after the departure of the kth customer A k X k-1 X k X k-1 X k D k-1 I k Z k Z k D k D D k-1 k Y k Y k Arrivals During Y k M/G/1 Queueing System! Pollaczek-Khinchin (PK) Formula where " X is the number of customers in the system " 1/µ is the average service time " # 2 is the variance of the service time distribution "! is the Poisson arrival rate " "=!/µ is the traffic intensity

M/G/1 Example! Consider the queueing systems M/M/1 and M/D/1! Compare their average number in the system and average system delay when the arrival is Poisson with rate! and the mean service time is 1/µ.! For the M/M/1 system, # 2 =1/µ 2, therefore! For the M/D/1 system, # 2 =0, therefore