## Design, Dimensioning, and Optimization of 4G/5G Wireless Communication Networks

View this Special IssueResearch Article | Open Access

Sławomir Hanczewski, Adam Kaliszan, Maciej Stasiak, "Convolution Model of a Queueing System with the cFIFO Service Discipline", *Mobile Information Systems*, vol. 2016, Article ID 2185714, 15 pages, 2016. https://doi.org/10.1155/2016/2185714

# Convolution Model of a Queueing System with the cFIFO Service Discipline

**Academic Editor:**Ioannis Moscholios

#### Abstract

This article presents an approximate convolution model of a multiservice queueing system with the continuous FIFO (cFIFO) service discipline. The model makes it possible to service calls sequentially with variable bit rate, determined by unoccupied (free) resources of the multiservice server. As compared to the FIFO discipline, the cFIFO queue utilizes the resources of a multiservice server more effectively. The assumption in the model is that the queueing system is offered a mixture of independent multiservice Bernoulli-Poisson-Pascal (BPP) call streams. The article also discusses the results of modelling a number of queueing systems to which different, non-Poissonian, call streams are offered. To verify the accuracy of the model, the results of the analytical calculations are compared with the results of simulation experiments for a number of selected queueing systems. The study has confirmed the accuracy of all adopted theoretical assumptions for the proposed analytical model.

#### 1. Introduction

In recent years there has been a rapid increase in development of networks, in particular of mobile networks. While the evidence indicates that increasing competition is bringing down the cost of mobile services and user equipment, this implicates that more and more of the percentage of total network traffic is generated by mobile devices. According to the report [1], there were 7.3 billion active mobile devices in 2015, while the identified global trends are reporting much increased smartphone-generated traffic with each smartphone generating 1.4 Gb per month on average. It is also worth noticing that the ever-increasing number of users is using wireless broadband access networks. The authors of the report estimate that the number of active devices operating on LTE networks will reach 4.3 billion in 2021.

The dynamic development of telecommunications (mobile) networks, the growing number of offered online services with strictly defined Quality of Service (QoS) parameters, and the ever-increasing number of network users cause network operators to introduce a number of different traffic management mechanisms that increase the effectiveness of the network. Good examples of the above mechanisms are threshold compression mechanisms [2–4] and thresholdless compression mechanism that allow elastic and adaptive traffic to be supported [5, 6]. In the case of systems with real time services it is necessary to ensure an appropriate delay level and/or deviation from true periodicity of a presumed periodic signal (*jitter*). As a result, analyses of these systems are not feasible unless appropriate models of queueing systems that would take into account the multiservice nature of traffic are constructed. A great number of works, including [7–9], discuss a solution that is based on a single-service Erlang C model with the accompanying assumption that all call classes have similar characteristics and service conditions. Another solution proposes an application of a recurrent formula that describes the occupancy distribution in full-availability multiservice systems with elastic traffic [5]. The model presented in [6] has an accompanying assumption that all streams of serviced traffic—if there is an instance of a lack of free resources in the system—can be compressed without limit, which eventually leads to lossless traffic service. Here, the level of compression can be interpreted as a measure of delay for all call classes. These solutions, however, do not allow us to determine delay parameters for individual call classes. In [10–13] new queueing models are proposed that make it possible to determine individual delay characteristics for particular call classes offered to the system with a queueing discipline called state-dependent FIFO (SD-FIFO). This discipline ensures access to a multiservice server to be available for all classes of calls, while the division of the resources of a server between individual call classes depends on the number of all calls that are currently in the queueing system and corresponds to the balanced fairness algorithm [14, 15] for the division of resources in multiservice systems.

This article proposes an approximate model of a queueing system with continuous FIFO (cFIFO) service discipline for a mixture of multiservice Bernoulli-Poisson-Pascal (BPP) call streams. In general, the cFIFO discipline assumes that calls that are in the queue are serviced according to the FIFO discipline. However, when the server has a lower amount of resources than demanded by the first call in the queue, then it can start servicing this call with a lower bitrate than the one demanded by the call. In the proposed model, a convolution algorithm is used to determine the occupancy distribution in the queueing system. Appropriate convolution algorithms to model full-availability multiservice systems with losses and multiservice access network systems (the so-called multiservice tree network) are proposed in [16–18]. Papers [19, 20] propose convolution algorithms to model multiservice full-availability systems with resource reservation. Paper [21] develops and discusses a two-dimensional convolution model for a multiservice overflow system. Yet another work [22] proposes a generalized version of a convolution algorithm that makes it possible to model different multiservice state-dependent systems. The advantage of convolution algorithms is that they provide a possibility to model approximately systems with call streams with different distributions, provided they are not mutually dependent. Hence, in order to present the capabilities of the proposed model, this article also presents the results of a comparison of the analytical model with the data obtained by the digital simulation experiments for traffic streams other than BPP call streams. It is worthwhile to stress that convolution algorithms have never been used in the analysis of queueing systems before.

