## 1. Introduction

Channel holding time (CHT) is of paramount importance for the analysis and performance evaluation of mobile cellular networks. This time variable allows one to derive other key system parameters such as channel occupancy time, new call blocking probability, and handoff call dropping probability. CHT depends on cellular shape, cell size, user’s mobility patterns, used handoff scheme, and traffic flow characteristics. Traffic flow characteristics are associated with unencumbered service time (UST), while the overall effects of cellular shape, users’ mobility, and handoff scheme are related to cell dwell time (CDT).

For convenience and analytical/computational tractability, the teletraffic analysis of mobile cellular networks has been commonly performed under the unrealistic assumption that CDT and/or CHT follow the negative exponential distribution (Lin et al., 1994; Hong & Rappaport, 1986). However, a plenty of evidences showed that these assumptions are not longer valid (Wang & Fan, 2007; Christensen et al., 2004, Fang, 2001, 2005; Orlik & Rappaport, 1998; Fang & Chlamtac, 1999; Fang et al., 1999; Alfa & Li, 2002; Rahman & Alfa, 2009; Soong & Barria, 2000; Yeo & Jun, 2002; Pattaramalai, et al., 2007). Recent papers have concluded that in order to capture the overall effects of users’ mobility, one needs suitable models for CDT distribution (Lin, 1994; Hong & Rappaport, 1986). In specific, the use of general distributions for modeling this time variable has been highlighted. In this research direction, some authors have used Erlang, gamma, uniform, deterministic, hyper-Erlang, sum of hyper-exponentials, log-normal, Pareto, and Weibull distributions to model the pdf of CDT; see (Wang & Fan, 2007; Fang, 2001, 2005; Orlik & Rappaport, 1998; Fang & Chlamtac, 1999; Fang et al., 1997, 1999; Rahman & Alfa, 2009; Pattaramalai et al., 2007, 2009; Hidata et al., 2002; Thajchayapong & Toguz, 2005; Khan & Zeghlache, 1997; Zeng et al. 2002; Kim & Choi, 2009) and the references therein. Fang in (Fang, 2001)) emphasizes the use of phase-type (PH) distributions for modeling CDT. The reason is twofold. First, PH distributions provide accurate description of the distributions of different time variables in wireless cellular networks, while retaining the underlying Markovian properties of the distribution. Markovian properties are essential in generating tractable queuing models for cellular networks. Second, there have been major advances in fitting PH distributions to real data. Among the PH probability distributions, the use of either Coxian or Hyper-Erlang distributions are of particular interest because their universality property (i.e, they can be used to approximate any non negative distribution arbitrarily close) (Soong & Barria, 2000; Fang, 2001).

Due to the discrepancy and the wide variety of proposed models, it appears mandatory to investigate the implications of the cell dwell time distribution on channel holding time characteristicsin mobile wireless networks. This is the topic of research of the present chapter. Let us describe the related work reported in this research direction.

### 1.1. Previously related work

In (Fang, 2001; Zeng et al. 2002), it is observed that, depending on the variance of CDT, the mean channel holding time for new calls (CHTn) can be greater than the mean channel holding time for handoff calls (CHTh). However, in these works, it is neither explained nor discussed the physical reasons for this observed behavior. This phenomenon (which is addressed in Section 3.1) and the lack of related published numerical results have motivated the present chapter.

Most of the previously published papers that have developed mathematical models for the performance analysis of mobile cellular systems considering general probability distribution for cell dwell time have either only presented numerical results for the Erlang (Wang & Fan, 2007; Fang et al., 1999; Rahman & Alfa, 2009; Kim & Choi, 2009) or Gamma distributions with shape parameter greater than one[1] - (Yeo & Jun, 2002; Fang, 2005), and/or only for the CHTh[2] - (Fang, 2001; Fang & Chlamtac, 1999), or have not presented numerical results at all (Fang, 2005; Alfa & Li, 2002; Soong & Barria, 2000). Thus, numerical results both for values of the coefficient of variation (CoV) of CDT greater than one and/or for the CHTn have been largely ignored. Exceptions of this are the papers (Orlik & Rappaport, 1998; Fang et al., b, 1997; Pattaramalai, et al., 2009).

On the other hand, probability distribution of CHT has been determined under the assumption of the staged distributions sum of hyper-exponentials, Erlang, and hyper-Erlang for the CDT (Orlik & Rappaport, 1998; Soong & Barria, 2000). However, to the best of the authors’ knowledge, probability distribution of CHT in mobile cellular networks with neither hyper-exponential nor Coxian distributed CDT has been previously reported in the literature.

