1

Sum-Rate Analysis for High Altitude Platform (HAP) Drones
with Tethered Balloon Relay

arXiv:1712.06583v1 [cs.IT] 18 Dec 2017

Sudheesh P. G, Mohammad Mozaffari, Student Member, IEEE, Maurizio Magarini, Member, IEEE, Walid
Saad, Senior Member, IEEE, and P. Muthuchidambaranathan

Abstract—High altitude platform (HAP) drones can provide
broadband wireless connectivity to ground users in rural areas
by establishing line-of-sight (LoS) links and exploiting effective
beamforming techniques. However, at high altitudes, acquiring
the channel state information (CSI) for HAPs, which is a key
component to perform beamforming, is challenging. In this
paper, by exploiting an interference alignment (IA) technique,
a novel method for achieving the maximum sum-rate in HAPbased communications without CSI is proposed. In particular,
to realize IA, a multiple-antenna tethered balloon is used as a
relay between multiple HAP drones and ground stations (GSs).
Here, a multiple-input multiple-output X network system is
considered. The capacity of the considered M ×N X network with
a tethered balloon relay is derived in closed-form. Simulation
results corroborate the theoretical findings and show that the
proposed approach yields the maximum sum-rate in multiple
HAPs-GSs communications in absence of CSI. The results also
show the existence of an optimal balloon’s altitude for which the
sum-rate is maximized.

I. I NTRODUCTION
Satellite and terrestrial cellular communications are the two
most widely used communication systems for providing global
connectivity to mobile ground users. While satellites can
deliver wireless service to users in remote areas, their spectral
efficiency is limited by their large footprints [1]. Meanwhile,
terrestrial communication systems cannot guarantee a reliable
service for users in remote, rural areas, due to the lack of
infrastructure nodes, such as base stations (BSs). High altitude
platform (HAP) drones can substantially extend the coverage
of terrestrial networks by establishing line-of-sight (LoS) links
and adjusting their altitude [2], [3].
To exploit the spatial dimension and enhance spectral efficiency, HAPs will typically rely on highly directive antennas to
communicate with ground stations [4]. In a single HAP system
with multiple antennas at the transmitter, a spatial multiplexing
gain cannot be typically achieved due to a high correlation
between parallel paths [5]. However, the deployment of multiple spatially separated HAPs can be a promising solution
to exploit spatial multiplexing and boost spectral efficiency.
In particular, by using a large number of antennas at the
HAPs, one can provide a precise beamforming which is a key
requirement for spatial multiplexing. To this end, channel state
information (CSI) at the transmitter (CSIT) is required [5],
Sudheesh P. G and P. Muthuchidambaranathan are with Department of Electronics and Communication Engineering, National Institute of Technology,
Tiruchirappalli, India (e-mail: pgsudheesh@gmail.com; muthuc@nitt.edu).
Mohammad Mozaffari and Walid Saad are with Wireless@VT, Electrical
and Computer Engineering Department, Virginia Tech, Blacksburg, VA,
24061, USA, (e-mail: mmozaff@vt.edu; walids@vt.edu).
Maurizio Magarini is with Dipartimento di Elettronica, Informazione e
Bioingegneria Politecnico di Milano, 20133 Milano, Italy (e-mail: maurizio.magarini@polimi.it).