The article is organized as follows. Section 2 provides a description of a multiservice queueing system with the SD-FIFO discipline. Section 3 presents a description of a multiservice queueing system with the continuous discipline cFIFO. Section 4 proposes an analytical convolution model of the considered queueing system. A number of exemplary results of a comparison of the analytical model with the results of a simulation are provided and discussed in Section 5. The last section, Section 6, is a short summary of the article.

#### 2. Queueing System with the SD-FIFO Discipline

The multiservice queueing system with the SD-FIFO queue service discipline is described in [10, 11, 13]. The analytical model developed for this system allows the most important parameters that characterize the queueing system (such as the average time for a call to be in the queue or the average queue length) to be determined for each call class offered to the system. An exemplary SD-FIFO system is shown in Figure 1. The system is composed of a multiservice server with the capacity allocation units (AUs) (allocation unit is the capacity unit for broadband systems [23, 24]. Typically, it is defined as the greatest common divisor of the maximum bitrates of individual calls [25] or equivalent bandwidths [26] determined for offered call streams. A method for a determination of allocation units in telecommunications systems is provided in the Appendix) and a buffer with the capacity AUs. The notion of the multiservice server is understood to be a server with the capacity that enables concurrent service to a number of calls with different demanded bitrates to be effected. The system (Figure 1) is offered 3 call classes. In the considered SD-FIFO queueing system, virtual queues for each call class serviced by the system are created within one buffer. Depending on the number of calls that are currently in the system, appropriate resources of the server are allocated to individual classes. The amount of these resources varies every time the number of calls of individual classes that are currently in the system is changed (resulting by a termination of service or admittance of a new call). According to this discipline, in each occupancy state of the system, appropriate resources of the service are allocated to calls of all classes. It should be stressed that the method for resource allocation in the server for individual call classes is compatible with the balanced fairness algorithm, well known from the literature of the subject [14, 15].

The system is composed of a multirate server with the capacity AUs and a buffer with the capacity AUs. The system is offered a set of traffic classes of the type (), where(i)Er denotes the Erlang traffic (Poisson call stream),(ii)En denotes the Engset traffic (Bernoulli call stream),(iii)Pa denotes the Pascal traffic (Pascal call stream),(iv) denotes a type of traffic, with non-Poissonian call stream whose properties have been determined either empirically or in a simulation.The cardinality of the set , that is, the number of traffic classes offered to the system, is equal to . Any randomly chosen traffic class demands AUs for service. This notation method will be also used in the remaining part of the article to identify the parameters related to the considered call classes. We notice that the mixture of Erlang, Engset, and Pascal traffics is denoted as BPP traffic.

Let us consider now a queueing system with the SD-FIFO discipline composed of a multiservice server with the capacity AUs and a buffer with the capacity AUs. The system is offered a set of traffic classes of the Erlang-type call classes. Individual call classes are described by the following parameters:(i) is the intensity of a call stream of class ,(ii) is the intensity of a service stream of class ,(iii) is the number of AUs demanded by a call of class ,(iv) is the traffic intensity of class :

The occupancy distribution in the system with the SD-FIFO discipline is determined on the basis of the following recurrence [11]: where is the probability that there are occupied AUs in the system with server capacity AUs and buffer capacity AUs. To the system a set of traffic classes is offered.

It is possible to determine on the basis of the distribution the parameter , that is, the number of calls of class serviced in the server in state : The parameter determines then the resources occupied in the server by a given traffic class and is dependent on the occupancy distribution . This means that the resources occupied in the server by a given call class are dependent on the occupancy state of resources by all traffic classes in the system (serviced in the server and those in the queue).

Thus, Formula (3) defines the queueing service discipline in the system determined by distribution (2), that is, that from a virtual queue of class ; in a given occupancy state of system , this number of calls of this class is taken for service that satisfies (3). Because of these dependencies, this discipline is labeled SD-FIFO [11].

On the basis of distribution (2) one can determine the appropriate queueing characteristics for the system under consideration [11], for example, average waiting time, average queue length, or blocking probability.

#### 3. Queueing System with Continuous cFIFO Discipline

Let us consider the multiservice queueing system presented in Figure 2. Calls that are in the queue are serviced according to the proposed cFIFO discipline. In the case when the multiservice server has a lower throughput than the bitrate required by the first call that is in the queue to be serviced, then the server can start servicing this call with a lower bitrate than the demanded AUs. Figure 2 shows a call, marked with the hyphenated line, that demands 3 AUs. The server only has two free AUs; hence—in line with the cFIFO service discipline—2 AUs of the considered call will be serviced in the server, while one AU will be waiting for service in the buffer. Only one call in the system can be serviced with a lower number of AUs than that demanded for service. Such a concept ensures the maximum usage of the resources of the server providing at the same time a simple algorithm for buffer service.

