This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Access Points in the Air: Modeling and Optimization of Fixed-Wing UAV Network

Junyu Liu, Min Sheng, Ruiling Lyu, Yan Shi, Jiandong Li
State Key Laboratory of ISN, Xidian University, Xi’an, Shaanxi, 710071, China
Abstract

Fixed-wing unmanned aerial vehicles (UAVs) are of great potential to serve as aerial access points (APs) owing to better aerodynamic performance and longer flight endurance. However, the inherent hovering feature of fixed-wing UAVs may result in discontinuity of connections and frequent handover of ground users (GUs). In this work, we model and evaluate the performance of a fixed-wing UAV network, where UAV APs provide coverage to GUs with millimeter wave backhaul. Firstly, it reveals that network spatial throughput (ST) is independent of the hover radius under real-time closest-UAV association, while linearly decreases with the hover radius if GUs are associated with the UAVs, whose hover center is the closest. Secondly, network ST is shown to be greatly degraded with the over-deployment of UAV APs due to the growing air-to-ground interference under excessive overlap of UAV cells. Finally, aiming to alleviate the interference, a projection area equivalence (PAE) rule is designed to tune the UAV beamwidth. Especially, network ST can be sustainably increased with growing UAV density and independent of UAV flight altitude if UAV beamwidth inversely grows with the square of UAV density under PAE.

I Introduction

Benefiting from rapid event response and flexible deployment, unmanned aerial vehicle (UAV) networks have experienced explosive development and evolution in recent years [1, 2]. Among the applications of great potential, low-altitude UAVs are capable of serving as the access points (APs) to provide service to ground users (GUs) in short-distance line-of-sight wireless channels [3, 4]. Moreover, the deployment of UAV APs is flexible to provide instant coverage in a variety of scenarios including public safety, dense crowds, internet of things (IoT) applications and emergency scenarios. Thanks to better aerodynamic performance, greater payload versatility and longer flight endurance, fixed-wing UAVs are promising in serving as the APs in the air [5]. Moreover, millimeter wave (mmWave) by nature can be used to convey wireless backhaul for UAV APs, since there are few obstacles for the air-to-ground channel. In particular, it was reported that up to 10 Gbps peak rates could be reached using multi-user multiple-input-multiple-output (MU-MIMO) mmWave for long-range transmissions [6].

Despite the great potential, fixed-wing UAVs have to hover in the air when providing coverage to a given ground area. With the increase of hover radius, it is difficult for fixed-wing UAVs to provide steady and continuous service to GUs. Moreover, mmWave beams are basically oriented to the hovering UAVs to guarantee the backhaul capacity. Therefore, increasing the UAV hover radius may potentially enlarge the mmWave beamwidth, thereby degrading the mainlobe gain and backhaul capacity. Hence, it is crucial to investigate the impact of UAV hover radius and how to effectively enhance the performance of the fixed-wing UAV network.

I-A Related Work

The research on UAV networks, where UAVs serve as APs, has attracting considerable attention from both academia and industry [7, 8, 9, 10, 11, 4]. The UAV trajectory optimization problem was considered in [8], where UAVs provided service to cell edge users to offload the data traffic of terrestrial cellular network. Specifically, UAV trajectory and user scheduling were alternately optimized to maximize the sum rate of UAV-served edge users with the rate constraint of all the users. In [9], a distributed UAV deployment algorithm was proposed to minimize the distance from UAV APs to GUs so as to improve the user coverage. Moreover, the UAV AP deployment problem has been extended into heterogeneous networks [10]. Especially, a latency-aware approach has been proposed, which jointly optimizes the location and association coverage of UAVs, to minimize the total average latency of GUs considering the existing terrestrial base stations (TBSs). If frequency resources are shared by UAV APs and TBSs, potential cross-layer interference will be generated, which may degrade the performance of the UAV network. The performance of a UAV integrated terrestrial cellular network is evaluated in [11, 4]. Especially, the impact of UAV deployment density on the performance of the coexisting system was captured in [11] and an effective interference avoidance scheme was proposed to mitigate the air-to-ground interference and significantly improve the spectrum efficiency of the coexisting system in [4].

Due to flexible deployment and low expenditure, rotary-wing UAVs are applied in most of the available research [7, 8, 9, 10, 11, 4], where rotary-wing UAVs were assumed to move according to the preset trajectory and steadily stay in a given position. Accordingly, rotary-wing UAV systems have been designed and developed by a number of companies and operators [12, 13]. However, since most of off-the-shelf rotary-wing UAVs are battery driven, they are energy-inefficient and excessive power is required especially when they are taking off, climbing and changing the flight attitude. Therefore, the flight endurance of the rotary-wing UAVs is significantly limited (dozens of minutes for battery-driven ones). Moreover, limited wireless backhaul will bottleneck the performance of the UAV system supposing that high data rates are required by the connected GUs. To address the issues of limited flight endurance and backhaul, AT&A was reported to develop a tethered UAV system, named Flying COW [13]. Particularly, each rotary-wing UAV, which was connected to the ground by a thin tether, is capable of providing the long-term evolution (LTE) coverage to GUs. The tether was used for conveying highly secure backhaul via fiber and supplying power, which allows for longer flight time.

Although Flying COW is designed to ultimately provide coverage to an area up to 40 square miles, the tether connection will inevitably limit the deployment and mobility of UAVs. On the other hand, researchers and engineers are increasingly focusing on the development of the fixed-wing UAVs [14, 15, 5, 16]. Compared with the rotary-wing UAVs, the deployment of fixed-wing UAVs is more flexible. For instance, fixed-wing UAVs could offer better aerodynamic performance, which makes them well suited for the application in higher flight altitude, greater payload and longer ranges and flight endurance. In [16], an all-weather fixed-wing emergency communication system was developed for the quick network recovery in the emergency communication scenarios. According to the aerodynamic principle, nevertheless, fixed-wing UAVs cannot steadily stay in a given position and have to hover in the air. If the hover radius is large, rotary-wing UAVs would fail to provide continuous service to GUs. Moreover, although sufficient backhaul can be provided through mmWave link [17, 18], the high mobility of the hovering fixed-wing UAVs may seriously deteriorate mmWave link capacity since the mmWave beam is highly directional. Worse still, the acquisition of real-time perfect channel state information (CSI) in the highly dynamic network is challenging and remains to be an open problem [19]. Especially, the acquired CSI would be easily outdated at the moment of decision due to the on-the-move feature of fixed-wing UAVs. Therefore, the performance and deployment of fixed-wing UAV network remains to be further investigated, which motivates this work.

I-B Contribution and Outcome

In this paper, we model a downlink UAV network, where fixed-wing UAV APs provide service to GUs and mmWave is applied to provide wireless backhaul. In particular, the impact of key parameters, including UAV flight altitude, deployment density, hover radius, backhaul limitation and user association rules, etc., on the performance of UAV network has been evaluated. On this basis, we further investigate how to effectively improve the UAV network performance through adjusting UAV beamwidth. The main conclusions of this work are summarized as follows.

  • Real-time VS semi-real-time user association rules. We evaluate the performance of two typical user association rules, namely, real-time user association (RTNA), where GUs always connect the closest UAVs, and semi-real-time user association (Semi-RTNA), where GUs connect to the UAVs whose hover centers are the closest, in terms of network spatial throughput (ST). It is shown that network ST is greatly degraded by increasing the UAV hover radius under Semi-RTNA. The reason is that desired signal power is likely to be reduced if the associated UAVs hover apart under large hover radius. Worse still, the increase of hover radius would notably increase the mmWave mainlobe beamwidth, thereby reducing the mainlobe gain and backhaul capacity. On the contrary, the impact of hover radius on the performance of RTNA is shown to be minor even though more handover overhead is introduced. Especially, if ignoring the backhaul limitation, network ST is shown to be independent of the hover radius under RTNA.

  • Optimization of UAV beamwidth. The impact of UAV beamwidth on network ST is further investigated. Particularly, reducing the UAV beamwidth first increases (due to alleviating overlap-cell interference) and then decreases (due to limiting UAV coverage) network ST. Moreover, the optimal beamwidth is shown to be dependent on the activated UAV density λa\lambda_{\mathrm{a}}. Accordingly, we propose a projection area equivalence policy to optimize the UAV beamwidth. In particular, network ST could be sustainably increased with UAV density if UAV half-beamwidth inversely grows with 1λa\frac{1}{\sqrt{\lambda_{\mathrm{a}}}}. More importantly, the optimized network ST is proved to be independent of the UAV flight altitude in the backhaul-unlimited case. In other words, aerial spatial resources can be fully exploited by the proposed policy.

We organize the remaining parts of this paper as follows. System model is given in Section II followed by the performance analysis of two typical user association rules in UAV network with respect to user coverage probability and network ST in Section III. In Section IV, we further investigate the impact of UAV antenna beamwidth on the performance of directional-antenna UAV network and optimize network ST through adjusting UAV beamwidth. Finally, conclusions are drawn in Section V. The main parameter notations used in the paper are summarized in Table I.

Table I: Summary of Notations
Symbol Meaning Symbol Meaning
UAVi\mathrm{UAV}_{i}, GUj\mathrm{GU}_{j} iith TBS, jjth UAV Φ\Phi UAV half-beamwidth
ΠUAV\Pi_{\mathrm{UAV}}, ΠGU\Pi_{\mathrm{GU}} UAV and GU sets G(ϕ,φ)G\left(\phi,\varphi\right) antenna mainlobe gain
Π~UAV\mathrm{\tilde{\Pi}_{UAV}} activated UAV set RpR_{\mathrm{p}} projection radius
λ\lambda, λGU\lambda_{\mathrm{GU}} UAV and GU densities ppp_{\mathrm{p}} projection probability
λa\lambda_{\mathrm{a}} activated UAV density α\alpha pathloss exponent
qaq_{\mathrm{a}} UAV activated probability RmR_{\mathrm{m}} backhaul coverage radius
η\eta η=λ/λGU\eta=\lambda/\lambda_{\mathrm{GU}} CtC_{\mathrm{t}} normalized backhaul capacity
μ\mu μ\mu=3.5 CbC_{\mathrm{b}} backhaul capacity for each UAV
RhR_{\mathrm{h}} UAV hover radius ϑ\vartheta mmWave mainlobe beamwidth
θ\theta UAV hover angle Gm(ϑ)G_{\mathrm{m}}\left(\vartheta\right) mmWave mainlobe gain
hUAVh_{\mathrm{UAV}}, hGUh_{\mathrm{GU}} UAV and GU altitudes ε\varepsilon mmWave orientation error
hLh_{\mathrm{L}} lower bound of hUAVh_{\mathrm{UAV}} ε¯\bar{\varepsilon} mean of ε\varepsilon
hUh_{\mathrm{U}} upper bound of hUAVh_{\mathrm{UAV}} F|ε|(ϑ)F_{\left|\varepsilon\right|}\left(\vartheta\right) mmWave orientation error
Δh\Delta h Δh=hUAVhGU\Delta h=h_{\mathrm{UAV}}-h_{\mathrm{GU}} r0r_{0} 2D distance from UAV0\mathrm{UAV}_{0} to GU0\mathrm{GU}_{0}
h~L\tilde{h}_{\mathrm{L}} h~L=hLhGU\tilde{h}_{\mathrm{L}}=h_{\mathrm{L}}-h_{\mathrm{GU}} CC projection scaling parameter
h~U\tilde{h}_{\mathrm{U}} h~U=hUhGU\tilde{h}_{\mathrm{U}}=h_{\mathrm{U}}-h_{\mathrm{GU}} PP UAV transmit power

II System Model

II-A Network Model

