Efficient 3D Positioning of UAVs and User Association Based on Hybrid PSO-K -Means Clustering Algorithm in Future Wireless Networks

Unmanned aerial vehicles (UAVs) play an important role in the future of 5G and 6G communication networks. UAV-assisted communication ofers the benefts of improved network capacity and coverage. A typical communication setup is for UAVs to connect users to the core network via a backhaul channel. Some of the challenges in such a setup include user-UAV association and management of the backhaul channel. Tese two challenges are greatly impacted by the positioning of the UAVs in the network. In this article, we address these challenges by considering a joint UAV placement and user association problem under data rate, signal to interference and noise ratio, and bandwidth constraints. To overcome this problem, a hybrid PSO -K -means clustering algorithm is used in two stages. In the frst stage, we use a K -means algorithm to cluster users and determine their horizontal locations. In the second stage, we use particle swarm optimization (PSO) to fnd the efcient 3D position of UAVs to maximize various network designs, namely, the network-centric approach and the user-centric approach. Te performance of the proposed solution is verifed using simulation results.


Introduction
Massive network capacity, low latency, and high data rates are some of the signifcant demands and goals of future wireless networks.Tese necessitate a practical transformation and efcient improvements in the current wireless communication network.Terefore, to meet these requirements, many technologies as promising candidates are considered for future wireless networks, including network densifcation, massive multiple-input and multiple-output (MIMO), and millimeter-wave (mmWave) communication.By operating millimeter-wave (mmWave) with large bandwidth where the frequency ranges from 30 to 300 GHz, we can enhance the data rate and capacity performance to a gigabit/s level [1][2][3][4].
Recently, it has become difcult to predict user trafc patterns within certain regions.Tis requires the rapid and massive deployment of ground stations.However, this is impossible due to expenditures, either capital or operational.In these scenarios, UAVs as base stations (BSs) have become a promising solution in future wireless networks.UAV-BSs can help a terrestrial BS network provide high data rate coverage at any time and place it is required [5][6][7].Also, UAV-BSs can deal with the issue of temporary coverage in rural and poorly populated areas and when terrestrial wireless infrastructure has been harmed by a natural catastrophe [8][9][10].
Along with the motivations mentioned above, there are numerous challenges in integrating the mmWave band with UAV-BSs.First of all, shorter wavelengths, for instance, cause smaller objects (i.e., humans) to become radio propagation obstructions in the line-of-sight (LoS) [11,12].Terefore, it is critical to take the blockage caused by the human body into consideration when analyzing the performance or when mmWave-BS deployment is being planned.On the contrary, the path loss (PL) in higher frequencies increases dramatically with the increasing distance between a transmitter and a receiver [13,14].At mmWave bands, it is another challenge.Consequently, there is a trade-of between locating a UAV at the maximum altitude (where the LoS improves) and maintaining a low PL (where PL increases with the increasing distance).Secondly, in contrast to the ground base stations (GBSs), UAVs have wireless backhaul between them and the core network.Consequently, as the data rates of access links between users and UAV-BS have been substantially increased, the capacity of the backhaul link becomes relatively limited [15,16].Terefore, there is a need to consider the constraints of the wireless backhaul link as one of the limitations in the design and deployment of the UAV-BSs.
Tirdly, a few articles have focused on the user-UAV association and UAV positioning issues [17][18][19][20][21].In the context of a single UAV, the authors of [17] examine UAV 3D location to increase the number of served users.Te association problem between the users and the UAV was addressed by using the signal-to-noise ratio (SNR).In [18], the problem of efcient 3D positioning of the UAV was examined.Two metrics were used: user-centric and network-centric to maximize the sum rate and the number of served users.Also, they imposed constraints including backhaul data rates and bandwidth and addressed the association problem by maximum path loss.Nevertheless, in the previous works [17,18], the authors employed exhaustive search algorithms to solve the association problem.Te algorithms are computationally expensive and thus impractical to employ.In the context of multiple UAVs, the authors in [19] seek to determine the least number of UAVs needed to service a set of users with high data rate requirements while placing multiple UAVs in 3D.Te users' association with UAV depends on the best signal-to-interference-plus-noise ratio (SINR) value.It is worth noting that the authors used a particle swarm optimization (PSO) algorithm to solve this problem.Nevertheless, they still sufer from high computational costs.Te authors in [20] developed and studied the association problem with UAVhubs using a simple greedy algorithm by taking into account the following constraints: the backhaul data rate of the link between the core network and the mother-UAV hub; the maximum bandwidth of each UAV-hub available for smallcell base stations; the maximum number of links that every UAV-hub can support; and the minimum SINR.In [21], authors employed an unsupervised learning-based k-means clustering technique for deploying, followed by the association of the small-cell base stations (SCBSs) by taking into account the same constraints that were mentioned in [20].Nonetheless, one of the limitations of this work is the fact that the UAVs are all assumed to be the same height, which might not be practically possible or advisable.
Finally, quite a few schemes [22][23][24] have been proposed to address the resource allocation problem in UAV-assisted communication.In [22], the authors employ two sleep scheduling policies for massive machine-type communication devices, which are defned by three diferent multiple access protocols.Moreover, they develop closed-form formulas for the massive machine-type communication devices' peak age of information (AoI), which are formulated as the optimization objective under the constraints of energy harvesting power, status update rate, and stability conditions.For the underwater data collection to be efciently completed, a heterogeneous autonomous underwater vehicle (AUV) auxiliary information collection system was proposed in [23].Te AUV trajectory with low time complexity was obtained by PSO.Additionally, to iteratively balance the trade-of between energy efciency and system queue backlog, a two-stage joint optimization algorithm based on Lyapunov optimization was created.By jointly optimizing the UAV trajectory and radio resource allocation for multiple access techniques, the authors of [24] characterized the capacity region of a UAV for multiple users.Te authors demonstrated that nonorthogonal multiple access signifcantly outperforms orthogonal multiple access in terms of rate regions for multiple users, while frequency division multiple access achieves higher rate regions than time-division multiple access.
To address these issues, we investigate the efcient deployment of a mmWave-UAV-BS by taking into account the properties of mmWave communication, where a human body may be able to block the LoS link by assuming UAVs are hovering in the quasi-stationary form at a specifc height and the height of the mmWave-UAV-BS is comparable to the altitude of the BSs mounted on the building walls [25].In our optimization problem, we utilize the terrestrial and airto-ground (A2G) channel models [26].In addition, we investigate optimizing efcient 3D positioning of UAVs and users-UAVs association to maximize the number of served users for various network designs by considering multiple communication-related constraints.To cope with this optimization problem, UAV positions are determined using a hybrid PSO-K-means clustering algorithm in two stages.In addition, a greedy algorithm with low complexity is proposed to solve the user-UAV association problem.
Tese paper's contributions can be summarized as follows: (i) To account for human body blockage, we adopt a terrestrial mmWave channel model for air-toground mmWave communication.(ii) To solve the user-UAV association problem, we propose an efcient greedy algorithm for the association problem by using multiple communication-related factors, including the maximum backhaul data rate, the maximum bandwidth of each UAV, and the maximum number of links that every UAV can support.(iii) We propose a joint user-UAV of data rate assignment optimization problem.Te rest of this article is structured as follows.In Section 2, the system model is presented.In Section 3, the optimization problem is formulated.In Section 4, the proposed approach is described.After that, the simulation results are presented in Section 5. Finally, the conclusions are shown in Section 6.