Let us consider now the operation of a queueing system with the system with the parameters AUs and AUs, presented in Figure 3(a), as an example. This system can be analyzed either at the microstate level (Figure 3(b)) or at the macrostate level (Figure 3(c)). Microstate of the multiservice service process is defined by the number of calls of individual classes serviced in the server and waiting in the buffer [11]. Macrostate is then defined by the total number of AUs that are occupied by calls that are in the system, that is, those that are being serviced in the server and those waiting in the buffer [11]. The macrostate does not take into consideration the distribution of occupied AUs between individual classes of calls. The system presented in Figure 3(a) is offered two call classes that demand and AUs, respectively. Each microstate of the process is represented by the ordered set = , where the values and define the number of calls of classes and that are serviced in the server and the values and define then the number of calls of classes and that are waiting in the buffer. In Figure 3(b), the microstate determines such an occupancy state of the considered queueing system in which a call of class 2 is serviced in the server with a lower number of AUs than that demanded; that is, one AU (“0.5 of the call”) is serviced in the server and one AU (“0.5 of the call”) is waiting in the buffer. Let us consider two possible ways of service termination for a call of class 2 in the microstate under consideration. In the first case, after the termination of service for a call of class 1, the server starts to service a call of class 2 with the demanded number of = 2 AUs (in Figure 3(b) it is the transition from microstate to microstate ). In the second case, because of long service time for the call of class 1, the considered call of class 2 is serviced with AU throughout the service (in Figure 3(b) it is the transition from microstate to microstate ).

**(a)**

**(b)**

**(c)**

#### 4. Model of a Queueing System with cFIFO Discipline

The model of a queueing system proposed in the article has been developed at the macrostate level. This means that the occupancy state is described by just one parameter, namely, the total number of occupied AU in the system (i.e., together in the server and buffer). Each macrostate is the sum of such microstates that satisfies the following condition: A description of a system at the macrostate level greatly simplifies its analysis in that it limits the analysis to a lower number of states. In the system presented in Figure 3(a), and defined by the parameters AUs, AUs, the number of microstates is equal to 12 (Figure 3(b)), whereas the number of macrostates is equal to 5 (Figure 3(c)). In Figure 3(b), the contours of the microstates that belong to one macrostate are marked by identical line. For example, macrostate includes the following microstates: , , and .

##### 4.1. Convolution Algorithm

Multiservice telecommunications systems can be generally modelled by algorithms that analyze either the dependencies between microstates or macrostates. The one approach is characterized by a large computational complexity [27], whereas the other approach is effective in modelling systems that are offered Poisson call streams. Papers [16, 17] propose convolution algorithms that make it possible to model systems with multiservice traffic on the basis of mutually independent occupancy distributions. These distributions are determined for all call streams offered to the system and can be determined on the basis of appropriate models and theoretical single-service systems. They can also directly result from conducted measurements in real systems or from simulation experiments for a number of selected classes of calls.

Let us consider now the operation of a convolution algorithm with a multiservice system that is composed of a server with the capacity and a buffer with zero capacity () as an example [28, 29]. In the literature of the subject this system is frequently labelled as the full-availability group (FAG) [30]. Our assumption is that the system is offered a set of traffic classes of the type () and any randomly chosen traffic class demands AUs for service.

The input data for the convolution algorithm are the occupancy distributions for single classes , where . In the adopted notation for the distribution the first expression in the lower index defines the capacity of the server, while the second element defines the capacity of the buffer (in the system under investigation, the buffer with zero capacity is considered). The distributions are determined with the assumption that the system with the capacity is offered only one class of calls. For Erlang, Engset, or Pascal traffic classes, to determine this distribution one can use appropriate theoretical distributions [31]; for example, in the case of class of Erlang traffic, the distribution for a single class can be determined on the basis of the Erlang distribution:where and is the traffic intensity of class , whereas the parameter denotes the maximum number of calls that can be serviced in the server with the capacity :

In the case of Engset and Pascal traffic, the distributions of single classes can be determined on the basis of the following formulas:where the parameters and define the average intensity of calls offered by a single free Engset and Pascal traffic source. The parameters and are the total number of Engset and Pascal traffic sources. The relationship between the average traffic offered to the system by Engset traffic stream and Pascal traffic stream and the parameters and is as follows:

The occupancy distribution , , of a single class , determining the occupancy probabilities for all macrostates (such that AUs) in the server with the capacity , is related to distribution (5) by the following dependence:

In the case of other traffic streams (of the type ), distributions of single classes can be determined on the basis of measurements taken in real systems or on the basis of relevant simulation experiments.

By having the single distributions for all classes offered to the system at our disposal we can now determine the aggregated distribution which is the result of the convolution operation of all distributions for single classes offered to the system: The distribution is not a normalized distribution. After an appropriate normalization process, the normalized distribution can be written in the following form: where is the normalization constant:

Let us note that, for example, as a result of the convolution operation of two normalized distributions with the length , a distribution with the length is obtained. From the point of view of the considered system with the capacity of AUs, all states will never occur. Therefore, such states, in the distribution with the length must be removed and distribution must be normalized. This process is discussed in detail in [16, 17, 19], among others.

In (12), the symbol denotes the convolution operation which is defined in the following way: where the symbol denotes the normalized distribution of a single class or the aggregated distribution of a number of classes that belong to a given set (), whereas the symbol denotes the normalized distribution of a single class or the aggregated normalized distribution of a number of classes that belong to a given set ( and ). The symbol defines the aggregated nonnormalized distribution that is the result of a convolution of the distributions from the sets and .

