I Introduction
Deploying unmanned aerial vehicles (UAVs) for wireless networking is a flexible and costeffective approach to providing ondemand communications [12]. However, the limited bandwidth of radio frequency bands and limited energy of UAVs confine the applicability of UAVenabled wireless networks. This challenge can be addressed by equipping UAVs with visible light communication (VLC) capabilities. Indeed, VLC has recently attracted attention due to its large licensefree bandwidth and high energy efficiency. For instance, a VLC system that uses lightemitting diodes (LEDs) to transmit wireless signals can provide both illumination and communication services. Moreover, the altitude of the UAVs ensures the line of sight channel for VLC. Therefore, using VLC can be a promising approach to provide energyefficient UAV communications with sufficiently available bandwidth. However, deploying VLCenabled UAVs also faces many challenges that range from illumination interference detection and prediction to UAV deployment optimization and energy efficiency.
The existing literature such as in [7, 10, 3, 4, 13] has studied a number of problems related to UAV deployment. In [7], the authors proposed an optimization framework in UAVassisted backscatter networks to identify the tradeoff between UAV’s altitude, number of backscatter devices and backscatter coefficients. The authors in [10]
derived the average coverage probability and the system sumrate as a function of the UAV altitude and the number of users. However, the works in
[7] and [10] only consider the altitude of UAVs without optimizing their locations. In [3], the authors studied the optimal UAVs’ locations based on the prediction of human behavior so as to optimize the qualityofexperience of wireless devices. The authors in [4]studied the use of neural network based learning algorithm to optimize the performance of the UAV based wireless networks. The work in
[13] optimized the locations of UAVs in twouser broadcast channel and proved that the UAVs must hover around the user who has a large rate requirement. However, the works in [3, 4, 13] ignored the energy efficiency of UAVs in optimizing the deployment of UAVs. Moreover, all of the existing works such as in [7, 10, 3, 4, 13] are over limited capacity radio frequency bands which may not allow the UAVs to meet the high data rate demands of ground users. Instead, VLCenabled UAVs can be considered to provide high speed communications [14]. However, nighttime illumination such as vehicle lights, street lights and building lights causes strong interference to VLC link. Therefore, it is necessary to analyze the illumination distribution of the service areas so as to optimize the deployment of VLCenabled UAV.The main contribution of this work is a novel framework for dynamically optimizing the locations of VLCenabled UAVs based on accurate predictions of the illumination distribution of a given area. To the best of our knowledge, this is the first work that studies the use of the predictions of the illumination distribution to provide a powerefficient deployment of VLCenabled UAVs. Our key contributions include:

We consider a VLCenabled UAV network, in which the UAVs must find their optimal locations by predicting the distribution of ambient lighting so as to provide illumination as well as communication services to ground users. This problem is formulated as an optimization problem whose goal is to minimize the total transmit power of UAVs under both illumination and communication constraints.

To solve this optimization problem, we first use a Gaussian mixture model (GMM) to fit the illumination distribution, which can be obtained using an expectationmaximization (EM) algorithm. Based on derived illumination distribution, a gated recurrent units (GRUs) based prediction algorithm is proposed to predict future illumination distribution. The proposed prediction algorithm can model the temporal characteristics of the longterm historical illumination distribution thus enabling the UAVs to predict future illumination distributions.

Given the predicted illumination distribution, the original optimization problem is proved to be a convex problem which is then solved using duality.

