Next Article in Journal
Design and Implementation of Sigma-Delta ADC Filter
Previous Article in Journal
Performance of Millimeter Wave Dense Cellular Network Using Stretched Exponential Path Loss Model
Previous Article in Special Issue
Applications of Multi-Agent Systems in Unmanned Surface Vessels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ability-Restricted Indoor Reconnaissance Task Planning for Multiple UAVs

1
School of Automation, Beijing Institute of Technology, Beijing 100081, China
2
Peng Cheng Laboratory, Shenzhen 518055, China
*
Authors to whom correspondence should be addressed.
Electronics 2022, 11(24), 4227; https://doi.org/10.3390/electronics11244227
Submission received: 14 November 2022 / Revised: 10 December 2022 / Accepted: 11 December 2022 / Published: 19 December 2022
(This article belongs to the Special Issue Recent Advances in Intelligent Unmanned Systems)

Abstract

:
For indoor multi-task planning problems of small unmanned aerial vehicles (UAVs) with different abilities, task assignment and path planning play a crucial role. The multi-dimensional requirements of reconnaissance tasks bring great difficulties to the task execution of multi-UAV cooperation. Meanwhile, the complex internal environment of buildings has a great impact on the path planning of UAVs. In this paper, the ability-restricted indoor reconnaissance task-planning (ARIRTP) problem is solved by a bi-level problem-solving framework. In the upper level, an iterative search algorithm is used to solve the task assignment problem. According to the characteristics of the problem, a solution-space compression mechanism (SSCM) is proposed to exclude solutions that do not satisfy the task requirements. In the lower level, based on a topological map, the nearest neighbor (NN) algorithm is used to quickly construct the path sequence of a UAV. Finally, the genetic algorithm (GA) and simulated annealing (SA) algorithm are applied to the upper level of the framework as iterative search algorithms, which produces two hybrid algorithms named the GA-NN and SA-NN, respectively. ARIRTP instances of different scales are designed to verify the effectiveness of the SSCM and the performance of the GA-NN and SA-NN methods. It is demonstrated that the SSCM can significantly compress the solution space and effectively improve the performance of the algorithms. The proposed bi-level problem-solving framework provides a methodology for the cooperation of multi-UAV to perform reconnaissance tasks in indoor environments. The experimental results show that the GA-NN and SA-NN methods can quickly and efficiently solve the ARIRTP problem. The performance of the GA-NN method is similar to that of the SA-NN method. The GA-NN method runs slightly faster. In large-scale instances, the performance of the SA-NN method is slightly better than that of the GA-NN method.

1. Introduction

Due to the advantages of flexibility and mobility, unmanned aerial vehicles (UAVs) are increasingly used in various environments [1,2,3,4,5,6,7,8,9,10,11]. In mountainous environments, UAVs are widely used in searches, reconnaissance, surveillance, fault detection, data collection, target recognition, and classification [12,13,14,15,16]. In suburban environments, UAVs are widely used in agriculture such as for watering, sowing, and spraying pesticides [17]. In urban environments, UAVs are widely used in searches, reconnaissance, and other tasks. Recently, there has been an increasing demand for UAVs to solve complex indoor tasks such as target searches [18], surveillance, search and rescue [19], package delivery to large indoor facilities, and task-scheduling systems [20].
Compared with outdoor environments, such as urban environments and mountainous environments, indoor environments have more restrictions on UAVs [21]. In outdoor environments, there may be no-fly areas and dangerous areas because of bad weather or electromagnetic interference [22], which require UAVs to avoid those areas or pass through those areas in the shortest possible time. Obstacle avoidance is an important problem for UAVs in path planning. In contrast, such constraints in indoor environments may be different. Specifically, for mobility, UAVs cannot move forward smoothly in an indoor environment because of the narrow indoor space and safe flight distance of UAVs. In addition, compared with outdoor open environments, the connection relationships between indoor rooms are complex and have a serious impact on the path planning of UAVs.
In existing research on UAVs, the task assignment and path planning of UAVs are often studied separately [23,24]. In [23], an improved algorithm based on the simulated annealing algorithm (SA), namely the swap-and-judge simulated annealing (SJSA) algorithm is used to solve the multi-task planning problem of UAVs. However, in this problem, the task requirements and UAV abilities are single. The impact of environmental factors on UAV path planning is not specifically considered. In [25], a bi-level ant colony optimization (BACO) algorithm was proposed to solve the capacitated electric vehicle routing problem (CEVRP). The upper level determines the visiting sequence of the tasks, and the lower level adjusts the task sequence generated by the upper level according to the electricity constraint. In addition, most of the existing research on indoor path planning is to find a path that meets the requirements so that the UAV can smoothly reach the destination from the starting point [21,26]. Most studies aim to find an optimal path for one UAV. There are also some studies focusing on multi-UAV cooperative path planning but the process of task assignment is not considered [27,28]. In recent studies, most reconnaissance tasks have single requirements and the abilities of UAVs are also single. When the multi-task assignments of multiple UAVs are carried out, these problems can be modeled as the classical vehicle routing problem (VRP). However, when reconnaissance tasks have multi-dimensional requirements and the abilities of UAVs are heterogeneous, it may be necessary to use multiple UAVs with different abilities to complete a single reconnaissance task.
Motivated by the above observations, in order to balance the multi-task assignment and path planning of multiple UAVs, this paper considers the multi-dimensional requirements of reconnaissance tasks. Meanwhile, the influence of the connection relationships between rooms on the path planning of UAVs is considered based on the indoor topological map. This paper focuses on solving the ability-restricted indoor reconnaissance task-planning (ARIRTP) problem. The main contributions of this paper are as follows:
  • The ARIRTP problem is modeled as a combinatorial optimization problem in which the requirements of the reconnaissance tasks and the abilities of UAVs have multiple dimensions. A topological map is used to present the connection relationships between indoor rooms and help UAVs to conduct global path planning;
  • A bi-level problem-solving framework is proposed for the ARIRTP problem. The upper level uses an iterative search algorithm to solve the task assignment problem of UAVs. According to the characteristics of the problem, a solution space compression mechanism is proposed to make the generated task assignment schemes meet the requirements of the reconnaissance tasks. The lower level uses the nearest neighbor (NN) algorithm to quickly construct the path sequence of each UAV based on the topology information of buildings;
  • Two hybrid algorithms are proposed by applying the genetic algorithm (GA) and SA algorithm as the iterative search algorithms to the upper level of the bi-level problem-solving framework, respectively.
The remainder of this paper is organized as follows. The problem description and the mathematical model are described in Section 2. The details of the proposed algorithms are described in Section 3. The discussions of the design of the experiments and the experimental results are given in Section 4. Finally, Section 5 provides the conclusion and perspectives.

2. Problem Formulation

In this section, first, the problem description of ARIRTP is given. Then, the environment modeling process is described. Finally, the problem model is presented.

2.1. Problem Description

A group of UAVs equipped with various sensors is required to conduct a set of reconnaissance tasks with different requirements distributed in a building. In an indoor environment, multiple rooms are connected by doors. A series of reconnaissance tasks with multiple dimensions and different requirements are scattered indoors. Based on the requirements of the reconnaissance tasks, we judge whether the UAVs are able to complete the whole task or part of the task. If the UAVs cannot independently meet the requirements of all dimensions of a certain reconnaissance task, multiple UAVs need to be combined to make their common capabilities meet the reconnaissance requirements of the whole task.
Figure 1 shows the situation of multiple reconnaissance task planning for multiple UAVs in an indoor environment. Three UAVs are required to take off from the entrance point (start point) and enter the building to perform seven reconnaissance tasks. If one UAV can meet the requirements of all dimensions of a task, it can complete the task independently. For example, UAV U 1 independently performs reconnaissance tasks T 7 and T 6 . However, when the UAV cannot meet all the requirements of the task alone, it needs to complete the corresponding reconnaissance task cooperatively in a certain combination mode. For instance, UAV U 2 and UAV U 3 cooperate to conduct reconnaissance task T 1 because they cannot complete the task alone. When UAV U 2 and UAV U 3 pass through door D 4 , the two UAVs may collide. Now, UAV U 2 can choose to wait for UAV U 3 to pass through door D 4 first according to the priority of the task to avoid a collision. Topological maps represent the environment using graphs, where the vertices and edges represent the rooms and connection relationships between the rooms, respectively [29,30,31,32].
Due to ability restrictions, a variety of task assignment situations could occur. One UAV performs one reconnaissance task, one UAV performs multiple reconnaissance tasks, or multiple UAVs perform one reconnaissance task.

2.2. Environment Modeling

The contents of the indoor environment are modeled as follows: the index sets of rooms, doors, task points, and UAVs are denoted by R n , D n , T n , and U n , respectively. The coordinate matrices of the doors and reconnaissance task points are described by P D and P T , where x and y represent the horizontal and vertical coordinates of the different doors or task points.
R n = { R 1 , R 2 , , R N R } .
D n = { D 1 , D 2 , , D N D } .
T n = { T 1 , T 2 , , T N T } .
U n = { U 1 , U 2 , , U N U } .
P D = x D 1 y D 1 x D 2 y D 2 x D N D y D N D .
P T = x T 1 y T 1 x T 2 y T 2 x T N T y T N T .
The vector R indicates the room the task points are in. r j represents the index number of the room that task j belongs to.
R = ( r 1 , r 2 , , r j , , r N T ) , r j R .
The connectivity matrix between rooms is described by A. If two rooms share a door, the two rooms are connected.
A = a 1 , 1 a 1 , 2 a 1 , N R a i , j a N R , 1 a N R , 2 a N R , N R .
a i , j = 0 if i = j , i , j R , g if room i and room j are directly connected , ( g D ) , if room i and room j are indirectly connected .
The distance matrix between two doors is described by B. If the two doors are in the same room, the distance between the doors can be approximately represented by the Euclidean distance of the positions of the doors (represented by d i s t a n c e ( D h , D g ) ). If the two doors are not in the same room, the distance between the doors can be obtained based on the Dijkstra algorithm (represented by d i s t a n c e D i j ( D h , D g ) ). a r g m i n ( d i s t a n c e D i j ( D h , D g ) ) represents the coordinate matrix of the doors connecting door h and door g.
B = b 1 , 1 b 1 , 2 b 1 , N D b h , g b N D , 1 b N D , 2 b N D , N D .
b h , g = 0 if h = g , h , g D , d i s t a n c e ( D h , D g ) if door h and door g are in one same room , d i s t a n c e D i j ( D h , D g ) if door h and door g are not in one same room .
a r g m i n ( d i s t a n c e D i j ( D h , D g ) ) = x D h y D h x D g y D g .

2.3. Problem Model

In this section, the ARIRTP problem model is given in detail. First, the related assumptions are made. Then, the constraints and decision variables in the ARIRTP problem are represented. Finally, the objective function is described. The list of relevant parameters and variables is shown in Table 1.

2.3.1. Assumptions

The assumptions about the relevant situations in the ARIRTP problem are made from three aspects: the UAVs, reconnaissance tasks, and indoor environment [33,34,35,36].
  • UAVs
    (1) The curvature constraint of the UAVs is not considered because small-rotor UAVs are usually used in indoor environments and their flight speed is slow.
    (2) The collision avoidance problem of UAVs is ignored. If the trajectories of the UAVs have intersections, the UAV with a low task priority waits for the UAV with a high task priority to pass through the place where a collision may occur. Then, the UAV with the low task priority starts to perform its task again.
    (3) All the UAVs fly at the same constant speed.
    (4) The communication between the UAVs is ignored.
  • Reconnaissance tasks
    (1) The reconnaissance tasks do not require the UAVs to perform the tasks synchronously.
    (2) The UAVs are not required to reach and leave the task points at the same time.
    (3) The reconnaissance time can be ignored.
    (4) If two task points are distributed in the same room, the distance between these two task points can be approximately expressed by the Euclidean distance. If two task points are distributed in different rooms that can be connected directly, the distance between the two task points can be approximately expressed by the sum of d j , h and d h , q . If two task points are distributed in different rooms that can be connected indirectly, the distance between the two task points can be approximately expressed by the sum of d j , h , d h , g , and d g , q . If two task points are distributed in two completely disconnected rooms, the distance between the two task points can be defined as infinity.
  • Indoor environments
    (1) The information on the indoor map is known a priori.
    (2) Two connected rooms share only one door.
    (3) The rooms inside the building are fully connected, that is, a room can reach any other room through a certain topological relationship.
    (4) The no-fly zone and interference zone for the UAVs are not considered in indoor environments.
    (5) The obstacles in indoor environments will not affect the path planning of the UAVs because they can adjust their height to avoid obstacles and the time cost caused by adjusting the height can be ignored.

2.3.2. Decision Variables

Decision variables include two parts: X and P. X = [ x i , j ] N U × N T represents a task assignment scheme. x i , j { 0 , 1 } means reconnaissance task j is assigned to UAV i or not. P = { P 1 ; P 2 ; ; P i ; ; P N U } denotes the task sequence set of each UAV. P i is shown in Formula (14). p i 0 and p i m i + 1 represent the starting point and ending point, that is, the point with coordinates (0,0).
X = x 1 , 1 x 1 , 2 x 1 , N T x i , j x N U , 1 x N U , 2 x N U , N T .
P i = ( p i 0 , p i 1 , , p i m i , p i m i + 1 ) , i I U .

2.3.3. Constraints

In order to ensure the feasibility of the planned scheme, the following constraints need to be met.
• Endurance constraint
l i e i , i I U .
Formula (15) stipulates that the tour length of each UAV must be less than its maximum endurance.
• Ability demand constraint
i I U : u i n x i j θ j n , j I T , n I S .
Formula (16) ensures that for each task, at least one UAV assigned to this task can meet the d-th requirement.

2.3.4. Cost Function

In this paper, the task assignment scheme and the path planning results are evaluated by constructing a specific cost function. The mathematical definition of the ARIRTP problem is shown below:
F ( X , P ) = i = 1 N U l i .
l i = d p i 0 , p i 1 + d p i 1 , p i 2 + + d p i m i 1 , p i m i + d p i m i , p i m i + 1 , m i N T , i I U .
Formula (17) represents the total tour length of all UAVs and the objective of the ARIRTP problem is to minimize the total tour length of the UAVs while performing the tasks. Formula (18) represents the tour length of each UAV.
The ARIRTP problem model is shown below:
m i n F ( X , P ) s . t . ( 11 ) , ( 12 ) , x i , j { 0 , 1 }

3. Algorithm Design

In this section, in order to better solve the ARIRTP problem modeled in the previous section, we propose a bi-level problem-solving framework. The upper level uses an iterative search algorithm to solve the ability-restricted task assignment problem of UAVs. The results of the task assignment of the upper level are input to the lower level. Then, the lower level uses the nearest neighbor algorithm to quickly construct the task sequence of each UAV based on the topology information of buildings. Meanwhile, the upper and lower levels jointly participate in the decoding process to obtain a complete solution to evaluate the task assignment scheme generated at the upper level.

3.1. Bi-Level Problem-Solving Framework

