Next Article in Journal
Achromatic and Athermal Design of Aerial Catadioptric Optical Systems by Efficient Optimization of Materials
Next Article in Special Issue
Movie Scene Event Extraction with Graph Attention Network Based on Argument Correlation Information
Previous Article in Journal
An Innovative Cloud-Fog-Based Smart Grid Scheme for Efficient Resource Utilization
Previous Article in Special Issue
Indoor Positioning Design for Mobile Phones via Integrating a Single Microphone Sensor and an H2 Estimator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Swarm Intelligence Algorithm with a Parasitism-Relation-Based Structure for Mobile Robot Path Planning

1
School of Information and Communication Engineering, Communication University of China, No.1 Dingfuzhuang East Street, Chaoyang District, Beijing 100024, China
2
State Key Laboratory of Media Convergence of Communication, Communication University of China, Beijing 100024, China
3
Key Laboratory of Acoustic Visual Technology and Intelligent Control System, Ministry of Culture and Tourism, Beijing 100024, China
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(4), 1751; https://doi.org/10.3390/s23041751
Submission received: 29 December 2022 / Revised: 19 January 2023 / Accepted: 2 February 2023 / Published: 4 February 2023
(This article belongs to the Special Issue Smart Mobile and Sensing Applications)

Abstract

:
A multi-swarm-evolutionary structure based on the parasitic relationship in the biosphere is proposed in this paper and, according to the conception, the Para-PSO-ABC algorithm (ParaPA), combined with merits of the modified particle swarm optimization (MPSO) and artificial bee colony algorithm (ABC), is conducted with the multimodal routing strategy to enhance the safety and the cost issue for the mobile robot path planning problem. The evolution is divided into three stages, where the first is the independent evolutionary stage, with the same evolution strategies for each swarm. The second is the fusion stage, in which individuals are evolved hierarchically in the parasitism structure. Finally, in the interaction stage, a multi-swarm-elite strategy is used to filter the information through a predefined cross function among swarms. Meanwhile, the segment obstacle-avoiding strategy is proposed to accelerate the searching speed with two fitness functions. The best path is selected according to the performance on the safety and consumption issues. The introduced algorithm is examined with different obstacle allocations and simulated in the real routing environment compared with some typical algorithms. The results verify the productiveness of the parasitism-relation-based structure and the stage-based evolution strategy in path planning.

1. Introduction

With the development of control technology, mobile robotics has been developed in multiple fields, such as rescue, military, industry, etc. As one of the essential parts of the mobile robotics field, the path planning problem establishes an effective path for the robot to reach the target and finish the task without any collisions based on the specific environment. According to the preknowledge of the working environment, the path planning approach can be classified into two parts, which are the global path planning with static information and the local planning based on the sensor information. The former one needs to find a suitable route according to previous map information. While the local aspect should have the decision capacity upon real-time information and find where the problem is, such as the local obstacle distribution, so as to optimize a solution from the current node to a sub-target until the mission is completed. Successful planning should meet both optimization criteria in terms of time and traveling distance, etc.
Recently, various methods of intelligent planning have been studied to find the most productive solution. All the methods can be classified into traditional path-planning methods based on environment modeling, search-based method, and artificial intelligence algorithm (seen in Figure 1). The traditional path planning algorithms are built with the previously defined map, and most of them require environment information in advance to guide robot movement [1] or make a mobility prediction [2]. It is normally utilized to solve a global problem, such as the potential field category, which produces an artificial field based on the motion environment. The movement of the robot can be guided by descent direction, such as gravity, to avoid the repulsive fields (obstacle) from the start to the target. However, it cannot guarantee the path is the global best, even for a certain search range. Hence, the optimization algorithm is normally utilized to optimize the path generated by the artificial potential field (APF), so as to increase the effectiveness of the hybrid algorithm [1]. Another category is the search-based method, such as Dijkstra and A* algorithm, whose complexities increase with the dimensions of problems, resulting in lower effectiveness [3]. Further, a dynamic environmental problem is hard to deal with for those methods that may increase the cost of the planning. Moreover, the environment-modeling-based method, such as the Voronoi Diagram [4], Visibility Graph [5], and Cell Decomposition [6], decomposes the environment into several regions and transforms the complex workspace into a simple map search problem. This type of algorithm has a strong ability to guarantee safety, but the local optima cannot be avoided.
The artificial intelligence algorithm is probably the most commonly adopted approach for mobile robot path planning, which transforms the best path planning into a constrained optimization with one or several objective functions. The optimization algorithm is utilized to solve the nonlinear and multi-constrained models for an optimum or approximate solution, such as particle swarm optimization (PSO), ant colony optimization (ACO), artificial bee colony algorithm (ABC), genetic algorithm (GA), artificial neural network [7,8,9,10,11,12], etc. They have a great ability to handle the uncertainty condition in the complex environment and are flexible for global or local routing problems. However, the performance of this kind of method is unstable due to defects in the algorithm, such as the premature problem and the balance between exploration and exploitation, which may cause inferior work efficiency with a redundant computation. Therefore, the algorithm should be modified according to its merits. For instance, to address the slow convergence problem in the GA, Hong et al. introduced a co-evaluation strategy to provide some margin of error during evolution. Similarly, Ref. [13] proposed a parallel strategy in ACO to address the premature problem. The hybrid algorithm is another productive and flexible method that can merge the advantages of each algorithm to find the best path. Ref. [14] introduced a hybrid algorithm that combined biogeography-based optimization (BBO) and PSO based on the Voronoi diagram to deal with static path planning. Normally, a traditional path planning algorithm could generate a more stable static path. Ref. [15] utilized the merit of the potential field method, combined with a bacterial evolutionary algorithm, to reduce the disadvantages of the intelligence algorithm when it is applied in an environment with a dynamic condition. The traditional map-based method is applied with the evolution algorithm is also an effective method. Ref. [16] used fuzzy logic to enhance the searchability for a dynamic path optimization based on the A* and GA. In addition, the hybrid method can extract some special evolution mechanisms to compensate for the defect in the iteration process. Ref. [17] pointed out that the parameters in the bat algorithm (BA) can be optimized by PSO in multi-objective optimization. In this category, however, the robot always moves precisely along the predetermined path [18]. In the case of an accident, the original path would be affected and redesigned with a local search method. Hence, the falling into a local optimum and computational complexities problems of the intelligence algorithm is trickier, so as to reduce the robustness and working safety. The convergence recourse should be regrouped and allocated to balance exploration and exploitation. Most swarm intelligence algorithms are improved by the evolutionary approach, such as the whale algorithm [19] and the bat algorithm [20], which changes the form of evolution. However, such modification cannot solve the search balance problem. The other method is to construct evolutionary structures, such as ABC [21], artificial fish algorithm [22], or wolf colony algorithm [23]. Although these algorithms can effectively improve the efficiency of population utilization by allocating responsibilities to the populations, they still have inevitable problems, such as the weak local exploitation ability in the ABC. The reasons are that it cannot control the evolutionary direction and the nectar replacement mechanism, resulting in a lower convergence accuracy.
In this paper, a conception of the dual-community-evolutionary structure inspired by the parasitism relation is proposed to balance the exploration and exploitation in optimization. Then, the parasitism-relation-based algorithm, composed of the PSO and ABC algorithm (ParaPA), is applied tp the path planning problem. In the paradigm algorithm, the memory swarm in PSO is utilized to prevent the loss of optima at the superior level, while the evolution approach in ABC is tailored to maximize the convergence-resource usage at the bottom level. The swarm intelligence algorithm is a search optimization based on probability where a more uniform distribution at an earlier stage obtains a better optimal solution easier. Hence, during the evolution, the initialization phase uses a chaos-based logistic map to create a chaotic status in the first stage. Then, the personal best particles in PSO are selected in the superior population, whereas the global best is produced by the nectar selection method in the ABC process. Meanwhile, the environment is divided into several segments to decrease the dimension of the variable, which improves the adaptability of the algorithm. In addition, multi-swarm evolution with the elitist-based information changing strategy is conducted to guarantee the algorithm diversity directly in the upper layer through the cross function. Finally, the proposed algorithm is examined in some path-planning environments and compared with other path-planning for verifying its effectiveness and safety.
Notations: 
· represents Euclidean norm. R n is the n-dimensional real number space. ⊊ is the non-true subset operator.

