Downlink Coverage Analysis for a Finite 3D Wireless Network of Unmanned Aerial Vehicles

# Downlink Coverage Analysis for a Finite 3D Wireless Network of Unmanned Aerial Vehicles

## Abstract

In this paper, we consider a finite network of unmanned aerial vehicles (UAVs) serving a given region. Modeling this network as a uniform binomial point process (BPP), we derive the downlink coverage probability of a reference receiver located at an arbitrary position on the ground assuming Nakagami- fading for all wireless links. The reference receiver is assumed to connect to its closest transmitting node as is usually the case in cellular systems. After deriving the distribution of distances from the reference receiver to the serving and interfering nodes, we derive an exact expression for downlink coverage probability in terms of the derivative of Laplace transform of interference power distribution. In the downlink of this system, it is not unusual to encounter scenarios in which the line-of-sight (LOS) component is significantly stronger than the reflected multipath components. To emulate such scenarios, we also derive the coverage probability in the absence of fading from the results of Nakagami- fading by taking the limit . Using asymptotic expansion of incomplete gamma function, we concretely show that this limit reduces to a redundant condition. Consequently, we derive an accurate coverage probability approximation for this case using dominant interferer-based approach in which the effect of dominant interferer is exactly captured and the residual interference from other interferers is carefully approximated. We then derive the bounds of the approximate coverage probability using Berry-Esseen theorem. Our analyses reveal several useful trends in coverage probability as a function of height of the transmitting nodes and the location of reference receiver on the ground.

## 1Introduction

With significant advancements in the drone technology, like increased payload capacity, longer average flight time, better power management techniques, and the capability to harvest solar energy, unmanned aerial vehicles (UAVs) can serve a multitude of purposes such as surveillance, localization and communication, making them a flexible solution to augment and enhance the capabilities of the current cellular systems. They provide an especially attractive solution to provide connectivity in the wake of disasters and accidents, which may completely cripple the terrestrial networks due to damaged equipment and/or loss of power [2]. In general, UAVs provide a realistic solution in scenarios where there is a temporary need for network resources. These could include first responder situations, such as the one discussed above, or even usual civilian scenarios, such as football games or concerts. In order to provide short-term connectivity in such scenarios, temporary deployment of UAVs may be faster and more cost-effective compared to the temporary installation of conventional base stations. They are also currently being investigated as a possible candidate for providing ubiquitous connectivity in remote areas that lack traditional cellular infrastructure. While there is no doubt about the deployment flexibility and general benefits of UAVs, their performance in terms of the coverage and capacity provided to the terrestrial users is not quite well understood. This is especially true for a realistic use case of finite UAV networks, where we have a given number of UAVs serving users in a given region (such as a city). In this paper, we use tools from stochastic geometry to derive downlink signal-to-interference ratio distribution for this setup, which immediately provides useful insights into the coverage performance of the resulting three-dimensional network. Several intermediate results derived in this paper are also of general interest for the analysis of finite wireless networks.

### 1.1Motivation and Related Work

The improvements in payload capacity and prolonged flight times have enabled the commercial use of UAVs, especially for communication purposes. UAV networks differ significantly from conventional wireless networks in terms of the mobility, energy constraints, as well as the propagation conditions. This has stimulated interest in the design of application-oriented protocols for the effective utilization of aerial networks [4]. For instance, a cluster-based protocol, which improves the resilience to frequent link failures resulting from the motion of UAVs, has been proposed in [5]. The flexibility offered by the mobility of UAVs has motivated a lot of algorithmic research efforts towards finding efficient trajectories and deployment strategies aimed at optimizing different network resources [7]. For instance, an algorithm to optimize transmit power and frequency spectrum for autonomous self-deployment was proposed in [7]. An adaptive algorithm for adjusting the UAV heading was proposed in [8] to improve the uplink performance and minimize mutual interference. An approach to optimize the altitude of UAVs to maximize coverage on the ground was proposed in [9]. The performance of UAVs acting as relays between terrestrial users and base stations was investigated in [10]. The problem of efficient placement of UAVs with slightly different objectives was studied in [11].

Another direction of research, which is somewhat complementary to the one discussed above, is to develop techniques for the realistic system-level analysis of UAV networks. As is the case in terrestrial networks, such as cellular networks, these techniques can then be used to compare the performance of different deployment strategies and to benchmark their performance against standard baselines. In the case of UAV networks, the system-level performance has mostly been studied through field tests and simulations [14]. For instance, in [15] the outage time and average goodput were compared for different routing algorithms using real-world experiments. While field-tests or simulations can provide initial insights into the behavior of the network, these methods are usually not scalable when the number of simulation parameters is large. One way of reducing the dimensionality of such problems is to endow the locations of the nodes with a distribution, which additionally allows the use of powerful tools from stochastic geometry to derive easy-to-use expressions for key performance metrics. While stochastic geometry has already emerged as a preferred tool for the analysis of ad hoc and cellular networks [17], its potential has not yet been exploited for the analysis of UAV networks. One relevant prior art is [18], which studies the co-existence of a device-to-device (D2D) communication network with a single UAV. In this work, we will develop the first comprehensive model aimed at the downlink analysis of a finite multi-UAV network using tools from stochastic geometry.

