Sum Rate Optimization for Multiple Access in Multi-FD-UAV-Assisted NOMA-Enabled Backscatter Communication Network

: With the rapid development of the Internet of Things (IoT) network, research on low-power and energy-saving devices has attracted extensive attention from both academia and the industry. Although the backscatter devices (BDs) that utilize the environmental power to activate circuits and transmit signals are a promising technology to be deployed as IoT nodes, it is challenging to design a ﬂexible data backhaul scheme for massive BDs. Therefore, in this paper, we consider an unmanned-aerial-vehicle (UAV)-assisted backscatter communication network, where BDs are served by multiple full-duplex (FD) UAVs with the non-orthogonal multiple access (NOMA) schemes and modulate their signals on the downlink signals, which are generated by the UAVs to serve the coexisting regular user equipments (UEs). To maximize the sum rate of the considered system, we construct an optimization problem to optimize the reﬂection coefﬁcient of BDs, the downlink and the backhaul transmission power, and the trajectory of UAVs jointly. Since the formulated problem is a non-convex optimization problem and is difﬁcult to solve directly, we decouple the original problem into three sub-problems and solve them with the successive convex approximation (SCA) method, thereby addressing the original problem by a block coordinate descent (BCD)-based iterative algorithm. The simulation results show that, compared with the benchmark schemes, the proposed algorithm can obtain the highest system sum rate and utilize limited time-frequency resources more efﬁciently.


Introduction
As a typical scenario of the future wireless communication networks, the internet of things (IoT) is changing the way of future life with the emergence of new applications, such as power sensing, smart factories, and smart cities [1][2][3].With the rapid development of the IoT industry, a massive scale of devices should be densely arranged in vast areas to fulfill the sensing requirements.Under the small size of the external structure, the energy storage and power supply of sensing devices are limited.To improve the device's lifetime and save the cost of extensive device deployment, low-power communication transmission technologies have a broad vision of application and received widespread attention.As typical self-sustainable and low-power devices, backscatter devices (BDs) can transmit information by reflecting the existing environmental signals, which can be regarded as a promising solution to enable battery-free communications [4,5].

