Rate Maximization in a UAV Based Full-Duplex Multi-User Communication Network Using Multi-Objective Optimization

: In this paper, we study an unmanned-aerial-vehicle (UAV) based full-duplex (FD) multi-user communication network, where a UAV is deployed as a multiple-input–multiple-output (MIMO) FD base station (BS) to serve multiple FD users on the ground. We propose a multi-objective optimization framework which considers two desirable objective functions, namely sum uplink (UL) rate maximization and sum downlink (DL) rate maximization while providing quality-of-service to all the users in the communication network. A novel resource allocation multi-objective-optimization-problem (MOOP) is designed which optimizes the downlink beamformer, the beamwidth angle, and the 3D position of the UAV, and also the UL power of the FD users. The formulated MOOP is a non-convex problem which is generally intractable. To handle the MOOP, a weighted Tchebycheff method is proposed, which converts the problem to the single-objective-optimization-problem (SOOP). Further, an alternative optimization approach is used, where SOOP is converted in to multiple sub-problems and optimization variables are operated alternatively. The numerical results show a trade-off region between sum UL and sum DL rate, and also validate that the considered FD system provides substantial improvement over traditional HD systems.


Introduction
In the past few years, Unmanned Aerial Vehicle (UAV) communication systems have become a focus of researchers for their unique characteristics that produce interesting communications problems. In particular, UAVs are expected to be integral in various applications in 5G wireless systems. Due to the fact that UAVs are easily operable, highly maneuverable, and have increasing payload capability, it is expected that the use of UAVs will increase in the coming years. There are numerous applications of UAVs including emergency search and rescue, communication relaying, package delivery, emergency broadband service, and infrastructure inspection [1].
In the telecommunication domain, a common use of the UAV is a mobile base station (BS) to form a wireless network [2]. In [3], the authors worked on multi-UAV communication systems to serve ground users and primarily focused on designing the UAV trajectories and transmit power to improve throughput. UAVs can also be utilized to work as a relay. The authors in [4] worked on UAV-aided relay systems to minimize outage of the system by optimizing trajectory and transmit power of the UAV. In [5], authors studied joint optimization of altitude and beamwidth of UAV to improve throughput in three different UAV enabled multi-user communication models, where the UAV was assumed to have directional antennas with adjustable beamwidth angles. The results show that the optimal altitude and beamwidth angle critically depend on the communication model. Furthermore, the antenna beamwidth optimization was also considered in [6].
Concurrently, the increase in wireless data traffic forced the limits of communication systems in terms of reliability and throughput [7]. In band full duplex (IBFD) communication allows bi-directional communication at the same frequency and time, which doubles the system capacity and improve spectrum efficiency compared to the traditional half duplex (HD) systems [8][9][10]. The IBFD systems are of great interest among researchers, and there has been some work in the literature where IBFD systems are utilized with UAVs to propose and solve new problems. For example, in [11], the UAV was deployed as a full duplex (FD) relay to minimize the outage probability of the relaying system. In [12], spectrum sharing planning was studied in a scenario where FD UAV was considered as a relay in device-to-device cellular systems. In [13], FD UAV-based small cell wireless system was considered, where the UAV was deployed as a base station to serve downlink (DL) and uplink (UL) users in a FD manner. The UAV trajectory, user scheduling, and UL user transmit power are alternately optimized to maximize the total capacity of the system.
In addition to the spectral efficiency, beamforming with MIMO system has gained wide research attention for the design of a next generation wireless network. Beamforming is a technique used for directional transmission and reception at a particular device/user, rather than spreading the wireless signal in all directions, thereby reducing interference and increasing the wireless link gain. In [14], transmit and receive beamforming optimization for a IBFD cloud radio access network (CRAN) was considered. Similarly, beamforming optimization for power over wireless FD MIMO systems was discussed in [15], where considerable performance gains were achieved over the HD scheme when the time splitting parameter and beamformer are jointly optimized. The UAV enabled FD relay was considered in [16], where beamforming and power allocation were jointly optimized to maximize the instantaneous data rate of the system. In [17], a sum-rate maximization problem for FD multi-user MIMO system was considered, where the authors proposed a beamformer design for the UL and DL FD link. Similarly, in [18], a power efficient resource allocation algorithm was considered in a multi-user FD network. The authors proposed a multi-objective optimization framework to jointly optimize the beamformer at the MIMO FD base station, and the UL transmit power from the HD users.
In recent works, effectiveness of FD in a UAV network was shown in [19], where the authors considered maximizing the DL sum rate while maintaining the minimum quality-ofservice (QoS) for the UL users, by optimizing the resources and location of the UAV. In [20], FD was considered in a UAV-assisted wireless powered Internet-of-Things (IoT) network, where the FD UAV collects data from the target IoT device and simultaneously charges remaining devices utilizing fly-hover-communicate protocol. Moreover, a FD UAV relay network was considered in [21], to jointly optimize the HD UL and DL users scheduling, and the UAV trajectory. Finally, in [22], a UAV was considered to charge nodes using Wireless powered communication via backscatter and/or harvest-then-transmit in the IoT network. The UAV trajectory with other resource allocations were optimized to maximize the energy efficiency of the network. The UAV serves all the nodes by sequentially visiting each node in respective time slot which leads to increase in energy consumption due to propulsion. Noting that the propulsion energy consumption during hover is usually less than the frequent maneuverings [23]. In this work, we consider optimal UAV deployment with resource allocation to serve multiple users on the ground.
To the best of our knowledge, UAV coverage area optimization with resource allocation was never considered in an MIMO UAV FD network. Coverage area optimization has been considered for traditional HD single antenna system in the literature [5,24], where antenna beamwidth and altitude of the UAV were optimized to maximize the system performance. For FD system, adjustable beamwidth antenna was also considered in [19] for directional communication from UAV. Here, we consider coverage optimization by including adjustable beamwidth antennas in an MIMO FD UAV BS and formulate bi-direction FD rates where interference terms are shown to be a function of antenna beamwidth. In the following, we propose a novel multi-objective-optimization-problem (MOOP) to maximize the sum UL-DL rate while providing QoS for a multi-user MIMO FD UAV network. The proposed formulation takes into account the coverage area (beamwidth and altitude), beamformer, horizontal location of the UAV, and UL transmit power of the FD users. The optimization problem is non-convex and non-tractable. MOOP was also considered in [20], where a deep learning approach is used to solve the conflicting system objectives. Contrary, we first adopt weighted Tchebycheff method to convert MOOP to a single-objective-optimization-problem (SOOP). SOOP are singular weighted/fractional order objective systems which are discussed in [25,26]. Second, we divide SOOP into multiple sub-problems and provide an independent solution for each sub-problem. Lastly, we propose a joint optimization algorithm to solve the proposed problem.
The main contribution of the paper are summarized as follows: • We consider adjustable beamwidth antennas on the FD MIMO UAV in a multi-user aerial network where the users have FD capability and propose a novel optimization problem to optimize the UAV's coverage area with resource allocation; • We propose a computationally-efficient solution for transmit beamforming from FD MIMO UAV towards the single antenna FD users by utilizing a zero-forcing (ZF) approach; • We found that the UL power optimization in a multi-user FD aerial network leads to a non-convexity in the optimization problem due to co-channel interference; • We proposed a joint optimization solution for the optimization of DL beamformer, coverage area and location of the UAV, and UL power of the users in a multi-user FD aerial network when the objective is to maximize the sum bi-directional FD rate; • Numerical results show that the proposed solution to the MOOP achieves better performance compared to the baseline FD algorithms. Similarly, FD functionality in the proposed system model provides notable performance improvement over a HD system.
The rest of the paper is organized as follows. In Section 2, the system model of a UAVbased FD multi-user communication network is presented, including the channel model and the relevant performance metrics. The proposed MOOP formulation is described in Section 3, where the sum UL and DL rates are maximized. The optimization solution, proposed joint resource allocation algorithm, and complexity and convergence analysis of the algorithm are investigated in Section 4. The numerical results based performance analysis of the proposed algorithm is presented in Section 5. Finally, the paper is concluded in Section 6.
Notation: The bold face lower case and capital case letters are used to denote vectors and matrices, respectively. X H , Trace(X), Rank(X), and X −1 indicate the Hermitian transpose, trace, rank, and inverse of the matrix X, respectively. X 0 shows that X is a positive definite matrix. ||.|| denotes the Euclidean vector norm. C M×M and I M denotes the space of M × M matrices with complex entries and M × M identity matrix, respectively. E(.) denotes the expectation with respect to random variable x, and CN (µ, σ 2 ) denotes the circularly-symmetric complex Gaussian distribution with mean µ and variance σ 2 .