While infinite homogeneous Poisson Point Process (PPP) has become a canonical model for the spatial locations of terrestrial base stations [17], it is not quite suitable for UAV networks, especially when a given number (likely small) of UAVs is deployed to cover a given finite region. For such scenarios, a simple yet reasonable model for the spatial distribution of UAVs is the homogeneous binomial point process (BPP) [19]. More sophisticated models incorporating inter-point interaction are usually far less tractable. While BPP has not yet been used for the analysis of UAVs, it has received significant attention for the analysis of terrestrial networks with a given number of nodes. Until recently, however, the analysis was focused on ad hoc networks, in which a given number of nodes were assumed to be distributed uniformly at random in a circular region with the reference receiver located at the center of the circle [21]. The outage probability for this reference receiver is then derived assuming that it is served by a reference transmitter located at a fixed distance (not a part of the BPP). In order to model cellular systems meaningfully with this setup, two key generalizations are required: (i) reference receiver can lie anywhere in the region, and (ii) serving base station for the reference receiver will be chosen from the BPP itself. For the latter, it is reasonable to assume that the reference receiver is served by the closest base station from the BPP. The exact analysis of this finite cellular network setup was done very recently in [26]. An approximate analysis of a related setup also appears in [27]. Building on the distance distributions derived in [26], we will perform downlink analysis for an arbitrarily located user on the ground that is served by a finite network of UAVs. In addition to providing the first such system-level analysis of a finite UAV network, several intermediate results provide constructs that are more generally applicable to the analysis of finite wireless networks. With this brief discussion, we now provide a precise summary of our contributions.

### 1.2Contributions

Modeling of finite three-dimensional network. We develop a general framework for the downlink coverage analysis of finite three-dimensional networks under a fairly general channel fading model. In particular, we consider a finite network of a given number of UAVs whose locations are modeled as a uniform BPP in a plane at a fixed altitude above the ground. The reference receiver is assumed to be located at some arbitrary position on the ground. As noted above, we assume that the reference receiver connects to its closest transmitting UAV node. We then derive the distribution of distances from the reference receiver to the serving and interfering nodes. As the air-to-ground channel models are still under active investigation, we assume Nakagami- fading, which allows us to control the severity of multi-path fading (to emulate variety of scenarios) while retaining analytical tractability. In some deployment scenarios, the LOS component may be significantly stronger than the reflected multipath components, which means it may be reasonable to study system performance in the absence of small-scale fading. Such scenarios, which will henceforth be referred to as “no-fading” environments, are also studied in detail.

Coverage probability. We derive an exact expression for coverage probability of the reference receiver located at some arbitrary position on the ground for Nakagami- fading channel in terms of the Laplace transform of interference power distribution. While this problem has been studied in [25] for terrestrial networks, the complexity and form of the final expression forbids any further analysis or simplifications. In this paper, we develop an approach that explores the possibility of deriving the results for no-fading case from the results of Nakagami- fading by applying the limit . We use well-established mathematical results, especially asymptotic expansion of incomplete gamma function, in the evaluation of this non-trivial limit. Quite interestingly, we discover that this limit renders a redundant condition for coverage probability, thereby not yielding an explicit expression for the no-fading case. As a result, we derive a simple yet accurate approximation for coverage probability using dominant interferer-based approach in which the effect of dominant interferer is accurately captured and the aggregate interference from the rest of the interferers is carefully approximated. We then derive the bounds of the approximate coverage probability using Berry-Esseen theorem (BET) [29].

Performance analysis. We analyze the trends in coverage probability for different system parameters such as the altitude of the UAVs, the location of the reference receiver, the path-loss exponent of the channel, and the Nakagami- fading parameter. We demonstrate that the coverage probability of the reference receiver degrades as the altitude of the transmitting nodes increases when the area over which the transmitters are scattered remains unchanged. We also observe that the coverage probability of the reference receiver increases as the path-loss exponent of the channel increases. These observations offer useful guidelines for the system design. Several intermediate mathematical results are of general interest to the analysis of finite networks.

## 2System Model

