  This paper considers a work-conserving FIFO single-server queue with multiple
batch Markovian arrival streams governed by a continuous-time finite-state
Markov chain. A particular feature of this queue is that service time
distributions of customers may be different for different arrival streams.
After briefly discussing the actual waiting time distributions of customers
from respective arrival streams, we derive a formula for the vector generating
function of the time-average joint queue length distribution in terms of the
virtual waiting time distribution. Further assuming the discrete phase-type
batch size distributions, we develop a numerically feasible procedure to
compute the joint queue length distribution. Some numerical examples are
provided also.