The convolution operation (15) enables us to determine the average number of serviced calls of class in the server with the capacity AUs with zero buffer that is in the occupancy state AUs. A determination of this parameter should take into consideration the following reasoning [30]. First, the aggregated nonnormalized distributions of all classes, except class , are to be determined:

Then, on the basis of the distributions and the value of the parameter can be determined, that is, the average number of calls of class in occupancy state of the server:

##### 4.2. Model of the cFIFO Queueing System for Poisson Traffic Streams

Let us consider now the queueing system with the cFIFO service discipline the operation of which is presented in Section 3. Our assumption is that this system is offered Erlang traffic streams for which call arrival intensities of new calls are described by Poisson distributions with the parameters for each traffic class, respectively. Service time is described with exponential distributions with the parameters . Individual call classes demand, respectively, AUs for service. In the case of Erlang traffic, it is assumed that the number of traffic sources is infinite. This means that the call arrival intensity for new call arrivals of a given class is independent of the number of calls of this class that are currently in the system: The intensity of offered Erlang traffic of class can be determined on the basis of Formula (1):By expressing the call intensity for a given class in allocation units per time unit (the product ), we can express traffic offered to the system in the following form:

Paper [30] analyzes the dependence between the occupancy distribution in a multiservice full-availability system and the occupancy distribution in a multiservice state-dependent system, with the assumption of identical capacity and traffic offered to both systems. A similar approach will be now applied to the queueing model proposed in the article to determine the dependence between the occupancy distribution in the server with zero buffer and the distribution in the considered queueing system: where defines the relation between the corresponding probabilities and , that is, the probability that the queueing system is in state AUs and the occupancy probability AUs in the server with zero buffer , with the assumption that traffic offered to both systems is identical. To determine the parameter we can use the recurrence dependencies derived for multiservice SD-FIFO queueing systems. The distribution in the cFIFO system is then approximated by distribution (2) that, with (20) taken into consideration, will be rewritten in the following way:

Assuming, then, that the method for a determination of the occupancy distribution in the queueing system and in the system with zero buffer (when the system is offered Erlang call classes) is known, it is possible to determine the value of the parameter by dividing both sides of (22) by the probabilities : Equation (23), with (21) taken into consideration, can be easily transformed into an iterative formula for the parameter :

The Markov process in the server with zero buffer to which a mixture of Erlang traffic classes is offered is a reversible process [28, 29]. For the neighboring states and , because of class , it is possible then to write local balance equations: Now, taking into consideration (25), (24) can be rewritten in the following form:

The convolution operation makes a determination of the occupancy distribution and the average number of calls of individual classes serviced in the server with zero buffer (Formulas (16) and (17)) possible. These results enable us then to assess the value of the parameter (Formula (26)) and eventually to evaluate the occupancy distribution in the queueing system on the basis of (21).

##### 4.3. Commentary

In Section 4.2, to find the general dependence between the cFIFO queueing system described by the parameters () and the system with zero buffer () the occupancy distribution in the SD-FIFO system is used. The possibility of an approximation of the cFIFO system by the SD-FIFO system results from significant similarities in both queue service disciplines involved. In the cFIFO discipline considered in the article only one call can be partly serviced in the server and partly waiting in the queue. The SD-FIFO discipline, in turn, allows even a number of calls to be partly serviced simultaneously. In both instances, however, resources of the server are used maximally and it is just this fact that determines the probability of occupancy distributions for systems with cFIFO and SD-FIFO queues. Figures 4 and 5 show the occupancy distribution in a queueing system with the structural parameters AUs and AUs, respectively. The occupancy distributions in the SD-FIFO system and the cFIFO system obtained in a simulation were then compared. The systems were offered three Erlang traffic streams with the following demands: AU, AUs, and AUs. Traffic was offered in the following proportions: . Figure 4 shows the occupancy distributions for the case of small system loads ( Erl/AU) and Figure 5 for large system loads ( Erl/AU), where is the average traffic offered per AU of the system:

The simulation results are shown with 95% confidence interval determined on the basis of the Student distribution for 5 series, 100,000 calls each. The results for the presented comparison indicate a very good convergence of occupancy distributions in the SD-FIFO and cFIFO queueing systems. The simulation experiments for other systems differentiated by their capacity, number, and demands of offered traffic classes conducted by the authors earlier confirm strong convergence of occupancy distributions of both queueing systems.

It should be stressed that queueing characteristics (such as the average number of calls of particular classes in the queue in corresponding occupancy states of the system) are not characterized by just as good convergence as the occupancy distributions. This would be the first reason why the SD-FIFO model cannot be directly applied to determine queueing characteristics for the cFIFO system. The other reason is the fact that to determine the parameter (Formula (26)) the occupancy distribution of the SD-FIFO queueing system, derived for Erlang traffic, is used. In the case of other traffic streams, we do not observe such a good convergence between distributions for the SD-FIFO and cFIFO systems any longer. However, Formula (26) turns out to be very universal and versatile and approximates cFIFO queueing systems that service traffic with any, mutually independent, call streams very well. Appropriate numerical examples are presented in Section 5.

##### 4.4. Characteristics of cFIFO System with Erlang Traffic

