Modification of the asymptotic analysis method under heavy load condition on the example of the study of the retrial queueing system M vertical bar M vertical bar 1

In the paper, three methods of the asymptotic analysis modification are offered for retrial queueing systems under heavy load condition as mathematical models of telecommunication systems. The characteristic feature of these models is existence of repeated attempts of calls to get service after a random time. There are a large number of papers devoted to retrial queues researching, but mainly problems are solved by numerical methods and simulation. Analytical results are obtained only for systems with Poisson input process. Earlier we proposed the method of asymptotic analysis (first and second order) for retrial queueing systems researching under heavy load condition [13, 14]. It has been shown that the first order asymptotics has the form of gamma distributions for all types of the system, but its application range is quite narrow. The asymptotics of second order has allowed to extend the application range of the results by 4 times, but the asymptotic formula has not unified form and it needs time-taking rather laborious formula derivations for each system. Therefore, we suggest ways to modify the method of the asymptotic analysis under heavy load condition to increase the method accuracy. The method is described on the example of the simplest retrial queue M vertical bar M vertical bar 1 in order to compare the results with the known exact distribution. The system of differential Kolmogorov equations for Markov process {k(t), i(t)} is composed, then it is rewritten for characteristic functions in steady state. The parameter rho=lambda/mu is the load rate of the system. The principle of asymptotic analysis method under heavy load condition (where rho up arrow 1) is briefly described. Asymptotic characteristic functions of the first and second order are presented. Note that during the derivation of the first order asymptotic formula the function F-0(w) was obtained. We propose to study "intermediate" asymptotics h*(u), which has the form of a weighted sum of two gamma distribution characteristic functions: h*(u) = q Gamma(1)(u)+(1-q)Gamma(0)(u), where Gamma(1) (u)(1- ju/beta)(-alpha), Gamma(0)(u)=(1-ju/beta)(-alpha+1), q=1/2-rho, alpha= mu/sigma+1 beta=1-rho We suggest to call this distribution as hypergamma distribution (like a hyperexponential distribution). The numerical analysis shows that hypergamma asymptotics is more accurately than the first order asymptotics, but it is worth than the second order asymptotics. However, the derivation of the second order asymptotic formula for more complex retrial queues is difficult, whereas "intermediate" asymptotics does not require additional calculations, because its shape is the same for all retrial queues and hypergamma distribution parameters are obtained in the first order asymptotic formula theorem. Then a modification of hypergamma asymptotics of shifting the argument to the parameter delta is proposed: F(x) = q (Gamma) over tilde1(x+delta)+(1-q)(Gamma) over tilde (0) (x+delta), where delta is found from the condition Gamma(1)(delta). (1-rho)(2). The numerical analysis shows that "modified" distribution significantly improves results of hypergamma asymptotics. In the paper conclusion a combination of "modified" distribution and the second asymptotics pm((2)) 1/2(pm+p((2)) is considered. As a result of numerical analysis, it was observed that this modification is more accurately than the second order asymptotics and can be applied for load rate rho >= 0.7 and even rho >= 0.5 for large values of delay parameter s. In summary, it can be concluded that for practical tasks which does not require high accuracy of approximation, it is useful to apply the "intermediate" asymptotics results because of its simplicity of construction and rather broad area of applicability (rho >= 0.8), and for applications requiring higher accuracy of approximation, it is advisable to apply the proposed modifications 1 and 2. In the future, the modifications can be applied to the study of various retrial queueing systems including with non Poisson incoming process.