[6]. However, in HAP drone systems, acquiring precise CSIT
is challenging due to the high altitudes and the movement of
the drones. Consequently, in HAPs-to-ground stations (GSs)
communications1 exploiting spatial multiplexing, which can
yield a maximum possible sum-rate, is also challenging.
One practical approach to achieve a maximum possible sumrate for HAP drones-to-GSs communications is via the use of
interference alignment (IA) schemes [6]–[8]. In particular, at
a high SNR regime, which is a typical case in HAP communications, IA can achieve a maximum sum-rate by restricting
interference beams to a smaller subspace that does not overlap
with the desired signal space [6]. Unlike terrestrial wireless
systems in which the BSs’ positions are fixed, acquiring CSIT
to implement IA in HAPs-GSs communication systems is challenging due to imperfect HAP drone stabilization. As a result,
the lack of exact CSI at HAP drones can yield a significant
degradation of the sum-rate performance. Nevertheless, with
the use of relays, it is possible to achieve maximum sum-rate
when CSIT is not available [9].
Unlike the time domain realization of IA, the performance
of IA in the spatial domain is limited due to the restrictions
in designing precoding matrices [6]. Therefore, the maximum
possible sum-rate can be achieved by implementing relayassisted IA in the time domain. In this case, one can use
popular relaying mechanisms such as amplify and forward
(AF) or decode and forward (DF) [9]. An AF relaying scheme
with multiple relays achieves an upper bound for the DoF of a
generalized X network [9]. Hence, by adopting an AF relaying
scheme such as the one in [9], multiple relays can be used to
achieve maximum possible sum-rate of the system where HAP
drones have knowledge of CSI. Similar results can be achieved
with a DF-based relaying scheme with single relay in a twouser X channel [10]. However, the previous works in [9] and
[10] did not investigate the use of a DF relaying mechanism
for IA in an HAP drones system.
The main contribution of this paper is a novel framework
for maximizing the sum-rate of a relay-aided HAP drones
wireless system when the CSI is not available. In particular,
to achieve the maximum sum-rate, we propose a DF scheme
involving M HAP drones, N ground receivers, and one relay
with (M − 1) × (N − 1) antennas. In this scenario, we show
that it is possible to achieve the maximum possible sum-rate
by exploiting the IA scheme. Moreover, we derive a closed1 An HAP drones-GSs wireless system in which each HAP drone carries
a dedicated symbol for each GS can be modeled as an X network. The X
network houses all possible channel models such as the interference channel,
the Y channel, and the Z channel. Beyond offering a generalized structure, an
X network offers a maximum capacity as compared to other channel models
such as the interference channel [6].

2

form analytical expression for the capacity of an M × N X
channel with tethered balloon relay. Simulation results verify
our analytical results and show that a significant sum-rate gain
can be achieved by using the proposed scheme.
II. S YSTEM M ODEL
Consider a geographical area with N GSs (or receivers)
and a tethered balloon attached to a control station, as shown
in Fig. 1. This control station provides the power required to
operate the tethered balloon. Meanwhile, the GSs receive data
from M HAPs, which are located at altitudes within the range
of 17-22 km [1]. Each HAP and receiver houses A antennas
while the tethered balloon has (M − 1) × (N − 1) antennas.
Each GS receives data from each HAP, forming an M × N
X network. Unlike [4], which uses a frequency duplexing
technique to avoid the interference in HAP communications,
our proposed model uses a relay that operates in half duplex
mode in the same frequency band.
A. Channel Model
For terrestrial communications, the channel is typically
modeled as Rayleigh in urban areas and Rician in suburban scenarios. However, in air-to-ground communications,
the channel has different characteristics [11], [12]. In urban
environments, the air-to-ground channel experiences Rician
fading due to the presence of LoS links. In suburban areas, a
Rayleigh fading is experienced due to the presence of reflected
signals which are stronger than LoS signals [11].
Here, we adopt a Rician channel model in which both
LoS and non-LoS (NLoS) paths are considered. Therefore,
the channel gain matrix can be represented as [13]:
r
r
κ
1
H LoS +
H NLoS ,
(1)
H=
1+κ
1+κ
where H LoS and H NLoS represent, respectively, the channel
matrices for LoS and non-LoS
communication. The Rician
σ2
2
2
factor κ is given by, κ = σ2LoS [14], where σLoS
and σNLoS
NLoS
are the power of LoS path and NLoS path, respectively. For
our model, we consider A antennas at each HAP and GS.
The role of HAP and GS as transmitter or receiver can be
reversed, using the reciprocity property. The static MIMO
channel, excluding the path loss, is given by [14]:
 