2. Preliminary Knowledge and Analysis of ABC Algorithm

2.1. Basic Conception of Artificial Bee Colony Algorithm

The artificial bee colony algorithm is derived from the honey-harvesting behavior of honey bees. The bee colony is divided into several groups with different tasks and shares the group information to find the optimal solution. The major assignments can be divided into two swarms, namely employed bees and onlookers [21]. The former forages for food sources and searches only in the local region. The number of food sources equals the onlooker bees population. Once the nectar collection is completed, associated bees can become scout bees and repeat the food search within the entire space and bring the new location into the colony. The procedure for mass fundamentals is described below:
A.
Employed bee
In charge of the food exploitation by the defined equation as follows:
v i j = x i j + ϕ x i j x k j
where i k , j indicates the dimension index from 1 , 2 , D , and ϕ is a random factor selected from 1 , 1 . X i = x i 1 , x i 2 , , x i D represents the current food source location, D is the dimensions number of X i , and v i = v i 1 , v i 2 , , v i D is the location of a new food source searched for by employed bee. Note that if the fitness of the new sources is higher, the memories of employed bees will keep the new position. Otherwise, the previous one is kept to the next iteration until the nectar is substituted.
B.
On looker
A bee is looking for suitable food sources by roulette wheel selection, as shown in Equation (2):
p i = f i t i i = 1 S N f i t i
where f i t i is the fitness value of X i , and p i the is the probability of selection. Once the target is locked, onlookers are transformed into employed bees and exploit food sources.
C.
Scout bee
If a food source cannot be further improved, new nectar would be produced randomly in the searching space, using Equation (3):
x i j = x min , j + r a n d o m 0 , 1 x max , j x min , j
where x max , j and x min , j are the maximum and minimum search bounds.

2.2. The Pros and Cons Analysis of ABC

In ABC, the group evolution strategy can maximize the utilization of convergence resources. For instance, the setting of the nectar harvesting mechanism allows the algorithm to retain a dynamic search ability even in the later period. However, because of this kind of mechanism, a potential optimum may be disturbed, resulting in the algorithm being unable to achieve better local convergence accuracy in the limited number of iterations. To address the question, Li et al. pointed out that the successful experience generated in the evolution can be used to guide the next foraging behavior [24]. The best value is kept throughout the whole process. While Ref. [25] found that the nectar collection is a random process, if different nectars have the same objective value, the search phase of onlookers might be invalid in the evolution. Therefore, they proposed a neighborhood selection method to improve the updated format. In addition, Zhou et al. introduced a multi-elite strategy to increase the guidance in position updating. Moreover, the employed bee and onlooker bee adapt two different search equations [26]. Nevertheless, this approach does not solve the problem of search range overlapping among elites. Regardless of either approach, the fundament is to prevent the loss of optimal solutions with iterations at the cost of the exploration capacity of the algorithm.
The diversity and local convergence accuracy in the later stage is a pair of contradictory factors. A better solution is to explore the search space as much as possible in the initialization and establish an appropriate promotion method in the middle stages of evolution so as to avoid the limitations on diversity caused by the homogeneous potential optimum in the later stages. While the ABC algorithm can deal with the problem by collecting and reusing the converging resources through the division of functional responsibilities into different swarms. These advantages can be utilized, and this article conducts a novel approach incorporated with PSO to address the shortcomings of the ABC algorithm in a dual-community structure. The community can be regarded as a swarm with a specific evolution pattern.

3. The Proposed Parasitism-Relation Structure and ParaPA Algorithm

It is worth noting that the improvement on certain parts of the algorithm always comes at the expense of other performance, and this approach can only achieve better application results with specific requirements. However, in most of the modifications, the efficiency of the algorithm does not increase, and the extra parameter setting in the mechanism also asks for higher requirements during the initialization process. While the structure of dual communities can compensate for the defects of each algorithm by establishing the parasitic relationship to integrate the advantages of algorithms. The problem is how to ensure the information interchange and rebuild the evolutionary mechanism based on the relationship to balance the exploration and exploitation ability during the whole process. Hence, the major purpose of the structure is to address the search balance problem.

3.1. The Parasitism-Relation Structure

The superiority of the PSO is that it can utilize the best personal and global memories to increase evolution efficiency. This paper establishes the relationship between particle swarms in PSO and bee colonies by the parasitism phenomenon in the biosphere. In previous studies, some researchers have used the symbiosis phenomenon to construct the evolutionary process (symbiotic organisms search algorithm, called the SOS algorithm) [27,28]. Although the same conception of symbiosis is utilized, the understanding of symbiosis is quite different. In essence, the SOS algorithm divides the evolution process into three stages, named mutualism, commensalism, and parasitism, which correspond to the initialization stages for producing diversity, the middle stage for evolving toward the best position, and the last stage for preventing populations from stagnating by generating a random sample, respectively. In terms of modes and the formulation established in SOS, the essence is a DE algorithm with a phased evolution. As the concept of symbiosis described in Ref. [29], to mimic the symbiosis phenomenon, a more appropriate approach is to focus on the swarm relationship rather than individuals. The survival of any swarm is determined by whether it can converge to the final evolution. While in a parasitic relationship, the host has a superior resource, and the inferior party can only attach to the dominant position to keep evolution and competition for a chance to survive. The organism, however, has limited nutrition, which means the exploited position should keep changing throughout the whole process. When the convergence resources are exhausted, the attached particle also loses the right to survive and thus enters a chaotic state. The interpretation of this relation in mathematics is (Note that the multi-path planning is related to the multimodal optimization, which is adapted in this paper.):
Definition 1.
Global best. If x * S , for x S where S is the search space with S R n , can have f x f x * , then the x * is regarded as the global best decision variable in S and f x * is named the global best value.
Definition 2.
Multimodal Function. f * is the best value on S based on f, if there are different decision variables x 1 , x 2 , , x m S where f x i ( i = 1 , 2 , m ) are the global best or local best value, the function f x can be named as the multi-peaks function or multimodal function.
Based on Definition 2, the process of finding the best decision variables in S based on the multimodal function f can be named multimodal optimization.
Definition 3.
Parasitism relationship in multimodal optimization. S is regarded as the living space for the decision variable x. Suppose there are m optima in S, recorded as Φ 1 , Φ 2 , , Φ m , S 1 S , when the iteration is t and the S 1 can be recorded as S 1 t , then P , 0 < P < 1 , satisfy
lim t i = 1 m P Φ i S 1 t P
P is the convergence probability, which is defined in fuzzy logic. Further, for S 2 S , S 1 S 2 , for 0 < δ < 1 , can satisfy
lim t i = 1 m P Φ i S 2 t < δ
Then called the S 1 and S 2 is the parasitism relationship in the living space S.
Definition 4.
Suppose Φ i is the i-th optimum in the living space S, if β S , δ > 0 , if Φ i β < δ , then
P Φ i S = 1 .
If a threshold value is determined, let 0 < θ 1 < 1 , θ 2 and 0 < θ 2 < θ 1 , when β S and Φ i β < θ 1 , then
P Φ i S = Φ i β θ 1 θ 2 .
Further, for β S , it always has Φ i β > θ 1 , then
P Φ i S = 0
To summarize, as
P Φ i S = 1 , T o t a l l y c o n v e r g e n c e 0 , T o t a l l y n o t c o n v e r g e n c e Φ i β θ 1 θ 2 , F u z z y c o n v e r g e n c e
The commensalism and mutualism relationship is described in Appendix A.
In the ParaPA algorithm, the bee colony is regarded as the host part to control the better convergence resources, while the survival of particles in PSO is living in the bottom population, such as the structure, as shown in Figure 2. The life of particles in PSO swarm leeches onto the onlookers and personal best swarms. In other words, the bottom particles cannot find the evolution direction without the instruction from the pbest and bees swarms. If the populations are divided into different hierarchies according to their objective function fitness, the result can be illustrated as in Figure 3. The S 4 level represents the bottom swarm whose evolution is dependent on S 3 . Meanwhile, the S 3 level is constructed from S 4 , and they hold exclusive evolution strategies, such as the multi-swarm-elites strategy in pbest swarm and look limitation strategy in the bee colony, which also aims to enhance the algorithm performance in exploitation and diversity, respectively. S 1 and S 2 belong to the superior level, which is mainly responsible for exchanging the best information in its population with other populations and also conveying the feedback to S 3 . Hence, the major function of a superior level is to guarantee the information changes among the different populations to avoid evolutionary problems from a monotonous population.