In this Chapter, the statistical relationships among residual cell dwell time (CDTr), CDT, and CHT for new and handoff calls are revisited and discussed. In particular, under the assumption that UST is exponentially distributed and CDT is phase-type distributed, a novel algebraic set of general equations that examine the relationships both between CDT and CDTr and between CDT and channel holding times are obtained. Also, the condition upon which the mean CHTn is greater than the mean CHTh is derived. Additionally, novel mathematical expressions for determining the parameters of the resulting CHT distribution as functions of the parameters of the CDT distribution are derived for hyper-exponentially or Coxian distributed CDT.

## 2. System model

A homogeneous multi-cellular system with omni-directional antennas located at the centre of each cell is assumed; that is, the underlying processes and parameters for all cells within the cellular network are the same, so that all cells are statistically identical.As mobile user moves through the coverage area of a cellular network, several variables can be defined: cell dwell time, residual cell dwell time, channel holding time, among others. These time variables are defined in the next section.

### 2.1. Definition of time interval variables

In this section the different time interval variables involved in the analytical model of a mobile cellular network are defined.

First, the *unencumbered service time* per call *x*_{s} (also known as the *requested call holding time* (Alfa and Li, 2002) or *call holding duration* (Rahman & Alfa, 2009)) is the amount of time that the call would remain in progress if it experiences no forced termination. It has been widely accepted in the literature that the unencumbered service time can adequately be modeled by a negative exponentially distributed random variable (RV) (Lin et al., 1994; Hong & Rappaport, 1986). The RV used to represent this time is **X**_{s} and its mean value is

Now, c*ell dwell time* or *cellresidence timex*_{d}^{(j)} is defined as the time interval that a mobile station (MS) spends in the *j-*th (for *j* = 0, 1, …) handed off cell irrespective of whether it is engaged in a call (or session) or not. The random variables (RVs) used to represent this time are **X**_{d}^{(j)} (for *j* = 0, 1, …) and are assumed to be independent and identically generally phase-type distributed. For homogeneous cellular systems, this assumption has been widely accepted in the literature (Lin et al., 1994; Hong & Rappaport, 1986; Orlik & Rappaport, 1998; Fang & Chlamtac, 1999; Alfa & Li, 2002; Rahman & Alfa, 2009).

In this Chapter, cell dwell time is modeled as a general phase-type distributed RVwith the probability distribution function (pdf)

The *residual cell dwell timex*_{r} is defined as the time interval between the instant that a new call is initiated and the instant that the user is handed off to another cell. Notice that residual cell dwell time is only defined for new calls. The RV used to represent this time is **X**_{r}. Thus, the probability density function (pdf) of **X**_{r}, **X**_{d} using the excess life theorem (Lin et al., 1994)

where *E*[**X**_{d}] and _{d}.

Finally, we define *channel holding time* as the amount of time that a call holds a channel in a particular cell.In this Chapter we distinguish between channel holding times for handed off (CHTh) and channel holding time for new calls (CHTn). CHTh (CHTn) is represented by the random variable

## 3. Mathematical analysis

### 3.1. Relationship between X_{d} and X_{r}

The relationship between the probability distributions of CDT and CDTr is determined by the residual life theorem. In Table I some particular typically considered CDT distributions and the corresponding CDTr distributions obtained by applying the residual life theorem are shown.

The functional relationship between the moments of the residual cell dwell time and the cell residual time was obtained in (Kleinrock, 1975) applying the Laplace transform to the residual life theorem. That is,

This equation can be rewritten as

The *n-th* moment of the residual cell dwell time in terms of the moments of the cell dwell time can be obtained by deriving *n* times (equation 3) with negative argument and substituting *s*=0. Then (Kleinrock, 1975),

The mean residual cell dwell time as function of the moments of cell dwell time can be obtained as (Kleinrock, 1975)

*CoV*{**X**_{d}} represents the coefficient of variation of CDT, the condition for which the mean CDTr is greater than the mean CDT

In this way, the relationship between mean CDT and mean CDTr only depends on the value of the *CoV* of CDT. Thus, the mean CDTr is greater than the mean CDT (i.e., *E*{**X**_{r}} > *E*{**X**_{d}}) when the *CoV* of CDT is greater than one. This behavior (i.e., *E*{**X**_{r}} > *E*{**X**_{d}}) may seem to be counterintuitive due to the fact that, for a particular realization and by definition, CDTr cannot be greater than CDT[3] -
. This occurs because in such conditions there is a high variability on the cell dwell times in different cells and it is more probable to start new calls on cells where users spent more time. Then, residual cell dwell times tend to be greater than the mean CDT. This phenomenon that may seem to be counterintuitive is now explained and mathematically formulated in this Chapter.

