I Introduction
Due to the stationary locations and high cost of traditional macro base stations (BSs), it is extremely arduous to provide ubiquitous coverage via terrestrial cellular networks, especially for critical applications, e.g. disaster rescue, firefighting, reconnaissance, etc [2, 3, 4]. Therefore, an effective solution is urgently needed for the next generation of wireless networks. Under these circumstances, unmanned aerial vehicles (UAVs) become increasingly popular, owing to their flexibility and autonomy [5]. In terms of the frequency band used for UAV networks, two main candidates have been proposed: sub6 GHz and millimeter wave (mmWave) [2, 4, 6]. Benefited by mature wireless techniques and tolerance to blockages, the existing sub6 GHz in traditional networks can be reused effortlessly in a few denseobstacle environments [2]. However, in most of the cases, UAVs are able to establish lineofsight (LOS) connections to users by adjusting their locations and altitudes [4]. Under this condition, mmWave bands become the best choice for UAVaided networks due to its larger available bandwidth and consequent higher data rate transmission in LOS than sub6GHz [7]. Moreover, with the aid of large antenna scales deployed at mmWave devices, sharp directional beams can be generated to increase network capacity [8] and mitigate the Doppler spread [4].
In order to analyze the average performance of wireless networks, stochastic geometry has rekindled the strong interest of academia [9]. More specifically, the locations of transceivers are distributed according to different random spatial point processes, which are capable of statistically depicting the nature of wireless systems [10]. One popular point pattern is Poisson point process (PPP), which has been widely adopted in recent works for modeling network nodes such as BSs and multiple pieces of user equipment [8, 10, 11]
. The main advantage of this pattern is that all nodes are uniformly distributed in a plane. However, in practical scenarios, most wireless systems have a clustered property. For example, in devicetodevice (D2D) communications or cognitive networks, users with the same demands are frequently gathered to form a cluster and the PPP model fails to accurately describe this property. As a result, another point pattern named Poisson cluster process (PCP), which includes multiple clusters, becomes increasingly attractive in the literature
[12, 13].Ia StateoftheArt and Motivation
The evaluation of the performance of UAV networks plays a vital role in proposing new relative protocols and designing efficient network structures. Accordingly, AlHourani et al. [14] provided a theoretical approach to obtain the optimal altitude of UAVs, with an aim to maximize the coverage on the ground. However, this primary work was based on deterministic locations. Unfortunately, acquiring the accurate location information is a difficult and costly task. Therefore, researchers turned to analyze randomly deployed spatial models for UAV networks. Mozaffari et al. [15] investigated the average coverage probability and throughput in a UAVaided network with underlaid D2D communications, where the D2D users were modeled as a homogeneous PPP and the downlink users were uniformly distributed in a finite area forming a variant of PPP, known as binomial point process (BPP). However, this article only considered a singleUAV scenario. Chetlur and Dhillon [16] extended this work to a multiUAV case, where the locations of UAVs obey a BPP. With the aid of this framework, the coverage performance of downlink transmissions was analysed. Recently, Zhang and Zhang in [17] introduced a threedimensional (3D) PPP incorporating the adjustable altitude, in order to model the locations of UAVs. Analytical expressions for the coverage probability of drone smallcells were derived.
In terms of the frequency band utilized for UAV networks, the aforementioned research contributions [14, 15, 16, 17] focused on the sub6 GHz scenarios. However, when considering mmWave scenarios, blockage environment, antenna patterns, and fading channel model should be carefully modified. Instead of using the ray tracing method as discussed in [18], Bai et al. [19] leveraged the random shape theory to propose a mathematical blockage model for high frequencies. This model was approximated by a fixed lineofsight (LOS) disc in [8], which had acceptable accuracy with high calculation efficiency. By invoking the altitude information of devices, the conventional twodimensional (2D) blockage model was extended to the 3D case with the aid of Rayleigh distributed buildings [20]. In addition to the widely used omnidirectional antenna pattern in sub6 GHz, Zhu et al. [21] introduced a 3D sectorized beamforming pattern for mmWave antennas to depict the sharp directional beam. Furthermore, due to the huge difference between LOS and nonLOS (NLOS) mmWave transmissions, the Nakagami fading channel was preferable in mmWaveenabled networks. This smallscale fading model was considered in [11, 22, 21, 8, 13].
In UAVaided networks, for the communications between terrestrial BSs and UAVs, the locations of transceivers can be modeled as two independent PPPs to capture the randomness of the networks. However, for communications between UAVs and users, in most of the cases one UAV is dedicated to serve a cluster of users with the same requirements. To this end, PCP becomes a popular pattern for modeling these scenarios. Furthermore, Ganti and Haenggi [12] first proposed PCP to model the node locations in clustered wireless ad hoc networks. In this work, two specific models, i.e., Matern cluster processes (MCPs) and Thomas cluster processes (TCPs), were provided. Then, several networks were studied with these two types of PCPs, e.g. cognitive networks [23], heterogeneous networks [24], and devicetodevice (D2D) networks [13]. To the best of our knowledge, the research in the context of UAV networks by utilizing PCP is still in its infancy.
When considering the altitude of UAVs and the cluster property of users, a unified 3D spatial framework including both uplink and downlink phases is needed, which can be effortless adjusted to fit most practical applications. Additionally, the evaluation of average coverage performance in UAV networks, especially for mmWave scenarios, is still at the very early stage. These two factors are key motivations for this article.
IB Contributions and Organization
We propose a unified spatial framework for UAVaided networks, where a UAV first collects messages from one user in an uplink phase and then it retransmits the desired information to the requiring BS in a downlink phase. The transceivers in the downlink phase and uplink phase are distributed with the aid of a PPP and a PCP, respectively. Particularly, two commonly used PCPs, namely TCP and MCP, are investigated, and mmWave communications are considered as well. In summary, the main contributions of this paper are as follow:

