We use antagonistic stochastic games and fluctuation analysis to
examine a single-server queue with bulk input and secondary work during server's multiple vacations. When the buffer contents become exhausted the server leaves the system to perform some diagnostic service of a minimum of L jobs clustered in packets of random sizes (event A). The server is not supposed to stay longer than T units of time (event B). The server returns to the system when A or B occurs, whichever comes first. On the other hand, he may not break service of a packet in a middle even if A or B occurs. Furthermore, the server waits for batches of customers to arrive if upon his return the queue is still empty. We obtain a compact and explicit form functional for the queueing process in equilibrium.
1. Introduction
In this paper we use a game-theoretic analysis for a single-server queue with maintenance. When the queue becomes empty, the server leaves the system. It is usually associated with vacations. However, in our case, the server is vacant from the system not just to perform some unspecified maintenance (as it is the case in most systems with server vacations), but to render a random quantity of jobs. This maintenance can be associated with a semiroutine diagnostic, as computer servers render from time to time. When the server has a break, it is a right time to perform a background work.
In our case, the server will be facing random amounts of various diagnostic procedures and he will have to do some minimum quantity of them. In addition, his vacation from the main system should also be limited by some reasonable time limit. Consequently, the server is scheduled to be vacant until one of the two events will take place first: the total amount of jobs performed is L or the time he spends in the maintenance facility is T. The server's absence policy is however not rigid, due to one more factor. Jobs are included in packets of random quantities (as it is the case in computer networks and telecommunications), so that by no means can the server break a packet with unfinished jobs and return to the system even if the total amount of jobs crosses L or his time spent with maintenance expires beyond T. The server alone cannot break a single job in a middle if his time in the maintenance facility is expired. The server will ultimately return to the primary system after he is done with a packet. This policy is more reasonable than a rigid one, but the system modeling is obviously more complex.
Consequently, the return time, the amount of done jobs, and most importantly, the number of customers so far accumulated in the buffer become a nontrivial problem that is investigated by means of fluctuation analysis. More specifically, it seems plausible to introduce a multivariate random walk process, with two antagonistic components (the number of jobs and the cumulative time), say players A and B. In light of stochastic game theory, the game will be over when a player wins. However, the information on player's victory is not available in the real time, but at one of the epochs of time designated for (random) observations. It means that the end of the game does not occur at the moment when the defeat takes place, but with a random delay. Thus, the information collected about the game can be crude or fine dependent on the intensity of the observation process. Now, due to delays, and because the game will actually continue even after one of the players is defeated, a delayed information can cause the winner to endure damages exceeding his threshold. For example, if server's time in the maintenance facility is expired, but he is not done with a packet (an analog of a delayed observation), he remains in the facility until he completes the packet. Meanwhile, the total quantity of jobs he processed by then also exceeds L. Then, he returns to the system with both, the amount of jobs done greater than L and the time spent at the facility greater than T. (Here the game is pronounced to be over, at the opportune observation epoch, with both players winning.)
We need to notice that there is one more player C involved in the game, though passively. This player receives strikes from both players A and B, but never responds and his defeat is imminent. It occurs either during the game between players A and B, or thereafter when either one of the players or their coalition attacks C until C is defeated. In our case, player C is represented by the flow of customers arriving to the main system during server's work at the maintenance facility and possibly during server's waiting time in the system.
In a nutshell, the game aspect of this analysis lies in two active players exerting random damages to each other at random times. The first phase of the game continues until the casualties to one of the players exceed a certain threshold (L for player A and T for player B) upon one of the legitimate “observation” times (successive completion of servicing packets). When this takes place, the first phase of the game is over. There is a third player (C) who suffers casualties (from both players) but does not actively participate in neither phase of the game. At the beginning of phase two, the total amount of casualties to player C is observed and the entire game is over if there is at least one. Otherwise, the game enters its second phase and continues until the player C suffers one or more hits. (Each strike carries out multiples hits.) Player C is represented by the input flow of customers to the system and the “total casualties” to player C is the amount of customers accumulated in the system during the two phases of the game.
In the past analyses of antagonistic games [1–13], the focus was on warfare [2–6, 8, 12], economics [1, 5, 7, 9, 11], and ecology [10]. For the first time (to the best of our knowledge), stochastic games are directly applied to queueing in this form.
Description of the System
The input to the system is a marked Poisson process with position independent marking and the FIFO discipline of lined up units. The server takes only one unit at a time and have it departed once being processed. (More on this below.)
When the queue becomes exhausted, the server leaves the system (for special vacations) beginning with his diagnostic checkup. Suppose that the server first observes a batch of X1 jobs that he needs to complete. Each job needs a random amount of time beginning with δ1. So, the X1 jobs the server work on sequentially willbe done in Δ1:=δ1+⋯+δX1 units of time.
After this is finished, the server turns to another batch of jobs, X2 and so on. We assume that Xi≥1, a.s. The total amount of jobs to be rendered is limited but some minimum of them, say L, must be rendered before the server returns to the system. In addition, the total amount of time needed to complete these jobs is also limited by a fixed positive number T. We need to clarify this constraints. In general the total quantity of all jobs is limited by L. However, the server must not break a batch to comply with this rule. In other words, if for some k, he will be serving k batches with
Ak:=X1+⋯+Xk≥L,
while
Ak-1=X1+⋯+Xk-1<L,k=1,2,…,X0=0.
Analogously, the time the server spends on diagnostic is limited by T, but never at expense of breaking a job or even a packet (i.e., a batch), so that
τk:=Δ1+⋯+Δk≥T,
but
τk-1=Δ1+⋯+Δk-1<T,k=1,2,…,Δ0=0.
Therefore the server returns to the system when one of the two events take place (whichever comes first): the cumulative amount of jobs exceeds L or the total time spent with the jobs exceeds T after completion of a last batch of jobs (any of these taking place for the first time).
If τn is as in (1.3), we obviously have
τn=δ1+⋯+δX1+δX1+1+⋯+δX2+⋯+δXn,
as the time spent for X1+⋯+Xn jobs.
We need to find how many jobs were actually done upon server's return to the system, as well as the return time, and most importantly, how many customers were so far accumulated in the buffer. We first define the random indices
ν1=inf{n:An=X1+⋯+Xn≥L},ν2=inf{n:τn≥T},ν:=min{ν1,ν2}.
With ν so defined, we see that τν will be the return time (also referred to as the first passage time or the exit time), while Aν and Bν are the total quantity of jobs done during the maintenance time and the total number of customers fill in the buffer, respectively, upon server's return to the system. Notice that if the server, upon his return, finds no customers in the system, he will rest and wait for a first batch to come. We will discuss this new phase in due course.
Consequently, we will be interested in the joint functional
Φν(u,z,θ):=EuAνzBνe-θτν,
in which both Aν andτν are “active components,” while Bν is a “passive component”. The reason for these names is that Aν and τν both actively involved in server's return and their values are related to L and T, respectively, while Bν is not and its value is only assumed on server's return.
Aν is the amount of jobs rendered by the server upon his return to the system, Bν is the number of customers by the end of the maintenance period, and τν is the first passage time, that is, the time when the server returns to the system and thus exits the maintenance mode.
The Game
In one of the setting for stochastic antagonistic games [4, 14–17], there is a bivariate marked point process
𝒢1:=∑k=0∞(Xk,Yk)ετk
representing incremental casualties Xk and Yk to players A and B observed at time τk,k=0,1,…. The players actions in reality take place at different times, but they are observed jointly upon times 𝒯={τ0,τ1,…}. The first passage time takes place at one of the observed epochs from 𝒯 when the casualties to one of the players become “unbearable.”
In our case, however, the two players will be represented by incremental marks Xk and τk-τk-1 instead, while marks Yk will be continually observed and associated with the casualties to the third player C whose status will be raised upon moving on to phase two of the game. (One of the authors of this paper has been involved in multiphase games [15, 18–20] directly related to this model.)
Without loss of generality, we will be using the same notation for the random measure 𝒢1, however meaning Xk's and τk-τk-1's as active components representing players A and B, while Yk is meant to be a passive component representing player C.
In the second phase the game turns to
𝒢2:=∑j=0∞Zjεtj,
with Z0=Yν and t0,t1,… are the times when the winner of phase one exerts hits to player C in increments of Z0,Z1,… until the total casualties Cm:=Z0+⋯+Zm to C at tm cross level N. The total game is over.
The Queue on Phase II
When the server returns to the system at τν, he finds Bν customers in the buffer, and unless Bν≥N, the server does not resume his service and waits for more units to arrive at times t1,t2,…. Once, Cm≥N, a new busy period begins. In this paper, we assume for simplicity, N=1 and leave a more general case of N for an upcoming work (in its association with the N-policy). Nevertheless, the use of fluctuation analysis on the second phase with N=1 is almost mandatory, because we are not able to discern analytically whether or not Bν>0.
The Layout of the Paper
The introductory section is followed by the formalism of our model and the use of fluctuation analysis. The named functional Φν(u,z,θ):=EuAνzBνe-θτν is obtained in a closed form with the focus on the marginal functional EzBνe-θτν (in Section 3) giving the number of units accumulated in the buffer upon exit from the secondary work (phase I). The tractability of EzBνe-θτν is then demonstrated on a special case of the marked point process 𝒢1 in phase I (with no restriction on the queueing system), allowing us in Section 4 to arrive at fully explicit formulas. Section 5 deals with phase II of the game (waiting time after server's return) and the output functional Φ̃μ(z,θ):=EzCμe-θtμ, where Cμ is the number of customers in the system at time tμ (the exit from phase II and the beginning of a busy period). Section 6 offers a Kendall-like compact formula for the queue in the steady state in the form of the probability generating function. The paper concludes with a summary of the results obtained and the discussion of the upcoming work on continuous time parameter queueing process and an extended global control under the revival of rendered jobs during the maintenance period.
A Relationship to the Existing Literature
The present model generalizes the very old classical MX/G/1/∞ queue [21, page 301], but it goes beyond it. It also modifies N- and T-policy (cf. Tian and Zhang [22, pages 27-28 and 32], resp.), multiple vacations systems (cf. Jain et al. [23, page 357]) all in one, let alone the inclusion of the secondary work, and not just vacations.
We use the method of fluctuation analysis [4–6, 14–19, 24, 25], some of which previously applied to antagonistic games [4, 5, 15, 18, 19] and described the system (at least in part) as a game. To the best of our knowledge such an approach is novel and only fragments of this model are related to the past work. A lion portion of this paper is related to sequential games [15, 18–20]. The latter includes two or more phases of a game, mostly a conflict between two players. For example article [15] introduced a game of two players A and B, in which phase I involved a series of one-sided attacks on player A (such as economic sanctions, small ambushes, or threats). They continued until damage to player A became serious enough to respond. The latter turned the game to phase II with mutual attacks, in which attacks to player A became more vigorous than during phase I. In another paper modeling a sequential game [19], phase I runs a mutual conflict with weak attacks on each other that were “tolerated” until cumulative damages to one of the players were crossing some thresholds. Once this took place, the second phase began with more aggressive and dangerous attacks leading to the destruction of one of the players. The paper dealt with the exit time from the game as well as the cumulative damages to either player at the exit. In various variants of such or a similar game, the authors treated discrete- or continues-valued components (cf. [4]). In some cases, like in the present variant of the conflict, the components are mixed. The variety of the nature of the component suggests different types of inverse operators involved to solve the problem. Article [25] by the second author offers an introduction to such different cases. Other work on sequential games can be found in [26–30].
The notion of T-policy as stated in Tian and Zhang [22] and earlier in Heyman [31] is to call the server back as soon as his vacation time hits T. A more plausible scenario of T-policy is to call the server back after his vacation time hits T, but not at expense of breaking his vacation in the middle, especially if it is a multiple vacation case. Furthermore, if the server (as it is in our case) performs explicit and not virtual jobs, then it makes sense to restrict the exit from his vacations by not only having him to end a job, but also a packet of jobs. So, formally we do not generalize [31], but modify it to a seemingly more realistic scenario. As far as N-policy alone (i.e., if T=∞), we also modify it due to the server's obligation not to break a random batch of jobs until it is finished. A classic and most often cited to paper on N-policy and multiple vacations is by Lee et al. [32], which is somewhat related to our present model.
While a majority of work on N-policy with (or without) multiple vacations use the “supplementary variable” technique by finding the probability generating function (pgf) of the continuous time parameter queueing process (cf. Lee et al. [32]), we work on the embedded process upon departures. However, the results we obtain here will be used later on to have the pgf of the continuous time parameter queue using techniques of time-sensitive functionals, which offers benefits of no mandatory assumption on the service time distribution to be absolutely continuous (used in the supplementary variables method).
As far as the techniques utilized in this paper, we relate to fluctuations of stochastic processes [24, 33–35] widely used in queueing [36, 37], economics [14, 17, 38–40], physics [41–44], astronomy [42], and biology [42]. Pertinent cases of fluctuation phenomena in single- and multivariate marked point processes are treated by the second author and his collaborators in [14, 16, 17, 24, 25]. The reader is also advised to see more general work in [33–35].
Finally our present paper is also related to Abolnikov et al. [36] about “hybrid queues” which comes closest of all to our present work. The authors there considered a queue also with maintenance jobs similar to ours, but without a time restriction, and most importantly, without explicit treatment of the process during the maintenance period.
2. Formalism of the Model and the Main Past Results
To formulate a result previously obtained in Dshalalow and Liew [16] we utilize a generic bivariate marked point process. Let
(𝒜,ℬ,𝒯):=∑i≥0(Xi,Yi)ετi,(Xi,Yi)∈{0,1,…},τi≥0
(εais the Dirac mass) be a delayed marked renewal process on a probability space (Ω,ℱ(Ω),P). We assume that (𝒜,ℬ,𝒯) is with position dependent marking such that
γ0(u,z,θ):=EuX0zY0e-θτ0,∥u∥≤1,∥z∥≤1,Re(θ)≥0,γ(u,z,θ):=EuX1zY1e-θ(τ1-τ0),∥u∥≤1,∥z∥≤1,Re(θ)≥0.
Given some fixed thresholds L=1,2,…, and T>0, we define
ν1:=inf{n≥0:An=X0+⋯+Xn≥L},ν2:=inf{n≥0:τn=Δ0+⋯+Δn>T},ν:=min{ν1,ν2},
the random index such that the cumulative marks Aν and τν will exceed L or T, respectively. Notice that Bν is the only one of the three random components that assumes a value at τν, making ℬ=(Y0,Y1,…) a passive component. Since 𝒜=(X0,X1,…) and 𝒯=(Δ0,Δ1,…) are directly related to the control levels L and T, they are referred to as active components of (𝒜,ℬ,𝒯). The following functional:
Φν:=Φν(u,z,θ)=EuAνzBνe-θτν
will describe the status of the process uponits exit from the rectangle [0,L)×[0,T] at τν.
We will make use of the Laplace-Carson transform
ℒ𝒞p(·)(s):=s∫p=0∞e-sp(·)dp,Re(s)>0,
whose inverse obviously is
ℒ𝒞s-1(·)(p)=𝔏s-1(·1s)(p),
where 𝔏-1 denotes the inverse of the Laplace transform.
Next, we introduce another (discrete) transformation
Dp{f(p)}(x):=∑p=0∞xpf(p)(1-x),∥x∥<1,
where f is an integrable function defined on set ℕ0. The inverse operator below can restore f, if we apply it for all k:
𝒟xk(Dp{f(p)}(x))=f(k),k=0,1,…,
where the inverse 𝒟k is
k↦𝒟xkφ(x,y)={limx→01k!∂k∂xk[11-xφ(x,y)],k≥0,0,k<0,
if applied to a function φ(x,y) analytic at zero in the first variable.
Theorem 2.1 by Dshalalow and Liew [16] given below establishes an explicit formula for Φν.
Theorem 2.1.
The functional Φν of (2.4) satisfies the following formula:
Φν=Φν(u,z,θ)=ℒ𝒞y-1[γ0(xu,z,θ+y)1-γ(xu,z,θ+y)𝒟xL-1{γ0(xu,z,θ+y)1-γ(xu,z,θ+y)γ0(u,z,θ)-γ0(xu,z,θ+y)+γ0(xu,z,θ+y)1-γ(xu,z,θ+y)[γ(u,z,θ)-γ(xu,z,θ+y)]γ0(xu,z,θ+y)1-γ(xu,z,θ+y)}](T).
In the sequel, we will drop the delay in (𝒜,ℬ,𝒯) by setting
X0=Y0=τ0=0,
thereby making
γ0(u,z,θ)=1.
The latter will reduce (2.10) to
Φν=ℒ𝒞y-1[𝒟xL-1{γ(u,z,θ)-γ(xu,z,θ+y)+1-11-γ(xu,z,θ+y)}](T)=1-[1-γ(u,z,θ)]ℒ𝒞y-1[𝒟xL-1{11-γ(xu,z,θ+y)}](T).
In particular, for u=1 we obtain the marginal functional
EzBνe-θτν=1-[1-γ(1,z,θ)]ℒ𝒞y-1[𝒟xL-1{11-γ(x,z,θ+y)}](T).
We will now focus on γ(u,z,θ). With the conventional formula for the double expectation we have
γ(u,z,θ)=E[uX1zY1e-θΔ1]=E[uX1E[zY1e-θ(δ1+⋯+δX1)∣X1]].
Now, we can regard Y1 as the sum of customers that enter the system during the times δ1,…,δX1, that is,
Y1=y1+⋯+yX1
with
yj={numberofcustomersaccumulatedin(δj-1,δj]},j=1,…,X1.
Since (y1,δ1),(y2,δ2)… are conditionally (given σ-algebra σ(X1)) independent and jointly identically distributed random vectors, we have
E[zy1e-θδ1⋯zyX1e-θδX1∣σ(X1)]=[g(z,θ)]X1,
where
g(z,θ):=Ezy1e-θδ1
is the common Laplace-Stieltjes transform (LST) of the r.v.'s δ1,δ2,…. Using (2.18)-(2.19) we have
γ(u,z,θ)=E[uX1[g(z,θ)]X1]=G(ug(z,θ)),
where
G(v):=EvX1,X1≥1,
is the probability generating function (pgf) of r.v. X1.
Finally, substituting (2.20) in (2.14), we get
EuAνzBνe-θτν=1-[1-G(ug(z,θ))]ℒ𝒞y-1[𝒟xL-1{11-G(xug(z,θ+y))}](T),
giving the joint marginal transform of r.v.'s Bν (the number of units upon server's return to the system) and τν (the time of server's return). We note that the closest functional to (2.22) was the one considered in Abolnikov et al. [36], but without ℒ𝒞y-1, as no restriction in T was there imposed. Consequently, the vacation model there was simpler from the fluctuations stand. In the rest of the comprehensive literature on vacations (cf. the monograph [22]), there is none to deal with secondary jobs and fluctuation analysis that will bear fruits as we will see it from the upcoming sections.
3. EuAνzBνe−θτν Adopted to the Queueing System
In this section we return to the queueing system and adopt the generic bivariate marked point process (𝒜,ℬ,𝒯) of the previous section, now interpreting the component 𝒜 as the maintenance routine taking place during server's vacationing from the main system as observed over 𝒯. The component ℬ describes the accumulation process of the incoming units in the system during the server's absence. Considering that the input to the system is marked Poisson (also referred to as compound Poisson), with position independent marking and intensity λ,
ℐ=∑i=1∞Uiεti,
where the marks Ui's are iid r.v.'s with the common pgf and the mean
a(z):=EzUi,a:=EUi,i=1,2,….
Assuming that the times δ1,δ2,… needed to process successive jobs during the maintenance period are iid r.v.'s with the common LST
δ(θ):=Ee-θδ1.
we easily conclude that the joint functional g(z,θ) of the number of customers entering the main facility during processing any job (left behind in (2.19)) will then be
g(z,θ)=δ(θ+λ-λa(z)),
equation (3.4) is due to the usual double expectation routine occurring in queueing.
Consequently, with (3.4), formula (2.22) for the joint transform of the customers accumulated during the entire maintenance process and the server's return time will be
EuAνzBνe-θτν=1-[1-G(uδ(λ-λa(z)+θ))]×ℒ𝒞y-1[𝒟xL-1{11-G(xuδ(y+λ-λa(z)+θ))}](T),L≥1.
Specifically, for L=1, EuAνzBνe-θτν degenerates to
EuAνzBνe-θτν=G(uδ(λ-λa(z)+θ)),L=1.
4. A Special Case of the Maintenance Process
First, we will avoid a relative simplicity with L=1, as in this case the operator 𝒟L-1 will degenerate. So, we assume that L≥2. Next, make two assumptions on the maintenance process and go with the special case.
There is the LST of the processing time of a job
Ee-θδ1=δ(θ)=δδ+θ.
The distribution of the batch sizes X1,X2,… is geometric with parameter p and pgf
G(u)=pu1-qu,p+q=1.
Given the above assumptions, we calculate γ(x,z,y) of (2.20) as needed for (3.5). After straightforward algebra we have
γ(xu,z,y+θ)=G(xuδ(y+λ-λa(z)+θ))=pxu(δ/(δ+λ-λa(z)+y+θ))1-qxu((δ/δ+λ-λa(z)+y+θ))=pxuδδ+λ-λa(z)+y+θ-δxuq,1-γ(xu,z,y+θ)=δ+λ-λa(z)+y+θ-δxuq-δxupδ+λ-λa(z)+y+θ-δxuq=δ+λ-λa(z)+y+θ-xuδδ+λ-λa(z)+y+θ-δxuq,11-γ(xu,z,y+θ)=11-G(xuδ(y+λ-λa(z)+θ))=δ+λ-λa(z)+y+θ-δxuqδ+λ-λa(z)+y+θ-xuδ=q+pδ+λ-λa(z)+y+θδ+λ-λa(z)+y+θ-xuδ.
We first apply the operator 𝒟xL-1 to (4.5) considering that it is linear and has constant functions as its fixed points (cf. [4, 14, 25]):
𝒟xL-1{11-γ(xu,z,y+θ)}=𝒟xL-1{11-G(xuδ(y+λ-λa(z)+θ))}=𝒟xL-1{q+pδ+λ-λa(z)+y+θδ+λ-λa(z)+y+θ-xuδ}=q+𝒟xL-1{pδ+λ-λa(z)+y+θδ+λ-λa(z)+y+θ-xuδ}=q+p𝒟xL-1{11-xuδ/(δ+λ-λa(z)+y+θ)}.
Furthermore, using the property that
𝒟xk{11-bx}=1-bk+11-b
(cf. [4]) we obtain
𝒟xL-1{11-γ(xu,z,y+θ)}=𝒟xL-1{11-G(xuδ(y+λ-λa(z)+θ))}=q+p1-[δu/(δ+λ-λa(z)+y+θ)]L1-δu/(δ+λ-λa(z)+y+θ).
We now prepare 𝒟xL-1{1/(1-γ(xu,z,y+θ))}=𝒟xL-1{1/(1-G(xuδ(y+λ-λa(z)+θ)))} for the further application of ℒ𝒞y-1, and after a simple algebra and partial fractions decomposition we arrive at
𝒟xL-1{11-γ(xu,z,y+θ)}=𝒟xL-1{11-G(xuδ(y+λ-λa(z)+θ))}=1+pδuδ-δu+λ-λa(z)+y+θ-p(δu)L(δ-δu+λ-λa(z)+y+θ)[δ+λ-λa(z)+y+θ]L-1.
To apply the Laplace-Carson inverse we need to divide 𝒟xL-1{1/(1-γ(xu,z,y+θ))} by y first and then use the Laplace inverse. Using again the partial fraction decomposition, after the division by y, we have the expression
1y𝒟xL-1{11-γ(xu,z,y+θ)}=1y+pδuδ-δu+λ-λa(z)+θ(1y-1δ-δu+λ-λa(z)+θ+y)-p(δu)Lδ-δu+λ-λa(z)+θ×[1y(δ+λ-λa(z)+θ+y)L-1-1(δ-δu+λ-λa(z)+θ+y)(δ+λ-λa(z)+θ+y)L-1].
Applying the Laplace inverse to the first two terms of (4.10) and using the standard tables for the Laplace transform we have
ℒy-1(1y+pδuδ-δu+λ-λa(z)+θ(1y-1δ-δu+λ-λa(z)+θ+y))(T)=1+pδuδ-δu+λ-λa(z)+θ(1-e-(δ-δu+λ-λa(z)+θ)T).
To process the third term we first notice that
ℒy-1{1[δ+λ-λa(z)+θ+y]L-1}(T)=e-(δ+λ-λa(z)+θ)TTL-2(L-2)!
and, in particular, for δ=0 and L=2,
ℒy-1{1λ-λa(z)+θ+y}(T)=e-(λ-λa(z)+θ)T.
Then we will use the property of the Laplace transform applied to the convolution of two functions h and g:
ℒp(h*g(p))(y)=ĥ(y)ĝ(y).
Consequently, the inverse of the product of two Laplace images is the convolution of their inverses. So,
ℒy-1(1(δ-δu+λ-λa(z)+θ+y)(δ+λ-λa(z)+θ+y)L-1)(T)=∫u=0Te-(δ+λ-λa(z)+θ)vvL-2(L-2)!e-(δ-δu+λ-λa(z)+θ)(T-v)dv=e-(δ-δu+λ-λa(z)+θ)T∫v=0Te-δvvL-2(L-2)!dv=e-(δ-δu+λ-λa(z)+θ)T(δu)L-1[1-∑j=0L-2e-δuT(δuT)jj!].
From (4.10)–(4.13) and (4.15)-(4.16) we obtain
ℒ𝒞y-1[𝒟xL-1{11-G(xuδ(y+λ-λa(z)+θ))}](T)=1-pδuH(δu,z,θ)Δ(u,z,θ),
where
H(δu,z,θ)=(δuh(δ,z,θ))L-1Σ(δ,z,θ)+e-Δ(u,z,θ)T[1-Σ(δu,1,0)]-1,Σ(δ,z,θ)=1-∑j=0L-2e-h(δ,z,θ)T[h(δ,z,θ)T]jj!,h(δ,z,θ)=δ+λ-λa(z)+θ,Δ(u,z,θ)=δ-δu+λ-λa(z)+θ.
From (3.5), using (4.4) for x=1 and (4.17)–(4.21), after simple algebra, we finally have
EuAνzBνe-θτν=1-[1-G(uδ(λ-λa(z)+θ))]×ℒ𝒞y-1[𝒟xL-1{11-G(xuδ(y+λ-λa(z)+θ))}](T)=1-Δ(u,z,θ)h(δ,z,θ)-δuq(1-pδuH(δu,z,θ)Δ(u,z,θ))=pδu1+H(δu,z,θ)h(δ,z,θ)-δuq,L≥2.
In particular,
EzBνe-θτν=Φν(1,z,θ)=pδ1+H(δ,z,θ)h(δp,z,θ),L≥2,
where
H(δ,z,θ)=(δh(δ,z,θ))L-1Σ(δ,z,θ)+e-h(0,z,θ)T[1-Σ(δ,1,0)]-1,Σ(δ,1,0)=1-∑j=0L-2e-δT[δT]jj!.
From (3.6) and (4.3), for L=1 we get
EuAνzBνe-θτν=puδδ+λ-λa(z)+θ-δuq,L=1.
In particular,
EzBνe-θτν=pδpδ+λ-λa(z)+θ,L=1.
In a nutshell we have the following.
Theorem 4.1.
Given assumptions (4.1) and (4.2), the number of customers Bν that enter the system during the server's maintenance work (i.e., between his departure from the system and his return at τν) along with the time τν and the number of jobs Aν rendered by τν in the form EuAνzBνe-θτν satisfies formulas (4.18)–(4.22) and (4.26). The corresponding marginal functional EzBνe-θτν satisfies formulas (4.19)-(4.20), (4.23)–(4.25), (4.27).
5. Phase 2: The Waiting Period
Now, when the server returns to the system, he finds Bν customers in the buffer, which does not have to be positive. In case if Bν=0, the server does not leave the system as he did in the conclusion of a busy period, and so he waits for any batch of customers. Consequently, his service finally resumes when a group of units enters the system. However, it sounds like most likely customers will wait for the server rather than other way around. The problem is not trivial as it is for a regular MX/G/1/∞ queue, with the server waiting for the first batch of customers unconditionally. We turn again to fluctuation analysis considering a generic delayed marked renewal process
(𝒵,T):=∑k≥0Zkεtk,
with generally position dependent marking. The process 𝒵 runs its course until its cumulative mark
Cn=Z0+⋯+Zn
crosses some integer N at tn. In other words, with
μ=inf{n:Cn≥N}Cμ≥N at the exit time tμ. We describe the process (𝒵,T) by
Γ0(z,θ):=E[zZ0e-θt0],∥z∥≤1,Re(θ)≥0,Γ(z,θ):=E[zZ1e-θ(t1-t0)],∥z∥≤1,Re(θ)≥0.
We are interested in the joint functional
Φ̃μ(1,z,θ):=E[zCμe-θtμ].
The latter is known [16] to obey the formula
Φ̃μ(1,z,θ)=Γ0(z,θ)-[1-Γ(z,θ)]𝒟xN-1{Γ0(xz,θ)1-Γ(xz,θ)},
with
Γ0(z,θ)=Φν(1,z,θ),
when bearing in mind a period when the server returns to the system, with the history of
t0=τν,Z0=Bν.
Note that the server would then wait until N customers or more will accumulate upon one of the arrivals. (Recall that the input is bulk.) In the event, Z0≥N, Γ0(z,θ) will contain factor zN (with the rest being analytic in z) and thus make the second term of (5.6) vanishes, as
𝒟xN-1(xNφ(1,xz,θ))=𝒟x0(xφ(1,xz,θ))=0.
In this case, Φ̃μ(z,θ) would equal Γ0(z,θ). Thus, formula (5.6) will take care of the case when Z0≥N by itself.
Obviously,
Γ(z,θ)=a(z)λλ+θ
in our case of the marked Poisson process with position independent marking. Consequently, we have the marginal functional
Φ̃μ(1,z,θ)=E[zCμe-θtμ]=Φν(1,z,θ)-[1-λλ+θa(z)]𝒟xN-1{Φν(1,xz,θ)1-λ/(λ+θ)a(xz)}.
Furthermore, we assume that the threshold N=1 (although in some upcoming versions we would like to retain the generality of N) thereby further reducing (5.11) to
Φ̃μ(1,z,θ)=Φν(1,z,θ)-[1-λλ+θa(z)]𝒟x0{Φν(1,xz,θ)}=Φν(1,z,θ)-[1-λλ+θa(z)]Φν(1,0,θ),
where
Φν(1,z,θ)=δp1+H(δ,z,θ)h(δp,z,θ).
To find Φν(1,0,θ) explicitly we turn to formulas (4.18)–(4.21):
Φν(1,0,θ)=δp1+H(δ,0,θ)δp+λ+θ,H(δ,0,θ)=(δδ+λ+θ)L-1[1-∑j=0L-2e-(δ+λ+θ)T[(δ+λ+θ)T]jj!]-[1-∑j=0L-2e-(δ+λ+θ)T(δT)jj!].
The above can be summarized as Theorem 5.1. For convenience, we bring all formulas together in one place.
Theorem 5.1.
The number of units Cμaccumulated in the buffer upon the beginning of a busy period tμ, jointly with tμ, for L≥2, is given by the transform
E[zCμe-θtμ]=Φ̃μ(1,z,θ)=Φν(1,z,θ)-[1-λλ+θa(z)]Φν(1,0,θ),
where
Φν(1,z,θ)=δp1+H(δ,z,θ)h(δp,z,θ),H(δ,z,θ)=(δh(δ,z,θ))L-1Σ(δ,z,θ)+e-h(0,z,θ)T[1-Σ(δ,1,0)]-1,Σ(δ,z,θ)=1-∑j=0L-2e-h(δ,z,θ)T[h(δ,z,θ)T]jj!,h(δ,z,θ)=δ+λ-λa(z)+θ,Φν(1,0,θ)=δp1+H(δ,0,θ)δp+λ+θ,H(δ,0,θ)=(δδ+λ+θ)L-1[1-∑j=0L-2e-(δ+λ+θ)T[(δ+λ+θ)T]jj!]-[1-∑j=0L-2e-(δ+λ+θ)T(δT)jj!],Σ(δ,1,0)=1-∑j=0L-2e-δT(δT)jj!.
The marginal functional EzCμ can be obtained from (5.15)-(5.16) by letting θ=0:EzCμ=Φ̃μ(1,z,0)=Φν(1,z,0)-[1-a(z)]Φν(1,0,0),
where, in particular,
Φν(1,0,0)=δp1+H(δ,0,0)δp+λ=pδpδ+λ{(δδ+λ)L-1Σ(δ,0,0)+e-λT(1-Σ(δ,1,0))},H(δ,0,0)=(δδ+λ)L-1[1-∑j=0L-2e-(δ+λ)T[(δ+λ)T]jj!]-[1-∑j=0L-2e-(δ+λ)T(δT)jj!].
For L=1,Φ̃μ(1,z,θ)=pδpδ+λ-λa(z)+θ-[1-λλ+θa(z)]pδpδ+λ+θ,L=1.
For the sequel, we will need ECμ and for thus, we need EBν=Φν′(1). After straightforward computations, we arrive at the following.
Proposition 5.2.
The mean number of units α:=ECμ at the beginning of a busy period satisfies the following formula:
α=ECμ=EBν+aΦν(1,0,0),L≥2,
where
EBν=λapδ+λa[T+(L-1δ-T)Σ(δ,1,0)-Te-δT(δT)L-2(L-2)!],L≥2,α=ECμ=λapδ+apδpδ+λ,L=1,
with Φν(1,0,0) satisfying (5.18) and a being defined in (3.2).
6. The System on a Busy Period and in Equilibrium
In Section 5 (Theorem 5.1) we calculated the pgf of the quantity of the units accumulated in the buffer at the beginning of a busy period. The time during which the system gained in Cμ customers consisted of two phases. The first phase was when the server performed a maintenance completing Aμ amount of jobs. Phase 2 involved a waiting period in the event the buffer was empty.
The rest of the analysis (as regards the embedded process) is almost identical to that of the usual bulk input MX/G/1/∞ queue. We therefore restrict ourselves to the main formalism and present the final formulas.
The input to the system is bulk Poisson ℐ=∑i=1∞Uiεti with position independent marking and a(z)=EzUi, a=EUi,i≥1, as per (3.1)-(3.2), and intensity λ of its point process ∑i=1∞εti.
Service is general and independent including independence of the input that is, service times σ1,σ2,… are iid r.v. with a common LST
β(θ)=Ee-σ1,Re(θ)≥0,b:=Eσ1.
Q(t) is the right continuou squeueing process (the number of customers at time t≥0).
T1,T2,… are successive departures of individual units.
{Qn:=Q(Tn);n=0,1,…} is the embedded Markov chain upon departures.
The transition probability matrix P=(pij) of {Qn} is a delta-2 matrix (cf. [36]) similar to that of MX/G/1/∞ queue, with only zero row different, with no impact on ergodicity condition ρ:=λab<1 (necessary and sufficient). The zero row contains the transition probabilities over the two-phase maintenance with waiting and service period (zero-service cycle), so that the queue length upon the end of such cycle is Cμ+ℐ([tμ,tμ+σ1)) and its pgf, after a straightforward calculation, is EzCμ+ℐ([tμ,tμ+σ1))=β(λ-λa(z))EzCμ. So, the pgf of the zero-row in the MX/G/1/∞ system is
P0(z):=E[zQ1∣Q0=0]=β(λ-λa(z))α(z)z-1,
where
α(z):=EzCμ
(of Theorem 5.1) will replace a(z) in the MX/G/1/∞ system [21].
Consequently, if Kendall's formula (most often referred to as Pollaczek-Khinchin formula) in the MX/G/1/∞ system is
P(z)=p0β(λ-λa(z))a(z)-1z-β(λ-λa(z)),
(cf. Medhi [21]) all we need to do is to replace a(z) with α(z) to get
P(z)=p0β(λ-λa(z))α(z)-1z-β(λ-λa(z)).
Also, in Kendall's formula, we have
p0=1-ρa.
So, we also replace a with
α=ECμ
of Proposition 5.2. Here P(z) is the pgf of the embedded queueing process in equilibrium.
We can summarize the above as the follows.
Theorem 6.1.
The embedded queueing process {Qn} in the MX/G/1/∞ type queue with a two-phase maintenance is ergodic given ρ<1. Under this condition the pgf of the invariant probability measure p=(p0,p1,…) satisfies the Kendall-like formula
P(z)=1-ραβ(λ-λa(z))α(z)-1z-β(λ-λa(z)),
with α and α(z) of Proposition 5.2 and Theorem 5.1, respectively.
Summary
In conclusion, we studied an exhaustive type MX/G/1/∞ queue with a single server performing secondary maintenance jobs during his vacationing periods. Unlike most of queues with vacations, in which the vacationing time is “anonymous,” our server takes on real jobs. In addition, we introduce a so-called “modified T-Policy” that nonrigidly restricts the server to a time T to be vacant from the system. However, if the server is in a middle of a batch of jobs, he must finish them first before coming back. Thus, the return will take place at the first opportune time. Furthermore, the server must also complete a minimum of L jobs before returning to the system. Therefore, his discharge from maintenance will take place at some epoch of time (first passage time) when the server first completes a minimum of jobs (without interrupting an individual batch) or being vacant during T units of time (with the same restriction), whichever of these two occurs first. To find the number of units Bν present in the system (if any) upon server's return, we use game-theoretic modeling (of two players A and B whoever wins first) and fluctuation analysis for multivariate random walk processes.
Now, coming back, if no unit is waiting in the buffer, then the server waits himself for new units to come. We use again fluctuation analysis to find the total quantity of units Cν upon server's return to work and complete the paper with Kendall-like formula for the stationary distribution of units present in the system upon departure epochs.
In the upcoming work with this model, we plan to use time-sensitive functionals to interpolate the probabilities on departures to those for any time (time dependent process). This requires a different approach and would expand the paper. We would also like to consider a multivalued process which would also revive the number of jobs done during the maintenance periods in order to more comprehensively optimize the values of L and T. Again all these would run out of the scope of this paper that aims at introducing a game-theoretic approach and fluctuations to find the pgf of the equilibrium queue distribution. We notice that a supplementary variables approach would be very difficult if not impossible for our model. In addition, it also restricts the service time distributions to those absolutely continuous.
Acknowledgments
The authors are grateful to anonymous referees and the associate editor for their valuable suggestions that led to a considerable improvement of the paper. This research is supported by the U.S Army Grant no. W911NF-07-1-0121.
BasarT. S.OlsderG. J.1982160Orlando, Fla, USAAcademic Pressxii+430Mathematics in Science and EngineeringMR649247BreitnerM. H.PeschH. J.GrimmW.Complex differential games of pursuit-evasion type with state constraints. I. Necessary conditions for optimal open-loop strategies1993783419441MR124043010.1007/BF00939876ZBL0796.90078DavidovitzA.ShinarJ.Two-target game model of an air combat with fire-and-forget all-aspect missiles198963213316510.1007/BF00939571MR1026862ZBL0662.90103DshalalowJ. H.Random walk analysis in antagonistic stochastic games2008264738783MR243107810.1080/07362990802128354ZBL1151.91330DshalalowJ. H.KeH.-J.Multilayers in a modulated stochastic game2009353255356510.1016/j.jmaa.2008.12.026MR2508958ZBL1173.91308DshalalowJ. H.TreerattrakoonA.Antagonistic games with an initial phase2009932772862-s2.0-0036054811FishburnP. C.Noncooperative stochastic dominance games197871516110.1007/BF01763120MR497568IsaacsR.1999New York, NY, USADoverKonstantinovR. V.PolovinkinE. S.Mathematical simulation of a dynamic game in the enterprise competition problem2004405720725PerryJ. C.RoitbergB. D.Games among cannibals: competition to cannibalize and parent-offspring conflict lead to increased sibling cannibalism2005186152315332-s2.0-000014772610.1111/j.1420-9101.2005.00941.xShashikhinV. N.Antagonistic game with interval payoff functions2004404556564ShimaT.Capture conditions in a pursuit-evasion game between players with biproper dynamics20051263503528MR216480410.1007/s10957-005-5495-3SobelM. J.Noncooperative stochastic games19714219301935MR030956310.1214/aoms/1177693059ZBL0229.90059DshalalowJ. H.On exit times of multivariate random walk with some applications to finance2005635–7e569e5772-s2.0-000212780110.1016/j.na.2005.01.083ZBL1159.60346DshalalowJ. H.HuangW.A stochastic game with a two-phase conflict2008chapter 18Cambridge, UKCambridge Academic Publishers201209DshalalowJ. H.LiewA.On exit times of a multivariate random walk and its embedding in a quasi Poisson process200624245147410.1080/07362990500522478MR2204722ZBL1105.60065DshalalowJ. H.LiewA.On fluctuations of a multivariate random walk with some applications to stock options trading and hedging2006449-1093194410.1016/j.mcm.2006.02.022MR2253752ZBL1133.91416DshalalowJ. H.HuangW.On noncooperative hybrid stochastic games20082380381110.1016/j.nahs.2007.11.001MR2431714DshalalowJ. H.HuangW.Sequential antagonistic games with initial phase (jointly with Weijun Huang)2010chapter 2New York, NY, USANova Science1536HuangW.DshalalowJ. H.Tandem antagonistic games20097112e259e2702-s2.0-5674916431410.1016/j.na.2008.10.051MedhiJ.1991Boston, Mass, USAAcademic Pressxiv+444MR1118150TianN.ZhangZ. G.2006New York, NY, USASpringerxii+385International Series in Operations Research & Management ScienceMR2248264JainJ. L.MohantyS. G.BöhmW.2007Boca Raton, Fla, USAChapman & Hall/CRCxii+461Statistics: Textbooks and MonographsMR2257773AgarwalR. P.DshalalowJ. H.On multivariate delayed recurrent processes20051533549MR2144193ZBL1078.60036DshalalowJ. H.On the level crossing of multidimensional delayed renewal processes1997104355361MR148801810.1155/S1048953397000415BrandtsJ.SolàcC.Reference points and negative reciprocity in simple sequential games2001362138157MR185056710.1006/game.2000.0818ZBL1037.91505ExmanI.Solving sequential games with Boltzmann-learned tactics496Proceedings of the 1st Workshop on Parallel Problem Solving from Nature (PPSN '90)October 1990Dortmund, GermanySpringer216220Lecture Notes In Computer ScienceRadzikT.SzajowskiK.Sequential games with random priority199094361377MR1108811ZBL0745.62080SiegristK.SteeleJ.Sequential games200138410061017MR187655510.1239/jap/1011994188ZBL1002.60082WenQ.A folk theorem for repeated sequential games200269249351210.1111/1467-937X.00214MR1920843ZBL1030.91008HeymanD.The T-policy for the M/G/1 queue19772377577810.1287/mnsc.23.7.775ZBL0357.60022LeeH. W.LeeS. S.ParkJ. O.ChaeK. C.Analysis of the MX/G/1-policy and multiple vacations1994312476496MR127480310.2307/3215040AlsmeyerG.On generalized renewal measures and certain first passage times199220312291247MR117526110.1214/aop/1176989690ZBL0759.60088BinghamN. H.ShanbhagD. N.RaoC. R.Random walk and fluctuation theory200119Amsterdam, The NetherlandsElsevier Science171213Handbook of StatisticsMR1861724ZBL0982.60038TakácsL.RotaG.-C.On fluctuations of sums of random variables19782New York, NY, USAAcademic Press4593Advances in Mathematics, Supplementary StudiesMR517253ZBL0447.60053AbolnikovL.DshalalowJ. H.TreerattrakoonA.On a dual hybrid queueing systems2008219610910.1016/j.nahs.2007.05.002MR2381040TakácsL.On fluctuation problems in the theory of queues197683548583MR042359310.2307/1426144ZBL0357.60021KyprianouA. E.PistoriusM. R.Perpetual options and Canadization through fluctuation theory20031331077109810.1214/aoap/1060202835MR1994045ZBL1039.60044MellanderE.VredinA.WarneA.Stochastic trends and economic fluctuations in a small open economy19927436939410.1002/jae.3950070405MuzyJ. F.DelourJ.BacryE.Modelling fluctuations of financial time series: from cascade process to stochastic volatility model20001735375482-s2.0-003164095110.1007/s100510070131Garrido, Ed.L.Fluctuations and stochastic phenomena in condensed matter268Proceedings of the Sitges Conference on Statistical Mechanics1987Barcelona, SpainSpringerLectures Notes in PhysicsHida, Ed.T.Mathematical approach to fluctuations: astronomy, biology and quantum dynamicsProceedings of the Iias WorkshopMay 1992Kyoto, JapanWorld Scientific10.1287/mnsc.23.7.775HorsthemkeW.KondepudiD.1984New York, NY, USASpringerRednerS.2001Cambridge, UKCambridge University Pressx+312MR1851867