The proposed framework is shown in Figure 2. In the upper level, an iterative search method is used to solve the task assignment problem. Common iterative search algorithms include global search algorithms and local search algorithms or exact algorithms and meta-heuristic algorithms. According to the characteristics of the ARIRTP problem, the encoding and decoding mechanisms are designed and the knowledge contained in the problem is mined. Meanwhile, a solution space compression mechanism suitable for the characteristics of the problem is proposed. After the solution space is compressed, the iterative search algorithms will use the compressed solution space to generate new solutions. All the newly generated solutions can meet the requirements of the tasks. Then, the task list assigned by the upper level for each UAV will be input to the lower level. In the lower level, the nearest neighbor algorithm is used to quickly construct the task-visiting sequences for the UAVs. In addition, the lower level also participates in the decoding process of the upper level and constructs an executable path that conforms to the topological constraints. A complete task-planning scheme is obtained through decoding by both the upper and lower levels. The termination condition will be determined according to the specific iteration search algorithm used in the upper level.

3.1.1. Task Assignment

• Encoding and decoding
Based on the encoding scheme, an individual can be described using Formula (20) and the initial population is formed by randomly generating multiple individuals.
S = ( s 1 , s 2 , , s j , , s N T ) , s j { 1 , 2 , , 2 N u 1 } .
C = c 1 , 1 c 1 , 2 c 1 , N U c k , i c 2 N u 1 , 1 c 2 N u 1 , 2 c 2 N u 1 , N U
s j represents the index number of the combination mode of the UAVs that perform task j. The constant matrix C represents all the combination modes of the UAVs. The vector C ( k , : ) is the k-th row of C and represents the k-th combination mode and c k , i { 0 , 1 } . The data in C ( k , : ) are removed and expressed in binary form, recorded as b k . b 1 =1 and b k + 1 = b k +1, 0 < k < 2 N u 1 , k N .
The decoding process is shown in Table 2. For ease of understanding, an example of three UAVs and ten task points is shown in Figure 3. When the number of UAVs is determined to be 3, the corresponding constant matrix C can be determined. Referring to matrix C, the solution vector S can be converted into a matrix X and then the task list L i of each UAV can be extracted.
• Solution space compression mechanism
The solution space compression mechanism is a method proposed to eliminate UAV combinations that cannot meet the task requirements. When there are more solutions that do not meet the task requirements in the instances, the effect of solution space compression will be more obvious. The algorithm flow of solution space compression is shown in Algorithm 1. First, the set t j = { 1 , 2 , , 2 N u 1 } is constructed, representing all the combination modes of the UAVs. Then, we judge whether the combination mode can satisfy the requirements of each task. According to Formula (16), in one combination mode, at least one UAV assigned to this task can meet the requirements; this combination mode will be retained; otherwise, it will be abandoned.
Algorithm 1 Solution space compression mechanism
1. Input set t j , j I T
2. For j = 1 : N T              \ \ Traverse each task
3.    For k = 1 : 2 N U 1     \ \ Traverse each combination mode of UAVs
4.       If /   u i n θ j n , n I S   \ \ The UAVs in the combination mode k all cannot satisfy the
5.                                            \ \ requirements of task j
6.           t j = t j ∖ { k }
7.       End if
8.    End for
9. End for
10. Output t j , j I T              \ \ The compressed task list
The specific process and results of solution space compression in the example of three UAVs and ten tasks are shown in Figure 4. Since there are two options for whether the UAVs are assigned to perform tasks, the number of the whole assignment scheme is ( 2 3 1 ) 10 . The size of the compressed solution space becomes 0.00007 times the original.

3.1.2. Path Planning

In the lower level, based on the task list of each UAV obtained in the upper level, the path planning is carried out using the nearest neighbor algorithm to obtain the path sequence of each UAV based on an indoor topology map.
• Nearest neighbor algorithm based on topological map
In order to meet the constraints of indoor environments, we use a topology map for path planning. This is necessary to determine not only the sequences of the task points but also the sequence of the path points that can be executed by each UAV. The algorithm flow of the nearest neighbor algorithm based on a topological map is shown in Algorithm 2.
• Path sequence generation method
According to the connectivity of the indoor rooms, it needs to be judged whether any two adjacent task points in the UAV task sequence meet the constraints of the environment. According to the specific indoor environment, the following three situations are summarized: (1) if two connected task points are in the same room, they do not need to be processed; (2) if they are in two directly connected rooms, we need to add the door connecting the two rooms to the sequence, and (3) if they are in two indirectly connected rooms, we will find a connecting sequence of doors based on Formula (12) and add it to the sequence.
Remark 1.
The task list indicates the set of tasks assigned to each UAV after the task assignment. The task sequence represents the visiting sequence vector of the task points of the UAVs. The path sequence indicates the visiting waypoint vector of the UAVs.
Algorithm 2 Nearest neighbor algorithm based on topological map
1. Input: task list L, start point p s t a r t , end point p e n d
2. P i = 0 , P i = 0 , d = 0 , the final path sequence list P i , total tour length l t o t a l = 0
3. For i = 1 : N U
4.     P i = [ P i ; p s t a r t ] , P i = [ P i ; p s t a r t ]
5.    While L i
6.        p s i z e = s i z e ( P i ) , m i = s i z e ( L i )        \ \ Function s i z e ( ) is used to find the length of L i or P i
7.       For j = 1 : m i
8.           d ( j ) = d i s t a n c e ( P i ( p s i z e ) , L i ( j ) )
9.       End for
10.        p c u r r = L i ( a r g m i n ( d ) )        \ \   p c u r r is temporary point
11.        P i = [ P i ; p c u r r ] , P i = [ P i ; p c u r r ] , L i = L i { p c u r r }
12.        m i = m i 1 , p s i z e = p s i z e + 1
13.       If A ( R ( P i ( p s i z e ) ) , R ( P i ( p s i z e 1 ) ) ) ! = 0 & A ( R ( P i ( p s i z e ) ) , R ( P i ( p s i z e 1 ) ) ) ! =
14.                                      \ \ Two points in two directly connected rooms
15.           P i = [ P i ; P D ( ( A ( R ( P i ( p s i z e ) ) , R ( P i ( p s i z e 1 ) ) ) ) , : ) ]
16.       Elseif A ( R ( P i ( p s i z e ) ) , R ( P i ( p s i z e 1 ) ) ) =        \ \ Two points in two indirectly
17.                                                                          \ \ connected rooms
18.           P D i j = a r g m i n ( d i s t a n c e D i j ( D h , D g ) )   \ \ Based on Formulas (11) and (12)
19.                                \ \   P D i j is a coordinate matrix of the doors (connecting the two points)
20.           P i = [ P i ; P D i j ]
21.       End if
22.    End while
23.     P i = [ P i ; p e n d ] , P i = [ P i ; p e n d ]
24.     p s i z e = s i z e ( P i )
25.       If A ( R ( P i ( p s i z e ) ) , R ( P i ( p s i z e 1 ) ) ) ! = 0 & A ( R ( P i ( p s i z e ) ) , R ( P i ( p s i z e 1 ) ) ) ! =
26.           P i = [ P i ; P D ( ( A ( R ( P i ( p s i z e ) ) , R ( P i ( p s i z e 1 ) ) ) ) , : ) ]
27.       Elseif A ( R ( P i ( p s i z e ) ) , R ( P i ( p s i z e 1 ) ) ) =
28.           P D i j = a r g m i n ( d i s t a n c e D i j ( D h , D g ) )
29.           P i = [ P i ; P D i j ]
30.       End if
31.     l i = l e n g t h ( P i ) , l t o t a l = l t o t a l + l i        \ \ Sequence length function l e n g t h ( ) is formula (14)
32.     P i = P i                                         \ \   P i is the visiting waypoint matrix of the i-th UAV
33.     d = 0 , P i = 0 , P i = 0 , l t o t a l = 0
34. End for
35. P = { P 1 ; P 2 ; ; P N U }                          \ \   P is the visiting waypoint matrix of all UAVs
36. Output P , l t o t a l

3.2. Iterative Search Algorithm

The iterative search method is mainly composed of three parts: initializing solutions, evaluating solutions, and generating new solutions according to certain rules. The iterative search algorithm includes the global search algorithm and local search algorithm. As a typical global search algorithm, the GA has been widely used since it was proposed. After years of development, many studies still concentrate on improving the GA at different stages, such as the initialization stage, evolution stage, and individual evaluation, among others [37]. As another promising option, the SA [38] was one of the winning algorithms for the electric vehicle routing problem in the IEEE WCCI2020 competition [39], which shows its efficiency in solving the VRP and its variants or other combinatorial optimization problems [40,41].
Therefore, we apply the GA and SA as iterative search algorithms to the proposed bi-level problem-solving framework, which produces two hybrid algorithms named the GA-NN and SA-NN, respectively.

3.2.1. Genetic Algorithm

• Initialization stage
First, we load the reconnaissance tasks, UAVs, and indoor environment parameters. Then, the population is initialized. For the population initialization process, we use the compressed solution space to randomly generate individuals, which ensures that the UAVs can meet the requirements of each reconnaissance task in the task assignment scheme.
•Evolutionary stage
(1) Fitness evaluation
In order to ensure the algorithm searches in the expected direction, a solution needs to be properly evaluated. In this paper, the fitness function is realized by the following formula:
h p g = 1 ( f p g f m i n g ) / ( f m a x g f m i n g + ε ) .
h p g represents the fitness value of the hth individual in the gth generation and h p g ( 0 , 1 ) . f p g represents the objective function value of the hth individual in the gth generation. f m i n g and f m a x g represent the minimum and maximum objective function values in the gth generation, respectively. We set the parameter ε = 0.001 in order to avoid f m i n g f m a x g + ε =0.
(2) Operator design
Selection operator
In this paper, the objective is to minimize the total tour length of UAVs. When the objective value is smaller, the fitness evaluation value is closer to 1. So, we prefer to retain individuals with a high fitness value. For each individual, we generate a random number. If the fitness evaluation value of the individual is bigger than this random number, the individual is retained. Otherwise, the individual is abandoned, that is, if the following formula is satisfied, we retain the hth individual in the gth generation.
h p g > r a n d , r a n d ( 0 , 1 ) .
Crossover operator
After the selection operation, the selected individuals form a new population, that is, the offspring population. However, the genes of each individual in the offspring population do not change. Now, in order to generate new individuals, it is necessary to change the genes of the individual. In this paper, the uniform crossover operator is used. First, two individuals are randomly selected from the population. Then, each gene of the first individual is traversed and the random number r a n d j ( 0 , 1 ) is generated and compared with the crossover probability p c to determine whether or not to carry out the crossover operation. If r a n d j > p c , the jth gene of the first individual exchanges with the jth gene of the second individual.
Mutation operator
In this paper, the uniform mutation operator is used. First, each gene of the individual is traversed and the random number r a n d j ( 0 , 1 ) is generated and compared with the mutation probability p m . Then, it is determined whether or not to carry out a mutation operation in each gene location. If r a n d j > p m , the jth gene of the individual mutates and the mutation value are randomly selected from t j .
Recombination operator
We set the population size as N P . After the selection operation, some individuals are abandoned. In order to ensure that the population size is N P , the crossover and mutation operations are required to be made to generate new individuals until the population size reaches N P .
To solve the ARIRTP problem, we propose a hybrid algorithm based on the GA and the nearest neighbor algorithm. The algorithm flow is shown in Algorithm 3. The GA-NN method mainly includes three stages: the solution space compression stage, initialization stage, and evolutionary stage. In the initialization stage, the compressed solution space is used to generate new individuals. When evaluating individuals, the GA and the NN algorithm jointly participate in the decoding process. The GA is responsible for decoding the task assignment results of the upper level. The NN algorithm is responsible for decoding the path planning results of the lower level.
Algorithm 3 GA-NN
1. Input: instances
2. Set iteration number g e n = 1 , maximum number of iterations M a x g e n , population size N P
3. Selected population P s e l = , population size P s i z e = 0
4. Solution space compression              \ \ According to Section 3.1.1
5. Initialize current population P           \ \ Task assignment
6. While g e n M a x g e n                          \ \ Main loop
7.    For i = 1 : N P
8.        r = r a n d                                   \ \ Rand number r a n d ( 0 , 1 )
9.        f ( i ) = F ( P ( i ) )   \ \ Objective value of individual i, according to Equations (17) and (18) and Algorithm 2
10.        h ( i ) = f i t n e s s ( f ( i ) )                 \ \ Fitness value of individual i, according to Equation (22)
11.       If h ( i ) > r
12.           P s e l = P s e l { P ( i ) }
13.       End if
14.    End for
15.     P = , P = P s e l , P s e l =
16.     P s i z e = s i z e ( P )
17.    While P s i z e < N P
18.        p 1 = r a n d P s i z e , p 2 = r a n d P s i z e        \ \ Two individuals are randomly selected
19.        ( p 1 , p 2 ) = c r o s s o v e r ( p 1 , p 2 )                 \ \ According to Section 3.2.1
20.        p 1 = m u t a t e ( p 1 ) , p 2 = m u t a t e ( p 2 )        \ \ According to Section 3.2.1
21.        P = P { p 1 , p 2 }
22.        P s i z e = s i z e ( P )
23.    End while
24.    If P s i z e > N P
25.        P = P ( 1 : N P , : )                    \ \ Population size is kept as N P
26.    End if
27.     g e n = g e n + 1
28. End while
29. Output the best individual

3.2.2. Simulated Annealing Algorithm

The SA is different from the GA and other swarm intelligence optimization algorithms that solve problems using the force of the swarm. After performing several neighborhood operations on a current solution, the SA finally obtains the optimal solution that it can search for. The characteristic of the SA in the search process is to accept a solution worse than the current solution with a certain probability. The SA tries to accept a new solution that is slightly worse than the local optimal solution as the current solution and then searches for the current solution. Once the SA finds a better solution than the local optimal solution, it means that the SA has escaped from the local optimal solution.
In this paper, since the value range of each gene in an individual will be different after the solution space is compressed, the traditional neighborhood search strategies, the exchange, reversal, and insertion operations, are not applicable. So, we use the mutation operator designed in Section 3.2.1 to generate new solutions. Meanwhile, the mutation probability should be correspondingly increased.
So, we propose a hybrid algorithm based on the SA algorithm and the nearest neighbor algorithm. The algorithm flow is shown in Algorithm 4.
Algorithm 4 SA-NN
1. Input: instances
2. Set initial temperature T, cooling factor α , stop iteration temperature T s t o p
3. Solution space compression              \ \ According to Section 3.1.1
4. Generate an initialization solution p 0 , set the current solution p c u r r = p 0
5. Set a best solution p b e s t = p 0
6. h ( p b e s t ) = f i t n e s s ( F ( p 0 ) )        \ \ Objective value of individual i, according to Equations (17) and (18) and
7.                                            \ \ Algorithm 2. Fitness evaluation according to Equation (22)
8. Whereas T > T s t o p
9.   For i = 1 : L \ \ Metropolis chain length L
10.     Generate a new solution p n e w from p c u r r based on the neighborhood search rules
11.      If h ( p n e w ) > h ( p c u r r )
12.        p c u r r = p n e w
13.      Else
14.        r = U [0,1] \ \ Generate a random number
15.       If r < e x p [ ( h ( p n e w ) h ( p c u r r ) ) / T ]
16.          p c u r r = p n e w
17.       End if
18.      End if
19.      If h ( p c u r r ) < h ( p b e s t )
20.        p b e s t = p c u r r
21.        h ( p b e s t ) = h ( p c u r r )
22.      End if
23.   End for
24.    T = α T
25. End while
26. Output the best individual