Consider a downlink UAV network, where fixed-wing UAV APs provide service to GUs. For the two-dimension (2D) locations, UAVs and GUs are distributed as two independent homogeneous Poisson Point Processes (PPPs) ΠUAV={UAVi}\Pi_{\mathrm{UAV}}=\left\{\mathrm{UAV}_{i}\right\} (i)\left(i\in\mathbb{N}\right) with density λ\lambda and ΠGU={GUj}\Pi_{\mathrm{GU}}=\left\{\mathrm{GU}_{j}\right\} (j)\left(j\in\mathbb{N}\right) with density λGU\lambda_{\mathrm{GU}}, respectively. For practical concern, we assume that fixed-wing UAVs hover in the air with identical hover radius RhR_{\mathrm{h}} and different flight altitude hUAVh_{\mathrm{UAV}}, which follows uniform distribution hUAV𝒰(hL,hU)h_{\mathrm{UAV}}\sim\mathcal{U}\left(h_{\mathrm{L}},h_{\mathrm{U}}\right). Supposing that GUs are of identical height hGUh_{\mathrm{GU}}, the vertical distance Δh(=hUAVhGU)\Delta h\left(=h_{\mathrm{UAV}}-h_{\mathrm{GU}}\right) from UAVs to GUs follows uniform distribution Δh𝒰(h~L,h~U)\Delta h\sim\mathcal{U}\left(\tilde{h}_{\mathrm{L}},\tilde{h}_{\mathrm{U}}\right), where h~L=hLhGU\tilde{h}_{\mathrm{L}}=h_{\mathrm{L}}-h_{\mathrm{GU}} and h~U=hUhGU\tilde{h}_{\mathrm{U}}=h_{\mathrm{U}}-h_{\mathrm{GU}}. It is assumed that multiple GUs can simultaneously connect to one UAV AP for service. For scheduling fairness, each UAV AP randomly and independently serves the connected GUs in a time-division manner. Therefore, the GUs in one UAV cell have the equal chance to be served. To improve frequency reuse, spectrum is reused by different UAV cells. As a result, the neighboring UAV APs may potentially generate inter-cell interference to the intended GU. If not properly handled, the inter-cell interference may significantly degrade the performance of the downlink UAV network. Besides, saturated data model is adopted such that each GU always has data request from the connected UAVs.

Refer to caption
Figure 1: Illustration of a fixed-wing UAV network. RTNA and Semi-RTNA are used for user association.

II-B User Association Model

We adopt the following two user association rules to balance the UAV coverage performance and handover overhead.

  • Real-time nearest association (RTNA). One GU is associated with the UAV, which is horizontally closest to the GU111The flight altitude of UAV APs may be varying due to air buoyancy and attitude control to ensure that the hovering UAVs could cover a given area. Therefore, each GU is assumed to connect to horizontally closest UAV for stronger average received power and better coverage.. It is shown in Fig. 1 that GU1\mathrm{GU}_{1} is located closer to the hover center of UAV1\mathrm{UAV}_{1} compared to that of UAV5\mathrm{UAV}_{5}. However, the hovering UAV5\mathrm{UAV}_{5} is closer to GU1\mathrm{GU}_{1} in the current time slot. Therefore, GU1\mathrm{GU}_{1} is associated with UAV5\mathrm{UAV}_{5} under RTNA. In the coming time slots, when UAV1\mathrm{UAV}_{1} hovers closer, GU1\mathrm{GU}_{1} will be handover to UAV1\mathrm{UAV}_{1}.

  • Semi-real-time nearest association (Semi-RTNA). One GU is associated with the UAV, the hover center of which is horizontally closest to the GU. It is shown in Fig. 1 that, even if GU2\mathrm{GU}_{2} is closer to UAV2\mathrm{UAV}_{2} in the current time slot, GU2\mathrm{GU}_{2} is associated with UAV3\mathrm{UAV}_{3} under Semi-RTNA to reduce the handover overhead.

It is assumed that real-time location information of UAVs could be acquired by GUs for performing RTNA and Semi-RTNA. Accordingly, GUs could calculate the distance to the UAVs, which could provide coverage, to determine the UAV to connect to. Then, the handover of user association is initiated by the GU, which sends the handover request to the associated UAV AP and the UAV AP to connect to. Note that the handover frequency is dependent on the location of GUs, number of fixed-wing UAV APs near the GUs, hover radius and velocity of UAV APs. If not properly handled, the frequent handover will result in considerable handover overhead, thereby degrading the performance of the UAV network.

II-C Antenna and Channel Model

It is shown in Fig. 1 that each UAV with constant transmit power PP is equipped with one directional antenna. The azimuth and elevation beamwidths are 2Φ2\Phi with Φ(0,π2)\Phi\in\left(0,\frac{\pi}{2}\right). Denoting ϕ\phi and φ\varphi as the azimuth and elevation angles, respectively, the antenna gain in the direction (ϕ,φ)\left(\phi,\varphi\right) (corresponding to the antenna mainlobe) is modeled as [20]

G(ϕ,φ)={G0/Φ2,ΦϕΦ,ΦφΦ,0,otherwise,G\left(\phi,\varphi\right)=\begin{cases}G_{\mathrm{0}}/\Phi^{2},&-\Phi\leq\phi\leq\Phi,-\Phi\leq\varphi\leq\Phi,\\ 0,&\mathrm{otherwise},\end{cases} (1)

where G0=2.2846G_{0}=2.2846. Accordingly, the projection of each UAV on the ground is a disk region with radius Rp=ΔhtanΦR_{\mathrm{p}}=\Delta h\tan\Phi. If one GU is out of the projection area of any UAV, no service will be provided. Note that, although increasing Φ\Phi could enhance the UAV coverage, it will result in greater overlap of the projection areas, which may introduce excessive inter-cell interference.

The channel from UAVs to GUs is assumed to consist of pathloss dαd^{-\alpha} and Rayleigh fading h𝒞𝒩(0,1)h\sim\mathcal{CN}\left(0,1\right) [21, 22], where dd and α\alpha denote the distance from one UAV to its associated GU and pathloss exponent, respectively.

II-D Backhaul Model

We consider that mmWave is used to convey wireless backhaul for fixed-wing UAVs owing to 1) higher frequency bands, which is non-overlapping with those in UAV-GU channel, and 2) greater bandwidth and backhaul capacity [6, 23]. For instance, it was tested in [24] that approximately 10Gbps transmission rate with <20% outage could be provided through 200m mmWave transmissions over Ka band (28GHz), V band (60GHz) and E band (73GHz) even if 3 blockages exist between mmWave transmitters and receivers. Since the number of ground gateways is limited, we assume that mmWave backhaul capacity is shared by the UAVs within the region of radius RmR_{\mathrm{m}} [23, 25, 26]. Specifically, each UAV occupies equal proportion of backhaul resources. According to [25, 26], the backhaul capacity for each UAV is given by

Cb=Gm(ϑ)F|ε|(ϑ)CtλaπRm2,bits/(sHz)C_{\mathrm{b}}=\frac{G_{\mathrm{m}}\left(\vartheta\right)F_{\left|\varepsilon\right|}\left(\vartheta\right)C_{\mathrm{t}}}{\lambda_{\mathrm{a}}\pi R_{\mathrm{m}}^{2}},\>\mathrm{bits}/\left(\mathrm{s\cdot Hz}\right) (2)

where CtC_{\mathrm{t}} denotes the normalized backhaul capacity in bits/(sHz)\mathrm{bits}/\left(\mathrm{s\cdot Hz}\right) and λa\lambda_{\mathrm{a}} denotes the density of activated UAVs222Note that UAVs would keep inactivated with no backhaul requirement if no GUs are connected.. Therefore, the denominator λaπRm2\lambda_{\mathrm{a}}\pi R_{\mathrm{m}}^{2} denotes the number of activated UAVs within the region.

In (2), Gm(ϑ)=2πϑG_{\mathrm{m}}\left(\vartheta\right)=\frac{2\pi}{\vartheta} denotes the mmWave mainlobe gain and F|ε|(ϑ)F_{\left|\varepsilon\right|}\left(\vartheta\right) denotes the mmWave orientation effective function [23], which is dependent on mmWave mainlobe beamwidth ϑ\vartheta and orientation error ε\varepsilon. Given |ε|ϑ2\left|\varepsilon\right|\leq\frac{\vartheta}{2}, backhaul could be effectively conveyed through mmWave. Otherwise, UAV fails to obtain backhaul due to the orientation error. Supposing that the mmWave mainlobe could cover the UAV hover region, we have ϑ=2arctan(RhhUAV)\vartheta=2\arctan\left(\frac{R_{\mathrm{h}}}{h_{\mathrm{UAV}}}\right) under large hUAVh_{\mathrm{UAV}}. Accordingly, increasing the UAV hover radius may alleviate the impact of orientation error, while reduces the mmWave mainlobe gain.

II-E Performance metric

In this work, we use coverage probability (CP) and network ST to evaluate the performance of the UAV network. In particular, one GU is in coverage when 1) the GU is in the projection area of UAVs and 2) the data transmission from the associated UAV to the GU is successful. Without loss of generality, we denote the link, which consists of UAV0\mathrm{UAV}_{0} and GU0\mathrm{GU}_{0}, as the typical link. Accordingly, the CP of the typical link is defined by

𝖢𝖯\displaystyle\mathsf{CP} ={UAV0GU0<Rp,𝖲𝖨𝖱>τ}\displaystyle=\mathbb{P}\left\{\left\|\mathrm{UAV}_{0}-\mathrm{GU}_{0}\right\|<R_{\mathrm{p}},\mathsf{SIR}>\tau\right\}
=pp{𝖲𝖨𝖱>τ|Ep}.\displaystyle=p_{\mathrm{p}}\mathbb{P}\left\{\mathsf{SIR}>\tau\left|E_{\mathrm{p}}\right.\right\}. (3)

In (3), EpE_{\mathrm{p}} denotes the event that GU0\mathrm{GU}_{0} is within the projection area of UAV0\mathrm{UAV}_{0} and we denote the corresponding projection probability as pp={UAV0GU0<Rp}p_{\mathrm{p}}=\mathbb{P}\left\{\left\|\mathrm{UAV}_{0}-\mathrm{GU}_{0}\right\|<R_{\mathrm{p}}\right\}. If EpE_{\mathrm{p}} occurs, {𝖲𝖨𝖱>τ|Ep}\mathbb{P}\left\{\mathsf{SIR}>\tau\left|E_{\mathrm{p}}\right.\right\} denotes the transmission success probability, where 𝖲𝖨𝖱\mathsf{SIR} denotes the signal-to-interference ratio (SIR) at GU0\mathrm{GU}_{0} and τ\tau is the SIR threshold.

Based on (3), network ST is defined by

𝖲𝖳\displaystyle\mathsf{ST} =λa𝖢𝖯min{log2(1+τ),Cb}[bits/(sHzm2)].\displaystyle=\lambda_{\mathrm{a}}\mathsf{CP}\min\left\{\log_{2}\left(1+\tau\right),C_{\mathrm{b}}\right\}\>\left[\mathrm{bits/(s\cdot Hz\cdot m^{2})}\right]. (4)

With the constraint of CbC_{\mathrm{b}}, 𝖲𝖳\mathsf{ST} captures how many bits could be successfully conveyed over unit time, frequency and area by the UAV network under limited backhaul.

III User Association in UAV Network

In this section, we evaluate the performance of the UAV network under two user association rules, i.e., Semi-RTNA and RTNA. To capture the impact of hover radius, we consider that the half-beamwidth Φ\Phi in (1) equals π2\frac{\pi}{2}. In consequence, the radius of the projection area of each UAV approaches infinity, i.e., RpR_{\mathrm{p}}\rightarrow\infty, and UAV antenna gain degenerates into 1. Besides, GUs are always in the projection area of the UAVs, i.e., pp=1p_{\mathrm{p}}=1 in (3).

Following the definitions of CP and network ST in Section II-E, the key to the analysis is to calculate the distribution of SIR at GU0\mathrm{GU}_{0}. However, the SIR distributions under Semi-RTNA and RTNA rules are different, which will be discussed in the following.

III-A Semi-RTNA Rule

With Semi-RTNA, the SIR at GU0\mathrm{GU}_{0} can be expressed as

𝖲𝖨𝖱GU0SR=\displaystyle\mathsf{SIR}_{\mathrm{GU}_{0}}^{\mathrm{SR}}= PHUAV0d0αISR,\displaystyle\frac{PH_{\mathrm{UAV}_{0}}d_{0}^{-\alpha}}{I^{\mathrm{SR}}}, (5)