We consider a network of transmitting devices (UAVs) uniformly distributed in a finite area forming a BPP. While the devices can be strategically placed to optimize the network utility, in the absence of exact traffic patterns, these optimal locations are not known, which justifies the BPP assumption. As discussed in the previous section, BPP is a finite-network analogue of a PPP, which is a popular model for infinite networks. This is a simple yet reasonable first step towards comprehensive understanding of these networks. The locations of the devices are uniformly distributed in a disk of radius centered at , as depicted in Fig. ?. For simplicity, we assume that all the devices are positioned at the same height . The locations of the devices are denoted by . The distance of the node from is denoted by for . The sequence of distances is unordered, which means that the indices are assigned arbitrarily to the nodes. For this setup, we perform downlink coverage analysis for a reference receiver located at some arbitrary position on the ground, at a distance of from the origin . For brevity, the reference receiver will henceforth be referred to as only the receiver. Since the point process is invariant to the orientation of the axes, we can assume, without loss of generality, that the receiver lies on the -axis, i.e., the location of the receiver is . The distance between the receiver and the projection of the location of transmitting device onto the ground plane is denoted by , as shown in Figure 2. The angular separation between the receiver position and the projection of the device location is denoted by , which is uniformly distributed in the range . Note that and are related by the cosine rule. The receiver is assumed to connect to its closest transmitter from . The unordered set of distances from the receiver to the transmitters is denoted by . The ordered set of distances is denoted by , where is the distance between the receiver and the closest transmitter to the receiver. From this set, the serving distance is denoted by . The distance to the closest interfering node is denoted by and the unordered set of distances between the receiver and remaining interferers is denoted by .

We assume that all the nodes transmit at the same power. For notational simplicity, we assume that the thermal noise is negligible as compared to the interference experienced at the receiver and is hence ignored. Owing to the limited knowledge of air-to-ground channel models for low altitude platforms, we choose Nakagami- fading, which is a generalized model that mimics various fading environments. We denote the Nakagami- fading parameter for the serving link and interfering links by and , respectively. We restrict the values of to integers for analytical tractability. The signal-to-interference ratio (SIR) at the receiver is

where is the path-loss exponent, and and are channel fading gains corresponding to serving and interfering links. The channel gains and follow a gamma distribution with the probability density function (PDF) [30]

## 3Coverage Probability

In this section, we derive the coverage probability for the setup introduced in the previous section. Using properties of gamma function, we will express it in terms of the derivative of Laplace transform of interference power distribution. We also attempt to obtain the coverage probability for no-fading channel from the results of Nakagami- fading by taking the limits and . Quite interestingly, we discover that this limit renders a redundant condition for coverage probability, thereby not yielding an explicit expression for the no-fading case. Hence, we provide an alternate approach to compute the coverage probability in this case by approximating the aggregate interference from devices other than the dominant interferer with a normal random variable and capturing the effect of dominant interferer exactly.

### 3.1Relevant Distance Distributions

We begin our analysis by characterizing the distribution of distances between the receiver and the transmitters in this subsection. While somewhat similar expositions about distance distributions can be found in [26] in the context of terrestrial networks, the distributions corresponding to the dominant interferer-based approach are unique to this paper because this paper is the first one to apply that approach to the analysis of finite cellular networks.

For a receiver located at the origin , this piece-wise expression for the PDF reduces to a simple expression, which is given in the following Corollary.

If the receiver is located at the origin, the above result reduces to a simple expression, which is given in the next Corollary.

We now derive the distribution of distance between the receiver and interferers, conditioned on the serving distance . This distribution will be useful in characterizing the interference experienced at the receiver.

For the case where the receiver is located at the origin, the conditional PDF of distances between the receiver and the interferers is given in the following Corollary.

We now characterize the joint distribution of serving distance and the distance between the receiver and the dominant interferer . This distribution holds the key to the derivation of approximate coverage probability using the dominant interferer approach (for the no-fading case) as discussed in the sequel. Note that the dominant interferer in the absence of fading is the second closest transmitter to the receiver. Since it is easy to visualize and understand these cases in two dimensions, we consider the projections of the transmitting nodes onto the ground plane. We will then use these results to derive the distributions in the actual three-dimensional setup. Let the distance from the receiver to the projections of locations of serving transmitter and dominant interferer be denoted by and , respectively. Using the same argument as presented in Lemma ?, the piece-wise nature of joint distribution of and can be attributed to the following three cases: (i) both the disks and are entirely contained in the disk , (ii) the disk is entirely contained in , while overlaps only partially with , and (iii) both the disks and overlap only partially with . Note that the converse of case (ii) is not valid since .

For the case where the receiver is at the origin, this result reduces to a simple expression, which is given in the next Corollary.

We now determine the distribution of the set of distances between the receiver and the interferers , conditioned on the serving distance , and the distance of the dominant interferer from the receiver . Let the distances from the receiver to the projections of remaining interferers be denoted by . The piece-wise distribution of these distances depends on the following four cases: (i) all the disks , , and are completely contained in , (ii) the disks and are completely contained in , while the disks overlap only partially with , (iii) the disk is completely contained in the disk , while and partially overlap with , and (iv) all the disks , , and overlap partially with . Note that other combinations are not valid because .

