Abstract

We consider an M/M/1 retrial queue with working vacations, vacation interruption, Bernoulli feedback, and N-policy simultaneously. During the working vacation period, customers can be served at a lower rate. Using the matrix-analytic method, we get the necessary and sufficient condition for the system to be stable. Furthermore, the stationary probability distribution and some performance measures are also derived. Moreover, we prove the conditional stochastic decomposition for the queue length in the orbit. Finally, we present some numerical examples and use the parabolic method to search the optimum value of service rate in working vacation period.

1. Introduction

In the queueing theory, vacation queues and retrial queues have been intensive research topics; we can find general models in Tian and Zhang [1] and Artalejo and Gómez-Corral [2]. In 2002, Servi and Finn [3] first introduced working vacation policy and studied an M/M/1/WV queue. Their work is motivated and illustrated by the analysis of a WDM optical access network using multiple wavelengths which can be reconfigured. The study of queueing system with working vacations can also provide the theory and analysis method to design the optimal lower speed period. Wu and Takagi [4] extended the M/M/1/WV queue to an M/G/1/WV queue. Using the matrix-analytic method, Baba [5] considered a GI/M/1 queue with working vacations. Krishnamoorthy and Sreenivasan [6, 7] analyzed an M/M/2 queue with working vacations.

Furthermore, during the working vacation period, if there are customers at a service completion instant, the server can stop the vacation and come back to the normal working level. This policy is called vacation interruption. In some practical situations, the server can take service in the vacation period and must come back to work at times. For example, when the number of customers exceeds the special value and if the server continues to take the vacation, the costs of waiting customers and providing service in the vacation period will be large. In 2007, Li and Tian [8] first introduced vacation interruption policy and studied an M/M/1 queue. Next, Li et al. [9] analyzed the GI/M/1 queue. Using the method of a supplementary variable, Zhang and Hou [10] considered an M/G/1 queue with working vacations and vacation interruption. Sreenivasan et al. [11] studied an MAP/PH/1 queue with working vacations, vacation interruption, and N-policy.

Retrial queueing systems are described by the feature that the arriving customers who find the server busy join the retrial orbit to try their requests again. Retrial queues are widely and successfully used as mathematical models of several computer systems and telecommunication networks. For example, peripherals in computer systems may make retrials to receive service from a central processor. Choi et al. [12] analyzed an M/M/1 retrial queue with general retrial times. Martin and Gómez-Corral [13] considered an M/G/1 retrial queue with liner control policy. Lillo [14] investigated a G/M/1 retrial queue. Sherman and Kharoufeh [15] studied an M/M/1 retrial queue with unreliable server.

One additional feature which has been widely discussed in retrial queueing systems is the Bernoulli feedback of customers. The phenomena of feedback in retrial queueing systems occurred in many practical situations. For example, the retrial queue with feedback can be used to model the Automative Repeat Request protocol in a high frequency communication network. Falin [16] studied an M/M/1 retrial queue with feedback. Kumar et al. [17] investigated an M/G/1 retrial queue with feedback and starting failures. Ke and Chang [18] considered a modified vacation policy for the M/G/1 retrial queue with balking and feedback. Kumar et al. [19] discussed an M/M/1 retrial queue with feedback and collisions. Kumar et al. [20] analyzed an M/G/1 retrial queue with feedback and negative customers.

Do [21] first studied an M/M/1 retrial queue with working vacations. Zhang and Xu [22] considered an M/M/1/WV queue with N-policy. On the basis of the model in [21, 22], we study an M/M/1 retrial queue with working vacations, vacation interruption, Bernoulli feedback, and N-policy simultaneously. Compared with the model in [23], we will see that the infinitesimal generator is different, and in our model, .

This paper is organized as follows. In Section 2, we introduce the model and obtain the infinitesimal generator. In Section 3, we derive the stability condition and the stationary probability distribution. Some important and interesting performance measures are also given. In Section 4, we prove the conditional stochastic decomposition for the queue length given that the server is busy and there are at least customers in the orbit. Some numerical results are presented in Section 5. Using the parabolic method, a cost minimization problem is also analyzed. Finally, Section 6 concludes this paper.

2. Quasi Birth and Death (QBD) Process Model

In this paper, we consider an M/M/1 retrial queue with working vacation interruption and feedback under N-policy. The detailed description of this model is given as follows.

The interarrival times of customers are exponentially distributed with parameter . Upon the arrival of customers, if the server is busy, customers are forced to wait in the orbit of infinite size. If the server is not occupied, arriving customers get service immediately.

Request retrials from the orbit follow a Poisson process with rate . Upon the arrival of requests, if the server is busy, the retrial customers come back to the orbit. If the server is free, on the other hand, requests get service immediately.

The service time in regular busy period is governed by an exponential distribution with parameter . The service time in working vacation period follows an exponential distribution with parameter .

The server begins a working vacation each time when the system becomes empty, and the vacation time follows an exponential distribution with parameter . When a vacation ends, if there are at least customers in the orbit, the server switches to the normal working level. Otherwise, the server begins another vacation.

In the working vacation period, if there are at least customers in the orbit at a service completion instant, the server will stop the vacation and come back to the normal busy period, which means that vacation interruption happens. If the number of customers in the orbit is less than , the server will continue the vacation.

When a customer completes his/her service, s/he may leave the system with probability or join the retrial group for another service with probability .

We assume that interarrival times, interretrial times, service times, and vacation times are mutually independent. Let be the number of customers in the orbit at time , and let be the state of server at time . There are four possible states of the server as follows:

Clearly, is a Markov process with state space

