1. Introduction
The explosive number of devices demands data traffic growth and new radio spectrum resources in future 6G systems and the IoT context. Key enabling technologies are (i) the underutilized millimeter wave (mmWave) band (between 30 GHz and 300 GHz), which could be valuable for its wide bandwidth and higher spectral efficiency, and (ii) NOMA, which could be valuable for simultaneously supporting multiple users on the same time-frequency resources. The joint use of the large spectrum available in mmWaves together with massive multiple-input-multiple-output (MIMO) strategies allows ultra-high data rates to be guaranteed through spatial directional transmissions compensating for the high propagation loss of mmWaves communications [
1]. This directional nature of mmWaves transmissions needs the support of one radio frequency (RF) chain for each user on the same time-frequency resource. Therefore, the hardware complexity and costs of the mmWave MIMO system increase with an increase in the number of users. Hybrid architectures, which combine phase shifters based on analog precoding and digital precoding, reduce the costs practically by selecting a reduced number of RF chains. However, even if hybrid beamforming structures are implemented, the user’s channels are highly correlated in mmWave communications, and thus, the users cannot be separated by linear operations. Such a correlation facilitates integration with NOMA technology. NOMA can simultaneously serve multiple users on the same time-frequency resource by converting their channel gains into multiplexed gains in the power domain and by using successive interference cancellation (SIC) at the receiver to remove intra-channel interference with a decoding order based on the channel conditions [
2,
3,
4].
In the literature, different papers analyze the spectral efficiency (SE) and the energy efficiency (EE) maximization problem mainly for a downlink mmWave MIMO with hybrid beamforming [
5,
6]. In [
3], a power allocation (PA) problem to optimize EE is considered for an uplink NOMA-assisted mmWave MIMO system under users’ quality-of-service (QoS) and quality-of-experience requirements [
7]. In [
4], power allocation and beamforming are jointly considered to maximize the sum rate of a pair of users in a mmWave NOMA system by using an analog beamforming structure indeed of a hybrid mmWave beamforming structure.
Simultaneous wireless information and power transfer (SWIPT) techniques are integrated in mmWave massive MIMO-NOMA systems to maximize the energy efficiency in [
8], and consequently, each user can extract both information and energy from the received RF signals by using a power splitting receiver.
In the 6G vision, low Earth orbit (LEO) satellites and diverse aerial platforms, such as UAVs, are considered to support IoT development in remote areas and in emergency situations thanks to the mobility, flexibility and good channel conditions of UAVs. UAVs have been developed for their monitoring capabilities, implemented by on board sensors, for services in agriculture or security border controls. In addition, UAVs are used as a flying base station (BS) to provide ubiquitous connectivity and effectively increase the coverage and throughput of wireless systems through the optimization of UAVs’ positions and trajectories [
9].
Recently, UAVs have been considered to act as relays for cooperative communications due to the high probability to establish LoS links with the ground terminals. UAV-assisted relaying systems operate according to the two classical types of transmission protocols, namely decode-and-forwarding (DF) and amplify-and-forwarding (AF) [
10].
Several papers consider UAV communications combined with orthogonal multiple access (OMA) to maximize the energy efficiency or throughput by optimizing the source/relay transmit power and the UAV speed and trajectory design as, e.g., in [
11]. The paper [
12] focuses on a UAV full duplex (FD) relay with joint beamforming and power allocation to optimize the instantaneous data rate when the UAV flight follows a circular trajectory. FD relaying allows for a relay node to simultaneously transmit and receive in the same band, unlike half-duplex (HD) mode [
10,
13]. Therefore, a natural choice is to combine FD relays with NOMA to enhance spectral efficiency.
A UAV-supported clustered-NOMA system for the 6G-enabled IoT is detailed in [
14], where the numerous terminals are partitioned into clusters and the UAV provides services to the clusters by using wireless-powered communication (WPC) to optimize the uplink average achievable sum rate of all terminals by designing the UAVs’ trajectory.
In [
15], a multiple-UAV-aided NOMA scheme is proposed to improve spectral and energy efficiency for cellular uplinks. In particular, half of users are partitioned in clusters served by multiple UAV relays, and the other ones communicate with the BS directly. A location-based user pairing (UP) scheme associates the clustered users with the multiple-UAV-aided NOMA to minimize the resource allocation problem.
However, all the aforementioned works focus on UAV-relay-aided NOMA without considering mmWave communications. The severe signal power attenuation in the mmWave band impacts UAV connectivity performance, especially when very long communications distances exist between the ground users and the associated UAV which serves them.
Therefore, this paper addresses the design of a UAV-enabled FD relaying network in the mmWave band to aid an uplink NOMA cellular system.
In detail, a multiple-UAV-relaying network supports NOMA technology, and hybrid mmWave beamforming is considered at the base station (BS), which can estimate the DoA information of UAVs to improve the overall sum rate of the system. Indeed, the DoAs are unknown and the MUltiple SIgnal Classification (MUSIC) method [
16] considered in this paper estimates directly the DoAs at the hybrid mmWave beamformer.
In [
17], a beamspace MUSIC algorithm is used to estimate path directions for mmWave channel estimation problem showing that the hybrid precoding structure can avoid the spectrum ambiguity and maximize the number of resolvable path directions. The mmWave channel estimation problem is also considered in [
18] by using MUSIC for the hybrid analog/digital beamformer with 2D co-prime arrays where the directions can be uniquely estimated by finding the common peaks of the 2 decomposed subarrays.
Deep-learning approaches are considered to evaluate the angle-of-arrival (AoA) information in the uplink of an mmWave communication system based on MUSIC to enhance classification accuracy in [
19]. The above papers use the MUSIC algorithm to derive the DoA information in mmWave hybrid beamforming, but they do not deploy an UAV-relay network, as in this paper.
To the best of our knowledge, this paper is the first contribution that considers jointly (i) the UL mmWave communications, (ii) the hybrid beamforming with DoA estimation based on the MUSIC method and (iii) multiple UAVs acting as an aerial BS to relay NOMA transmissions and, consequently, achieve better data rate for users who can suffer severe channel conditions. Moreover, a novel maximization design of the overall sum-rate is proposed for the uplink mmWave transmission of the multiple-UAV-relay-aided NOMA 6G system.
The contributions of this paper can be summarized as follows:
We consider a multiple-UAV-aided NOMA network where each UAV acts as a mobile FD relay in mmWave UL for 6G cellular systems;
We propose the use of the MUSIC algorithm at the hybrid beamforming to detect the DoA estimations for each UAV and improve the performance of the UAV-aided NOMA cellular system;
We propose an optimization procedure in order to maximize the average achievable sum rate of the UL mmWave UAV relay network by taking into accounts LoS obstruction, channel time-varying condition, the DoA information at the hybrid beamforming due to the different spatial directions of UAVs, as well as the requirements in terms of quality of service (QoS) for each user.
The paper is organized as follows.
Section 2 presents the proposed system model, channel model, the mmWave hybrid beamforming, the MUSIC algorithm and the problem formulation to optimize the global UL sum rate of the UAV-network-aided NOMA cellular system. In
Section 3, the numerical results providing a comparison with the NOMA cellular system without the use of multi-UAV relays are shown. Finally,
Section 4 concludes the paper and outlines future research activities.
2. System Model and Problem Formulation
A UAV-enabled full-duplex relaying system is considered to aid an uplink mmWave NOMA cellular system consisting in a BS, with
N users and
K UAVs, as shown in
Figure 1. The
N users are randomly distributed in the cell of radius
R, and the BS is located at the center of the cell.
We assume that the BS is unable to deliver the superimposed signals to the NOMA users in a far subarea of the cell because the link between users and BS is negligible due to severe blockage.
Each UAV acts as a DF relay to help data transmission between the BS and users and operates in FD mode.
We suppose that
K UAVs are flying at height
in such a way that all links from UAVs to BS are LOS channels. The elevation angles between the BS and the UAVs are denoted as
, and
is the elevation angle between the
n-th user (named as user equipment UE
) and UAV
, as shown in
Figure 1.
2.1. Path Loss Model
The links from ground users and UAV are LOS or no LOS (NLOS) channels due to the presence of buildings, vegetation etc., which can obstruct the signals propagation. In detail, the channel gain between the
k-th UAV and the BS is, according to [
20,
21,
22],
and if the
is in LOS with the
user, the channel gain is defined as
where
and
represent the distances between UAV
and BS and between UAV
and UE
, respectively.
denotes the additional attenuation factor of the LOS channel at the reference distance
m, and
is the path loss exponent at the air channel. UAVs are in
with probability
, where
B and
C are constants related to the environment, whereas UAVs are in NLOS according to the complementary probability
.
2.2. MmWave Hybrid Beamforming
In mmWave communications, large phase arrays are usually adopted to overcome the high propagation losses, and in combining with NOMA, mmWave beamforming is used to increase beam gain and serve multiple users. Usually, hybrid analog/digital beamforming is adopted in mmWave NOMA communications, where the precoding is performed in hybrid mode by combining the digital baseband precoding with an analog RF beamforming driven by a limited number of RF chains. This hybrid analog/digital beamforming is a cost-effective solution due to the use of massive antennas with limited RF chains. It can be easily implemented through the use of analog phase shifters together with the abilities of digital precoding, which allow the beams to be directed towards the desired user and remove inter-user interferences [
23,
24].
In particular, we consider the hybrid beamforming architecture at the BS with a number
of RF chains exploiting NOMA in each RF chain and spatial division multiple access (SDMA) between RF chains [
8], as shown in
Figure 2.
From the angle domain perspective, the knowledge of DoAs plays a fundamental role.
In this paper, the UAVs’ angle information are discovered by using the MUSIC algorithm implemented at the hybrid mmWave beamforming. We assume that the number of UAVs does not exceed the number of available RF chains, i.e., .
The received signal at the BS after RF and digital beamforming combining can be expressed as [
19,
24]
where
is the channel matrix (detailed in Equation (
6)),
is the number of receiving antennas at the
,
denotes the transmit signal,
is the zero-mean independent and identically distributed (i.i.d.) Gaussian white noise vector with power
and
,
are the analog and digital beamforming matrices (see Equations (
4) and (
5)), respectively, with
.
In partially-connected hybrid mmWave beamforming, the antenna array with
elements can be organized into groups, called subarrays, and each subarray connected to one RF chain processes each received signal by a phase shifter. Then, all of them are added up as shown in
Figure 2.
Therefore, the analog matrix
of Equation (
3) is a diagonal matrix where
is the vector of subarray
j with
;
is chosen from a uniform distribution in the range of
. The RF signal passes through
parallel RF chains. It is down-converted and then the digital beamforming operation follows. The zero-forcing scheme can be used as a digital beamformer as
where
considers the actual channel matrix
filtered by the analog beamforming matrix
.
2.3. Channel Model
We consider a ray-tracing channel model widely used in mm-Wave communications with a limited number of L scattering paths, as, due to the spatial sparsity in the mm-wave channel, it is expected that the propagation paths are along a small number of directions [
19,
25,
26].
Therefore, each column of the channel matrix
is defined as
where
is the complex gain of the i-path due to small-scale fading,
is the path loss between the
and the UAV
and
is the angle of arrival of the
path of the
UAV.
Without loss of generality, we assume a uniform linear array (ULA) at the
for simplicity, and the steering vector can be expressed as:
where
is the signal wavelength, and
l is the distance among antenna elements, where
.
In the following, we assume that in Equation (
6) the variations of the channel are only caused by the path gains
, and the path angles remain unchanged according to the mmWave channel measurements in [
27].
2.4. MUSIC Algorithm
The DoA information is estimated by MUSIC spectral estimation on the filtered version of the received signal, as shown in
Figure 2, [
19,
25,
26]. The MUSIC algorithm developed by Schmidt [
16] is an eigenstructure-based DOA-finding method and, similar to other parametric algorithms such as ESPRIT, has demonstrated a superior resolution with respect to the non-parametric methods.
By considering a partially connected hybrid mmWave beamforming, after the analog beamforming and the analog-to-digital conversion, the baseband signal is [
19,
28]:
By performing eigenvalue decomposition on the covariance matrix
of the output vector
of the virtual array, we have:
where
is the
matrix of the signal eigenvectors corresponding to the
K largest eigenvalues,
is the
noise subspace matrix with eigenvectors corresponding to the smallest
singular values, and
is the
diagonal matrix containing the
K largest eigenvalues
of
. It is clear that the signal subspace
and the noise subspace
are orthogonal.
The MUSIC algorithm utilizes the orthogonality between the two complementary spaces to estimate the spatial signal. Therefore, the DoA estimation
consists of finding the values of
, whereby the filtered vector
is related to the signal subspace of
, where
is the array manifold vector of the virtual array.
By using the definition of a pseudo-spectrum of the MUSIC algorithm, the estimated DoA of the emitter direction can be calculated by maximizing the function
which provides high resolution of angle separation.
2.5. Problem Formulation
Considering the NOMA method in each beam, intra-beam superposition coding at the transmitter and SIC at the receiver are performed. In the case of uplink mmWave NOMA, the users begin to transmit uplink signals
and
at the same time and in the same frequency band. For the proposed multiple-UAV relay network shown in
Figure 1, each UAV decodes the mixed signals of the two users in its beam (called, for example, UE
and UE
) and then transmits the superimposed signal
to the BS according to
where
and
are the decoded signals of UE
and UE
, respectively, and
and
cannot exceed the maximum transmission power of the UAV.
In the conventional NOMA with a single-antenna, usually the information of the user with a lower channel gain is decoded first to maximize the sum rate. In contrast, in mmWave NOMA, the decoding order depends on both channel gain and beamforming gain. Without loss of generality, assuming that UE
has a better channel condition with respect to UE
in the area covered by UAV
[
4,
29], the achievable rates are
and the achievable rate from
is
where
,
and
are the transmission power of UE
, UE
and UAV
respectively,
,
are the channel response vector between the UE
, UE
and UAV
,
is the channel response vector between the UAV
and the BS,
and
represent jointly the analog and digital precoding matrices at the UAV
and at the BS, whereas
and
are the power of the zero-mean additive Gaussian white noise at the UAV
and BS, respectively.
Therefore, the available rate received at the BS considering the UAV
acting as a relay is
To maximize the uplink average achievable sum rate of all terminals of the multiple-UAV-aided NOMA mmWave system by dynamically tracking the DoAs of multiple UAVs, the optimization problem can be formulated as
where
denotes the minimal data rate constraint for each user and the last constraint indicates the transmitted power constraint with P being the maximum total transmitted power.
This optimization problem is very complicated to be solved directly because the problem is non-convex and may not be converted to a convex problem with simple manipulations. Consequently, in order to validate the effectiveness of the proposed optimization problem and the accuracy of the derived analytical model, we resort to numerical simulations.
3. Simulation Results
In this section, we evaluate the performance of the proposed network of multiple UAV relays supporting the uplink mmWave NOMA system. In the considered system, UAV relays are located in different positions with different spatial directions with respect to the BS, and hybrid mmWave beamforming with MUSIC technique is adopted to derive the different DoAs in order to optimize the achievable uplink sum-rate. In detail, we simulate three different scenarios, i.e., an UL mmWave NOMA without a UAV network, a single UAV and a network of multiple UAVs supporting the UL mmWave NOMA. The cellular network has an area of m with the BS located in the center. A total of 32 users are randomly and uniformly distributed in this area.
The BS is equipped with an ULA of antennas and RF chains to simultaneously serve a number K of UAVs with . The UAVs deployed in the cell have different DoAs at the BS uniformly distributed between , and each UAV covers a pair of users who have poor connections with the BS due to a severe blockage condition. The UAV has receiving antennas which transmits towards the BS with one antenna.
We choose to associate the considered blocked user pair not to the closest UAV in term of distance, but to the one offering the best communication performance, as especially in a very dense urban scenario, the nearest available UAV may be in NLOS visibility. An accurate DoA estimation at the BS can monitor the fluctuations of the UAV motion and consequently improve the performance of the overall uplink sum rate.
The main parameters for an urban scenario are summarized in
Table 1 [
20,
27].
In
Figure 3, the achievable sum rate is shown in terms of the varying signal-to-noise ratio (SNR) for the multi-UAV network with respect to a single UAV scheme in the case of an LOS environment. The multi-UAV scheme outperforms the single UAV scheme even for low SNR values, as the channels from users to each UAV are in LOS in these simulations realizations and consequently demand low power to transmit towards the UAV.
In the same figure, the case of a multi-UAV scheme without DoA information acquisition is highlighted to validate the performance of the MUSIC method. Indeed, the DoA estimation error is an important problem due to the UAV mobility and dynamic channel variations. However, a more accurate DoA estimation occurs when the UAV is far from the BS [
30], and in particular the use of the MUSIC technique allows both the DoA and received powers to be estimated more accurately with respect to other methods [
19,
24,
31].
Figure 4 shows the achievable sum rate results with respect to SNR when multiple UAVs are considered for both the case of LOS and NLOS environments. The MUSIC method at the BS to acquire the DoA information of the UAVs is used. In these simulation results, the heights of UAVs are not considered too high, because even if the links have a higher probability to be in LOS, the impact of the increased distance between the pair of users and the serving UAV decreases the overall link budget, and the links can have a worse SNR. The results are compared with the performance of the UL NOMA mmWave system without the multi-UAV network, i.e, the users can communicate directly by NOMA mmWave link to the BS. The NOMA mmWave link to BS is considered only for evaluating the performance comparison as, practically, the user pair covered by each UAV has a severe blockage and cannot communicate with the BS. The achievable sum rate of the multiple-UAV network achieves a considerably better performance than the direct NOMA mmWave connection, as shown in
Figure 4, providing a improvement of about 6 dB for high SNR values.