The following Corollary gives the conditional PDF of distances of interferers from the receiver for the special case where the receiver is located at the origin.

Now that we have determined the necessary distance distributions, we will proceed to derive the coverage probability for the receiver under Nakagami- fading in the following subsection.

### 3.2Coverage Probability under Nakagami-m Fading Channels

The coverage probability is formally defined as the probability with which the SIR at the receiver exceeds a pre-determined threshold necessary for a successful communication. Recall that the SIR at the receiver in Nakagami- fading is given by

where is the path-loss exponent, and and are channel fading gains with parameters and , respectively. We denote the set of gains for the interfering links by and the set of distances of the interferers from the receiver by . We first calculate the coverage probability, conditioned on , as a derivative of conditional Laplace transform of interference power distribution, which is given in the following Lemma.

For the case where the receiver is located at the origin, this result reduces to a simple expression, which is given in the next Corollary.

Using the Laplace transform of conditional interference power distribution, we derive the coverage probability in the following theorem.

In this subsection, we attempt to derive the coverage probability for a no-fading environment from the results of Nakagami- fading by applying the limits and . This approach has mostly been overlooked in the literature due to the complexity of the results of Nakagami- fading. A partial attempt to compute this limit was made in [31] where the limit is applied only on the interfering links but not on the desired link. We show that it is more challenging to evaluate the limit for the desired link . Our analysis will rely on an asymptotic expansion of incomplete gamma function.

From our approach, we observe that it is convenient to take this limit in one of the intermediate steps in the derivation of under Nakagami- fading. Substituting in , we get the conditional coverage probability in Nakagami- fading as follows:

where (a) follows from switching the order of differentiation and expectation by applying Dominated Convergence Theorem (DCT) [32]. Therefore, the coverage probability conditioned on and can now be written as

The results obtained upon applying the limits on the conditional coverage probability are given in the following Theorem.

We notice that the condition under which the coverage probability converges to 1 is , i.e., . This is nothing but the condition of in no-fading channels. Therefore, quite interestingly, we obtain only a redundant condition when the coverage probability in no-fading channels is evaluated from the results of Nakagami- fading. To the best of our understanding, this insight has not been reported before in the context of the limiting case of Nakagami- fading. In the next subsection, we provide an accurate approximation to compute coverage probability in the absence of fading.

### 3.4Dominant Interferer Approach

In this method, we capture the effect of dominant interferer exactly and approximate the aggregate interference from rest of the interferers to a Gaussian random variable. The at the receiver in the absence of fading is

Let . Since the distance of the interferers from the receiver conditioned on and are i.i.d., the terms that constitute the sum are also conditionally i.i.d. Therefore, by central limit theorem (CLT), the sum of i.i.d. random variables can be approximated by a normal random variable, whose mean and variance are given by the following Lemmas.

While these integrals can not be reduced to closed-form, it is easy to evaluate them numerically. However, for a receiver located at the origin, these expressions can be simplified to closed-form expressions given in the following Corollaries.

Using these results, we derive an accurate approximation for coverage probability in the following Theorem.

### 3.5Bounds of Coverage Probability Approximation

In this subsection, using BET, we analyze the tightness of the coverage probability approximation proposed in the previous subsection. BET gives a bound on the maximal deviation of the normal distribution from the true distribution in terms of the moments of the distribution. By BET, for a sequence of random variables , with , , and sample mean , the error between the actual distribution of the random variable and the standard normal distribution is bounded by , i.e.,

where is a constant. The best known estimate of is [37]. We now rewrite the expression for interference in a form that will allow us to conveniently apply BET. Let and where is the mean interference from Lemma ?. Therefore, we have and , which is given by Lemma ?. The third moment of absolute value of , , is computed in the following Lemma.

While it is difficult to get a simple closed-form expression for the above result, it can be easily evaluated numerically. Using this result, along with the first and second moments of , we derive the bounds of coverage probability in the following Theorem.

While these bounds are loose for very small values of , we observe that the error between the normal distribution and true distribution decreases at the rate . In other words, as the number of transmitting devices grows to a large value, the approximate coverage probability converges to the actual value. However, in the next section, we show that this approximation is surprisingly accurate even for a small number of nodes in the network.

## 4Results and Discussion

In this section, we validate our analytical results by comparing the theoretical coverage probabilities with the simulation results for a finite network of UAVs. In addition to providing useful design insights, we also discuss the applicability of our proposed analytic approaches to an urban setting in which the visibility of some UAVs is obstructed/blocked by buildings.

### 4.1Numerical Results