where d0d_{0} denotes the distance from UAV0\mathrm{UAV}_{0} to GU0\mathrm{GU}_{0} and ISR=UAViΠ~UAV\UAV0PHUAVidiαI^{\mathrm{SR}}=\underset{\tiny{\mathrm{UAV}_{i}\in\mathrm{\tilde{\Pi}_{UAV}}\backslash\mathrm{UAV}_{0}}}{\sum}PH_{\mathrm{UAV}_{i}}d_{i}^{-\alpha} denotes the interference stemming from other activated UAV cells. HUAViH_{\mathrm{UAV}_{i}} denotes the channel power gain due to small-scale fading and Π~UAV\mathrm{\tilde{\Pi}_{UAV}} denotes the set of activated UAVs with density λa=qaλ\lambda_{\mathrm{a}}=q_{\mathrm{a}}\lambda. Due to UAV association, the UAV activation probability is given by [27]

qa\displaystyle q_{\mathrm{a}} =1(1+(μη)1)μ,\displaystyle=1-\left(1+\left(\mu\eta\right)^{-1}\right)^{-\mu}, (6)

where η=λ/λGU\eta=\lambda/\lambda_{\mathrm{GU}} denotes the ratio of UAV density to GU density and μ\mu=3.5. Moreover, we suppose that the mmWave orientation error ε\varepsilon in (2) follows exponential distribution[23, 28]. Therefore, the mmWave effective function is given by the following lemma333Note that the results in Lemma 1 can be feasibly extended to the cases, where other distributions are applied to model the orientation error..

Lemma 1.

Supposing that the absolute mmWave beam orientation error |ε|\left|\varepsilon\right| follows exponential distribution truncated to [0,π]\left[0,\pi\right], the mmWave effective function F|ε|(ϑ)F_{\left|\varepsilon\right|}\left(\vartheta\right) is given by

F|ε|(ϑ)\displaystyle F_{\left|\varepsilon\right|}\left(\vartheta\right) =1exp(ϑ2ε¯)1exp(πε¯),ϑ(0,π2)\displaystyle=\frac{1-\mathrm{exp}\left(-\frac{\vartheta}{2\bar{\varepsilon}}\right)}{1-\mathrm{exp}\left(-\frac{\pi}{\bar{\varepsilon}}\right)},\>\vartheta\in\left(0,\frac{\pi}{2}\right) (7)

where ε¯\bar{\varepsilon} denotes the mean of ε\varepsilon (prior to the truncation).

Proof: Please refer to Appendix A-A.∎

On average, the mmWave effective function F|ε|(ϑ)F_{\left|\varepsilon\right|}\left(\vartheta\right) in (7) could capture the impact of mean orientation error ε¯\bar{\varepsilon} and mmWave mainlobe beamwidth ϑ\vartheta on the performance of mmWave transmission. Specifically, F|ε|(ϑ)F_{\left|\varepsilon\right|}\left(\vartheta\right) is shown to be reduced by either the increase of ε¯\bar{\varepsilon} or the decrease of ϑ\vartheta, which complies with intuition. According to Lemma 1, it can be shown that Gm(ϑ)F|ε|(ϑ)G_{\mathrm{m}}\left(\vartheta\right)F_{\left|\varepsilon\right|}\left(\vartheta\right) in (2) is a decreasing function of the mmWave mainlobe beamwidth ϑ\vartheta. Since ϑ=2arctan(RhhUAV)\vartheta=2\arctan\left(\frac{R_{\mathrm{h}}}{h_{\mathrm{UAV}}}\right), increasing the UAV hover radius would potentially degrade the mmWave backhaul capacity.

Aided by (5), (6) and Lemma 1, CP and network ST under Semi-RTNA could be obtained in the following proposition.

Proposition 1.

Each UAV equipped with Φ=π2\Phi=\frac{\pi}{2} directional antenna, network ST of the fixed-wing UAV network under Semi-RTNA rule is given by

𝖲𝖳π/2SR\displaystyle\mathsf{ST}_{\pi/2}^{\mathrm{SR}} =λa𝖢𝖯π/2SRmin{log2(1+τ),Cb}.\displaystyle=\lambda_{\mathrm{a}}\mathsf{CP}_{\pi/2}^{\mathrm{SR}}\min\left\{\log_{2}\left(1+\tau\right),C_{\mathrm{b}}\right\}. (8)

In (8), 𝖢𝖯π/2SR\mathsf{CP}_{\pi/2}^{\mathrm{SR}} is given by

𝖢𝖯π/2SR\displaystyle\mathsf{CP}_{\pi/2}^{\mathrm{SR}} =𝔼r0,θ,Δh[exp(2πλaτd0αl2αω1(α,τd0αlα)α2)],\displaystyle=\mathbb{E}_{r_{0},\theta,\Delta h}\left[\exp\left(-2\pi\lambda_{\mathrm{a}}\frac{\tau d_{0}^{\alpha}l^{2-\alpha}\omega_{1}\left(\alpha,\tau d_{0}^{\alpha}l^{-\alpha}\right)}{\alpha-2}\right)\right], (9)

where θ\theta denotes the hover angle, d0=r02+Rh22r0Rhcosθ+Δh2d_{0}=\sqrt{r_{0}^{2}+R_{\mathrm{h}}^{2}-2r_{0}R_{\mathrm{h}}\cos\theta+\Delta h^{2}} and

l={Δh,r0Rh(r0Rh)2+Δh2.r0Rhl=\begin{cases}\Delta h,&r_{0}\leq R_{\mathrm{h}}\\ \sqrt{\left(r_{0}-R_{\mathrm{h}}\right)^{2}+\Delta h^{2}}.&r_{0}\geq R_{\mathrm{h}}\end{cases}

If denoting F12(,;;){}_{2}F_{1}\left(\cdotp,\cdotp;\cdotp;\cdotp\right) as the standard Gaussian hypergeometric function, we have ω1(x,y)=F12(1,12x;22x;y)\omega_{1}\left(x,y\right)={}_{2}F_{1}\left(1,1-\frac{2}{x};2-\frac{2}{x};-y\right).

The probability density functions (PDFs) of r0r_{0}, θ\theta and Δh\Delta h are, respectively, given by

fr0(x)\displaystyle f_{r_{0}}\left(x\right) =2πλxexp(πλx2),x(0,)\displaystyle=2\pi\lambda x\exp\left(-\pi\lambda x^{2}\right),\>x\in\left(0,\infty\right) (10)
fθ(x)=\displaystyle f_{\theta}\left(x\right)= 12π,x[0,2π]\displaystyle\frac{1}{2\pi},\>x\in\left[0,2\pi\right] (11)
fΔh(x)\displaystyle f_{\Delta h}\left(x\right) =1h~Uh~L.x[h~L,h~U]\displaystyle=\frac{1}{\tilde{h}_{\mathrm{U}}-\tilde{h}_{\mathrm{L}}}.\>x\in\left[\tilde{h}_{\mathrm{L}},\tilde{h}_{\mathrm{U}}\right] (12)

Proof: Please refer to Appendix A-B.∎

It is observed from Proposition 1 that the CP and ST of UAV APs are dependent on the flight altitude of UAVs, UAV hover radius and mmWave backhaul constraint, etc. Moreover, it is shown from (9) that the expression of CP is in complicated form, which is due to the varying flight altitude of UAVs. To shed light on the impact of hover radius on the performance of UAV network, we analyze the upper bound of CP given a fixed flight altitude in the following corollary.

Corollary 1.

Given a fixed UAV flight altitude, the conditional CP of GU0\mathrm{GU}_{0} is upper bounded by

𝖢𝖯^π/2|ΔhSR=\displaystyle\hat{\mathsf{CP}}_{\pi/2\left|\Delta h\right.}^{\mathrm{SR}}= exp(πλδ(α,τ,qa)(Rh2+Δh2))1+δ(α,τ,qa)\displaystyle\frac{\exp\left(-\pi\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\left(R_{\mathrm{h}}^{2}+\Delta h^{2}\right)\right)}{1+\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}
×\displaystyle\times (1+δ(α,τ,qa)πλRh(1+δ(α,τ,qa))32(1+Erf(δ(α,τ,qa)πλRh1+δ(α,τ,qa)))),\displaystyle\left(1+\frac{\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\pi\sqrt{\lambda}R_{\mathrm{h}}}{\left(1+\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\right)^{\frac{3}{2}}}\left(1+\mathrm{Erf}\left(\frac{\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\sqrt{\pi\lambda}R_{\mathrm{h}}}{\sqrt{1+\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}}\right)\right)\right), (13)

where δ(α,τ,qa)=2qaτω1(α,τ)α2\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)=\frac{2q_{\mathrm{a}}\tau\omega_{1}\left(\alpha,\tau\right)}{\alpha-2} and Erf()\mathrm{Erf}\left(\cdot\right) denotes the standard error function.

Proof: Please refer to Appendix A-C.∎

Following Corollary 1, it can be shown that 𝖢𝖯^π/2|ΔhSRRh<0\frac{\partial\hat{\mathsf{CP}}_{\pi/2\left|\Delta h\right.}^{\mathrm{SR}}}{\partial R_{\mathrm{h}}}<0. This indicates that increasing the UAV hover radius would deteriorate the performance of UAV network under Semi-RTNA rule444Note that deconditioning Δh\Delta h through calculating the expectation of Δh\Delta h will not disprove the above results.. On the one hand, when the associated UAV is departing from the GU, the transmission distance will be increased, thereby reducing the desired signal power. On the other hand, when the neighboring UAVs are approaching the intended GU, inter-cell interference becomes severe and dominates the performance of air-to-ground transmission.

To verify the above results, we plot Fig. 2, which shows CP and network ST as a function of UAV hover radius RhR_{\mathrm{h}} under different UAV deployment densities λ\lambda and backhaul capacity constraints CtC_{\mathrm{t}}. It can be seen that either CP or network ST decreases with RhR_{\mathrm{h}} under different UAV densities. In particular, network ST is more significantly degraded given smaller CtC_{\mathrm{t}}, e.g., Ct=10C_{\mathrm{t}}=10bit/(sHz)\left(\mathrm{s}\cdot\mathrm{Hz}\right), and greater λ\lambda, e.g., λ>60/km2\lambda>60/\mathrm{km}^{2}. As discussed, the reason is that the hovering feature of UAVs would result in the connection discontinuity of GUs. In consequence, the UAV-GU transmission is more vulnerable to air-to-ground interference especially when more UAVs are deployed. Moreover, the mmWave mainlobe beamwidth has to be adaptively increased with the UAV hovering range so as to provide sufficient backhaul, which nevertheless reduces the mmWave mainlobe gain and limits the backhaul link capacity.

Refer to caption
(a) CP v.s. RhR_{\mathrm{h}}.
Refer to caption
(b) ST v.s. RhR_{\mathrm{h}}.
Figure 2: CP and network ST under different hover radius RhR_{\mathrm{h}}. For system parameters, we set P=30P=30dBm, Δh𝒰(90m,110m)\Delta h\sim\mathcal{U}\left(90\mathrm{m},110\mathrm{m}\right), τ=1\tau=1, α=3.5\alpha=3.5 and λGU=1×104\lambda_{\mathrm{GU}}=1\times 10^{4} user/km2\mathrm{km}^{2}. For backhaul parameters, we set Rm=500R_{\mathrm{m}}=500m. In this and following figures, numerical and simulation results are denoted by lines and markers, respectively.

III-B RTNA Rule

If frequent handover is allowed under RTNA, the performance of UAV network could intuitively be improved since GUs always connect to the closest hovering UAVs. In this light, we analyze CP and network ST under RTNA in the following proposition.

Proposition 2.

Each UAV equipped with Φ=π2\Phi=\frac{\pi}{2} directional antenna, network ST of the fixed-wing UAV network under RTNA rule is given by

𝖲𝖳π/2R\displaystyle\mathsf{ST}_{\pi/2}^{\mathrm{R}} =λa𝖢𝖯π/2Rmin{log2(1+τ),Cb}.\displaystyle=\lambda_{\mathrm{a}}\mathsf{CP}_{\pi/2}^{\mathrm{R}}\min\left\{\log_{2}\left(1+\tau\right),C_{\mathrm{b}}\right\}. (14)

In (14), 𝖢𝖯π/2R\mathsf{CP}_{\pi/2}^{\mathrm{R}} is given by