System Model
In our model, we assumed a downlink wireless heterogeneous network (HetNet) area that consists of three major nodes, namely, users, UAVs, and a core network gateway, as shown in Figure 1.We assume a crowded temporary event.Terefore, we will have ofoaded users due to the increase in their numbers or since their data rate exceeds the determined capacity.
Let I be the set of all ofoaded users where the total number of elements of the set is defned by N where i ∈ I � 1, 2, . . ., N { }.Te location of the i th user being r i � x i , y i  .Tese users, who are considered human blockers, are designed as cylinder shapes.It is worth mentioning that we assume all users are blockers to each other, where the mobile phone is located at the altitude h i .On the other hand, UAVs are denoted as J where M denotes the total number of elements of the set, where j ∈ J � 1, 2, . . ., M { }.We suppose that the control information, including backhaul data rate, available bandwidth, and SINR, is shared by all UAVs.In the examined model, each user is associated with their own UAV.In addition, the wireless backhaul link between each UAV and the core network gateway is a mmWave link.In our assumption, the proposed system model is assumed to be static for the time of operation.
To communicate between the UAVs and users, we adopted the A2G standard linear model.Te A2G model has two components: line-of-sight (LoS) and non-line-of-sight (NLoS) communication links at mmWave frequencies.For a given j th UAV and i th user, the standard linear PL models for two components are modeled as follows [26]: where α and β are the foating intercept and the line slope parameters, respectively.Te 3D distance between the UAV and user D ij is expressed by To account for the blockage in the human body, the probability of LoS was adopted [27], P LoS , and given as Also, P NLoS � 1 − P LoS , where λ and g B are the density and the diameter of human blockers, respectively.Te parameters denote h j , h i , and h B represent the height of UAV, user, and human blocker, respectively, and the 2D distance between user and UAV s ij is given as Ten, the average PL between i th user and j th UAV, in decibel (dB), is expressed by We assume P j is the transmit power for each UAV.Tus, when the j th UAV transmits a frame, the user i th receives the frame with power P j g ij L ij , where g ij denotes multipath fading and is assumed to follow an exponential distribution with a mean of μ.Ten, the SINR value is given as follows: where L ij � 10 −L dB ij /10 and σ 2 is the power of additive white Gaussian noise.