Overview and Motivation
In Ref. [6], the authors provide an overview of the fundamentals, applications and challenges of backscatter communications (BackCom).Ref. [7] provided a literature survey of the BackCom system and focused on the key signal processing techniques applied in this system.In existing works about BackCom, the authors in Refs.[8][9][10] developed theoretical analysis of the BackCom system to provide an accurate assessment of the system performance.Ref. [8] performed ergodic rate analysis and derived the upper bound of the backscatter link.Ref. [9] provided an outage performance analysis of BackComs.Ref. [10] conducted a coverage probability analysis by considering the interference of two communication links in the BackCom system.Then, the throughput and the energy efficiency of the system were both considered to be improved by optimizing the schedule strategy of BDs.In Refs.[11][12][13], the scheduling method of BDs was optimized to maximize the overall throughput of the BackCom system.In Ref. [14], a low-complexity method was proposed to jointly optimize the transmit power and the reflection coefficients of BDs.The authors in Ref. [15] considered the fairness between the BDs and optimized the resource allocation scheme to maximize the max-min throughput of the hybrid backscattering and harvest-then-transmit system.In particular, a deep reinforcement learning (DRL)-aided resource allocation method for both the BDs and UEs was investigated in Ref. [16] to improve the overall throughput of the system.In Ref. [17], a Gini threshold-based method was proposed to maximize the energy efficiency through balancing the traffic load of each BD.In Ref. [18], the authors considered the fairness of BDs and jointly optimized the transmission power and reflection coefficients to maximize the max-min energy efficiency of BDs.In addition, some of the literature investigated the code and decoding scheme of the BackCom system.Ref. [19] proposed an energy-efficient code scheme for the BackCom to minimize energy consumption.In Ref. [20], multiple readers were assumed to work cooperatively to collect the signals from BDs, and a joint receiving algorithm was designed to reduce the interference.
In the BackCom system, to assist the BDs in completing the signal transmission, the carrier emitters (CEs) and readers must be deployed to generate and receive the signals, respectively.Since BDs need a forward haul process under a passive method to transmit signals, the energy of the received signals at the reader side is usually low [7].Therefore, in a wide region where BDs are distributed, there is also a vital requirement for dense deployment of CEs and readers to ensure successful data reception and decoding, which result in an unfordable burden to low-power IoT systems.To address this issue, integrating unmanned-aerial-vehicles (UAVs) with BackCom is an attractive technology that can reduce the cost of implementing the system [21].With high flexibility and motility, the UAV can move close to the distributed BDs to efficiently complete energy supply or signal collection, making them well-suited for deployment as CEs or readers in the BackCom system.In Ref. [22], the fixed CEs on the ground and a flying reader in UAV are considered to serve the BDs.In this paper, the trajectory of the UAV, the transmitting power of CE, and the orthogonal multiple access (OMA) scheduling of BDs are jointly optimized to improve the energy efficiency of the system.In Ref. [23], UAV acts as a relay to collect the data from BDs and upload them to the base station.The authors consider the energy constraint in the UAV and optimize its locations for data collection to maximize energy efficiency.In Ref. [24], UAV is treated as a flying BD to forward the signal to the base station.The trajectory of UAV, as well as the time-splitting method between the backscattering and common transmitting, are jointly optimized to maximize the throughput.In Refs.[25][26][27][28], the UAV was operated as both the CE and the reader.In Ref. [25], the authors consider the security problem and optimize the scheduling of BDs to maximize the secrecy rate.In Ref. [26], the UAV trajectory, resource allocation, and reflection coefficient of BDs are jointly optimized to maximize the ergodic capacity of the UAV-aided BackCom system.In Ref. [27], the mobility of UAV was explored to implement the over-the-air computation for the BackCom network, where the UAV flies over the coverage area obtain the sum channel gain from BDs in several samples and utilize the samples to assist the decoding of the signals transmitted by BDs.Ref. [28] proposed a deep reinforcement learning-based algorithm to design the trajectory of UAV to improve energy efficiency.
To further improve the efficiency of data acquisition and accelerate the aggregation of data from a large number of BDs in a short time, some of the literature investigated the integration of non-orthogonal multiple access (NOMA) with the BackCom.In Ref. [29], the grouping method and the adjusting strategy of each BD was solved to maximize the sum rate in a NOMA-aided BackCom system.In Ref. [30], the fairness between BDs was considered, and the scheduling method of BDs with the NOMA scheme was optimized to maximize the max-min throughput of the system.The authors in Ref. [31] considered the power constraint of both the BD and transmitter, then proposed a power allocation scheme to maximize the capacity of this system.In Ref. [32], the authors optimized the transmitting power and the reflection coefficient for NOMA-aided BackCom by taking into account the imperfect successive interference cancellation decoder to improve the sum rate.The work in Ref. [33] investigates the BackCom system under NOMA and fullduplex (FD) transmission, where the resource allocation scheme for the fixed base station is optimized to maximize the throughput of the common users and the BDs.In particular, the studies in Refs.[34][35][36] focus on merging UAV and NOMA jointly to assist backscatter communications.According to the operation and internal structure of the devices, we can classify the backscatter devices into three categories, that is, the monostatic backscatter, the bistatic backscatter, and the ambient backscatter [6,7].In contrast to the other two BDs, the ambient backscatters modulate the signal on the available ambient sources to further reduce the cost of the BackCom system [37], which are mainly considered in this paper.
The analysis and optimization of the existing work have investigated the scenario where a single UAV is deployed to assist the BackCom.However, considering the wide distributed region and the limited energy supply capacity of BDs, a joint operation of multiple UAVs is essential for ensuring the performance of the BackCom system, which has not yet been researched as far as we know.