3.2. Evolutionary Process of ParaPA

3.2.1. Independent Evolution Stage

A chaotic distribution has been proven to achieve a better statistical property and have faster convergence to the algorithm [30]. To prosper the diversity of the food sources initialization, in this paper, the Logistic map, which is widely used in the chaos-based initialization, is given as follows:
x n + 1 = μ x n 1 x n
where n is the number of chaotic variables and x 0 0 , 0.25 , 0.5 , 0.75 , 1 . The chaotic control parameter μ 3.57 , 4 . We use μ = 4 to produce the chaotic system.
During the initiation, the food resources are selected based on objective fitness, and then all the particles in PSO will transform to onlooker bees to choose the nectar as the global best. Moreover, personal memory is also kept to create more possibilities during the iteration.

3.2.2. The Fusion Stage

Based on the independent evolution, the novel position update incorporated with the bee colony is introduced as follows:
V i d t + 1 = ω V i d t + r 1 p b e s t i d t X i d t + r 2 o n l o o k e r i d t X i d t
X i d t + 1 = X i d t + V i d t + 1
where V i d is the d dimensional velocity of the current particle after t iterations, X i d is the d dimensional location of the current particle after t iterations, and the inertial weight w is utilized for the exploration ability, changing from 0.1 , 0.9 (as shown in Equation (7)) through the evolutionary process, which is a self-regulating method.
ω = ω max ω max ω min g M a x g e n
where g is the current iterative time, and Maxgen is the total number of iterations. r 1 and r 2 are two constants, which are taken randomly from 0 to 2 in this paper. The global best is replaced by the best onlooker bee, which is selected by the roulette wheel, as shown in Equation (2).
The limitation of nectar collection is set to three times the number of the population, which means that when all individuals visit the current solution more than three times, the nectar should be abandoned and the onlooker bee on the current nectar would be respawned inside the solution space according to Equation (3). Meanwhile, the original onlooker bees, regarded as the host, compete with the parasitic individuals during the evolution. The host would be directly replaced once a better position appears and the visited number is also restarted. It is worth noting that there is no memory preservation for the onlooker. Originally, the attached individuals would fall into a chaos status when the onlooker is regenerated. However, due to the record of the personal best position, the parasitic individual can still move toward the best position in its memory after losing the global guidance so that the algorithm can keep the exploitation capacity in the middle stage of evolution and avoid the waste of iterations.

3.2.3. Interaction Stage with the Multi-Swarm Elite Strategy

In order to enhance the diversity of the algorithm, this paper adopts a multiple swarm parallel strategy. Each swarm is an independent biosphere, and its parasitic relationships are bound in the biosphere, which is not influenced by other environments. To avoid overlapping regions of the host during the convergence, each population takes out the best memory to join the mixing pool after each evolution (as illustrated in Figure 2). The quality of all elites in the mixing pool is judged by the cross function instead of the previous objective function. Furthermore, the worst individuals are usually discarded via evolution, but the amount of information carried by them can create substantial value in multimodal situations. For the consideration of diversity, the worst individuals in the personal best memory are also mixed with the individuals selected from the mixed pool who have better cross function values and then re-screened by the related objective function. Finally, the individuals with the worst fitness in each swarm would be replaced by those chosen memories in the next iteration (the pseudocode of multi-swarm elites selection is shown in Algorithm 1).
Algorithm 1 Multi-swarm-elite selection strategy.
Input: Select elites from the best personal memories in each sub-swarm based on the objective functions and add them to the mixed pool.
Steps: for each sub-swarm
              1. Select several worst personal memories in each sub-swarm (suppose the number is w).
              2. Choose the best individuals from the mixed pool according to the cross fitness. Then, mix the worst w individuals into with them and remove other elites in the pool.
              3. for each individual in the mixed pool
                         Evaluated it by the cross function
                              if current individual is better
                                      Record it in the pool-output group (In this paper, the population of pool-output group is equal to 1)
                              end
                   end
              4. Replace the worst individuals in sub-swarm by the pool-output group into next iterations
             end
Output: Each sub-swarm with elites who has the best cross fitness. To note that the replaced individuals always keep the best memories.

4. Problem Formulation and the Strategies for Path Planning

In this paper, several static obstacles are listed in the environment along with some premises and assumptions, shown as follows:
(1)
Global path planning is the main target in this paper, which means all obstacles are known before algorithm execution.
(2)
The environment of path planning is built in a 2D workspace. The path planning will consider obstacle avoidance without height.
(3)
If the planning path can avoid the obstacles successfully under environmental constraints, it means the algorithm has the ability to build a safe path in the static condition. Hence, the physical characteristics of the robot are not considered in this paper.
(4)
In the static condition, the robot speed is constant.

4.1. Workspace Formulation

Normally, the path planning question has two or three decision variables. Although it can enrich the diversity in the evolution process, the convergence efficiency of the algorithm will be significantly reduced in the face of complex obstacle situations. Therefore, a new coordinate is applied to connect the start and target positions, as shown in Figure 4. Then the new x -axis is divided into N segments averagely, which means the x dimension is fixed in the new coordinate. While the y dimension is optimized along the parallel line L i ( i = 1 , 2 , , N 1 ), which is vertical to the new x -axis. Hence, the planning path can be presented as ( S 1 , S 2 , S 3 , , S n ), where n = 1 , 2 , , N and the start and target point are fixed. The corresponding transformation formula is shown in Equation (8).
x y = cos ϕ sin ϕ sin ϕ cos ϕ x y + S t a r t x S t a r t y
where the x , y and S t a r t x , S t a r t y belong to the original cartesian coordinate frame. ϕ is the rotation angle to the new frame. Moreover, the map boundary is changed. For instance, in Figure 4, the new y -axis intersects with the original y-axis. For any point on the parallel line L i , the boundary of y in the new frame is determined as follows:
y max = x i sin ϕ , x i < x m a n d y i < y m u l y y i cos ϕ , x i > x m a n d y i > y m
y min = y i cos ϕ , x i < x m a n d y i < y m u l x x i sin ϕ , x i > x m a n d y i > y m
where u l y and u l x are the upper limits of the y-axis and x-axis, respectively. Further, the segment points are located on the x -axis, such as P i ( x i , y i ) , with respect to the original frame. For any points that go over the boundary during the evolution, the whole path ( S 1 , S 2 , S 3 , , S n ) should be regenerated in the solution space. To satisfy the safety issue, each segment S n should not intersect with any obstacles, and all the objective functions are calculated on the original cartesian coordinate frame.