𝖢𝖯π/2R\displaystyle\mathsf{CP}_{\pi/2}^{\mathrm{R}} =Erf(δ(α,τ,qa)h~U)Erf(δ(α,τ,qa)h~L)2λδ(α,τ,qa)(1+δ(α,τ,qa))(h~Uh~L).\displaystyle=\frac{\mathrm{Erf}\left(\sqrt{\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}\tilde{h}_{\mathrm{U}}\right)-\mathrm{Erf}\left(\sqrt{\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}\tilde{h}_{\mathrm{L}}\right)}{2\sqrt{\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}\left(1+\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\right)\left(\tilde{h}_{\mathrm{U}}-\tilde{h}_{\mathrm{L}}\right)}. (15)

Proof: Please refer to Appendix A-D.∎

It is shown in Proposition 2 that the CP under RTNA is independent of the UAV hover radius. This is because, when hover radius is increased or decreased, GUs could always be associated with the closest UAV. In consequence, the problem raised by the hover of UAVs in Semi-RTNA rule will be alleviated, i.e., dominant interference due to approaching neighboring UAVs and reduced desired signal power due to departing associated UAVs.

Even though the results in Proposition 2 are still in complicated forms, we make a comparison of the two rules with the aid of Corollary 1. Deconditioning Δh\Delta h in (13), we have

𝔼Δh[𝖢𝖯^π/2|ΔhSR]\displaystyle\mathbb{E}_{\Delta h}\left[\hat{\mathsf{CP}}_{\pi/2\left|\Delta h\right.}^{\mathrm{SR}}\right] <(a)𝔼Δh[exp(πλδ(α,τ,qa)Δh2)1+δ(α,τ,qa)]\displaystyle\overset{\left(\mathrm{a}\right)}{<}\mathbb{E}_{\Delta h}\left[\frac{\exp\left(-\pi\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\Delta h^{2}\right)}{1+\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}\right]
=𝖢𝖯π/2R,\displaystyle=\mathsf{CP}_{\pi/2}^{\mathrm{R}}, (16)

where (a) follows by setting Rh=0R_{\mathrm{h}}=0 because 𝖢𝖯^π/2|ΔhSR\hat{\mathsf{CP}}_{\pi/2\left|\Delta h\right.}^{\mathrm{SR}} is a decreasing function of RhR_{\mathrm{h}}. Since 𝖢𝖯π/2SR<𝔼Δh[𝖢𝖯^π/2|ΔhSR]\mathsf{CP}_{\pi/2}^{\mathrm{SR}}<\mathbb{E}_{\Delta h}\left[\hat{\mathsf{CP}}_{\pi/2\left|\Delta h\right.}^{\mathrm{SR}}\right] (see Corollary 1) and the impact of mmWave backhaul is identical to RTNA and Semi-RTNA, (16) indicates that RTNA outperforms Semi-RTNA in terms of CP and network ST.

In Fig. 3, we plot CP and network ST with varying UAV deployment density λ\lambda under Semi-RTNA and RTNA. It is shown that RTNA outperforms Semi-RTNA in terms of CP and network ST. As indicated by Proposition 2, the CP is shown to be independent of the UAV hover radius under RTNA in Fig. 3a. Moreover, when λ\lambda is small, the network ST obtained by RTNA overlaps under different hover radii, which is shown in Fig. 3b. When λ\lambda increases, limited wireless backhaul will be shared by more UAVs. In consequence, it can be seen in Fig. 3b that network ST is more significantly degraded (solid lines) in the backhaul-limited region. Moreover, it can be seen that a greater hover radius would lead to a more rapid decrease of network ST. As discussed earlier, the reason is that the mmWave link capacity will be decreased with the high mobility of UAVs due to reduced mainlobe gain. Consequently, the UAV network performance is more likely to be dominated by limited backhaul.

In addition, compared to the terrestrial network, in which the coverage of terrestrial APs is limited, the coverage of UAVs is greatly enhanced due to higher deployment altitude, which increases the overlapping area of the UAV cells. As a result, it is shown in Fig. 3 that network ST begins to diminish with UAV density λ\lambda even when λ\lambda is small, e.g., λ>40/km2\lambda>40/\mathrm{km}^{2}, due to the inter-cell interference caused by overlapping UAV cells. This can be analytically verified through Proposition 2. For this reason, it is crucial to investigate how to alleviate the overwhelming interference in the UAV network, which is discussed in Section IV.

Refer to caption
(a) CP v.s. λ\lambda.
Refer to caption
(b) ST v.s. λ\lambda.
Figure 3: CP and Network ST under different UAV density λ\lambda. Except for RhR_{\mathrm{h}}, system parameters are identically set as those in Fig. 2.

IV Analysis and Optimization of Directional-Antenna UAV Network

The adjustment of directional-antenna UAV beamwidth is of great potential to avoid the overlap of UAV cells and mitigate the inter-cell interference. In this light, we first evaluate the performance of the directional-antenna UAV network under the two user association rules in the following.

IV-A Performance Analysis

When Semi-RTNA is applied, the SIR at GU0\mathrm{GU}_{0} can be expressed as

𝖲𝖨𝖱^GU0SR=\displaystyle\mathsf{\hat{SIR}}_{\mathrm{GU}_{0}}^{\mathrm{SR}}= PHUAV0G(ϕ,φ)d0αI^SR,\displaystyle\frac{PH_{\mathrm{UAV}_{0}}G\left(\phi,\varphi\right)d_{0}^{-\alpha}}{\hat{I}^{\mathrm{SR}}}, (17)

where G(ϕ,φ)G\left(\phi,\varphi\right) is given by (1) in Section II-C, I^SR=UAViΠ~UAV\UAV0PHUAViG(ϕ,φ)diα\hat{I}^{\mathrm{SR}}=\underset{\tiny{\mathrm{UAV}_{i}\in\mathrm{\tilde{\Pi}_{UAV}}\backslash\mathrm{UAV}_{0}}}{\sum}PH_{\mathrm{UAV}_{i}}G\left(\phi,\varphi\right)d_{i}^{-\alpha} denotes the interference stemming from the overlapping activated UAV cells. According to (3), in addition to transmission success, CP is as well dependent on whether GUs are within the projection area of the associated UAVs. Since UAV0GU0=r02+Rh22r0Rhcosθ\left\|\mathrm{UAV}_{0}-\mathrm{GU}_{0}\right\|=\sqrt{r_{0}^{2}+R_{\mathrm{h}}^{2}-2r_{0}R_{\mathrm{h}}\cos\theta} is the 2D distance from UAV0\mathrm{UAV}_{0} to GU0\mathrm{GU}_{0}, the projection probability ppp_{\mathrm{p}} is given by (18) to facilitate the analysis.

pp\displaystyle p_{\mathrm{p}} =(UAV0GU0<Rp)\displaystyle=\mathbb{P}\left(\left\|\mathrm{UAV}_{0}-\mathrm{GU}_{0}\right\|<R_{\mathrm{p}}\right)
=UAV0GU0<Rpr0λexp(πλr02)dθdr0.\displaystyle=\underset{\left\|\mathrm{UAV}_{0}-\mathrm{GU}_{0}\right\|<R_{\mathrm{p}}}{\iint}r_{0}\lambda\exp\left(-\pi\lambda r_{0}^{2}\right)d\theta dr_{0}. (18)

Following (17) and (18), CP and network ST are given in Proposition 3.

Proposition 3.

Each UAV equipped with directional-antenna Φ(0,π2]\Phi\in\left(0,\frac{\pi}{2}\right], network ST of the fixed-wing UAV network under Semi-RTNA rule is given by

𝖲𝖳dirSR\displaystyle\mathsf{ST}_{\mathrm{dir}}^{\mathrm{SR}} =λa𝖢𝖯dirSRmin{log2(1+τ),Cb}.\displaystyle=\lambda_{\mathrm{a}}\mathsf{CP}_{\mathrm{dir}}^{\mathrm{SR}}\min\left\{\log_{2}\left(1+\tau\right),C_{\mathrm{b}}\right\}. (19)

In (19), 𝖢𝖯dirSR\mathsf{CP}_{\mathrm{dir}}^{\mathrm{SR}} is given by

𝖢𝖯dirSR\displaystyle\mathsf{CP}_{\mathrm{dir}}^{\mathrm{SR}} =𝔼r0,θ,Δh[pp𝖢𝖯dir|r0,θ,ΔhSR|Ep].\displaystyle=\mathbb{E}_{r_{0},\theta,\Delta h}\left[p_{\mathrm{p}}\mathsf{CP}_{\mathrm{dir}\left|r_{0},\theta,\Delta h\right.}^{\mathrm{SR}}\left|E_{\mathrm{p}}\right.\right]. (20)

In (20), 𝖢𝖯dir|r0,θ,ΔhSR\mathsf{CP}_{\mathrm{dir}\left|r_{0},\theta,\Delta h\right.}^{\mathrm{SR}} is given by

𝖢𝖯dir|r0,θ,ΔhSR=\displaystyle\mathsf{CP}_{\mathrm{dir}\left|r_{0},\theta,\Delta h\right.}^{\mathrm{SR}}= {exp[πppλa(R^p2ω2(α,R^pατd0α)d^02ω2(α,d^0ατd0α))],r0Rhexp[πppλa(R^p2ω2(α,R^pατd0α)Δh2ω2(α,Δhατd0α))],otherwise\displaystyle\begin{cases}\exp\left[-\pi p_{\mathrm{p}}\lambda_{\mathrm{a}}\left(\hat{R}_{\mathrm{p}}^{2}\omega_{2}\left(\alpha,\frac{\hat{R}_{\mathrm{p}}^{\alpha}}{\tau d_{0}^{\alpha}}\right)-\hat{d}_{0}^{2}\omega_{2}\left(\alpha,\frac{\hat{d}_{0}^{\alpha}}{\tau d_{0}^{\alpha}}\right)\right)\right],&r_{0}\geq R_{\mathrm{h}}\\ \exp\left[-\pi p_{\mathrm{p}}\lambda_{\mathrm{a}}\left(\hat{R}_{\mathrm{p}}^{2}\omega_{2}\left(\alpha,\frac{\hat{R}_{\mathrm{p}}^{\alpha}}{\tau d_{0}^{\alpha}}\right)-\Delta h^{2}\omega_{2}\left(\alpha,\frac{\Delta h^{\alpha}}{\tau d_{0}^{\alpha}}\right)\right)\right],&\mathrm{otherwise}\end{cases} (21)

where R^p=Rp2+Δh2\hat{R}_{\mathrm{p}}=\sqrt{R_{\mathrm{p}}^{2}+\Delta h^{2}}, ω2(x,y)=F12(1,2x;1+2x;y)\omega_{2}\left(x,y\right)={}_{2}F_{1}\left(1,\frac{2}{x};1+\frac{2}{x};-y\right), d^0=(r0Rh)2+Δh2\hat{d}_{0}=\sqrt{\left(r_{0}-R_{\mathrm{h}}\right)^{2}+\Delta h^{2}} and d0=r02+Rh22r0Rhcosθ+Δh2d_{0}=\sqrt{r_{0}^{2}+R_{\mathrm{h}}^{2}-2r_{0}R_{\mathrm{h}}\cos\theta+\Delta h^{2}}.

Proof: Please refer to Appendix A-E.∎

Based on the results in Proposition 3, the impact of UAV projection radius RpR_{\mathrm{p}}, which is dependent on UAV flight altitude and beamwidth, on the performance of UAV network under Semi-RTNA can be revealed. To better illustrate the impact, we give Lemma 2 in the following.

Lemma 2.

Define Ψ(x)=x2F12(1,2z;1+2z;bxz)\Psi\left(x\right)=x^{2}{}_{2}F_{1}\left(1,\frac{2}{z};1+\frac{2}{z};-bx^{z}\right) (x0)\left(x\geq 0\right), where b(>0)b\left(>0\right) and z(>0)z\left(>0\right) are constant. Ψ(x)\Psi\left(x\right) is an increasing function of xx.

Proof: Please refer to Appendix A-F.∎

Following the projection probability in (18), more GUs can be covered by UAVs under a greater projection radius RpR_{\mathrm{p}} (or equivalently UAV beamwidth) by increasing the projection probability. According to Lemma 2, R^p2ω2(α,R^pατd0α)\hat{R}_{\mathrm{p}}^{2}\omega_{2}\left(\alpha,\frac{\hat{R}_{\mathrm{p}}^{\alpha}}{\tau d_{0}^{\alpha}}\right) in (21) is an increasing function of R^p=Rp2+Δh2\hat{R}_{\mathrm{p}}=\sqrt{R_{\mathrm{p}}^{2}+\Delta h^{2}}. Therefore, it indicates from Proposition 3 that increasing RpR_{\mathrm{p}} may as well reduce the coverage probability by increasing R^p\hat{R}_{\mathrm{p}} in (21). Intuitively, this is due to the degrading inter-cell interference among different UAV cells. For this reason, there is a tradeoff on adjusting the UAV projection radius in practical UAV network.

Based on Proposition 3, we further analyze the network ST of the directional-antenna UAV network under RTNA.

Proposition 4.

Each UAV equipped with directional-antenna Φ(0,π2]\Phi\in\left(0,\frac{\pi}{2}\right], network ST of the fixed-wing UAV network under RTNA rule is given by

𝖲𝖳dirR\displaystyle\mathsf{ST}_{\mathrm{dir}}^{\mathrm{R}} =λa𝖢𝖯dirRmin{log2(1+τ),Cb}.\displaystyle=\lambda_{\mathrm{a}}\mathsf{CP}_{\mathrm{dir}}^{\mathrm{R}}\min\left\{\log_{2}\left(1+\tau\right),C_{\mathrm{b}}\right\}. (22)

In (22), 𝖢𝖯dirR\mathsf{CP}_{\mathrm{dir}}^{\mathrm{R}} is given by

𝖢𝖯dirR\displaystyle\mathsf{CP}_{\mathrm{dir}}^{\mathrm{R}} =𝔼Δh,r0<Rp[p^pexp(πp^pλa(R^p2ω2(α,R^pατdˇ0α)dˇ02ω2(α,1τ)))],\displaystyle=\mathbb{E}_{\Delta h,r_{0}<R_{\mathrm{p}}}\left[\hat{p}_{\mathrm{p}}\exp\left(-\pi\hat{p}_{\mathrm{p}}\lambda_{\mathrm{a}}\left(\hat{R}_{\mathrm{p}}^{2}\omega_{2}\left(\alpha,\frac{\hat{R}_{\mathrm{p}}^{\alpha}}{\tau\check{d}_{0}^{\alpha}}\right)-\check{d}_{0}^{2}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)\right)\right)\right], (23)