Paper Contribution
In this paper, we consider multiple full-duplex UAVs cooperate to serve the BDs with the NOMA scheme, while transmitting the downlink signal to normal users.We aim to provide an analytical framework for the multi-FD-UAV-assisted NOMA-enhanced backscatter communication system, and further consider the interference between the backscattering communication from the BDs and the active communication from the UEs.The main contributions of this paper are summarized as follows: • We consider a multi-FD-UAV-assisted NOMA-enabled backscatter communication network, where multiple FD-UAVs are deployed to cooperate to serve the ground BDs with the NOMA scheme while communicating with the downlink UE.Then, the UAVs forward back the data of BDs to a nearby BS through the OMA scheme.To maximize the sum rate of BDs under the considered system, we formulate a sum rate maximization problem to jointly optimize the reflection coefficient of BDs, the downlink and the backhaul transmission power, and the trajectory of UAVs with constraints such as the successive interference cancellation (SIC) decoding constraints and the backhaul capacity constraint.

•
Due to the non-convexity of the constructed problem, we split the original problem into three sub-problems, that is, the BD reflection coefficient optimization problem, the UAV transmission power optimization problem, and the UAV trajectory optimization problem.Then we apply the successive convex approximation (SCA) method to solve each sub-problem.Based on the solution of each sub-problem, we propose a block coordinate descent (BCD)-based iterative algorithm to solve the original problem.Furthermore, the convergence and complexity of the proposed algorithm are also proven.

•
To verify the effectiveness of the proposed algorithm, we provide extensive numerical results.The numerical results show that the proposed algorithm outperforms the two benchmark schemes, and the fast convergence of the proposed algorithm is also verified.

Paper Organization
The rest of the paper is organized as follows.Section 2 presents the system model of the multi-FD-UAV-assisted NOMA-enabled backscatter communication system and formulates the optimization problem.Section 3 provides the details about the proposed optimization algorithm.In Section 4, the numerical simulation results are presented.Then, Section 5 concludes this paper.

System Model
As shown in Figure 1 To illustrate the considered network clearly, we adopt the 3D Cartesian coordination to describe the positions of the UAVs, UEs, BDs, and BS.Let S u (t) = (x u (t), y u (t), h uav ), y d , 0) and W = (x, y, h BS ) denote the coordinate of the u-th UAV, the i-th UE, the d-th BD and the BS at the time t for ∀u ∈ U , ∀i ∈ I and ∀d ∈ D, respectively.For analytical convenience, we denote T as the total service time of the UAVs, and we discretize it into Q slots, each with a duration as ∆t = T/Q.Therefore, the trajectory of the UAVs can be further expressed as S u (q) = (x u (q), y u (q), h uav ) for ∀q ∈ {1, 2, • • • , Q}.In the considered system, each UAV provides service for one UE and a set of BDs, namely a cell.