4. Experimental Results

In this section, first, the data sets and parameter settings are described. Then, two groups of experiments are conducted to illustrate the feasibility and superiority of the algorithms. In the first group of experiments, small-scale instances with different degrees of sparsity solution space were used to evaluate the effectiveness of the solution space compression mechanism. Then, in the second group of experiments, the performances of the GA-NN, SA-NN, brute force algorithm, BACO [25], and SJSA [23] methods were compared. In the modeled problems in [23,25], one UAV or vehicle can perform multiple tasks. In the problem modeled in this paper, there is not only a case where a UAV can perform multiple tasks but also a case where multiple UAVs can perform a task cooperatively. All the compared algorithms were implemented in MATLAB R2018b on a workstation (Intel(R) Core (TM) i7-8700 CPU @ 3.20GHz 3.19 GHz, 16.00 GB of RAM) and run independently on each instance 30 times.
In the brute force algorithm, we traversed all the task assignment schemes in the upper level. In the lower level, we used the same nearest neighbor algorithm used in the GA-NN and SA-NN methods to construct the paths.

4.1. Data Sets

The test instances included the locations of the task points, requirements of the reconnaissance tasks, abilities of the UAVs, and information about the buildings. The normal distribution function was used to generate the task requirements and UAV abilities. The locations of the reconnaissance tasks were randomly generated in the selected scene. The number of UAVs was chosen from {3, 4}. The number of tasks was chosen from {10, 15, 20}. The solution space included four different degrees of sparsity. We designed a total of 24 different instances. Here, we specified that if the number of feasible solutions in the instances was less than 0.001 times the size of the original solution space, this type of instance was a sparse type; otherwise, it was a dense type.

4.2. Parameter Settings

In order to fairly compare the running times of the two algorithms, the parameters were adjusted appropriately so that the evaluation times of the individuals were the same when the two algorithms were run once, respectively. We carried out experiments on Instances 1 and 4, respectively (the details of Instances 1 and 4 are shown in Table 3). Through the experimental results, it was found that the GA-NN and SA-NN methods had the best performance when the relevant parameters were set as follows: the population size was set to 200, and 500 generations were given as the training-stopping criterion. The crossover rate was 0.8 and the mutation rate was 0.15. The stop iteration temperature T s t o p = 0.001 . The cooling factor α = 0.99. The initial temperature T = 285 . The metropolis chain length L = 80 .

4.3. Performance Evaluation of Solution Space Compression Mechanism

In the first group of experiments, we compared the GA-NN method with the brute force algorithm to verify the accuracy of the GA-NN method and compared the GA-NN method and the GA-NN method without the SSCM to verify the effectiveness of the SSCM. In the first group of experiments, the number of UAVs was 3, the number of reconnaissance tasks was 10, and the solution space had 4 different densities.
According to theoretical analysis, the performance of the algorithm using the SSCM was better than that of the algorithm without the SSCM. According to the number of feasible solutions, the solution space of the instances was divided into sparse and dense. Therefore, when conducting the experiments on instances with sparse solution spaces, the performance of the algorithm was significantly improved after using the SSCM. Compared to sparse solution spaces, when conducting experiments on instances with dense solution spaces, the effect was not obvious.
In order to test the performance of the solution space compression mechanism, the brute force algorithm, GA-NN method without the SSCM, and GA-NN method were tested on different instances with different solution spaces. The experimental results are shown in Table 3. i j k represents the number of UAVs, reconnaissance tasks, and sparse type, respectively. ‘–’ represents data that cannot be calculated. ‘Concentration’ represents the proportion of feasible solutions to the whole solution space. ‘(0.00002 times)’ represents that the number of feasible solutions is 0.00002 times the size of the original solution space. The minimum values of the objective value were counted because the objective function was to minimize the total tour length of the UAVs.
From the experimental data, we can draw the following conclusions:
(1) When the solution space was sparse, the brute force algorithm quickly found the optimal solution, but when the solution space became dense, the running time of the brute force algorithm became immeasurable.
(2) It was found that the optimal objective value obtained using the GA-NN method was consistent with the optimal objective value obtained using the brute force algorithm, whereas the gap between the optimal objective value obtained using the GA-NN method without the SSCM and the real optimal objective value became larger with the gradual increase in the density of the solution space.
(3) It was found that the running time of the GA-NN method without the SSCM was faster than that of the GA-NN method alone. The reason is that during the operation of the GA-NN method without the SSCM, the newly generated solutions were not always feasible. If the solution was not feasible, the path of the UAVs was not planned so it took less time.

4.4. ARIRTP Problem-Solving Results

In the second group of experiments, we designed different scale instances with four different density solution spaces to compare the performance of the GA-NN and SA-NN methods with the brute force algorithm, BACO, and SJSA. There were a total of 24 instances. The experimental results are shown in Table 4. In order to save space, we selected 12 groups of instances to show their evolution curves. Among them, ‘3-10-1’ represents three UAVs, ten tasks, and the solution space of the first sparse type. The routes of the UAVs in Instance 3-10-1 are shown in Figure 5.
From the experimental data, we can draw the following conclusions:
(1) With the increase in the number of reconnaissance tasks or the gradual density of the solution space, the brute force algorithm failed.
(2) As can be seen from the experimental data, in most instances, the objective function values obtained using the SA-NN and GA-NN methods after the convergence of the algorithms were the same, but the running time of the SA-NN method was longer than that of the GA-NN method.
(3) However, in some large-scale instances, the SA-NN method performed better than the GA-NN method.
(4) The GA-NN method was superior to the SJSA in terms of both the algorithm convergence time and optimal objective function values obtained after the algorithm convergence. After the algorithm convergence, the optimal objective function values obtained using the GA-NN method and the BACO in most instances were the same. However, in some large-scale instances, the objective function value obtained using the GA-NN method was better than that obtained using the BACO. Moreover, the GA-NN method converged faster than the BACO.
(5) The SA-NN method converged faster than the SJSA. After the algorithm convergence, in most instances, the optimal objective function values obtained using the SA-NN method and the SJSA were the same. However, in some large-scale instances, the objective function value obtained using the SA-NN method was better than that obtained using the SJSA.

4.5. Performance Evaluation of the GA-NN and SA-NN Methods

The objective value evolution curves of the GA-NN and SA-NN methods for Instances 3 and 4 are shown in Figure 6 and Figure 7, respectively. However, it can be seen from the figures that the GA-NN method converged faster than the SA-NN method. In this section, we used the Mann–Whitney U test and Kolmogorov–Smirnov test with 95% confidence to test whether there were significant differences between the GA-NN and SA-NN methods. The statistical hypothesis test results are shown in Table 5. As far as the Mann–Whitney U test and Kolmogorov–Smirnov test are concerned, the performance of the GA-NN method was no different from that of the SA-NN method. However, in some large-scale instances, the SA-NN method performed better than the GA-NN method.
Remark 2.
Instances ( / § / ): †, §, and ≈ represent that the performance of the SA-NN method is superior to, inferior to, and no different from the GA-NN method, respectively.

5. Conclusions and Perspectives

This paper focuses on the ability-restricted indoor reconnaissance task-planning problem for multiple UAVs. First, the ARIRTP problem is modeled as a combinatorial optimization problem in which the multi-dimensional requirements of the reconnaissance tasks, multi-dimensional abilities of UAVs, and impact of the connectivity between indoor rooms on the path planning of UAVs are considered. Second, a bi-level problem-solving framework is proposed. The upper level uses an iterative search algorithm to solve the task assignment problem of UAVs. According to the characteristics of the problem, a solution space compression mechanism is proposed to make the generated task assignment schemes meet the requirements of the reconnaissance tasks. The lower level uses the nearest neighbor algorithm to quickly construct the path sequence of UAVs. Third, a hybrid algorithm based on the GA and the nearest neighbor algorithm and a hybrid algorithm based on the SA and the nearest neighbor algorithm are proposed to solve the ARIRTP problem under the bi-level problem-solving framework. The experimental data show that the proposed algorithms, the GA-NN and SA-NN, can solve the ARIRTP problem quickly and accurately. The performance of the GA-NN method is no different from that of the SA-NN method; however, the GA-NN method runs slightly faster. In large-scale instances, the performance of the SA-NN method is slightly better than that of the GA-NN method. Furthermore, the GA-NN and SA-NN methods perform better than the SJSA, and the GA-NN method also performs better than the BACO.
This paper also has some limitations. We assume that the indoor environment has full coverage for communication and that UAVs can communicate with each other. However, due to the complex indoor environment and large number of walls, the communication between UAVs has a certain attenuation.
In future work, we will consider the communication between UAVs in an indoor environment. The connection relationship of the communication topology between UAV teams is helpful for them to cooperate in reconnaissance, rounding up, and other tasks. Initially, we plan to add a maximum communication distance as a constraint to ensure smooth communication between UAV teams. Subsequently, the attenuation of the indoor environment to communication may be considered to make the communication model closer to reality. In indoor environments, when UAVs pass through a door, there are certain risks. In the future, we will consider introducing a risk coefficient to measure the risk degrees of the paths of UAVs.

Author Contributions

Conceptualization, R.Z.; methodology, B.X.; software, R.Z.; validation, R.Z.; formal analysis, R.Z.; investigation, R.Z.; resources, R.Z.; data curation, R.Z.; writing—original draft preparation, R.Z.; writing—review and editing, R.Z.; visualization, R.Z.; supervision, Y.D.; funding acquisition, Q.W. and L.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Outstanding Youth Talent Support Program 61822304, NSFC, under Grant 61873033, the Basic Science Center Programs of the NSFC under Grant 62088101, the Beijing Advanced Innovation Center for Intelligent Robots and Systems, the Shanghai Municipal Science and Technology Major Project (2021SHZDZX0100), the Shanghai Municipal Commission of Science and Technology Project (19511132101), the National Natural Science Fund of China under Grant 62003044, and the National Science Fund for Distinguished Young Scholars of China under Grant 62025301.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ARIRTPAbility-restricted indoor reconnaissance task planning
SSCMSolution space compression mechanism
GAGenetic algorithm
SASimulated annealing algorithm
NNNearest neighbor algorithm
GA-NNGenetic algorithm–nearest neighbor algorithm
SA-NNSimulated annealing algorithm–nearest neighbor algorithm
UAVsUnmanned aerial vehicles
SJSASwap-and-judge simulated annealing algorithm
BACOBi-level ant colony optimization
CEVRPCapacitated electric vehicle routing problem

References

  1. Allam, A.; Nemra, A.; Tadjine, M. Parametric and implicit features-based UAV-UGVs time-varying formation tracking: Dynamic approach. Unmanned Syst. 2022, 10, 109–128. [Google Scholar] [CrossRef]
  2. Zammit, C.; Kampen, E.J. Comparison between A* and RRT algorithms for 3D UAV path planning. Unmanned Syst. 2022, 10, 129–146. [Google Scholar] [CrossRef]
  3. Khan, S.; Tufail, M.; Khan, M.T. A novel framework for multiple ground target detection, recognition and inspection in precision agriculture applications using a UAV. Unmanned Syst. 2022, 10, 45–56. [Google Scholar] [CrossRef]
  4. Li, Q.; Hua, Y.; Dong, X.; Yu, J.; Ren, Z. Time-varying formation tracking control for unmanned aerial vehicles with the leader’s unknown input and obstacle avoidance: Theories and applications. Electronics 2022, 11, 2334. [Google Scholar] [CrossRef]
  5. Zhang, J.; Chen, Y.; Yang, Q.; Lu, Y.; Shi, G.; Wang, S.; Hu, J. Dynamic task allocation of multiple UAVs based on improved A-QCDPSO. Electronics 2022, 11, 1028. [Google Scholar] [CrossRef]
  6. Chung, S.H.; Sah, B.; Lee, J. Optimization for drone and drone-truck combined operations: A review of the state of the art and future directions. Comput. Oper. Res. 2020, 123, 105004. [Google Scholar] [CrossRef]
  7. Macrina, G.; Pugliese, L.D.; Guerriero, F.; Laporte, G. Drone-aided routing: A literature review. Transp. Res. Part Emerg. Technol. 2020, 120, 102762. [Google Scholar] [CrossRef]
  8. Tamke, F.; Buscher, U. A branch-and-cut algorithm for the vehicle routing problem with drones. Transp. Res. Part Methodol. 2021, 144, 174–203. [Google Scholar] [CrossRef]
  9. Liu, Y.; Liu, Z.; Shi, J.M.; Wu, G.; Pedrycz, W. Two-echelon routing problem for parcel delivery by cooperated truck and drone. IEEE Trans. Syst. Man, Cybern. Syst. 2021, 51, 7450–7465. [Google Scholar] [CrossRef]
  10. Vasquez, S.A.; Angulo, G.; Klapp, M.A. An exact solution method for the TSP with drone based on decomposition. Comput. Oper. Res. 2021, 127, 105127. [Google Scholar] [CrossRef]
  11. Luo, Z.; Pan, B.; Zhang, Z.; Liu, Z.; Lim, A. The multi-visit traveling salesman problem with multi-drones. Transp. Res. Part Emerg. Technol. 2021, 128, 103172. [Google Scholar] [CrossRef]
  12. Cui, Y.J.; Dong, W.H.; Hu, D.X.; Liu, H.B. The application of improved harmony search algorithm to multi-UAV task assignment. Electronics 2022, 11, 1171. [Google Scholar] [CrossRef]
  13. Khan, A.; Zhang, J.; Ahmad, S.; Memon, S.; Qureshi, H.A.; Ishfaq, M. Dynamic positioning and energy-efficient path planning for disaster scenarios in 5G-assisted multi-UAV environments. Electronics 2022, 11, 2197. [Google Scholar] [CrossRef]
  14. Rudys, S.; Ragulis, P.; Laučys, A.; Bručas, D.; Pomarnacki, R.; Plonis, D. Investigation of UAV detection by different solid-state marine radars. Electronics 2022, 11, 2502. [Google Scholar] [CrossRef]
  15. Mostafa, S.A.; Mustapha, A.; Gunasekaran, S.S.; Ahmad, M.S.; Mohammed, M.A.; Parwekar, P.; Kadry, S. An agent architecture for autonomous UAV flight control in object classification and recognition missions. Soft Comput. 2021. [Google Scholar] [CrossRef]
  16. Greenwood, W.W.; Lynch, J.P.; Zekkos, D. Applications of UAVs in civil infrastructure. J. Infrastruct. Syst. 2019, 25, 04019002. [Google Scholar] [CrossRef]
  17. Sun, F.; Wang, X.; Zhang, R. Task scheduling system for UAV operations in agricultural plant protection environment. J. Ambient. Intell. Humaniz. Comput. 2020. [Google Scholar] [CrossRef]
  18. Shi, T.; Wang, H.; Cui, W.; Ren, L. Indoor space target searching based on EEG and EOG for UAV. Soft Comput. 2019, 23, 11199–11215. [Google Scholar] [CrossRef]
  19. Sampedro, C.; Rodriguez-Ramos, A.; Bavle, H.; Carrio, A.; de la Puente, P.; Campoy, P. A fully-autonomous aerial robot for search and rescue applications in indoor environments using learning-based techniques. J. Intell. Robot. Syst. 2019, 95, 601–627. [Google Scholar] [CrossRef]
  20. Khosiawan, Y.; Park, Y.; Moon, I.; Nilakantan, J.M.; Nielsen, I. Task scheduling system for UAV operations in indoor environment. Neural Comput. Appl. 2019, 31, 5431–5459. [Google Scholar] [CrossRef]
  21. González de Santos, L.M.; Frías Nores, E.; Martínez Sánchez, J.; González Jorge, H. Indoor path-planning algorithm for UAV-based contact inspection. Sensors 2021, 21, 642. [Google Scholar] [CrossRef] [PubMed]
  22. Bouzid, Y.; Bestaoui, Y.; Siguerdidjane, H. Guidance-control system of a quadrotor for optimal coverage in cluttered environment with a limited onboard energy. J. Intell. Robot. Syst. 2019, 95, 707–730. [Google Scholar] [CrossRef]
  23. Huo, L.; Zhu, J.; Wu, G.; Li, Z. A novel simulated annealing based strategy for balanced UAV task assignment and path planning. Sensors 2020, 20, 4769. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, H.X.; Nan, Y.; Yang, Y. Multi-UAV reconnaissance task assignment for heterogeneous targets based on modified symbiotic organisms search algorithm. Sensors 2019, 19, 734. [Google Scholar] [CrossRef] [Green Version]
  25. Jia, Y.H.; Mei, Y.; Zhang, M.J. A bilevel ant colony optimization algorithm for capacitated electric vehicle routing problem. IEEE Trans. Cybern. 2022, 52, 10855–10868. [Google Scholar] [CrossRef]
  26. Gomez, C.; Fehr, M.; Millane, A.; Hernandez, A.C.; Nieto, J.; Barber, R.; Siegwart, R. Hybrid topological and 3D dense mapping through autonomous exploration for large indoor environments. In Proceedings of the 2020 IEEE International Conference on Robotics and Automation (ICRA), Paris, France, 31 May–31 August 2020; pp. 9673–9679. [Google Scholar]
  27. Ruan, W.Y.; Duan, H.B. Multi-UAV obstacle avoidance control via multi-objective social learning pigeon-inspired optimization. Front. Inf. Technol. Electron. Eng. 2020, 21, 740–748. [Google Scholar] [CrossRef]
  28. Dai, H.; Lu, W.; Li, X.; Yang, J.; Meng, D.; Liu, Y.; Liang, B. Cooperative planning of multi-agent systems based on task-oriented knowledge fusion with graph neural networks. Front. Inf. Technol. Electron. Eng. 2022, 23, 1069–1076. [Google Scholar] [CrossRef]
  29. Dhiman, N.K.; Deodhare, D.; Khemani, D. Where am I? Creating spatial awareness in unmanned ground robots using SLAM: A survey. Sadhana 2015, 40, 1385–1433. [Google Scholar] [CrossRef]
  30. Mccammon, S.; Hollinger, G.A. Topological path planning for autonomous information gathering. Auton. Robot. 2021, 45, 821–842. [Google Scholar] [CrossRef]
  31. He, S.W.; Kilgour, D.M.; Hipel, K.W. A three-level hierarchical graph model for conflict resolution. IEEE Trans. Syst. Man, Cybern. Syst. 2021, 51, 1424–1433. [Google Scholar] [CrossRef]
  32. Zhu, Z.M.; Kilgour, D.M.; Hipel, K.W. A new approach to coalition analysis within the graph model. IEEE Trans. Syst. Man, Cybern. Syst. 2020, 50, 2231–2241. [Google Scholar] [CrossRef]
  33. Zhang, H.; Xin, B.; Dou, L.; Chen, J.; Hirota, K. A review of cooperative path planning of an unmanned aerial vehicle group. Front. Inf. Technol. Electron. Eng. 2020, 21, 1671–1694. [Google Scholar] [CrossRef]
  34. Ding, Y.L.; Xin, B.; Chen, J. A review of recent advances in coordination between unmanned aerial and ground vehicles. Unmanned Syst. 2021, 09, 97–117. [Google Scholar] [CrossRef]
  35. Chen, J.; Ding, Y.; Xin, B.; Yang, Q.; Fang, H. A unifying framework for human-agent collaborative systems-part I: Element and relation analysis. IEEE Trans. Cybern. 2022, 52, 138–151. [Google Scholar] [CrossRef]
  36. Chen, J.; Ding, Y.; Xin, B.; Yang, Q.; Fang, H. A unifying framework for human-agent collaborative systems-part II: Design procedure and application. IEEE Trans. Cybern. 2021, 52, 11990–12002. [Google Scholar]
  37. Tao, Y.; Wen, Y.; Gao, H.; Wang, T.; Wan, J.; Lan, J. A path-planning method for wall surface inspection robot based on improved genetic algorithm. Electronics 2022, 11, 1192. [Google Scholar] [CrossRef]
  38. Xiao, S.; Tan, X.; Wang, J. A simulated annealing algorithm and grid map-based UAV coverage path planning method for 3D reconstruction. Electronics 2021, 10, 853. [Google Scholar] [CrossRef]
  39. Available online: https://mavrovouniotis.github.io/EVRPcompetition2020/ (accessed on 19 July 2020).
  40. Wang, Y.Y.; Jiao, X.H. Multi-objective energy management for PHEV using Pontryagin’s minimum principle and particle swarm optimization online. Sci.-China-Inf. Sci. 2021, 64, 119204. [Google Scholar] [CrossRef]
  41. Chen, L.; Xin, B.; Chen, J. Interactive multiobjective evolutionary algorithm based on decomposition and compression. Sci.-China-Inf. Sci. 2021, 64, 202201. [Google Scholar] [CrossRef]