System Model
In this section, we present the UAV based multi-user FD communication network. We then introduce a multi-user MIMO FD communication model and performance metrics for evaluating our algorithms.

System Model
Consider a UAV-based FD communication network as depicted in Figure 1. The system consists of a UAV-based FD radio BS and N FD users. The UAV BS is equipped with M T antennas, while each user n ∈ {1, 2, . . . , N} is equipped with a single antenna for simultaneous UL and DL transmission in the same frequency band. An example of such a feasible solution is a circulator-based FD radio prototype at both BS and users which can achieve bi-directional transmission on the same antenna while providing low hardware complexity [27]. A UAV is deployed with the altitude H and the horizontal location Q = (Q x , Q y ) and the location of users are denoted by O n = (O x n , O y n ). The UAV is equipped with adjustable beamwidth antennas to modify the coverage area for the users. Since the MIMO antennas on the commercial UAVs are usually close to each other, we assume that the azimuth and the elevation half-power beamwidth of all the M T antennas are equal and defined by 2Φ ∈ (0, π) radians. The antenna gain in the direction (θ, Ψ) can be determined by: Here, θ and Ψ are the azimuth angle and the elevation angle, respectively. In (1), the former interval corresponds to the actual beamwidth of the antenna, and latter is the range of θ and Ψ values outside the beamwidth. Additionally, G o = 30,000 2 2 × ( π 180 ) 2 ≈ 2.2846, and g is the channel gain outside the beamwidth of the antenna [5]. To simplify, we set g = 0.
In contrast, each user n is equipped with an omni-directional antenna. Therefore, radius of the ground coverage area of the UAV is defined by r c = H tan(Φ).

