In this article we define a class of distributions called bilateral phase type (BPH), and study its closure and computational properties. The class of BPH distributions is closed under convolution, negative convolution, and mixtures. The one-sided version of BPH, called generalized phase type (GPH), is also defined. The class of GPH distributions is strictly larger than the class of phase-type distributions introduced by Neuts, and is closed under convolution, negative convolution with nonnegativity condition, mixtures, and formation of cohere'nt systems. We give computational schemes to compute the resulting distributions from the above operations and extend them to analyze queueing processes. In particular, we present efficient algorithms to compute the steady-state and transient waiting times in GPH/GPH/ 1 queues and a simple algorithm to compute the steady-state waiting time in M/GPH/ 1 queues.
INTRODUCTIONNeuts [ 171, using the properties of exponential, generalized Erlang, and hyperexponential distributions as a base, developed a phase-type probability distribution function. Through a sequence of several articles, he has established its closure and computational properties, which make its use very attractive in applied probability modeling (see [20] for a recent account of these results). The purpose of this article is to present two classes of distributions, called bilateral phase type and generalized phase type, and to establish their closure and computational properties. The class of generalized phase-type distributions is strictly larger than the class of phase-type distributions.The phase-type distribution of Neuts can be interpreted as a random sum of independent exponential random variables via the uniformization procedure [9; 18,p.78;23;24]. The stopping time for this sum has an associated discrete phase-type distribution. A natural extension of the phase-type distribution is to consider the sum of independent and identically distributed (iid) exponential random variables (rv) over an arbitrary random number. We will call such a sum a generalized phase-type random variable. Let (g(n))," be the probability distribution function of a random number L , and define where (En): are iid exponential rv's with mean l / A and Eo = 0 with probability (wp) 1. The survival function Fx of X is then *Present address: