Two-Layer Optimization Algorithm for Multi-UAV Conflict Resolution considering Individual Fairness

In order to solve the unfair individual payment costs problem in the low-altitude unmanned aerial vehicle (UAV) conflict resolution process, a multi-UAV conflict resolution algorithm based on the cooperative game concept “coalition complaint value” is proposed. Firstly, based on the low-altitude multi-UAV conflict scene characteristics, according to the “coalition complaint value” concept, the UAV conflict resolution payment matrix is established. Secondly, combined with the advantages of the artificial potential field (APF) method and the genetic algorithm (GA), a hybrid solution strategy for conflict resolution based on APF-GA is proposed. The final simulation results show that the APF-GA hybrid solution strategy has the best efficiency by combining the three evaluation indicators of calculation time, feasibility, and system efficiency. The reliability of the proposed algorithm is verified based on the Monte Carlo algorithm. The solution strategy based on the cooperative game “coalition complaint value” can improve individual fairness to a certain extent. At the same time, it can achieve the rapid planning goal with priority drones at the expense of a small amount of overall benefits.


Introduction
With the development of aviation technology, UAVs have gradually entered the low-altitude fusion airspace due to their strong maneuverability, convenient operation, and low cost. And they have been widely used in many fields such as military, agriculture, logistics, and fire control [1]. Currently, isolated airspace for UAVs and strict approval and control are generally adopted to reduce operational risks. However, as the scale of UAVs continues to increase, it brings huge challenges to UAV air traffic management. Therefore, it is significant to study the low-altitude multi-UAV conflict resolution technology under free flight for the UAV safe operation and the low-altitude airspace resources efficient utilization. At the same time, after entering the low-altitude fusion airspace, the UAV can assist in the management of public safety (PS) by establishing elastic and reliable public safety networks (PSNs), having characteristics such as fast deployment, adaptive operation, coverage guarantees, and low latency [2].
The problem of multi-UAV conflict resolution is a typical complex system, which has the characteristics of large scale, nonlinear, multiconstrained, and high dimensional. It needs to find effective methods to simplify and solve. In recent years, many intelligent optimization algorithms and optimizers have been proposed to solve some practical and complex problems. Currently, the more mature methods for multi-UAV conflict resolution include intelligent algorithms, potential field method, and mathematical optimization methods. Based on the classic ant colony optimization (ACO) algorithm; Tang et al. [3] proposed a multi-UAV conflict resolution based on the speed adjustment strategy and the heading adjustment strategy, and it improved the computational efficiency. Perez-Carabaza et al. [4] optimized the ACO based on the target probability and spatial attributes to improve the convergence speed. Liu et al. [5] combined the APF method with the ACO. Firstly, the APF method was used to quickly search for the initial conflict resolution path, and then, the ACO was used for global optimization; it was both time-effective and feasible. Kim and Yoon [6] proposed a low-altitude unmanned airspace drone cooperative collision avoidance method based on satisfaction game. By combining self-preferences and cooperative preferences to make individual heading decisions, the satisfaction framework provides a strategy to avoid conflicts, and it can increase throughput while reducing unnecessary conflicts. Based on the mathematical optimization algorithm, Yang et al. [7] proposed a two-layer optimization of multi-UAV conflict resolution. First, the stochastic parallel gradient descent method was used to search for the initial solution, and, then the mixed integer linear programming model was used to accurately plan the path. Cai and Zhang [8] proposed a mixed-integer nonlinear programming conflict resolution method. It stipulates that the aircraft is allowed to change the speed and altitude to maintain the interval. Hernández-Romero et al. [9] proposed an aircraft conflict resolution method considering the wind forecast uncertainty influence, and they expressed the problem as a constrained nonlinear programming problem.
The conflict resolution method based on APF has the rapid planning characteristics. However, it often falls into the local optimum and cannot get the global optimum solution. It often generates paths with state mutations, which is not suitable for UAV performance constraints. As a new global optimization search algorithm, the GA algorithm is widely used in nonlinear system optimization problems due to its simple and universal, strong robustness, and suitability for parallel processing. However, due to the complex coding and random search of the initial population, the computational complexity is often high. In this paper, combining the advantages of APF and GA, a hybrid multi-UAV conflict resolution method is proposed. APF is used to search the initial solution quickly. GA is used to search the global optimization, so as to obtain the resolution method, which satisfies the UAV flight constraints.
At the same time, the current research on multi-UAV conflict resolution is limited to the optimal overall efficiency and ignores the fairness issues and priority of individual UAV resolution. In low-altitude airspace, the urgency of UAV missions is different, so it is necessary to consider the factors of different UAV priorities and individual fairness in the resolution process. Based on the cooperative game theory, this paper applies the concept of "coalition complaint value" to the conflict resolution field and studies the individual conflict resolution methods considering priority and fairness.
The contributions of this paper mainly include the following: (1) A multi-UAV conflict resolution strategy with better comprehensive performance was proposed (2) A welfare distribution strategy that considers individual fairness to ensure individual fairness in the UAV conflict resolution process was proposed (3) A strategy that can ensure the rapid arrival of UAVs with high priority was proposed

Conflict Resolution Strategy Based on the "Coalition Complaint Value" of Cooperative Game
Cooperative game is an important part of game theory. It is widely used in management, resource allocation, business, and other fields. Cooperative game has many important distribution solutions, such as nucleolus and Shapley value.
In the process of studying the cooperative game solution, the coalition complaint value is an important concept. It measures the coalition dissatisfaction with the distribution plan. The larger the coalition complaint value, the more dissatisfied the alliance is with the income distribution plan. Scholars have proposed some solutions to cooperative game based on the coalition complaint value. Tae et al. [10] and   International Journal of Aerospace Engineering Baïou and Barahona [11] first based on the coalition complaint value and obtained the prenucleolus and nucleolus of the cooperative game through the lexicographic method. In recent years, scholars have expanded the coalition complaint value and proposed other types of complaint values. Mu and Han [12] put forward the optimistic complaint value and pessimistic complaint value of the coalition. Sakawa and Nishizaki [13] defined the personal complaint value by summing the coalition complaint value where the player belongs. Kong et al. [14] defined the individual jealousy complaint value by considering both the coalition to which the player belongs and its complementary coalition. When multi-UAVs are engaged in conflict resolution in the low-altitude airspace, individual UAVs tend to choose the strategy beneficial to their own side for conflict resolution. In order to balance UAV individual fairness, based on the characteristics of low-altitude UAV conflict resolution and cooperative game, the "alliance complaint value" idea in cooperative game is introduced. The minimum "alliance complaint value" is used as a negotiation strategy between conflicting parties. The "alliance complaint value" essence in the cooperative game is to minimize the maximum coalition dissatisfaction value. The solution is unique. The following defines the "alliance complaint value" in conflict resolution.
In the cooperative game two-tuple hN, υi, the player set is N = fa 1 , a 2 ,⋯,a P g, and υ is the payment value [15,16]. In the grand coalition structure, S = fs 1 , s 2 ,⋯,s K g is the conflict resolution strategies set. For each resolution strategy s k ðk = 1, 2,⋯,KÞ, the conflict UAV cost is denoted as υðs k , a i Þ ðk = 1, 2,⋯,K, i = 1, 2,⋯,PÞ. Referring to the excess value Step L … Step 1 … Step L 3 International Journal of Aerospace Engineering definition of cooperative game, the new excess value definition is eðs k , a i Þ = υðs k , a i Þ. It represents the dissatisfaction level of the individual paying the price in the grand coalition. The closer the excess value is to 0, the smaller the individual dissatisfaction degree. The excess value eðs k , a i Þ of all conflicting individuals in the resolution strategy s k is arranged in order from large to small. The sequence is marked as εðs k Þ. Definition 1. There is a unique configuration γ in all the cooperative game grand coalition set N. And it meets the validity configuration, such that where γ is defined as the minimum "coalition complaint value" in the cooperative game. ≺ les is defined as "lexographically less." If y≺ les z, where y and z are both k-dimensional vectors, then l ∈ f1, ⋯, kg exists, such that y 1 = z 1 , ⋯, y l−1 = z l−1 , y l < z l is true. It should be noted that the cardinality of x and γ are not important. When the cardinality is changed to one dimension, the size of the two can be directly compared. It is easy to see that the nonempty minimum "coalition complaint value" exists in every cooperative game.

UAV Conflict Resolution Grand Coalition Payment
Matrix. For a grand coalition with the number of UAVs N = P, if there are K conflict resolution strategies, the grand coalition payment matrix V is expressed as In the formula, υðs k , a i Þðk = 1, 2,⋯,K, i = 1, 2,⋯,PÞ represents the resolution payment value of the ith UAV under the kth resolution strategy.  For the kth ðk = 1, 2,⋯,KÞ resolution strategy, arrange the resolution payment value in descending order to obtain the sequential payment matrix V ′ : where max υðs k , a i Þ is the maximum value of the kth row in matrix V, min υðs k , a i Þ is the minimum value of the kth row in matrix V. According to Definition 1, the conflict resolution "coalition complaint value" can be determined.

UAV.
In conflict UAV group ðA, υÞ, if there is a highpriority individuali * performing emergency tasks such as rescue and disaster relief, the UAV i * payment cost υðs k , a i * Þ is regarded as the largest in all conflict resolution solutions, that is, max υðs k , a i Þ = αυðs k , a i * Þðk = 1, 2,⋯,K, i = 1, 2,⋯,AÞ , α is the gain coefficient. Its value must satisfy the following conditions.
The UAV sequential payment matrix V ′ is expressed as That is to say, the UAV conflict resolution strategy must be produced in the scheme that has the priority UAV i * to pay the least cost, so as to ensure that UAV i * task is processed first and the task is completed first.
Here is a toy example to illustrate the selection process of the minimum "coalition complaint value" in the cooperative game.
Suppose there is 5 × 6 grand coalition payment matrix V, which means there are 5 conflict resolution strategies and 6 conflict individuals.     : ð7Þ According to the minimum "coalition complaint value" concept in the cooperative game, the 5th conflict resolution strategy is the optimal solution, and it can ensure the individuals fairness to the greatest extent.
If the 4th conflicting individual has the allocation priority, α is 8, and the sequential payment matrix V ′ becomes At this point, the optimal solution is the 2nd conflict resolution strategy.

Search Algorithm for Optimal Resolution
Strategy Based on APF-GA According to the ATC (air traffic control) safety regulations and the UAV actual flight situation, the UAV conflict resolution problem has been reasonably simplified.
(1) Considering that the UAV is flying at a fixed altitude during actual flight, except during take-off and landing, and taking into account the UAV fuel consumption, the UAV conflict resolution is simplified into two-dimensional plane (2) It is assumed that the UAV speed does not change significantly during the actual cruising flight, so it is assumed that the UAV speed remains unchanged. When the UAV executes the release strategy, three maneuvering strategies can be adopted, that is, maintaining the original heading, yawing 30°to the left and yawing 30°to the right (3) Assuming that UAVs are equipped with ADS-B surveillance equipment, the UAV position, speed, and altitude information at each time point can be obtained. Meanwhile, the UAV target position is known Assuming that the UAVs in the airspace are all agents, all drones can obtain local and surrounding drone status data through ADS-B. They can communicate with ground station, as shown in Figure 1. The low-altitude UAV system is a distributed multi-intelligence system. The conflict between UAV and other UAVs is a kind of dynamic  International Journal of Aerospace Engineering cooperative task solving problem. The process of proposing a conflict resolution mechanism is as follows.
Step 1: all UAVs with associated conflicts form a conflict group.
Step 2: the UAV in the conflicting group upload their respective status and attribute data to the ground station.
Step 3: the ground station forms a conflict resolution plan based on the negotiation strategy.
Step 4: plans are issued to each UAV and executed.
3.1. Artificial Potential Field Method. The APF method constructs an abstract artificial potential field through simulation to carry out trajectory planning. The APF is composed of a gravitational field, repulsive field, and vortex field. The gravitational field is generated by the target point position and acts on the UAV. The potential field vector direction is from the UAV to the target point, and the potential energy is inversely proportional to the distance to the target point. The repulsive field is generated by obstacles. The potential field vector direction points from the obstacle to the UAV, and the potential energy is inversely proportional to the distance to the obstacle. The vortex field is tangent to the repulsion field [17,18]. As shown in Figure 2, the position coordinate of UAV i is p i = ðx i , y i Þði = 1, 2, 3, ⋯, NÞ, and N is the UAV number in the area. For UAV i, the target position is p gi = ðx gi , y gi Þ. In order to ensure that the UAV is close to the target point, the gravity is defined as the negative gradient of the gravitational field. Then, the gravitational potential field function U a ðp i Þ and gravity F a ðp i Þ are expressed as In the formula, k ai is the gravitational constant. Its value is greater than 0.
In order to realize the collision problem when multi-UAVs are gathered and run, the repulsion field U r ðp i , p j Þ will be excited between the UAV i and UAV j. Define the repulsion force as the negative gradient of the repulsion field. Then, the repulsion potential field function U r ðp i , p j Þ and the repulsion force F r ðp i , p j Þ are

10
International Journal of Aerospace Engineering In the formula, k ri is the repulsive force constants. Its value is greater than 0. ρ 0 represents the range of the UAV's repulsive force field. When the UAV i is in the range of the UAV F's repulsive force field, it will be subjected to a repulsive force.
In order to prevent UAVs from randomly turning to generate a large number of redundant paths in the conflict resolution process, it is necessary to introduce rules to ensure that UAVs turn in unison, so a vortex field is introduced, and it is tangent to the repulsion field.
Superimpose the gravity, repulsion, and vortex forces experienced by the UAV to obtain the UAV operating equation in the conflict resolution process: Both the repulsive force and the vortex force are taken between 0 and 1 to control the UAV to avoid collisions. The gravity is normalized so that the force received by the UAV is kept at the same order of magnitude; so as to avoid the large difference between the forces, it may lead to the conflict resolution failure or the generation of redundant paths.
Normalize the UAV motion equation p i · to get the speed direction. Assuming the speed is k i , the UAV speed can be expressed as Using the UAV resolution speed v i , iterate gradually according to the time step, and obtain the multi-UAV conflict resolution trajectory.

Genetic
Algorithm. GA is a global optimization probability search algorithm formed by simulating the genetic and evolution process of living beings in the natural environment. It repeatedly acts on the entire evolutionary group through operation operators such as selection, crossover, and mutation. Finally, the optimal solution or approximate optimal solution of the problem is obtained [19][20][21].
3.2.1. Coding Method. Divide the path from the initial position to the target position of the N UAVs into L steps, respectively. Through analysis, the UAV has three steering strategies. If the UAV reaches point B from point A at time t, then, the reachable positions at time t + 1 are C 1 , C 2 , and C 3 . If the UAV current position and heading information matrix is denoted as A k ½x, y, θ, the next step UAV information matrix is denoted as A k+1 ½x + vt cos θ′, y + vt sin θ′, θ′. In the formula, θ′ can take θ + ðπ/6Þ,θ and θ − ðπ/6Þ three kinds of adjustment strategies.
If the three UAV maneuvering strategies are coded as -1 (representing a 30°left turn), 0 (representing a heading maintenance), and +1 (representing a right turn 30°), the Step number

12
International Journal of Aerospace Engineering UAV system conflict resolution strategy can be encoded into The specific encoding structure is shown in Table 1.

Objective Function and Constraints.
The UAV uses appropriate maneuvers to resolve conflicts in accordance with the prescribed coding strategy. At the end of the resolution process, the UAV reaches the shortest distance from the destination. Therefore, the fitness function for conflict resolution is set as In the formula, Delay i represents the ith UAV delay distance from the destination at the end of conflict resolution. The smaller the delay distance is, the better the conflict resolution effect. Denote the position of UAV i at the kth step as A k i , then, the UAV i delay distance Delay i is expressed as the Euclidean distance between the position A L i ðx L i , y L i Þ after L steps and the destination position B i ðx g , y g Þ; it is expressed as In order to avoid the UAV collision in the conflict resolution process, it is required that the distance between UAV iðx i , y i Þ andjðx j , y j Þ maintains a minimum safety interval at any time: In the formula, d represents the minimum safe interval among UAVs.

Selection, Crossover, and Mutation.
According to the fitness of each individual in the group, the roulette method is used to select the fittest survival among the individuals in the group. According to the crossover probability P C , two parent individuals are randomly selected from the group for single-point crossover to generate new individuals. Taking into account that the one-point crossover operator is less likely to destroy individuals with higher fitness values, this method is adopted to ensure the superiority of increasing the population. As shown in Figure 3. In the meantime perform gene flip mutation with probability P m .

Optimal Retention.
In order to retain the excellent individuals that have been obtained, the optimal retention is added to the selection operation. That is, the good individuals that have appeared in the previous generation are directly inherited to the next generation without any crossover and mutation.

Conditions for Genetic Termination.
When the number of inherited generations exceeds a predetermined value, the inheritance is terminated, and the individual with the greatest fitness in the final chromosome is output. Step number 14 International Journal of Aerospace Engineering

Multi-UAV Conflict Resolution Process Based on APF-GA Algorithm
Firstly, based on the APF method principle, a preliminary search is made for the initial path of multi-UAV conflict resolution. After the iteration is terminated, the initial path is used to initialize the population information. Finally, the global optimum is searched through constant iteration. The whole calculation process is shown in Figure 4. The specific implementation steps are described as follows: Step 1: initialize the parameters. Set the initial conflict position p 0 i = ðx 0 i , y 0 i Þ and target point p gi = ðx gi , y gi Þ for the UAVs in the area.
Step 2: force analysis of the potential field. Calculate the gravitational force F a ðp i , p gi Þ, repulsion force F r ðp i , p j Þ, and vortex force F v ðp i , p j Þ of each UAV and calculate the motion equation of each UAV according to formula (12).
Step 3: update location. According to the calculated UAV speed direction and given speed size, step length, and other information, iteratively update the position step by step.
Step 4: judge the result. To monitor in real time whether the UAV i has reached its destination p gi = ðx gi , y gi Þ, if it has reached the destination, the conflict resolution process is completed; if it has not been reached, it turns to Step 3 to update the position.
Step 5: the initial path generated by the APF method is approximated and encoded and set as the initial value of the GA search.
Step 6: apply the concept of "coalition complaint value" in cooperative game to the conflict resolution process, calculate the fitness value of each body in the group, and iterate the population information through the selection, crossover, mutation, and other operations in the GA algorithm, so as to obtain the optimal calculation result of the global area and corresponding conflict resolution path.
Step 7: if the specified termination condition is reached, the final conflict resolution path will be output; if not, turn to Step 6 to iterate until the requirements are met. Assuming that the conflict circle radius is 100 m, the initial paths are linear, as shown in Figure 5. By changing course, UAVs can reach their respective destinations while ensuring a safe interval between them. As an extreme case, classical conflict scenario can better test the merits and demerits of conflict resolution algorithm. For N conflict individuals, there are N × ðN − 1Þ /2 distance curves in the conflict resolution process. Figure 6, the UAVs maintain a certain minimum safety interval during the operation, and they converge at the intersection point at a certain moment, and the initial path is a straight line. The entire system conflict resolution is completed by adjusting the direction.

Scenario 2: UAV Flow Conflict Scene. As shown in
In this calculation example, it is assumed that each cluster contains 4 individual UAVs. The parameters required for simulation in the three conflict scenarios are shown in Table 2.

Algorithm Validation.
Taking the 2, 6, and 10 aircraftto-flight conflict scenario in scenario 1 as an example, two algorithms, APF and APF-GA, are used to solve the solution strategy. The results are shown in Figures 8 and 9.
In order to visually demonstrate the effectiveness of the algorithm proposed in this paper, the convergence curves of 2 UAVs, 6 UAVs, and 10 UAVs are shown in Figure 10. After analyzing the simulation phenomenon, it can be seen that the algorithms can converge quickly without oscillation. Figure 11 shows the running distance of the convergence curve of 2, 6, and 10 UAVs, respectively. There is ð2 × 1Þ/2 = 1 distance curve during the 2 UAVs' conflict resolution process, ð6 × 5Þ/2 = 15 distance curve exists during 6 UAVs' conflict resolution process, and ð10 × 9Þ/2 = 45 distance curve exists during 10 UAVs' conflict resolution process. According to the simulation data, it can be seen that the distance during the UAV operation is greater than the safety interval d = 20 m, indicating that the solved conflict resolution path meets the safety constraints. It should be noted that due to the symmetry of conflict scenes, curves in Figure 11 overlap and the number of lines remains unchanged.  In order to evaluate the effectiveness of the multi-UAV conflict resolution algorithm, three indicators of calculation time, feasibility, and system efficiency are now proposed to measure the resolution effectiveness.
In calculation time, the conflict resolution process needs to meet real-time requirements, so the shorter the calculation time is, the stronger the real-time performance is, and the better the resolution efficiency is. The calculation time expressed as In the formula, m represents the simulations number for conflict resolution, and T j represents the time required for the jth UAV to resolve.
For feasibility, it is used to measure whether the conflict resolution path satisfies the UAV flight. Generally speaking,   International Journal of Aerospace Engineering if the single heading change is too large, the path point is regarded as nonflyable. The greater the feasibility value is, the safer the conflict resolution process is. The feasibility is expressed as In the formula, f j represents the number of nonflyable UAV in the UAV j conflict resolution path.
For system efficiency, it is used to measure the delay degree in the conflict resolution process. If all UAVs run along the original path, the system efficiency is 100%. When the UAVs have conflicts, the system efficiency is reduced accordingly; it is expressed as In the formula, S gj and S j , respectively, represent the preset step length and delay of the jth UAV.
The UAV number is selected from 2 to 28. In the antiflight conflict scene, the solution strategy performance is verified by comparing calculation time, feasibility, and system efficiency. The reference algorithm includes APF, GA method, stochastic parallel gradient descent (SPGD), and a typical nonlinear optimization solver Snopt. As shown in Figure 12, in terms of computing time indicators CT, the APF-GA algorithm can maintain a calculation speed similar to APF, which is far lower than the GA algorithm, SPGD, and Snopt. In terms of feasibility indicators F, APF generates more path points with state mutations, which is not conducive to UAV flight. But APF-GA can absorb the GA advantages and can restrict the path points in a targeted manner to generate a flyable path. In terms of conflict resolution efficiency SE, with the increase of conflicting individuals, the efficiency of APF-GA lies between APF, GA. SPGD, and Snopt, and the efficiency declines slowly. Therefore, APF-GA can quickly generate a flyable path and can maintain a high conflict resolution efficiency. On the whole, the performance is better.

Reliability Analysis of Conflict Resolution Strategy
Based on Monte Carlo Algorithm. In the above simulation analysis, the conflicts of multi-UAVs are all successfully resolved. However, due to the complicated environment during the actual operation of the UAV, the above simulation is only based on the conflict situation where the position, speed, and heading are all fixed, and the universality of the algo-rithm cannot be guaranteed. Therefore, a large number of random UAV conflicts are generated by the Monte Carlo method to obtain the conflict probability of multi-UAVs. Similar to the above, this paper is still based on the simulation of small UAVs.
In order to evaluate the performance of the algorithm in this paper, some Monte Carlo simulation parameters need to be set. The performance of the algorithm can be evaluated by the collision probability P ct .
In the formula, N t is the number of UAV distance intervals in the area, which is N × ðN − 1Þ/2, and N c is the number of intervals that do not meet the safety interval during the conflict resolution process.
In order to simplify the simulation conditions, the UAVs in the Monte Carlo simulation are kept flying in a fixed square area of200 m × 200 m. The random situation of the four UAVs is shown in Figure 13. The initial position p i ðx i , y i Þ and the target position p gi ðx gi , y gi Þ of the UAV are randomly distributed in the square area A squ , and the range of the horizontal and vertical coordinates is ½−200 m, 200 m. Other simulation conditions are consistent with Section 5.1.
Taking the number of UAVs of 2-50 as the object, through the Monte Carlo simulation experiment, the trend of conflict probability P ct under different numbers of individuals is obtained. As shown in Figure 14, the conflict probability decreases with the increase of conflicting individuals. But the overall level is at a relatively high level, meeting the reliability index.

The Impact of Cooperative Game "Coalition Complaint
Value" and Traditional Cooperative Solution Strategy on Individual Fairness. Taking the UAV flow conflict scenario in scenario 2 as an example, APF-GA can quickly solve the resolution strategy. Based on the two solution strategies of cooperative game "coalition complaint value" and traditional cooperation (take the total costs paid by the individual as an indicator), two solution strategies are obtained. The individual UAV payment costs under the two solution strategies are shown in Table 3.
In order to more intuitively compare the two solution strategies impact of cooperative game "coalition complaint value" and traditional cooperation on individual fairness, draw the Lorentz curve, as shown in Figure 15.
Taking the UAV conflict scenario in the square area in scenario 3 as an example, APF-GA can quickly solve the resolution strategy. In the same way, the individual UAV  Figure 16.
According to Figures 15 and 16, the solution strategy based on the cooperative game "coalition complaint value" is fair. The curve is closer to the fairness baseline, and it indicates that the solution strategy of the cooperative game "coalition complaint value" can improve individual fairness to a certain extent. 5.5. UAV Priority Allocation. In order to more clearly compare the different priority allocation strategies impact on the cost of individual resolution payment, 4 aircrafts in scenario 1 are selected to simulate the flight conflict scene. The UAV priority deployment rule based on cooperative game "coalition complaint value" proposed in this paper is compared with the deployment method without considering priority. If UAV1 is performing emergency tasks such as rescue and disaster resolution, the priority is the highest (α takes 100), and the resolution effect is shown in Figure 17.
The individual payment costs under the two deployment strategies are shown in Table 4. Under the deployment strategy without considering priority, each UAV individual conflict status is symmetrical and equal, and the payment cost is basically the same. The priority deployment based on the cooperative game gain coefficient can reduce the UAV1 payment cost and promote the high-priority UAV to reach the destination quickly, but it will sacrifice the conflict system overall interests. Compared with the nonpriority deployment, the priority deployment increases 29.06%, but the payment cost of UAV1 decreases 50.26%. Therefore, it can be analyzed that priority deployment based on "coalition complaint value" can realize the UAV rapid planning with priority to reach the target at the premise of sacrificing a small amount of overall benefits.

Conclusion
Based on the cooperative game "coalition complaint value" concept, this paper proposes a low-altitude multi-UAV conflict resolution algorithm, and it realizes the conflict resolution under the premise of ensuring the individual resolution payment costs fairness.
(1) Considering the three conflict resolution evaluation indexes of calculation time, feasibility, and system efficiency, the multi-UAV conflict resolution algorithm based on APF-GA has better performance and can quickly generate a resolution path. And, it satisfies UAV performance constraints (2) The reliability of the proposed algorithm is verified based on the Monte Carlo algorithm. The conflict probability decreases with the increase of conflicting individuals, but the overall level is at a relatively high level, which meets the reliability index (3) The conflict resolution strategy based on the cooperative game "coalition complaint value" can improve individual fairness to a certain extent (4) Priority deployment based on "coalition complaint value" can achieve the rapid planning goal with priority UAV at a small amount of overall benefit expense Abbreviations N = fa 1 , a 2 ,⋯,a P g: Player set υðs k , a i Þðk = 1, 2,⋯,K , i = 1, 2,⋯,PÞ: Payment value S = fs 1 , s 2 ,⋯,s K g: Conflict resolution strategies set s k ðk = 1, 2,⋯,KÞ: Resolution strategy eðs k , a i Þ: New excess value εðs k Þ: The sequence when the resolution strategy s k is arranged in order from large to small γ: The minimum "coalition complaint value" in the cooperative game ≺ les : Lexographically less V: The grand coalition payment matrix α: Gain coefficient V ′: The sequential payment matrix i * : High priority individuals p i = ðx i , y i Þði = 1, 2, 3 , ⋯, NÞ: Position coordinate of UAV i p gi = ðx gi , y gi Þ: Target position of UAV i k ai : Gravitational constant U a ðp i Þ: Gravitational potential field function F a ðp i Þ: Gravity U r ðp i , p j Þ: Repulsion potential field function F r ðp i , p j Þ: Repulsion force k ri : Repulsive force constants ρ 0 : The range of the UAV's repulsive force field F v ðp i , p j Þ: Vortex force p i · : UAV operating equation k i : UAV speed size v i : UAV resolution speed L: Step number t: Time A k ½x, y, θ: UAV position and heading information matrix at time k A k+1 ½x + vt cos θ ′ , y + vt sin θ ′ , θ ′ : UAV position and heading information matrix at time k + 1 X = ½x 1 , x 2 ,⋯,x N×L : Conflict resolution path coding vector θ: Angle Delay i : The ith UAV delay distance from the destination at the end of conflict resolution f : Fitness function for conflict resolution A k i : The position of UAV i at the kth step d: The minimum safe interval P C : Crossover probability P m : Mutation probability N max : Iterations 20 International Journal of Aerospace Engineering

CL:
Calculation time m: The simulations number for conflict resolution T j : The time required for the jth UAV to resolve. F: Feasibility f j : The number of nonflyable UAV in the UAV j resolution path.

SE:
System efficiency S gj : The preset step length of the j-th UAV. S j : The delayof the jth UAV. P ct : Collision probability N t : The number of UAV distance intervals in the area N c : The number of intervals that do not meet the safety interval.

Data Availability
No data were used to support this study.

Conflicts of Interest
The authors declare that they have no conflicts of interest.