T

1
1
dT
dR
 


 ej2π λ sin(θA )   ej2π λ sin(θD ) 

, (2)


H̄ LoS =
..
..
·

.
.
 


dR
dT
j2π λ (M −1)
j2π λ (N −1)
sin(θA )
sin(θD )
e
e
where dR and dT are the antenna spacing at the receiver
and transmitter (L  dR , dT ), and λ is the wavelength.
Also, θA and θD represent, respectively, the angle-of-arrival
at the receiver and angle-of-departure at the transmitters. We
also note that the NLoS MIMO channel follows a Rayleigh
distribution.
III. DF RELAY AIDED I NTERFERENCE A LIGNMENT FOR
SYSTEMS WITHOUT CSIT
Knowledge of CSIT is a key requirement for designing a
precoder at the HAPs. However, at high altitudes acquiring a

Fig. 1: System model.
precise CSIT is challenging due to the difficulty in stabilizing
the aerial platform that is affected by wind and other natural
factors. In such a case, to achieve IA, a feasible solution
is to accommodate a tethered balloon relay to perform DF
operation. A DF scheme will require at least M antennas
to decode M transmitted symbols. However, the DF relay
must accommodate (M − 1) × (N − 1) antennas to perform
IA. By exploiting temporal domain characteristics, we can
transmit M N number of symbols in M + N − 1 time slots
for transmitter and receiver.
A. Feasibility of DF in tethered balloon relay
For IA, the tethered balloon must have (M − 1) × (N − 1)
antennas. To tightly pack these antennas in a tethered balloon
relay, the HAPs must be placed sufficiently apart so that the
links from HAPs to GSs become uncorrelated. We deduce the
minimal separation required between HAPs, when they are
located at 18 km above earth, operating at 48 GHz [1]. From
the concept of LoS MIMO [5], we have:
Lλ
,
(3)
dHAP dGS =
β
where L is the distance between an HAP and the GS, β represents the degrees of freedom. Also, dHAP and dGS , respectively,
represent the inter-HAP and inter-GS distances. To determine
the HAP separation distance, we choose the following parameters: dGS = 0.1 m, L = 18 km, λ = 0.00625 m, and β = 1.
Now, by using (3), we find that the HAPs must be spaced
1125 m apart to ensure that the channels will be uncorrelated.
In this case, each HAP drone-GS wireless communication
channel will be full rank which allows sending data over
multiple paths.
The single tethered balloon relay performs the DF operation
and operates in half-duplex mode as shown in [9]. The
communication between HAPs and GSs is carried out in two
phases, direct transmission and relay aided transmission. In the
first phase, GSs receive signal directly from the HAP drones
while the relay remain silent. In this case, HAP drones transmit
data to one specific GS during one time slot. In the second
phase, tethered balloon relay is active and transmits data to the

3

GS after precoding. Therefore, the GS receives signals from
the relay and the first transmitter.
B. Capacity of Rician X network
The asymptotic sum-rate of a network as a function of the
signal-to-noise-ratio (SNR) can be expressed as [6], [8]:
C = β · log(γ) + O(log(γ)),
(4)
where γ is the SNR value at a given receiver.
The DoF for an X network with M transmitters and N
NA
receivers each with A antennas, is equal to β = MM+N
−1 .
However, we consider a relay-aided system that uses a DF
relaying mechanism. We denote the channel matrices between
the relay and HAP i by H i , and between GS j and the relay
by Gj . Since the Rician factor of the HAP-to-relay link is
greater than that of the relay-to-GS link, we model H i and
Gj with different κ values. In this case, the small-scale fading
matrices H̄ i and Ḡj are obtained by replacing κ in (1) by κui
and κlj , respectively. After adding the path loss components
ψj
i
to H̄ i and Ḡj , we get H i = (dα
2 H̄ i and Gj = (d
2 Ḡj ,
Ri )
jR )
where αi and ψj represent, respectively, the channel gains in
H i and Gj at a 1 m reference distance. Also, dRi and djR
are the link distances in HAP i-tethered balloon and tethered
balloon-GS j communications. In order to find capacity of the
considered relay aided HAPs-GSs system, we first calculate
the SNR of each stream. Here, we have two phases: HAPstethered balloon communications, and tethered balloon-GSs
communications. We begin with the first phase, and the same
steps can be used to find the sum-rate of the second phase. The
zero-forcing detection (ZF) SNR of the k th parallel channel,
γk , is given by [13]:
Γ
γk =
,
(5)
−1
[W ]k,k
where Γ is the transmit power per symbol, and W = H H H.
Now, given the knowledge of ZF SNR, we derive the explicit
expression for ZF-capacity as follows.