4.2. Design of Objective Function

The quality of a path is examined by the fitness value. In the previous study, only feasible paths were evaluated, and information on infeasible paths was neglected [31], which makes it difficult to take the valid yet easily ignored information into account for global planning. Therefore, in this paper, different objective functions are used for the assessment of feasible and infeasible paths, respectively. According to the previous definition, when any segments are overlapped with obstacles, the current path is regarded as an infeasible solution. For the evaluation of infeasible paths, the penalty for worse individuals should be increased so that the evolution of the infeasible solutions can move toward the local best, which has the least intersecting segments. While for the feasible path, when the speed is constant, the total path length should be considered primarily. Since the idealized robot model is adopted in this paper, the kinetics effects are not considered, but the path should be smoothed as much as possible during planning. Based on the aforementioned hypothesis, a path is composed of N segments and N + 1 nodes; then the evolutionary variable can be represented as path i p 1 , p 2 , , p N + 1 where p 1 should be the fixed start point, and p N + 1 is the target point. Hence, the length of the path can be calculated as
f 1 = n = 1 N x n + 1 x n 2 + y n + 1 y n 2 = n = 1 N d x 1 , x N + 1 N 2 + y n + 1 y n 2
where x n and y n are the original coordinates of node p n ( x n , y n ) . d ( x 1 , x N + 1 ) is the distance from the start point to the target along the x -axis and y n is the transformed coordinates.
The smooth function can be defined as follows:
f 2 = π j = 2 N cos 1 x j x j 1 x j + 1 x j + y j y j 1 y j + 1 y j x j x j 1 2 + y j y j 1 2 × x j + 1 x j 2 + y j + 1 y j 2
which represents the summation of the angles between every connected segment. It starts from the second point and calculates the angle between the first and second segments, et cetera. The small f 2 value means a small direction change in each turn, representing a better path.
The overall objective function for the feasible function is expressed as follows:
F f e a s i b l e = 1 f 1 + k f 2 + C
where k is the weight of smoothing and a higher k can obtain a smoother path, but the diversity might be affected. C is the feasible and practical reward parameter which is a positive number that is not greater than the maximum path length in the current environment.
For the infeasible condition, except for the length of the path, the ratio of infeasible segments and the ratio of the infeasible path over the total length are involved in the estimation as follows:
f 3 = N inf N
f 4 = i = 1 N inf d o b s i f 1
where the N i n f is the number of infeasible segments and d o b s i is the overlapping distance between each infeasible segment and the obstacles. The objective function for the infeasible function is expressed as follows:
F i n f e a s i b l e = 1 w f 1 + 1 f 3 + f 4
where w is the weight to adjust the influence of total infeasible path length.

4.3. Design of the Cross Function

The main purpose of the multimodal strategy is to increase the efficiency of path planning, which can provide several collision-free paths from the start position to the target. Once a path falls into an emergency condition, the potential solution can change to another option immediately. The target and initial point for the robot are already known before the real application. The motion of the robot is from its current position to reach the next subsequent segment, and the process will continue until it reaches its goal position. Hence, segments of the robot should not intersect with obstacles or any other potential solution. For the design of the cross function, the objective is to generate a constraint that minimizes the arrival time for each planned path, and meanwhile, each segment cannot overlap with other paths as much as possible. Based on the above analysis of the algorithm, the best individuals in each sub-swarm are first mixed, and the cross fitness of each population elite is calculated based on the following equation, Equation (17), assuming that the number of elites is E.
c r o s s f u n c t i o n = e = 1 E n = 1 N C p n e = e = 1 E n = 1 N C p 1 1 C p N 1 C p 1 E C p N E
where the index n is the n-th segment in the current evaluated path, the index e is the e-th elite in the other sub-swarms, and C p n e is the cross value, which represents the cross fitness of individual i and interactive elite e. The cross function should be applied to each sub-swarm, and the cross evaluation should go through the sub-swarms except for the current swarm. Moreover, only the elite with the lowest number of segment intersections will be selected and mixed again with the worst individuals in each group memory. After that, all individuals are examined by the cross function again, and the elite who has the greatest difference from others in the mixed group is selected in the next iteration. It is worth noting that all mixing processes are evaluated only by cross fitness to avoid the influence of other factors.

4.4. Multimodal Path Planning Strategy

Multimodal path planning is designed to find multiple paths in a single run. When an emergency occurs on a path, the potential solution can be immediately changed to another option, thus addressing the inefficiencies of traditional path planning. Different from the scheme of a bug algorithm robot walking around obstacles when encountering obstacles and walking along a straight line without encountering obstacles [32], the multimodal path planning strategy adopts the segmented method as described in Section 4.1. Each node position of the segmented ( P 1 , P 2 , . . . , P N ) corresponds to a dimension of the variable; then the paraPA algorithm is used to find multiple paths. The optimal value with the best fitness that each swarm finally converges to is a global optimal path, while coevolution between multiple swarms ensures the diversity of the paths.
During the demonstration, we found that normally, there are one or two segments that are infeasible, causing the failure of the whole path, which is especially common among the superior populations at the late stage of convergence. To address this problem, this paper proposed a multi-path-based reverse planning strategy, which can be regarded as a complementary strategy for the replacement strategy to nectar sources in the bee colony. If the replaced nectar source only has one or two infeasible segments, the reverse planning strategy is triggered. Specific details of this strategy is shown in Algorithm 2:
Algorithm 2 Multi-path-based reverse planning strategy.
Input: Paths with segments that are infeasible.
Steps: for each infeasible path
              1. Find the index of the non-viable section (suppose S i is the infeasible segment).
              2. Go through the p b e s t group to check.
                      if there is a feasible solution in the same section ( S j )
                             Examine the segments from S j back to the start segment.
                                  if segments from start to S j in p b e s t p are all feasible
                                           Record the corresponding node in this p b e s t p and replace the corresponding positions in nectar.
                                   end
                        else
                              Trigger the random generation process of nectar.
                        end
                  end
Output: New feasible paths.

5. Experiments and Analysis

5.1. Environments and Comparison of Algorithms