### 3.2. Channel holding time distribution for handed off and new calls

Channel holding times for handed off and new calls (denoted by X_{C}^{(h)} and X_{C}^{(N)}, respectively) are given by the minimum between UST and CDT or CDTr, respectively. The CDF of the CHTh and CHTn are, respectively, given by

Due to the fact that the Laplace transform of the pdf of both UST and CDTr are rational functions, the Laplace transform of the pdf of CHTn can be obtained using the Residue Theorem as follows (Wang & Fan, 2007)

where **X**_{r}) by cell dwell time (**X**_{d}).

Under the condition that UST is general phase type (PH) distributed, the authors of (Alfa & Li, 2002) prove that the CDT is PH distributed if and only if the CHTn is PH distributed or the CHTh is PH distributed.

The probability distributions of CHTn and CHTh for different staged probability distributions of CDT assuming that the UST is exponentially distributed are shown in Table II. The first entry of this table is a well known result[4] - . In (Soong & Barria, 2000), it was shown that when CDT has Erlang or hyper-Erlang distribution, channel holding times have the uniform Coxian and hyper-uniform Coxian distribution, respectively.Uniform Coxian is a special case of the Coxian distribution where all the phases have the same parameter (Perros & Khoshgoftaar, 1989). The hyper-uniform Coxian distribution is a mixture of uniform Coxian distributions.

Next, it is shown that when the UST is exponentially distributed and CDT has hyper-exponential distribution of order *n*, the distribution of CHTh has also a hyper-exponential distribution of order *n*. Similarly, when CDT has Coxian distribution of order n, the distribution of CHTn has also a Coxian distribution of order *n*.

#### 3.2.1. Case 1: Hyper-exponentially distributedcell dwell time

Considering that CDT has a hyper-exponential pdf of order *n* given by

For exponentially distributed UST and using (4), the CDF of the CHTh can be expressed as follows

This expression corresponds to a hyper-exponential distribution of order n with phase parameters *μ + η _{i}* and probabilities

*P*

_{i}of choosing stage

*i*(for i = 1, …, n).

As the CDTr is hyper-exponentially distributed when CDT has hyper-exponential distribution, the CHTn is also hyper-exponentially distributed. In this case, the probability of choosing stage *i* (for *i* = 1, …, *n*) is given by

#### 3.2.2. Case 2: Coxian distributedcell dwell time

Considering that cell dwell time has an *m*-th order Coxian distribution (which diagram of phases is shown in Fig. 1) with Laplace transform of its pdf given by

where

(1- α_{i}) represents the probability of passing from the *i*-th phase to the (*i*+1)-th one.

For exponentially distributed UST and using (9), the Laplace transforms of the pdf of CHTh and CHTn are given by

Replacing (13) into (15), it can be written as

where

for i = 1, …, m. Then, CHTh has also a Coxian distribution of order m but with parameters (*μ + η _{i}*), for i = 1, …, m.

On the other hand, the Laplace transform of the residual cell dwell time can be shown to be given by

where

for *j* = 1, …, *m*(*m*+1)/2. Substituting (19) into (16), Laplace transform of CHTn can be written as

where

Equation (23) corresponds to the Laplace transform of a generalized Coxian pdf.

The above analytical results show that CHTh (CHTn) has the same probability distribution as CDT (CDTr) but with different parameters of the phases, probabilities of reaching the absorbing state after each phase, and probabilities of choosing each stage. The detailed derivation of the last entry of Tables I and II (i.e., when cell dwell time has generalized Coxian distribution) is addressed in (Corral-Ruiz et al., a, 2010).

### 3.3. Relationship between X c ( h ) and X c ( N )

Using (15) and (16) it is straightforward to show that the mean values of CHTn and CHTh are, respectively, given by

At this point, it is important to mention that authors in (Fang, 2001; Zeng et al., 2002) stated that, depending on the variance of CDT, the mean CHTn can be greater than the mean CHTh. However, it was neither explained nor discussed the physical reasons for this observed behavior. This behavior occurs because the residual cell dwell times tend to increase as the variance of cell dwell time increases, as it was explained above.

Using (25) and (26), the condition for which the mean CHTn is greater that the mean CHTh, that is,

can be easily found. This condition is given by

Thus, the relationship between the mean new and handoff call channel holding times is determined by the mean values of both CDT and UST and by the Laplace transform of the pdf of CDT evaluated at the inverse of the mean UST.

Finally, in a similar way, the squared coefficient of variation for CHTn and CHTh can be shown to be given, respectively, by

It can be shown that the n-th moments for new and handoff call channel holding times are given, respectively, by

## 4. Numericalresults and discussion

In this section, numerical results on how the distribution of cell dwell time (CDT) affects the characteristics of channel holding time (CHT) are presented. We use different distributions to model CDT, say, negative-exponential, constant (deterministic), Pareto with shape parameter α in the range (1, 2] (i.e., when infinite variance is considered), Pareto with α>2 (i.e., when finite variance is considered), log-normal, gamma, hyper-Erlang of order (2,2), hyper-exponential of order 2, and Coxian of order 2. Three different mobility scenarios for the numerical evaluation are assumed: *E*{**X**_{d}}=5*E*{**X**_{s}} (low mobility), *E*{**X**_{d}}=*E*{**X**_{s}} (moderate mobility), and *E*{**X**_{d}}=0.2*E*{**X**_{s}} (high mobility). In the plots of this section we use *E*{**X**_{s}}=180 s. In our numerical results, the effect of CoV and skewness of CDT on CHT characteristics is investigated. In the plots presented in this section, “HC” and “NC” stand for channel holding time for handoff calls (CHTh) and channel holding time for new calls (CHTn), respectively.

### 4.1. Cell dwell time distribution completely characterized by its mean value

Fig. 2 plots the mean value of both CHTn and CHTh versus the mean value of CDT when it is modeled by negative-exponential (EX), constant, and Pareto with 1<2 distributions. It is important to remark that all of these distributions are completely characterized by their respective mean values. As expected, Fig. 2 shows that, for the case when CDT is exponentially distributed, mean CHTn is equal to mean CHTh. An interesting observation on the results shown in Fig. 2 is that, irrespective of the mean value of CDT, there exists a significant difference between the mean value of CHTn when CDT is modeled as exponential distributed RV and the corresponding case when it is modeled by a heavy-tailed Pareto distributed RV (this behavior is especially true for the case when α=1.1). Notice, however, that this difference is negligible for the case when α=2 and high mobility scenarios (say, *E*{**X**_{d}}<50 s) are considered. Similar behaviors are observed if mean CHTh is considered. Consequently, for high mobility scenarios where CDT can be statistical characterized by a Pareto distribution with shape parameter close to 2, the exponential distribution represents a suitable model for the CDT distribution. Fig. 2 also shows that, for

a given value of the mean CDT and considering the case when CDT is Pareto distributed with α=1.1 (α=2), mean CHTn always is greater (lower) than mean CHTh. This behavior can be explained by the combined effect of the following two facts. First, as α comes closer to 1 (2), the probability that CDT takes higher values increases (decreases). This fact contributes to increase (reduce) the mean CHTh. Second, in general, new calls are more probable to start on cells where users spent more time and, as α comes closer to 1, this probability increases. This fact contributes to increase mean CHTn relative to the mean CHTh. Then, the combined effect is dominated by the first (second) fact as α comes closer to 2 (1). This leads us to the behavior explained above and illustrated in Fig.2. It may be interesting to derive the condition upon which the mean CHTn is greater than the mean CHTh when CDT is heavy-tailed Pareto distributed. This represents a topic of our current research.

### 4.2. Cell dwell time distribution completely characterized by itsfirst two moments

Fig. 3 plots the mean value of both CHTn and CHTh versus the CoV of CDT when it is modeled by Pareto with shape parameter α>2, lognormal, and Gamma distributions; all of them with mean value equal to 180 s. It is important to remark that all of these distributions are completely characterized by their respective first two moments. Fig. 3 shows that both mean CHTn and mean CHTh are highly sensitive to the type of distribution of CDT; this fact is especially true for *CoV*>2. Notice that, for the particular case when *CoV*=0, the mean values of both CHTn and CHTh are identical to the corresponding values for the case when CDT is deterministic with mean value equals 180 s, as expected. Fig. 3 also shows that, for values of CoV of CDT greater than 1 (1.2), mean CHTn is greater that mean CHTh when CDT is Gamma (log-normal) distributed. On the other hand, when CDT is Pareto distributed and irrespective of the value of its CoV, CHTh always is greater that mean CHTn. This behavior is mainly due to the heavy-tailed characteristics of the Pareto distribution.

### 4.3. Cell dwell time distribution completely characterized by itsfirst three moments

Figs. 4, 5, and 6 (7, 8, and 9) plot the mean value (CoV) of both CHTn and CHTh versus both the CoV and skewness of CDT when it is modeled by hyper-Erlang (2,2), hyper-exponential of order 2, and Coxian of order 2 distributions, respectively. It is important to remark that all of these distributions are completely characterized by their respective first three moments. Results of (Johnson & Taaffe, 1989; Telek & Heindl, 2003) are used to calculate the parameters of these distributions as function of their first three moments. In Figs. 4 to 9, two different values for the mean CDT are considered: 36 s (high mobility scenario) and 900 s (low mobility scenario). From Figs. 2, 5 and 6 the following interesting observation can be extracted. Notice that, for the case when CDT is modeled by either hyper-exponential or Coxian distributions and irrespective of the mean value of CDT, the particular scenario where skewness and *CoV* of CDT are, respectively, equal to 2 and 1, corresponds to the case when CDT is exponential distributed (in the exponential case mean CHTn and mean CHTh are identical).

On the other hand, Fig. 4 shows that the case when hyper-Erlang distribution with skewness equals 2 and *CoV* equals 1 is used to model CDT does not strictly correspond to the exponential distribution; however, the exponential model represents a suitable approximation for CDT in this particular case. From Figs. 4 to 9, it is observed that the qualitative behavior of mean and *CoV* of both CHTn and CHTh is very similar for all the phase-type distributions under study. The small quantitative difference among them is due to moments higher than the third one. Analyzing the impact of moments of CDT higher than the third one on channel holding time characteristics represents a topic of our current research.

From Fig. 10 is observed that the difference among the mean values of CHTn and CHTh is strongly sensitive to the CoV of the CDT, while it is practically insensitive to the skewness of the CDT. This difference is higher for the case when the CDT is modeled as hyper-exponential distributed RV compared with the case when it is modeled as hyper-Erlang distributed RV. Also, it is observed that this difference remains almost constant for the entire range of values of the CoV of the CDT.

Finally, in Fig. 11 the mean channel holding time for new and handoff calls considering the gamma, hyper-Erlang (2,2), hyper-exponential of order 2, and Coxian of order 2 distributions for the cell dwell time are shown for different values of the coefficient of variation. The numerical results shown in Fig. 11 are obtained by equaling the first three moments of the different distributions to those of the gamma distribution. From Fig. 11, it is observed that for the hyper-exponentialand Coxian distributions practically the same results are obtained for the mean channel holding time for both new and handoff calls. The differences among the other distributions are due to the fact that they differ on the higher order moments. To show this, the forth standardized moment (i.e., excess kurtosis) of the different distributions is shown in Fig. 12 for different values of the coefficient of variation, equaling the first three moments of the different distributions to those of the gamma distribution. From Fig. 12, it is observed that the hyper-exponentialand Coxian distributions practically have the same value of excess kurtosis but this differs for that of the gamma and hyper-Erlang distributions. The gamma distribution shows the more different value of the excess kurtosis and, therefore, for this distribution the more different values of the mean channel holding times in Fig. 11 are obtained. Then, it could be necessary tocapture more than three moments, even though the lower ordermoments dominate in importance. Similar conclusion was drawn in (Gross &Juttijudata, 1997).

## 5. Conclusions

In this Chapter, under the assumption that unencumbered service time is exponentially distributed, a set of novel general-algebraic equations that examines the relationships between cell dwell time and residual cell dwell time as well as between cell dwell time and new and handoff channel holding times was derived. This work includes relevant new analytical results and insights into the dependence of channel holding time characteristics on the cell dwell time probability distribution. For instance, we found that when cell dwell time is Coxian or hyper-exponentially distributed, channel holding times are also Coxian or hyper-exponentially distributed, respectively. Also, our analytical results showed that the mean and coefficient of variation of the new and handoff call channel times depend on Laplace transform and first derivative of the Laplace transform of the probability density function of cell dwell time evaluated at the inverse of the mean unencumbered service time as well as on the mean of both cell dwell time and unencumbered service time. Additionally, we derive the condition upon which the mean new call channel holding time is greater than the mean handoff call channel holding time. Similarly, the condition upon which the mean residual cell dwell time is greater than the mean cell dwell time was also derived. To the best authors’ knowledge, this phenomenon that may seem to be counterintuitive has been explained and mathematically formulated in this Chapter. We believe that the study presented here is important for planning, designing, dimensioning, and optimizing of mobile cellular networks.