where p^p=1exp(πλRp2)\hat{p}_{\mathrm{p}}=1-\exp\left(\pi\lambda R_{\mathrm{p}}^{2}\right) and dˇ0=r02+Δh2\check{d}_{0}=\sqrt{r_{0}^{2}+\Delta h^{2}}.

Proof: Please refer to Appendix A-G.∎

Similarly as the results under Semi-RTNA in Proposition 3, there is also a tradeoff on setting the projection radius RpR_{\mathrm{p}} under RTNA. According to (23), increasing RpR_{\mathrm{p}} may enhance the UAV coverage, while degrades the transmission success probability by introducing more inter-cell interference. Based on Propositions 3 and 4, we make a comparison of Semi-RTNA and RTNA with directional-antenna equipped by each UAV in the following.

Refer to caption
(a) 𝖲𝖳\mathsf{ST} v.s. Φ\Phi.
Refer to caption
(b) 𝖲𝖳\mathsf{ST} v.s. λ\lambda.
Figure 4: Network ST under Semi-RTNA and RTNA. System parameters are identically set as those in Fig. 2. In addition, we set Ct=10C_{\mathrm{t}}=10bit/(sHz)\left(\mathrm{s}\cdot\mathrm{Hz}\right).

Fig. 4 plots network ST under Semi-RTNA and RTNA. In particular, Fig. 4a plots network ST with varying UAV beamwidth under different UAV deployment densities. It is observed that increasing the UAV beamwidth could first enhance network ST since a greater air-to-ground coverage could be provided. Nevertheless, network ST will be degraded due to the dominant air-to-ground interference if UAV beamwidth further increases. Meanwhile, it is shown in Figs. 4a and 4b that a greater network ST can always be obtained under RTNA. Moreover, the optimal UAV beamwidth, which maximizes the network ST, inversely grows with the UAV deployment density. This indicates that the optimal UAV beamwidth is critically dependent on the UAV deployment density.

It is shown in Fig. 4 that appropriately decreasing the UAV beamwidth could reduce the overlap of UAV projection areas and therefore is of great potential in mitigating the inter-cell interference and enhancing network ST. For this reason, we study how to optimize the UAV beamwidth in the following.

IV-B Performance Optimization

In this part, we intend to maximize the network ST through tuning the UAV beamwidth. To this end, we first formulate the optimization problem as follows.

Φ\displaystyle\Phi^{*} =argmaxΦ𝖲𝖳dir(Φ)\displaystyle=\arg\underset{\Phi}{\max}\>\mathsf{ST}_{\mathrm{dir}}\left(\Phi\right) (24)
s.t.\displaystyle\mathrm{s.t.} Φ(0,π2].\displaystyle\Phi\in\left(0,\frac{\pi}{2}\right].

From Propositions 3 and 4, it is shown that the UAV half-beamwidth Φ\Phi has a complicated impact of network ST under either Semi-RTNA or RTNA rule. As a consequence, it is difficult to obtain a closed-form solution for the optimal half-beamwidth Φ\Phi^{*}. Even if the optimization problem in (24) is of single variable, numerically solving it may consume much time in practical applications. On this account, we propose a simple but effective PAE policy to tune the UAV beamwidth in the following.

It is shown in Fig. 4a that the optimal UAV beamwidth, which maximizes network ST, always decreases with the UAV deployment density. Moreover, it can be readily shown from (10) that the average 2D UAV-GU distance r0¯\bar{r_{0}} follows r0¯1λ\bar{r_{0}}\propto\frac{1}{\sqrt{\lambda}}. Therefore, the concept of the PAE policy is that the UAV projection radius RpR_{\mathrm{p}} should be inversely proportional to the square of activated UAV density λa\lambda_{\mathrm{a}}555Note that no coverage will be provided by the inactivated UAVs, which are connected by no GUs.. Especially, we have

Rp=\displaystyle R_{\mathrm{p}}= Cλa,\displaystyle\frac{C}{\sqrt{\lambda_{\mathrm{a}}}}, (25)

where C(>0)C\left(>0\right) is defined as the projection scaling parameter. A greater CC will result in a greater projection radius and equivalently UAV beamwidth. Therefore, supposing that the density of UAV APs is small and inter-cell interference is moderate, increasing CC would greatly increase the coverage of UAV APs and improve the UAV network performance. Following (25), the half-beamwidth of UAV is tuned as Φ=arctan(CΔhλa)\Phi^{*}=\arctan\left(\frac{C}{\Delta h\sqrt{\lambda_{\mathrm{a}}}}\right). Note that the projection radius of each UAV could be identically set as (25) through adjusting the beamwidth even if UAVs are in different flight altitudes. Moreover, perfect CSI is not required by the PAE policy. Therefore, the proposed policy can be applied in practical fixed-wing UAV networks, where it is difficult to acquire perfect CSI due to the high mobility of UAVs.

Since RTNA outperforms Semi-RTNA in terms of network ST, we then verify the efficiency of the PAE policy in improving network ST under RTNA. In particular, we study the scaling behavior of network ST with the growing UAV deployment density.

Theorem 1.

When each UAV adjusts the half-beamwidth according to Φ=arctan(CΔhλa)\Phi^{*}=\arctan\left(\frac{C}{\Delta h\sqrt{\lambda_{\mathrm{a}}}}\right), the scaling behavior of network ST under RTNA rule is given by

limλ\displaystyle\underset{\lambda\rightarrow\infty}{\mathrm{lim}} 𝖲𝖳dirR=λap^pδ1(p^p,C,α,τ)min{log2(1+τ),Cb},\displaystyle\mathsf{ST}_{\mathrm{dir}}^{\mathrm{R}}=\lambda_{\mathrm{a}}\hat{p}_{\mathrm{p}}^{\dagger}\delta_{1}\left(\hat{p}_{\mathrm{p}}^{\dagger},C,\alpha,\tau\right)\min\left\{\log_{2}\left(1+\tau\right),C_{\mathrm{b}}\right\}, (26)

where p^p=1exp(πC2)\hat{p}_{\mathrm{p}}^{\dagger}=1-\exp\left(-\pi C^{2}\right) and δ1(p^p,C,α,τ)=exp(p^pπC2ω2(α,1τ))exp(πC2)1p^pω2(α,1τ)\delta_{1}\left(\hat{p}_{\mathrm{p}}^{\dagger},C,\alpha,\tau\right)=\frac{\exp\left(-\hat{p}_{\mathrm{p}}^{\dagger}\pi C^{2}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)\right)-\exp\left(\pi C^{2}\right)}{1-\hat{p}_{\mathrm{p}}^{\dagger}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)}.

Proof: Please refer to Appendix A-H.∎

According to Theorem 1, we further study the network ST scaling behavior in the following two cases.

1) Backhaul unlimited case. This case occurs when sufficient backhaul could be provided through mmWave and GU density is limited. Accordingly, (26) in Theorem 1 degenerates into

limλ𝖲𝖳dirR\displaystyle\underset{\lambda\rightarrow\infty}{\mathrm{lim}}\mathsf{ST}_{\mathrm{dir}}^{\mathrm{R}} =(a)λGUp^pδ1(p^p,C,α,τ)log2(1+τ),\displaystyle\overset{\left(\mathrm{a}\right)}{=}\lambda_{\mathrm{GU}}\hat{p}_{\mathrm{p}}^{\dagger}\delta_{1}\left(\hat{p}_{\mathrm{p}}^{\dagger},C,\alpha,\tau\right)\log_{2}\left(1+\tau\right), (27)

where (a) follows because λa=qaλ\lambda_{\mathrm{a}}=q_{\mathrm{a}}\lambda and limλqa=λGU/λ\underset{\lambda\rightarrow\infty}{\mathrm{lim}}q_{\mathrm{a}}=\lambda_{\mathrm{GU}}/\lambda according to (6). It is shown in (27) that network ST will linearly increase with the GU density λGU\lambda_{\mathrm{GU}}. This indicates that greater network ST could be obtained if more GUs request service from UAVs. In other words, the proposed PAE policy could sustainably improve spectrum reuse as long as sufficient wireless backhaul is provided. This is fundamentally different from the cases in Figs. 3 and 4b, where network ST is degraded by over-deployment of UAVs supposing that constant-beamwidth directional-antenna is equipped by each UAV.

More importantly, it is shown in (27) that identical network ST can be obtained even if UAVs are in different altitudes with the application of the proposed PAE policy. For this reason, the optimization of beamwidth could help compensate for the loss of network ST, which is due to the increase flight altitude.

2) Backhaul limited case. This case occurs when either limited backhaul could be provided or GU density is sufficiently large. According to Theorem 1, (26) degenerates into

limλ𝖲𝖳dirR\displaystyle\underset{\lambda\rightarrow\infty}{\mathrm{lim}}\mathsf{ST}_{\mathrm{dir}}^{\mathrm{R}} =(a)p^pδ1(p^p,C,α,τ)Gm(ϑ)F|ε|(ϑ)CtπRm2.\displaystyle\overset{\left(\mathrm{a}\right)}{=}\frac{\hat{p}_{\mathrm{p}}^{\dagger}\delta_{1}\left(\hat{p}_{\mathrm{p}}^{\dagger},C,\alpha,\tau\right)G_{\mathrm{m}}\left(\vartheta\right)F_{\left|\varepsilon\right|}\left(\vartheta\right)C_{\mathrm{t}}}{\pi R_{\mathrm{m}}^{2}}. (28)

In this case, it is observed from (28) that network ST is limited by the mmWave link capacity and mmWave effective function, which is dependent on the UAV flight altitude and hover radius (see Section II-D). According to Lemma 1, Gm(ϑ)F|ε|(ϑ)G_{\mathrm{m}}\left(\vartheta\right)F_{\left|\varepsilon\right|}\left(\vartheta\right) is a decreasing function of ϑ=2arctan(RhhUAV)\vartheta=2\arctan\left(\frac{R_{\mathrm{h}}}{h_{\mathrm{UAV}}}\right). Therefore, the convergence value of network ST can be improved by increasing the flight altitude under the PAE policy in the backhaul-limited case.