By having the knowledge of the coefficients at hand we are in position to determine the occupancy distribution in a queueing system with the cFIFO discipline. This distribution can also provide a basis for a determination of appropriate QoS parameters, such as the blocking probability and the average number of occupied AUs in queue . The phenomenon of blocking for calls of class in the queueing system occurs when the queue lacks free AUs to store a new call of class . Therefore, The average queue length (for calls of all classes) is a mean value for the resources (expressed in AUs) occupied by calls of all classes that are waiting in the queue:

Let us consider now a method for a determination of the average number of busy AUs in a queue occupied by calls of class . This article proposes an approximate determination of this parameter on the basis of the following reasoning. First, it is the average number of AUs occupied by calls of class that are waiting in the queue that has to be determined: where is the average number of calls of class waiting in the queue in state of the system. The product is the average number of AUs occupied in the buffer by waiting calls of class in state . In states such that , the queue is empty; therefore the product is equal to zero.

In states such that , the total number of occupied AUs in the queue (by calls of all classes) can be determined as follows: Between the neighboring states and there is a difference in the number of busy AUs occupied by calls of class waiting in the queue, which, in the article, is denoted by the symbol . We can thus write where Let us assume that the parameter is proportional to the intensity in which the buffer is occupied by calls of a given class . We can therefore determine the value of the parameter in the following way: where Equation (35) is the result of a particular service discipline in the cFIFO queue. If in system state a new call of class arrives, then in the occupancy states of the server that are determined by the condition a number of AUs of the considered call of class will be immediately serviced by the server, while the remaining part will be directed to the queue. If , then all AUs of a new call of class will be immediately placed in the queue. Equations (34) and (35) can be written as just one equation in the following way: Ultimately, (32), taking into account (36), can be written as follows: where . A determination of all values will make it possible, on the basis of (30), to evaluate and assess all queue lengths for calls of individual classes.

##### 4.5. Model of a Queueing System with cFIFO Discipline and State-Dependent Call Streams

To determine the occupancy distribution in a queueing system with the cFIFO queueing discipline and offered BPP traffic (or traffic with any distribution of the call stream) we will use the general recurrence dependencies derived for a system with Erlang traffic (Section 4.4). To determine the occupancy distribution in the cFIFO queueing system with BPP traffic we use then dependencies (21) and (25) that, in their generalized form, can be rewritten as follows: where the parameter is determined on the basis of a convolution operation (Formulas (16) and (17)). The occupancy distribution for the server with zero buffer in Formula (38) is determined on the basis of a convolution algorithm for which the input data are the distributions of single classes, determined for Erlang traffic by Formula (5) and for Engset and Pascal traffic by Formulas (7) and (9) as well as (8) and (10), respectively.

It should be emphasized that the application of the convolution algorithm in the proposed model makes it possible to determine characteristics of queueing systems for call streams other than BPP streams. Distributions of single classes that are the input data for the convolution algorithm can be then determined empirically on the basis of measurements or by simulation experiments.

On the basis of the occupancy distribution , according to (28) and (29), it is possible to determine the blocking probability and the average queue length for calls of all classes.

Let us consider now the method for a determination of the average queue length for calls of a single class . We will first determine the average number of occupied resources in the buffer in state of the system by calls of class . By the adoption of the same initial assumptions (30)–(33) that were earlier adopted for the Erlang call stream, (32) can be rewritten in such a way as to include any type traffic : Equation (36) will be also rewritten in our adopted notation to include the dependence between the call stream and the state: where determines the call arrival intensity of new calls of any traffic class of type in state of the system.

In the case of Engset traffic, the parameter in state depends on the total number of traffic sources and on the number of currently serviced calls in the server and the number of calls waiting in the buffer : where is the intensity of the call arrival process for calls of class , determined with the assumption that all sources are not occupied: By approximating the number of occupied resources in the considered queueing system with nonzero buffer by calls of class in state by the average number of occupied resources in state in the server with zero buffer, Formula (42) can be rewritten as follows: Using the analogous simplifying assumptions as for Engset traffic we are in position to determine the intensity of the call arrival process for Pascal calls in state of the system in the following way: The call arrival intensities that have been described by (44) for calls of Engset class, and (45) for Pascal class, make it possible to determine, on the basis of (40) and (41), the average number of occupied resources of a queue by calls of individual classes.

Independently of the considered call streams, the proposed model can be written in the form of the following method, henceforth called the cFIFO method.

*cFIFO Method*(1)Determination of occupancy distributions for individual classes: in multiservice server with the capacity and with zero buffer.(2)Determination of nonnormalized aggregated occupancy distributions of all classes, except a class (): , .(3)Determination of the average number of serviced calls in state () for each class ().(4)Determination, by convolution algorithm, of the occupancy distribution for multiservice server with the capacity and with zero buffer.(5)Determination of transformation coefficients .(6)Determination of the occupancy distribution for queueing system with the cFIFO discipline.(7)Determination of the average queue length and the blocking probability () for calls of individual classes.(8)Determination of the average number of occupied AUs by calls of class () waiting in the queue in such states of the system that .(9)Determination of the average number of occupied AUs by calls of class , (), waiting in the queue.