We propose a novel 3D analytical framework for UAVaided networks by invoking a tractable 3D blockage model and a 3D sectorized antenna pattern, in which the altitude of all transceivers and the cluster property of users are jointly considered. Additionally, both uplink and downlink phases are assessed in order to model the entire transmission process of UAVaided communications.

We investigate the distributions of communication distances in the downlink phase. Based on these distributions, closedform expressions for the coverage probability are derived to enhance the evaluation efficiency. Regarding the uplink phase, we derive a general expression for coverage probabilities. Moreover, this expression is able to characterize TCP as well as MCP. Note that the number of users in each cluster can be a constant or to follow a Poisson distribution depending on different network requirements. We theoretically demonstrate that these deployment scenarios perform similarly when the number of interferers is large.

We provide a practical cooperation transmission to model the behaviours of UAVs from the uplink phase to the downlink phase. Three typical applications of UAVs are examined with the aid of the unified framework, and thus validating the flexibility of our work. In order to study the performance of the entire communication process, we also derive tractable expressions for the system coverage probability.

We show that, 1) when the geographical size of clusters and the interfering number of users are small, the intercluster interference can be ignored to simplify the analysis; 2) the effect of thermal noise is negligible in the proposed UAV networks. In terms of mmWave communications, NLOS transmissions can be ignored as well, especially in a lowdensity blockage environment; 3) there exists an optimal altitude of UAVs for achieving the maximum system coverage probability; and 4) a large antenna scale is able to enlarge its main beam gain and narrow the beamwidth for compensating path loss, thereby enhancing coverage performance.
The rest of this paper is organized as follows: Section II proposes the network model including the PPP distributed downlink phase and the PCP distributed uplink phase. Section III first describes the coverage performance in the proposed UAVaided networks. Then, three practical applications are studied with the aid of the unified framework. Section IV provides the validating and numerical results. Section V presents the conclusions.
Ii Network Model
We introduce a unified spatial framework for UAVaided networks by integrating mmWave communications, which describes a general case where a typical BS is located far from the served users, and hence, UAVs are deployed to establish the desired communication link. Due to the mobility of UAVs, we assume that the serving UAV is able to fly from the user region to the BS region. Therefore, one corresponding UAV receives the message from a transmitting user and then it flies toward the typical BS to forward this message^{1}^{1}1In this paper, we focus on the coverage performance of the UAVaided mmWave networks. The analysis of delay caused by the flight process will included in our future work..In addition to the flight assumption, more practical cooperation transmission scenarios are discussed at the end of this section. The entire transmission process can be divided into two phases: 1) Uplink Phase, where the transmitting user uploads the desired message to the corresponding UAV; and 2) Downlink Phase, where the typical BS downloads the required message from the corresponding UAV. These two phases are processed in different time slots. The details are presented in the following.
Iia Spatial Distribution
Different spatial distributions are assumed for the two phases. For ease of understanding, we first introduce the downlink and then the uplink phase.
IiA1 Downlink Phase
In this phase, we consider multiple transceivers. To capture the randomness of networks, we assume that macro BSs and UAVs are modeled as two independent PPPs with density and , denoted by and , respectively. Regarding the altitude of each device, UAVs are hovering at a height of , while the receiving antennas of macro BSs are located at an altitude , as illustrated in Fig. 1.
IiA2 Uplink Phase
For the uplink phase, since users with similar requests frequently stay together, we utilize PCPs to model this cluster property [12]. In one PCP, parent points are distributed following a homogeneous PPP with density . Around each parent point at , daughter points (denoted by ) are independently and identically distributed (i.i.d.) forming a cluster. The number of the daughter points can be constant,
, or a Poisson random variable with mean
for different purposes. Based on this spatial model, we assume that parent and daughter points represent locations of UAVs and users, respectively.Note that TCP and MCP are two widely used patterns of PCPs. The main difference between them is the distribution of daughter points in each cluster. For TCP, these points are modeled as a symmetrical normal distribution around the central parent point, with a standard deviation
. For MCP, the daughter points are uniformly distributed in a disc with radius . Therefore, the two density functions for the distance between a user at to the above UAV at can be expressed as(1)  
(2) 
where the superscript represents TCP and represents MCP. The is the unit step function, which is given by
(3) 
Regarding the altitude in the uplink phase, we assume that all users are located at the ground with an altitude and each UAV hovers above one cluster center at an altitude to serve the intracluster users.
IiB Blockage Model
The authors in [20] have proposed a tractable blockage model. The extracted theoretical expressions are derived based on a three dimensional obstacle environment, where the average height and the density of obstacles are adjustable. Therefore, this model is suitable for various practical scenarios, e.g., the suburban with lowdensity and lowaltitude obstacles, the urban with dense and high buildings, and so forth. In this model, the density of obstacles and the ratio of the obstacle area to the total area are represented by m and , respectively. The height of each obstacle is modeled as a Rayleigh distribution with a scale parameter . When the horizontal transmission distance is
, then the probability density function (PDF) of a LOS link is given by
[25](4) 
where and is the floor function. The subscript and represent transmitters and receivers, respectively^{2}^{2}2In the downlink phase, the transmitters are UAVs and the receivers are BSs, while in the uplink phase, and represent users and UAVs, respectively. We omit this explanation in the rest of this paper.. Intuitively, the PDF of a NLOS link is . As a result, the path loss law in the UAV network can be expressed as
(5) 
where is a Bernoulli random variable with parameter of success probability . The parameters and are the considered intercept and path loss exponent, respectively. Additionally, , where represents LOS links and denotes NLOS links.
IiC Directional Beamforming
We consider the typical sectorized antenna pattern, as mentioned in [26] at all transceivers to accomplish the sharp directional beamforming (see Fig. 2). The numbers of antenna elements at macro BSs, UAVs and users are , and , respectively. Four main characteristics of antennas are introduced in this pattern, which are the halfpower beamwidth in the azimuth plane , the counterpart in the elevation plane , the main beam gain , the side lobe gain (). Based on the hierarchical beam search scheme [27], the location information of devices can be obtained at the transmitter. Then, it adjusts the antenna direction toward the corresponding receiver in order to achieve the maximum beamforming gain .
Regarding one interfering transmission, the angles deviating from the boresight direction in the azimuth plane and the elevation plane are and , respectively. We assume that is uniform in the range and is uniformly distributed over the range for all interfering transmissions [21]. The directivity gain at one receiver with the interfering transmitter located at can be expressed as follows:
(6) 
where denotes the directional antenna gain. Therefore, has four patterns as shown in Table. I. Each pattern has the value with the probability , where and .
IiD Signal Model
IiD1 Downlink Phase
We assume that the number of BSs is greater than that of UAVs, namely , such that all UAVs in the downlink phase are active at the considered time slot. The typical BS at is randomly selected and it is fixed at the origin of the downlink plane. To achieve the best quality of service, the corresponding UAV, which is located at , is the closest UAV to the typical BS. Since the interfering signals are offered by the remaining UAVs, the received SINR at the typical BS can be expressed as
(7) 
where is the transmit power for each UAV, is the small scale gain for Nakagami fading channel with parameter such that is a normalized Gamma variable, and is the thermal noise.
IiD2 Uplink Phase
Since the downlink and uplink phases are analyzed in two independent planes, the corresponding UAV can belong to any of the clusters in the uplink plane. We assume that this UAV is located at the origin of the uplink plane and its location is denoted by . In one time slot, we randomly select an intracluster user at to be the transmitting user. In contrast to the downlink phase, the interference in the uplink phase is originated by both intracluster and intercluster users. Then, the received SINR at the corresponding UAV is given at the top of next page.
(8) 
In (8), is the transmit power for each user.
Iii Performance Evaluation for UAVaided Networks
In this section, we first evaluate the performance of UAVaided networks with the aid of the proposed unified framework. Then, we provide several practical scenarios which can be modeled by the proposed framework after minor adjustments.
Iiia Performance of the Downlink Phase
In this part, we analyze the coverage performance of downlink transmissions based on the distribution of distances^{3}^{3}3In the remainder of this paper, unless other specified, the distance stands for the horizontal transmission distance. and Laplace transform of interference.
IiiA1 Distance Distributions in Poisson Point Processes
Assuming that the distance realizations of UAVs in the downlink phase form a set , where and is the number of UAVs. The subscript is the parameter of ranked distances, namely . Note that the corresponding UAV is the closest node to the reference BS. Therefore, the PDF of the nearest communication distance is given by [10]
(9) 
In addition to the corresponding UAV, the rest of UAVs which are located further than are interfering transmitters. Due to the Slivnyak’s theorem, the density of interfering UAVs is still in the area , where represents an annulus with the inner radius and outer radius [28]. Moreover, each interferer is distributed independently such that the index can be dropped from . Therefore, the PDF of the interfering distance can be derived via the probability generating function of PPP [29].
IiiA2 Laplace Transform of Interference in Poisson Point Processes
For simplifying the notation, we define the Laplace transform of interference as follows:
(10) 
where is the expectation function, is the transform parameter, and represents the received power of interference, which is given by
(11) 
Next, we provide closedform expressions for the Laplace transform of interference in the following lemma and corollary.
Lemma 1.
Since the corresponding UAV at is the closest node to the typical BS with a communication distance , all interfering UAVs are located in the area , namely . A closedform expression for the Laplace transform of interference in the downlink phase conditional on is given by
(12) 
where
(13) 
and
(14) 
when , can be expressed as
(15) 
with being Gauss hypergeometric function.
For most of mmWave frequencies, the path loss exponent of LOS transmissions equals to two [30, 31, 32]. Therefore, when , is changed to
(16) 
and .
Proof.
See Appendix A. ∎
IiiA3 Coverage Probability for the Downlink Phase
In order to insure the quality of the downlink communicating service, a targeted rate is predecided in most of the cases. Then, the corresponding desired SINR threshold is given by , where is the bandwidth for each downlink resource block. Therefore, the coverage probability is defined as the proportion of received SINR that exceeds , which can be expressed as
(17) 
where denotes probability. With the aid of the aforementioned distance distribution and Laplace transform of interference, we present the coverage probability in the following theorem.
Theorem 1.
When choosing the nearest UAV at as the corresponding UAV, then the serving communication distance is . Accordingly, a closedform expression for the coverage probability at the typical BS is given by
(18) 
where
(19) 
and . The set of nonnegative integers , and is a GaussChebyshev node, which equals . When , the equality holds. We are able to change the value of to balance the complexity and efficiency [10].
Proof.
See Appendix B. ∎
Remark 1.
Since we randomly select the typical BS, the coverage provability in the downlink phase is independent of the distribution of the macro BSs .
Remark 2.
Although the range of is infinite due to , we are able to choose the first three values, namely , for enhancing the computation efficiency. The reason is that when the communication distance increases, the received power of both signals and interference decreases, which have a negligible impact on the received SINR.
Remark 3.
Note that the point process of UAVs is same with that of macro BSs. When considering the uplink transmissions from the BSs to UAVs in the proposed downlink phase, the expression of coverage probability in Theorem 1 is also valid.
IiiB Performance of the Uplink Phase
Due to the cluster property of PCP, the analysis of the uplink phase is more challenging compared with the downlink phase. Similarly, we first present the distance distribution and Laplace transform of the interference in the uplink phase. Then, the coverage performance at the corresponding UAV is analyzed.
IiiB1 Distance Distribution in the Corresponding Cluster
The communication distances in the same cluster with the corresponding UAV form a group and , where is the number of users in the considered cluster. By utilizing the fair selection strategy, the subscript can be removed from and . Regarding TCP, the PDF of the distance between one user in the corresponding cluster to the corresponding UAV is given by [13]
(20) 
Then, for MCP, the PDF of that distance can be expressed as [10]
(21) 
IiiB2 Distance Distributions in Other Clusters
In other clusters, the distances between the users with a cluster center and the corresponding UAV compose a set and . As all users are i.i.d, we are able to drop the subscript from and . For TCP, the PDF of this distance is conditional on the distance , which is given by [13]
(22) 
where is the first kind for the modified Bessel function with order zero. Then, for MCP, the PDF of this distance is [33]
(23) 
Remark 4.
In fact, eqs. (22) and (IIIB2) are the PDF of the distance between a daughter point to the origin. The origin can be anywhere on the plane. For example, (22) and (IIIB2) are also valid for D2D communications, where one daughter point is fixed at the origin instead of one parent point as mentioned in the proposed system.
IiiB3 Laplace Transform of Interference in Poisson Cluster Processes
There exist two kinds of interference in the uplink phase. One is intracluster interference and the other is intercluster interference . As a result, Laplace transform of interference can be defined as
(24) 
where
(25)  
(26) 
In several special networks, e.g., D2D communications and nonorthogonal multiple access (NOMA) networks, the number of users in one cluster should be fixed as they are paired. For other networks, the number of users should be random across different clusters to enhance the generality. Therefore, we consider two cases in this paper: 1) Fixed Case: the number of users in each cluster is fixed as ; and 2) Random Case: the number of users in each cluster follows Poisson distribution with the mean .
Lemma 2.
When the distance between an intracluster interferer and the corresponding UAV is , then the Laplace transform of the interference in the fixed case is given by
(27) 
where
(28) 
and , , and .
On the other hand, if we consider the random case, the Laplace transform of the intracluster interference can be expressed as
(29) 
Proof.
See Appendix C. ∎
In terms of the intercluster interference, the Laplace transform of this interference can be expressed in the following lemma and corollary.
Lemma 3.
When the distance between an intercluster interferer and the corresponding UAV is , for the fixed case, the Laplace transform of intercluster interference can be expressed as
(30) 
and for the random case, the corresponding Laplace transform of intercluster interference is changed to
(31) 
where
(32) 
Proof.
See Appendix D. ∎
Remark 5.
Note that when is large, , . Since and , if the number of users in each cluster is large and , the difference between Laplace transform of interference in the fixed case and the random case is negligible.
IiiB4 Coverage Probability for the Uplink Phase
We predecide a targeted rate for the uplink phase such that the corresponding SINR threshold is given by and is the bandwidth for each uplink resource block. Then, the coverage probability for the uplink phase is defined as
(33) 
Note that in the corresponding cluster, the transmitting user is randomly selected. The communication distance between the corresponding UAV and the transmitting user is . Based on eqs. (8) and (33), the analytical expression for the uplink coverage probability can be expressed in the following theorem and corollary.
Theorem 2.
When the SINR threshold for the uplink phase is , then the coverage probability at the corresponding UAV is given by
(34) 
where
(35) 
and , .
Proof.
In terms of coverage probability, there exist two differences between the downlink phase and uplink phase. One is the distribution of communication distance and the other is the received interference power. Based on Theorem 1, we replace (9) by (20) and (21). Additionally, we utilize instead of . After that, we obtain Theorem 2. ∎
By analyzing Theorem 1 and Theorem 2, we are able to abstract a general expression of the coverage probability with an SINR threshold for all point processes, including PPP and PCP.
Proposition 1.
Assuming that Nakagami fading channel is utilized in the considered wireless networks, then the coverage probability for both PPP and PCP can be expressed as
(36) 
where the index represents the general case, is the signaltonoiseratio (SNR) coverage probability part, is the Laplace transform of interference, and is the PDF of communication distance between the serving transmitter and the served receiver.
Remark 6.
When the Nakagami parameter , namely Rayleigh fading channel, then the equality holds in (1).
IiiC System Coverage Probability
Without loss of generality, we assume that the probability of the message being successfully transmitted from the uplink to the downlink phase is . In the proposed system, we assume that the corresponding UAV is able to successfully flies from the uplink plane to the downlink plane, such that the messages can be guaranteed to reach the downlink phase, namely . Therefore, we are capable to deduce the system coverage probability for the proposed UAV networks.
Proposition 2.
Assuming that the SINR thresholds for the downlink phase and uplink phase are and , respectively, then the system coverage probability is given by
(37) 
IiiD Practical Applications
Next, we provide three typical UAV applications: UAVaided ubiquitous coverage, UAVaided information dissemination and data collection, and UAVaided relaying [2]. The proposed unified framework is able to model them after minor modifications. As it can be seen from Fig. 3, for the ubiquitous coverage case, since the corresponding UAV aims to offload the data traffic of the central BS, all transceivers should be modeled in one plane. For the information dissemination and data collection case, the proposed framework can be utilized directly. For the relaying case, since two corresponding UAVs act as two relays, the success probability is mainly decided by their communication conditions. In order to make the analysis more complete, we introduce a simplified cooperative transmission between two relays in the following part, which can be extended to other complicated scenarios. Moreover, the time slot for this transmission is located between the uplink and downlink phases. Special modifications for modeling these three applications are summarized in Table II.
Applications  Uplink Phase  Downlink Phase  Cooperative Transmission  

