Abstract
In this paper, we consider a ground terminal (GT) to an unmanned aerial vehicle (UAV) wireless communication system where data from GTs are collected by an unmanned aerial vehicle. We propose to use the ground terminalUAV (GU) region for the energy consumption model. In particular, to fulfill the data collection task with a minimum energy both of the GTs and UAV, an algorithm that combines optimal trajectory design and resource allocation scheme is proposed which is supposed to solve the optimization problem approximately. We initialize the UAV’s trajectory firstly. Then, the optimal UAV trajectory and GT’s resource allocation are obtained by using the successive convex optimization and Lagrange duality. Moreover, we come up with an efficient algorithm aimed to find an approximate solution by jointly optimizing trajectory and resource allocation. Numerical results show that the proposed solution is efficient. Compared with the benchmark scheme which did not adopt optimizing trajectory, the solution we propose engenders significant performance in energy efficiency.
1. Introduction
With series of features such as low cost, long duration, high flexibility, and high adaptability, extensive research endeavour has been rendered to exploring the application of UAV. In [1], the authors made a survey on autonomous cargo pickup with an UAV helicopter by designing systems including selftracking, payload pickup, and deployment of cargo. In [2], the UAVs are employed to airborne maritime surveillance, which commits to monitoring of the marine environment, safeguarding national security, sovereignty, and as sea rescue. Anwar et al. [3] studied the framework of detecting and classifying sounds of amateur aircraft out in noisy environment by leveraging machine learning to maintain public safety. They also [4] proposed a disasterresilient architecture including centralized control and edge computing for disaster areas, which can reduce delay effectively. Moreover, UAV networks are developing rapidly embracing the wireless technologies such as IEEE 802.11n, IEEE 802.15.4, and 3G/LTE [5–9]. However, UAVs will also face many challenges in its application. Researchers are also trying to solve these problems. In particular, it can not only be hired for delivery of goods, but also serve as lowaltitude aircraft to improve the coverage and rate of wireless networks in different scenarios. First, placement/deployment optimization is a major challenge for quasistatic drones [10–14]. In [10], coverage was optimized by adjusting the height of the platform. In [11], the authors characterize the performance and tradeoffs of UAV with underlaid devicetodevice communication. Compared with terrestrial base stations, UAV has more advantages as airborne static base stations, for its excellent capability to have lineofsight communication links (LOSs) to the ground terminal, which will achieve the maximum communication coverage of the ground terminal. Thus, in [12], the authors propose an efficient deployment method based on circle packing theory, using minimum transmit power and maximizing the total coverage area. BorYaliniz et al. [13] emphasized the performances of dronecell placement problem by designing it as a 3D placement that leads to maximum revenue of the network. Second, UAV’s energy efficiency and resource management are crucial factors in its communication systems. There are many differences between traditional wireless communication system and UAV communication. Energy efficiency maximization of traditional wireless communication systems is to save energy of communication, but in the UAV communication system, the energy for communication is much lower than that of the UAV itself. So, the UAV’s energy consumption is the main factor that needs to be considered toward high efficiency and energy saving [14–16]. Inspired by this, trajectory optimization of energyefficient UAV communication has been studied, in which sequential convex optimization techniques are applied to deal with the nonconvex trajectory optimization problems [17], despite its research object being single UAV and a ground user, but the frame can extend to a joint trajectory and communication design in multiUAVs wireless network. Motivated by this, the joint optimization problem of user scheduling, transmit power, and UAV trajectories was studied in the multiUAV system [18]. In IoT networks, an efficient mobile scheme of uplink data collection was studied [19]. Wu and Zhang [20] consider the user communication delay requirements, maximizing the system total throughput via joint resource allocation and UAV trajectory. Furthermore, energy saving has been recognized as a vital target in designing future wireless communication system [21]. In particular, UAV communication systems are considered. In [22], the author proposed a throughputdelay tradeoff in a new cyclical multiple access scheme in UAV communications. In the previous work, we have investigated an energy tradeoff between a UAV and single ground user [23], with results being extendable to UAVenabled multiuser communication networks.
As we discussed previously, energy efficiency is a factor that cannot be ignored for UAV due to the limited onboard power. He et al. [24] investigated the throughput by jointly optimizing the UAV’s flying altitude and antenna beamwidth in case of downlink multicasting, downlink broadcasting, and uplink multiple access models without considering energy consumption. Lyu et al. [25] optimize the deployment of the UAV to provide wireless coverage for ground terminal which does not study the issue of energy efficiency. Zeng et al. [26] aim to design trajectories to minimize the mission completion time and also hardly consider energy consumption. Zeng and Zhang [17] study energy efficiency by trajectory optimization, but only a UAV and a user included. For a UAV communication system, besides coverage, UAV’s energy efficiency and resource management is necessary to investigate, especially in multiuser UAV communication.
In this paper, a flexible UAV is deployed to collect data from a group of ground terminals at known location in ground terminals to a UAV (G2U) system. Intuitively, the study needs to jointly consider the uplink transmission energy of GTs and the UAV’s propulsion energy consumption. To obtain the most fundamental insights, we focus on a G2U wireless communication system, where a group of ground terminals are collected data by a UAV, as shown in Figure 1. Specially, an efficient algorithm that combines optimal trajectory design and resource allocation scheme is proposed to approximately solve the optimization problem. For optimal trajectory design and resource allocation, we gain the optimal GTs transmit power and UAV trajectory, respectively. Numerical results evaluate the tradeoff involved in the proposed designs. A special case where a G2U communication scheme considering two ordinary UAV trajectories with circular and straight trajectories has been studied in [23].
The rest of this paper is organized as follows. Section 2 introduces the system model, and problem formulation for a highmobile UAV is deployed regularly to total ground nodes for data collection. In Section 3, we first formulate the optimization problem to describe the Pareto boundary of the energy region. An efficient iterative algorithm is proposed by combining optimal trajectory design and resource allocation scheme; then the numerical results are solved in Section 4. Finally, we conclude the paper in Section 5.
2. System Model and Problem Formulation
2.1. System Model
We consider a wireless system consisting of K ground nodes, denoted as , where the coordinate of is denoted as . We assume that each of the ground node periodically generates data of size for each period of duration , which need to be regularly uploaded to the fusion center for data processing. This could correspond to the practical periodic sensing applications. We assume that there is no established communication link between the ground nodes and fusion center due to their long separation distance as well as the high cost in deploying fixed communication infrastructure such as relays. Instead, we assume that a highmobile UAV is despatched regularly to fly over the ground nodes for data collection and then brings it back to the fusion center for offline processing.
We assume that the UAV flying at a fixed altitude H around the ground terminals for data collection for a finite time horizon T. Note that T is the time required for the UAV to complete the data collection task for each period . In practice, we usually have (e.g., dozens of minutes versus several days) since the data are only generated intermittently. We denote the timevarying coordinate of the UAV’s horizontal location as , . Thus, represents the UAV trajectory, which is to be optimized. We assume that the initial and final locations of the UAV are predetermined, which are denoted as and , respectively. Furthermore, we assume that the maximum UAV speed and acceleration is and , respectively. Thus, the UAV trajectory must satisfy and , where and denotes the instantaneous velocity and acceleration, respectively.
We assume that the channels between the UAV and the ground nodes are dominated by LoS links. Furthermore, we assume that the Doppler effect due to the UAV mobility can be perfectly compensated. Thus, the channel coefficient from the ground terminal to the UAV at time t can then be expressed aswhere denotes the channel gain at the reference distance meter, and represents the distance between ground node and the UAV. Denote by B the total bandwidth assigned for the uplink data collection system, which is assumed to be dynamically allocated to the K ground terminals for orthogonal multiple access. Specifically, denote by the fractional of the spectrum bandwidth that is assigned to at time t. Note that, in practical implementation, is determined by the number of subbands allocated for normalized by the total number of frequency subbands, which has discrete values in general. However, we assume that the number of subbands is sufficiently large so that can be approximated to have continuous values between 0 and 1 for simplicity. Denote by the transmit power of at time t. Then, the maximum transmission rate for the uplink from to the UAV can be expressed aswhere is the noise power spectrum density in Watts/Hz, and represents the reference SNR. It is observed from (2) that the data collection rate for is a function of , as well as the UAV trajectory . The aggregated transmission throughput over the time horizon T for can be presented aswhere the total rate of each ground terminal should be larger than the data collection task : .
There are two parts of energy consumption in the UAVenabled data collection system:(1)Ground Terminal Energy Consumption. The total energy consumption for ground terminals can be written as(2)UAV Energy Consumption. The UAV energy consumption of the considered data collection system consists of two main components. The first one is the energy consumption related to the communication functions, such as circuitry and signal processing. The other component is the propulsion energy consumption, which is required for the UAV to remain aloft as well as for supporting its mobility. In practice, the communicationrelated energy is much smaller than the propulsion energy and hence is neglected in this paper.
For fixedwing UAVs, the total propulsion energy as a function of the trajectory can be modelled as [21]where and are the parameters depending on the aircraft’s weight, wing area, air density, etc, m is the mass of the aircraft including all it payload, and denotes the gravitational acceleration with nominal value , and and are the instantaneous UAV velocity and acceleration vectors, respectively, which are related to its trajectory as
The energy consumption for UAV can be expressed as
Note that the achievable data rate for each ground terminal is a function of total operating time T, spectrum, and power allocation , as well as UAV trajectory . Denote by the set of strategies in terms of the resource allocations and UAV trajectories such that the data collection requirement can be fulfilled. The set can be expressed as
Consider the data collection system which both ground terminals and UAV would consume energy to support the data transmission and UAV flying. There are two parts of energy consumption: (1) ground terminal energy and (2) UAV energy . It is seen that the optimal strategies for minimizing the transmitting power of ground terminals and UAV energy consumption for flying are in general different. For data transmission, UAV would be slowly flying when it is near to the ground terminals. On the other hand, for UAV’ flying, the steady level flight would be a better choice for energy decreasing. It thus motivates our investigation of the following question: what are the optimal UAV flying trajectory strategies for data transmission and UAV flight energy consumption? To answer this question, we propose to use ground terminalUAV (GU) region for the energy consumption model to characterize all the ground terminal transmit power and UAV energy consumption pairs under a given data collection rate constraint. We define the feasible energy region to be the set of all energy pairs of ground terminal and UAV that can be achieved using the wireless resource and trajectory vector that satisfy the rate constraint. It can be presented as
The outer boundary of this region is called the Pareto boundary, because it consists of operating points for which it is impossible to decrease one of the energies, without simultaneously increasing one of the other energies. More precisely, we define the Pareto optimality of an operating point as follows:
Definition 1. An energy pair is Pareto optimal if there is no other energy pair with and (the inequality is componentwise) can meet the rate constraints for the data collection system.
The Pareto boundary of the energy region characterizes that the minimum energy consumption for UAV and ground terminal for data collection task. It can be adopted to evaluate the tradeoff of the energy consumption in the data collection system. It is an interesting topic to investigate the UAV trajectory strategy and resource allocation scheme for UAV and ground terminals to finish the data collection task cooperatively.
3. Multiple Ground Terminals with UAV Arbitrary Flight
In this section, we address the data collection system in a general scenario where there are multiple ground terminals which is served by a flexible flying UAV. We first formulate the optimization problem to describe the Pareto boundary of the energy region. Then, we investigate the optimal solution of the optimization problem for this scenario, from which we obtain an upper bound for the achievable energy consumption pairs in the GU region. Then, we propose an alternating iterative method to derive the optimal resource allocation scheme and optimal trajectory design strategy.
3.1. Problem Formulation for Optimal Resource Allocation and Trajectory Design
In order to characterize the tradeoff of the energy region effectively, we adopt the strategy to minimize the ground terminal energy consumption with fixed UAV energy consumption. By traversing , the Pareto boundary of GU region can be obtained. In this case, the design objective for data collection system is to minimize the energy consumption for ground terminals with the fixed UAV energy consumption. The optimization problem can be formulated as
The constraint is the data collection requirement for each ground node, which denotes the total amount of data to be collected for . is the energy consumption constraint for UAV, with denoting the total available energy consumption for UAV. represent the UAV’s initial/final location and velocity constraints, respectively. It is noted that the amount of available energy for UAV is larger than the minimum energy consumption for UAV flying from the initial location to final location. Without considering the UAV’s energy consumption, the minimum transmit power for ground terminal is the situation that the UAV would be sequentially located above each ground terminal during the data collection time T. However, the optimal solution to the above problem is difficult to be directly solved.
3.2. Linear Discrete StateSpace Approximation
Note that the problem (P2) is difficult to be directly solved since the trajectory is a continuous variable respect to the optimization problem and its first and secondorder derivatives. To tackle the problem more tractable, we utilize the linear discrete statespace approximation to reformulate the optimization problem (P2).
For ease of exposition, the time horizon T is discreteized into equally spaced time slots with step , where is chosen to be sufficiently small so that the UAV location can be assumed to be approximately constant. Therefore, the UAV trajectory can be represented by the sequence . Moreover, and are the instantaneous UAV velocity and acceleration vectors, respectively. Considering the relationship between and , and are first and secondorder derivatives, respectively. Then, applying Taylor approximation, the velocity and acceleration are related to its trajectory as
Correspondingly, the spectrum and power allocation can be expressed as for the discretization.
For the problem (P2), the objective is to minimize the total ground node transmit power while satisfying the data collection requirements of all the ground terminals, subject to a fixed UAV energy consumption. It should jointly optimize the operation period N with fixed step size , the UAV trajectory , and spectrum and power allocation . The problem (P2) can be reformulated as
The constraint , , are the discrete mode for the data collection requirement for each ground node, energy consumption constraint for UAV, and spectrum allocation constraint for ground terminals, respectively. to represent the UAV’s initial/final location and velocity constraints, respectively. It is noted that the amount of available energy for UAV is larger than the minimum energy consumption for UAV flying from the initial location to final location.
Problem (13) is nonconvex; hence, it cannot be directly solved by standard convex optimization techniques. In the following sections, we propose an efficient algorithm to find an approximate solution to (13) via iteratively optimizing trajectory and resource optimizations.
3.3. Combining Optimal Trajectory Design and Resource Allocation Scheme for Multiuser UAV Network
In this section, an efficient algorithm that combines optimal trajectory design and resource allocation scheme is proposed to approximately solve the original optimization problem. The sequential convex optimization method is applied to meet the optimized goal by iteratively obtaining the optimal trajectory and allocating the wireless resource to the ground terminals.
The original optimization problem can be classified into two subproblems: (1) wireless resource allocation issue for the multiple ground nodes and (2) the UAV trajectory optimization problem. To make the optimization problem more trackable, we propose to adopt the iteratively optimization method which assumes that one issue is determined when the other issue is considered to be optimized, and vice versa.
3.3.1. Minimizing Energy Consumption for Constrained UAV Flight
Firstly, we should find a feasible solution for the UAV flight trajectory from the initial point to the final location by the predefined available energy . In other words, we should examine whether the UAV can finish the flight under the initial/finial location, velocity constraints within the operation time slots without considering any rate constraints of ground terminals. The optimization problem can be formulated as
It is noted that the constraints are convex sets. However, the objective function is nonconvex. The problem cannot be directly solved with the standard convex optimization techniques. We should adopt some relaxation method to construct the standard convex form and solve the problem (P3.1) efficiently.
The UAV energy consumption of the objective function can be upperbounded by the following equation:where denotes the amount of the UAV’s kinetic energy variation, which would be determined as a constant value by the initial and final velocity constraints . It is noted that the energy consumption upper bound is tight for the circular flight mode, in which case . Therefore, the objective function can be approximately formulated as
Considering the function is still a nonconvex set for the variable , we propose to introduce a slack variable to reformulate the objective function as
It can be shown that we must have , for all n $ at the optimal solution for (P3.1.2), since otherwise one can always increase and decrease when the UAV near the ground terminal to decrease the objective function value and obtain a strictly lower energy consumption. Until now, the objective is a convex with respect to , but with a new nonconvex constraint . Similarly, to tackle the nonconvexity of the constraint, a local convex approximation is adopted to tackle this issue. For any given local point , we have the following expression by applying the firstorder Taylor expansion:where the equality holds at the point . The function is a linear function. Then, the new constraint is a convex set, which can be expressed as
Then, the inequality in firstorder Taylor expansion shows that the new constructed convex constraint always implies the nonconvex constraint , but the reverse is not true in general.
Then, the optimization problem () can be reformulated as
Based on the previous discussions, the optimal solution of (P3.1.3) is the upper bound to that of problem (P3.1). Fortunately, the optimization problem (P3.1.3) is a standard convex optimization problem with the convex objective function with all convex constraints. It can be solved efficiently via the bisection method or the standard Dinkelbach’s algorithm.
3.3.2. Optimal Resource Allocation with Fixed Trajectory
For the wireless resource allocation issue, it is proposed to minimize the ground nodes transmit power under the fixed UAV trajectory assumption. Then, the optimization problem for wireless resource allocation can be formulated as
This subproblem may correspond to the practical scenario when the UAV’s trajectory is predetermined due to other tasks (e.g., surveillance) rather than data collection. In this case, the optimal solution of the subproblem (P3.2) can minimize the transmit power of the ground nodes and extend the life time of the node.
Since the trajectory is fixed, the constraint is a convex set combined with respect to . It can be verified that the subproblem (P3.2) satisfied the Slater’s condition. Using the dual Lagrangian method, the optimal power allocation can be expressed aswhere , is the Lagrange multiplier (or dual variable) associated with the constraint.
; is the optimal frequency band allocation factor. This is one instance of the waterfilling solution to the optimal power allocation.
Substitute the optimal in the equation, the optimal frequency band allocation factor can be obtained as equation (22):where is the Lagrange multiplier (or dual variable) associated with the constraint ; can be any real number between 0 and 1, which can be calculated by satisfying the constraints and . Generally, the optimal solution for is binary solution if the amount of time slot n is large enough. Otherwise, is adopted to finetune the resource allocation solution to meet the constraint.
It is proposed to adopt an iterative subgradient method (Algorithm 1) to optimize the dual variable : the subgradient for of ground node k is given by and for of each time slot is given by . The update procedure of Lagrangian dual variables are performed aswhere are the step size for the iterative search process. Several step size rules have been proven to guarantee convergence under some general conditions.
The optimal solution includes two steps: (1) each of subproblem in slot n can be solved with fixed Lagrangian multipliers using the combinatorial method, (2) the convergency updates of Lagrangian multipliers are performed to meet the constraints. For the optimal dual variables optimization, it can also adopt the ellipsoid method to find the optimal solution. Our proposed procedure above is convenient for the transceiver processing for both UAV and ground nodes.
3.3.3. UAV Trajectory Optimization with Fixed Resource Allocation
Since the ground node has been allocated with fixed frequency band and power resource, we should adopt a dual problem to achieve the optimization goal by optimizing the UAV trajectory. In this section, we consider the issue to maximize the total rate of ground nodes by optimizing the UAV’s trajectory with fixed frequency band and power allocation. It is noted that maximizing the total rate of ground nodes can minimize the total energy consumption of ground nodes with the rate constraints. They are dual problems with each other. The optimization subproblem (P3.3) can be formulated as
It is noted that the constraints are convex sets. However, the constraint and are nonconvex, and the objective function is nonconvex neither. The problem cannot be directly solved with the standard convex optimization techniques. We should adopt some relaxation method to construct the standard convex form and solve the problem (P3.3) efficiently.
Firstly, to tackle the nonconvexity of the objective function, for any local point , define function for each ground node aswhere
Note that is a concave function with respect to . Then the upper contour set is convex set. Moreover, as mentioned in [21], holds true for any given . The lower bound equation is the first order Taylor expansion of a convex differential function which is its global underestimator. Both and have an identical gradient.
Therefore, the optimization problem (P3.3) can be approximately solved by maximizing the lower bound aswhere the constraint of is reformulated as a convex set .
Secondly, considering the constraint set is still a nonconvex set for the variable . Just like the approximation method as in minimizing energy consumption for constrained UAV flight, we propose to introduce a slack variable to reformulate the constraint as
Then, all the constraints are convex set. Until now, we reformulate the optimization problem for any given local point , and it can be written as:
Based on the previous discussions, the optimal solution of (P3.3.2) is the lower bound to that of problem (P3.3). Fortunately, the optimization problem (P3.3.2) is a standard convex optimization problem with the convex objective function with all convex constraints. It can be solved efficiently via the bisection method or the standard Dinkelbach’s algorithm. Then, the nonconvex problem (P3.3) can be solved by iteratively optimization (P3.3.2) with local point updated in each iteration. It can be summarized as the following table (Algorithm 2).