Simulation results show that the proposed approach can achieve up to 22.1% reduction in terms of transmit power compared to a conventional optimal UAV deployment without considering illumination distribution. Furthermore, simulations indicate that UAVs should hover over areas with strong illumination.
The rest of this paper is organized as follows. The system model and the problem formulation are described in Section II. The use of the GMM to model illumination distribution, GRUbased algorithm for illumination prediction, and the optimization of the UAV deployment are proposed in Section III. In Section IV, the numerical results are presented and discussed. Finally, conclusions are drawn in Section V.
Ii System Model and Problem Formulation
Consider a wireless network composed of a set of D VLCenabled UAVs that serve a set of U ground users distributed over a geographical area . The UAVs provide downlink transmission and illumination simultaneously as shown in Fig. 1. Hereinafter, we use aerial cell to refer to the service area of each UAV. We assume that each UAV only serves the users located in its aerial cell .
Given a UAV located at and a ground user located at , the channel gain of the VLC link between UAV and user can be given by [9]:
(1) 
where is the detector area and is the distance between UAV and ground user . is the Lambert order with being the transmitter semiangle (at half power). and represent the angle of incidence and irradiance, respectively. As such, . is the receiver field of vision (FOV) semiangle. The gain of optical concentrator is defined as:
(2) 
where represents the refractive index.
For static user located at , the channel capacity at time can be given by:
(3) 
where is the illumination target, is the transmit power of UAV for user at time , and
represents the standard deviation of the additive white Gaussian noise. We define the illumination distribution of the target area as
that will be specified in Section III. Note that the illumination distribution is the luminance of the given area and it is also the interference over the VLC transmission links caused by illumination due to human activity. Here, is the interference over the VLC link between the UAV and the user located at .Due to the limited energy of UAVs, their deployment must be optimized to minimize the transmit power while satisfying the data rate and illumination requirements of users. As done in [8], we do not consider the mobility energy consumption of the UAVs.
Iia Problem Formulation
To formalize the deployment problem, we must first determine the minimum transmit power that each UAV uses to meet the data rate and illumination requirements of its associated users. To satisfy the data rate constraint of each user located at , the required power of UAV at time is:
(4) 
A UAV can successfully satisfy all the users’ requirements once the user that has the maximum power requirement is satisfied. Therefore, the minimum transmit power of UAV satisfying the data rate requirements of its associated users is given by:
(5) 
where is the set of users associated with UAV .
Given this system model, our goal is to find an effective deployment of UAVs that meets the data rate requirements of each user while minimizing the transmit power of the UAVs. This problem involves predicting the illumination and adjusting the locations as well as the transmit powers of the UAVs. The optimization problem is formulated as follows:
(6)  
(6a)  
(6b) 
where denotes the illumination demand and is the illumination of UAV at time [6]. (6a) indicates that the UAV needs to provide illumination to meet the illumination threshold of each user . (6b) indicates that the transmit power of UAV should satisfy the data rate requirements of its associated users from (5). Here, we note that ambient illumination increases the VLC interference while reducing the illuminance threshold. Since the illumination distribution is timevarying, it is necessary to predict the illumination distribution of the target area to deploy the UAVs at the beginning of each time interval. Hence, we introduce a machine learning algorithm to predict the illumination distribution of the service area. Based on the prediction, UAVs are deployed according to the solution of (6) which remain unchanged during each prediction period.
Iii Machine Learning for Illumination Prediction and UAV Deployment
Here, we first introduce the use of GMM to model the illumination distribution. Then, we propose a GRUbased machine learning algorithm for predicting the illumination distribution. The proposed learning algorithm enables the UAVs to analyze the relationship among historical illumination distributions, to predict the future illumination distribution. Based on the prediction, the UAVs can be optimally deployed.
Iiia GMM Fitting of Illumination Distribution
Since illumination is caused by human activities such as business and industrial operation, the illumination distribution of a given area is centered on several small areas. Therefore, we assume that the illumination distribution of area at time follows a GMM given by:
(7) 
where is the GMM amplitude and is the number of Gaussian components with with being the weight of each component. For each Gaussian component , is the central coordinate while and represent the standard deviation on the and axes, respectively. The parameters related to Gaussian component at time
can be represented as a vector
. Given (7), the spatial features of illumination distribution at time can be represented by a vector, which can be given by:(8)  
The optimal vector that is used to model the illumination distribution can be determined using the EM algorithm [1].
IiiB Illumination Distribution Prediction
Next, we introduce the use of GRUs [5]
for the prediction of the illumination distribution. GRUs are extensions of conventional recurrent neural networks (RNNs)
[2]. GRUs can effectively solve the gradient vanishing and the gradient exploding problem in longterm memory RNNs. Due to interconnected neurons at hidden layers and their internal gating mechanisms, GRUs can model the temporal characteristics of the longterm illumination distribution. In addition, GRUs can dynamically update the model based on the current illumination distribution due to the variablelength recurrent structure, hence, GRUs enable the UAVs to predict future illumination distribution.
A GRUbased prediction algorithm consists of three components: a) input, b) output, and c) GRU model. The key components of our GRUbased prediction approach are:

Input: The input of the GRUbased prediction algorithm is a matrix that represents the time series data of illumination distribution, which can be given by where is the length of the time series.

Output: The output of the GRUbased prediction algorithm is a vector , that represents the illumination distribution at time slot .

GRU model: A GRU model is used to approximate the function between the input and output , thus building a relationship between historical illumination distribution and future illumination distribution. A GRU model is essentially a dynamic neural network that consists of an input layer, a hidden layer, and an output layer. The hidden states of the units of the in hidden layer at time are used to store information related to the illumination distribution from time slot to . For each time , the hidden states of the GRU are updated based on the input and . Next, we introduce how to update the hidden state of hidden unit given a new illumination distribution .
At each time slot , the hidden state is determined by two gates: reset gate and update gate . First, the reset gate is used to determine the historical illumination distribution information retained in the candidate hidden state , which can be given by:
(9) where
is the logistic sigmoid function and
is element of a vector. and represent the weight matrices of reset gate, where is the length of the input and is the number of the units in hidden layer. Based on the value of the reset gate , the candidate hidden state that is used to combine the input illumination distribution with the previous memory is given by:(10) where is a reset gate vector at time and is an elementwise multiplication. For example, given two vectors and , . and represent the hidden state weight matrices.
Similarly, the update gate is used to decide the size of the information stored in the candidate hidden state to update the hidden state , which can be given by:
(11) where and represent the weight matrices of the update gate. The actual hidden state of hidden unit is updated by:
(12) The proposed GRU model iteratively updates the hidden states to store the input until the hidden state of the current time is computed. The output layer of the GRU model will predict the illumination distribution at time based on the hidden state :
(13) where is the output weight matrix. In fact, (13) is used to build the relationship between output and the hidden state that stores the information of input .
To build this relationship, a batch gradient descent approach is used to train the weight matrices which are initially generated randomly via a uniform distribution. The update rule of the gradient descent approach is:
(14) where is the learning rate, and . and
are the gradients of the loss function
which is defined as:(15) where is the illumination distribution prediction and is the actual illumination distribution at time . The specific process of using the GRUbased prediction algorithm to predict the illumination distribution for each UAV is summarized in Algorithm 1.
IiiC Optimization of the UAV Deployment
Once the illumination distribution is predicted, the UAVs can determine their optimal deployment at the beginning of each time interval by solving the optimization problem defined in (6). As analyzed in Section II, a UAV only needs to consider the user with the maximum power requirement since, by doing so, the requirements of all other users will be automatically satisfied. Therefore, substituting (1) and (5) into (6), we have:
(16)  
(16a)  
(16b) 
where , and .
Since the service area of each UAV does not overlap with the service areas of other UAVs, we ignore the interference caused by other UAVs. In consequence, problem (16) can be decoupled into multiple subproblems. For each UAV , the location optimization subproblem can be formulated as:
(17)  
(17a) 
where .
Due to its convex objective functions and constraints, problem (17) is a convex problem, which can be optimally solved by using the dual method [15]. The Lagrange function of problem (17) will be:
(18) 
where is the dual variable associated with constraint in (17a).
Given , and , the value of can be determined by the gradient method [15]. The updating procedure is:
(24) 
where is a dynamic stepsize. With regards to the optimality, we formulate the original optimization problem as a convex problem and, hence, it can always converge to the optimal solution.
The proposed algorithm used to solve problem in (6) is summarized in Algorithm 2. The complexity of the proposed algorithm lies in training a GRUbased prediction model and updating UAV location . A ground server can be deployed to train a GRUbased prediction model using Algorithm 1. Based on the trained model, the complexity of predicting the illumination distribution at the next time slot is . In addition, the complexity of calculating of each UAV is , where is the number of iterations of UAV until (17) convergence and is the number of users covered by UAV . Therefore, the proposed algorithm can run independently on each UAV due to the linear algorithm complexity.
Iv Simulation Results and Analysis
For our simulations, a m m square area is considered with uniformly distributed users and UAV. The downlink rate requirement of each user is generated randomly and uniformly over [0.5,1.5] Mbps. Other parameters are listed in Table I. The time series illumination data used to train GRU based prediction algorithm is a dataset of average radiance composite nighttime remote sensing images, obtained from the Earth observations group (EOG) at NOAA/NCEI [11].
Parameters  Value  Parameters  Value 

