Fast simulation of Markov chains with small transition probabilities

被引:32
作者
Juneja, S [1 ]
Shahabuddin, P
机构
[1] Indian Inst Technol, New Delhi 110016, India
[2] Columbia Univ, New York, NY 10027 USA
关键词
simulation; Markov chains; reliability models; steady-state measures; importance-sampling;
D O I
10.1287/mnsc.47.4.547.9827
中图分类号
C93 [管理学];
学科分类号
12 ; 1201 ; 1202 ; 120202 ;
摘要
Consider: a-finite-state Markov chain where the transition probabilities differ by orders of magnitude. This Markov chain has an "attractor state," i.e., from any state of the Markov chain there exists:a sample path of significant probability-to the-attractor state. There also exists a "rare set," which is accessible from the attractor state only by sample paths of very. small probability. The problem is to estimate the probability that starting from the attractor state, the Maykov-chain hits the rare set before returning to the attractor state. Examples of this setting arise in the case of reliability models with highly reliable components as well as in the case of-queueing networks with low traffic. Importance-sampling is a commonly used simulation technique for the fast estimation of rare-event probabilities. It-involves simulating the Markov chain under a new probability-measure that emphasizes the most likely paths to therare set.-Previous research focused on developing importance-sampling schemes fora special case of Markov chains that did not include "high-probability cycles." We show, through examples that the Markov chains,used to model many commonly encountered systems do have high-probability cycles, and existing importance-sampling schemes can lead to infinite variance in simulating such systems, We then develop the insight that in the:presence. of high-probability cycles care should be taken in allocating:the new transition probabilities so that the variance accumulated over these cycles does not increase without bounds. Based on this observation we develop two importance-sampling techniques that have the bounded; relative error property, i.e., the simulation run-length required to estimate the, rare-event probability to a fixed degree of accuracy remains bounded the event of interest becomes more rare.
引用
收藏
页码:547 / 562
页数:16
相关论文
共 25 条
[1]  
BLUM A, 1994, 24 ANN INT S FAULT T, P137
[2]  
BLUM A, 1993, 219S IBM
[3]  
CARRASCO JA, 1992, P 5 INT C MOD TECHN, P351
[4]   EFFECTIVE BANDWIDTH AND FAST SIMULATION OF ATM INTREE NETWORKS [J].
CHANG, CS ;
HEIDELBERGER, P ;
JUNEJA, S ;
SHAHABUDDIN, P .
PERFORMANCE EVALUATION, 1994, 20 (1-3) :45-65
[5]  
CONWAY AE, 1987, P 17 INT S FAULT TOL, P230
[6]   SIMULATING STABLE STOCHASTIC SYSTEMS .3. REGENERATIVE PROCESSES AND DISCRETE-EVENT SIMULATIONS [J].
CRANE, MA ;
IGLEHART, DL .
OPERATIONS RESEARCH, 1975, 23 (01) :33-45
[7]   IMPORTANCE SAMPLING FOR STOCHASTIC SIMULATIONS [J].
GLYNN, PW ;
IGLEHART, DL .
MANAGEMENT SCIENCE, 1989, 35 (11) :1367-1392
[8]   MODELING AND ANALYSIS OF COMPUTER-SYSTEM AVAILABILITY [J].
GOYAL, A ;
LAVENBERG, SS .
IBM JOURNAL OF RESEARCH AND DEVELOPMENT, 1987, 31 (06) :651-664
[9]   A UNIFIED FRAMEWORK FOR SIMULATING MARKOVIAN MODELS OF HIGHLY DEPENDABLE SYSTEMS [J].
GOYAL, A ;
SHAHABUDDIN, P ;
HEIDELBERGER, P ;
NICOLA, VF ;
GLYNN, PW .
IEEE TRANSACTIONS ON COMPUTERS, 1992, 41 (01) :36-51
[10]  
GOYAL A, 1986, P 16 INT S FAULT TOL, P84