The present invention belongs to the technical field of data acquisition and optimization for UAV uplink communication, relates to a design solution for high efficiency green communication of a single unmanned aerial vehicle (UAV) and ground sensors, and particularly relates to a method for jointly optimizing the flight trajectory, sensor wake-up scheduling and time slot by UAV during data acquisition to achieve the purpose of maximizing the energy efficiency of the system.
UAV is an emerging technology that has been widely used in military, public and civil fields due to its high maneuverability and low cost. With the development of Internet and Internet of Things technology in the future, the UAV can assist in satisfying communication needs of mass connection and high information rate. Application scenarios of UAV assisted information transmission and data acquisition have been extensively studied. Traditionally, the information of the wireless sensor network is transmitted to the data center through multi-hop. Each sensor node not only sends own data, but also forwards data from other nodes, resulting in that the energy of some sensor nodes may be quickly exhausted. Network connection is intermittent and is difficult to guarantee the fairness of users. The UAV is used as an auxiliary mobile node, and each sensor can directly send information to the UAV, thereby ensuring the fairness of the users. In addition, the general Line of Sight (LOS) channel conditions between the UAV and sensors make the information transmission rates faster, and high maneuverability of the UAV is also more suitable for wireless communication in complex environments.
Although the UAV communication system is widely used, the limited on-board energy of the UAV fundamentally limits the endurance and communication time of the UAV. Therefore, it is important to maximize energy efficiency in UAV green communication. The energy efficiency is defined as the transmission bit of information per unit of energy consumption, which can directly increase the amount of information that the UAV can communicate before recall, giving consideration to communication service quality and UAV energy consumption. The present invention appropriately designs the parameters in the system for the purpose of maximizing energy efficiency.
The purpose of the present invention is to solve the problem of high energy efficiency green communication in UAV data acquisition systems. In a network of uplink communication of a single UAV and I ground sensors, the UAV receives data periodically. The specific solution is shown in schematic diagram 1. The UAV flight trajectory W, the sensor wake-up scheduling S and the time slot t are jointly optimized to ensure that the transmission information amount and energy consumption of the sensor satisfy system requirements, while maximizing the energy efficiency EE of the system.
To achieve the above purpose, the present invention adopts the following technical solution:
A design method of a high energy efficiency UAV green data acquisition system comprises the following steps:
Step 1, constructing a system optimization objective:
(1) Serving a set of I ground sensors which are randomly distributed through time division multiple access (TDMA) by a UAV. The UAV and the sensors are respectively provided with a single antenna, and the sensors do not interfere with each other during serving.
(2) Flying, by the UAV, at a fixed altitude H with a maximum flight speed of Vm and a total cycle of T, and discretizing the cycle T into N time slots by a time discretization method, with the length of each time slot of
The coordinate of the UAV is w[n]=[x(n),y(n)]T∈R2×1 in time slot n, wherein x(n), y(n) are the x-coordinate and y-coordinate of the UAV respectively, and R2×1 is a two-dimensional vector space. For a sensor set SI={1, 2, . . . , I} of random distribution, the coordinate of sensor i is fixed as Li=[xi,yi]T ∈R2×1,i ∈SI, each sensor supports the energy of Ei, i∈SI and the data amount to be transmitted is Bi,i∈SI. If communication between the UAV and ground is line of sight (LoS) link communication, channel quality only depends on a distance between the UAV and the sensor, and power gain in unit reference distance is expressed as ρ0 Then the channel power gain hi[n] of the sensor i in the time slot n conforms to a free space path loss model, i.e.,
is the distance between the UAV and the sensor i in the three-dimensional space.
(3) Assuming that the UAV serves only one sensor in one time slot, and defining a binary variable Si[n]∈{0,1} to represent wake-up scheduling of the sensor. Si[n]=1 indicates that the UAV establishes communication with the sensor i in the time slot n. Si[n]=0 indicates that the UAV does not establish communication with the sensor i in the time slot n. Then, the information transmission rate Rui[n] between the UAV and the sensor i in the time slot n can be expressed as:
wherein σ2 is additive white gaussian noise (AWGN) at a receiving end of the UAV, and PA is the transmission power of a ground sensor during communication. The unit of Rui[n] is bps/Hz. The total information amount
For a rotary-wing UAV, when parameters are constant, the propulsion power P(V) of the UAV is mainly related to flight speed V, which can be expressed as:
The propulsion power is composed of three parts: blade profile power, parasite power and induced power. Due to time discretization, the speed of the time slot n can be approximately expressed as
and Δn is defined as the flight distance of the time slot n. Then the propulsion power Pprop[n] of the time slot n can be approximated expressed by the following formula:
In the formula, P0 and Pi are the blade profile power and the induced power respectively in a hovering state; Ω is blade angular velocity; r is rotor radius; d0 represents fuselage drag ratio; ρ is air density; s is rotor solidity; A is rotor disc area; v0 is mean rotor induced velocity; the above parameters are constants. The total propulsion energy E consumed by the UAV in one cycle of serving can be expressed as:
According to the definition of energy efficiency, the system optimization objective can be represented as:
Step 2, constructing an optimization problem according to the energy efficiency formula in step 1, wherein an optimization objective is maximization of EE({W},{t},{S}), and constraints comprise UAV trajectory constraints, sensor wake-up scheduling constraints, sensor energy constraint and data amount constraint to construct the following optimization problem:
In the optimization problem, formulas (7b)-7(d) are trajectory constraints, Vm is the maximum speed of the UAV and the UAV returns to an initial position after flying by a cycle. Formulas (7e) and (7f) are the sensor wake-up scheduling constraints. Formula (7g) is the data amount constraint of the sensor and Bi is the data amount to be transmitted by the sensor i. Formula (7h) is the sensor energy constraint and Ei is maximum energy supported by the sensor i in each cycle.
The above optimization problem is a non-convex optimization problem. An original problem (7) is decomposed into two approximate concave-convex fractional sub-problems based on a block coordinate descent method and a successive convex approximation technique, and Dinkelbach algorithm is used to obtain a suboptimal solution.
Step 3, decomposing an original problem into two sub-problems according to a block coordinate descent method. For the two sub-problems, approximately converting two non-convex problems into two convex optimization problems by a successive convex approximation technique and designing algorithm solving, as follows:
(1) Optimization Sub-Problem of Wake-Up Scheduling S and Time Slot t
fixing UAV trajectory W so that the sub-problem is the non-convex optimization problem of wake-up scheduling S and time slot t; firstly, for a binary variable S, slacking S to a continuous variable within a range [0,1]; then, introducing an auxiliary variable z[n] to satisfy
using z[n] to replace the third term of the propulsion power Pprop[n] in formula (4) to obtain the UAV propulsion power PpropA[n] under the sub-problem; introducing an auxiliary variable R_t[i] to satisfy
after introducing the auxiliary variables, applying the successive convex approximation technique for non-convex constraints, converting hyperbolic constraints into SOCP and approximating the original non-convex sub-problem as a convex problem, which can be expressed as:
In the sub-problem (8), PpropA[n] is the propulsion power after the auxiliary variable z[n] is introduced, and is a convex function of t and z[n]; R_tlb[i] is the lower bound of first-order taylor expansion of the auxiliary variable R_t[i]2, and is a linear function of R_t[i]; (1/t)lb is the lower bound of first-order taylor expansion of 1/t, and has a linear relationship with t. The constraints of the sub-problem (8) are convex constraints; the optimization objective (8a) is a standard concave-convex fractional programming problem with concave numerator over convex denominator; and the existing Dinkelbach algorithm and convex optimization tool CVX can be used for calculation. It should be noted that because the constraint range is reduced by the successive convex approximation technique, the optimal solution of the convex problem after approximation is the lower bound of the optimal solution of an original sub-problem.
(2) Optimization Sub-Problem of UAV Trajectory W
fixing wake-up scheduling S and time slot t so that the sub-problem is a non-convex optimization problem of the UAV trajectory W; introducing an auxiliary variable y[n] to satisfy
using y[n] to replace the third term of the propulsion power Pprop[n] in formula (4) to obtain the UAV propulsion power PpropB[n] under the sub-problem; after introducing the auxiliary variables, applying the successive convex approximation technique for non-convex constraints, and approximating the original non-convex sub-problem as a convex problem, which can be expressed as:
In the sub-problem (9), PpropB[n] is the propulsion power after the auxiliary variable y[n] is introduced, and is a convex function of w[n]; Rui,lb[n] is the lower bound of first-order taylor expansion of the information transmission rate Rui[n] on ∥w[n]−Li∥, and is a concave function of w[n]. The solving method of the sub-problem (9) is the same as that of the sub-problem (8); and the Dinkelbach algorithm and convex optimization tool CVX can be used for solving. The optimal solution of the convex problem after approximation is the lower bound of the optimal solution of the original sub-problem.
(4) Overall Iterative Algorithm Design
Based on the above results, the present invention proposes an overall iterative algorithm. In each iteration, by solving the sub-problem (8) and the sub-problem (9), alternately optimizing wake-up scheduling S, the time slot t and the UAV trajectory W. The solution obtained in each iteration is used as the input of next iteration. The termination condition for iteration is that the increase of optimization values of one iteration and the previous iteration is less than a set threshold. A specific algorithm process is as follows:
4.1) Setting an iteration termination threshold ε, an initial trajectory w0 and an iteration index r=0.
4.2) In the r+1 iteration, using the trajectory wr obtained from the r iteration to solve the sub-problem (8) to obtain the optimization result of the sub-problem (8) of the r+1 iteration, namely, wake-up scheduling Sr+1 and time slot tr+1.
4.3) Solving the sub-problem (9) by the given wr, Sr+1 and tr+1, to obtain the optimization result of the sub-problem (9) of the r+1 iteration, namely trajectory wr+1.
4.4) If the increase of an optimization target value is greater than a threshold ε, then updating the iteration index r=r+1; skipping back to step 4.2) for the next iteration; and if the increase of the target value is less than the threshold ε, terminating the iteration.
The present invention has the beneficial effects that: the UAV flight trajectory, the sensor wake-up scheduling and the time slot are jointly optimized to illustrate how to realize the energy saving communication of UAV data acquisition with energy efficiency as an index, so as to provide a reference value method for maximization of the energy efficiency.
The present invention is described below in detail in combination with the drawings and embodiments.
It is assumed that a UAV serves 6 ground sensors which are randomly distributed. The UAV flies at a fixed altitude H=100 m with a maximum flight speed Vm=50 m/s. One cycle T is fixedly divided into N=60 time slots. The coordinates of the sensors are expressed with a matrix as L=[−1100,500;−425,400;600,1100;200,200;800,−400;−700,−600]T. Additive white gaussian noise (AWGN) at the receiving end of the UAV is σ2=−110 dBm, and power gain of reference distance is p0=−60 dB. The transmission power of the ground sensors is PA=0.1 W. If the UAV flies above the sensors, the channel power gain is
In the case, the information transmission rate Rui[n] in formula (1) is the maximum, and the maximum is Rui[n]=9.9672 bps/Hz.
For the parameters in formula (3), the parameter values of the classic rotary-wing UAV is taken in embodiment 1: blade angular velocity Ω=300 r/s; rotor radius r=0.4 m; fuselage drag ratio d0=0.6; air density ρ=1.225 kg/m3; rotor solidity s=0.05; rotor disc area A=0.503 m2; mean rotor induced velocity v0=4.03 m/s. Then, in formula (3), the velocity that minimizes the propulsion power P(V) is Vmin=10.0125 m/s.
In this scenario, the present invention assumes that each sensor needs to transmit the same amount of data and has the same energy constraints. Namely, Bi=B and Ei=E. The above parameters are substituted into the optimization problem (7) for solving, to obtain the trajectory design for maximizing energy efficiency proposed in the present invention, as shown in
According to the design scenario of embodiment 1, in order to demonstrate the superiority of the present invention, this section proposes two other benchmark solutions and compares the performance. Solution 1: energy efficiency maximization solution (the present invention). Solution 2: flying—hovering solution. Solution 3: energy efficiency maximization solution under fixed circular trajectories.
The above embodiments only express the implementation of the present invention, and shall not be interpreted as a limitation to the scope of the patent for the present invention. It should be noted that, for those skilled in the art, several variations and improvements can also be made without departing from the concept of the present invention, all of which belong to the protection scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
202110134201.0 | Jan 2021 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2021/099267 | 6/10/2021 | WO |