Problem Formulation
As we assumed, the network is in a download case where the data download from the ground core network getaway to N users by M UAVs.Furthermore, we suppose that the users in Cartesian coordinates are distributed randomly.Tere is a need to estimate the required number of UAVs such that all users can be handled.In order to accomplish this, each UAV can handle the maximum number of links up to L j from a Mobile Information Systems total number of N users.Consequently, the required M can be expressed by where ⌈.⌉ is the ceiling function.Let Q be a user-UAV association matrix.Te entries of this matrix are denoted as q i,j , where the row and column indicate the user and UAV indexes, respectively.Tus, q i,j is given as Assuming the desired data rate of users from the UAV is defned by r ij , this implies that the required bandwidth of users from the UAV is given as Te objective is to obtain efcient 3D positioning with a predefned efcient number of UAVs to maximize various network design parameters.Tese parameters include the sum rate and the number of served users for various rate demands of the overall network.To this end, the optimization problem is formulated as follows: Subject to: SINR ij q ij ≥ SINR min , ∀i, j, ifq ij ≠ 0, (10d) Constraint (10a) ensures that the total data rates of all UAVs do not exceed the data rate allowed by the backhaul channel R j .Tis is used as a QoS requirement.Constraint (10b) forces the total bandwidth allocated by each UAV not to exceed the allowable bandwidth for that particular UAV B j .Constraint (10c) limits on the number of users associated with each UAV not to exceed a certain number L j .One reason to justify such a requirement is to limit any possible interference, thus providing a better QoS.Constraint (10d) ensures that the SINR between a user and their UAV does not fall below a certain threshold given by SINR min .Tis implies that users with severely degraded channels are not included in the association process.Constraints (10e)-(10g) specify the limits of the spatial dimension of the UAVs.Constraint (10h) indicates the binary nature of the association problem.
For the cost function in (10), we consider the following scenarios: (i) Scenario 1: α ij � r ij in this scenario, the setup is called the user-centric approach, and the purpose is to maximize the total sum rate of the network (ii) Scenario 2: α ij � 1 in this scenario, the setup is called the network-centric approach, and the purpose is to maximize the number of served users of the network.

Hybrid PSO-K-Means Clustering Algorithm
Te objective function in ( 10) is a binary optimization problem.Such problems are shown to be NP-hard.Tis implies the high computational complexity of providing an optimal solution.Tis requires the development of efcient solutions that have low computational complexity.A twostage solution is presented to deal with high computational complexity for efcient UAV positions.In the frst stage, we clustered the ground users and determined their 2D location by using the k-means clustering algorithm.In the second stage, 3D UAV positions are determined using the PSO algorithm to maximize various network designs, whether they are a user-centric approach or a network-centric approach.