#### 5. Numerical Results

In order to verify the proposed analytical method the results of the calculations were compared with the results provided by the simulation experiments. For this purpose, a dedicated simulation program for a cFIFO queueing system evaluation was developed. The simulator uses an experiment with steady system time in which the process interaction method is used. The classes with finite number of sources were implemented according to the principles described in [32]. The program was written in the C++ language in the Qt environment. The software architecture and process interaction method is described in [33]. The obtained results are presented in graphs as the function of traffic offered to one AU in the server:

Each simulation experiment to determine the characteristics of the system under investigation for particular values involved 10 series. The prerequisite condition for a single series to be completed was call service for 10,000,000 calls or a loss of 10,000 calls. The adopted length of a single series in the simulation made it possible to determine confidence intervals at the 95% level. The results obtained on the basis of the analytical model are shown in the graphs with lines, whereas the results of the simulation are presented with appropriate symbols. In most cases, the values for the confidence intervals in the graphs are so small that they do not exceed the value of the symbol that denotes the result of a simulation. The proportions of offered traffic in all presented cases were . For each of the considered queueing systems, the results for the blocking probability and the average queue lengths for individual classes, expressed in the number of occupied AUs, will be presented.

Figures 6–8 show the results for three queueing systems that service Erlang traffic. The assumption in the systems presented in Figures 6 and 7 was that the intensities of service streams of particular call classes were, respectively, equal to , , , whereas in the case of the results shown in Figure 8 they were, respectively, , , .

The next stage in the process of verification of the proposed analytical model was to test its accuracy for systems to which a mixture of Erlang, Engset, and Pascal traffic was offered (Figure 9). In this case, the intensities of service streams of particular call classes were, respectively, equal to , , . The number of traffic sources for Engset and Pascal classes is equal to , , respectively.

In the first case of considered types of non-Poissonian offered traffic (denoted by symbol NE), each call stream is described by normal distribution and each service stream is described by exponential distribution. It is assumed that the mean value is equal to the inverse of offered traffic intensity of a given call class () and variance is equal to . So, we consider the system with normally distributed call stream . The intensities of service streams are equal to (similar to previously considered queueing systems). The results for this queueing system are presented on Figure 10.

Figure 11 shows the results for a queueing system that, for each serviced traffic class, the call stream and service stream are described by normal distribution. This type of traffic is denoted by symbol NN. Similar to previous case, each call stream is described by the distribution . Each service stream is given by the distribution , where in the experiment was admitted .

The results presented on Figure 12 were worked out for queueing system in which call streams are described by uniform distribution and service streams are described by exponential distribution. This type of traffic is denoted by symbol UE. It is assumed that for each traffic class the uniform distribution is given by parameters and (where is the traffic intensity of class ). The intensities of service streams are equal to , , .

The results of the study showed in Figures 6–12 confirm high accuracy of the proposed approximate analytical model of the queueing system with the cFIFO discipline. The diversified accuracy of the results stems from the dependence between the method for a determination of the average number of calls for individual classes in the system adopted in the study (independently for the server and queue) and the type of offered traffic. More accurate results were obtained for BPP streams, while those for non-Poisson streams were less accurate. This results from the adoption of the method of the approximation of the distribution in the system with the cFIFO discipline by the distribution in the SD-FIFO system that can be only precisely determined for Poisson streams. During the study no dependence between the accuracy of the model and the number of offered classes of calls and their demands was observed.

#### 6. Conclusions

This article proposes an approximate analytical model of a queueing system that services multiservice traffic with the cFIFO queueing discipline. The model is based on a convolution algorithm in which to determine certain characteristics, such as the average number of serviced calls, the formulas derived on the basis of accurate Markovian models for multiservice systems with zero and nonzero buffer (with the SD-FIFO service discipline for the queue) are used. The introduction of these assumptions makes it possible to construct a general model of a multiservice queueing system with the cFIFO discipline. The proposed model is characterized by high accuracy, which makes an analysis of cFIFO systems for call streams with any, mutually independent, probabilistic distributions possible. This high accuracy stems from the fact that any errors that may possibly result from the adopted approach will be minimized by the convolution operation.

#### Appendix

#### Resource Allocation Units

The application of multirate analytical models to analyze present-day broadband telecommunications systems is possible with a proper bandwidth discretization. Discretization allows the capacity of a system and bit rates demanded by calls to be expressed in the so-called allocation unit (AU).

The discretization process itself consists in replacing a variable bit rate (VBR) of call streams with a specific constant bit rate (CBR) called the equivalent bandwidth. There are many methods and algorithms specific methods of determining the equivalent bandwidth for given network types and services in the literature of the subject, for example, [24, 26, 34–37]. The value of equivalent bandwidth depends on various parameters, including the admissible packet-loss ratio, admissible delay, link flow capacity (bit rate), average and peak bit rates, packet stream type (e.g., self-similar streams), and network type.

The procedure of discretization process for system services call classes any type of , and can be written in the form of the following algorithm.