We simulate a finite network of UAVs with , uniformly distributed in a disk of radius km. We evaluate the coverage probability for different system parameters and compare them with the theoretical results obtained in Theorems ? and ?. We find that our theoretical results match exactly with the simulations as shown in Fig. ?. The key factors that affect the coverage probability are: (i) channel fading parameter , (ii) path-loss exponent , (iii) the height of UAVs , and (iv) the distance of the receiver from the origin . We study the impact of each parameter on the coverage probability in the rest of this Subsection.

Impact of fading. We compute the coverage probability of the receiver as a function of SIR threshold for and . Note that is nothing but the no-fading scenario. The other parameters of the simulation are km, km, and . As expected, the variance of SIR decreases with the increase in . In other words, SIR starts concentrating as we move from Rayleigh fading () case to the no-fading () case.

Impact of path-loss exponent. We study the impact of path-loss exponent on coverage in Figure 5, where we plot the coverage probability as a function of threshold for different values of . The simulations are run for km, , and km. It can be observed that the coverage probability degrades with a decrease in the path-loss exponent. While reducing increases the received power of the desired signal, it also increases the interference power, thereby degrading the overall SIR and hence the coverage probability.

Impact of height. We compare the coverage probability of the receiver as a function of threshold for different values of the height of UAVs (2, 4, 6, and 8 km) . The other simulation parameters were fixed at , , and km. It can be observed from Fig. ? that the coverage probability deteriorates as the height of the UAVs increases. An increase in the height of UAVs increases the distance between the receiver and the transmitters. Intuitively, when viewed from a receiver that moves away from the transmitters, the separation between the serving and interfering nodes tends to diminish. This worsens the SIR and hence the coverage.

Impact of receiver distance from the origin. The impact of receiver distance from the origin on coverage probability can be studied from Figure 7 where we plot coverage as a function of for different values of . The other simulation parameters were km, , dB, and . It can be observed that the coverage probability varies significantly with the location of the receiver, which highlights the importance of assuming arbitrarily located receiver.

### 4.2Applicability of the Proposed Results to a Relevant Urban Model

As discussed already, the main technical objective of this paper was to develop a comprehensive framework for the analysis of a reasonable canonical model for finite UAV networks. Before concluding this Section, we demonstrate that this canonical model can be enriched to some extent while retaining its tractability. In particular, we incorporate the effect of shadowing/blocking due to high-rises that will be prominent in urban regions, especially at higher frequencies. The system setup is illustrated in Fig. ?. The blockages may result in attenuation of the received signal and hence affect the coverage probability. We begin with the spatial model considered in this paper, where transmitting devices are uniformly distributed in a disk at an altitude above the ground. Due to the presence of buildings, all the UAVs may not be visible at the receiver. We denote the number of visible UAVs by . In our analytical treatment, we assume that the LOS paths from the receiver to each UAV are blocked independently of each other. In other words, we ignore the correlation in blocking introduced by the spatial distribution of buildings in the area. As a result, can be modeled as a binomial random variable. This independence assumption will be validated through a numerical comparison in Figure 9.

For simplicity, we limit our discussion to direct path propagation and ignore multi-path fading. We assume that if a UAV is hidden behind a building, its signal is attenuated by a fixed factor . While we can easily extend this discussion to a more general setup, this simple scenario is sufficient to fix the key ideas. For this setup, the received signal power from the transmitter, located at a distance from the receiver is , where if the UAV is visible at the receiver and otherwise. Recall that the distribution of is given by Lemma ?. Now conditioned on , we get two independent BPPs: (i) a BPP formed by visible UAVs, and (ii) a BPP formed by blocked UAVs. Conditional on , the coverage analysis can be performed following the proposed approach. We do not go into the mathematical details due to lack of space. Note that if , we get only one BPP (of visible UAVs), which reduces this setup to that of Theorem ?, which will be used for numerical comparisons below.

The key approximation made in the above analysis is the independent blocking assumption. We validate this assumption numerically. In particular, we simulate a urban scenario in MATLAB with 5 UAVs uniformly distributed over a circular area of radius of 10 km. We assume a uniform distribution of 50 buildings that are of dimensions m m m. We choose the simple case of , i.e., we receive signals from only those UAVs that are visible. For a receiver at a distance of 1 km from the origin and path-loss exponent of , we obtain the coverage probability from simulations. For each realization, we numerically obtain , which is used instead of in Theorem ? to obtain the conditional coverage probability analytically. Monte-Carlo simulations are used to average over . Note that the purpose of this comparison is to show that the independent blocking assumption is reasonable. This is quite evident in the results presented in Figure 9, where the simulation results are the ones obtained from actual numerical experiments without any assumptions, and the analytical result is obtained under independent blocking assumption. This discussion shows that the canonical setup introduced in this paper can be extended in many meaningful ways to study various aspects of UAV networks.

## 5Conclusion