UAV-BD Link Transmission Model
In this work, the BDs access the UAV by reflecting the transmitted downlink signal with the power-domain NOMA scheme.For the associations between the BDs and UAVs, each BD can only associate with one UAV in one slot.Let α d,m [q] denote the association indicator of BDs, which illustrates whether the d-th BD accesses the U-th UAV at the q-th slot.Specifically, the association indicator α d,m [q] = 1 when the d-th BD accesses the u-th UAV; otherwise α d,m [q] = 0.
Without loss of generality, we assume that information is shared between the UAVs (e.g., the UAV location information and UE data).Therefore, the UAV is able to eliminate the downlink UE signals sent by another UAV.Furthermore, the associated UAV can decode the data that the BD modulates on the transmitted signal of other UAVs.As a result, the received signal at the d-th BD is obtained as where P u [q] and s u [q] represent the downlink transmission power and the u-th UAV data symbols at the q-th slot, respectively.z d denotes the additive white Gaussian noise (AWGN) with zero mean and variance σ 2 d at the BD.In this work, all the channel links are assumed to be Rayleigh channels dependent on the path loss [41].Considering the buildings or the trees would block the line-of-sight (LoS) path, our channel hypothesis combines the LoS and non-line-of-sight (NLoS) paths, that is, we assume that all channels in the same cell are the LoS channel similar to Refs.[22,39], while the channels between different cells will be dominated by NLoS path loss.Therefore, the channel from the u-th UAV to the d-th BD at the q-th slot h u-b u,d [q] in (1) can be expressed as where f c1 and c denote the carrier frequency for the UAV downlink and speed light in vacuum, respectively.The coefficients η represent the additional attenuation.The distance from the u-th UAV to the d-th Based on the received signals, the d-th BD activates the circuit and sends its own data s d .Then the signal at the u-th UAV is given by where β d [q] denotes the reflection coefficient of the d-th BD and z u is the zero mean AWGN noise with variance σ 2 u at the UAV.Since the BDs access the UAVs by the NOMA scheme, the UAVs adopt the SIC method to decode the data of the BDs.In order to obtain a better decoding result, the decoding order of the SIC is from UEs with better channels to those with poorer channel conditions.To clarify the decoding order, we define a channel strength comparison indicator as λ i,j,u [q] to demonstrate the strength relationship between two UAV-BD channels.Specifically, λ i,j,u [q] = 0 represents that the channel strength of the BD i is weaker than that of BD j when accessing the UAV u at slot q.Then the UAV first decodes the BD j data by considering the signal of BD i as interference.Otherwise, λ i,j,u [q] = 1.
Therefore, the signal-to-interference and noise ratio (SINR) of the BD d at the UAV u is given by where Îu [q] is Îu Generally, to ensure the successful decoding of the SIC method, the SINR of accessing BDs must be greater than a threshold, which can be denoted as where κ is defined as the given SINR threshold.
Then the sum rate at the UAV u can be obtained as Specifically, according to Ref. [42], we can simplify

UAV-UE Link Transmission Model
For the ground UEs, since each UE is allowed to connect one UAV, we denote θ u,i [q] as the association indicator to describe the service state between the i-th UE and the u-th UAV at the q-th slot.The association indicator θ u,i [q] = 1 when the i-th UE is served by u-th UAV; otherwise, θ u,i [q] = 0. Therefore, the signal received at the i-th UE is given by where P u [q] and s u [q] represent the downlink transmission power and the data symbols of the u-th UAV at the q-th slot, respectively, and h u-u u,i [q] is denoted as the average channel gain from the u-th UAV to the i-th UE at the q-th slot, that is, where The second term on the right side of ( 9) is the interference of BD to the UE, and z ue,i denotes the additive AWGN at the UE with zero mean and variance σ 2 ue .Therefore, the sum rate for the i-th UE is obtained as where h b-u d,i is the average channel gain from the b-th BD to the i-th UE, which can be expressed as where is the distance from the d-th BD to the i-th UE.

UAV-BS Link Transmission Model
After decoding the signal from BDs, the UAVs will send them to the adjacent BS through the backhaul link.Due to the small number of UAVs, the backhaul link of UAV-BS is transmitted with the OMA scheme to enhance the backhaul capacity.Let P B u [q] denote the transmission power of the UAV u at the slot q in the backhaul link.The sum rate at the BS for the UAV u can be obtained as where σ 2 B is a variance of the AWGN noise at the BS and h u-B u is the average channel gain from the u-th UAV to the BS, that is, where f c2 denotes the carrier frequency for the UAV-BS backhaul link and d u-B u [q] is the distance from the U-th UAV to the BS, that is,