Finally, we verify the efficiency of the proposed beamwidth optimization policy. In particular, Fig. 5 plots network ST as a function of UAV density under RTNA. Compared to the results in Fig. 4b, network ST could linearly increase and then converge with the growing UAV density when RTNA is applied under the PAE policy. Moreover, even if UAVs are deployed over different altitudes, identical network ST could be obtained supposing that the backhaul capacity is sufficiently large, i.e., CtC_{\mathrm{t}}\rightarrow\infty, which verifies the validity of the results in (27). Otherwise, when backhaul capacity is limited (see the lines under Ct=10bits/(sHz)C_{\mathrm{t}}=10\mathrm{bits/\left(\mathrm{s}\cdot\mathrm{Hz}\right)} or Ct=100bits/(sHz)C_{\mathrm{t}}=100\mathrm{bits/\left(\mathrm{s}\cdot\mathrm{Hz}\right)} in Fig. 5), increasing the flight altitude could enhance network ST under large UAV deployment density.

In practice, the UAV flight altitude has to be raised in specific scenarios such as emergency and cooperative engagement, which is detrimental to the spatial reuse of available spectrum resources. However, the proposed simple but effective PAE policy could significantly enhance the spatial reuse in higher flight altitude considering both backhaul unlimited and backhaul-unlimited cases. Moreover, it is worth noting that the PAE policy could be feasibly applied in the multi-antenna system [29, 30], where UAV APs are equipped with multiple antennas. The principle of PAE and multi-antenna techniques such as coordinated beamforming is to enhance the desired signal power and alleviate the interference power received at the intended receiver. The difference is that received power distribution will be different and accordingly the UAV beamwidth should be further optimized, which is dependent on the applied multi-antenna techniques.

Refer to caption
Figure 5: Network ST v.s. UAV density under PAE policy with projection scaling parameter C=1C=1. Note that RTNA is applied for user association. Except for Δh\Delta h, other system parameters are identically set as those in Fig. 2.

V Conclusion

In this paper, we have modeled and evaluated the performance of fixed-wing UAV network, where UAVs serve as APs to provide air-to-ground coverage to GUs with mmWave backhaul. It was shown that the hovering feature of UAV would result in unstable connections of GUs and reduce the mmWave mainlobe gain, which degrades CP and network ST. Even though the impact of UAV hovering is minor supposing that GUs are associated with the closest UAVs in a real-time manner, excessive handover overhead will be introduced. More importantly, it was shown that the over-deployment of UAVs would result in a rapid decrease of network ST since coverage of aerial APs is large. To enhance the performance of the fixed-wing UAV network, we propose a PAE policy to adaptively adjust UAV beamwidth according to the UAV deployment density. Notably, network was proved to increase with the UAV density and independent of the UAV flight altitude under PAE. Therefore, the results of this paper are helpful for the design, deployment and beamwidth optimization of fixed-wing UAV network.

Appendix A *

A-A Proof for Lemma 1

According to [23], we assume that the absolute mmWave beam orientation error |ε|\left|\varepsilon\right| follows exponential distribution, which is truncated to [0,π]\left[0,\pi\right]. Therefore, the PDF of |ε|\left|\varepsilon\right| is given by

f|ε|(x)\displaystyle f_{\left|\varepsilon\right|}\left(x\right) ={exp(xε¯)(1exp(πε¯))ε¯,x[0,π]0.otherwise\displaystyle=\begin{cases}\frac{\mathrm{exp}\left(-\frac{x}{\bar{\varepsilon}}\right)}{\left(1-\mathrm{exp}\left(-\frac{\pi}{\bar{\varepsilon}}\right)\right)\bar{\varepsilon}},&x\in\left[0,\pi\right]\\ 0.&\mathrm{otherwise}\end{cases} (29)

Since mmWave backhaul is effective only if |ε|ϑ2\left|\varepsilon\right|\leq\frac{\vartheta}{2}, the proof can be completed by computing 0ϑ/2f|ε|(x)dx\intop_{0}^{\vartheta/2}f_{\left|\varepsilon\right|}\left(x\right)\mathrm{d}x.

A-B Proof for Proposition 1

Since each UAV independently hovers with identical radius RhR_{\mathrm{h}}, hover angle θ\theta follows uniform distribution over [0,2π]\left[0,2\pi\right], i.e., θ𝒰(0,2π)\theta\sim\mathcal{U}\left(0,2\pi\right). In consequence, the 2D locations of UAVs still follow PPP with density λ\lambda [31]. Meanwhile, when the half-beamwidth of UAV directional antenna equals π2\frac{\pi}{2}, GUs are always in the projection area of the associated UAVs, i.e., pp=1p_{\mathrm{p}}=1. Therefore, the CP under the Semi-RTNA is given by

𝖢𝖯π/2SR=\displaystyle\mathsf{CP}_{\pi/2}^{\mathrm{SR}}= {PHUAV0d0αISR>τ}\displaystyle\mathbb{P}\left\{\frac{PH_{\mathrm{UAV}_{0}}d_{0}^{-\alpha}}{I^{\mathrm{SR}}}>\tau\right\}
=\displaystyle= {HUAV0>τISRPd0α}\displaystyle\mathbb{P}\left\{H_{\mathrm{UAV}_{0}}>\frac{\tau I^{\mathrm{SR}}}{Pd_{0}^{-\alpha}}\right\}
=\displaystyle= 𝔼r0,θ,Δh[ISR(s1)],\displaystyle\mathbb{E}_{r_{0},\theta,\Delta h}\left[\mathcal{L_{\mathit{I}^{\mathrm{SR}}}}\left(s_{1}\right)\right], (30)

where ISR(s1)=exp(s1ISR)\mathcal{L_{\mathit{I}^{\mathrm{SR}}}}\left(s_{1}\right)=\exp\left(-s_{1}I^{\mathrm{SR}}\right) denotes the Laplace Transform of ISRI^{\mathrm{SR}} evaluated at s1=τPd0αs_{1}=\frac{\tau}{Pd_{0}^{-\alpha}}. Note that d0=r02+Rh22Rhr0cosθ+Δh2d_{0}=\sqrt{r_{0}^{2}+R_{\mathrm{h}}^{2}-2R_{\mathrm{h}}r_{0}\cos\theta+\Delta h^{2}} follows according to Law of Cosines [32]. In particular, ISR(s1)\mathcal{L_{\mathit{I}^{\mathrm{SR}}}}\left(s_{1}\right) can be expressed as

ISR(s1)=(a)\displaystyle\mathcal{L_{\mathit{I}^{\mathrm{SR}}}}\left(s_{1}\right)\overset{\left(\mathrm{a}\right)}{=} exp[2πλalx(111+s1Pxα)dx]\displaystyle\exp\left[-2\pi\lambda_{\mathrm{a}}\int_{l}^{\infty}x\left(1-\frac{1}{1+s_{1}Px^{-\alpha}}\right)\mathrm{d}x\right]
=\displaystyle= exp[2πλaτd0αl2αω1(α,τd0αlα)α2],\displaystyle\exp\left[-2\pi\lambda_{\mathrm{a}}\frac{\tau d_{0}^{\alpha}l^{2-\alpha}\omega_{1}\left(\alpha,\tau d_{0}^{\alpha}l^{-\alpha}\right)}{\alpha-2}\right], (31)

where (a) follows owing to the probability generating functional (PFGL) of PPP [33]. Given r0Rhr_{0}\leq R_{\mathrm{h}}, the hovering interfering UAVs can be sufficiently close to GU0\mathrm{GU}_{0} such that l=Δhl=\Delta h. Otherwise, when r0Rhr_{0}\geq R_{\mathrm{h}}, we have l=(r0Rh)2+Δh2l=\sqrt{\left(r_{0}-R_{\mathrm{h}}\right)^{2}+\Delta h^{2}}.

In addition, the PDF of r0r_{0} can be obtained using contact distribution [33], i.e., fr0(x)=2πλxexp(πλx2),x(0,)f_{r_{0}}\left(x\right)=2\pi\lambda x\exp\left(-\pi\lambda x^{2}\right),x\in\left(0,\infty\right). Hence, CP and network ST can be obtained through substituting (31) into (30) and deconditioning r0r_{0}, θ\theta and Δh\Delta h.

A-C Proof for Corollary 1

According to (9) in Proposition 1, when UAV flight altitude is fixed, the conditional CP is given by

𝖢𝖯π/2|ΔhSR\displaystyle\mathsf{CP}_{\pi/2\left|\Delta h\right.}^{\mathrm{SR}} =𝔼r0,θ[exp(2πλaτd0αl2αω1(α,τd0αlα)α2)]\displaystyle=\mathbb{E}_{r_{0},\theta}\left[\exp\left(-2\pi\lambda_{\mathrm{a}}\frac{\tau d_{0}^{\alpha}l^{2-\alpha}\omega_{1}\left(\alpha,\tau d_{0}^{\alpha}l^{-\alpha}\right)}{\alpha-2}\right)\right]
<(a)𝔼r0,θ[exp(2πλaτd0αd^02αω1(α,τd0αlα)α2)]\displaystyle\overset{\left(\mathrm{a}\right)}{<}\mathbb{E}_{r_{0},\theta}\left[\exp\left(-2\pi\lambda_{\mathrm{a}}\frac{\tau d_{0}^{\alpha}\hat{d}_{0}^{2-\alpha}\omega_{1}\left(\alpha,\tau d_{0}^{\alpha}l^{-\alpha}\right)}{\alpha-2}\right)\right]
<(b)𝔼r0,θ[exp(2πλaτd^02ω1(α,τd0αlα)α2)]\displaystyle\overset{\left(\mathrm{b}\right)}{<}\mathbb{E}_{r_{0},\theta}\left[\exp\left(-2\pi\lambda_{\mathrm{a}}\frac{\tau\hat{d}_{0}^{2}\omega_{1}\left(\alpha,\tau d_{0}^{\alpha}l^{-\alpha}\right)}{\alpha-2}\right)\right]
<(c)𝔼r0[exp(πλδ(α,τ,qa)d^02)]\displaystyle\overset{\left(\mathrm{c}\right)}{<}\mathbb{E}_{r_{0}}\left[\exp\left(-\pi\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\hat{d}_{0}^{2}\right)\right]
=exp(πλδ(α,τ,qa)Δh2)𝔼r0[exp(πλδ(α,τ,qa)(r0Rh)2)]\displaystyle=\exp\left(-\pi\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\Delta h^{2}\right)\mathbb{E}_{r_{0}}\left[\exp\left(-\pi\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\left(r_{0}-R_{\mathrm{h}}\right)^{2}\right)\right]
=𝖢𝖯^π/2|ΔhSR(λ),\displaystyle=\hat{\mathsf{CP}}_{\pi/2\left|\Delta h\right.}^{\mathrm{SR}}\left(\lambda\right), (32)

where (a) follows because we substitute ll by d^0=(r0Rh)2+Δh2\hat{d}_{0}=\sqrt{\left(r_{0}-R_{\mathrm{h}}\right)^{2}+\Delta h^{2}} and d^0>l\hat{d}_{0}>l, (b) follows because d^0<d0=r02+Rh22r0Rhcosθ+Δh2\hat{d}_{0}<d_{0}=\sqrt{r_{0}^{2}+R_{\mathrm{h}}^{2}-2r_{0}R_{\mathrm{h}}\cos\theta+\Delta h^{2}}, and (c) follows because we substitute ll by d0d_{0}, l<d0l<d_{0} and ω1(x,y)\omega_{1}\left(x,y\right) is a decreasing function of xx (see Lemma 1 in [34]).

A-D Proof for Proposition 2

Similarly to the proof in Appendix A-B, 𝖢𝖯π/2R\mathsf{CP}_{\pi/2}^{\mathrm{R}} under RTNA is given by

𝖢𝖯π/2R=\displaystyle\mathsf{CP}_{\pi/2}^{\mathrm{R}}= 𝔼r0,Δh[IR(s2)],\displaystyle\mathbb{E}_{r_{0},\Delta h}\left[\mathcal{L_{\mathit{I}^{\mathrm{R}}}}\left(s_{2}\right)\right], (33)

where s2=τPdˇ0αs_{2}=\frac{\tau}{P\check{d}_{0}^{-\alpha}} and dˇ0=r02+Δh2\check{d}_{0}=\sqrt{r_{0}^{2}+\Delta h^{2}}. The Laplace Transform IR(s2)\mathcal{L_{\mathit{I}^{\mathrm{R}}}}\left(s_{2}\right) can be obtained by

IR(s2)=\displaystyle\mathcal{L_{\mathit{I}^{\mathrm{R}}}}\left(s_{2}\right)= 𝔼IR[es2IR]\displaystyle\mathbb{E}_{\mathit{I}^{\mathrm{R}}}\left[e^{-s_{2}\mathit{I}^{\mathrm{R}}}\right]
=\displaystyle= 𝔼IR[11+s2Pdˇ0α]\displaystyle\mathbb{E_{\mathrm{\mathit{I}^{\mathrm{R}}}}}\left[\prod\frac{1}{1+s_{2}P\check{d}_{0}^{-\alpha}}\right]
=(a)\displaystyle\overset{\left(\mathrm{a}\right)}{=} exp(2πλadˇ0x(111+s2Pxα)dx)\displaystyle\exp\left(-2\pi\lambda_{\mathrm{a}}\int_{\check{d}_{0}}^{\infty}x\left(1-\frac{1}{1+s_{2}Px^{-\alpha}}\right)\mathrm{d}x\right)
=\displaystyle= exp(2πλaτdˇ02F12(1,12α,22α,τ)α2)\displaystyle\exp\left(-2\pi\lambda_{\mathrm{a}}\frac{\tau\check{d}_{0}^{2}{}_{2}F_{1}\left(1,1-\frac{2}{\alpha},2-\frac{2}{\alpha},-\tau\right)}{\alpha-2}\right)
=(b)\displaystyle\overset{\left(\mathrm{b}\right)}{=} exp(πλδ(α,τ,qa)r02)exp(πλδ(α,τ,qa)Δh2),\displaystyle\exp\left(-\pi\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)r_{0}^{2}\right)\exp\left(-\pi\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\Delta h^{2}\right), (34)