In this paper, we have presented a comprehensive downlink coverage analysis for a finite three-dimensional wireless network formed by UAVs. Modeling the network of UAVs as a BPP, we characterized the distribution of distances from the receiver to the serving and interfering nodes. We first derived an exact expression for coverage probability for the reference receiver under independent Nakagami- fading channels in terms of the derivatives of the Laplace transform of interference power distribution. Using asymptotic expansion of incomplete gamma function, we showed that the coverage probability for an important special case of no-fading cannot be obtained explicitly as the limiting case of Nakagami- fading by taking the limit . As a result, we developed an alternate approach to compute the approximate coverage probability in which the interference from dominant interferer is modeled exactly and the distribution of residual interference from other interferers is approximated to a normal distribution by CLT. We then obtained the bounds of the approximate coverage probability using Berry-Esseen theorem, which quantifies the rate of convergence of the normal approximation to the true distribution. Our analysis revealed several useful performance trends in terms of the heights of the UAVs and channel propagation characteristics. We also discussed a possible extension of the proposed canonical model to a simple urban scenario in which the the UAVs are shadowed by high-rises.

This work has numerous extensions. The mathematical tools developed in the paper can be applied to the analysis of more general three-dimensional finite networks. The setup studied in this paper can also be extended to study the co-existence of UAV networks and terrestrial cellular networks. From modeling perspective, a useful direction of work wold be to develop more sophisticated but tractable three-dimensional spatial models that account for the curvature of the earth. In terms of performance evaluation, the proposed framework can be used to study other useful metrics besides coverage, such as throughput and energy efficiency.

### @.1Proof of Lemma

The cumulative distribution function (CDF) of each element of the sequence is

where is the CDF of the distance between the receiver and the projection of the location of transmitter onto the ground plane. The projections of the locations of the transmitters form a two dimensional BPP on the ground plane. Using the approach presented in [38], the conditional CDF of is computed as the the area of intersection of the disks and , divided by the area of the disk . Depending on the range of , there are two possible cases: (i) the disk is entirely contained in the disk , and (ii) partially overlaps with . Therefore, we obtain a piece-wise conditional CDF of as given below:

with

Substituting in , we obtain the conditional CDF of . This completes the proof.

### @.2Proof of Lemma

Since the receiver connects to the closest transmitter, the serving distance is by definition , where the distribution of is given in Lemma ?. The conditional CDF of can therefore be computed as

where (a) follows from the i.i.d nature of the set of distances . Differentiating the above expression w.r.t. , the PDF of the serving distance is obtained as

The PDF of can be obtained by substituting the results from Lemmas ? and ? in the above equation. This completes the proof.

### @.3Proof of Lemma

The joint density function of the ordered subset conditioned on the serving distance , and is

where (a) follows from the joint density function for the order statistics of a sample of size drawn from the distribution of , and (b) follows from the result derived in Lemma ?. Following the same argument presented in Lemma 3 in [26], we can say that indicates all possible permutations of the elements in the ordered set . Hence, by the joint density function for ordered set, the unordered set of distances are i.i.d. with PDF .

### @.4Proof of Lemma

For a sequence of i.i.d. random variables, , with each element characterized by PDF and CDF , the order statistics are random variables defined by sorting the realizations of the sequence in the increasing order. By order statistics [39], the joint PDF of the smallest two random variables is

In our case, the serving distance and the distance of the dominant interferer from the receiver are the smallest distances in the set . The joint PDF of the two distances is obtained by substituting the results from Lemmas ? and ? in . This completes the proof.

### @.5Proof of Lemma

The Laplace transform of interference power distribution conditioned on serving distance can be derived as follows:

where (a) follows from the independence of channel gains and the distances of interferers from the receiver, (b) follows from rewriting the expectation of product as the product of expectation owing to i.i.d. channel gains , (c) follows from the moment generating function (MGF) of gamma random variable , and (d) follows from conditionally i.i.d. distances of interferers from the receiver. Now, by applying the definition of mean and using the conditional PDF of from Lemma ?, we get the final result. This completes the proof.

### @.6Proof of Asymptotic Expansion of Incomplete Gamma Function

In order to prove the desired result, we first asymptotically estimate the following auxiliary integral as ,

where are taken to be real for simplicity. Using Taylor series expansion, we have

From multinomial expansion, it follows that

where

So, we can now write

Substituting in ,

With the substitution in the following integral, we get

We know by definition that

Here, the fractional powers of are to be understood as having principal values. Substituting in ,

Now we define,

Substituting and in ,

Using this key intermediate result, we will now prove our main result by deriving an asymptotic estimate for the incomplete gamma function as .

where (a) follows from the substitution . For , we observe that

For , this limit results in three sub-cases depending on relative to :

Using these limits in and substituting the resulting expression in , we obtain the asymptotic expansion of incomplete gamma function for the special case. Now, considering , we get

where,