In recent years, researchers have been using the population intelligence approach for path planning [33]. In order to examine the effectiveness of the ParaPA algorithm on different obstacle situations, three different types of scenarios are conducted in this paper to examine the robustness, efficiency, and convergence performance of the algorithm. The first is 20 × 20 maps, which are used to detect the sensitivity of algorithms to different types of obstacles. Second, rectangles are randomly generated in 50 × 50 maps with various sizes as obstacles, which is applied to test the adapted capacity of the algorithm under the complex environment. The third category, an actual scenario simulation, named Small Cultural Complex, is adopted as the application environment built-in 100 × 100 maps. The Small Cultural Complex is a new cultural industry model, which achieves different cultural function requirements by transforming the inside construction of the building, and during the transformation, mobile robots would take charge of the functional equipment transport, such as stage props or lighting. In this paper, the exhibition, sports activity, and theatrics functions are modeled, respectively. Meanwhile, to further show the advantages of the multi-swarm strategy, this paper uses PSO, CPSO, WPSO, ABC, and ACO algorithms for comparison experiments. Furthermore, MSPO and MABC are applied to test the effect of the proposed structure, also to evaluate the algorithm in terms of planning effectiveness, stability, etc.

5.2. Experimental Settings

For all algorithms, the data are calculated from 1000 runs, and the maximum evolution times are 200. For parameters setting in the PSO algorithm, the population size for each swarm is N P = 50 , and the maximum and minimum evolution steps are V m a x = 2 and V m i n = 2 , respectively. The acceleration factors are accepted as a constant 1.49445 in this paper. Specifically, for the WPSO, the inertia weight w is linearly changed from 0.1 to 0.9. Otherwise, the w is adapted as a constant 0.7. While for CPSO, the chaos-based initialization is utilized with a logistic chaos map to enrich the algorithm diversity where μ = 4 is adapted. For the parameter setting of the ABC algorithm, the look limitation is adapted as 5 × N P , which means all individuals in the swarm should visit the nectar more than five times, then the nectar can be replaced. In ABC and MABC, the size of onlookers is equal to the size of employed bees, but in the ParaPA algorithm, this number is equal to N P N P 5 5 . Moreover, the number of elites in each swarm is taken as 10 while 3 of them will be selected in the cross-verification process. Moreover, the number of elites in each swarm is taken as 10, while 3 of them will be selected in the cross-verification process. The start points for 20 × 20 and 50 × 50 maps are set at point S 0.5 , 0.5 and the target points are E 20 19.5 , 19.5 and E 50 49.5 , 49.5 , respectively. While for 100 × 100 maps, the start points are set to S 100 3.5 , 3.5 , and the target points are E 100 99.5 , 99.5 .

5.3. Results Representation and Analysis

5.3.1. Scenario 1: Path Planning on 20 × 20 Maps with Different Types of Obstacles

In the first case, different obstacles are constructed in a 20 × 20 map, as shown in Figure 5a–g to examine the sensitivity of the algorithm on obstacle shapes. The length of each chromosome is set to 12 for all algorithms. Results are shown in Table 1, and all the best values are highlighted in bold. First, by comparing with the particle swarm algorithm in the original coordinate, it can be clearly verified that the performance of the algorithm with the new coordinate system is significantly improved, but the problem is also ParaPA. As can be seen in the test of Map 4, the mean and variance under the original coordinates are the best, which means that the new coordinate system is inferior to dealing with the S-shaped trajectory planning. Because one dimension is fixed in the new coordinate system, the changeable range for each individual is less than the original coordinate. Hence, when in an environment with several continuous transverse obstacles, such as Map 4, the path generated by the algorithm is relatively monotonous, resulting in a worse convergence effect. It is worth noting that the best value of optima, as well as the worst individual in Table 1, are not in the original coordinate system, indicating that the algorithm in the new coordinate system could be superior during the evolution once the algorithm can find a feasible sample as soon as possible. From the results of other maps, the ParaPA algorithm has the best performance in 1000 runs, and the final convergence achieved the shortest path length in most scenarios. While the algorithm with the best stability is ACO, its evolutionary process is easily stagnant, resulting in unsatisfactory final convergence results. In comparison, the multi-swarm strategies, such as MABC and ParaPA algorithm, can increase the path diversity in a complex environment, and the proposed structure can be better exploited locally for better convergence results.
Table 2 shows the success rate. Note that the multiple population strategy proposed in this paper, combined with the parasitic relationship of two communities, can increase the probability of a successful search for the special continuous transverse obstacles in Map 4. In most cases, the algorithm can successfully find a suitable path. The results, however, are different in terms of the algorithm’s efficiencies combined with data from Table 1. The structure of ParaPA can utilize the population differences in multiple swarms and transform them into the diversity of path samples so as to improve the effectiveness of the algorithm, and the diversity of the algorithm is expressed by the individual disparity in the population in the following:
D I G = 1 N P i = 1 N P X i G 1 N P j = 1 N P X j G 2
where D I is the individual diversity [34], G means generation, and N P is the number of populations.

5.3.2. Scenario 2: Path Planning on 50 × 50 Maps with Randomly Generated Rectangles as Obstacles

While the merits of diversity are not obvious in the simple condition. To verify the convergence and searching capacity of different algorithms, more complex environments are built with randomly generated rectangles in 50 × 50 maps, as shown in Figure 6. In terms of the success rate from Table 2, ParaPA, ABC, and MABC are the only algorithms that can guarantee the finding of a feasible path, while Table 3 shows the convergence performance of the algorithm in the 50 × 50 environments. It is worth noting that the most stable algorithm is MABC, while the ParaPA algorithm makes some stability sacrifices in order to obtain a better search breadth in return for more accurate local convergence. This can be seen from the best value of optima in each map, indicating that the ParaPA algorithm has a stronger local searching ability. The standard deviation reflects the stability of the algorithm’s performance in a scenario, while the value of D I is an indicator of the diversity of the algorithm. In Table 4, the above algorithms with better performance are selected to compare the diversity. The standard deviations of CPSO and WPSO are larger compared with others. Meanwhile, from the performance of D I , the algorithm does not construct more feasible paths during chaos search, resulting in poorer sample diversity. The proposed structure performs chaos at the bottom to create more feasible paths, and the superior individuals are responsible for fine-tuning the optima. Combining the optima value in Table 3 and the D I performance in Table 4 shows the capacity of the ParaPA algorithm both in local convergence accuracy and global exploration.
The convergence curves in 50 × 50 maps for the compared algorithms are shown in Figure 7 to observe the convergence variation during an evolution where the maximum number of evolutions is set to 4000. The frequency and magnitude of the changes in the path length are small enough to be ignored after approximately 300 iterations. Hence only the first 400 iterations are taken into account, as the results shown in Figure 7a. It can be seen that all algorithms complete the planning for a feasible path within 60 iterations, after which the best value is planned locally around it, and the changing magnitude is reduced. Note that the best path is recorded based on the parasitism swarm, indicating that the best position has been locked by the inferior population at around 400 iterations. After that, what needs to be performed for the algorithm is to maintain the diversity to create more possibilities for path convergence. The D I is measured in 4000 generations to observe the tendency, as shown in Figure 7a. It can be seen that it is difficult to maintain a stable diversity for CPSO as well as WPSO, which means that their convergence tendency will fall into a unified position. ParaPA and MABC are more stable in terms of diversity performance.

5.3.3. Scenario 3: Path Planning for Real Application on 100 × 100 Maps

In practical map applications, as illustrated in Figure 8 [35], this paper chooses algorithms with a better performance in the above test to make a comparison.
The results in Table 5 and Table 6 show that the ABC series algorithm is able to create more sample diversity but hardly converts the recourses to the final convergence accuracy. For the ParaPA algorithm with a parasitic relationship, the bottom particles can guarantee the search ability, but the “nutrients” they can pass to the superior layer would be reduced. It can be solved by building multiple population interactions to increase the diversity in the upper layer population directly. From the test on Maps 2 and 3, MABC has a better performance in terms of stability. In other words, MABC does not have any internal mechanism to refine the local search, so the higher stability can be regarded as the outcome at the sacrifice of the exploitation ability as the optima result in Table 5. Consequently, the proposed algorithm has a better performance on most of the measurements, but how to choose an algorithm should consider the specific requirements of the real application. Furthermore, in order to show the results of multimodal path planning more intuitively, we take map (a) in Figure 8 as an example and show the multiple paths planned by the single-run algorithm in Figure 9. In a single run, it generates the three paths with the best fitness values, as shown in Figure 9.

6. Conclusions

Primarily, this paper introduces a dual-community-based evolutionary algorithm model that mimicked parasitic relationships in the biosphere. In the established structure, the bottom particles guarantee the algorithm diversity, while the superior community controls the better convergence resources and is responsible for completing the local convergence. Meanwhile, a multiple population strategy is conducted that is utilized to build information interaction channels among superior populations and directly shares the best information contained in superior populations. Finally, based on the above algorithm structure, ParaPA is proposed to solve the path planning problem. Meanwhile, we design a cross function to filter the high-level interactive information through the proposed multiple swarm framework so as to ensure path diversity in the superior population. From the comparison of the algorithm diversity and the average length of paths, the proposed approach is able to produce more path possibilities by using the structure, achieving a better solution in regard to path length.
Compared with traditional path planning strategies, such as A* algorithm, because of the characteristics of ABC and PSO algorithms, the ParaABC algorithm can effectively solve high-dimensional path planning problems. In the process of algorithm operation, once the target point is unreachable, the proposed algorithm can discard the current path directly or re-plan the infeasible section according to prior knowledge, thus saving performance loss. In addition, the ParaABC algorithm can plan multiple optimal paths through a single operation. Compared with the popular method of neural networks in recent years, the proposed method only needs to adjust a few parameters and has low equipment requirements. However, although the performance of ParaABC is improved compared with similar algorithms, it still faces the problems of diversity disappearance and premature convergence.
In the future, the three symbiosis-relation-based evolutionary structures introduced in this paper could be studied for different applications. Under the mutualism relationship, the convergence probabilities of multiple populations are equivalent, similar to the ring topology, which is suitable for multimodal or multi-objective optimization problems. For the commensal relationship, the interspecies relationship is inclined to a competition that has a continuous dynamic optimization capability. While in the parasitic relationship, where the convergence will tend to be homogeneous during the evolution, the avoidance of the premature problem is required. Note that the symbiosis relationship construction establishes an information bridge among multiple populations rather than a specific individual. Meanwhile, the development of symbiosis-relation-based evolutionary conception for different applications through the adapted framework is an interesting field that may be undertaken in the future.

Author Contributions

Conceptualization, H.R. and X.S.; Methodology, L.G. and X.S.; Software, L.G. and X.S.; Validation, H.R. and L.G.; Formal analysis, H.R. and X.S.; Investigation, H.R., L.G., X.S., and M.L.; Resources, H.R. and W.J.; Data curation, H.R., L.G., X.S., and M.L.; Writing—original draft, H.R., L.G., and X.S.; Writing—review and editing, H.R. and L.G.; Visualization, L.G., X.S., and M.L.; Supervision, H.R. and W.J.; Project administration, H.R.; Funding acquisition, H.R. and W.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Funding grant number 2018YFB1403703.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

This manuscript has not been published and is not under consideration for publication elsewhere. The authors have no conflicts of interest to disclose.

Appendix A

According to the description in [29], the commensalism and mutualism relationship should be defined as follows. All the relationships are defined under the multimodal optimization condition.
Commensalism relationship: The weaker population should give resources to the superior population while it is not affected by the interaction.
Mutualism relationship: The evolutional process will not be interrupted by others among two or more populations and can obtain benefits from group interaction. In other words, each evolution is independent while the information changing among different populations can contribute to their convergence.
We define the relationships in mathematics as:
Definition A1.
Suppose there are m optima in S, recorded as Φ 1 , Φ 2 , , Φ m . There exists S 1 S , S 2 S , S 1 S 2 . In the t-th iteration, S 1 and S 2 are recorded as S 1 t and S 2 t . If the following function is true
lim t i = 1 m P Φ i S 1 t > lim t i = 1 m P Φ i S 2 t
then the S 1 and S 2 is the commensalism relationship in the living space S.
Definition A2.
Suppose there are m optima in S, recorded as Φ 1 , Φ 2 , , Φ m . S 1 S , when the iteration is t and S 1 can be recorded as S 1 t , then P 1 t , 0 < P 1 t < 1 , satisfies
lim t i = 1 m P Φ i S 1 t = P 1 t
and for S 2 S , S 1 S 2 , P 2 t , 0 < p 2 t < 1 satisfies
lim t i = 1 m P Φ i S 2 t = P 2 t
For δ > 0 , always have
lim t P 1 t P 2 t < δ
called the S 1 and S 2 is the mutualism relationship in the living space S.

References

  1. Bounini, F.; Gingras, D.; Pollart, H.; Gruyer, D. Modified artificial potential field method for online path planning applications. In Proceedings of the 2017 IEEE Intelligent Vehicles Symposium (IV), Los Angeles, CA, USA, 11 June 2017. [Google Scholar]
  2. Yan, M.; Li, S.; Chan, C.A.; Shen, Y.; Yu, Y. Mobility Prediction Using a Weighted Markov Model Based on Mobile User Classification. Sensors 2021, 21, 1740. [Google Scholar] [CrossRef] [PubMed]
  3. Chang, J.R.; Jheng, Y.H.; Chang, C.H.; Lo, C.H. An Efficient Algorithm for Vehicle Guidance Combining Dijkstra and A* Algorithm with Fuzzy Inference Theory. J. Internet. Technol. 2015, 16, 189–200. [Google Scholar]
  4. Xiong, C.; Chen, D.; Lu, D. Path planning of multiple autonomous marine vehicles for adaptive sampling using Voronoi-based ant colony optimization. Robot. Auton. Syst. 2019, 115, 90–103. [Google Scholar] [CrossRef]
  5. Rashid, A.T.; Ali, A.A.; Frasca, M. Path planning with obstacle avoidance based on visibility binary tree algorithm. Robot. Auton. Syst. 2013, 61, 1440–1449. [Google Scholar] [CrossRef]
  6. Gonzalez, R.; Kloetzer, M.; Mahulea, D. Comparative study of trajectories resulted from cell decomposition path planning approaches. In Proceedings of the 2017 21st International Conference on System Theory, Control and Computing (ICSTCC), Sinaia, Romania, 19 October 2017. [Google Scholar]
  7. Song, B.; Wang, Z.; Zou, L. On global smooth path planning for mobile robots using a novel multimodal delayed PSO algorithm. Cogn. Comput. 2017, 9, 5–17. [Google Scholar] [CrossRef]
  8. Yen, C.T.; Cheng, M.F. A study of fuzzy control with ant colony algorithm used in mobile robot for shortest path planning and obstacle avoidance. Microsyst. Technol. 2018, 24, 125–135. [Google Scholar] [CrossRef]
  9. Rajput, U.; Kumari, M. Mobile robot path planning with modified ant colony optimisation. Int. J. Bio Inspired Comput. 2017, 9, 106–113. [Google Scholar] [CrossRef]
  10. Lamini, C.; Benhlima, S.; Elbekri, A. Genetic algorithm based approach for autonomous mobile robot path planning. Procedia Comput. Sci. 2018, 127, 180–189. [Google Scholar] [CrossRef]
  11. Zhang, Y.; Li, S.; Guo, H. A type of biased consensus-based distributed neural network for path planning. Nonlinear Dynam. 2017, 89, 1803–1815. [Google Scholar] [CrossRef]
  12. Yan, M.; Yuan, H.; Xu, J.; Yu, Y.; Jin, L. Task allocation and route planning of multiple UAVs in a marine environment based on an improved particle swarm optimization algorithm. Eurasip. J. Adv. Sig. Pr. 2021, 94, 2021. [Google Scholar] [CrossRef]
  13. Yan, F. Autonomous vehicle routing problem solution based on artificial potential field with parallel ant colony optimization (ACO) algorithm. Pattern Recogn. Lett. 2018, 116, 195–199. [Google Scholar] [CrossRef]
  14. Mo, H.; Xu, L. Research of biogeography particle swarm optimization for robot path planning. Neurocomputing 2015, 148, 91–99. [Google Scholar] [CrossRef]
  15. Montiel, O.; Orozco-Rosas, U.; Sepúlveda, R. Path planning for mobile robots using bacterial potential field for avoiding static and dynamic obstacles. Expert. Syst. Appl. 2015, 42, 5177–5191. [Google Scholar] [CrossRef]
  16. Oleiwi, B.K.; Al-Jarrah, R.; Roth, H. Multi objective optimization of trajectory planning of non-holonomic mobile robot in dynamic environment using enhanced GA by fuzzy motion control and A*. In Proceedings of the International Conference on Neural Networks and Artificial Intelligence 2014, Brest, Belarus, 3–6 June 2014. [Google Scholar]
  17. Ajeil, F.H.; Ibraheem, I.K.; Sahib, M.A. Multi-objective path planning of an autonomous mobile robot using hybrid PSO-MFB optimization algorithm. Appl. Soft. Comput. 2020, 89, 1–13. [Google Scholar] [CrossRef]
  18. Zafar, M.N.; Mohanta, J.C. Methodology for path planning and optimization of mobile robots: A review. Procedia Comput. Sci. 2018, 133, 141–152. [Google Scholar] [CrossRef]
  19. Gharehchopogh, F.S.; Gholizadeh, H. A comprehensive survey: Whale Optimization Algorithm and its applications. Swarm. Evol. Comput. 2019, 48, 1–24. [Google Scholar] [CrossRef]
  20. Yang, X.S.; He, X. Bat algorithm: Literature review and applications. Int. J. Bio Inspired Comput. 2013, 5, 141–149. [Google Scholar] [CrossRef]
  21. Contreras-Cruz, M.A.; Ayala-Ramirez, V.; Hernandez-Belmonte, U.H. Mobile robot path planning using artificial bee colony and evolutionary programming. Appl. Soft. Comput. 2015, 30, 319–328. [Google Scholar] [CrossRef]
  22. Neshat, M.; Adeli, A.; Sepidnam, G. A review of artificial fish swarm optimization methods and applications. Int. J. Smart. Sens. Int. 2012, 5, 107–148. [Google Scholar] [CrossRef]
  23. Qu, C.; Gai, W.; Zhang, J. A novel hybrid grey wolf optimizer algorithm for unmanned aerial vehicle (UAV) path planning. Knowl-Based. Syst. 2020, 194, 1–14. [Google Scholar] [CrossRef]
  24. Li, X.; Yang, G. Artificial bee colony algorithm with memory. Appl. Soft. Comput. 2016, 41, 362–372. [Google Scholar] [CrossRef]
  25. Wang, H.; Wang, W.; Xiao, S.; Cui, Z.; Xu, M.; Zhou, X. Improving artificial Bee colony algorithm using a new neighborhood selection mechanism. Inform. Sci. 2020, 527, 227–240. [Google Scholar] [CrossRef]
  26. Zhou, X.; Lu, J.; Huang, J.; Zhong, M.; Wang, M. Enhancing Artificial Bee Colony Algorithm with Multi-elite Guidance. Inform. Sci. 2020, 543, 242–258. [Google Scholar] [CrossRef]
  27. Cheng, M.Y.; Prayogo, D. Symbiotic organisms search: A new metaheuristic optimization algorithm. Comput. Struct. 2014, 139, 98–112. [Google Scholar] [CrossRef]
  28. Ezugwu, A.E.; Prayogo, D. Symbiotic organisms search algorithm: Theory, recent advances and applications. Expert. Syst. Appl. 2019, 119, 184–209. [Google Scholar] [CrossRef]
  29. Ren, H.; Shen, X.; Jia, X. A novel dual-biological-community swarm intelligence algorithm with a commensal evolution strategy for multimodal problems. J. Supercomput. 2021, 77, 10850–10895. [Google Scholar] [CrossRef]
  30. Tharwat, A.; Elhoseny, M.; Hassanien, A.E. Intelligent Bézier curve-based path planning model using Chaotic Particle Swarm Optimization algorithm. Cluster Comput. 2019, 22, 1–22. [Google Scholar] [CrossRef]
  31. Tuncer, A.; Yildirim, M. Dynamic path planning of mobile robots with improved genetic algorithm. Cluster Comput. 2012, 38, 1564–1572. [Google Scholar] [CrossRef]
  32. Melo W., D.; Jorge, D.; Marques, V. Low-cost thermal explorer robot using a hybrid neural networks and intelligent bug algorithm model. Int. J. Comput. Appl. Technol. 2021, 65, 245–252. [Google Scholar] [CrossRef]
  33. Nie, Z.; Yang, X.; Gao, S. Research on autonomous moving robot path planning based on improved particle swarm optimization. In Proceedings of the 2016 IEEE Congress on Evolutionary Computation (CEC), Vancouver, BC, Canada, 24 July 2016. [Google Scholar]
  34. Tang, L.; Dong, Y.; Liu, J. Differential evolution with an individual-dependent mechanism. IEEE. Trans. Evolut. Comput. 2014, 19, 560–574. [Google Scholar] [CrossRef] [Green Version]
  35. Liu, Z.; Tian, H.; Wang, Y.; Wu, Z. Talking about the Construction Requirements of Multi-functional Basic-level Cultural Service Complex. Eniertainment Technol. 2020, 10, 66–72. [Google Scholar]