K-Means Clustering Algorithm. Te k-means clustering
algorithm is an unsupervised algorithm that was proposed by Stuart Lloyd of Bell Labs in 1957 [28].Ṫhe algorithm attempts to divide I users into M predefned clusters, with each user being assigned to only one group.K-means is also a measure of intracluster similarity and intercluster dissimilarity that is maximized through iterative hill-climbing.In our context, the k-means clustering algorithm deals with the 2D placement of the ground users.Algorithm 1 illustrates the k-means clustering method's pseudo-code, which was inspired by the k-means algorithm [29].
As we mentioned before, there is a need to divide a set of users I into M clusters, where each user belongs to the cluster with the closest mean.According to Algorithm 1, selecting the number of clusters (i.e., M) is the key stage in our algorithm before randomly initializing clusters centroids (i.e., step 5).Te algorithm performs two tasks throughout each iteration: (1) stage of cluster assignment and (2) stage of centroids moving (step 6).In the stage of assigning clusters, each point is examined by the algorithm, which selects the 4 Mobile Information Systems nearest centroid and assigns the point to it.In the moving centroids stage, the algorithm determines the cluster's mean point and shifts the centroids there.Tis two-step procedure is repeated until the algorithm converges.When the assignments stop changing, the algorithm converges.Te temporal complexity of the k-means clustering technique is O(NM), where N is the number of users and M denotes the number of clusters.Based on the suggested algorithm, we assume that each cluster will only have one UAV covered.
After it has fnished clustering the users, it uses the PSO algorithm to determine the location of the UAV in three dimensions to maximize the overall network total sum rate and the number of served users based on the applied approach.
4.2.Efcient 3D Positioning of UAVs.One of the most popular intelligence approaches and evolutionary global optimization techniques is the particle swarm optimization (PSO) algorithm.Te algorithm was frst proposed in [30] and was inspired by social animal group behavior including bird focks and fsh swarms.Tese sociable creatures are well-known for cooperating to maximize their accessibility to food by exchanging information among themselves.Te fnding sparked the idea of random agents collaborating to fnd better solutions [30].In PSO, each individual agent, also called a particle, changes its search velocity based on its best performance as well as the best performance of its group.In this stage, we apply the PSO algorithm to fnd an efcient 3D position of the UAV to maximize the objective function.Specifc steps are illustrated as follows: (a) Initialization.Te maximum number of iterations is determined in this frst stage, along with initializing the population position and setting the maximum particle speed.Te position information is used to determine the search space, and initialization is done by selecting random points within the speed range and search space.Te initial fying speed of each particle is then chosen at random, with M designating the size of the particle swarm.Te beginning population is set as the initial position of the UAV because the objective of this approach is to fnd the most efcient UAV 3D position.(b) Function of Fitness.Te objective function must satisfy multiple communication-related factors, including ( 1)-( 4) constraints in order to maximize the objective function, whether it is the total sum rate or the number of served users, in a user-centric approach or a network-centric approach, respectively.Consequently, the objective function can be used as the evaluation function as follow max x j ,y j ,h j ,q ij  M j�1  N i�1 α ij q ij .(c) Update Position and Speed.Following that, the ftness value of each particle is computed using the ftness function in an iterative manner to improve the candidate solutions based on two best values, i.e., the best local experience of each particle ω(i) best , and the best global experience ω(i) globalbest .Te best local location for each particle ω(i) best and the best global location ω(i) globalbest are updated in each iteration, and based on them, the particle positions and velocities are determined [19].How much the position can change is indicated by the speed V(i) value, which is provided by where rand(V size ) is uniformly distributed random variables with independent distributions in the range [0, 1], c 1 and c 2 are positive parameters known as acceleration coefcients, that determine the maximum step size between iterations, and w is the inertia weight.Additionally, each particle position ω(i) is updated as follows: Te computational cost of the algorithm is determined by the number of iterations (T max ), candidate solutions (N pop ), and iterations used to update each particle's velocity and position (ω).As a result, O(ωT max N pop ) gives the PSO algorithm's worst-case complexity.In comparison to other meta-heuristic algorithms (i.e., Artifcial Bees Colony (ABC), Genetic Algorithm (GA), and Ant Colony), the PSO algorithm is less complicated and takes less time to run.
Having determined the location of the UAVs via the PSO-K-means deployment scheme, we turn our attention towards the UAV-User association.Te association problem is an NP-hard binary problem that is computationally intensive.To solve this problem, we propose a greedy approach, which consists of two stages that are described below.
(i) SINR Ranking.Te purpose of this step is for UAVs to estimate the SINR values at the user side.Tis is necessary since the SINR values are not known beforehand and are sensitive to channel efects.To produce this estimate, the UAV initializes communication with a predetermined signal to the users.Each user calculates the SINR value at their side.Next, a user orders the M SINR values corresponding to M UAV and chooses the user with the maximum SINR value.We note that the user excludes all UAVs with corresponding SINR values that violate the minimum SINR values set in the constraint (A4).Upon determining the UAV with maximum SINR value, say the j th UAV, it sends confrmation feedback to that UAV.Combining all Mobile Information Systems user feedback to the UAVs, an association matrix Q can be constructed.Te Q matrix is a binary N × M matrix in which the (i, j) th entry corresponds to the association between the i th user and j th UAV.For i th user, if the j th UAV is the UAV with the max SINR then Q(i, j) � 1, otherwise it is set at 0. (ii) Q-Matrix Pruning.
We note that the association matrix Q established in the previous step is only concerned with ensuring that the SINR requirement is met.In order to satisfy the other constraints simultaneously with maximizing the objective function in (10), the association matrix Q needs to be pruned.For pruning the Q matrix, each UAV attempts to fnd the users who are willing to provide service from the set of users who send feedback requests to UAVs.In other words, for the j th UAV, the goal is to prune the j th column (i.e., column Q(: , m)).
To achieve this, we propose a greedy association algorithm shown in Algorithm 2. In this algorithm, the association is locally decided at each UAV, where the decision process is sequential in nature.We consider two cases; the user-centric scenario and the network-centric scenario objectives.
(i) Scenario 1 (User-Centric Scenario).In this scenario, the goal is to maximize the sum of data rates served by the UAV.Tus, the objective function in (10) is such that α ij � r ij .To maximize the objective, the proposed algorithm associates users to each UAV in a sequential greedy fashion.Te idea is to select the users with the maximum data rates and keep adding them to an association vector as long as the constraints (i.e., (A 1 ), (A 2 ) and (A 3 )) are not violated.
Let I m denote the association vector of the j th UAV.Te vector I m is an (N m × 1) vector with N n denoting the number of non-zero entries in the j th column of Q denoted as Q(: , m).Furthermore, let 1 q,m denote the vector of user indexes corresponding to nonzero entries in the column Q(: , m).Moreover, let r → users denote the set of data rate required corresponding to the user 1 q,m .We are now ready to explain the steps of the proposed algorithm.Te frst step is to fnd the maximum rate in the vector r → users .
Assume that it corresponds to the index i th , then the frst user to be associated is the i th user.Next, the association vector I m is updated to have a non-zero entry in the index i th .Te constraints are then modifed and checked for binary validity.Te association process continues until at least the constraints are violated.(ii) Scenario 2 (Network-Centric Scenario).
Tis is similar to the previous scenario, except for the fact that the goal here is to maximize the number of users associated with each UAV.To achieve this, the α ij in the cost function is modifed such that α ij � 1. Tis combined with the rate constraint implies that users should be added in a way that gives more priority to users with low data rates.Tis is in contrast to the user-centric approach explained earlier.Tus, we can employ the same greedy algorithm described where above but with the following modifcation instead of the descending ordering used earlier.Te entity of the r → users vector is selected in ascending order.