Expanding the summation and ignoring the terms for , we obtain the same result for as follows:

We also know that [40]

Therefore, from and ,

### @.7Proof of Lemma

Since and , the conditional PDF of is given by

Now, by definition,

We know that is the smallest interference that could be caused by any of the transmitting nodes at the receiver. Since will always be smaller than the mean interference , . Therefore, we split the integral and change the limits accordingly. Substituting in and using the conditional distribution of from Lemma ?, we obtain the final expression.

### @.8Proof of Theorem

The coverage probability is given by

Let , where , and be the conditional CDF of . The coverage probability can now be written as

By BET, the CDF term in the integrand is related to the CDF of standard normal distribution as follows:

where . Substituting in and integrating the resulting inequality, we get the bounds of coverage probability.

### References

1. V. V. Chetlur and H. S. Dhillon, “Downlink coverage probability in a finite network of unmanned aerial vehicle (UAV) base stations,” in Proc., IEEE SPAWC, Jul. 2016.
2. A. Merwaday, A. Tuncer, A. Kumbhar, and I. Guvenc, “Improved throughput coverage in natural disasters: Unmanned aerial base stations for public-safety communications,” IEEE Veh. Technology Magazine, vol. 11, no. 4, pp. 53–60, Dec. 2016.
3. S. Kandeepan, K. Gomez, L. Reynaud, and T. Rasheed, “Aerial-terrestrial communications: terrestrial cooperation and energy-efficient transmissions to aerial base stations,” IEEE Trans. on Aerospace and Electronic Systems, vol. 50, no. 4, pp. 2715–2735, Oct. 2014.
4. I. Bekmezci, O. K. Sahingoz, and S. Temel, “Flying ad-hoc networks (FANETs),” Ad Hoc Netw., vol. 11, no. 3, pp. 1254–1270, May 2013.
5. B. Fu and L. A. DaSilva, “A mesh in the sky: A routing protocol for airborne networks,” in Proc., IEEE MILCOM, Oct. 2007.
6. S. Hayat, E. Yanmaz, and R. Muzaffar, “Survey on unmanned aerial vehicle networks for civil applications: A communications viewpoint,” IEEE Commun. Surveys and Tutorials, vol. 18, no. 4, Fourthquarter 2016.
7. H. Claussen, “Autonomous self-deployment of wireless access networks,” Bell Labs Tech. Journal, vol. 14, no. 1, pp. 55–71, Spring 2009.
8. F. Jiang and A. Swindlehurst, “Optimization of UAV heading for the ground-to-air uplink,” IEEE Journal on Sel. Areas in Commun., vol. 30, no. 5, pp. 993–1005, Jun. 2012.
9. A. Al-Hourani, S. Kandeepan, and S. Lardner, “Optimal LAP altitude for maximum coverage,” IEEE Wireless Commun. Letters, vol. 3, no. 6, pp. 569–572, Dec. 2014.
10. P. Zhan, K. Yu, and A. L. Swindlehurst, “Wireless relay communications with unmanned aerial vehicles: Performance and optimization,” IEEE Trans. on Aerospace and Electronic Systems, vol. 47, no. 3, pp. 2068–2085, Jul. 2011.
11. M. Mozaffari, W. Saad, M. Bennis, and M. Debbah, “Efficient deployment of multiple unmanned aerial vehicles for optimal wireless coverage,” IEEE Commun. Letters, vol. 20, no. 8, Aug. 2016.
12. R. I. Bor-Yaliniz, A. El-Keyi, and H. Yanikomeroglu, “Efficient 3-D placement of an aerial base station in next generation cellular networks,” in Proc., IEEE Intl. Conf. on Commun. (ICC), May 2016.
13. D. T. Ho, E. I. Grøtli, P. B. Sujit, T. A. Johansen, and J. B. Sousa, “Cluster-based communication topology selection and UAV path planning in wireless sensor networks,” in Proc., IEEE Intl. Conf. on Unmanned Aircraft Systems (ICUAS), May 2013.
14. W. Guo, C. Devine, and S. Wang, “Performance analysis of micro unmanned airborne communication relays for cellular networks,” in Proc., IEEE Intl. Symposium on Commun. Systems, Networks, and Digital Sig. Proc. (CSNDSP), Jul. 2014.
15. S. Rosati, K. Krużelecki, G. Heitz, D. Floreano, and B. Rimoldi, “Dynamic routing for flying ad hoc networks,” IEEE Trans. on Veh. Technology, vol. 65, no. 3, Mar. 2016.
16. R. C. Palat, A. Annamalau, and J. R. Reed, “Cooperative relaying for ad-hoc ground networks using swarm UAVs,” in Proc., IEEE MILCOM, Oct. 2005.
17. J. G. Andrews, A. K. Gupta, and H. S. Dhillon, “A primer on cellular network analysis using stochastic geometry,” arXiv preprint, Oct. 2016, available online: arxiv.org/abs/1604.03183.
18. M. Mozaffari, W. Saad, M. Bennis, and M. Debbah, “Unmanned aerial vehicle with underlaid device-to-device communications: Performance and tradeoffs,” IEEE Trans. on Wireless Commun., vol. 15, no. 6, pp. 3949–3963, Jun. 2016.
19. M. Haenggi, Stochastic Geometry for Wireless Networks.1em plus 0.5em minus 0.4emCambridge University Press, 2013.
20. S. Srinivasa and M. Haenggi, “Distance distributions in finite uniformly random networks: Theory and applications,” IEEE Trans. on Veh. Technology, vol. 59, no. 2, pp. 940–949, Feb. 2010.
21. J. Chen, M. Ding, and Q. T. Zhang, “Interference statistics and performance analysis of mimo ad hoc networks in binomial fields,” IEEE Trans. on Veh. Technology, vol. 61, no. 5, pp. 2033–2043, Jun. 2012.
22. K. Venugopal, M. C. Valenti, and R. W. Heath, “Interference in finite-sized highly dense millimeter wave networks,” in Proc., Information Theory and its Applications Workshop (ITA Workshop), Feb. 2015.
23. Z. Gong and M. Haenggi, “Interference and outage in mobile random networks: Expectation, distribution, and correlation,” IEEE Trans. Mobile Computing, vol. 13, no. 2, pp. 337–349, Feb. 2014.
24. D. Torrieri and M. C. Valenti, “The outage probability of a finite ad hoc network in Nakagami fading,” IEEE Trans. on Commun., vol. 60, no. 11, pp. 3509–3518, Nov. 2012.
25. J. Guo, S. Durrani, and X. Zhou, “Outage probability in arbitrarily-shaped finite wireless networks,” IEEE Trans. on Commun., vol. 62, no. 2, pp. 699–712, Feb. 2014.
26. M. Afshang and H. S. Dhillon, “Fundamentals of modeling finite wireless networks using binomial point process,” submitted to IEEE Trans. on Wireless Commun., 2016, available online: arxiv.org/abs/1606.04405.
27. S. Banani, A. Eckford, and R. Adve, “Analyzing the impact of access point density on the performance of finite-area networks,” IEEE Trans. on Commun., vol. 63, no. 12, pp. 5143–5161, 2015.
28. Z. Khalid and S. Durrani, “Distance distributions in regular polygons,” IEEE Trans on Veh. Technology, vol. 62, no. 5, pp. 2363–2368, 2013.
29. A. C. Berry, “The accuracy of the gaussian approximation to the sum of independent variates,” Trans. of the American Mathematical Society, vol. 49, no. 1, pp. 122–136, 1941.
30. D. Wackerly, W. Mendenhall, and R. Scheaffer, Mathematical Statistics with Applications.1em plus 0.5em minus 0.4emDuxbury, 2002.
31. R. Tanbourgi, H. S. Dhillon, J. G. Andrews, and F. K. Jondral, “Dual-branch MRC receivers under spatial interference correlation and Nakagami fading,” IEEE Trans. on Commun., vol. 62, no. 6, pp. 1830–1844, Jun. 2014.
32. P. Billingsley, Probability and Measure.1em plus 0.5em minus 0.4emWiley, 1995.
33. M. Abramowitz and I. Stegun, Handbook of Mathematical Functions. 1em plus 0.5em minus 0.4emDover Publications, 1965.
34. F. Tricomi, “Asymptotische eigenschaften der unvollständigen gammafunktion,” Mathematische Zeitschrift, vol. 53, pp. 136–148, 1950.
35. N. M. Temme, “The asymptotic expansion of the incomplete gamma functions,” SIAM J. Math. Anal., vol. 10, no. 4, pp. 757–766, 1979.
36. W. Gautschi, “The incomplete gamma functions since Tricomi,” in In Tricomi’s Ideas and Contemporary Applied Mathematics, Atti dei Convegni Lincei, n. 147, Accademia Nazionale dei Lincei, 1998, pp. 203–237.
37. I. Shevtsova, “On the absolute constants in the Berry-Esseen type inequalities for identically distributed summands,” arXiv preprint arXiv:1111.6554, 2011.
38. A. M. Mathai, An introduction to geometrical probability: Distributional aspects with applications.1em plus 0.5em minus 0.4emCRC Press, 1999, vol. 1.
39. M. Ahsanullah and V. Nevzorov, Order Statistics: Examples and Exercises.1em plus 0.5em minus 0.4emNova Science Publishers, 2005.
40. W. Feller, An introduction to probability theory and its applications: Volume I.1em plus 0.5em minus 0.4emJohn Wiley & Sons London-New York-Sydney-Toronto, 1968, vol. 3.