Problem Formulation
In this work, we aim to maximize the sum rate of the BDs by jointly optimizing the transmission power of UAVs, the reflection coefficient of BDs, and the trajectory of UAVs with the NOMA scheme.Let denote the reflection coefficient of BDs, the downlink transmission power of UAVs and the backhaul transmission power of UAVs, respectively.Therefore, the sum rate optimization with the considered multi-FD-UAV-assisted NOMAenabled backscatter communication network over Q slots is formulated as (P1) max S,P,P B ,β s.t. ( 6), ( 16) ||S ||S where the constraint ( 17) is to ensure that the capacity of the backhaul link is greater than the data transmitted by the BDs to prevent data from accumulating and overflowing at the UAVs.To satisfy the quality of service (QoS) of the UE, we have the minimum sum rate constraint for UEs as shown in (18), where R UE min denote the minimum sum rate for each UE over all the slots.The reflection coefficient limitation of BDs is shown in (19).The inequalities (20) and ( 21) represent the speed constraint and the trajectory constraint, where V max is the maximum flying velocity.Given the minimum distance between two UAVs d min , the inequality (22) ensures that collisions between UAVs do not occur.The inequality ( 23)-( 25) are the power limitation of the UAVs, where P max represents the maximum transmission power of the UAV.

Proposed Algorithm for Sum Rate Maximization
Since the formulated Problem (P1) is intractable to solve due to the non-convexity of the objective function and the constraints, in this section, we propose a BCD-based iterative algorithm to jointly optimize the reflection coefficient, transmission power, and trajectory of the UAV.Specifically, we split Problem (P1) into three independent sub-problems, that is, the reflection coefficient optimization problem, the transmission power of the UAV problem, and the trajectory optimization problem.Then, we use the SCA method to convert each sub-problem into a convex problem and solve it.Furthermore, the complexity and the convergence of the proposed algorithm are given in this section.

Reflection Coefficients Optimization
Given the transmission power and the trajectory of the UAV, the reflection coefficients optimization problem can be obtained as s.t. ( 6), ( 17), ( 18), (19).
However, the formulated Problem (P2) is still non-convex due to the fractional forms in the logarithm function.Therefore, we will simplify the logarithmic form and adopt the SCA method to solve the problem.
For the non-convexity of R u [q], based on (8), we first split it into R u [q] = R u [q] − Ru [q], where R u [q] and Ru [q] are denoted as which are both concave with respect to (w.r.t.) β d [q].Therefore, to transform Ru [q] to a convex function, we apply the first-order Taylor expansion (FTE) to Ru [q] and we can obtain its upper bound as where β r k [q] is the local point at the r-th SCA iteration and A r u [q] is denoted as Therefore, the non-convexity of the objective function R u [q] is addressed.For the backhaul capacity constraint (17), since R B u is convex w.r.t.β d [q], we need to transform R u [q] to a convex function.Therefore, based on (30), we can obtain the lower bound of R u [q] at the r-th iteration as where B r u [q] is denoted as Therefore, the constraint ( 17) is converted to be convex as For the constraint of QoS of the UE, following the similar steps of the transformation of the objective function R u [q], we can divide R u i [q] into two parts as , which can be specifically expressed as Similarly, according to FTE, Ru i [q] is upper-bounded by where C r i [q] is given by Therefore, the constraint ( 18) is converted to be convex as As for the SINR constraint (6), we can make it convex by the simple algebraic operation, which is expressed as Consequently, with the SCA method, we can recast Problem (P2) at the r-th iteration as 34), ( 34), ( 39), (40).(42) Since Problem (P2.1) is convex w.r.t. the reflection coefficient β, it can be solved directly by the standard convex tools.

