In this paper we present a technique for generating random variates from an empirical distribution using the matrix exponential representation of the distribution. In our experience, a matrix exponential representation of an empirical distribution produces random variates with an excellent fit with the empirical distribution. This technique is particularly important when the empirical data is very bursty, i.e., has a high variance. In this paper we discuss how to find the matrix exponential representation of an empirical distribution and we present our technique for generating random variates from the empirical distribution using its matrix exponential representation. We show how the matrix exponential representation of an empirical distribution is found through an example and then we show that matrix exponential random variates are an excellent fit with the empirical data through an χ2 goodness-of-fit test.
[1]
C. Loan,et al.
Nineteen Dubious Ways to Compute the Exponential of a Matrix
,
1978
.
[2]
Marcel F. Neuts,et al.
Matrix-Geometric Solutions in Stochastic Models
,
1981
.
[3]
T. Stieltjes.
Recherches sur les fractions continues
,
1995
.
[4]
W. Gragg,et al.
On the partial realization problem
,
1983
.
[5]
Edward F. Brown.
A distribution-free random number generator via a matrix-exponential representation
,
1992,
SAC '92.
[6]
Kenneth D. Lawrence,et al.
Modern Statistical Systems and GPSS Simulation
,
1992
.
[7]
Edward Dudewicz,et al.
The Generalized Bootstrap
,
1992
.
[8]
Lester Lipsky,et al.
Queueing Theory: A Linear Algebraic Approach
,
1992
.