where δ(α,τ,qa)=2qaτω1(α,τ)α2\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)=\frac{2q_{\mathrm{a}}\tau\omega_{1}\left(\alpha,\tau\right)}{\alpha-2}. In (34), (a) follows due to probability generating functional (PFGL) of PPP [33] and (b) follows due to dˇ0=r02+Δh2\check{d}_{0}=\sqrt{r_{0}^{2}+\Delta h^{2}} since GUs always connect to the closest UAVs for service under RTNA. In this case, the CP of GUs is independent of the hover angle.

With (33), (34) and the uniform distribution Δh𝒰(h~L,h~U)\Delta h\sim\mathcal{U}\left(\tilde{h}_{\mathrm{L}},\tilde{h}_{\mathrm{U}}\right), we obtain the 𝖢𝖯π/2R\mathsf{CP}_{\pi/2}^{\mathrm{R}} by deconditioning Δh\Delta h as

𝖢𝖯π/2R=\displaystyle\mathsf{CP}_{\pi/2}^{\mathrm{R}}= 𝔼Δh[exp(πλδ(α,τ,qa)Δh2)1+δ(α,τ,qa)]\displaystyle\mathbb{E}_{\Delta h}\left[\frac{\exp\left(-\pi\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\Delta h^{2}\right)}{1+\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}\right]
=\displaystyle= Erf(δ(α,τ,qa)h~U)Erf(δ(α,τ,qa)h~L)2λδ(α,τ,qa)(1+δ(α,τ,qa))(h~Uh~L).\displaystyle\frac{\mathrm{Erf}\left(\sqrt{\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}\tilde{h}_{\mathrm{U}}\right)-\mathrm{Erf}\left(\sqrt{\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}\tilde{h}_{\mathrm{L}}\right)}{2\sqrt{\lambda\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)}\left(1+\delta\left(\alpha,\tau,q_{\mathrm{a}}\right)\right)\left(\tilde{h}_{\mathrm{U}}-\tilde{h}_{\mathrm{L}}\right)}. (35)

A-E Proof for Proposition 3

When directional-antenna is equipped by each UAV under Semi-RTNA, 𝖢𝖯dir|r0,θ,ΔhSR\mathsf{CP}_{\mathrm{dir}\left|r_{0},\theta,\Delta h\right.}^{\mathrm{SR}} is given by

𝖢𝖯dir|r0,θ,ΔhSR=pp\displaystyle\mathsf{CP}_{\mathrm{dir}\left|r_{0},\theta,\Delta h\right.}^{\mathrm{SR}}=p_{\mathrm{p}} 𝔼r0,θ,Δh[^SR(s3)],\displaystyle\mathbb{E}_{r_{0},\theta,\Delta h}\left[\mathcal{L_{\hat{I}^{\mathrm{SR}}}}\left(s_{3}\right)\right], (36)

where ^SR(s3)\mathcal{L_{\hat{I}^{\mathrm{SR}}}}\left(s_{3}\right) denotes the Laplace Transform of I^SR\hat{I}^{\mathrm{SR}} at s3=τPG(ϕ,φ)d0αs_{3}=\frac{\tau}{PG\left(\phi,\varphi\right)d_{0}^{-\alpha}}. With the projection probability ppp_{\mathrm{p}} given by (18), ^SR(s3)\mathcal{L_{\hat{I}^{\mathrm{SR}}}}\left(s_{3}\right) can be obtained as

^SR(s3)=\displaystyle\mathcal{L_{\hat{I}^{\mathrm{SR}}}}\left(s_{3}\right)= 𝔼I^SR[exp(s3I^SR)]\displaystyle\mathbb{E}_{\hat{I}^{\mathrm{SR}}}\left[\exp\left(-s_{3}\hat{I}^{\mathrm{SR}}\right)\right]
=\displaystyle= 𝔼I^SR[11+s3PG(ϕ,φ)d0α]\displaystyle\mathbb{E_{\hat{\mathit{I}}^{\mathrm{SR}}}}\left[\prod\frac{1}{1+s_{3}PG\left(\phi,\varphi\right)d_{0}^{-\alpha}}\right]
=(a)\displaystyle\overset{\left(\mathrm{a}\right)}{=} exp(2πppλal^R^px(111+s3PG(ϕ,φ)xα)dx)\displaystyle\exp\left(-2\pi p_{\mathrm{p}}\lambda_{\mathrm{a}}\int_{\hat{l}}^{\hat{R}_{\mathrm{p}}}x\left(1-\frac{1}{1+s_{3}PG\left(\phi,\varphi\right)x^{-\alpha}}\right)\mathrm{d}x\right)
=\displaystyle= exp(πppλa(R^p2ω2(α,R^pατd0α)l^2ω2(α,l^ατd0α))),\displaystyle\exp\left(-\pi p_{\mathrm{p}}\lambda_{\mathrm{a}}\left(\hat{R}_{\mathrm{p}}^{2}\omega_{2}\left(\alpha,\frac{\hat{R}_{\mathrm{p}}^{\alpha}}{\tau d_{0}^{\alpha}}\right)-\hat{l}^{2}\omega_{2}\left(\alpha,\frac{\hat{l}^{\alpha}}{\tau d_{0}^{\alpha}}\right)\right)\right), (37)

where (a) follows since the inter-cell interference is limited due to the projection of directional-antenna UAVs. Specifically, interference only stems from the set of UAVs with density ppλap_{\mathrm{p}}\lambda_{\mathrm{a}}, which are distributed within the projection area around GU0\mathrm{GU}_{0} with radius RpR_{\mathrm{p}}.

In (37), when r0RhRpr_{0}\leq R_{\mathrm{h}}\leq R_{\mathrm{p}} holds, the distance from the interfering UAVs to GU0\mathrm{GU}_{0} is greater than l^=Δh\hat{l}=\Delta h. Otherwise, when r0>Rhr_{0}>R_{\mathrm{h}}, the distance from the closest interfering UAVs to GU0\mathrm{GU}_{0} is l^=(r0Rh)2+Δh2\hat{l}=\sqrt{\left(r_{0}-R_{\mathrm{h}}\right)^{2}+\Delta h^{2}}.

Integrating (36) and (37), we complete the proof.

A-F Proof for Lemma 2

The proof is completed by showing the derivative of Ψ(x)\Psi\left(x\right) in terms of xx is greater than 0. Specifically, the derivative is given by

dΨ(x)dx\frac{\mathrm{d}\Psi\left(x\right)}{\mathrm{d}x}= 2x1+bxz.\displaystyle\frac{2x}{1+bx^{z}}. (38)

Since b>0b>0, dΨ(x)dx>0\frac{\mathrm{d}\Psi\left(x\right)}{\mathrm{d}x}>0 holds.

A-G Proof for Proposition 4

When directional-antenna is equipped by each UAV under RTNA, we have p^p=1exp(πλRp2)\hat{p}_{\mathrm{p}}=1-\exp\left(\pi\lambda R_{\mathrm{p}}^{2}\right). Therefore, 𝖢𝖯dirR\mathsf{CP}_{\mathrm{dir}}^{\mathrm{R}} is given by

𝖢𝖯dirR\displaystyle\mathsf{CP}_{\mathrm{dir}}^{\mathrm{R}} =𝔼r0,Δh[p^p^R(s4)]\displaystyle=\mathbb{E}_{r_{0},\Delta h}\left[\hat{p}_{\mathrm{p}}\mathcal{L_{\hat{I}^{\mathrm{R}}}}\left(s_{4}\right)\right]
=𝔼r0,Δh,I^R[p^pexp(s4I^R)]\displaystyle=\mathbb{E}_{r_{0},\Delta h,\hat{I}^{\mathrm{R}}}\left[\hat{p}_{\mathrm{p}}\exp\left(-s_{4}\hat{I}^{\mathrm{R}}\right)\right]
=𝔼r0,Δh[p^pexp(2πp^pλadˇ0R^px(111+s4PG(ϕ,φ)xα)dx)]\displaystyle=\mathbb{E}_{r_{0},\Delta h}\left[\hat{p}_{\mathrm{p}}\exp\left(-2\pi\hat{p}_{\mathrm{p}}\lambda_{\mathrm{a}}\int_{\check{d}_{0}}^{\hat{R}_{\mathrm{p}}}x\left(1-\frac{1}{1+s_{4}PG\left(\phi,\varphi\right)x^{-\alpha}}\right)\mathrm{d}x\right)\right]
=𝔼r0,Δh[p^pexp[πp^pλa(R^p2ω2(α,R^pατdˇ0α)dˇ02ω2(α,1τ))]],\displaystyle=\mathbb{E}_{r_{0},\Delta h}\left[\hat{p}_{\mathrm{p}}\exp\left[-\pi\hat{p}_{\mathrm{p}}\lambda_{\mathrm{a}}\left(\hat{R}_{\mathrm{p}}^{2}\omega_{2}\left(\alpha,\frac{\hat{R}_{\mathrm{p}}^{\alpha}}{\tau\check{d}_{0}^{\alpha}}\right)-\check{d}_{0}^{2}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)\right)\right]\right], (39)

where ^R(s4)\mathcal{L_{\hat{I}^{\mathrm{R}}}}\left(s_{4}\right) denotes the Laplace Transform of I^R\hat{I}^{\mathrm{R}} at s4=τPG(ϕ,φ)dˇ0αs_{4}=\frac{\tau}{PG\left(\phi,\varphi\right)\check{d}_{0}^{-\alpha}}, dˇ0=r02+Δh2\check{d}_{0}=\sqrt{r_{0}^{2}+\Delta h^{2}} and I^R=UAViΠ~UAV\UAV0PHUAViG(ϕ,φ)diα\hat{I}^{\mathrm{R}}=\underset{\tiny{\mathrm{UAV}_{i}\in\mathrm{\tilde{\Pi}_{UAV}}\backslash\mathrm{UAV}_{0}}}{\sum}PH_{\mathrm{UAV}_{i}}G\left(\phi,\varphi\right)d_{i}^{-\alpha} is the interference stemming from the overlapping activated UAV cells. The remaining parts of the proof are identical to those in Appendix A-D and thus omitted due to space limitation.

A-H Proof for Theorem 1

According to Proposition 4, given a fixed flight altitude, the conditional CP of GUs is given by