1 cm  1.5  
10 m  1  
16  
358  3  
0.01  100 
Fig. 2 shows how the EM algorithm approximates the illumination distribution of the target area using a GMM model. Fig. 2(a) shows the actual illumination distribution and Fig. 2(b) shows the illumination distribution of the target area modeled by GMM. In Fig. 2(b), we can see that EM algorithm uses 3 Gaussian components to model the actual illumination distribution. This is due to the fact that the area in Fig. 2(a) has only three lighting places.
Fig. 3 evaluates the prediction accuracy of our GRUbased algorithm. 268 data samples are used to train the GRUbased model and the remaining 90 data samples are used to test the accuracy of the model. In Fig. 3
, we can see that the proposed model can accurately predict the illumination distribution after 100 epochs. Over 5000 independent runs, the average rootmeansquare error (RMSE) of training data prediction and test data prediction are 3.36 and 3.65, respectively. This is due to the fact that the GRUbased model can build a relationship between the prediction and the historical illumination distribution. Therefore, the GRUbased algorithm can predict the illumination distribution accurately.
Fig. 4 shows how the transmit power used to meet the users’ data rate and illumination requirements changes as the number of users varies. In this figure, we can see that the proposed algorithm achieves up to 22.1% gain in terms of transmit power reduction compared to a conventional optimal UAV deployment without considering the illumination distribution. This is due to the fact that the power required by the users is related to the illumination of the service area. From Fig. 4, we can see that the proposed algorithm is closer to the UAV deployment optimization using actual illumination distribution and the gap between the two schemes is less than 3.4%. This is because the proposed prediction algorithm can accurately predict the illumination distribution so as to optimize UAV deployment. Fig. 4 also shows that, as the number of users increases, the performance gain of the proposed deployment becomes less significant. This is because when enough users are considered, the users will be uniformly distributed in the square and the optimal position of the UAV will be fixed.
In Fig. 5, we show an example of how the proposed algorithm can optimize the deployment of the UAV. From Fig. 5, we can see that the optimal location of the UAV without considering the illumination distribution is the center of all the users and the optimal location of the UAV obtained by the proposed algorithm are shifted to the area with strong illumination. This is due to the fact that the illumination increases the interference for VLC link and, hence, the users located in a bright area need more transmit power compared to those located in a dark area. Under the collective effect of all the users in the service area, the optimal UAV locations move towards the area with strong illumination to minimize the total transmit power.
V Conclusion
In this paper, we have developed a novel UAV deployment framework for dynamically optimizing the locations of UAVs in a VLCenabled UAV based network. We have formulated an optimization problem that seeks to minimize the transmit power while meeting the illumination and communication requirements of each user. To solve this problem, we have developed a GRUbased prediction algorithm, which can model the longterm historical illumination distribution and predict the future illumination distribution. We have then shown that the location optimization problem is convex and the optimal solution is obtained by using the dual method. Simulation results have shown that the proposed approach yields significant power reduction compared to conventional approaches.
References
 [1] (2012Jan.) A componentwise EM algorithm for mixtures. Journal of Computational and Graphical Statistics. no.4 pp. 697–712. External Links: 1201.5913, Link Cited by: §IIIA.
 [2] (2019to appear,) Artificial neural networksbased machine learning for wireless networks: a tutorial. IEEE Communications Surveys Tutorials (). External Links: Document, ISSN 1553877X Cited by: §IIIB.
 [3] (201705) Caching in the sky: proactive deployment of cacheenabled unmanned aerial vehicles for optimized qualityofexperience. IEEE Journal on Selected Areas in Communications 35 (5), pp. 1046–1061. External Links: Document, ISSN 07338716 Cited by: §I.
 [4] (2019Mar.) Liquid state machine learning for resource and cache management in LTEU unmanned aerial vehicle (UAV) networks. IEEE Transactions on Wireless Communications 18 (3), pp. 1504–1517. External Links: Document, ISSN 15361276 Cited by: §I.