Figure 1. An example of three UAVs performing seven reconnaissance tasks in an indoor environment.
Figure 1. An example of three UAVs performing seven reconnaissance tasks in an indoor environment.
Electronics 11 04227 g001
Figure 2. The bi-level problem-solving framework.
Figure 2. The bi-level problem-solving framework.
Electronics 11 04227 g002
Figure 3. The decoding process of the upper level of the example of three UAVs and ten tasks.
Figure 3. The decoding process of the upper level of the example of three UAVs and ten tasks.
Electronics 11 04227 g003
Figure 4. Example of solution space compression process of three UAVs and ten tasks.
Figure 4. Example of solution space compression process of three UAVs and ten tasks.
Electronics 11 04227 g004
Figure 5. The routes of the UAVs in Instance 3-10-1.
Figure 5. The routes of the UAVs in Instance 3-10-1.
Electronics 11 04227 g005
Figure 6. Objective value evolution curve for Instance 3.
Figure 6. Objective value evolution curve for Instance 3.
Electronics 11 04227 g006
Figure 7. Objective value evolution curve for Instance 4.
Figure 7. Objective value evolution curve for Instance 4.
Electronics 11 04227 g007
Table 1. List of parameters and variables.
Table 1. List of parameters and variables.
SymbolsDescription
N U Total number of UAVs
N T Total number of tasks
N S Total number of sensor types
iIndex of UAVs
j , q Index of tasks
h , g Index of doors
nIndex of sensor types
kIndex of the combination mode of UAVs
l i Tour length of UAV i
e i Maximum endurance of UAV i
θ j n Threshold of the demand on the n-th
sensor type in task j
u i n Ability value of the n-th sensor type for UAV i
L i Task list of UAV i
P i Task sequence of UAV i
m i Total number of tasks assigned to UAV i
d j , q Distance between the position of task j and task q
d j , h Distance between the position of task j and door h
d h , g Distance between the position of door h and door g
d g , q Distance between the position of door g and task q
X = [ x i , j ] N U × N T Assignment scheme
x i , j { 0 , 1 } Task j is assigned to UAV i or not
L = { L 1 ; L 2 ; ; L N U }Task list of total UAVs
P = { P 1 ; P 2 ; ; P N U }Task sequence set of total UAVs
I U ={1, 2, ⋯, N U }Index set of the UAV
I T ={1, 2, ⋯, N T }Index set of the task
I S ={1, 2, ⋯, N S }Index set of the sensor type
Table 2. The decoding process of the upper level.
Table 2. The decoding process of the upper level.
Decoding of the upper level
1. Vector S
2. Construct combination mode matrix C        \ \ According to Formula (17)
3. Transfer S to X referring to C  \ \ Obtain a complete task assignment scheme
4. Obtain L i based on X                    \ \ Obtain the task lists
5.    Set L i =
6.    For i = 1 : N U
7.       For j = 1 : N T
8.           x i , j = C ( s j , i )
9.          If X ( i , j ) = = 1
10.              L i = L i { j }
11.          End if
12.       End for
13.    End for
14. Output L
Table 3. Verification experiment of solution space compression mechanism.
Table 3. Verification experiment of solution space compression mechanism.
InstanceScaleConcentration Brute Force
Algorithm
GA-NN without
SSCM
GA-NN
i j k tmintmintmin
13-10-1Sparse
(0.00002 times)
1.439 s 2444 0.688   s 27620.937 s 2444
23-10-2Sparse
(0.00014 times)
43.892 s 2036 0.706   s 25851.022 s 2036
33-10-3Dense
(0.003 times)
451.579 s 1928 0.716   s 27761.185 s 1928
43-10-4Dense
(0.019 times)
0.709   s 27621.209 s 1478
Table 4. Experimental results.
Table 4. Experimental results.
InstanceScaleConcentration Brute Force
Algorithm
SA-NN GA-NN BACO SJSA
i j k tmintmintmintmintmin
13-10-1Sparse
(0.00002 times)
1.439 s 2444 8.178 s 2444 0.937 s 2444 0.771   s 2444 9.101 s 2444
23-10-2Sparse
(0.00014 times)
43.892 s 2036 7.269 s 2038 1.022 s 2036 0.836   s 2038 7.026 s 2036
33-10-3Dense
(0.003 times)
451.579 s 1928 7.676 s 1928 1.185 s 1928 1.089   s 1928 8.755 s 1928
43-10-4Dense
(0.019 times)
7.624 s 1478 1.209 s 1478 1.038   s 1478 8.034 s 1478
53-15-1Sparse
(0.00000011 times)
321.018 s 3620 11.582 s 3620 0.945   s 3620 1.912 s 3620 13.135 s 3620
63-15-2Sparse
(0.00038 times)
12.310 s 2878 1.144   s 2878 1.384 s 2878 14.347 s 2878
73-15-3Dense
(0.0031 times)
11.724 s2144 1.135   s 2128 1.209 s 2128 13.078 s 2128
83-15-4Dense
(0.0111 times)
8.901 s 252.05 1.874   s 252.05 1.921 s 252.05 11.668 s 252.05
93-20-1Sparse
(0.0000000006 times)
17.188 s 4206 1.767   s 4206 1.861 s 4206 20.181 s 4206
103-20-2Sparse
(0.000002 times)
15.059 s 3434 1.957   s 3434 2.762 s 3434 17.796 s 3434
113-20-3Sparse
(0.00036 times)
16.316 s4028 2.624   s 2960 3.118 s 2960 17.817 s 2960
123-20-4Dense
(0.0028 times)
10.665 s 308.8 2.737   s 308.08 3.564 s 308.8 15.101 s 308.8
134-10-1Sparse
(0.00002 times)
9.791 s 2444 1.245   s 2444 3.570 s 2444 13.268 s 2444
144-10-2Sparse
(0.00089 times)
8.426 s 2037 1.536   s 2037 2.854 s 2037 14.786 s 2037
154-10-3Dense
(0.0014 times)
8.942 s 2170 1.785   s 2170 3.859 s 2170 14.135 s 2170
164-10-4Dense
(0.0601 times)
8.908 s 1477 1.733   s 1477 2.616 s 1477 12.078 s 1477
174-15-1Sparse
(0.00000012 times)
13.933 s 3620 2.579   s 3620 3.277 s 3620 15.347 s3639
184-15-2Sparse
(0.0000022 times)
15.607 s 3576 3.320   s 3576 3.645 s 3576 17.455 s3693
194-15-3Sparse
(0.000048 times)
13.645 s 2730 3.428   s 2730 3.716 s 2730 16.078 s2983
204-15-4Dense
(0.1501 times)
9.734 s 1675 2.721   s 1675 3.233 s 1675 13.263 s1836
214-20-1Sparse
(0.0000000007 times)
15.777 s 3611 3.495   s 42065.196 s431019.181 s4341
224-20-2Sparse
(0.0000014 times)
14.654 s 3710 5.080   s 37156.401 s 3710 19.822 s3857
234-20-3Sparse
(0.0000338 times)
12.500 s 2826 4.025   s 40286.005 s415219.796 s3169
244-20-4Dense
(0.087 times)
12.305 s 1895 5.868   s 20667.163 s220815.817 s1930
Table 5. Statistical hypothesis test results.
Table 5. Statistical hypothesis test results.
123456789101112
Mann–Whitney U test0/0/300/0/300/0/300/0/300/0/300/0/300/0/300/0/300/0/300/0/300/28/20/0/30
Kolmogorov–Smirnov test0/0/300/0/300/0/300/0/300/0/300/0/300/0/300/0/300/0/300/0/300/27/30/0/30
Instances131415161718192021222324
Mann–Whitney U test0/0/300/0/300/0/300/0/300/0/300/0/300/0/300/0/3028/0/227/0/330/0/024/0/6
Kolmogorov–Smirnov test0/0/300/0/300/0/300/0/300/0/300/0/300/0/300/0/3030/0/029/0/130/0/022/0/8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, R.; Dou, L.; Wang, Q.; Xin, B.; Ding, Y. Ability-Restricted Indoor Reconnaissance Task Planning for Multiple UAVs. Electronics 2022, 11, 4227. https://doi.org/10.3390/electronics11244227

AMA Style

Zhang R, Dou L, Wang Q, Xin B, Ding Y. Ability-Restricted Indoor Reconnaissance Task Planning for Multiple UAVs. Electronics. 2022; 11(24):4227. https://doi.org/10.3390/electronics11244227

Chicago/Turabian Style

Zhang, Ruowei, Lihua Dou, Qing Wang, Bin Xin, and Yulong Ding. 2022. "Ability-Restricted Indoor Reconnaissance Task Planning for Multiple UAVs" Electronics 11, no. 24: 4227. https://doi.org/10.3390/electronics11244227

APA Style

Zhang, R., Dou, L., Wang, Q., Xin, B., & Ding, Y. (2022). Ability-Restricted Indoor Reconnaissance Task Planning for Multiple UAVs. Electronics, 11(24), 4227. https://doi.org/10.3390/electronics11244227

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop