CC BY
123
Accelerated Life Testing Design Using Geometric Process for Generalized Rayleigh Distribution with Complete Data
Kamal Ullah, Inthekhab Alam, Showkat Ahmad Lone*
•
Department of Statistics & Operations Research, Aligarh Muslim University, Aligarh. Email: [email protected]
Abstract
The log-linear function between life and stress which is just a simple re-parameterization of the original parameter of the life distribution is used to obtain the estimates of original parameters in many of the studies concerning Accelerated life testing (ALT). But from the statistical point of view, it is preferable to work with the original parameters instead of developing inferences for the parameters of the log-linear link function. In this study we introduce the geometric process for the analysis of accelerated life testing with Generalized Rayleigh Distribution for constant stress. Assuming that the lifetimes of units under increasing stress levels form a geometric process, the maximum likelihood estimation approach is used for the estimation of parameters. The confidence intervals (CIs) of the model parameters are derived. A Simulation study is also performed to check the statistical properties of estimates of the parameters and the confidence intervals.
Keywords: Geometric process, Generalized Rayleigh Distribution, Maximum Likelihood Estimator, Fisher Information Matrix, Asymptotic Confidence Interval, Simulation Study.
1. Introduction
Accelerated life testing is the process of testing a product by subjecting it to conditions (stress, strain, temperatures, voltage, vibration rate, pressure etc.) in excess of its normal service parameters in an effort to uncover faults and potential modes of failure in a short amount of time. By analyzing the product's response to such tests, statisticians can make predictions about the service life and maintenance intervals of a product.
In general, ALT deals with three types of stress patterns: constant stress, step stress and Progressive stress. In the former case, each unit is run at a pre-specified constant stress level which does not vary with time. This means that every item is subjected to only one stress level until the item fails or the test is stopped for other reasons. In use, most products such as semiconductors and microelectronics, capacitors, lamps ...etc, run at a constant stress. This type of stress is widely used and preferred because the stress is constant in most applications, it is much easier to apply and quantify constant stress and models for constant stress are available, widely publicized and empirically verified.
There is a lot of literature on constant stress accelerated life testing, for example, Ahmad et al. [1], Islam and Ahmad [2], Ahmad and Islam [3], Ahmad et al.[4] and Ahmad [5] discuss the optimal constant stress accelerated life test designs under periodic inspection and Type-I censoring. Yang [6] proposed an optimal design of 4-level constant stress ALT plans considering different censoring
times. Pan et al. [7] proposed a Bivariate constant stress accelerated degradation test model by assuming that the copula parameter is a function of the stress level that can be described by the logistic function. Wilkins and johns [8] considered constant stress accelerated life test based on Weibull distribution with constant shape and a log-linear link between scale the stress factor which is terminated by a Type-II censoring regime at one of the stress levels.
The concept of geometric process in accelerated life testing was first introduced by Lam [9] in the problems of repair replacement. Lam [10] studied the geometric process model for a multistate system and concluded a replacement policy to minimize the long run average cost per unit time. Since then a lot of studies in maintenance problems and system reliability have been shown that a GP model is a good and simple model for analysis of data with a single trend or multiple trends, for example, Lam and Zhang [11], Lam [12] and Zhang [13]. Huang [14] introduced the GP model for the analysis of constant stress ALT with complete and censored exponential samples. Kamal et al. [15] extended the GP model for the analysis of complete Weibull failure data in constant stress ALT. Zhou et al. [16] implement the GP in ALT based on the progressive Type-I hybrid censored Rayleigh failure data. Kamal et al. [17] used the geometric process for the analysis of constant stress accelerated life testing for Pareto Distribution with complete data. S. Saxena [18] introduces the Rayleigh geometric process model for the analysis of accelerated life testing under constant stress. Sadia Anwar et al. [19] presented the mathematical model of accelerated life testing for Marshall-Olkin extended exponential distribution using geometric process and extended her work using type I censored data [20]. Recently Kamal [21] presented an application of the geometric process in accelerated life testing analysis on type-I censored Weibull failure data.
In the present study, the GP model is implemented in the analysis of ALT for the Generalized Rayleigh life distribution under constant stress with complete data. Maximum likelihood (ML) estimates of parameters and their asymptotic confidence intervals (CIs) are obtained. The performance of the estimates is evaluated by a simulation study.
2. The Model and Test Procedure
2.1. The Geometric Process
A geometric process describes a stochastic process {xn,n = 1,2,...}, where there exists a real-valued X> 0 such that {¿n-1 x„,n = 1,2,...} forms a renewal process. It can be shown that if {Xn,n = 1,2,...} is a GP and the probability density function of Xx is f (x) with mean ¡1 and
variance a2 then the probability density function of Xn will be X 1 f (X 1 x) with E(Xn ) = —d
X
2
and var(XK ) = ~a—r . Thus X , ¡1 and a2 are three important parameters of GP.
X( )
2.2. The Generalized Rayleigh Distribution
The probability density function (pdf) of a generalized Rayleigh distribution is given by
f(xla,¡)=<!(1 - e ^ ) x > 0,a> 0,¡> 0 [0, elsewhere
where, a > 0 is the shape parameter and 5 > 0 is the scale parameter of the distribution.
Generalized Rayleigh distribution is a member of the family of Burr distributions which was appeared since 1942. It is known also Burr type X distribution. The cumulative distribution function (cdf) of generalized Rayleigh distribution is
F(x la, 0) = e PX ) ,x>0 ,a> 0 ,p > 0 (0 , elsewhere
The Hazard function of the Generalized Rayleigh distribution takes the following form
5(x la,0) = 1 -(l - e-0x2 )" The failure rate (or hazard rate) for the Generalized Rayleigh distribution is given by
h(x l a, 0) =
2 I a. 2 \a-1
2a0xex (l - e ^ ) 1 -(l-e0 )a
The two-parameter Generalized Rayleigh distribution was first proposed by (Raqab and Kundu; 2003) [22] and is denoted by GR(a, 0) ■ It is observed that the hazard function of a Generalized Rayleigh distribution can be either bathtub type or increasing function, depending on the shape
parameter a . For a < 1, the hazard function is bathtub type and for a > — , it has an increasing
hazard function. Surles and Padgett (2001) [23] showed that the two-parameter GR distribution can be used quite effectively in modelling strength data and also modelling general lifetime data.
2.3. Assumptions and test procedure
1. Under any constant stress, the time to failure of test unit follows Generalized Rayleigh distribution.
a
2. The Generalized Rayleigh shape parameter ^ is constant, i.e. independent of stress.
3. Let the sequence of random variables X0,Xx ,..XS denote the lifetimes under each stress level,
X ix k = 12 s)
where 0denotes lifetime of an item under the design stress. We assume ^ k' ' ''"' ' is a geometric process with ratio X > 0 .
4. Suppose that an ALT under zk, k = 1, 2,..., s, arithmetically increasing stress levels is performed. A
random sample of Ni' * =1'2'."'n, identical items are placed under each stress level and start to operate at the same time. Whenever an item fails, it is removed from the test and its observed
failure time Xki is recorded.
5. The scale parameter is a log-linear function of stress that is log0 = a + bSt, here a and b are unknown parameters depending on the nature of the product and the test method.
Theorem: If the stress level in an ALT is increasing with a constant difference then under each stress level the lifetimes of items forms a GP. That is, If Sk+1 - Sk is constant for k = 1'2'...,s -1, then X,k = 1,2,...,s) forms a GP. Proof: From assumption (5), we get
(B ^
'og(t ,
V Pk y
This shows that the increased stress levels form an arithmetic sequence with a constant difference
AS.
Now the above equation can be written as
= b(Sk+1 - Sk ) = bAS
All = = ¿(say) (1)
pk
It is clear from (1) that
P = Xpk 1 = £pk 2 =...= Xp
The lifetime PDF of an item at the kth stress level is fXi (x) = 2aP2 xe )2 (i _ e -(pkx)2 jT
= 2a(Xp)2 xe -(xp)2 (1 _ e -(xp? J11 (2)
This implies that
fXk (x) = XkfXQ (Xx) (3)
Now, from the definition of GP and from expression (3) it is clear that, if density function of X 0
is fx0 (x), then the pdf of X^ will be given by ^ fX {fix) , k = 1,2,...,s. Therefore, it is clear
that lifetimes under a sequence of arithmetically increasing stress levels form a GP with ratio X . Now, the pdf of a lifetime of an item at the k th stress level is
fXk (x/a,PX) = = 2a(lkp)2xeJ* (1 _ e^^ (4)
It is clear from above expression that if lifetimes of items under a sequence of increasing stress level form a geometric process with ratio X and if the life distribution at design stress level is
P ,th
generalized Rayleigh with characteristic H , then the life distribution at k stress level will also be generalized Rayleigh with characteristic life pX .
2.4. Maximum likelihood Estimation
The likelihood function for constant stress ALT for complete case generalized Rayleigh distribution failure data using GP for s stress levels is given by:
L(X a, p) = nft 2a(Xp) x, e^ F (i _ e^ f (5)
k=1 i=1
The log likelihood of (5) can be written a
log 2a l2k log X 1 2 log p 1 log % _ (Xpxh )2 l
l X,a, p) = ZZ
s n k=1 i =1
(a _ 1) log( 1 _ e^
Partial derivatives of above equation with respect to X and p are:
§ =£i {f_2kXi_1 pxty 12kXk_1 (a_OK) 1} (6)
&¡_ dß
= t t 1-2 - )2 + ßa- ^ )
k=1 i=1
ß
A D
(7)
Where
A= e ^ r and D= fl - e ^?
From equations (6) and (7), it is observed that these equations are non-linear. Therefore, the closed forms of MLEs of X and ß do not exist. So, Newton-Raphson method must be used to solve these equations simultaneously to obtain the MLEs of X and ß .
3. Asymptotic Confidence Interval
Let I(X, ß) denotes the Fisher Information matrix, then observed Information matrix of I(X, ß) is given as
I (X, ß) =
Where
111 1 12
121 122
Ill =-
fd 2l^
y dX ;
= tt
k=1 i=1
2k + 2k (2k - l)Z 2X2(k-l)+(a-1) X
f 4ADZ4k2X2(k-1) - 2kAD(2k - 1)Z2X2(k-1) + 4k2A2Z4X2(2k-1
D2
I = -
1 22
f d2l
dß2
=tt
k=1 i=1
4ADß2Z4 - 2ADZ2 + 4A2ß2Z4
D2
— + 2Z2 +(a-1
I12 =-
f d2l >
ydXdßj
= 121 =-
f d2l > ydßdX ;
=tt
k=1 i=1
4ßkZ2 X
+ 4
a
-1 M ADkß3Z4 - ADkßZ2 + Z4A2ß3k
X
y X ;
D2
Where
Z = (Xkx. )
Now, the variance-covariance matrix can be written as
-1-1
var(X) " I" "
var(ß) _I 22 _
The 100(1 asymptotic confidence interval for X and P are then given respectively as
s n
And
X ± z ^ var (X)
111
p±z ¿yfva
4. Simulation Studies:
Simulation of data is the initial task for studying different properties of parameters. It is an attempt to model an assumed condition to study the behaviour of the function.
1. First, to perform the simulation study, first, a random sample is generated from Uniform distribution by using R software.
2. Now, we use inverse cdf method to transform the cdf at kth stress level in terms of u and get the
expression of Xki, k = 1,2,...,s; i = 1,2,...,n. 1
Xki =_
In (1 _ u)
pXk
2a
k = 1,2,...,s; i = 1,2,..., n.
Where X^ is obtained for n=20, 40 and 60.
3. The values of parameters and numbers of the stress levels are chosen to be a = 1,p = 2.8,X = 1.1 and s = 4 or 6.
4. By using optim() function, we obtain ML estimates, the mean squared error (MSE), relative absolute bias (RAB), relative error (RE) and lower and upper bound of 95% and 99% confidence intervals for different sample sizes n=20,40 and 60. The results obtained in the above simulation study are summarized in Table1 & 2.
Table 1: Simulation results of Generalized Rayleigh distribution using GP at a = 1, p = 2.8, X = 1.1 and s = 4 .
2
Sample Estimate Mean SE VMSE LCL UCL
20 (3 3.078 0.319 0.095 2.452 2.254 3.703 3.901
A 1.107 0.103 0.099 0.797 0.732 1.202 1.267
40 3.039 0.256 0.061 2.536 2.377 3.542 3.701
A 1.081 0.103 0.100 0.797 0.732 1.202 1.267
60 3.003 0.241 0.054 2.529 2.380 3.477 3.627
A 1.072 0.103 0.099 0.797 0.733 1.202 1.267
Table 2: Simulation results of Generalized Rayleigh distribution using GP at a = 1,/ = 2.8,1 = 1.1 and s = 4.
Sample Estimate Mean SE VMSE LCL UCL
20 (3 3.078 0.218 0.078 2.491 2.311 3.628 3.807
A 0.977 0.112 0.103 0.797 0.732 1.202 1.267
40 3.039 0.209 0.074 2.520 2.348 3.609 3.781
A 0.981 0.103 0.101 0.797 0.732 1.202 1.267
60 2.953 0.192 0.044 2.551 2.416 3.407 3.543
A 0.992 0.020 0.100 0.784 0.732 1.202 1.267
5. Conclusions
In this study, the geometric process is introduced for the analysis of accelerated life testing under constant stress when the life data are from a generalized Rayleigh distribution. It is a better choice for life testing because of its simplicity in nature. The Mean, SE and RMSE of the parameters are obtained. Based on the asymptotic normality, the 95% and 99% confidence intervals of the parameters are also obtained.
The results show in Table 1 and Table 2 that the estimated values of p and A are very close to true (or initial) value with very small SE and RMSE. As sample size increases, the value of SE and RMSE decreases and the confidence interval become narrower. For the Table 2, the maximum likelihood estimators have good statistical properties than the Table 1 for all sample sizes.
References
[1] Ahmad, R. K. (1994). Optimal design of accelerated life test plans under periodic inspection and type I censoring: the case of Rayleigh failure law. South African Statistical Journal, 28(2), 93101.
[2] Islam, A., & Ahmad, N. (1994). Optimal design of accelerated life tests for the Weibull distribution under periodic inspection and type I censoring. Microelectronics Reliability, 34(9), 1459-1468.
[3 Ahmad, N., & Islam, A. (1996). Optimal accelerated life test designs for Burr type XII distributions under periodic inspection and type I censoring. Naval Research Logistics, 43(8), 1049-1077.
[4] Ahmad, N., Islam, A., & Salam, A. (2006). Analysis of optimal accelerated life test plans for periodic inspection: The case of exponentiated Weibull failure model. International Journal of Quality & Reliability Management, 23(8), 1019-1046.
[5] Ahmad, N. (2010). Designing accelerated life tests for generalised exponential distribution with log-linear model. International Journal of Reliability and Safety, 4(2-3), 238-264.
[6] Yang, G. B. (1994). Optimum constant-stress accelerated life-test plans. IEEE Transactions on Reliability, 43(4), 575-581.
[7] Pan, Z., Balakrishnan*, N., & Sun, Q. (2011). Bivariate constant-stress accelerated degradation model and inference. Communications in Statistics—Simulation and Computation®, 40(2), 247257.
[8] Watkins, A. J., & John, A. M. (2008). On constant stress accelerated life tests terminated by Type II censoring at one of the stress levels. Journal of Statistical Planning and Inference, 138(3), 768-786.
[9] Lin, Y. L. Y. (1988). Geometric processes and replacement problem. Acta Mathematicae Applicatae Sinica, 4, 366-377.
[10] Yeh, L. (2005). A monotone process maintenance model for a multistate system. Journal of Applied Probability, 42(1), 1-14.
[11] Lam, Y., & Zhang, Y. L. (1996). Analysis of a two-component series system with a geometric process model. Naval Research Logistics (NRL), 43(4), 491-502.
[12] Yeh, L. (2005). A monotone process maintenance model for a multistate system. Journal of Applied Probability, 42(1), 1-14.
[13] Zhang, Y. L. (2008). A geometrical process repair model for a repairable system with delayed repair. Computers & Mathematics with Applications, 55(8), 1629-1643.
[14] Huang, S. (1911). Statistical inference in accelerated life testing with geometric process model (Doctoral dissertation, Sciences).
[15] Kamal, M., Zarrin, S., & Saxena, S. (2012). Weibull Geometric Process Model for the Analysis of Accelerated Life Testing with Complete Data. International Journal of Statistics and Applications, 2(5), 60-66.
[16] Zhou, K., Shi, Y. M., & Sun, T. Y. (2012). Reliability analysis for accelerated life-test with progressive hybrid censored data using geometric process.
[17] Kamal, M., Zarrin, S., & Saxena, S. (2012). Weibull Geometric Process Model for the Analysis of Accelerated Life Testing with Complete Data. International Journal of Statistics and Applications, 2(5), 60-66.
[18] Saxena, S., Zarrin, S., & Kamal, M. (2012). Optimum Step Stress Accelerated Life Testing for Rayleigh Distribution. International journal of statistics and applications, 2(6), 120-125.
[19] S. Anwar, M. Kamal and A. Islam, "Mathematical model of accelerated life testing using geometric process for Marshall-Olkin extended exponential distribution" , International Journal of Innovative Research in Science, Engineering and Technology, vol.2, no.12, pp. 7382-7390, 2013
[20] Anwar, S., Shahab, S., & Islam, A. U. (2014). Accelerated Life Testing Design Using Geometric Process For Marshall-Olkin Extended Exponential Distribution With Type I Censored Data. International Journal of Scientific & Technology Research, 3(1), 179-186.
[21] Kamal, M. (2013). Application Of Geometric Process in Accelerated Life Testing Analysis With Type-I Censored Weibull Failure Data. Reliability: Theory & Applications, 8(3).
[22] Kundu, D., & Raqab, M. Z. (2005). Generalized Rayleigh distribution: different methods of estimations. Computational statistics & data analysis, 49(1), 187-200.
[23] Surles, J. G., & Padgett, W. J. (2001). Inference for reliability and stress-strength for a scaled Burr Type X distribution. Lifetime Data Analysis, 7(2), 187-200.
Time-Dependent Analysis of a Single-Server Queuing Model with Discouraged Arrivals and Retention of
Reneging Customers
Rakesh Kumar* & Sapana Sharma •
Department of Mathematics, Shri Mata Vaishno Devi University, Katra Jammu and Kashmir, India-182320 Email: rakesh stat [email protected], [email protected] * Corresponding Author
Abstract
In this paper, a finite capacity Markovian single-server queuing system with discouraged arrivals, reneging, and retention of reneging customers is studied. The time-dependent probabilities of the queuing system are ob- tained by using a computational technique based on the 4th order Runge- Kutta method. With the help of the time-dependent probabilities, we develop some important measures of performance of the system, such as expected system size, expected reneging rate, and expected retention rate. The time-dependent behavior of the system size probabilities and the ex- pected system size is also studied. Further, the variations in the expected system size, the expected reneging rate, and the expected retention rate with respect to the probability of retaining a reneging customer are also studied. Finally, the effect of discouragement in the same model is ana- lyzed.
Keywords: time-dependent analysis, single server queuing system, discouraged arrivals, reneging, Runge-Kutta method, retention
Introduction
Queuing systems are used in the design and analysis of computer-communication networks, production systems, surface and air traffic systems, service systems etc. The study of queueing systems help to manage waiting lines and to construct an optimal system for balancing customer waiting time with the idle time of the server Gnedenko and Kovalenko (1989). The enormous literature in queuing theory is available where the customers always wait in the queue until their service is completed. But in many practical situations customers become impatient and leave the systems without getting service. Therefore, queuing systems with customers' impatience have attracted a lot of attention. The study of customers' impatience in queueing theory is started in the early 1950's. Haight (1959) studies a single-server queue in steady-state with a Poisson input and exponential holding time, for various reneging distributions. Ancker and Gafarian [(1963a), (1963b)] analyze an M/M/l/N queuing system with balking and reneging. In addition, the effect of reneging on an M/M/l/N queue is investigated in the works of Abou El-Ata (1991), Zhang et al. (2006), Al Seddy et al. (2009), and Wang and Chang (2002). Kovalenko (1961) discusses some queuing systems with restrictions.
Queuing systems with discouraged arrivals are widely studied due to their significant role in managing daily queueing situations. In many practical situations, the service facility possesses
defense mechanisms against long waiting lines. For instance, the congestion control mechanism prevents the formation of long queues in computer and communication systems by controlling the transmission rates of packets based on the queue length(of packets) at source or destination. Moreover, a long waiting line may force the servers to increase their rate of service as well as discourage prospective customers which results in balking. Hence, one should study queueing systems by taking into consideration the state-dependent nature of the system. In state-dependent queues the arrival and service rates depend on the number of customers in the system. The discouragement affects the arrival rate of the queueing system. Customers arrive in a Poisson fashion with rate that depends on the number of customers present in the system at that time i.e.
. Morse (1958) considers discouragement in which the arrival rate falls according to a negative exponential law. Natvig (1974), Van Doorn (1981), Sharma and Maheswar (1993), and Parthasarathy and Selvaraju (2001) have also studied the discouraged arrivals queuing systems. Ammar et. al (2012) derive the transient solution of an M/M/1/N queuing model with discouraged arrivals and reneging by employing matrix method. Abdul Rasheed and Manoharan (2016) study a Markovian queueing system with discouraged arrivals and self-regulatory servers. They discuss the steady-state behavior of the system. Rykov (2001) considers a multi-server controllable queuing systems with heterogeneous servers. He studies several monotonicity properties of optimal policies for this system. Koba and Kovalenko (2002) study retrial queuing systems which are used in the analysis of aircraft landing process. Efrosinin and Rykov (2008) study a multi-server heterogeneous queuing system and obtain its steady-state solution. They derive the waiting and sojourn time distributions. They also study the optimal control of the queuing system. Rykov (2013) generalizes the slow server problem to include additional cost structure. He finds that the optimal policy for the problem has a monotone property. Sani et al. (2017) perform the reliability analysis of a system subjected to deterioration before failure. They use system state transition probabilities to derive the Markov models of the system.
Queuing systems with customers' impatience have negative impact on the performance of the system, because it leads to the loss of potential customers. Kumar and Sharma (2012a) take this practically valid aspect into account and study an M/M/l/N queuing system with reneging and retention of reneging customers. Kumar (2013) obtains the transient solution of an M/M/c queue with balking, reneging and retention of reneging customers. Kumar and Sharma (2014) obtain the steady-state solution of a Markovian single server queueing system with discouraged arrivals and retention of reneging customers by using iterative method.
The steady-state results do not reveal the actual functioning of the system. Moreover, stationary results are mainly used within the system design process and it cannot give insight into the transient behavior of the system. That is why, we extend the work of Kumar and Sharma (2014) in the sense that the time-dependent analysis of the model is performed. The time-dependent numerical behavior is studied by using a numerical technique Runge-Kutta method.
1 Queuing Model Description
In this section, we describe the queueing model. The model is based on following assumptions:
1. We consider a single-server queueing system in which the customers arrive in a Poisson fashion with rate that depends on the number of customers present in the system at that time i.e. .
n+1
2. There is single server and the service time distribution is negative exponential with parameter n.
3. Arriving customers form a single waiting line based on the order of their arrivals and are served according to the first-come, first-served (FCFS) discipline.
4. The capacity of the system is finite (say N).
5. A queue gets developed when the number of customers exceeds the number of
servers, that is, when n > l. After joining the queue each customer will wait for a certain length of time T (say) for his service to begin. If it has not begun by then he may get renege with probability p and may remain in the queue for his service with probability q (= 1 — p) if certain customer retention strategy is used. This time T is a random variable which follows negative exponential distribution with parameter %. The reneging rate is given by
_ f0- 0 <n<1
= {(n — n>2
2 Mathematical Model
Let {X(t), t > 0} be the number of customers present in the system at time t. Let Pn(t) = P{X(t) = n},n = 0,1,... be the probability that there are n customers in the system at time t. We assume that there is no customer in the system at t = 0.
The differential-difference equations of the model are:
^^¡¡^)=—AP0(t)+^P1(t), (1)
(M + nfp)Pn+i(t),1 <n<N (2)
dPN(t) _ fX
dt
= g) PN-i(t) — &+(N — 1)Sp)PN(t), (3)
3 Transient analysis of the model
In this section, we perform the time-dependent analysis of a finite capacity single-server Markovian Queuing model with discouraged arrivals and retention of reneging customers using Runge-Kutta method of fourth order (R-K 4). The "ode45" function of MATLAB software is used to find the time-dependent numerical results corresponding to the differential-difference equations of the model.
We study the following performance measures in transient state:
1. Expected System Size (Ls(t))
Ls(t) = nPn(t)
2. Average Reneging Rate (Rr(t))
Rr(t) = Zn=1 (n — l)tpPn(t)
3.Average Retention Rate (Rr(t))
Rn(t) = !Z=i (n — 1)fqPn(t)
Now, we perform the time-dependent numerical analysis of the model with the help of a numerical example. We take N = 10, A = 2,^ = 3,% = 0.1, and p = 0.4. The results are presented in the form of Figures 1-5. Following are the main observations:
In Figure 1, the probabilities of number of customers in the system at different time points are plotted. We observe that the probability values P1(t), P2(t), ..., P10(t) increase gradually until they reach stable values except the probability curve P0(t) which decreases rapidly in the beginning and then attains steady-state with the passage of time.
Figure 2 shows the effect of the probability of retaining a reneging customer on the expected system size in transient state. One can observe that as the probability of retaining a reneging customer increases, the expected system size also increases. This establishes the role of probability of retention associated with any customer retention strategy.
In Figure 3, the change in average reneging rate with the change in probability of retaining
a reneging customer is shown. One can observe that there is a proportional decrease in average reneging rate with the increase in probability of retention, q.
The variation in average retention rate with probability of retention is shown in Figure 4. We can see that there is a proportional increase in RR(t) with increase in q, which justifies the functioning of the model.
In figure 5, the impact of discouraged arrivals on the performance of the system is shown. We compare two single server finite capacity Markovian queuing systems having retention of reneging customers with and without discouraged arrivals. One can see from Figure 5 that the expected system size is always lower in case of discouraged arrivals as compare to the queuing model without discouragement.
Time vs. Probabilities
1 r-r-1-1-1-1-1-1-1-
0.9 - \ 0.B - \ 0.7 - \
Time
Figure 1: The probability values for different time points are plotted for the case N = 10, A = 2,p = 3,% = 0.1, and p = 0.4
0 0.2 0.4 0.6 0.8 1
Probability of retention (q)
Figure 2: The expected system sizes versus probability of retention (q) are plotted for the case N = 10, A = 2,p = 3,$ = 0.1, t = 0.5, and q = 0.1,0.2, ...,0.9
Figure 3: Variation of average reneging rate with the variation in probability of retention for the case N = 10, A = 2,n = 3,% = 0.1, t = 0.5, and q = 0.1,0.2, .,0.9
Figure 4: Variation of average retention rate with the variation in probability of retention for the case N = 10, A = 2,/u = 3,{ = 0.1, t = 0.5, and q = 0.1,0.2, .,0.9
-Single server queuing system with discouraged arrivals and
retention of reneging customers
-Single server queuing system with reneging and retention of
reneging customers
1.6 -
1.4 -
S 1-2 -i/i
£ 1 -«
¿T 0.8
*u
€ o.6
Q.
UJ 0.4 0.2 -0
0 0.5 1 1.5 2 2.5
Time
Figure 5: The impact of discouragement on expected system size
4 Conclusions
The time-dependent analysis of a single-server queuing system with discouraged arrivals, reneging and retention of reneging customers is performed by using Runge Kutta method. The numerical results are computed with the help of MATLAB software. The effect of probability of retaining a reneging customer on various performance measures is studied. We also study the impact of discouraged arrivals on the system performance.
Acknowledgements
One of the authors (Dr. Rakesh Kumar) would like to thank the UGC, New Delhi, India for financial support given to him for this research work under the Major Research Project vide letter no. F.-43-434/2014(SR).
References
[1] Abdul Rasheed, K.V., Manoharan, M. (2016). Markovian queueing system with discouraged arrivals and self-regulatory servers. Advances in Operations Research, -Vol. 2016, Article ID 2456135, -P. 11 pages, doi:10.1155/2016/2456135.
[2] Abou-El-Ata, M.O. and Hariri, A.M.A. (1992). The M/M/c/N queue with balking and reneging. Computers and Operations Research, 19:713-716.
[3] Al-Seedy, R.O., El-Sherbiny, A.A., El-Shehawy, S.A, Ammar, S.I. (2009). Transient solution of the M/M/c queue with balking and reneging. Computers and Mathematics with Applications, 57:1280-1285.
[4] Ammar, S.I., El-Sherbiny, A.A., Al-Seedy, R.O. (2012). A matrix approach for the transient solution of an M/M/1/N queue with discouraged arrivals and reneging, International Journal of Computer Mathematics, 89:482-491.
[5] Ancker. Jr., C.J., Gafarian A.V. (1963a). Some queuing problems with balking and reneging I. Operations Research, 11:88-100
[6] Ancker. Jr., C.J., Gafarian A.V. (1963b). Some queuing problems with balking and reneging II. Operations Research, 11:928-937.
[7] Efrosinin, D., Rykov, V. (2008). On performance characteristics for queueing systems with heterogeneous servers. Automation and Remote Control, 69:61-75.
[8] Gnedenko, B.V., Kovalenko, I. N. (1989). Introduction to Queuing Theory. Birkhauser, 2nd Edition.
[9] Haight, F. A. (1959). Queuing with reneging. Metrika, 2:186-197.
[10] Koba, E. V., Kovalenko, I. N. (2002). Three Retrial Queueing Systems Representing Some Special Features of Aircraft Landing. Journal of Automation and Information Sciences, 34:DOI • 10.1615/JAutomatInfScien.v34.i4.10, 4 pages.
[11] Kovalenko, I. N. (1961). Some Queuing Problems with Restrictions. Theory of Probability and Its Applications, 6:204-208.
[12] Kumar, R. and Sharma, S. K. (2012a). An M/M/1/N queuing system with retention of reneged customers. Pakistan Journal of Statistic and Operation Research, 8:859-866.
[13] Kumar, R. and Sharma, S. K. (2014). A single-server Markovian queuing system with discouraged arrivals and retention of reneged customers. Yugoslav Journal of Operations Research, 24:119-216.
[14] Kumar, R. (2013). Economic analysis of an M/M/c/N queuing model with balking, reneging and retention of reneged customers. Opsearch, 50:383-403.
[15] Morse, P.M. Queues, inventories and maintenance. Wiley, New York, 1958.
[16] Natvig, B. (1974). On the transient state probabilities for a queueing model where potential customers are discouraged by queue length. Journal of Applied Probability, 11:345-354.
[17] Parthasarathy, P.R. and Selvaraju, N. (2001). Transient analysis of a queue where potential customers are discouraged by queue length. Mathematical Problems in Engineering, 7:433454.
[18] Rykov, V. (2001). Monotone Control of Queueing Systems with Heterogeneous Servers. Queueing Systems, 37391-403.
[19] Rykov, V. (2013). On a Slow Server Problem. Stochastic Orders in Reliability and Risk,
351-361.
[20] Sani, B., Gatawa, R. I. and Yusuf, I. (2017). Reliability Assessment of Deteriorating System. Reliability: Theory and Applications, 12:20-29.
[21] Sharma, O.P., Maheswar, M.V.R. (1993). Transient behavior of a simple queue with discouraged arrivals. Optimization, 27:283-291.
[22] Van Doorn, E. A. (1981). The transient state probabilities for a queueing model where potential customers are discouraged by queue length. Journal of Applied Probability, 18:499-506.
[23] Wang, K.H., Chang, Y.C. (2002). Cost analysis of a finite M/M/R queueing system with balking, reneging and server breakdowns. Mathematical Methods of Operations Research 56:169-180.
[24] Zhang, Y., Yue, D., Yue, W. (2006). Optimal performance analysis of an M/M/1/N queue system with balking, reneging and server vacation. International Journal of Pure and Applied Mathematics, 28:101-115.
