The purpose of this study is to present a subgeometric convergence formula for the stationary distribution of the finite-level M/G/1-type Markov chain when taking its infinite-level limit, where the upper boundary level goes to infinity. This study is carried out using the fundamental deviation matrix, which is a block-decomposition-friendly solution to the Poisson equation of the deviation matrix. The fundamental deviation matrix provides a difference formula for the respective stationary distributions of the finite-level chain and the corresponding infinite-level chain. The difference formula plays a crucial role in the derivation of the main result of this paper, and the main result is used, for example, to derive an asymptotic formula for the loss probability in the MAP/GI/1/N queue.

We consider two classes of irreducible Markovian arrival processes specified by the matrices C and D: the Markov-modulated Poisson process (MMPP) and the Markov-switched Poisson process (MSPP). The former exhibits a diagonal matrix D while the latter exhibits a diagonal matrix C. For these two classes we consider the following four statements: (I) the counting process is overdispersed; (II) the hazard rate of the event-stationary interarrival time is nonincreasing; (III) the squared coefficient of variation of the event-stationary process is greater than or equal to one; (IV) there is a stochastic order showing that the time-stationary interarrival time dominates the event-stationary interarrival time. For general MSPPs and order two MMPPs, we show that (I)–(IV) hold. Then for general MMPPs, it is easy to establish (I), while (II) is shown to be false by a counter-example. For general simple point processes, (III) follows from (IV). For MMPPs, we conjecture that (IV) and thus (III) hold. We also carry out some numerical experiments that fail to disprove this conjecture. Importantly, modelling folklore has often treated MMPPs as “bursty”, and implicitly assumed that (III) holds. However, to the best of our knowledge, proving this is still an open problem.

A multiclass c-server retrial queueing system in which customers arrive according to a class-dependent Markovian arrival process (MAP) is considered. Service and retrial times follow class-dependent phase-type (PH) distributions with the further assumption that PH distributions of retrial times are acyclic. A necessary and sufficient condition for ergodicity is obtained from criteria based on drifts. The infinite state space of the model is truncated with an appropriately chosen Lyapunov function. The truncated model is described as a multidimensional Markov chain, and a Kronecker representation of its generator matrix is numerically analyzed.

The Markovian arrival process generalizes the Poisson process by allowing for dependent and nonexponential interarrival times. We study the autocorrelation function of the two-state Markovian arrival process. Our findings show that the correlation structure of such a process has a very specific pattern, namely, it always converges geometrically to zero. Moreover, the signs of the autocorrelation coefficients are either constant or alternating.

In this paper we consider the problem of identifiability for the two-state Markovian arrival process (MAP2). In particular, we show that the MAP2 is not identifiable, providing the conditions under which two different sets of parameters induce identical stationary laws for the observable process.

An important class of queueing networks is characterized by the following feature: in contrast with ordinary units, a disaster may remove all work from the network. Applications of such networks include computer networks with virus infection, migration processes with mass exodus and serial production lines with catastrophes. In this paper, we deal with a two-stage tandem queue with blocking operating under the presence of a secondary flow of disasters. The arrival flows of units and disasters are general Markovian arrival processes. Using spectral analysis, we determine the stationary distribution at departure epochs. That distribution enables us to derive the distribution of the number of units which leave the network at a disaster epoch. We calculate the stationary distribution at an arbitrary time and, finally, we give numerical results and graphs for certain probabilistic descriptors of the network.

This paper presents an algorithmic procedure to calculate the delay distribution of a type k customer in a first-come-first-served (FCFS) discrete-time queueing system with multiple types of customers, where each type has different service requirements (the MMAP[K]/PH[K]/1 queue). First, we develop a procedure, using matrix analytical methods, to handle arrival processes that do not allow batch arrivals to occur. Next, we show that this technique can be generalized to arrival processes that do allow batch arrivals to occur. We end the paper by presenting some numerical examples.

We study a single-server queue in which the server takes a vacation whenever the system becomes empty. The service and vacation times and the arrival process are all assumed to be mutually independent. The successive service times and the vacation times each form independent, identically distributed sequences with general distributions. A new class of non-renewal arrival processes is introduced. As special cases, it includes the Markov-modulated Poisson process and the superposition of phase-type renewal processes.

Algorithmically tractable equations for the distributions of the waiting times at an arbitrary time and at arrivals, as well as for the queue length at an arbitrary time, at arrivals, and at departures are established. Some factorizations, which are known for the case of renewal input, are generalized to this new framework and new factorizations are obtained. The algorithmic implementation of these results is discussed.