Figure 1. Classification of mobile robot path planning approaches.
Figure 1. Classification of mobile robot path planning approaches.
Sensors 23 01751 g001
Figure 2. Evolutionary structure of the ParaPA algorithm.
Figure 2. Evolutionary structure of the ParaPA algorithm.
Sensors 23 01751 g002
Figure 3. The hierarchical illustration of the ParaPA algorithm with multi-swarm strategy.
Figure 3. The hierarchical illustration of the ParaPA algorithm with multi-swarm strategy.
Sensors 23 01751 g003
Figure 4. The coordinate transformation for path planning.
Figure 4. The coordinate transformation for path planning.
Sensors 23 01751 g004
Figure 5. Set of test 20 × 20 maps for the first scenario in which (ag) show the different distributions of obstacles on 20 × 20 maps.
Figure 5. Set of test 20 × 20 maps for the first scenario in which (ag) show the different distributions of obstacles on 20 × 20 maps.
Sensors 23 01751 g005
Figure 6. Set of test 50 × 50 maps for the second scenario in which (a,b) show the different distributions of obstacles on 50 × 50 maps.
Figure 6. Set of test 50 × 50 maps for the second scenario in which (a,b) show the different distributions of obstacles on 50 × 50 maps.
Sensors 23 01751 g006
Figure 7. Convergence comparison on 50 × 50 maps where (a) is the path length change curve during iteration and (b) shows the curves of the diversity of individuals ( D I ) within 4000 iterations.
Figure 7. Convergence comparison on 50 × 50 maps where (a) is the path length change curve during iteration and (b) shows the curves of the diversity of individuals ( D I ) within 4000 iterations.
Sensors 23 01751 g007
Figure 8. Real application on 100 × 100 maps of the Small Cultural Complex where (ac) are the layouts corresponding to the function of the exhibition, sports, and performance, respectively.
Figure 8. Real application on 100 × 100 maps of the Small Cultural Complex where (ac) are the layouts corresponding to the function of the exhibition, sports, and performance, respectively.
Sensors 23 01751 g008
Figure 9. Path diagram in which colorful lines represent three routes obtained by a signle run on a 100 × 100 map.
Figure 9. Path diagram in which colorful lines represent three routes obtained by a signle run on a 100 × 100 map.
Sensors 23 01751 g009
Table 1. Convergence performance comparison for scenario 1 on 20 × 20 maps.
Table 1. Convergence performance comparison for scenario 1 on 20 × 20 maps.
MAPs ParaPAPSO(ori)PSOCPSOWPSOMPSOACOABCMABC
Map1Mean28.531132.31130.05929.91930.07932.64438.00228.66728.583
Std0.1434.2452.7671.8511.7198.6400.0630.1700.042
Worst29.82558.06989.13040.22538.921212.31740.00031.32828.802
Optima28.45529.22828.45928.45928.45528.47538.00028.48028.480
Map2Mean27.02432.52628.04128.00728.00927.49638.03227.04427.033
Std0.0055.5961.1161.0771.0150.4920.2510.0180.006
Worst27.09853.12536.74934.03132.92630.31741.00027.17327.060
Optima27.02227.29527.02227.02227.02227.02537.00027.02227.022
Map3Mean28.26539.90531.82230.66830.54130.17138.00029.77728.670
Std0.4197.3006.3522.0841.8852.0560.0001.5420.571
Worst30.49584.38399.91741.31239.03542.67638.00037.59733.039
Optima27.84829.48227.88227.90027.91427.99438.00028.10228.170
Map4Mean57.44447.61772.30852.553NaN82.22181.61284.20673.934
Std14.6743.72314.3885.138NaN16.06917.74115.51413.658
Worst96.03769.024108.84660.659NaN114.122150.000113.900114.461
Optima41.73840.66546.06138.334NaN55.81450.00056.98948.293
Map5Mean28.92236.41131.51531.43031.53829.45538.00430.79529.431
Std1.7957.0862.3412.3702.2381.8450.0892.6351.627
Worst33.34369.18841.79142.17339.85536.34940.00039.65932.952
Optima27.29027.82827.29727.30227.29527.27538.00027.37627.395
Map6Mean28.48237.45529.82729.86829.78729.60038.00029.89929.117
Std0.3796.8851.4341.2761.0171.5980.0000.5320.119
Worst29.48763.96953.16937.79137.39973.79238.00032.87829.724
Optima27.98229.69628.04628.00128.04128.42138.00028.93629.111
Map7Mean30.91837.45233.14033.16932.95632.70038.03833.41430.018
Std1.9163.0873.1283.2383.0362.7300.2732.7950.388
Worst37.35354.30147.06946.97342.34429.08440.00040.82734.502
Optima29.13929.59429.15129.12729.13642.17538.00029.38529.285
1 The bold represents the best value.
Table 2. Success rate on various maps.
Table 2. Success rate on various maps.
MapsParaPAPSO(ori)PSOCPSOWPSOMPSOACOABCMABC
Maps 20 × 20Map110.9880.9980.9980.9990.998111
Map2111111111
Map310.4930.8710.940.9440.966111
Map40.49910.070.04100.040.830.0160.237
Map5111111111
Map610.3610.9440.9280.9260.995111
Map7110.9660.9690.9361111
Maps 50 × 50Map110.1310.6160.7940.6170.98111
Map210.3520.9580.7930.9470.9970.83611
Table 3. Convergence performance on 50 × 50 maps.
Table 3. Convergence performance on 50 × 50 maps.
MAPs ParaPAPSO(ori)PSOCPSOWPSOMPSOACOABCMABC
Map1Mean74.77694.88987.12596.96486.83380.842160.23382.59479.882
Std3.00521.59324.45839.28820.5847.85122.0513.9391.778
Worst103.862200.684295.924315.496219.431223.645266.000106.80686.220
Optima70.10772.95472.29973.70272.40370.300114.00074.30374.052
Map2Mean76.55891.57283.96397.12885.06278.558167.91680.01578.372
Std2.67114.28713.56638.19816.2074.34725.1262.8180.981
Worst83.936152.965255.782299.138238.817105.914320.00096.48087.981
Optima71.28771.95972.94674.96873.86771.361116.00075.30975.876
Table 4. Diversity of individuals ( D I ) on 50 × 50 maps.
Table 4. Diversity of individuals ( D I ) on 50 × 50 maps.
50 × 50 MapsParaPACPSOWPSOMPSOMABC
Map18.8975.0734.2796.7767.021
Map29.3387.5235.6187.0135.077
Table 5. Convergence performance on 100 × 100 maps.
Table 5. Convergence performance on 100 × 100 maps.
100 × 100 Maps ParaPACPSOMPSOABCMABC
Map1Mean147.727157.103155.820167.317161.270
Std5.7028.1307.9147.7174.940
Worst168.755203.790198.764206.522177.713
Optima141.006142.508142.262148.932148.287
Map2Mean145.878154.001153.771161.662157.390
Std9.35610.78610.3985.0134.021
Worst347.053263.723347.306188.500169.377
Optima139.750141.253141.306145.790145.765
Map3Mean166.445171.323176.595171.036169.858
Std18.0472.48612.65254.17580.481
Worst554.350175.755335.206255.509170.056
Optima165.459168.992165.795168.958169.838
Table 6. Diversity of individuals (DI) on 100 × 100 maps.
Table 6. Diversity of individuals (DI) on 100 × 100 maps.
100 × 100 MapsParaPACPSOMPSOABCMABC
Map17.1716.9739.24912.72813.409
Map25.6016.0348.10413.43014.664
Map34.0788.4249.12917.58310.043
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ren, H.; Gao, L.; Shen, X.; Li, M.; Jiang, W. A Novel Swarm Intelligence Algorithm with a Parasitism-Relation-Based Structure for Mobile Robot Path Planning. Sensors 2023, 23, 1751. https://doi.org/10.3390/s23041751

AMA Style

Ren H, Gao L, Shen X, Li M, Jiang W. A Novel Swarm Intelligence Algorithm with a Parasitism-Relation-Based Structure for Mobile Robot Path Planning. Sensors. 2023; 23(4):1751. https://doi.org/10.3390/s23041751

Chicago/Turabian Style

Ren, Hui, Luli Gao, Xiaochen Shen, Mengnan Li, and Wei Jiang. 2023. "A Novel Swarm Intelligence Algorithm with a Parasitism-Relation-Based Structure for Mobile Robot Path Planning" Sensors 23, no. 4: 1751. https://doi.org/10.3390/s23041751

APA Style

Ren, H., Gao, L., Shen, X., Li, M., & Jiang, W. (2023). A Novel Swarm Intelligence Algorithm with a Parasitism-Relation-Based Structure for Mobile Robot Path Planning. Sensors, 23(4), 1751. https://doi.org/10.3390/s23041751

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