Applied Mathematics

Volume 8, Issue 9 (September 2017)

ISSN Print: 2152-7385   ISSN Online: 2152-7393

Google-based Impact Factor: 0.58  Citations  

Inherent Numerical Instability in Computing Invariant Measures of Markov Chains

HTML  XML Download Download as PDF (Size: 361KB)  PP. 1367-1385  
DOI: 10.4236/am.2017.89101    1,099 Downloads   3,465 Views  Citations

ABSTRACT

Invariant measures of Markov chains in discrete or continuous time with a countable set of states are characterized by its steady state recurrence relations. Exemplarily, we consider transition matrices and Q-matrices with upper bandwidth n and lower bandwidth 1 where the invariant measures satisfy an (n + 1)-order linear difference equation. Markov chains of this type arise from applications to queueing problems and population dynamics. It is the purpose of this paper to point out that the forward use of this difference equation is subject to some hitherto unobserved aspects. By means of the concept of generalized continued fractions (GCFs), we prove that each invariant measure is a dominated solution of the difference equation such that forward computation becomes numerically unstable. Furthermore, the GCF-based approach provides a decoupled recursion in which the phenomenon of numerical instability does not appear. The procedure results in an iteration scheme for successively computing approximants of the desired invariant measure depending on some truncation level N. Increasing N leads to the desired solution. A comparison study of forward computation and the GCF-based approach is given for Q-matrices with upper bandwidth 1 and 2.

Share and Cite:

Baumann, H. and Hanschke, T. (2017) Inherent Numerical Instability in Computing Invariant Measures of Markov Chains. Applied Mathematics, 8, 1367-1385. doi: 10.4236/am.2017.89101.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.