[5]
(2014Oct.)
Learning phrase representations using RNN encoderdecoder for statistical machine translation.
In
Proc. of Conference on Empirical Methods in Natural Language Processing
, Doha, Qatar. External Links: 1406.1078 Cited by: §IIIB.  [6] (201404) Energyefficient brightness control and data transmission for visible light communication. IEEE Photonics Technology Letters 26 (8), pp. 781–784. External Links: Document, ISSN 10411135 Cited by: §IIA.
 [7] (201905) UAV data collection over NOMA backscatter networks: UAV altitude and trajectory optimization. In Proc. of IEEE International Conference on Communications, Shanghai, China. External Links: 1902.03061, Link Cited by: §I.
 [8] (2018Sep.) UAV offloading: spectrum trading contract design for UAVassisted cellular networks. IEEE Transactions on Wireless Communications 17 (9), pp. 6093–6107. External Links: Document, ISSN 15361276 Cited by: §II.
 [9] (2004Feb.) Fundamental analysis for visiblelight communication system using LED lights. IEEE Transactions on Consumer Electronics 50 (1), pp. 100–107. External Links: Document, ISSN 00983063 Cited by: §II.
 [10] (201606) Unmanned aerial vehicle with underlaid devicetodevice communications: performance and tradeoffs. IEEE Transactions on Wireless Communications 15 (6), pp. 3949–3963. External Links: Document, ISSN 15361276 Cited by: §I.
 [11] NOAA/NCEIearth observation groupdefense meteorological satellite progam, Boulder. Note: https://www.ngdc.noaa.gov/eog/viirs/download_dnb_composites.html/ Cited by: §IV.
 [12] (2019) A vision of 6G wireless systems: Applications, trends, technologies, and open research problems. IEEE Network. Cited by: §I.
 [13] (2018Sep.) Capacity characterization of UAVenabled twouser broadcast channel. IEEE Journal on Selected Areas in Communications 36 (9), pp. 1955–1971. External Links: Document, ISSN 07338716 Cited by: §I.
 [14] (201905) Power efficient visible light communication (VLC) with unmanned aerial vehicles (UAVs). IEEE Communications Letters, pp. 1–4. External Links: Document, ISSN 10897798 Cited by: §I.
 [15] (2018Aug.) Joint altitude, beamwidth, location and bandwidth optimization for UAVenabled communications. IEEE Commun. Lett. 22 (8), pp. 1716–1719. External Links: ISSN 10897798 Cited by: §IIIC, §IIIC.
Comments
There are no comments yet.