The present application relates to an unmanned aerial vehicle (UAV)-aided over-the-air computing system and a trajectory and power optimization method thereof, and in particular to a UAV-aided over-the-air computing system based on full-duplex relay and a trajectory and power optimization method.
Due to the advantages of strong mobility, flexible configuration and line-of-sight link, Unmanned Aerial Vehicle (UAV) is widely used in the field of wireless communication. UAV can also move to a place close enough to the sensor in the harsh field, which avoids long-distance information transmission, saves the power of the sensor and mitigates the influence of noise. In the air-to-ground transmission, the UAV is high in altitude and has usually a line-of-sight wireless transmission with the sensor, so the probability of channel depth fading is reduced. Because sensors and base stations can't communicate directly, using UAV as relay has become an important research direction of information collection in the Internet of Things based on UAV.
In 2021, in “UAV-assisted over-the-air computation” published by Min Fu et al., it is proposed to use high mobility and wireless line-of-sight transmission capability of UAV to assist the over-the-air computing system, so as to minimize the mean square error of over-the-air computing. In this system, UAV receives the information from sensors through the fusion of airborne base stations, and transmits and fuses the data of multiple sensors in a single time slot. However, in the wireless sensor network under complex environment, there is no direct communication between sensors and base stations.
The objective of the present application is to provide an unmanned aerial vehicle (UAV)-aided over-the-air computing system based on full-duplex relay which realizes direct communication between sensors and a base station and a trajectory and power optimization method thereof.
The UAV-aided over-the-air computing system includes one base station, one UAV and multiple sensors placed on the ground.
The base station receives information from the UAV; the multiple sensors transmit information to the UAV at the same time.
As a full-duplex relay, the UAV works in a Fusion and Forward (FF) mode, receiving fused information transmitted by the multiple sensors and transmitting the information to the base station at the same time.
The UAV collects and fuses data of the sensors in a way of over-the-air computing, and forwards the data to the base station in a way of full-duplex relay; the UAV flies according to an optimized flight trajectory.
A trajectory and power optimization method of the application includes following steps:
Optionally, in the S1, the optimization problem is established in the coordinate system, and an expression of the optimization problem is:
where
Optionally, in the S1, an expression of the denoising factor optimization sub-problem is:
When the denoising factor η[n] is optimized, it is necessary to give the transmitting power pk[n] of the sensors k, the UAV transmitting power P[n] and the UAV flight trajectory q[n].
Optionally, in the S1, an expression of the sensor transmitting power optimization sub-problem is:
When the transmitting power pk[n] of sensors k is optimized, it is necessary to give the denoising factor η[n], the UAV transmitting power P[n] and the UAV flight position q[n].
Optionally, an expression of the UAV transmitting power optimization sub-problem is:
When the UAV transmitting power P[n] is optimized, it is necessary to give the denoising factor η[n], the transmitting power pk[n] of the sensors k and the UAV flight position q[n].
Optionally, in the S1, the UAV flight trajectory optimization sub-problem is optimized by adopting a convex optimization method, and an expression of the UAV flight trajectory optimization sub-problem:
When the UAV flight position q[n] is optimized, it is necessary to give the denoising factor η[n], the transmitting power pk[n] of the sensors k and the UAV transmitting power P[n].
Optionally, the S2 is realized as follows:
Compared with the prior art, the application has the following remarkable effects: firstly, the application adopts the UAV as the full-duplex relay to receive and fuse all sensor data, estimates an interesting function, and simultaneously transmits the estimated function value of the time slots to the base station, thus realizing the minimum mean square error under the guarantee of communication rate; secondly, in the process of jointly optimizing the UAV trajectory and sensor power, the optimization problem is decomposed into four independent optimization sub-problems, and the UAV trajectory and power are optimized with low complexity algorithm.
The present application will be further described in detail with reference to the drawings and specific embodiments of the specification.
As shown in
A trajectory and power optimization method is realized as follows:
Horizontal coordinates of the sensors k are expressed as wk=[xk, yk]∈□1×2, where xk and yk represent an abscissa and an ordinate of the sensors, respectively.
The sensor group on the ground is represented as K□{1, 2, . . . K}, K>1, where K is the total number of the sensors.
During the flight of UAV, positions of the sensors on the ground are fixed, and the UAV has stored position information of the sensors. In addition, the UAV flies at a fixed altitude from the ground, denoted as H, a minimum flight altitude to ensure obstacle avoidance without frequent ascending and descending of the UAV. L is a fixed distance between a sending end and a receiving end of the UAV.
A time-varying trajectory of horizontal projection of the UAV is q(t)=[x(t), y(t)]∈□1×2, a starting position q[0]=[x0, y0] and an ending position q[T]=[xT, yT], where x0 and y0 represent an abscissa and an ordinate of the sensors at an initial position, respectively and xT and yT represent an abscissa and an ordinate of the sensors at ending position, respectively.
A time discretization method is adopted to deal with the continuous UAV trajectory. A duration T of a task is equally divided into N time slots: T=Nδ, where δ is an time step. An appropriate time step is selected so that a distance between the UAV and the sensors is approximately constant in each time slot, that is δVmax□H, where the Vmax is a maximum flight speed of the UAV. In time slots n, a movement constraint of the UAV in flight is expressed as:
∥q[n]−q[n−1]∥2≤Vmaxδ,n=1,2, . . . N (1),
q[0]=[x0,y0] (2),
q[N]=[x
N
,y
N] (3),
where the q[n] is a UAV flight position in the time slots n, and the q[n−1] is a UAV flight position in time slots n−1.
The target of UAV computing is fused data of all sensors on the ground, so a target function ƒ[n] of the UAV computing is expressed as:
where the ϕ represents a post-processing function of the UAV, the ψk presents a pre-processing function at the sensors k, the Zk[n] is data in the time slots n, and q[n] is the UAV flight position.
Pre-processed transmission signals of the sensor are sk[n]□ψk(Zk[n]), and assuming that the transmission signals are independent of each other, they are normalized by zero mean and unit variance, namely: E(sk[n])=0, E(sk[n]skH[n])=1, E(si[n]sj[n]H)=0, ∀i≠j. Therefore, after post-processing of averaging, a processing function received by the UAV is:
As the full-duplex relay, the UAV receives data from the sensors in each time slot and sends the data to the base station. A received signal y[n] of the UAV in the time slots n is:
where the sk[n] and the su[n] are transmitted signal of the sensors and transmitted signal at the sending end of the UAV, respectively, the βu is a self-interference cancellation coefficient, the bsk[n] and the bu[n] are transmitted precoding coefficients of sensors k and the UAV, the e[n] represents additive white Gaussian noise, and the hsk[n] and the hu[n] are a channel model of the sensors and the UAV and a channel model of the sending end to the receiving end of the UAV respectively.
A constraint of the transmitting power of the sensors k is:
E(|bsk[n]sk[n]|2)=|bsk[n]|2≤Pk (7),
An estimated average value {circumflex over (ƒ)}[n] of UAV transmission data is:
where η[n] is the denoising factor and K is the total number of the sensors.
Test performance is carried out with mean square error MSE[n], then:
where the Wk is the fixed horizontal position of the sensors, β0 is the channel gain per unit distance, βu is the self-interference cancellation coefficient, σ2 is the additive white Gaussian noise power, L is the distance from the sending end to the receiving end of the UAV, H is the lowest flight altitude at which the UAV does not need to ascend and descend in flight, and α is the path loss index, α≥2.
Therefore, the following optimization problem problem1 is established:
where the Gmin is the minimum information transmission rate between the UAV and the BS (BS: base station), and the B is a communication bandwidth; the pk[n] is the transmitting power of the sensors k in time slots n, the P[n] is the UAV transmitting power in time slots n, the UAV flight position in time slots n, and the η[n] is denoising factor in time slots n.
From the optimization problem problem1, it can be seen that optimization variables are highly coupled, so an iterative alternate optimization method is adopted to solve.
Sub-problem 1: when the denoising factor η[n] is optimized, it is necessary to give the transmitting power pk[n] of the sensors k, the UAV transmitting power P[n] and the UAV flight position q[n]. At this time, the sub-problem 1 is expressed as:
The optimization problem is decomposed into N sub-problems, and each sub-problem η[n] is optimized to minimize the mean square error of one time slot. Then the n-th sub-problem is expressed as:
Letting γ[n]=1/√{square root over (η[n])}, the problem represented by formula (13) is transformed into a convex quadratic problem, expressed as:
By setting a first derivative of an objective function of formula (14) to zero, the optimal solution is obtained:
Sub-problem 2: when the transmitting power pk[n] of the sensors k is optimized, it is necessary to give the denoising factor η[n], the UAV transmitting power P[n] and the UAV flight position q[n]. At this time, the sub-problem 2 is expressed as:
Because both
in the objective function are constants, the
and the
are ignored. The sub-problem 2 is decomposed into the following K sub-problems:
Since formula (17) is a typical convex linear constrained quadratic programming problem, the formula (17) is solved by a standard convex optimization method.
Sub-problem 3: when the UAV transmitting power P[n] is optimized, it is necessary to give the denoising factor η[n], the transmitting power pk[n] of the sensors k and the UAV flight position q[n]. At this time, the sub-problem 3 is expressed as:
For formula (18), since a constant term is ignored, so the formula (18) is solvable.
Sub-problem 4: when the UAV flight position q[n] is optimized, it is necessary to give the denoising factor η[n], the transmitting power pk[n] of the sensors k and the UAV transmitting power P[n]. At this time, the sub-problem 4 is expressed as:
By introducing a relaxation variable s={sk[n]=∥q[n]−wk∥22, ∀k, ∀n}, the sub-problem 4 is expressed as:
According to Taylor's formula, a global lower bound ĝklb[n] is obtained:
An inequality is obtained at the same time:
∥q[n]−wk∥22≥∥qr[n]−wk∥22+2(qr[n]−wk)T(q[n]−qr[n]) (25).
The sub-problem 4 is further transformed into:
Therefore, formula (26) transformed from the sub-problem 4 is a Quadratical Constraint Quadratic Programming (QCQP) problem and is solved by the standard convex optimization method. Through continuous iterative solution, the optimal power and UAV trajectory are finally obtained.
In order to minimize the time average mean square error of the system, the application adopts iterative optimization algorithm to solve each sub-problem step by step, and implementation steps are as follows:
Number | Date | Country | Kind |
---|---|---|---|
CN202111561590.1 | Dec 2021 | CN | national |
This application is a continuation of PCT/CN2022/105164, filed on Jul. 12, 2022, which claims priority to Chinese Patent Application No. 202111561590.1, filed on Dec. 16, 2021, the contents of which are hereby incorporated by reference.
Number | Date | Country | |
---|---|---|---|
Parent | PCT/CN2022/105164 | Jul 2022 | US |
Child | 18124751 | US |