Simulation Results and Discussion
We suppose the geographic area is 1 (km) 2 in which 100 users are uniformly distributed in the Cartesian coordinates.
Each UAV can accommodate 30 requests.According to (7), we have 4 UAVs.Furthermore, we assume that each UAV fies at a low altitude of [100-400] m.Now we allocate the data rates at random from the vector r → users and use the procedure described in Algorithm 2 to obtain the results.Table 1 lists the simulation parameters.Te results are simulated for two diferent scenarios: (i) Scenario 1: α ij � r ij where the priority is given to users who have the maximum data rate and minimum bandwidth.(ii) Scenario 2: α ij � 1 where the priority is given to users who have the minimum data rate and minimum bandwidth.
As mentioned before, the network selects the user based on α ij .Te network-centric approach attempts to serve as many users as feasible, regardless of the required data rate.When using a network-centric, most of the serviced users For every ground user i ∈ I, set ϵ i � arg max j∈M ‖r i − μ j ‖ 2 For each cluster j ∈ M, set μ j � ( i∈I,ϵ i �j r i / i∈I,ϵ i �j 1) (7) Output: (8) |M| Clusters.ALGORITHM 1:K-means clustering algorithm.6 Mobile Information Systems have a low data rate.5 G networks require a distinction between the applications and users; thus, ofering the service to users with a low data rate is not fair.From this point, we used a user-centric approach to identify user priorities, where the network selects the user to maximize the total sum rate and prioritize the users with the highest data rate.
In the two proposed scenarios, the results of Algorithm 2 are examined for both network designs.Figure 2(a) and 2(b) depict the users' distribution and the efcient 3D positioning of UAVs for user-centric and network-centric approaches, respectively.In both approaches, we observed that the UAVs move to the maximum potential altitude to maximize the Step 1 (2) for i � 1 to Ndo (3) Create a list of UAVs that satisfy A 4 constraint (4) end for (5) Step 2 (6) for j � 1 to Mdo (7) ifScenario 1then (8) choose users having maximum data rate (i.e., max.r ij ) ( 9) end (10) else ifScenario 2then (11) choose users having minimum data rate (i.e., min.r ij ) ( 12) end (13) Counters Initialize: (17) end if (18) end while (19) end for ALGORITHM 2: Users association with UAVs.Mobile Information Systems objective function.Also, as we can see, in the user-centric approach, the service will be provided to fewer users compared to the network-centric approach.Most service providers tend to choose the network-centric approach, which is considered the most favorable option in that it requires less payment for using the spectrum.Te license fee depends on the amount of bandwidth consumed per user across a geographic area.Table 2 illustrates the numerical performance of Figure 2 and compares the two proposed approaches (i.e., user-centric and network-centric) in terms of total associated users, total bandwidth consumption, total sum rate, and percentage of users in an outage.Overall, we note that network-centric is superior to user-centric in terms of the number of connected users and the percentage of users experiencing an outage from the service.Tis is due to the priority that was previously explained.
To have a better understanding of the proposed algorithm and limitations, we run a few experiments in the following to study the efects of various limitations due to constraints (A1), (A2), and (A2) on the objective function equation (11), i.e., total sum rate and number of served users of the overall system.Furthermore, we discuss the impact of the constraints on the number of associated users.

