Matrix Geometric Approach for Random Walks

Total Page:16

File Type:pdf, Size:1020Kb

Matrix Geometric Approach for Random Walks Matrix geometric approach for random walks Citation for published version (APA): Kapodistria, S., & Palmowski, Z. B. (2017). Matrix geometric approach for random walks: stability condition and equilibrium distribution. Stochastic Models, 33(4), 572-597. https://doi.org/10.1080/15326349.2017.1359096 Document license: CC BY-NC-ND DOI: 10.1080/15326349.2017.1359096 Document status and date: Published: 02/10/2017 Document Version: Publisher’s PDF, also known as Version of Record (includes final page, issue and volume numbers) Please check the document version of this publication: • A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website. • The final author version and the galley proof are versions of the publication after peer review. • The final published version features the final layout of the paper including the volume, issue and page numbers. Link to publication General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal. If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement: www.tue.nl/taverne Take down policy If you believe that this document breaches copyright please contact us at: [email protected] providing details and we will investigate your claim. Download date: 23. Sep. 2021 Stochastic Models ISSN: 1532-6349 (Print) 1532-4214 (Online) Journal homepage: http://www.tandfonline.com/loi/lstm20 Matrix geometric approach for random walks: Stability condition and equilibrium distribution Stella Kapodistria & Zbigniew Palmowski To cite this article: Stella Kapodistria & Zbigniew Palmowski (2017) Matrix geometric approach for random walks: Stability condition and equilibrium distribution, Stochastic Models, 33:4, 572-597, DOI: 10.1080/15326349.2017.1359096 To link to this article: https://doi.org/10.1080/15326349.2017.1359096 © 2017 The Author(s). Published with license by Taylor & Francis© Stella Kapodistria and Zbigniew Palmowski Published online: 31 Aug 2017. Submit your article to this journal Article views: 145 View related articles View Crossmark data Full Terms & Conditions of access and use can be found at http://www.tandfonline.com/action/journalInformation?journalCode=lstm20 Download by: [Eindhoven University of Technology] Date: 13 December 2017, At: 06:14 STOCHASTIC MODELS , VOL. , NO. , – https://doi.org/./.. Matrix geometric approach for random walks: Stability condition and equilibrium distribution Stella Kapodistriaa and Zbigniew Palmowskib aDepartment of Mathematics and Computer Science, Eindhoven University of Technology, Eindhoven, The Netherlands; bFaculty of Pure and Applied Mathematics, Wroclaw University of Science and Technology, Wroclaw, Poland ABSTRACT ARTICLE HISTORY In this paper, we analyze a sub-class of two-dimensional homo- Received April geneous nearest neighbor (simple) random walk restricted on the Accepted July lattice using the matrix geometric approach. In particular, we first KEYWORDS present an alternative approach for the calculation of the stability [30] Boundary value problem condition, extending the result of Neuts drift conditions and method; compensation connecting it with the result of Fayolle et al. which is based on approach; equilibrium Lyapunov functions.[13] Furthermore, we consider the sub-class distribution; matrix of random walks with equilibrium distributions given as series of geometric approach; random product forms and, for this class of random walks, we calculate walks; spectrum; stability the eigenvalues and the corresponding eigenvectors of the infi- condition R nite matrix appearing in the matrix geometric approach. This MATHEMATICS SUBJECT result is obtained by connecting and extending three existing CLASSIFICATION approaches available for such an analysis: the matrix geometric K; B approach, the compensation approach and the boundary value problem method. In this paper, we also present the spectral prop- erties of the infinite matrix R. 1. Introduction The objective of this work is to demonstrate how to obtain the stability condition and the equilibrium distribution of the state of a two-dimensional homogeneous nearest neighbor (simple) random walk restricted on the lattice using its underly- ing Quasi-Birth–Death (QBD) structure and the matrix geometric approach. This type of random walk can be modeled as a QBD process with the characteris- tic that both the levels and the phases are countably infinite. Then, based on the Downloaded by [Eindhoven University of Technology] at 06:14 13 December 2017 matrix geometric approach, if πn = (πn,0 πn,1 ···) denotes the vector of the equi- librium distribution at level n, n = 0, 1,...,itisknownthatπn+1 = πnR, n ≥ 1. Thisisaverywell-knownresult,butthecomplexityofthesolutionliesinthe calculation of the infinite-dimension matrix R. In this paper, we investigate how the matrix geometric approach can be extended to the case of countably infinite CONTACT Stella Kapodistria [email protected] Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box , MB Eindhoven, The Netherlands. Color versions of one or more of the figures in the article can be found online at www.tandfonline.com/lstm. Published with license by Taylor & Francis © Stella Kapodistria and Zbigniew Palmowski. This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License (http://creativecommons.org/licenses/by-nc-nd/./), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited, and is not altered, transformed, or built upon in any way. STOCHASTIC MODELS 573 phases and we discuss the challenges that arise with such an extension. Further- more, we propose an approach for the calculation of the eigenvalues and eigenvec- tors of matrix R, that complements the existing approaches for this type of random walks. Moreover, this approach can be numerically used for the approximation of the matrix R by considering spectral truncation instead of the usual state space truncation. 1.1. Literature overview ThemainbodyofliteratureonthetopicofcountablyinfinitephasesofQBDsis devoted to either stochastic processes in which the matrix R has a simple structure or to asymptotic results concerning the decay rate. Regarding QBDs with a simple structure, the literature is mainly devoted to stochastic processes for which the rate matrix has a simple property, see, e.g.,[11] and the references therein, or to stochastic process with a special structure in the allowed transitions, see, e.g.,[19,20] and the references therein. In Ref.[11],theauthors consider an infinite rate matrix R that can be written as the product of a vector col- umn times a row vector. This simple structure permits, under the assumptions of homogeneity and irreducibility, to show that the stationary distributions of these processes have a product form structure as a function of the level. Additionally, they apply their results in the case of the Cox(k)/MY /1 queue. In Refs.[19,20] and the references therein, the authors investigate stochastic processes with a special struc- ture regarding the allowed transitions satisfying the so-called successive lumpability property. This class of Markov chains is specified through its property of cal- culating the stationary probabilities of the system by successively computing the stationary probabilities of a propitiously constructed sequence of Markov chains. Each of the latter chains has a (typically much) smaller state space and this yields significant computational improvements. Such models are a special class of the GI/M/1 or the M/G/1 type of QBD processes. In Refs.[19,20],theauthorsdiscuss how the results for discrete-time Markov chains extend to semi-Markov processes and continuous-time Markov processes. Furthermore, in Ref.[20], the authors com- pare the successive lumping methodology developed with the lattice path count- ing approach[28] for the calculation of the rate matrix of a queueing model. These two methodologies are compared both in terms of applicability requirements and Downloaded by [Eindhoven University of Technology] at 06:14 13 December 2017 numerical complexity by analyzing their performance for some classical queueing models. Their main findings are (i) when both methods are applicable, the suc- cessive lumping-based algorithms outperform the lattice path counting algorithm, (ii) the successive lumping algorithms, contrary to the lattice path counting algo- rithm, include a method to compute the steady-state distribution using this rate matrix. Regarding the work devoted on the decay rate of QBDs with
Recommended publications
  • The Queueing Network Analyzer
    THE BELL SYSTEM TECHNICAL JOURNAL Vol. 62, No.9, November 1983 Printed in U.S.A. The Queueing Network Analyzer By W. WHITT* (Manuscript received March 11, 1983) This paper describes the Queueing Network Analyzer (QNA), a software package developed at Bell Laboratories to calculate approximate congestion measures for a network of queues. The first version of QNA analyzes open networks of multiserver nodes with the first-come, first-served discipline and no capacity constraints. An important feature is that the external arrival processes need not be Poisson and the service-time distributions need not be exponential. Treating other kinds of variability is important. For example, with packet-switchedcommunication networks we need to describe the conges­ tion resulting from bursty traffic and the nearly constant service times of packets. The general approach in QNA is to approximately characterize the arrival processes by two or three parameters and then analyze the individual nodes separately. The first version of QNA uses two parameters to characterize the arrival processes and service times, one to describe the rate and the other to describe the variability. The nodes are then analyzed as standard GI/G/m queues partially characterized by the first two moments of the interarrival­ time and service-time distributions. Congestion measures for the network as a whole are obtained by assuming as an approximation that the nodes are stochastically independent given the approximate flow parameters. I. INTRODUCTION AND SUMMARY Networks of queues have proven to be useful models to analyze the performance of complex systems such as computers, switching ma­ chines, communications networks, and production job shOpS.1-7 To facilitate the analysis of these models, several software packages have * Bell Laboratories.
    [Show full text]
  • Product-Form in Queueing Networks
    Product-form in queueing networks VRIJE UNIVERSITEIT Product-form in queueing networks ACADEMISCH PROEFSCHRIFT ter verkrijging van de graad van doctor aan de Vrije Universiteit te Amsterdam, op gezag van de rector magnificus dr. C. Datema, hoogleraar aan de faculteit der letteren, in het openbaar te verdedigen ten overstaan van de promotiecommissie van de faculteit der economische wetenschappen en econometrie op donderdag 21 mei 1992 te 15.30 uur in het hoofdgebouw van de universiteit, De Boelelaan 1105 door Richardus Johannes Boucherie geboren te Oost- en West-Souburg Thesis Publishers Amsterdam 1992 Promotoren: prof.dr. N.M. van Dijk prof.dr. H.C. Tijms Referenten: prof.dr. A. Hordijk prof.dr. P. Whittle Preface This monograph studies product-form distributions for queueing networks. The celebrated product-form distribution is a closed-form expression, that is an analytical formula, for the queue-length distribution at the stations of a queueing network. Based on this product-form distribution various so- lution techniques for queueing networks can be derived. For example, ag- gregation and decomposition results for product-form queueing networks yield Norton's theorem for queueing networks, and the arrival theorem implies the validity of mean value analysis for product-form queueing net- works. This monograph aims to characterize the class of queueing net- works that possess a product-form queue-length distribution. To this end, the transient behaviour of the queue-length distribution is discussed in Chapters 3 and 4, then in Chapters 5, 6 and 7 the equilibrium behaviour of the queue-length distribution is studied under the assumption that in each transition a single customer is allowed to route among the stations only, and finally, in Chapters 8, 9 and 10 the assumption that a single cus- tomer is allowed to route in a transition only is relaxed to allow customers to route in batches.
    [Show full text]
  • BUSF 40901-1/CMSC 34901-1: Stochastic Performance Modeling Winter 2014
    BUSF 40901-1/CMSC 34901-1: Stochastic Performance Modeling Winter 2014 Syllabus (January 15, 2014) Instructor: Varun Gupta Office: 331 Harper Center e-mail: [email protected] Phone: 773-702-7315 Office hours: by appointment Class Times: Wed, Fri { 10:10-11:30 am { Harper Center (3A) Final Exam (tentative): March 21, Friday { 8:00-11:00am { Harper Center (3A) Course Website: http://chalk.uchicago.edu Course Objectives This is an introductory course in queueing theory and performance modeling, with applications including but not limited to service operations (healthcare, call centers) and computer system resource management (from datacenter to kernel level). The aim of the course is two-fold: 1. Build insights into best practices for designing service systems (How many service stations should I provision? What speed? How should I separate/prioritize customers based on their service requirements?) 2. Build a basic toolbox for analyzing queueing systems in particular and stochastic processes in general. Tentative list of topics: Open/closed queueing networks; Operational laws; M=M=1 queue; Burke's theorem and reversibility; M=M=k queue; M=G=1 queue; G=M=1 queue; P h=P h=k queues and their solution using matrix-analytic methods; Arrival theorem and Mean Value Analysis; Analysis of scheduling policies (e.g., Last-Come-First Served; Processor Sharing); Jackson network and the BCMP theorem (product form networks); Asymptotic analysis (M=M=k queue in heavy/light traf- fic, Supermarket model in mean-field regime) Prerequisites Exposure to undergraduate probability (random variables, discrete and continuous probability dis- tributions, discrete time Markov chains) and calculus is required.
    [Show full text]
  • Inferring Queueing Network Models from High-Precision Location Tracking Data
    University of London Imperial College of Science, Technology and Medicine Department of Computing Inferring Queueing Network Models from High-precision Location Tracking Data Tzu-Ching Horng Submitted in part fulfilment of the requirements for the degree of Doctor of Philosophy in Computing of the University of London and the Diploma of Imperial College, May 2013 Abstract Stochastic performance models are widely used to analyse the performance and reliability of systems that involve the flow and processing of customers. However, traditional methods of constructing a performance model are typically manual, time-consuming, intrusive and labour- intensive. The limited amount and low quality of manually-collected data often lead to an inaccurate picture of customer flows and poor estimates of model parameters. Driven by ad- vances in wireless sensor technologies, recent real-time location systems (RTLSs) enable the automatic, continuous and unintrusive collection of high-precision location tracking data, in both indoor and outdoor environment. This high-quality data provides an ideal basis for the construction of high-fidelity performance models. This thesis presents a four-stage data processing pipeline which takes as input high-precision location tracking data and automatically constructs a queueing network performance model approximating the underlying system. The first two stages transform raw location traces into high-level “event logs” recording when and for how long a customer entity requests service from a server entity. The third stage infers the customer flow structure and extracts samples of time delays involved in the system; including service time, customer interarrival time and customer travelling time. The fourth stage parameterises the service process and customer arrival process of the final output queueing network model.
    [Show full text]
  • Markovian Queueing Networks
    Richard J. Boucherie Markovian queueing networks Lecture notes LNMB course MQSN September 5, 2020 Springer Contents Part I Solution concepts for Markovian networks of queues 1 Preliminaries .................................................. 3 1.1 Basic results for Markov chains . .3 1.2 Three solution concepts . 11 1.2.1 Reversibility . 12 1.2.2 Partial balance . 13 1.2.3 Kelly’s lemma . 13 2 Reversibility, Poisson flows and feedforward networks. 15 2.1 The birth-death process. 15 2.2 Detailed balance . 18 2.3 Erlang loss networks . 21 2.4 Reversibility . 23 2.5 Burke’s theorem and feedforward networks of MjMj1 queues . 25 2.6 Literature . 28 3 Partial balance and networks with Markovian routing . 29 3.1 Networks of MjMj1 queues . 29 3.2 Kelly-Whittle networks. 35 3.3 Partial balance . 39 3.4 State-dependent routing and blocking protocols . 44 3.5 Literature . 50 4 Kelly’s lemma and networks with fixed routes ..................... 51 4.1 The time-reversed process and Kelly’s Lemma . 51 4.2 Queue disciplines . 53 4.3 Networks with customer types and fixed routes . 59 4.4 Quasi-reversibility . 62 4.5 Networks of quasi-reversible queues with fixed routes . 68 4.6 Literature . 70 v Part I Solution concepts for Markovian networks of queues Chapter 1 Preliminaries This chapter reviews and discusses the basic assumptions and techniques that will be used in this monograph. Proofs of results given in this chapter are omitted, but can be found in standard textbooks on Markov chains and queueing theory, e.g. [?, ?, ?, ?, ?, ?, ?]. Results from these references are used in this chapter without reference except for cases where a specific result (e.g.
    [Show full text]
  • Stabilization of Branching Queueing Networks∗
    Stabilization of Branching Queueing Networks∗ Tomáš Brázdil1 and Stefan Kiefer2 1 Faculty of Informatics, Masaryk University, Czech Republic [email protected] 2 Department of Computer Science, University of Oxford, United Kingdom [email protected] Abstract Queueing networks are gaining attraction for the performance analysis of parallel computer sys- tems. A Jackson network is a set of interconnected servers, where the completion of a job at server i may result in the creation of a new job for server j. We propose to extend Jackson networks by “branching” and by “control” features. Both extensions are new and substantially expand the modelling power of Jackson networks. On the other hand, the extensions raise com- putational questions, particularly concerning the stability of the networks, i.e, the ergodicity of the underlying Markov chain. We show for our extended model that it is decidable in polyno- mial time if there exists a controller that achieves stability. Moreover, if such a controller exists, one can efficiently compute a static randomized controller which stabilizes the network in a very strong sense; in particular, all moments of the queue sizes are finite. 1998 ACM Subject Classification G.3 Probability and Statistics Keywords and phrases continuous-time Markov decision processes, infinite-state systems, per- formance analysis 1 Introduction Queueing theory plays a central role in the performance analysis of computer systems. In particular, queueing networks are gaining attraction as models of parallel systems. A queueing network is a set of processing units (called servers), each of which performs tasks (called jobs) of a certain type.
    [Show full text]
  • Introduction to Queueing Theory Review on Poisson Process
    Contents ELL 785–Computer Communication Networks Motivations Lecture 3 Discrete-time Markov processes Introduction to Queueing theory Review on Poisson process Continuous-time Markov processes Queueing systems 3-1 3-2 Circuit switching networks - I Circuit switching networks - II Traffic fluctuates as calls initiated & terminated Fluctuation in Trunk Occupancy Telephone calls come and go • Number of busy trunks People activity follow patterns: Mid-morning & mid-afternoon at All trunks busy, new call requests blocked • office, Evening at home, Summer vacation, etc. Outlier Days are extra busy (Mother’s Day, Christmas, ...), • disasters & other events cause surges in traffic Providing resources so Call requests always met is too expensive 1 active • Call requests met most of the time cost-effective 2 active • 3 active Switches concentrate traffic onto shared trunks: blocking of requests 4 active active will occur from time to time 5 active Trunk number Trunk 6 active active 7 active active Many Fewer lines trunks – minimize the number of trunks subject to a blocking probability 3-3 3-4 Packet switching networks - I Packet switching networks - II Statistical multiplexing Fluctuations in Packets in the System Dedicated lines involve not waiting for other users, but lines are • used inefficiently when user traffic is bursty (a) Dedicated lines A1 A2 Shared lines concentrate packets into shared line; packets buffered • (delayed) when line is not immediately available B1 B2 C1 C2 (a) Dedicated lines A1 A2 B1 B2 (b) Shared line A1 C1 B1 A2 B2 C2 C1 C2 A (b)
    [Show full text]
  • Elements of Probability and Queueing Theory
    APPENDIX A Elements of Probability and Queueing Theory The primary purpose of this appendix is to summarize for the reader's conve­ nience those parts of the elements of probability theory, Markov process, and queueing theory that are most relevant to the subject matter of and in consistent notations with this book. Mathematical terminologies are used only for the pur­ pose of introducing the language to the intended readers. No pretense at rigor is attempted. Except for some relatively new material and interpretations, the reader can find a more thorough introduction of this material in [Kleinrock 1975 Chapters 2 - 4] and [Cinlar 1975]. A.1 Elements of Probability Theory Let (Q, F, P) be a probability space, where F is a O'-field in Q and P is a prob- ability measure defined on F. A random variable X is areal (measurable) function defined on the probability space. For any roe Q, X(ro)e R is areal number. The distribution function of X is defined as F(x) = P(ro: X(ro) ~x). Function fex) is called the probability density of the random variable X if for any subset Ae R, we have P( Xe A) =i f(x) dx. If F(x) is differentiable at x, then fex) = dF(x)/dx. If X takes values only in a countable set of real numbers, X and F(x) are called discrete. A familiar dis­ crete distribution is the binomial distribution: 312 APPENDIX A r=O,I,2, ... ,n. If n=l, the distribution is called a Bernoulli distribution. For engineering pur­ poses, we can think of F(x) or f(x) as a function defining or characterizing the random variable, X.
    [Show full text]
  • Queueing Networks 1 QUEUEING NETWORKS
    J. Virtamo 38.3143 Queueing Theory / Queueing networks 1 QUEUEING NETWORKS A network consisting of several interconnected queues Network of queues • Examples Customers go form one queue to another in post office, bank, supermarket etc • Data packets traverse a network moving from a queue in a router to the queue in another • router History Burke’s theorem, Burke (1957), Reich (1957) • Jackson (1957, 1963): open queueing networks, product form solution • Gordon and Newell (1967): closed queueing networks • Baskett, Chandy, Muntz, Palacios (1975): generalizations of the types of queues • Reiser and Lavenberg (1980, 1982): mean value analysis, MVA • J. Virtamo 38.3143 Queueing Theory / Queueing networks 2 Jackson’s queueing network (open queueing network) Jackson’s open queueing network consists of M nodes (queues) with the following assumptions: Node i is a FIFO queue • – unlimited number of waiting places (infinite queue) Service time in the queue obeys the distribution Exp(µ ) • i – in each queue, the service time of the customer is drawn independent of the service times in other queues – note: in a packet network the sending time of a packet, in reality, is the same in all queues (or differs by a constant factor, the inverse of the line speed) – this dependence, however, does not markedly affect the behaviour of the system (so called Kleinrock’s independence assumption) Upon departure from queue i, the customer chooses the next queue j randomly with the • probability qi,j or exits the network with the probability qi,d (probabilistic routing) – the model can be extended to cover the case of predetermined routes (route pinning) The network is open to arrivals from outside of the network (source) • – from the source s customers arrive as a Poisson stream with intensity λ – fraction qs,i of them enter queue i (intensity λqs,i) J.
    [Show full text]
  • Stabilization of Branching Queueing Networks Brázdil Tomáš, Stefan Kiefer
    Stabilization of Branching Queueing Networks Brázdil Tomáš, Stefan Kiefer To cite this version: Brázdil Tomáš, Stefan Kiefer. Stabilization of Branching Queueing Networks. STACS’12 (29th Sym- posium on Theoretical Aspects of Computer Science), Feb 2012, Paris, France. pp.507-518. hal- 00678184 HAL Id: hal-00678184 https://hal.archives-ouvertes.fr/hal-00678184 Submitted on 3 Feb 2012 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Stabilization of Branching Queueing Networks∗ Tomáš Brázdil1 and Stefan Kiefer2 1 Faculty of Informatics, Masaryk University, Czech Republic [email protected] 2 Department of Computer Science, University of Oxford, United Kingdom [email protected] Abstract Queueing networks are gaining attraction for the performance analysis of parallel computer sys- tems. A Jackson network is a set of interconnected servers, where the completion of a job at server i may result in the creation of a new job for server j. We propose to extend Jackson networks by “branching” and by “control” features. Both extensions are new and substantially expand the modelling power of Jackson networks. On the other hand, the extensions raise com- putational questions, particularly concerning the stability of the networks, i.e, the ergodicity of the underlying Markov chain.
    [Show full text]
  • Studies in Heavy Traffic and in Production Systems
    Iowa State University Capstones, Theses and Retrospective Theses and Dissertations Dissertations 1999 Studies in heavy trafficnd a in production systems Zugeng Zheng Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/rtd Part of the Industrial Engineering Commons, Operational Research Commons, and the Statistics and Probability Commons Recommended Citation Zheng, Zugeng, "Studies in heavy trafficnd a in production systems " (1999). Retrospective Theses and Dissertations. 12499. https://lib.dr.iastate.edu/rtd/12499 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Retrospective Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, beginning at the upper left-hand comer and continuing from left to right in equal sections with small overiaps.
    [Show full text]
  • Networks of Queues
    Networks of queues c University of Bristol, 2012 1 Time reversal of Markov chains In studying the M=M=1 queue, we made use of reversibility in order to compute its invariant distribution. While this could have been done, in principle, by solving the global balance equations, reversibility made it much easier as we only needed to solve the local balance equations. This becomes even more important when trying to obtain the invariant distribution of Markov chains whose state space is multi-dimensional, as is the case for networks of queues. (For an M=M=1 queue, the state space is Z+, the set of non-negative integers, whereas for a network of d such queues, it would be d Z+.) In general, there is no good way of obtaining the invariant distribution unless the Markov chain is reversible. Unfortunately, not all Markov chains of interest are reversible! However, there is a class of queueing networks called Jackson networks which, while they aren't reversible, can be analysed using techniques related to time reversal. So we'll start with a study of those ideas, leading up to the development of Kelly's lemma. Let (Xt; t 2 R) be a continuous time Markov chain on some countable state space S, with transition rate matrix Q. Recall that any probability vector π on the state space S which solves the global balanace equations πQ = 0 P P (or, writing it out componentwise, j6=k πjqjk = −πkqkk = πk j6=k qkj for all k 2 S) is an invariant distribution of this Markov chain.
    [Show full text]