*Algorithm 1. * *Discretization*(1)Determination of the equivalent bandwidth for all offered call classes: .(2)Determination of the allocation unit (AU). The value of this parameter can be calculated as the greatest common divisor (GCD) of all equivalent bandwidths in a considered system: (3)Determination of the number of allocation units demanded by calls of class : (4)Determination of the capacity of the system (expressed in allocation units): where is the capacity of the telecommunications system under consideration.

#### Competing Interests

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

#### Acknowledgments

The present work is financed with the Ministry of Science and Higher Education resources for academic purposes in the year 2016 within the frame of own research project entitled “The Structure, Analysis and Designing of Modern Switching Systems and Telecommunications Networks.”

#### References

- Ericsson, “Ericsson mobility report,” Tech. Rep., Ericsson, 2016. View at: Google Scholar
- M. Sobieraj, M. Stasiak, J. Weissenberg, and P. Zwierzykowski, “Analytical model of the single threshold mechanism with hysteresis for multi-service networks,”
*IEICE Transactions on Communications*, vol. 95, no. 1, pp. 120–132, 2012. View at: Publisher Site | Google Scholar - J. S. Kaufman, “Blocking with retrials in a completely shared resource environment,”
*Performance Evaluation*, vol. 15, no. 2, pp. 99–116, 1992. View at: Publisher Site | Google Scholar - I. D. Moscholios, M. D. Logothetis, J. S. Vardakas, and A. C. Boucouvalas, “Performance metrics of a multirate resource sharing teletraffic model with finite sources under the threshold and bandwidth reservation policies,”
*IET Networks*, vol. 4, no. 3, pp. 195–208, 2015. View at: Publisher Site | Google Scholar - G. M. Stamatelos and V. N. Koukoulidis, “Reservation-based bandwidth allocation in a radio ATM network,”
*IEEE/ACM Transactions on Networking*, vol. 5, no. 3, pp. 420–428, 1997. View at: Publisher Site | Google Scholar - T. Bonald and J. Virtamo, “A recursive formula for multirate systems with elastic traffic,”
*IEEE Communications Letters*, vol. 9, no. 8, pp. 753–755, 2005. View at: Publisher Site | Google Scholar - J. W. Cohen, “The multiple phase service network with generalized processor sharing,”
*Acta Informatica*, vol. 12, no. 3, pp. 245–284, 1979. View at: Publisher Site | Google Scholar | MathSciNet - K. Lindberger, “Balancing quality of service, pricing and utilisation in multiservice networks with stream and elastic traffic,” in
*Proceedings of the 16th International Teletraffic Congress*, pp. 1127–1136, Edinburgh, UK, 1999. View at: Google Scholar - T. Bonald and J. W. Roberts, “Internet and the Erlang formula,”
*ACM SIGCOMM Computer Communication Review*, vol. 42, no. 1, pp. 23–30, 2012. View at: Publisher Site | Google Scholar - S. Hanczewski, M. Stasiak, and J. Weissenberg, “The queueing model of a multiservice system with dynamic resource sharing for each class of calls,” in
*Computer Networks*, A. Kwiecień, P. Gaj, and P. Stera, Eds., vol. 370 of*Communications in Computer and Information Science*, pp. 436–445, Springer, Berlin, Germany, 2013. View at: Publisher Site | Google Scholar - S. Hanczewski, M. Stasiak, and J. Weissenberg, “A queueing model of a multi-service system with state-dependent distribution of resources for each class of calls,”
*IEICE Transactions on Communications*, vol. 97, no. 8, pp. 1592–1605, 2014. View at: Google Scholar - S. Hanczewski, D. Kmiecik, M. Stasiak, and J. Weissenberg, “Multiservice queuing system with elastic traffic,” in
*Proceedings of the IEICE General Conference*, pp. S-46–S-47, The Institute of Electronics, Information and Communication Engineers, March 2016. View at: Google Scholar - M. Stasiak, “Queuing systems for the internet,”
*IEICE Transactions on Communications*, vol. 99, no. 6, pp. 1234–1242, 2016. View at: Publisher Site | Google Scholar - T. Bonald, A. Proutière, and J. Virtamo, “A queueing analysis of max-min fairness, proportional fairness and balanced fairness,”
*Queueing Systems*, vol. 53, no. 1, pp. 65–84, 2006. View at: Publisher Site | Google Scholar | MathSciNet - J.-P. Haddad and R. R. Mazumdar, “Congestion in large balanced multirate networks,”
*Queueing Systems*, vol. 74, no. 2-3, pp. 333–368, 2013. View at: Publisher Site | Google Scholar | MathSciNet - V. Iversen, “The exact evaluation of multi-service loss system with access control,” in
*Proceedings of the 17th Nordic Teletraffic Seminar*, pp. 56–61, Lund, Sweden, August 1987. View at: Google Scholar - K. Ross,
*Multiservice Loss Models for Broadband Telecommunication Network*, Springer, London, UK, 1995. - I. D. Moscholios, M. D. Logothetis, J. S. Vardakas, and A. C. Boucouvalas, “Congestion probabilities of elastic and adaptive calls in Erlang-Engset multirate loss models under the threshold and bandwidth reservation policies,”
*Computer Networks*, vol. 92, part 1, pp. 1–23, 2015. View at: Publisher Site | Google Scholar - M. Glabowski, A. Kaliszan, and M. Stasiak, “Asymmetric convolution algorithm for full-availability group with bandwidth reservation,” in
*Proceedings of the Asia-Pacific Conference on Communications (APCC '06)*, IEEE, Busan, South Korea, September 2006. View at: Publisher Site | Google Scholar - M. Głąbowski, A. Kaliszan, and M. Stasiak, “On the application of the asymmetric convolution algorithm in modeling of full-availability group with bandwidth reservation,” in
*Managing Traffic Performance in Converged Networks: 20th International Teletraffic Congress, ITC20 2007, Ottawa, Canada, June 17–21, 2007. Proceedings*, L. Mason, T. Drwiega, and J. Yan, Eds., vol. 4516 of*Lecture Notes in Computer Science*, pp. 878–889, Springer, Berlin, Germany, 2007. View at: Publisher Site | Google Scholar - M. Głąbowski, A. Kaliszan, and M. Stasiak, “Two-dimensional convolution algorithm for modelling multiservice networks with overflow traffic,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 852082, 18 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - A. Kaliszan, M. Głabowski, and M. Stasiak, “Generalised convolution algorithm for modelling state-dependent systems,”
*IET Circuits, Devices & Systems*, vol. 8, no. 5, pp. 378–386, 2014. View at: Publisher Site | Google Scholar - J. Roberts, Ed.,
*Performance Evaluation and Design of Multiservice Networks*, Final Report COST 224, Commission of the European Communities, Brussels, Belgium, 1992. - J. Roberts, V. Mocci, and I. Virtamo, Eds.,
*Broadband Network Teletraffic, Final Report of Action COST 242*, Commission of the European Communities, Springer, Berlin, Germany, 1996. - J. Y. Hui, “Resource allocation for broadband networks,”
*IEEE Journal on Selected Areas in Communications*, vol. 6, no. 9, pp. 1598–1608, 1988. View at: Publisher Site | Google Scholar - F. Kelly, “Notes on effective bandwidth,” Tech. Rep., University of Cambridge, 1996. View at: Google Scholar
- J. Karlsson, “Loss performance in trunk groups with different capacity demands,” in
*Proceedings of the 13th International Teletraffic Congress*, vol. Discussion, pp. 201–212, Copenhagen, Denmark, 1991. View at: Google Scholar - J. S. Kaufman, “Blocking in a shared resource environment,”
*IEEE Transactions on Communications*, vol. 29, no. 10, pp. 1474–1481, 1981. View at: Publisher Site | Google Scholar - J. Roberts, “A service system with heterogeneous user requirements—application to multi-service telecommunications systems,” in
*Proceedings of Performance of Data Communications Systems and Their Applications*, G. Pujolle, Ed., pp. 423–431, North Holland, Amsterdam, The Netherlands, 1981. View at: Google Scholar - M. Głąbowski, A. Kaliszan, and M. Stasiak, “Modeling product-form state-dependent systems with BPP traffic,”
*Performance Evaluation*, vol. 67, no. 3, pp. 174–197, 2010. View at: Publisher Site | Google Scholar - R. Syski,
*Introduction to Congestion Theory in Telephone Systems*, Studies in Telecommunication, North Holland, The Netherlands, 1986. - M. Głabowski and A. Kaliszan, “Simulator of full-availability group with bandwidth reservation and multi-rate bernoulli-poisson-pascal traffic streams,” in
*Proceedings of the International Conference on “Computer as a Tool” (EUROCON '07)*, pp. 2271–2277, Warszawa, Poland, September 2007. View at: Publisher Site | Google Scholar - S. Hanczewski and A. Kaliszan, “Simulation studies of queueing systems,” in
*Proceedings of the 10th International Symposium on Communication Systems, Networks and Digital Signal Processing (CSNDSP '16)*, Prague, Czech Republic, July 2016. View at: Google Scholar - T. Neame, M. Zuckerman, and R. Addie, “A paractical approach for multimedia traffic modeling,” in
*Proceedings of the 5th International Conference on Broadband Communications*, pp. 73–82, Hong Kong, November 1999. View at: Google Scholar - N. L. S. Fonseca, G. S. Mayor, and C. A. V. Neto, “On the equivalent bandwidth of self-similar sources,”
*ACM Transactions on Modeling and Computer Simulation*, vol. 10, no. 2, pp. 104–124, 2000. View at: Publisher Site | Google Scholar - S. Bodamer and J. Charzinski, “Evaluation of effective bandwidth schemes for self-similar traffic,” in
*Proceedings of the 13th ITC Specialist Seminar on IP Traffic Measurement, Modeling and Management*, pp. 21.1–21.10, Monterey, Calif, USA, 2000. View at: Google Scholar - J. Charzinski, “Internet traffic measurment and characterisation results,” in
*Proceedings of the 13th International Symposium on Services and Local Access (ISSLS '00)*, Stockholm, Sweden, June 2000. View at: Google Scholar

#### Copyright

Copyright © 2016 Sławomir Hanczewski et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.