Backhaul Data Rate Efect.
To study the backhaul data rate efect (i.e., (10a) constraint), Figure 3 illustrates the total sum data rate of the served users versus diferent ranges of R j , which is the peak data rate backhaul limit to the sum of the demanded data rate by the L j users for a single value of R j .For the user-centric and network-centric approaches, the ratio R j is kept the same by changing the backhaul data rate limit according to the demanded sum data rate of users.Additionally, by ofering more resources than necessary, other restrictions like those on bandwidth and the number of links are relaxed.Tis is done so that the impact of the backhaul data rate may be easily recognized.We can notice that the sum data rate increases with the increase in ratio R j until the ratio reaches R j � 0.8 Gbps.Ten, the sum data rate remains the same even with the increase in R j because, beyond R j � 0.8 Gbps, the algorithm has already associated all users.Tus, providing extra resources is unnecessary.
Figure 4 illustrates the efect of the R j constraint on the number of served users in the two approaches.As shown in the fgure, we severely limited the total number of served users due to the low backhaul data rate.Also, with the increased data rate backhaul per UAV in the two approaches, the total number of served users increases as well.At a certain data rate for backhaul (i.e., 1 Gbps), the increase in the overall number of users serviced is nearly constant.Tis is attributed to bandwidth resource exhaustion.Tus, the UAVs cannot serve more users.
In the user-centric approach, the rate of rising in the number of serviced users is practically constant (note the fxed slope of the blue dashed line), while in the networkcentric approach, it declines (take note of the red dashed lines' descending slope).Tis fxed slope is due to the proposed scenario where the users with a high data rate are served frst.Terefore, the backhaul capacity will increase, and users with a low data rate will receive service.However, in the network-centric approach, the slope is not constant but decreasing.In this scenario, the users with low data rate are served frst.Consequently, a few users with a high data rate obtain service, and in each step, the amount of increment in backhaul capacity is decreased by increasing the backhaul capacity.

Bandwidth Limit Efect.
To study the efect of bandwidth constraint (i.e., (10b) constraint), Figure 5 illustrates the total achievable sum rate versus diferent ranges of B j .Te two approaches show that the sum rate increases exponentially with the increase in bandwidth for each UAV, which is intuitive as users will be able to associate with more users due to increased bandwidth.We also note that the performance is almost similar in both approaches.Te total sum rate remains the same at a certain bandwidth value (i.e., 0.8 GHz).Tis is since the data rate backhaul is exhausted and all users have already associated.Peak Rate Backhaul (bps) User-centric approach Network-centric approach  8 Mobile Information Systems Figure 6 illustrates the efect of the B j constraint on the number of associated users in the two proposed approaches.By increasing the available bandwidth of each UAV, the total number of served users increases in the two approaches.Also, the network-centric approach can be noticed that the maximum number of served users is achievable.At a particular bandwidth value (i.e., 1 GHz), the increase in the total number of served users remains the same.Tis is attributed to the high data rate of backhaul exhaustion.Tus, the UAVs cannot serve more users.As for the slope, it bears the same interpretation that was previously illustrated in Figure 4.

Number of Links Limit Efect.
To study the efect of the number of links limit (i.e., (10c) constraint), Figure 7 illustrates the performance of the total sum rate when the number of connections accommodated by UAVs is raised from 5 to 40 links.Te backhaul data rate and bandwidth constraints are provided such that they do not afect the association of users.As a result, the infuence of the number of links can be seen entirely.Keep in mind that the number of UAVs remains the same (i.e., M � 4).With the increase of L j , the achieved sum rate by user-centric is greater than network-centric.Tis is due to the priority proposed in our greedy algorithms.Furthermore, in the two proposed approaches, the sum-rate value does not increase at a specifc value (i.e., L j � 30) because all the users are served at that point.User-centric approach Network-centric approach     and 8(e).Te UAVs adjust their positions dynamically, starting from their initial points (uncolored circles), which are chosen at random, and moving in the direction of the users they are serving in a few steps until arriving at their fnal efcient places (colorful stars).Te initial (uncolored circles) and fnal (colored circles) altitudes of UAVs are clearly illustrated in Figures 8(c) and 8(f), where these heights are plotted versus UAVs in X-coordinates.Te result demonstrates how UAVs modify their heights to alleviate interference, serve as many users as possible, and improve the LoS probability between the users and UAVs.

Conclusion
In summary, this paper presented the efcient 3D position of UAVs to provide mmWave backhaul connectivity to ground users with diferent rate requirements with the ground core network, using a hybrid PSO-K-means algorithm.Te user-UAV association is constrained by factors such as backhaul data rate, bandwidth available, the number of requests a UAV can accommodate, and the minimum SINR criterion.Furthermore, the objective function is done based on two approaches, namely, the user-centric approach that maximizes users' sum rate by giving priority to users having a maximum data rate and the network-centric approach that maximizes the total number of served users having a minimum data rate.Simultaneously, consuming less bandwidth in the two approaches is considered since the bandwidth resources will be fnite in the 5G network.Furthermore, we ran a few experiments to study the efects of various limitations due to constraints on peak data rate backhaul and bandwidth on the objective function.In addition, we discuss the efect of the constraints on the number of associated users.

2
Mobile Information Systems (iv) To solve the efcient 3D positions of UAVs problem that maximizes the objective function of the overall system, we propose a hybrid PSO-K-means clustering algorithm.(v) To solve the backhaul data rate assignment problem, we formulate the data assignment problem as a binary optimization problem.Te goal of this problem is to maximize the objective function by using multiple communication-related constraints.

( 1 )
Input: (2) Te placements of |I| ground users.(3) Te number of clusters |M|.(4) Start (5) Initialize the locations of the centroids μ 1 , μ 2 , . . ..μ |M| randomly (6) Repeat until convergence o o r d i n a t e ( m e t e r s ) X -c o o r d in a te ( m e te r s ) o o r d i n a t e ( m e t e r s ) X -c o o r d in a te ( m e te r s )

Figure 2 :
Figure 2: 3D location of UAVs and user-UAVs association in (a) user-centric approach and (b) network-centric approach.

Figure 3 :
Figure 3: Total sum rate of served users versus diferent backhaul data rates with restrictions B j � 5 GHz.

Figure 5 :
Figure 5: Total sum data rate of users versus available bandwidth limit with constraint R j � 5 Gbps.

Figure 4 :
Figure 4: Te total number of served users versus diferent backhaul data rates with restrictions B j � 5 GHz.

Figure 6 :
Figure 6: Te number of served users versus the bandwidth limit for the association algorithm with constraint R j � 5 Gbps.

Figure 7 :
Figure 7: Evaluation of total sum rate with diferent number of links constraint: B j � 5 GHz and R j � 5 Gbps. )

Table 2 :
Comparison of numerical results.