γ1 = 

Es
σ 2 NT
 ,
(H H H)−1 1,1

Es
where N
is the transmitted energy per symbol, and σ 2 is the
T
noise power. From [15], we know that [W −1 ]1,1 can be found
from the elements of H. That is,
 H

H
H
−1
H̃ ]h1 , (10)
(H H)−1 1,1 = hH
1 [I − H̃(H̃ H̃)

where, H = [h1 H̃]. Now, substituting (10) in (9) leads to:
H
H
Es H
h [I − H̃(H̃ H̃)−1 H̃ ]h1 .
(11)
γ1 = 2
σ NT 1
By substituting H i for H in (11), we get φi :
H
H
EHAP H
hi1 [I − H̃ i (H̃ i H̃ i )−1 H̃ i ]hi1 ,
(12)
φi = 2
σ NT
where EHAP is the transmitted power of each HAP. Similarly,
the SNR in GS j is:
H
H
EBL H
ϕj = 2
g j1 [I − G̃j (G̃j G̃j )−1 G̃j ]g j1 ,
(13)
σ NT
where EBL is the transmitted power of tethered balloon, associated with Gj . The overall system capacity can be defined
as in [16]:
CDF (γ) = β · min(C1 (γ), C2 (γ)),
(14)
PM
where C1 = i=1 log(1 + φi ) is the sum-rate between the
PN
HAPs-tethered balloon link and C2 = j=1 log(1 + ϕj ) is
the sum-rate between the tethered balloon-GSs link. Finally,
by substituting (12) and (13) in (14) this theorem is proved.

Using Theorem 1, it is possible to analyze the impact of κui ,
dRi , and M on the capacity of the Rician X-channel in
the HAP drones wireless system.
κlj ,

Remark 1. At high SNR, the system fails to achieve a higher
sum-rate. This is because, at larger κui and κlj , the columns of
H will be correlated. Hence, IA will fail to achieve a maximum
capacity for higher κui and κlj values.
Remark 2. When M = N and κui =κlj , the capacity of the
2

Theorem 1. The capacity of a relay aided multiple HAPs-GSs
communication is given by:
M
X

MN
C=
min
log2 1 + Γs hH
i1 Qi hi1 ,
M +N −1
i=1
N
X

log2 1 + Γs g H
j1 W j g j1



, (6)

j=1

where Qi and W j are given by:
H

H

Qi = [I − H̃ i (H̃ i H̃ i )−1 H̃ i ],
W j = [I −

H
H
G̃j (G̃j G̃j )−1 G̃j ],

(7)
(8)

with H̃ i , hi1 and G̃j , g j1 obtained from channel matrices H i
and Gj based on H i = [hi1 H̃ i ] and Gj = [g j1 G̃j ].
Proof. In the proposed model, we have two phases: HAPstethered balloon and tethered balloon-GSs communications.
We first find the SNR of HAP i-to-relay link, γi , and relayto-GS j, γj . Then, we proceed to calculate the sum-rate. In
general, the SNR of the first stream at the receiver is given
by:

(9)

E ×(d

)2

Ri )
= BL ψj jR . Also,
system is maximum when EHAP ×(d
αi
when the transmitted power EHAP = EBL , the maximum
capacity is obtained when the tethered balloon is at the center
of the HAPs-GSs link.

IV. N UMERICAL R ESULTS
Here, we evaluate the sum-rate performance of the relay
aided multiple HAPs-to-GSs communications. HAPs are located at an altitude of 18 km above earth and each GS receives
signals from all HAPs. We consider an X channel, consisting
of M = 2 and 3 HAPs, a tethered balloon, and 3 ground
stations. Fig. 2 shows the sum-rate of the Rician channel as
a function SNR for different Rician factors and number of
HAPs. As we can see from this figure, the sum-rate obtained
using simulations (the curves shown without a circle marker)
is in agreement with the analytical derivation. Without the use
of a tethered balloon, the sum-rate of the considered system
decreases as the maximum DoF cannot be achieved. While
exploiting IA, however, the system is equivalent to a system
in which CSI is known at the transmitters and, hence, a higher

4

5

u
l
dopt
RD = 15.5 km. Also, when κ = 30 dB and κ = 20 dB,
the optimal relay’s altitude is 14 km. In fact, as the difference
between κu and κl become smaller, dopt
RD converges towards
dSD
2 .

2 HAP (Simulation) , κ l = 15 dB
2 HAP (Theory), κ l = 15 dB

Sum−rate (bps/HZ)

2 HAP (Simulation) , κ l = 20 dB

4

2 HAP (Theory), κ l = 20 dB
3 HAP (Simulation) , κ l = 15 dB

3

3 HAP (Theory), κ l = 15 dB

V. C ONCLUSION

3 HAP (Simulation), κ l = 20 dB
3 HAP (Theory), κ l = 20 dB

2

2 HAP (non IA), κ l = 20 dB

1
0
0

10

20

30

SNR (dB)

Sum−rate (bps/HZ)

Fig. 2: Sum-rate vs. SNR.

5

Pt = 9.5 dB, κu = 30 dB, κl = 15 dB
Pt = 9.5 dB, κu = 30 dB, κl = 20 dB

4

Pt = 14.5 dB, κu = 30 dB, κl = 15 dB
Pt = 14.5 dB, κu = 30 dB, κl = 20 dB

R EFERENCES

3
2
1
0
0

2

4

6

8

10

12

14

In this paper, we have proposed an effective interference
alignment scheme for maximizing the sum-rate of HAPsground stations communications assisted by a tethered balloon relay. In particular, we have considered the half-duplex
relaying scheme using a tethered balloon relay to achieve
maximum DoF in HAPs-ground stations communications,
when the HAPs lack the knowledge of CSI. Our results have
shown that, using a tethered balloon relay for exploiting IA
provides a significant sum-rate gain in the HAP-based wireless
system that uses multiple interfering drones.

16

18

Altitude of tethered balloon (km)

Fig. 3: Sum-rate vs. tethered balloon’s altitude.
sum-rate is achieved. In this analysis, we consider κu = 30 dB
and dSR = 1 km and we vary κl and M . As expected, the
sum-rate increases by using the IA solution and decreasing the
Rician factor, κl . For instance, at SNR = 25 dB, the channel
with κl = 20 dB offers a sum-rate gain of 34.8%. Interestingly,
under the same SNR, the sum-rate increases by up to 1.7
times when κl = 15 dB. From Fig. 2, we can also see that,
as κ increases, the sum-rate degradation occurs at a higher
SNR. This is due to the fact that increasing κ increases the
correlation between channels, which results in a degradation
of the asymptotic performance.
From Fig. 3, we can see that sum-rate increases as the
tethered balloon moves away from the HAPs. Similarly, the
sum-rate decreases with κl , provided that the HAPs are located
sufficiently apart. This is due to the fact that exploiting IA
while deploying a tethered balloon relay allows achieving
a maximum DoF in the system. Fig. 3 shows that the optimal altitude of tethered balloon at which the HAPs-tethered
balloon-GSs link has the maximum sum-rate (dopt
RD ) does not
change by varying the transmit power. As we can see from
this figure, given dSD = 18 km, a channel with Pt = 14.5 dB,
κu = 30 dB and κl = 15 dB offers maximum sum-rate, when

[1] A. Mohammed, A. Mehmood, F. N. Pavlidou, and M. Mohorcic,
“The role of high-altitude platforms (HAPs) in the global wireless
connectivity,” Proceedings of the IEEE, vol. 99, no. 11, pp. 1939–1953,
Nov. 2011.
[2] M. Mozaffari, W. Saad, M. Bennis, and M. Debbah, “Unmanned aerial
vehicle with underlaid device-to-device communications: Performance
and tradeoffs,” IEEE Transactions on Wireless Communications, vol. 15,
no. 6, pp. 3949–3963, June 2016.
[3] M. Mozaffari, W. Saad, M. Bennis, and M. Debbah, “Efficient deployment of multiple unmanned aerial vehicles for optimal wireless
coverage,” IEEE Communications Letters, vol. 20, no. 8, pp. 1647–1650,
Aug. 2016.
[4] Y. Liu, D. Grace, and P. D. Mitchell, “Exploiting platform diversity
for GoS improvement for users with different high altitude platform
availability,” IEEE Transactions on Wireless Communications, vol. 8,
no. 1, pp. 196–203, Jan. 2009.
[5] A. Irish, F. Quitin, U. Madhow, and M. Rodwell, “Sidestepping the
Rayleigh limit for LoS spatial multiplexing: A distributed architecture
for long-range wireless fiber,” in Proc. of ITA, San Diego, CA, USA,
Feb. 2013, pp. 1–6.
[6] V. R. Cadambe and S. A. Jafar, “Interference alignment and the degrees
of freedom of wireless x networks,” IEEE Transactions on Information
Theory, vol. 55, no. 9, pp. 3893–3908, Sept. 2009.
[7] F. Pantisano, M. Bennis, W. Saad, M. Debbah, and M. Latva-aho,
“Interference alignment for cooperative femtocell networks: A gametheoretic approach,” IEEE Transactions on Mobile Computing, vol. 12,
no. 11, pp. 2233–2246, Nov. 2013.
[8] V. R. Cadambe and S. A. Jafar, “Interference alignment and degrees
of freedom of the K-user interference channel,” IEEE Transactions on
Information Theory, vol. 54, no. 8, pp. 3425–3441, Aug. 2008.
[9] Y. Tian and A. Yener, “Guiding blind transmitters: Degrees of freedom
optimal interference alignment using relays,” IEEE Transactions on
Information Theory, vol. 59, no. 8, pp. 4819–4832, Aug. 2013.
[10] H. Zebardast, A. Tajer, B. Maham, and M. Rezaee, “Relay X channels
without channel state information at the transmit sides: Degrees of
freedom,” in Proc. of WCNC, New Orleans, LA, March 2015, pp. 493–
497.
[11] A. Zajić, Mobile-to-mobile wireless channels. Artech House, 2012.
[12] M. Mozaffari, W. Saad, M. Bennis, and M. Debbah, “Optimal transport
theory for cell association in UAV-enabled cellular networks,” IEEE
Communications Letters, vol. 21, no. 9, pp. 2053–2056, Sep. 2017.
[13] A. Paulraj, R. Nabar, and D. Gore, Introduction to space-time wireless
communications. Cambridge university press, 2003.
[14] Y. S. Cho, J. Kim, W. Y. Yang, and C. G. Kang, MIMO-OFDM wireless
communications with MATLAB. Wiley & Sons, 2010.
[15] M. Kiessling and J. Speidel, “Analytical performance of mimo zeroforcing receivers in correlated rayleigh fading environments,” in Proc.
of SPAWC. IEEE, 2003, pp. 383–387.
[16] G. Levin and S. Loyka, “Amplify-and-forward versus decode-andforward relaying: Which is better?” 2012.

