Equilibrium customer and socially optimal balking strategies in a constant retrial queue with multiple vacations and -policy
Abstract
In this paper, equilibrium strategies and optimal balking strategies of customers in a constant retrial queue with multiple vacations and the -policy under two information levels, respectively, are investigated. We assume that there is no waiting area in front of the server and an arriving customer is served immediately if the server is idle; otherwise (the server is either busy or on a vacation) it has to leave the system to join a virtual retrial orbit waiting for retrials according to the FCFS rules. After a service completion, if the system is not empty, the server becomes idle, available for serving the next customer, either a new arrival or a retried customer from the virtual retrial orbit; otherwise (if the system is empty), the server starts a vacation. Upon the completion of a vacation, the server is reactivated only if it finds at least customers in the virtual orbit; otherwise, the server continues another vacation. We study this model at two levels of information, respectively. For each level of information, we obtain both equilibrium and optimal balking strategies of customers, and make corresponding numerical comparisons. Through Particle Swarm Optimization (PSO) algorithm, we explore the impact of parameters on the equilibrium and social optimal thresholds, and obtain the trend in changes, as a function of system parameters, for the optimal social welfare, which provides guiding significance for social planners. Finally, by comparing the social welfare under two information levels, we find that whether the system information should be disclosed to customers depends on how to maintain the growth of social welfare.
Keywords: Multiple vacations, Equilibrium strategies, Balking strategies, Particle Swarm Optimization algorithm, Information accuracy
1 Introduction
In many service and electronic commerce systems, there exists a new trend to study the behavior of customers in queuing models. In these models, customers can decide whether to join or balk, according to a natural tendency to maximize their personal utility. To this end, from the perspective of game-theory, the decentralized behavior of customers in the queuing system has attracted extensive attentions in recent decades. Generally, queuing systems are divided into the observable case and the unobservable case depending on whether customers can obtain the information about the system upon arrival. The observable case was first studied by Naor [16], who analyzed an queue model with a linear reward-cost structure, and obtained equilibrium and social optimal strategies. Subsequently, Naor’s study was extensively extended, see e.g. [6, 11, 18]. Specifically, Edelson and Hilderbrand [6] complemented the unobservable case to Naor’s model. Chen and Frank [2] generalized the model of Naor’s, who assumed that customers and servers use the same discount rate to maximize their expected discount utility. Afterward, some authors studied equilibrium strategies in various invisible models with many different characteristics. The monograph of Hassin and Haviv [10] summarized the main results of the subject under different levels of information.
The present paper aims to discuss equilibrium strategies and socially optimal balking strategies of customers in an constant retrial queue with multiple vacations and the -policy. Customers’ retrials are a common phenomenon in service systems and enterprise engineering. For example, arriving calls to a call center will be connected immediately if service staff is available, otherwise customers may have to retry for service after a random time. With the development of information technology, modern call centers may provide some levels of information to callers, e.g., the current number of customers waiting for service and/or expected waiting time, among other possibilities. Server vacation is another useful concept in modeling for situations, in which optimization of resources and/or reduction of cost are/is required. For vacation models, due to technical and cost (or other) reasons, the server might not be able to obtain the information about the current system capacity during the vacation, or it is impossible for the server to immediately return to work when the number of customers reaches a predetermined threshold, or the number of customers in the system is small at the end of the server vacation so that the server is reluctant to return to work, or return to work at the normal service rate. In addition, too frequent startups and changeovers on operations could lead to severe server wear and overhead on cost, the -policy is usually used to solve this predicament, such as batch traffic transfer systems, controlled manufacturing systems and possible others. In the literature, there are relatively fewer papers studying queueing models with -policy from the perspective of economic. Therefore, a model combining customer retrials, server multiple vacations, and the -policy is of practical interest and is our focus of this paper.
As for studies on equilibrium balking strategies of customers, Burnetas and Economou [1] considered queueing models with setup times under several information levels; Economou and Kanta [4] discussed balking strategies for an observable queue with breakdowns and repairs; Liu, Ma and Li [14] explored an observable queue under single vacation policy; Ma, Liu and Li [15] presented equilibrium balking behavior under a multiple vacation policy; Sun, Li and Cheng-Guo [19] investigated equilibrium strategies and optimal balking strategies for an unobservable queue with double adaptive working vacations. Customers’ equilibrium strategies for queue systems with retrials were also reported in the literature, for example, [12, 23] when balking is not allowed, and [5, 21, 13] when balking is allowed. Regarding models implemented with the -policy, Guo and Hassin [7, 8] investigated models at two information levels with homogeneous and heterogeneous customers, respectively; Guo and Li [9] addressed the same issue for systems, which are partially observable, such as the system capacity is observable, or the state of system is observable. Wang, Zhang and Huang [22] presented customers’ strategic behavior and the social optimal problem in a constant retrial queue with the -policy. Sun, Li and Tian [20] discussed equilibrium strategies and balking strategies with multiple vacations and the -policy.
However, different from the previously mentioned literature on the -policy, the present paper assumes that the system can be reactivated if and only if, upon the completion of a vacation, the server finds at least customers in the virtual orbit; otherwise, the server continues another vacation.
This paper studies equilibrium strategies and optimal balking strategies of customers in a queue with a constant retrial rate, multiple server vacations, and the -policy under two information levels (the observable case and the unobservable case). In this system, there is no waiting area in front of the server and an arriving customer will be serviced immediately if the state of server is idle; otherwise (when the state of the server is busy or on vacation, it has to leave the system to join a virtual retrial orbit waiting for retries according to the FCFS rule. After the completion of each service, the server will take a vacation if the system is empty, or becomes idle if there is at least one customer in the orbit. The idle server will serve the next customer, either a new arrival or a retried customer, whichever comes earlier. The server be reactivated, upon return from the vacation, if at least customers are presented in the virtual orbit; otherwise, the server will start another vacation. For each type of information level, we determine equilibrium strategies and optimal balking strategies of customers and social welfare. For the observable case, in order to ensure that the server can be reactivated, we derive the optimal balking threshold of customers in the vacation state, which must be greater than the optimal threshold in busy state, and also greater than . Therefore, there are three different queuing cases for the observable case, and we study the corresponding stationary distributions for the three queuing cases, and obtain the equilibrium social welfare per time unit. For the unobservable case, we derive the positive equilibrium arrival rate and optimal arrival rate, which are both unique. However, due to the complexity of equations involved, explicit expressions for the equilibrium balking thresholds of customers, socially optimal balking thresholds and optimal social welfare are not available in general. Hence, we use the Particle Swarm Optimization (PSO) algorithm to solve the complex analytic characteristics, by which the numerical optimal solution , and optimal social welfare and are obtained. By comparing the numerical results for the two different information levels, respectively, we conclude that the customers’ equilibrium behavior makes the system more congested than that under the socially optimal strategy, and whether the system information should be disclosed to customers depends on how to maintain the growth of the social welfare (i.e., potential demand arrivals). Obviously, in order to maximize the social welfare, which factor determines the level of information disclosure and when to disclose system information to customers are also crucial for the server or social planner. To conclude our main contributions made in this paper, we emphasize that to our best knowledge, a model, which combines features of retrials, multiple vacations, and the -policy, has not been considered in the literature for the purpose of customers’ equilibrium, and optimal balking strategies.
The remaining sections are organized as follows: In Section 2, we describe the model in detail. We derive the corresponding stationary distributions for the three queueing cases, equilibrium thresholds and social benefit per time unit in Section 3. Section 4 contributes to studies for the unobservable case, and we derive the equilibrium arrival rate and optimal arrival rate. Section 5 focuses on using numerical analysis to explore the theoretical findings in the previous sections, and compare the observable and unobservable cases of this model. Section 6 presents discussions and possible further studies.
2 Model description
Consider a single-server retrial queueing system with a constant retrial rate, multiple server vacations, and the exhaustive -policy. We assume that customers arrive to the system according to a Poisson process with rate , served by a single server with exponential service rate . There is no waiting area in front of the server. An arriving customer will be serviced immediately if the server is idle and leave the system immediately upon the completion of the service; otherwise (the server is either busy or on a vacation), it will join a virtual retrial orbit according to the first-come, first-served discipline (FCFS). In practice, a customer in the orbit can be viewed as a customer on the waiting list. After the completion of a service, the server becomes idle and immediately searches for the customer from the top of the waiting list. The time of the search is a random variable, exponentially distributed with rate . In the search process, if a new customer arrives, the search will be immediately interrupted and the server will return to serving the arriving customer; otherwise (no arrivals during the search process), the customer at the head of the waiting line will be served and will leave the system upon the completion of its service. After all customers in the system are served, or when the system becomes empty, the server will take a vacation of exponential amount of time with rate . During the vacation time, the server will be not available to serve customers. Upon the completion of a vacation, the server will continue to another (independent) vacation with the same parameter if there are fewer than customers in the system; otherwise, the server will return from vacations (to idle state) and immediately start the same search process as that for the case, described above, when the server becomes idle from busy. This type of queue systems is referred to as the M/M/1/MV queue. Inter-arrival times of customers, service times and the times of retrials are assumed to be mutually independent.
The state of the system at time can be represented by a random vector , where denotes the number of customers in the orbit, and denotes the state of the server at time :
Obviously, the stochastic process is a continuous-time Markov chain. The corresponding transition rate diagram is shown in Fig. 1. Moreover, the observable case means that arriving customers can observe all information about and , and the unobservable case implies that arriving customers can not observe any information about and .
The study of strategic customer behavior is important. In our case, we are interested in deciding whether an arriving customers would join or balk the system. Suppose that every customer receives the same reward of units for completing its service, which is used to quantify customer satisfaction or the added value of service. In addition, there is a waiting cost of units per time unit. The total waiting time for a customer is the continuous accumulation of the time when the customer reaches the system and until he leaves the system (including the service time). Customers are risk neutral and want to maximize their expected net benefit. Specifically, if a customer receives the reward of service is more than the expected waiting cost, then he will join the system. If a customer receives the reward of service equal to the expected waiting cost, the customer will be indifferent between joining and balking. Therefore, we only need to consider the reward with satisfying the following inequality:
(1) |
The above constrain assures that all customers, who find server being idle, always enter the system since his reward is more than the waiting cost during his expected service time (). We adopt a natural linear reward-cost structure, and is defined as the expected net benefit after the service completion, i.e., , where is mean sojourn time of the customer. Obviously, if the customer is balking, it will generate .
Under both levels of information levels, under the condition in (1), customers will be sure to enter the system if they find that the state of the system is idle upon arrivals. However, if an arriving customer finds a vacation or busy state, he has to decide whether to leave his contact details (enter retrial orbit) or leave for ever. We further assume that the arriving customers know the policy of the system, i.e., their decisions are irrevocable: the balking customers cannot retry and customers, who joined the system, cannot renege.
3 The observable case
As mentioned above, the observable case means that arriving customers can observe all information about and . Obviously, the information about (orbit length) is useful when the arriving customer finds the server being busy or on the vacation, but it is useless when the customer finds the server being idle upon arrivals, since in this case the customer will be served immediately regardless of the orbit length. Therefore, in all cases, the information about the system (both and ) is valuable for the arriving customer to make a better assessment on whether or not he should join the system. More specifically, if is idle, the customer will join the system for sure regardless of ; if is busy or on vacation, according to the FCFS discipline, the arriving customer knows his position in the orbit, which can help him in deciding whether entering the system is preferable. In the observable case, define to be the sojourn time of the marked customer, who joins the system at state (). For studying optimal balking strategies, we need to consider the expected (residual) net benefit of the marked customer, who is at the th position in the orbit and the state of the server is , after he receives the service. We denote the equilibrium balking threshold of customers and the integrated strategy at sate by and , respectively. In addition, we denote the socially optimal balking threshold of customers and the integrated strategy at state by and , respectively. To characterize and , we first give the following theorem.
Theorem 3.1
For the M/M/1/MV constant retrial queue with multiple vacations and -policy, when a marked customer is at th position in the orbit and the state of server is (), the mean (residual) sojourn time of the marked customer are given by, respectively,
(2) |
(3) |
(4) |
Proof. Consider a marked customer arrived to the system, who found that the server is busy or on vacation. Clearly, the mean overall sojourn time of the marked customer is not affected by customers who arrive after the marked customer by finding the server being busy or on vacation, but it is affected by the customers, who enter the system after the marked customer by finding the server being idle, since in this case, by our imposed condition (1)) they will join the system to receive the service immediately.
Since represents the mean residual service time of the customer, who is receiving the service, we have
(5) |
For , let be the probability of joining the virtual orbit for the arriving customer, who finds the server being busy and customers being in the orbit. Then, based on a first step argument and noticing that the mean time to the next event is , and the next event is an arrival or a service completion with probability or , respectively, we have
(6) |
When the server state is idle, we can similarly get
(7) |
For , we only need expressions for (see Remark 3.1), which is given by
(8) |
In terms of (7) and by solving (6) for , we get
(9) |
which leads to (3). Substituting (3) into (7) produces (4). Finally, substituting (4) into (8) gives (2), which completes the proof.
Remark 3.1
In the above theorem, we did not provide the expression for when , since for our purpose, we only need the expression when the server can be reactivated, or .
3.1 Equilibrium
We first study the equilibrium balking behavior of customers in the observable case, i.e., the customers can observe both information of and at time . As mentioned above, the condition (1) ensures that the customers who find the server is idle always enter the system. The customers who find a vacation or busy state have to decide whether to leave their contact details (enter retrial orbit) or leave for ever. Therefore, we only need to consider that the system is in the state of vacation or busy upon the customers arrivals.
From Theorem 3.1, it is easy to know that the sojourn time of the marked customer satisfies the following equation when he encounters the system state :
(10) |
Define to be the corresponding residual net benefit of the marked customer, and solve to get the equilibrium balking threshold:
(11) |
where the floor function is the largest integer smaller than .
Similarly, when the server state is , we have
(12) |
Define to be the corresponding the residual net benefit of the marking customer, and slove to get the balking threshold:
(13) |
Obviously, there are two possibilities: (i) , which implies ; and (ii) , which implies . Hence, we need to discuss the stationary distribution of the system in three cases: Case 1: ; Case 2: ; and Case 3: for the unobservable case. Our focus in this section is to characterize the integrated balking threshold strategy for these three cases.
Case 1: For , the corresponding transition rate diagram is showed in Fig. 2, and the state space of is given by:
(14) |
Define the stationary distribution as
The stationary distribution for this case (Caese 1) is given in Theorem 3.2.
Theorem 3.2
For the fully observable M/M/1/MV constant retrial queue with multiple vacations and the -policy, if , then the state space of is given by (14), and the stationary distribution is given by:
(15) |
(16) |
and
(17) |
where
(18) |
(19) |
(20) |
(21) |
and
(22) |
can be obtained by the normalization condition .
Proof. From Fig. 2, the corresponding balance equations of the stationary distribution are given as follows,
(23) | ||||
(24) | ||||
(25) | ||||
(26) | ||||
(27) | ||||
(28) | ||||
(29) | ||||
(30) | ||||
(31) |
We first consider the stationary distribution . From (23) and (24), we can obtain
(32) |
(33) |
Based on (26) and (33), we can get
(34) |
Therefore, we can get (15) from the above discussion.
Next, we consider the stationary distribution . From (28) and (30), we can obtain
(35) |
The solution of (35) can be given by the following homogeneous linear difference equation:
(36) |
The characteristic equation corresponding to (36) is
(37) |
which has two roots: 1 and . Let be the general solution of (36), where and are the coefficients that need to be determined. From (27) and (30), we can obtain
(38) |
which yields
(39) |
Therefore,
(40) |
where and are given by (39).
Now, let us continue to consider the stationary distribution . From (28) and (31), we can obtain
(41) |
The solutions of (3.1) can be obtained through solving the following system of nonhomogeneous linear difference equations:
(42) |
whose corresponding characteristic equation is given by
(43) |
Define as the general solution of (43), where is the general solution of the homogeneous version of (43), which is , and is a specific solution of (43).
We consider a specific solution of (43). Substituting it into (43), we can obtain
(44) |
Thus,
(45) |
where and are the coefficients that need to be determined. By considering (40), we get
(46) |
Therefore,
(47) |
where , and are given by (44) and (46), respectively. Specially, based on (28), (31), (15) and (47), we can obtain the stationary distribution of as follows:
(48) |
Based on (28), (30), (31), (15) and (48), we can get the stationary distribution of as follows:
(49) |
Continue our proof for the case of . In this case, the general solution of (36) is , where and are the coefficients that need to be determined. From (29), (30) and (49), we can obtain and
(50) |
Therefore,
(51) |
which leads to (16).
Finally, we consider the stationary distribution of . From (16), (30) and (31), we can easily get (17).
Based on Fig. 2 and Theorem 3.2, we know that the balking states of customers are and . For the social optimization, which will be considered later, we define to be the social benefit per time unit in Case 1: , or
(52) |
Indeed, the first summand of (52) is the effective arrival rate at the system times the reward , while the second summand is the mean number of customers in the system. Obviously, the equilibrium social benefit is .
Case 2: For , the corresponding transition rate diagram is showed in Fig. 3, and the state space of is given by
(53) |
The stationary distribution for this case (Caese 2) is given in Theorem 3.3.
Theorem 3.3
For the fully observable M/M/1/MV constant retrial queue with multiple vacations and the -policy, if , then the state space of is given by (53), and the stationary distribution is given by:
(54) |
(55) |
and
(56) |
where (), (), and are given by (18), (19), (20) and (21), respectively.
(57) |
can be obtained by the normalization condition .
Proof From Fig. 3, the corresponding balance equations of the stationary distribution are given as follows:
(58) | ||||
(59) | ||||
(60) | ||||
(61) | ||||
(62) | ||||
(63) | ||||
(64) | ||||
(65) | ||||
(66) | ||||
(67) | ||||
(68) |
We first consider the stationary distribution . From (58)–(61), the discussion is similar to the discussion for (23)–(26), which leads to (54).
We next consider the stationary distribution . From (62), (63) and (67), the discussion is similar to that for (40), and we can obtain
(69) |
where and are given by (39).
We now continue to consider the stationary distribution . From (54), (63) and (68), the discussion is similar to that for (47), and we can obtain
(70) |
where , and are given by (44) and (46), respectively. Specially, based on (54), (63) and (70), we can get the stationary distribution of as follows:
(71) |
Based on (54), (64), (68) and (71), we can get the stationary distribution of as follows
(72) |
Continue further to consider the stationary distribution of . Based on (54), (65) and (68), we can obtain that
(73) |
Recursively using (73), we can obtain that
(74) |
Thus,
(75) |
where
(76) |
Specially, based on (54), (65) and (68), we can get the stationary distribution of as follows:
(77) |
Therefore, taking into account all the above discussions, we can get (55).
Finally, we consider the stationary distribution of . From (55), (67) and (68), We can easily get (56).
Based on Fig 3 and Theorem 3.3, we know that the states at which customers will balk are and . Denote to be the social benefit per time unit in Case 2, or . Then,
(78) |
Obviously, the equilibrium social benefit is .
Case 3: For , the corresponding transition rate diagram is showed in Fig. 4, and the state space of is given by:
(79) |
The stationary distribution for this case (Caese 3) is given in Theorem 3.4.
Theorem 3.4
For the fully observable M/M/1/MV constant retrial queue with multiple vacations and the -policy, if , then the state space of is given by (79), and the stationary distribution is given by:
(80) |
(81) |
and
(82) |
where (, (, and are given by (18), (19), (20) and (21), respectively,
(83) |
can be obtained by the normalization condition .
From Fig. 4, we know that Case 3 and Case 2 have the same equilibrium equations. Hence, the stationary distribution of Case 3 is the same as that of Case 2, and therefore we omit the proof of Theorem 3.4. Based on Fig 4 and Theorem 3.4, we can obtain that the states, at which customers will balk are and . Denote the social benefit per time unit in Case 3 by . We then have
(84) |
Obviously, the equilibrium social benefit is .
3.2 Social optimization
Summarize the above discussion, we define as the social benefit per time unit, thus
(85) |
We define as the socially optimal social welfare. Obviously, it’s easy for us to get .
4 The unobservable case
In the observable case, we assume that the arriving customers can observe all information about and . Now, in the unobservable case, we assume that the arriving customers cannot observe any information about or . In this case, the customers join the system with probability (). So, the effective arrival rate is , the equilibrium mixed strategy of the customers is denoted by the equilibrium arrival rate ( is equilibrium joining probability of the customers), and the socially optimal mixed strategy is denoted by the optimal arrival rate , where is optimal joining probability of the customers. The corresponding transition rate diagram is showed in Fig. 5.
4.1 Equilibrium
In the unobservable case, the arriving customers can neither observe the state of the server nor the number of other customers in the orbit. In order to obtain the equilibrium arrival rate in this case, the stationary distribution needs to be determined first. From Fig. 5, it is easy to know that the process is a quasi-birth-and-death (QBD) process with Ω_un = { (0,n):n ≥0} ∪{ (1,n):n ≥0} ∪{ (2,n):n ≥1} . If , is defined as the stationary limit of the process . The stationary distribution is denoted by:
where
with the definition of
Let be defined as the infinitesimal generator of the process. Then, the stationary probability vector can be solved through the equations :
(86) | ||||
(87) | ||||
(88) | ||||
(89) | ||||
(90) | ||||
(91) | ||||
(92) |
From Fig. 5 and ordering the states in the state space lexicographically, we can write the infinite generator for the Markov process as a block-partitioned form as follows:
(93) |
where
Theorem 4.1
For the unobservable case, the M/M/1/MV constant retrial queue with multiple vacations and the -policy, given the arrival rate , the stationary distribution is given by
(94) |
(95) |
and
(96) |
where ,
(97) |
(98) |
and
(99) |
Proof In order to obtain the stationary distribution of the system, we first need to obtain the rate matrix , which is the minimum non-negative solution of the following matrix quadratic equation:
(100) |
By detailed calculations, we get the minimum non-negative solution of as follows:
(101) |
Using the matrix-geometric solution (see [17]), we have:
(102) |
and satisfies the following equation:
(103) |
where
(104) |
Substituting (104) into (103), we can get:
(105) |
By calculating (105), we can get
(106) |
(107) |
and
(108) |
where ,
(109) |
From (106)–(108), can be obtained. By (101), can be obtained as follows:
(110) |
where and see (98). Now, from (102), we can get:
(111) |
Hence, considering the above discussion, (94)–(96) can be obtained, and can be calculated by the following normalization condition:
(112) |
where the expression for is given in (99).
In Theorem 4.1, the stationary distribution under unobservable case was obtained by using the matric-analytic method, based on which we can get the mean queue length for the unobservable case, given by:
(113) |
By using Little’s law for the whole system, we can get mean sojourn time :
(114) |
We obtain the equilibrium arrival rate and socially optimal arrival rate by the following theorem.
Theorem 4.2
For the unobservable M/M/1/MV constant retrial queue with multiple vacations and the -policy, we have the following conclusions on the equilibrium arrival rate:
(i) It has no positive equilibrium arrival rate, if ;
(ii) It has one positive equilibrium arrival rate (if and only if ), if , where is the unique positive solution of ;
(iii)
(115) |
where , (with ) are the positive solutions of .
Proof From (114), we can obtain the expected net benefit of the marked customer:
(116) |
Since the second-order derivative of in is given by:
(117) |
If then is positive. In this case, is strictly convex of . If we denote the positive solution of equation by , and the positive solutions of by and (), then we have the following results:
(1) When , i.e., , is negative for every . Therefore, it has no positive equilibrium arrival rate, which leads to (i).
(2) When , i.e., , is negative for every . Therefor, it has one positive equilibrium arrival rate (if and only if ), which leads to (ii).
(3) When , i.e., . If , it has two positive equilibrium arrival rates , which leads to the first part of (115). If , it has two positive equilibrium arrival rates , which leads to the second part of (115). If , it has a unique positive equilibrium arrival rate , which leads to the third part of (115). If , obviously, it has no positive equilibrium arrival rate, which leads to the fourth part of (115).
4.2 Social optimization
In this section, we discuss the optimal balking behavior of customers in the unobservable case in terms of the following Theorem.
Theorem 4.3
For the unobservable M/M/1/MV constant retrial queue with multiple vacations and the -policy, the socially optimal mixed strategy is given by
(118) |
where () is the solution of , and .
Proof From (4.1), we can get the social welfare per time unite as follows:
(119) |
The second-order derivative of in is given by:
(120) |
If , then . So, is strictly concave of . If we denote the unique positive solution of equation by (), then we have the following results:
(1) When , obviously, the socially optimal mixed strategy of customers is unique, which is .
(2) When , the socially optimal mixed strategy of customers is unique, which is .
5 Numerical results
In this section, we explore the previous theoretical results through numerical experiments. One should note that, due to the complexity of equations (85) and (4.2), explicit expressions for the equilibrium balking threshold , , of customers, the socially optimal balking threshold , , and the optimal social welfare are not available in general. Hence, we use Particle Swarm Optimization (PSO) algorithm to numerically solve complex analytic characteristics in this section. The numerical optimal solution of and optimal social welfare and can be obtained by PSO algorithm. PSO algorithm was introduced by Kennedy and Eberhart [3] to solve continuous nonlinear optimization problems, and it has been widely used to solve global optimal solutions, since it does not require many constraints and objective functions. The key procedure of applying PSO algorithm to find the optimal solution (searching for socially optimal balking threshold ) is illustrated in Algorithm 1, where the velocity and position are generally provided by:
(121) |
and
(122) |
where is the number of particles, is the dimension, is a random number in , and are the learning factor and is the inertia factor.
5.1 Numerical results for the observable case
Based on a large number of numerical experiments with a series of parameter choices, we conclude that key qualitative properties are independent of the choice of parameters. To illustrate these properties, we present some exemplary results below. First, we explore the trend in changes for the socially optimal thresholds with respect to and in Fig 6 and Fig 7, respectively. They illustrate the following phenomena.
1. From Fig. 6, we can observe that and increase with , which illustrates that the social planner wants customers to actively join the system with the growth of .
2. In Fig. 6, it is clear when , when . The reason for this is that when , and when . Specially, when . It indicates that is the minimum threshold to ensure server activity. Therefore, when the social planner sets a larger value, the corresponding will be generated.
3. From Fig. 7, we can observe that both and decrease with , which illustrates that customers’ selfishness does not match the wishes of the social planner. Specially, when , which illustrates that the social planner does not want to accumulate many customers during the vocation.
Next, we explore the trend in changes for and with respect to and in Fig 8 and Fig 9, respectively. They reveal the following phenomena.
1. From Fig. 8 and Fig. 9, we can observe that both and increase with respect to and , respectively. It is obvious that the growth rate of is much faster than the growth rate of .
2. always holds as shown in Fig. 8 and Fig. 9, which illustrates that the customers’ individual behavior under the stable equilibrium can lead to system congestion more seriously.
Finally, Fig. 10 shows that the relationship between the optimal social welfare and , and the relationship between the optimal social welfare and . It reveals the following phenomena.
1. In Fig. 10 (a), the optimal social welfare decreases with . The reason is that when becomes larger, more customers will be accumulated, which leads to more waiting costs.
2. In Fig. 10 (b), the optimal social welfare increases with . When becomes bigger, it speeds up the operation of the system, which can then produce the social welfare more effectively. Moreover, when increases to a certain value, the social welfare reaches its maximum and remains stable afterwards.
5.2 Numerical results for the unobservable case
In the unobservable case, is unstable from the equilibrium point of view. Thus, we only study the stable one with , or in the following numerical results. First, we explore the impact of the parameters , , and on the equilibrium arrival rate and the optimal arrival rate in Fig. 11, respectively. They illustrate the following phenomena.
1. From Fig. 11, we can always see that (i.e. ) (i.e. with , , and , respectively.
2. In Fig. 11 (a), (i.e. ) decreases with , whereas (i.e. increases with . The reason for this is that a larger value of reduces the enthusiasm of customers to join the system. However, a larger value of makes the social planner encourages more customers to enter the system.
3. From Fig. 11 (b), (c) and (d), we can observe that the interests of customers and the social planner are coincident with each other with respect to , and . This is because that shortening the vacation time can reduce waiting costs, and increasing the retrial rate and service rate can accelerate switchover of the system back from vacations.
Moreover, Fig. 12 shows that (i.e. ) (i.e. as a function of . It is obvious that the value of has no significance on and when , but increases with . Thus, a higher value of will scare away customers unless is big enough.
If all customers follow the stable equilibrium mixed strategy , then their equilibrium social welfare per time unit can be achieved. Similarly, the optimal social welfare per time unit can also be obtained. Fig. 13 shows the trend in changes for the optimal social welfare per time unit with respect to and , respectively. Obviously, the case shown in Fig. 13 is similar to that shwon in Fig. 10.
5.3 The role of the information level on the equilibrium social welfare and optimal social welfare
An important issue in the strategic customer queuing model is the level of information that the social planner should provide to customers. Fig. 14 explores the trend in changes for the equilibrium social welfare of customers and optimal social welfare of customers under two levels of information. The properties shown here are presentative since the conclusions made are based on comprehensive numerical experiments with a broad choice of system parameters.
1. In Fig. 14 (a), it shows that the equilibrium strategy of customers is balking when is small, since smaller does not intend to activate the system. In this case, both and (i.e. ) have a zero-increase-decrease trend with . However, when , since the arrival rate is small and the number of people in the system is small. In this case, hiding the system information from customers helps to increase the number of customers entering the system, thereby increasing social welfare. Similarly, when , which shows that disclosing the system information can help reduce the system congestion, thus reducing waiting costs.
2. In Fig. 14 (b), it shows that the socially optimal mixed strategy of customers is balking when is small. Both and (i.e. ) keep growth, eventually becomes a constant. Similar to Fig. 14 (a), Fig. 14 (b) also shows that when , when . This also shows that the information level of the system has a serious impact on the social welfare. Therefore, the social planner should choose the strategy consistent with the system designer to achieve social optimum.
6 Conclusions and further research
In this paper, we studied equilibrium strategies and optimal balking strategies of customers in a constant retrial queue with multiple vacations and the -policy under two information levels (observable case and unobservable case), respectively. For each type of information levels, we determined equilibrium strategies and optimal balking strategies of customers, and the social welfare. For the observable case, in order to ensure that the server can be reactivated, we obtained that the optimal balking threshold of customers in the vacation state must be greater than the optimal threshold in busy state or greater than . Therefore, there are three different queuing cases for the observable case, and we obtained the corresponding stationary distributions for the three queuing cases, and determined the equilibrium social welfare per time unit. For the unobservable case, we obtained the positive equilibrium arrival rate and optimal arrival rate, which are unique. In Section 5, we explored the previous theoretical results through numerical experiments. However, due to the complexity of the involved equations, explicit expressions for the equilibrium balking thresholds of customers, socially optimal balking thresholds and optimal social welfare are not available in general. Hence, we use Particle Swarm Optimization (PSO) algorithm to solve the complex analytic characteristics. The numerical optimal solution and optimal social welfare and are obtained by PSO algorithm. By comparing the numerical results of the two information levels, we obtained that the customers’ behavior under the stable equilibrium makes the system more congested than that under the socially optimal one, and whether the system information should be disclosed to customers depends on how to maintain the growth of the social welfare (i.e., potential demand arrivals). Obviously, in order to maximize the social welfare, which factor determines the level of information disclosure and when to disclose the system information to customers are also crucial for the server or social planner. Fortunately, this paper achieved this goal. In the future, it is necessary for us to consider the almost observable case of this model, i.e., the state of the server can be observed, but the number of customers in the orbit cannot be observed.
Acknowledgements
This work was supported in partial by The National Natural Science Foundation of China (No. 61773014), the program of China Scholarships Council (No. 201906840070), and by The Natural Sciences and Engineering Research Council of Canada (NSERC).
Disclosure statement
None of the authors have any competing interests in the manuscript.
References
- [1] Burnetas A, Economou A. Equilibrium customer strategies in a single server markovian queue with setup times. Queueing Systems 2007;56(3-4):213–28.
- [2] Chen H, Frank MZ. State dependent pricing with a queue. Iie Transactions 2001;33(10):847–60.
- [3] Eberhart R, Kennedy J. A new optimizer using particle swarm theory. MHS’95 Proceedings of the Sixth International Symposium on Micro Machine and Human Science 1995;39–43.
- [4] Economou A, Kanta S. Equilibrium balking strategies in the observable single-server queue with breakdowns and repairs. Operations Research Letters 2008;36(6):696–9.
- [5] Economou A, Kanta S. Equilibrium customer strategies and social–profit maximization in the single-server constant retrial queue. Naval Research Logistics 2011;58(2):107–22.
- [6] Edelson NM, Hilderbrand DK. Congestion Tolls for Poisson Queuing Processes. Econometrica 1975;43(1):81–92.
- [7] Guo P, Hassin R. Strategic behavior and social optimization in markovian vacation queues. Operations research 2011;59(4):986–97.
- [8] Guo P, Hassin R. Strategic behavior and social optimization in markovian vacation queues: The case of heterogeneous customers. European Journal of Operational Research 2012;222(2):278–86.
- [9] Guo P, Li Q. Strategic behavior and social optimization in partially-observable markovian vacation queues. Operations Research Letters 2013;41(3):277–84.
- [10] Hassin R, Haviv M. To Queue or Not to Queue: Equilibrium Behavior in Queueing Systems. Boston: Kluwer Academic Publishers, 2003.
- [11] Johansen SG, Stidham S. Control of arrivals to a stochastic input–output system. Advances in Applied Probability 1980;12(4):972–99.
- [12] Kulkarni VG. A game theoretic model for two types of customers competing for service. Operations Research Letters 1983;2(3):119–22.
- [13] Kumar BK, Vijayalakshmi G, Krishnamoorthy A, Basha SS. A single server feedback retrial queue with collisions. Computers and Operations Research 2010;37(7):1247–55.
- [14] Liu W, Ma Y, Li J. Equilibrium threshold strategies in observable queueing systems under single vacation policy. Applied Mathematical Modelling 2012;36(12):6186–202.
- [15] Ma Y, Liu W, Li J. Equilibrium balking behavior in the geo/geo/1 queueing system with multiple vacations. Applied Mathematical Modelling 2013;37(6):3861–78.
- [16] Naor P. The Regulation of Queue Size by Levying Tolls. Econometrica 1969;37(1):15–24.
- [17] Neuts MF. Matrix-geometric solutions in stochastic models. Johns Hopkins University Press, Baltimore, Md, 1981.
- [18] Stidham S. Optimal control of admission to a queueing system. IEEE Transactions on Automatic Control 1985;30(8):705–13.
- [19] Sun W, Li S, Cheng-Guo E. Equilibrium and optimal balking strategies of customers in markovian queues with multiple vacations and n-policy. Applied Mathematical Modelling 2016;40(1):284–301.
- [20] Sun W, Li S, Tian N. Equilibrium and optimal balking strategies of customers in unobservable queues with double adaptive working vacations. Quality Technology and Quantitative Management 2017;14(1):94–113.
- [21] Wang J, Zhang F. Strategic joining in m/m/1 retrial queues. European Journal of Operational Research 2013;230(1):76–87.
- [22] Wang J, Zhang X, Huang P. Strategic behavior and social optimization in a constant retrial queue with the n-policy. European Journal of Operational Research 2017;256(3):841–9.
- [23] Zhang F, Wang J, Liu B. On the optimal and equilibrium retrial rates in an unreliable retrial queue with vacations. Journal of Industrial and Management Optimization 2012;8(4):861–75.