UAV Transmission Power Optimization
With the given reflection coefficients of the BD β and the trajectory of the UAV S, the optimization problem w.r.t. the transmission power of UAV can be formulated as (P3.1) max s.t. ( 6), ( 17), ( 18), ( 23), ( 24), (25).(44) For the above Problem (3.1), we can find the constraints ( 23), (24), and ( 25) are jointly convex w.r.t.P and P B .As for the constraint (6), we can rewrite it as , where Rup u [q] is obtained by applying FTE to Ru [q] with the given point P r u [r] at the r-th SCA iteration, which can be expressed as where E r u [q] is given by For the backhaul capacity constraint (17) and UE QoS constraint (18), by using the same method as optimizing β, we can obtain the transformed convex constraint as where F r u [q] and G r i [q] are given by Therefore, the transmission power optimization problem can be reformulated as s.t. ( 23), ( 24), ( 25), ( 45), ( 48), (49).
Since Problem (3.2) is a convex problem w.r.t.P and P B , the standard convex tool can be applied directly to solve it.
Note that due to the non-convexity of the objective function and constraints, Problem (3.1) is neither concave nor convex; thus, we will adopt the SAC technique to solve it.
According to the channel definition (2) and the split of the sum rate of UAV R u [q], we can re-express R u [q] and Ru [q] w.r.t. the trajectory of UAV S as where According to Ref. [43], R u [q] and Ru [q] are both convex w.r.t.||S j [q] − B d || for ∀j ∈ U , d ∈ D. Therefore, we introduce the slack variable with the following constraint as As a result, by applying FTE to ||S u [q] − B d || with the given point S r u [q] at the r-th SCA iteration, the constraint (59) can be transformed into which is convex w.r.t. S.
Similarly, for the UE QoS constraint (18), by introducing the slack variable π = {π u,i [q], ∀u ∈ U , i ∈ I, q ∈ Q} and the FTE, we can obtain the convex form of (18) as where 2 and can be converted to be convex by FTE as R u,lb i [q] is obtained by using FTE to R u i [q], which can be expressed as where K r i [q] is given by Finally, to address the non-convexity of the SINR constraint of the BD (6), we convert (6) to a convex function with the slack variable ω as where As a result, with the given point S r , we can reformulate Problem (3.1) as 20), ( 21), (60), (66), (67), (68), (71), (73) where ( 74) is obtained by applying FTE to ||S i [q] − S j [q]|| 2 .Since the reformulated Problem (3.2) is convex w.r.t. the UAV trajectory S, it can be solved efficiently by the standard convex tools.

Proposed Algorithm for Problem (P1)
Based on the solution of each sub-problem, we propose a BCD-based iterative optimization algorithm, which is shown in Algorithm 1. Specifically, the original problem is decomposed into three sub-problems corresponding to the problems in Sections 3.1-3.3.When optimizing each variable, other variables keep unchanged and the obtained optimal solution are the input of the next optimization problem.The whole iteration does not end until convergence.
It is clear from Algorithm 1 that the overall complexity of the proposed algorithm mainly comes from solving the three convex sub-optimization problems.The interior-point method was adopted to solve the standard convex optimization problem.Therefore, we can obtain the complexity of solving convex optimization Problem (P2.2),Problem (P3.2) and Problem (P4.2) as O((DQ) 3.5 ), O((2UQ) 3.5 ) and O((2QU + IUQ + UDQ) 3.5 ), respectively [44].Assuming that the number of iterations of the algorithm is I all , the overall complexity of the proposed algorithm is O I all (DQ) 3.5 + (2UQ) 3.5 + (2QU + IUQ + UDQ) 3.5 .