In summary, through adopting the relaxation method and sequential convex optimization technique, an efficient solution is proposed to solve the UAV trajectory optimization problem which is guaranteed that the optimal point is fulfill the Karush–Kuhn–Tucker conditions of the original nonconvex problem (P3.3).
3.3.4. Iterative Method for Optimal Resource Allocation and UAV Trajectory Design
Until now, we have solved the issues of UAV trajectory optimal design and resource allocation for ground terminals when one of them is fixed. In order to obtain the optimal solution for the original optimization problem (P3), we propose an alternate iterative solution to optimize the UAV trajectory and resource allocation, based on the methods mentioned in Sections 3.3.1.–3.3.3.
The proposed optimal solution can be summarized as follows (Algorithm 3).

Note that the solution for problem (P3.2) adopts an iterative subgradient method in Algorithm 1 whose complexity is where l is the iteration time. As Algorithm 2 requires to solve convex optimization problems, it has a polynomial complexity in the worst case each iteration [27]. Thus, the complexity of Algorithm 3 combining Algorithms 1 and 2 is the polynomial complexity because polynomial complexity multiplied by iterations is still polynomial complexity which is affordable.
4. Simulation Results and Analysis
In this part, we provide numerical results to prove the reliability of the proposed algorithm. We assume that ground users are randomly distributed in this system with an area of . Other reference values are as follows: UAV’s altitude is . The communication bandwidth is and noise power is assumed to be . Besides, we set that and , the maximum speed of UAVs are assumed as , and the maximum acceleration of UAVs are assumed as and = 100 Mbits (Mb).
Firstly, we evaluate the energy tradeoff between the consumption of ground terminals and UAV.
We consider a certain trajectory. Multiple ground terminal and UAV energy region are shown in Figure 2. We can see clearly from the figure that the energy consumption of the ground user decreases as the drone energy consumption increases with different data collect requirement of . More concretely, when increasing the UAV’s energy from 22 kJoule to 31 kJoule under the requirement of = 80 Mbits (Mb), the GTs energy consumption would be significantly reduced from 32 Joule to 27 Joule.
(a)
(b)
By comparing with the resource allocation optimization shown in Figure 3, the solution we propose that optimizes the resource allocation and trajectory jointly delivers significant performance in energy efficiency. Specifically, the GTs’ energy consumption of combined resource allocation optimization is 27 Joule when the UAV’s energy consumption is 35 kJoule, while that of resource allocation optimization only is 86 Joule. It can be observed that the total energy consumption for ground terminals achieved by the proposed algorithm dramatically decreases with the number of iterations and the algorithm converges in about 4–6 iterations from Figure 4. Thus, the solution we propose is able to converge fast and the computational complexity is suitable for practical implementation.
Next, we study the UAV trajectory optimization.
As shown in Figure 5, the UAV flies to the top of each GT eventually to fulfill the communication task. It is observed that the obtained trajectories tend to be converged for comparing the trajectory of fifth iteration and that of seventh iteration proposed by the algorithm. We also show the scheduling of every GT in Figure 6 where the UAV maintains communicating with GT all the time in order to save the energy.
(a)
(b)
We also collect the UAV trajectory optimization results. In multiuser scenario, we obtained the UAV trajectory under the constraints of limited flight energy and data collection rate, minimizing the energy consumption problem of user communication. The UAV trajectories under different UAV energy consumption are shown in Figure 7. With the increase of , the UAV hovers around more on the top of GTs so as to save that of the GTs but consume energy of the UAV. It is beneficial for energylimited GTs.
(a)
(b)
(c)
(d)
5. Conclusions
In this paper, we have studied the optimal resource allocation scheme and optimal trajectory design strategy for multiple ground terminals with UAV arbitrary flight. First, UAV’s energy consumption and GT’s energy consumption are derived, and we describe their Pareto optimal tradeoffs and adopt the strategy to minimize the ground terminal energy consumption with fixed UAV energy consumption and consider rate constraints of ground terminals. Second, through calculating minimum energy consumption for constrained UAV flight, we gain the initial UAV flight in the GU system. Next, we propose to get the optimal resource allocation under the fixed flight. Then, through maximizing the total rate of ground nodes by optimizing the UAVs trajectory with fixed frequency band and power allocation. Finally, we use an alternate iterative solution to optimize the UAV trajectory and resource allocation. Simulation results show the alternate solution we proposed can significantly improve the performance in energy efficiency compared with the scheme without optimizing trajectory. Through this paper, we consider a pointtopoint G2U communication scheme formerly [23], extending to a UAV to multiusers on designing energysaving UAV communication. We hope that there will be more development in the future wireless communication.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Authors’ Contributions
Dingcheng Yang and Qingmin Zhao are contributed equally to this work.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (grant nos. 61703197, 61561032, and 61461029), Graduate Student Innovation Special Funds of Nanchang University (grant no. CX2019077), and Key Research and Development Program of Jiangxi Province (grant no. 20182ABC28008).