Channel Model and Performance Metrics
In this work, we assume that the UL and DL channels between the UAV and user n have similar channel characteristics. Since a UAV to user channel is mostly dominated by the line-of-sight path, the UL and DL channels can be written as Here, h n and h n ∈ C (M T ×1) . Each element of h n is a random variable with zero mean and unit variance, and is the path loss exponent which is taken as = 2. The received signal at the FD BS and FD user n can be written as (3) and (4), respectively.
Here, S a n ∈ C and κ n ∈ C M T ×1 are the information bearing signal and transmit beamforming vector from the FD BS to user n, respectively. S b n ∈ C and p n are the UL transmit data and power sent from user n to the FD BS, respectively. g {m,n}= go d /2 {m,n} ∈ C denotes the channel between user m ∈ {1, 2, . . . , N} and user n, where g o is the channel gain at reference distance of 1 meter, d {m,n} is the distance between user m and n, and m = n. g {n,n} ∈ C and h D ∈ C M T ×M T are the self-interference (SI) channels at user n and the FD BS, respectively. π D ∼ CN (0, σ 2 D I M T ) and π n ∼ CN (0, σ 2 n ) denote the additive white Gaussian noise at the FD BS and user n, respectively. In (3), the term S a n represents the SI at the FD BS. Without loss of generality, it is assumed that E(|S a denotes the co-channel interference formed by UL transmissions from the users, the term g {n,n} S b n √ p n represents the SI at the user n, and the term shows the multi-user interference in the downlink channel. Additionally, we assume that the perfect channel state information for all users is known at the FD BS.
The FD BS applies a beamformer to the received signal y UL . The received signal after beamforming for the information transmitted from user n can be written as Here, ω n ∈ C M T ×1 is the receive beamforming vector at the FD BS to decode information transmitted from the user n. We adopt ZF beamforming for ω n for two main reasons: firstly, it requires low computation resources compared to the MMSE beamforming (Basic analysis and performance comparison of ZF and MMSE beamforming are given in [28]) and, therefore, it can be implemented more efficiently for UAV network that has limited computational resources, and secondly, the optimization problem is mathematically tractable due to the cancellation of interference terms in uplink rate expression as shown below. It is also reported that ZF beamforming approaches the performance of optimum MMSE beamforming when the noise is not a dominating factor or the number of antennas is sufficiently large in the system [29,30]. Additionally, the BS is usually equipped with a low noise amplifier, which facilitates in reduction in the noise at the BS. The receive beamforming vector for user n is calculated at the FD BS according to: Here, v n = [ . Additionally, to facilitate uplink signal detection, the number of antennas on the FD BS is assumed to be higher than the number of users.
The DL rate at user n for the information transmitted from FD BS can be written as where W is the system bandwidth. The UL rate at FD BS for the information transmitted from user n can be written as

Proposed MOOP Formulation
In this work, we study the maximization of two objective functions and formulate a MOOP. We aim to maximize the sum UL and DL rates while maintaining minimum QoS to the ground users by jointly optimizing beamwidth Φ, horizontal location Q and altitude H of the UAV, and UL power of the users. The proposed MOOP is formulated as follows: Here, P D and P B are the maximum DL power for the FD BS and the maximum UL power for the users, respectively. [H min , H max ] is the range of the UAV altitude H, which is determined by local authority regulations. R UL min n and R DL min n are the minimum required data rate for QoS in the UL and DL, respectively. [Φ min , Φ max ] is the half-beamwidth range, which is determined by an antenna beamwidth tuning technique. Constraint (9f) guarantees the placement of users in the coverage area of a UAV. Note that the MOOP helps to learn the trade-off between conflicting system objective functions using the concept of Pareto optimality. In the Pareto optimality or efficiency, a point is only Pareto optimal if there is no other point which improves at least one of the objective functions while keeping the same performance of the other objective functions.

Solution of the MOOP
The optimization problem (9) is non-convex due to non-convex objective functions and constraints. It is hard to find the optimal solution. We approach the proposed problem by first converting the MOOP to SOOP using the weighted Tchebycheff method. Then, we divide the SOOP in to three sub-problems and utilize alternating optimization algorithm to handle them iteratively. In the weighted Tchebycheff method, multiple objective functions are converted to a single objective by applying weights to each objective functions. The problem (9) can be reformulated as the following SOOP.
||κ n || 2 ≤ P D /N, ∀n = 1, 2, . . . , N, where α is weighting coefficient. The weighted Tchebycheff method guarantees to generate a set of Pareto optimal points, even for a non-convex MOOP. The problem (10) is still non-convex due to non-convex objective function (10a) and constraints (10b) and (10c). We divide the SOOP into three sub-problems that are described in the three subsequent subsections.

Optimal Downlink Beamformer
With fixed values for H, Φ, Q, and p n , problem (10) can be reformulated as the following sub-problem 1: Since, κ n is the only optimization variable in problem (11), the constraints that are not associated with κ n are not required. The problem (11) is still non-convex for the following two main reasons: (i) a multi-user interference term in the DL rate and (ii) an increasing function κ n of DL rate exists, while there is a decreasing function for the UL rate, which makes objective (11a) a difference of two concave functions. Therefore, we proceed with the ZF approach to simplify problem (11). We assume that for user n, DL beamformer vector κ n creates a null in the direction of the SI channel h D at the FD BS. Additionally, for user n, the projections of the channel vector h n on the DL beamformer vectors of remaining N − 1 users are 0. In other words, the following holds: h D κ n = 0, ∀n = 1, 2, . . . , N, and h n κ m = 0, ∀[m, n] = 1, 2, . . . , N, m = n.
In order to solve problem (11) efficiently, we reformulate the problem by adopting a semi-definite programming (SDP) method. To assist the SDP, we define K n = κ n κ H n , H D = h D h H D , and H n = h n h H n . The simplified SDP problem is given by: Trace(K n ) = P D /N, ∀n = 1, 2, . . . , N, The problem (12) is still non-convex because of the rank constraint (12g). Rank constrained problems are usually non-polynomial hard problems. Therefore, we relaxed the constraint (12g) by removing it from the SDP problem (12). The relaxed problem is given by: The problem (13) is known as a semi-definite relaxed (SDR) problem. Additionally, it is a convex problem and can be efficiently solved using convex problem solvers such as CVX [31]. Once solved, κ * n can be found using the eigenvalue decomposition of K n . The solution of problem (13) is optimal if the obtained matrix K n is Rank 1. Otherwise, it provides a lower bound for the problem (12). For the case when Rank(K n ) ≥ 1, ∀n = 1, 2, . . . , N, there are multiple methods in literature to improve the solution [32,33]. We follow the heuristic solution provided in [33], which implies that for user n, choose κ * n as the principle eigenvector of K n , which corresponds to the largest eigenvalue of K n .

Optimal Coverage (Altitude and Beamwidth)
With fixed κ * n , Q, Φ and p n , problem (10) can be reformulated as the following subproblem 2a: where A n = W log 2 and are the updated UL and DL rates after applying the ZF constraints of the DL beamformer, respectively. Similarly, with fixed κ * n , Q, H and p n , problem (10) can be reformulated as the following sub-problem 2b.
Note that the optimal values of H and Φ should satisfy all the respective KKT conditions. We now use a sub-gradient method to find the optimal values of all the Lagrange multipliers. The sub-gradient equations for λ, Ω, δ, U and J are given by: Here, β λ , β Ω , β δ , β U and β J are the positive gradient-search step sizes. λ n (t), Ω n (t) and δ n (t) are the values of the Lagrange multipliers at iteration t, and U n (o) and J n (o) are the values of the Lagrange multipliers at iteration o. Further, the optimal values of λ n , Ω n , δ n , U n and J n are used in their respective KKT conditions to find the optimal values of H and Φ. In the sub-gradient method, the optimal value of H should also satisfy constraint (14e) and Φ should also satisfy constraint (17d) in problem (14) and (17), respectively. It is to be noted that the optimal values of H and Φ are independently found using the sub-gradient method in iterative loops with iteration t (fixed Φ) and o (fixed H), respectively.

Optimal UAV Location and Uplink Power
For a fixed DL beamformer vector, altitude, and beamwidth, the problem (10) can be reformulated as a following sub-problem 3. The problem (32) is a non-convex optimization problem. Observing the updated DL rate of user n (i.e., B n ), the non-convexity in the problem (32) is mainly due to the UL power of N − 1 users. We use an augmented Lagrangian method (ALM) to solve the non-convex problem [35]. In ALM, constraints are augmented to the objective function to create a Lagrangian function. More specifically, a quadratic penalty term is added to the Lagrangian function, which minimizes the duality gap. Additionally, the augmented Lagrangian method is locally convex, when the penalty parameter is sufficiently large. ALM comprises of two main steps in iterations: (i) maximizing the augmented Lagrangian function, and (ii) updating the Lagrange multipliers and penalty parameters, until convergence. Note that the multipliers and penalty parameters are fixed in each iteration and updated between iterations. After applying ALM on the problem (32), the augmented Lagrangian function can be written as:

Proposed Joint Optimization Algorithm
In this work, we propose a joint optimization algorithm to solve the DL beamformer, beamwidth, altitude, and location of the UAV and also the UL power from the users in a UAV based multi-user FD communication network. Algorithm 1 illustrates the method and clarifies the joint optimization. It utilizes an alternating optimization approach which is shown in Figure 2. The algorithm comprises of three main parts as follows: Initialization, where the initial values of optimization variables, Lagrange multipliers and other variables are set; Iterative loop 1 (steps 2-21), where sub-problem 1 and sub-problem 2 (both a and b) are solved using the CVX program and the sub-gradient method, respectively; Iterative loop 2 (steps 13-18), where sub-problem 3 is solved using the ALM. Note that the loop 2 is a nested loop of the loop 1. Loop 2 iteratively runs until the convergence of ALM and loop 1 iteratively runs until all the optimization variables are converged or maximum number of iterationsv is achieved.   Calculate κ * n (v) using the eigen value decomposition of K * n (v)

7:
Else: 8: Choose κ * n (v) as the principle eigen vector of K * n (v), which corresponds to largest eigen value of K * n (v)

Complexity and Convergence Analysis
The computational complexity in Algorithm 1 is dominated by (i) step 3, where the DL beamformer of the UAV for each user n is calculated by using SDR problem (13) in a CVX program, (ii) steps 11-12, where the altitude and beamwidth of the UAV are optimized by using the sub-gradient method, and (iii) steps 14-19, where the horizontal location of the UAV and UL power of the users are independently optimized by using ALM. The complexity order for handling sub-problem 1 is O(max(N, M T ) 4 M 1/2 T log(1/ρ)) [32], for sub-problem 2a and 2b is O(N/ρ 2 ), and for sub-problem 3 is O(L 2 ), where ρ > 0 is the solution accuracy. Therefore, the overall complexity of the proposed joint optimization ). The convergence rate of the CVX program depends on the precision value, which can be manually-defined in the toolbox. For the sub-gradient method, the convergence rate is 1/ρ 2 . On the other hand, the convergence of ALM mainly relies on the adjustable penalty parameter (ξ). Generally, a faster convergence rate is obtained by a large penalty parameter. Additionally, the convergence rate tends to a constant value, which is proportional to the ratio (1/ξ), where ξ is greater than a threshold value ξ > 0.

Numerical Results
In this section, we present the simulation based results to validate our proposed optimization algorithm. The simulation parameters are given as follows. We consider that the MIMO FD UAV is equipped with M T = 8 antennas and deployed to serve N = 5 FD ground users which are located randomly in the circle with radius 100 meters. The minimum and maximum altitudes of the UAV are H min = 50 meters and H max = 200 meters. The QoS parameters are selected as R UL min n = R DL min n = 100 Kbps. The FD system bandwidth W is 10 MHz. The beamwidth range is between Φ min = 0.2 radians and Φ max = 1.2 radians. The UL reference channel gain g o = 1.42 × 10 −2 . The noise power values at the UAV and users are σ 2 D = σ 2 n = −169 dBm/Hz. For performance comparison, we consider following two baseline FD algorithms by fixing optimization variables in Algorithm 1: (i) fixed UAV location, where the location of the UAV is fixed and randomly chosen inside the radius r c , and (ii) fixed coverage (similar to the approach in [13]), where the altitude and beamwidth of the UAV are fixed and randomly chosen to satisfy constraint (9f).
We analyze the performance of the proposed system and validate its dependency on the maximum DL power in Figures 3 and 4. The maximum achievable sum FD rate for α = 0.5, P B = −10 dBm and different values of maximum DL power P D is shown in Figure 3. Similarly, the maximum achievable sum UL and DL rates are shown separately in Figure 4. It can be seen in Figure 3 that the sum FD rate is an increasing function with respect to the maximum DL power for all the considered algorithms. In Figure 4, the sum DL rate increases with P D while the sum UL rate is constant and does not depend on the P D due to the feasibility of constraint (13e) in problem (13). Moreover, even after the lower value of P B than P D and equal weights for UL and DL objectives, the sum UL rate is higher than the sum DL rate. This is mainly because of the co-channel interference which reduces the DL SINR at the users. The performance comparison of the algorithms shows that the proposed algorithm outperforms the other baseline algorithms because the proposed algorithm considers the joint optimization of DL beamfomer, beamwidth, altitude, and location of the UAV, as well as the UL power of the users. Additionally, we observed that the UAV's coverage optimization is more crucial than the location and it significantly improves the FD performance. It can be seen in Figure 3 that the fix UAV location algorithm has improved FD rate values than the fix coverage algorithm due to optimized altitude and beamwidth. 26    We analyze the trade-off between sum UL and sum DL rates of the proposed system in Figure 5. The trade-off region for P D =33.97 dBm is obtained by varying the weight coefficient α from 0.1 to 0.9. As expected, the sum DL rate decreases, while sum UL rate increases with α. From numerical results, it can be seen that the proposed algorithm provides improvement in the sum UL and DL rates compared to the baseline algorithms. To balance the rate performance on the UL and DL, α can be optimized through one dimensional search. In the case, where UL or DL is the only priority then α can be chosen as 1 or 0 in the proposed algorithm, respectively. Maximum Achievable Sum Downlink Rate (bps) 10 8 Proposed Fix-coverage Fix-UAV-location Figure 5. Trade-off region between sum UL and sum DL rates.
In addition, we also calculated the FD energy efficiency (EE) as where P D + ∑ N n=1 p n is the total FD power consumption. We examine the FD link power consumption in the proposed system by altering the maximum DL power. Figure 6 shows the maximum achievable FD energy efficiency for different values of P D . For all the considered algorithms, EE degrades with the increment in P D , mainly due to increased power consumption in the DL. It can be seen that the proposed algorithm has marginal EE improvement compared to the other baseline algorithms. To show the impact of the coverage area on sum UL and DL rates, we optimize sum UL and DL rates using fixed coverage algorithm for different values of coverage area radius r c and show the results in Figure 7. It is easy to see that both the sum rates decreases with increase in the coverage area and the effect is more prominent in the DL. Finally, we compare the proposed FD system with the similar HD system. In other words, we consider that the UAV and users are equipped with HD radio in Figure 1, where M T antennas are utilized for transmission and reception at separate and equal time instants. Therefore, we consider that R UL−HD n = 0.5R UL n , R DL−HD n = 0.5R DL n , and co-channel and self interference are neglected. Figure 8 shows the performance comparison between proposed FD and HD systems. For fair comparison, we consider that maximum UL and DL powers in the HD system are equal and defined by a single variable P in the Figure 8. Additionally, P D in the FD system is equal to P. From comparison, we found that on average, FD system provides rate improvement of 0.9399 Gbps and 1.0266 Gbps over HD system's UL and DL rates, respectively. Coverage Area Radius (m) Maximum Acheivable Rate (bps) 10 9 Proposed Proposed-HD-UL Proposed-HD-DL Figure 8. Maximum achievable sum rate versus maximum power.

Conclusions
In this paper, we have proposed a UAV based FD multiuser communication network in which a UAV is deployed as a MIMO FD BS to serve multiple FD users on the ground. We formulated a resource allocation MOOP which maximizes the achievable bi-directional sum rate of all users while providing QoS. The proposed MOOP is a non-convex problem which was converted to SOOP by using the Tchebycheff method to ease the design of a resource allocation algorithm. The SOOP was divided into three sub-problems with an independent solution of each sub-problem. Further, a joint optimization algorithm was proposed, which optimizes the resources, such as DL beamformer, UL power of users, and location and coverage of the UAV BS.
We evaluated the performance of the proposed algorithm through simulations and compared them with the performance of baseline FD algorithms. A trade-off region between sum UL and sum DL rate was unveiled, which can be used to prioritize UL or DL performance in an FD MOOP. It was found that the proposed algorithm significantly improves the FD rate performance compared to the other baseline algorithms while providing QoS to all the users. We also considered HD radio scheme in the proposed system model for comparison. Our results show that proposed FD scheme outperformed the HD scheme by 1.8616× in the UL and 2.0221× in the DL.