Simulation Results
In this section, we provide extensive numerical results to verify the efficiency of the proposed algorithm.We assume all the UEs and BDs are randomly distributed in the UAV service region.The location of the base station is (0,800,0).Unless otherwise stated, the remaining system parameters are listed in Table 1 [21][22][23].Furthermore, the initial transmission power of the UAV is evenly allocated for the downlink and the backhaul link.The BD reflection coefficient is set to 1.
In order to intuitively demonstrate the effect of the proposed algorithm, we display the optimized trajectory of the UAV in Figure 2. In Figure 2, the black squares and red squares represent the positions of BD and UE, respectively, and the trajectory of the UAV is plotted by the solid lines with different markers.The direction of the arrow represents the movement direction of the UAVs.It can be noted that in order to maximize the sum rate of the BDs, the UAV will traverse all the locations of the BD.Meanwhile, to keep the QoS of the UE, it will also approach the UE during the flight.
Figure 3 plots the optimized downlink transmission power of the UAV w.r.t.different slots.It can be noted that the downlink transmission power of UAV-1 is higher than that of UAV-2.This is because UAV-1 is closer to the BS, so it needs less backhaul transmission power to forward the data of BDs.Furthermore, as shown in Figure 3, the transmission power variation of the two UAVs is the opposite of satisfying the SINR constraint of the BD and the QoS constraint of the UE.When a UAV is close to the served BD, and another UAV is far away from the served BD, the farther away UAV will increase the transmission power to obtain a higher rate, while the closer UAV will decrease the transmission power to reduce interference to another BD and UE, thereby enhancing the sum rate with the constraint of the BD and UE satisfied.To illustrate the superiority of the proposed algorithm, we consider two benchmark schemes: (1) OMA scheme, that is, the UAV can access only one BD at each slot with the proposed algorithm; (2) fixed location, that is, the UAV stays hovering.Figure 4 demonstrates the sum rate w.r.t. the number of slots.As shown in Figure 4, more service time increases the system's sum rate.Furthermore, the performance of the proposed algorithm is better than that of the two fixed trajectory baseline algorithms, which shows the effectiveness of the proposed algorithm.Meanwhile, although the OMA scheme is also optimized with the same algorithm, its performance is still worse than that of the NOMA scheme, which shows the superiority of the NOMA scheme in the sum rate.Figure 5 shows the sum rate w.r.t. the number of BDs.In the simulation, we keep the total service slots Q and the initial trajectory of the UAV unchanged.With the increase in the number of BDs, the sum rate performance of the four schemes will be improved.However, too much BD makes the improvement of the sum rate inefficient.This is because the increase in space density of BD reduces the freedom of UAV trajectory optimization.Meanwhile, the superiority of the proposed optimization algorithm and the NOMA scheme are also intuitively shown in Figure 5.
We demonstrate the impact of the height of the UAV on the system sum rate in Figure 6.It can be noted that the sum rate of all three schemes decreases as the height of the UAV increases.This is due to the fact that the altitude of the UAV affects the quality of the channel; that is, as the UAV increases its flight altitude, the distance from the UAV to the BD, the user, and the BS also increases, resulting in greater path loss.Therefore, the increase in the UAV height reduces the sum rate of the system.In addition, we can see that the proposed algorithm still outperforms the other two schemes, which verifies the effectiveness of the proposed algorithm again.
Finally, we demonstrate the convergence of the proposed algorithm in Figure 7.It can be observed that the proposed algorithm can reach the optimal state within three iterations under different slots, which verifies the faster convergence of Algorithm 1.

Conclusions
In this work, we investigated a multi-FD-UAV-assisted, NOMA-enabled backscatter communication network, where multiple FD UAVs were deployed to serve the downlink UEs and also collect the data reflected by the BDs with the NOMA scheme, and a BS provided the backhaul link of the UAV.Then, we formulated an optimization problem to maximize the sum rate of the BDs by jointly optimizing the reflection coefficient of BDs, the downlink and the backhaul transmission power, and the trajectory of UAVs.To solve this non-convex problem, we decomposed it into three independent sub-problems and adopted the SCA method to solve them.Based on the transformed sub-problems, we proposed a BCD-based iterative algorithm to tackle the original problem.The simulation results validate the superiority and the fast convergence of the proposed algorithm and the efficiency of the developed NOMA system with limited time-frequency resources.

Figure 4 .
Figure 4.The sum rate w.r.t. the number of slots.

Figure 5 .Figure 6 .
Figure 5.The sum rate w.r.t. the number of BDs.

Figure 7 .
Figure 7.The sum rate w.r.t. the number of iterations.
which is convex w.r.t.P.As for the rest of the non-convex constraints, since Problem (3.1) w.r.t. the transmission power of UAV P and P B has the same structure as Problem (P2) w.r.t. the reflection coefficients β, we can use the same transformation steps to convert Problem (3.1) into a convex problem.Therefore, we can first obtain the convex lower bound of R u