limλ𝖢𝖯dir|ΔhR(λ)=\displaystyle\underset{\lambda\rightarrow\infty}{\mathrm{lim}}\mathsf{CP}_{\mathrm{dir}\left|\Delta h\right.}^{\mathrm{R}}\left(\lambda\right)= p^p𝔼r0<Rp[exp(πp^pλa(R^p2ω2(α,R^pατdˇ0α)dˇ02ω2(α,1τ)))]\displaystyle\hat{p}_{\mathrm{p}}\mathbb{E}_{r_{0}<R_{\mathrm{p}}}\left[\exp\left(-\pi\hat{p}_{\mathrm{p}}\lambda_{\mathrm{a}}\left(\hat{R}_{\mathrm{p}}^{2}\omega_{2}\left(\alpha,\frac{\hat{R}_{\mathrm{p}}^{\alpha}}{\tau\check{d}_{0}^{\alpha}}\right)-\check{d}_{0}^{2}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)\right)\right)\right]
=(a)\displaystyle\overset{\left(\mathrm{a}\right)}{=} p^p𝔼r0<Rp[exp(πp^pλaω2(α,1τ)(Rp2r02))]\displaystyle\hat{p}_{\mathrm{p}}\mathbb{E}_{r_{0}<R_{\mathrm{p}}}\left[\exp\left(-\pi\hat{p}_{\mathrm{p}}\lambda_{\mathrm{a}}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)\left(R_{\mathrm{p}}^{2}-r_{0}^{2}\right)\right)\right]
=\displaystyle= p^pexp(p^pπλaRp2ω2(α,1τ))exp(πλaRp2)1p^pω2(α,1τ),\displaystyle\hat{p}_{\mathrm{p}}\frac{\exp\left(-\hat{p}_{\mathrm{p}}\pi\lambda_{\mathrm{a}}R_{\mathrm{p}}^{2}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)\right)-\exp\left(-\pi\lambda_{\mathrm{a}}R_{\mathrm{p}}^{2}\right)}{1-\hat{p}_{\mathrm{p}}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)}, (40)

where (a) follows because limλRp0\underset{\lambda\rightarrow\infty}{\lim}R_{\mathrm{p}}\rightarrow 0 and limλr00\underset{\lambda\rightarrow\infty}{\lim}r_{0}\rightarrow 0. Substituting Rp=CλaR_{\mathrm{p}}=\frac{C}{\sqrt{\lambda_{\mathrm{a}}}} into (40), we have

limλ𝖢𝖯dir|ΔhR(λ)\displaystyle\underset{\lambda\rightarrow\infty}{\mathrm{lim}}\mathsf{CP}_{\mathrm{dir}\left|\Delta h\right.}^{\mathrm{R}}\left(\lambda\right) =p^pexp(p^pπC2ω2(α,1τ))exp(πC2)1p^pω2(α,1τ).\displaystyle=\hat{p}_{\mathrm{p}}^{\dagger}\frac{\exp\left(-\hat{p}_{\mathrm{p}}^{\dagger}\pi C^{2}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)\right)-\exp\left(\pi C^{2}\right)}{1-\hat{p}_{\mathrm{p}}^{\dagger}\omega_{2}\left(\alpha,\frac{1}{\tau}\right)}. (41)

In (41), limλ𝖢𝖯dir|ΔhR(λ)\underset{\lambda\rightarrow\infty}{\mathrm{lim}}\mathsf{CP}_{\mathrm{dir}\left|\Delta h\right.}^{\mathrm{R}}\left(\lambda\right) is independent of UAV flight altitude. Therefore, limλ𝖢𝖯dirR(λ)=limλ𝖢𝖯dir|ΔhR(λ)\underset{\lambda\rightarrow\infty}{\mathrm{lim}}\mathsf{CP}_{\mathrm{dir}}^{\mathrm{R}}\left(\lambda\right)=\underset{\lambda\rightarrow\infty}{\mathrm{lim}}\mathsf{CP}_{\mathrm{dir}\left|\Delta h\right.}^{\mathrm{R}}\left(\lambda\right). Hence, we complete the proof.

References

  • [1] S. Koulali, E. Sabir, T. Taleb, and M. Azizi, “A green strategic activity scheduling for UAV networks: A sub-modular game perspective,” IEEE Commun. Mag., vol. 54, no. 5, pp. 58–64, May. 2016.
  • [2] Y. Dong, M. Z. Hassan, J. Cheng, M. J. Hossain, and V. C. M. Leung, “An edge computing empowered radio access network with UAV-mounted FSO fronthaul and backhaul: Key challenges and approaches,” IEEE Wireless Commun., vol. 25, no. 3, pp. 154–160, Jun. 2018.
  • [3] N. Zhao, F. R. Yu, L. Fan, Y. Chen, J. Tang, A. Nallanathan, and V. C. M. Leung, “Caching unmanned aerial vehicle-enabled small-cell networks: Employing energy-efficient methods that store and retrieve popular content,” IEEE Veh. Technol. Mag., vol. 14, no. 1, pp. 71–79, Mar. 2019.
  • [4] J. Liu, M. Sheng, R. Lyu, and J. Li, “Performance analysis and optimization of UAV integrated terrestrial cellular network,” IEEE Internet Things Journal, vol. 6, no. 2, pp. 1841–1855, Apr. 2019.
  • [5] M. Erdelj, E. Natalizio, K. R. Chowdhury, and I. F. Akyildiz, “Help from the sky: Leveraging UAVs for disaster management,” IEEE Pervasive Comput., vol. 16, no. 1, pp. 24–32, Jan. 2017.
  • [6] Huawei, “5G: Huawei and Vodafone achieve 20 Gbps for single-user outdoor at E-band,” Huawei, Tech. Rep., 2016. [Online]. Available: http://www.huawei.com/en/press-events/news/2016/7/huawei-vodafone-5g-test
  • [7] M. Mozaffari, A. Taleb Zadeh Kasgari, W. Saad, M. Bennis, and M. Debbah, “Beyond 5G with UAVs: Foundations of a 3D wireless cellular network,” IEEE Trans. Wireless Commun., vol. 18, no. 1, pp. 357–372, Jan. 2019.
  • [8] F. Cheng, S. Zhang, Z. Li, Y. Chen, N. Zhao, F. R. Yu, and V. C. M. Leung, “UAV trajectory optimization for data offloading at the edge of multiple cells,” IEEE Trans. Veh. Technol., vol. 67, no. 7, pp. 6732–6736, Jul. 2018.
  • [9] A. V. Savkin and H. Huang, “Deployment of unmanned aerial vehicle base stations for optimal quality of coverage,” IEEE Wireless Commun. Lett., vol. 8, no. 1, pp. 321–324, Feb. 2019.
  • [10] X. Sun and N. Ansari, “Latency aware drone base station placement in heterogeneous networks,” in Proc. IEEE GLOBECOM, Singapore, Dec. 2017, pp. 1–6.
  • [11] C. Zhang and W. Zhang, “Spectrum sharing for drone networks,” IEEE J. Sel. Areas Commun., vol. 35, no. 1, pp. 136–144, Jan. 2017.
  • [12] KDDI, “Terra drone, KDDI launch new UAS services,” KDDI, Tech. Rep., 2019.
  • [13] AT&T, “When COWs fly: At&t sending LTE signals from drones,” AT&T, Tech. Rep., 2019. [Online]. Available: https://about.att.com/innovationblog/cows_fly
  • [14] Y. Zeng and R. Zhang, “Energy-efficient UAV communication with trajectory optimization,” IEEE Trans. Wireless Commun., vol. 16, no. 6, pp. 3747–3760, Jun. 2017.
  • [15] S. Jeong, O. Simeone, and J. Kang, “Mobile edge computing via a UAV-mounted Cloudlet: Optimization of bit allocation and path planning,” IEEE Trans. Veh. Technol., vol. 67, no. 3, pp. 2049–2063, March 2018.
  • [16] J. Zhou, Y. Su, and P. Li, “An emergency mobile communication system based on fixed-wing drone and satellite transmission access,” in Proc. IEEE DASC/PiCom/CBDCom/CyberSciTech, Fukuoka, Japan, Aug. 2019, pp. 465–470.
  • [17] R. He, B. Ai, G. L. Stüber, G. Wang, and Z. Zhong, “Geometrical-based modeling for millimeter-wave MIMO mobile-to-mobile channels,” IEEE Trans. Veh. Technol., vol. 67, no. 4, pp. 2848–2863, Apr. 2018.
  • [18] R. He, B. Ai, G. Wang, Z. Zhong, C. Schneider, D. A. Dupleich, R. S. Thomae, M. Boban, J. Luo, and Y. Zhang, “Propagation channels of 5G millimeter-wave vehicle-to-vehicle communications: Recent advances and future challenges,” IEEE Veh. Technol. Mag., vol. 15, no. 1, pp. 16–26, Mar. 2020.
  • [19] G. Wang, Q. Liu, R. He, F. Gao, and C. Tellambura, “Acquisition of channel state information in heterogeneous cloud radio access networks: challenges and research directions,” IEEE Wireless Commun., vol. 22, no. 3, pp. 100–107, Jun. 2015.
  • [20] Z. Yang, C. Pan, M. Shikh-Bahaei, W. Xu, M. Chen, M. Elkashlan, and A. Nallanathan, “Joint altitude, beamwidth, location, and bandwidth optimization for UAV-enabled communications,” IEEE Commun. Lett., vol. 22, no. 8, pp. 1716–1719, Aug. 2018.
  • [21] H. Zhu and J. Wang, “Chunk-based resource allocation in OFDMA systems - part I: chunk allocation,” IEEE Trans. Commun., vol. 57, no. 9, pp. 2734–2744, Sep. 2009.
  • [22] ——, “Chunk-based resource allocation in OFDMA systems - part II: Joint chunk, power and bit allocation,” IEEE Trans. Commun., vol. 60, no. 2, pp. 499–509, Feb. 2012.
  • [23] J. Wildman, P. H. J. Nardelli, M. Latva-aho, and S. Weber, “On the joint impact of beamwidth and orientation error on throughput in directional wireless poisson networks,” IEEE Trans. Wireless Commun., vol. 13, no. 12, pp. 7072–7085, Dec. 2014.
  • [24] W. Feng, Y. Li, D. Jin, L. Su, and S. Chen, “Millimetre-wave backhaul for 5G networks: Challenges and solutions,” IEEE Trans. Veh. Technol., vol. 16, no. 2, p. 892, Jun. 2016.
  • [25] B. Ejder, B. Mehdi, K. Marios, and D. Mérouane, “Cache-enabled small cell networks: modeling and tradeoffs,” EURASIP J. Wireless Commun. Netw., vol. 2015, no. 1, pp. 1–11, Feb. 2015.
  • [26] Z. Yan, S. Chen, Y. Ou, and H. Liu, “Energy efficiency analysis of cache-enabled two-tier hetnets under different spectrum deployment strategies,” IEEE Access, vol. 5, pp. 6791–6800, Mar. 2017.
  • [27] S. Lee and K. Huang, “Coverage and economy of cellular networks with many base stations,” IEEE Commun. Lett., vol. 16, no. 7, pp. 1038–1040, Jul. 2012.
  • [28] E. Turgut and M. C. Gursoy, “Energy efficiency in relay-assisted mmWave cellular networks,” in Proc. IEEE VTC-Fall, Montreal, Canada, Sep. 2016, pp. 1–5.
  • [29] C. Xing, S. Ma, and Y. Zhou, “Matrix-monotonic optimization for MIMO systems,” IEEE Trans. Signal Process., vol. 63, no. 2, pp. 334–348, Jan. 2015.
  • [30] Y. Dong, M. J. Hossain, J. Cheng, and V. C. M. Leung, “Cross-layer scheduling and beamforming in smart-grid powered cellular networks with heterogeneous energy coordination,” IEEE Trans. Commun., to be published.
  • [31] M. Haenggi and R. K. Ganti, Interference in large wireless networks.   Now Publishers Inc, 2009.
  • [32] F. Eriksson, “The law of sines for tetrahedra and n-simplices,” Geometriae Dedicata, vol. 7, no. 1, pp. 71–80, 1978.
  • [33] D. Stoyan, W. S. Kendall, J. Mecke, and L. Ruschendorf, Stochastic geometry and its applications.   Wiley Chichester, 1995, vol. 2.
  • [34] J. Liu, M. Sheng, L. Liu, and J. Li, “Effect of densification on cellular network performance with bounded pathloss model,” IEEE Commun. Lett., vol. 21, no. 2, pp. 346–349, Feb. 2017.