Using the lexicographical sequence for the states, the infinitesimal generator can be written as where Due to the block structure of matrix , is called a QBD process.

3. Stability Condition and Stationary Distribution

Theorem 1. The QBD process is positive recurrent if and only if .

Proof. The proof of this theorem is similar to the proof of Theorem  3.1 in [23]; we omit it here.

Theorem 2. If , the matrix equation has the minimal nonnegative solution where

Proof. The proof of this theorem is similar to the proof of Theorem  3.3 in [23]; we omit it here.
Under the stability condition, let () be the stationary limit of the process and denote

Note that if there is no customer in the orbit, the probability that the server is free in the normal service period is zero. Thus, .

Theorem 3. If , the stationary probability distribution of is given by for , and where where . Finally, can be obtained by the normalization condition.

Proof. Using the matrix-geometric solution method (see [24]), we have For ,

Subsisting into the above equation, we get (8). On the other hand, satisfies the equation where Thus, we obtain From (20) and (23), we get (10) by some computation. Taking (10) into (20), we get (9). In a similar way, (11) can be obtained from (25) and (28). Taking and into (29), together with (21), (24), and (26), we can derive , , , and . Then, , , and can be obtained from (19), (22), and (27). Let take in (23); using the expressions of , , , and , (13) can be derived. Finally, from the normalization condition , we can obtain .

Remark 4. From Theorem 3, we can see that, in our model, , (), which is different from the result in [23]. However, if we use the same technique to analyze the M/M/1 retrial queue with working vacations and feedback under N-policy but without vacation interruption, cannot be .

From Theorem 3, the probability that the server is busy is given by

The probability that the server is free is

The mean number of customers in the orbit is

The mean number of customers in the system is given by

Let be the waiting time of a customer in the orbit, using Little’s formula, . The expected sojourn time of a customer in the system is .

The system busy period is defined as the period that starts at an epoch when an arriving customer finds an empty system and ends at the departure epoch at which the system is empty. Using the theory of regenerative process, where is the amount of time in the state during a regenerative cycle. Obviously, . Thus, .

4. Conditional Stochastic Decomposition

Lemma 5. If , let be the conditional queue length of an retrial queue with feedback in the orbit given that the server is busy; then has probability generating function

Proof. Consider an M/M/1 retrial queue with feedback; let be the number of customers in the orbit at time , and then is a Markov process with state space . And the infinitesimal generator is given by where Following the steps we used before, the QBD process is positive recurrent if and only if , and the stationary probability distribution is where Thus,

For the model considered in this paper, we introduce a random variable . And is a conditional queue length given that the server is busy and there are at least customers in the orbit. Let be the probability that the server is busy and there are at least customers in the orbit. Obviously,

Theorem 6. If , the conditional queue length can be decomposed into the sum of two independent random variables: , where is defined in Lemma 5 and follows a geometric distribution with parameter . Additional queue length has a distribution

Proof. The proof of this theorem is similar to the proof of Theorem  5.2 in [23]; we omit it here.

5. Numerical Results

5.1. Sensitivity Analysis

In Figure 1, with the change of probability , the curves of probability (the server is busy) and (the server is free) are provided. Figure 2 illustrates the expected queue length with the change of probability at different retrial rate . In Figure 1, we find that increases as increases while decreases as increases. From Figure 2, we can see that decreases evidently with increasing value of . We can easily imagine that will increase dramatically with decreasing, as long as the stability condition in Theorem 1 holds.

From Figures 3 and 4, it is obvious that expected busy period and expected queue length both decrease evidently with service rate increasing. Thus, compared with ordinary vacation policy, working vacation policy can utilize the server and decrease the waiting jobs effectively. And it is easy to see that, if the other conditions are the same, the larger is, the smaller and become.

Under the stability condition, we vary the retrial rate from 3 to 5. Figures 5 and 6 illustrate the effect of on and , respectively. We can see that and both decrease with the rate increasing; this is due to the fact that the interretrial time becomes shorter. When the probability is small, and are sensitive to retrial rate ; this is because customers may join the retrial group for another service with probability .

5.2. Cost Analysis

Queueing managers are always interested in minimizing operating cost of unit time. In this section, we establish a cost function to search for the optimal service rate .

Define the following cost elements: cost per unit time for each customer present in the orbit; cost per unit time for service during a normal service period; cost per unit time for service in a working vacation period; fixed cost per unit time during a working vacation period.

We establish an expected operating cost function per unit time as

Assume that , , , and ; Figure 7 illustrates the curve of cost function with the change of . We can see that there is an optimal service rate to make the cost minimize. In order to solve the optimization problem (44), we can use the parabolic method in [25] to find the optimum value of , say . As is known to us, the unique optimum of a quadratic function agreeing with at 3-point pattern occurs at

Assume the stopping tolerance and with the information of Figure 7, we select the initial 3-point pattern , , and . After six iterations, Table 1 shows that the minimum expected operating cost per unit time converges to the solution with a value of 71.614654.

6. Conclusion

In this paper, we analyze an M/M/1 retrial queue with working vacation, interruption, and feedback under N-policy. Using the matrix-analytic method, the stationary probability distribution and some performance measures are obtained. The conditional stochastic decomposition is also given. We present several numerical examples to study the effect of some parameters. Finally, a cost optimization problem is studied.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the referees and editor for their valuable comments and suggestions. This work is supported by the National Nature Science Foundation of China (no. 11301306), the Post-Doctoral Research Project of China (no. 2013M540041), and the Scientific Research Starting Foundation for Doctors of Shandong University of Technology (no. 413012).