Ubiquitous coverage 

None,  

Keep same  Keep same  Flight process,  
Relaying  Keep same  Keep same 

IiiD1 Simplified Cooperative Transmission
In most of the cases, two relaying UAVs are connected through the air. Due to limited obstacles in the air, these links can be regarded as LOS transmissions. Assuming that the distance between the two relays is fixed, and for simplicity, the cochannel interferers can be ignored, then the received SNR at the receiver is given by
(38) 
In this cooperative transmission, represents the corresponding coverage probability, which can be calculated in the following lemma.
Lemma 4.
When the considered cooperative transmission distance is fixed as , the closedform coverage probability with an SNR threshold can be expressed as follows:
(39) 
where and is the upper incomplete gamma function.
Proof.
By using the cumulative distribution function (CDF) of the Gamma distribution and following similar proof as in
Theorem 1, we obtain Lemma 4. ∎IiiD2 Coverage Probability for the Relaying Case
The system coverage probability for the relaying case can be represented by Proposition 2 as well. The only difference is the probability . Therefore, the system coverage probability is expressed as follows:
Proposition 3.
Assuming that the SINR thresholds for the downlink phase, uplink phase, and cooperative transmission are , , and respectively, the system coverage probability for the relaying case with a fixed cooperative transmission distance is given by
(40) 
Remark 7.
Note that the LOS transmissions outperform NLOS transmissions in mmWave communications. When the altitude of UAVs increases, the PDF in (IIB) enlarges, which enhances the system coverage probability . However, the high altitude of UAVs also increases the communication distance such that the path loss in (IIB) raises. Therefore, can be maximized by selecting optimal altitude of UAVs.
Comments
There are no comments yet.