• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A Review on Representative Swarm Intelligence Algorithms for Solving Optimization Problems:Applications and Trends

    2021-10-25 01:41:04JunTangGangLiuandQingtaoPan
    IEEE/CAA Journal of Automatica Sinica 2021年10期

    Jun Tang,Gang Liu,and Qingtao Pan

    Abstract—Swarm intelligence algorithms are a subset of the artificial intelligence (AI) field,which is increasing popularity in resolving different optimization problems and has been widely utilized in various applications.In the past decades,numerous swarm intelligence algorithms have been developed,including ant colony optimization (ACO),particle swarm optimization (PSO),artificial fish swarm (AFS),bacterial foraging optimization(BFO),and artificial bee colony (ABC).This review tries to review the most representative swarm intelligence algorithms in chronological order by highlighting the functions and strengths from 127 research literatures.It provides an overview of the various swarm intelligence algorithms and their advanced developments,and briefly provides the description of their successful applications in optimization problems of engineering fields.Finally,opinions and perspectives on the trends and prospects in this relatively new research domain are represented to support future developments.

    I.INTRODUCTION

    IN nature,groups of thousands,millions,or trillions of individual elementary entities can self-organize into multifarious forms to fit a functional objective,purely through local and ordinary interactions.Examples at different biological levels include the self-assembly of crystals of bacterial flagella at the molecular scale,regular development of multicellular organisms at the cellular scale,and interconnected food searching or risk avoidance (for energysaving in ants) at the colony scale [1].Research on biological behavioral intelligence has been conducted for a long time;a study regarding Cape bees and“animal intelligence”was published in 1883 [2].

    Throughout nature,different organisms often profit from acting in swarms.Models of the simple collective behavior of individuals indicate that uncomplicated partial interactions are sufficient to produce and present group morphologies.Individuals constitute groups that can collectively process information and provide corporate decisions.By shared information,the group could make better decisions than a single individual;this phenomenon is called as“collective intelligence”[3].For example,bees appear to rely on decentralized decision-making for the critical choice of a new nest site.Although no individual evaluates the complete set of available information or directly compares the available options,the swarm efficiently integrates the resulting flow of information into a high-quality final decision,without central control.

    By studying the characteristics of individuals and their relationships with groups,algorithms for the corresponding mechanism,known as“swarm intelligence”have been formalized.These essentially comprise biologically inspired computations and have been deemed an emerging field and an integral part of artificial intelligence (AI) [4].Swarm intelligence algorithms can be organic combination of evolutionary computing,artificial neural networks (ANNs),and fuzzy systems in a series of computational intelligence(CI) approaches [4].They are known as a group of natureinspired methods that are used to deal with complex optimization problems for which mathematical or traditional approaches are invalid.This ineffectiveness is mainly caused by,e.g.,the overwhelming complexity of the procedure for mathematical reasoning,possible uncertainties,or a procedure that is stochastic by nature.Among the heuristic algorithms inspired by nature,in addition to the swarm-based intelligent optimization algorithms mainly introduced in this paper,some are also inspired by natural evolution,physical rules,and human behavior.The hierarchical distribution of these algorithms is shown in Fig.1.Specifically,we divide the heuristic algorithms inspired by nature into the above four categories.In each category,we have listed five classic algorithms.The first is the method of swarm intelligence inspired by the intelligent behavior between biological swarms,which will be discussed in this review.Secondly,algorithms based on the laws of cosmic physics are also an important branch,including the well-known simulated annealing algorithm,gravity search algorithm,sine cosine algorithm,artificial electric field algorithm,etc.Then,evolutionary algorithms based on evolutionary theory are produced by using the laws and processes of biological heredity and population evolution in nature.The most popular is the genetic algorithm that simulates Darwin’s evolution.In addition,classic algorithms include differential evolution,evolution strategy and evolution programming.Finally,there are algorithms based on human intelligent behavior,for example,the harmony search algorithm that simulates music creation by musicians,the teaching learning-based optimization algorithm imitating the influence of teachers on students,the political optimization algorithm that imitates the multi-stage process of politics,etc.

    Fig.1.The hierarchical distribution of nature-inspired algorithms.

    The term“swarm intelligence”was formally put forward by Beni and Wang in the background of studying a cellular robotics system [5].Currently,it is considered as a subdiscipline of CI,aiming to solve problems by modelling populations of agents that can self-organize and interact with each other.In view of well-organized principles,we have selected many of the most representative swarm intelligence algorithms and divided them into six categories,based on the number of algorithms used and the basic biology they are derived from:ant colony optimization (ACO),particle swarm optimization (PSO),artificial fish swarm (AFS),bacterial foraging optimization (BFO),artificial bee colony (ABC) and other swarm intelligence algorithms.As swarm intelligence algorithms possess features of self-organization,parallel operations,distributive operations,flexibility,and robustness,they have been gradually very widespread and are utilized in many events.For example,swarm intelligence algorithms are used for scheduling problems,robots,power systems,parameter optimization,system identification,image processing,signal processing,and other optimization problems of engineering fields.Hence,studying the swarm intelligence algorithms possesses extremely significant academic and practical value.

    From the number of publications based on Google Scholar data since 2000,we can obtain that abundant research has been conducted on swarm intelligence algorithms in the past two decades,as shown in Fig.2.Overall,from 2000 to 2005,the development of the various algorithms for swarm intelligence was relatively slow;this period can be regarded as the initial stage.From 2005 to 2012,the publication of various algorithms in Google Scholar contributed to a rapid growth trend,meaning that the swarm intelligence algorithms developed rapidly.After 2012,the related research was still increasing,but the trend was relatively flat.Specifically,among the five swarm intelligence algorithms,PSO has been the most researched;it developed particularly rapidly from 2005 to 2012.However,the growth rate slowed down after 2012,and particularly,the publications have continued to slightly decline since 2016.It is likely that the AI experts and scholars have gradually paid attention to machine learning,tree searching,reinforcement learning,and deep neural networks in recent years [6].The research on ACO is second only to that on PSO,and the overall research shows a steady growth trend.Unlike the other algorithms,whose growth rates slowed down after 2012,the ABC algorithm has become more pronounced since 2012,owing to the popular research topic of robot swarms.The AFS algorithm has had relatively fewer academic achievements over the years,but is the most stable.Academic publications regarding the BFO algorithm were particularly prominent in 2012,as more academic conferences were held on related topics.In addition to the above five types of algorithms,there are many other swarm intelligence algorithms,and the corresponding research is also increasing each year.

    Fig.2.Google trends indicator of swarm intelligence algorithms from 2000 to 2020.

    This review is intended to have reference significance for professional and technical researchers interested in the stateof-the-art in swarm intelligence algorithms and their trends in related fields.The remainder of this paper is organized as follows.Section II divides the swarm intelligence algorithms into six categories.Section III describes the achievements of the swarm intelligence schemes in various applications of related fields;several results from these applications are provided and are compared with the pipelines of commonly utilized approaches.Section IV summarizes major strengths and limitations in swarm intelligence algorithms,along with main trends that would be developed in the future.Section V concludes this review.

    II.SWARM INTELLIGENCE ALGORITHMS

    From flocks of birds to schools of fish to swarms of insects,biological aggregations displaying collective behavior emerge throughout nature.This self-organized,natural phenomenon has been the subject of intensive modelling for decades,as it is highly commonly employed in nature to efficiently solve problems and owing to its potential ability as a biomimetic control strategy for engineering systems [7].In recent years,swarm intelligence has been universally studied in various areas;consequently,many related approaches to collective behavior have emerged.A categorization of swarm intelligence algorithms with several representative works,are illustrated in Fig.3.In this section,we will briefly review each of these swarm intelligence algorithms,along with their most recent developments and applications.

    Essentially,swarm intelligence algorithms are iteratively based on stochastic search algorithms,where heuristic information is shared to conduct the search for the following iterations.Fig.4 shows a general framework for swarm intelligence algorithms [8].Before the initialization phase,the parameter values are required to be defined.The evolutionary process then begins with initialization and its corresponding strategies.Then,the termination condition used to stop the algorithm execution is set;it is either one condition or comprised of two separate conditions.The fitness function that is in charge of the evaluation of the search agents is evaluated;the fitness function can be a basic metric,or a composed one.The algorithm updates agents until the previous termination condition is achieved.Finally,the best search result is generated.For a specific swarm intelligence algorithm,the sequence of each process may be different,and some processes could be executed several times in a single iteration.

    A.Ant Colony Optimization

    In computer science and operational research,ACO is a probabilistic method for resolving computational problems which could be reduced to searching the optimal paths through graphs.It was first proposed by Marco Dorigo in 1992 based on his Ph.D.research [9],and the inspiration for this basic algorithm came from the behavior of ants finding paths in the procedure of food searching or risk avoidance.The ACO algorithm has characteristics of distributed computing,positive information feedback,and a heuristic search [10].Its basic idea is to imitate an ant’s dependence on pheromones,and to guide each ant’s action by using positive feedback among the ants.This algorithm can be utilized to deal with most optimization problems or can be transformed into optimization problems.Ants can leave a type of material called a pheromone in the paths they pass.They can sense the intensity of this material and can thereby guide their own direction of action in the process of foraging [11].The basic searching process of ACO is as follows [12].

    Fig.3.A categorization of the swarm intelligence algorithms and their representative works.

    Fig.4.General framework of the swarm intelligence algorithms.

    At first,mants are randomly placed in the location.The initial values of the pheromones on each path are equal,that is τij(0)=τ0.Then,thek-th(k=1,2,...,m) ant chooses the next location to be transferred to on account of the random proportion rule,and its selection probability is shown as follows [13]:

    In the above,τijis the pheromone of edge;ηij=1/dijis the heuristic factor of transferring from locationito locationj,wheredi jis the distance;andallowedkis the location set that thek-th ant is allowed to visit next.

    After thet-moment,all ants have accomplished a tour.The path length of each ant is calculated,and the shortest path length is saved.At the same time,the pheromones on each edges are updated.The pheromone update includes two parts[14]:i) the volatilization of pheromones,as in (2);and ii) the release of pheromones by ants on the edge they pass by (3)

    In the above,ρ is the volatilization coefficient of the pheromone,is the pheromone released by thek-th ant passing the edge(i,j).After the ants complete a cycle,they return to the original location and prepare for the next tour.

    The original ACO concept has been diversified to manage a wider group of numerical problems.However,several new problems have emerged,and additional studies have been required to improve the ACO performance,depending on different aspects of the ant behaviors.TheMAX–MINant system [15] is a typical representative method for directly improving an ant system.This algorithm modifies the approach to updating pheromones.After each iteration,only one ant can update the pheromone to obtain a better solution.To avoid the search stagnation,the pheromone on the path is limited to [MIN,MAX].In addition,the initial value of the pheromone takes its upper limit,helping to promote the algorithm’s search ability.Bullnheimeret al.[16] provided a novel rank-based version of an ant system.After each iteration,the paths of the ants are arranged from small to large.The algorithm gives different weights on account of the path length;the shorter the path length,the greater the weight.Huet al.[17] proposed a new ACO denoted“continuous orthogonal ACO”for solving continuous problems.The pheromone deposition mechanism of this algorithm enables ants to efficiently search for solutions.Guptaet al.[18]introduced the concept of depth for a recursive ACO.The depth determined the number of recursions,and each depth was based on an ordinary ant colony algorithm.Gaoet al.[19]introduced a K-means clustering algorithm to combine with ACO algorithm with three immigrant schemes for addressing the dynamic location routing problem.Hemmatianet al.[20]applied an elitist ACO algorithm to the multi-objective optimization of hybrid laminates,to obtain the minimum cost in the calculation process.

    B.Particle Swarm Optimization

    In groups of birds acting in a cooperative energy-saving manner,individuals in the group continuously optimize the search mode,in view of the experience of group members.Based on a study of swarm behavior,Kennedy and Eberhart put forward the PSO algorithm in 1995 [21].In PSO,a flock of particles retain motion in a constrained parameter space,interact with each other,and update their velocities and positions on account of their own and their neighbors’information,finding the global optimum [22],[23].The initial state of the algorithm is a group of random particles comprising a particle swarm;it also represents a random solution space.Each particle owns only two properties:velocity and position.Each particle searches for the optimal solution results in the corresponding space,and continuously updates its individual optimal value.Particle swarm members share their individual optimal values.The global optimal solution is the best of the current individual optimal values[24],[25].Based on the current individual optimal value and the current global optimal solution,the particles continuously adjust the velocities and positions.The above search process can be formulated as (4) and (5) [26],as follows:

    In the above,vtandxtare the velocity and position of particleiat timet,respectively.pi,tis the individual optimal value of particleifound before timet,andpg,trepresents the global optimal solution of the particle swarm found before timet.c1,c2,c3are the respective cognition coefficients,andr1,r2are the random parameters within [ 0,1].

    The first portion of (4) is the“memory term”and indicates the particle’s velocity vector in the current state.The second portion of (4) is a“self-cognition term”,which is a vector from the current value to the individual optimal value.It expresses the influence of the particle’s own experience.The third portion of (4) is a“global-cognition term”,i.e.,the vector from the particle’s current value to the current global optimal solution.It reflects the cooperation and information sharing between particles.

    Because of its simplicity,effectiveness,and low computational cost,PSO has obtained significant popularity,and many improvements have been proposed.Most researches of PSO improvement fall into categories which modify the model coefficients,consider the population structures,and alter the interaction modes [27].By combining PSO and a genetic algorithm,Ghamisi and Benediktsson [28] provided a new feature selection method.This method also adopts a support vector machine (SVM) as a fitness function,which improves the classification accuracy.To address the challenge of premature convergence to a local optimum in global optimization problems,Liet al.[29] introduced a competitive and cooperative PSO with a data-sharing mechanism.It allows particles to share each other’s individual optimal values,instead of simply reflecting the cooperation in the current global optimal solution.Furthermore,the comprehensive learning PSO proposed in [30] allows members in the swarm to share the historical individual optimal values of all other particles.This algorithm is particularly outstanding in solving multimodal problems.In[31],a parallel PSO algorithm with a special operator was proposed,and it has proven effective in solving the problem of PSO convergence at a suboptimal location.Since the PSO algorithm was proposed,its convergence has been an important aspect of the research.Via the depth analysis of basic PSO algorithm,Zeng and Cui proposed the stochastic PSO algorithm [32],which can ensure the probability of converging to be globally optimal.In contrast to the traditional PSO,which focused on solving continuous problems,Clerc [33] studied the reasonable adoption of a discrete PSO to the classic traveling salesman problem.The results represent that the discrete PSO algorithm could be easily combined with other algorithms to efficiently solve certain discrete problems.In [34],a PSO-based neural network identifier was constructed to model the unknown system dynamics,and proved to be more effective in resolving the Hamilton-Jacobi-Bellman equation by utilizing the estimated system states.Royet al.[35] proposed a PSO-based artificial neural network for enhanced forecasting of the reliability of software,considering the fault generation phenomenon with the fault complexity of different levels.Lvet al.[36] combined a surrogate-assisted PSO with Pareto active learning to solve the multi-objective optimization problem with high computational cost,and in this research PSO was used as a sampler to generate candidate solutions.

    C.Artificial Fish Swarm

    In nature,fish could discover more nutritious areas by individual searching or following other fish,and areas with more individuals are normally more nutritious.The stochastic population-based AFS optimization algorithm was initially proposed by Liet al.in 2002 [37].Its main idea is to imitate fish behaviors such as preying,swarming,and following,with a local search of individuals for generating a global optimum[38].The AFS algorithm possesses attractive features similar to those of a genetic algorithm,but it can reach a faster convergence speed and requires the adjustment of fewer parameters.The AFS algorithm does not include crossover and mutation procedures,so it can be executed more easily.It is also an optimizer in view of population;the system is initialized with a group of randomly generated potential solutions,and then iteratively executes a search for the optimum one [39].The AFS algorithm acts as a parallel and random search optimization method,with the advantages of a strong global search capacity and an excellent convergence speed.The fish group is represented by a set of points or solutions,with each agent representing a candidate solution.The feasible solution space constitutes the“waters”where the artificial fish moves and searches for the optimum.The survival,competition,and coordination mechanisms are introduced into the basic operation [40].The main challenge in the algorithm implementation is determining how to use simple and effective ways to construct and realize these behaviors.In addition,the PSO algorithm and the AFS algorithm are often used for comparison in research:the PSO algorithm is inspired by a flock of birds,and it uses some particles flying in the search space to find the best solution,that is,these particles continuously adjust themselves based on the original path according to the position of the optimal solution in the particle swarm;the AFS algorithm is mainly based on the foraging behavior of the fish school in the real environment with a bottom-up design idea,and it mainly uses the three operators of artificial fish school for foraging,clustering and tail-finding to construct the underlying behavior of the individual.

    The description of basic AFS processing in a global optimization problem for minimizingf(x)x∈Ω is as follows[41],[42].There exist a population of artificial fish includingnagents or pointsxi,i=1,2,...,nandxiis a floating-point encoding with the lower and upper boundsliandui.The search space ofxiis: Ω={x∈Rn:li≤xi≤ui,i=1,2,...,n}.Each pointxihas a closed neighborhood with a radiusv,which is called the“visual scope”,andv=δmaxj∈{1,2,...,n}(uj?lj).In addition,δ ∈(0,1) is a positive visual parameter which can be decreased by the iterative process.The“visual scope”of pointxiisnpinumber of points.The crowding factorCf=npi/Nindicates the degree of crowding.Ifnpi=0,the“visual scope”is empty,and pointxiwill move stochastically and look for better areas;ifCf≤θ,the“visual scope”is not crowded,and pointxican chase or swarm towards the best or central point of the“visual scope”.IfCf≥θ,the“visual scope”is congested,and the pointxihas a problem insofar as chasing.It will randomly choose and move toward another point.In the AFS algorithm,the swarm behavior is featured by the motion to the central point of the“visual scope”.

    During the use of the AFS algorithm,it has been developed and improved by numerous researchers,generating a series of novel variants.In [43],a binary version of the AFS is proposed for resolving 0–1 multidimensional knapsack problems.In this method,a binary string of 0/1 bits represents a point.Through copying the relevant bit from itself or from another specified point with identical probability,each bite of a trial point is produced.Zhanget al.[44] proposed a Paretoimproved AFS algorithm for resolving a multi-objective fuzzy disassembly line balancing problem,and it was useful in decreasing the existing uncertainty.To heighten the global search ability and the convergence speed,Zhu and Jiang [45]proposed a novel quantum AFS algorithm.The algorithm is on account of principles and concepts of quantum computing,e.g.,the quantum bit and quantum gate.As regular artificial fish can get stuck in blindness hunting,Gaoet al.[46]provided a knowledge-based AFS algorithm with crossover.In addition,the conventional regular AFS could only perform local or global searching by continuously initializing the step and visualization parameters.To remedy this,Yazdaniet al.[47] proposed two novel adaptive methods for utilizing fuzzy systems to adjust the step and visualization parameters during execution.The purpose is to adaptively adjust the capability for global and local searching.In the first one,the visualization and step parameters of all fish are adjusted uniformly,whereas in the second,a special fuzzy controller controls the step and visualization parameters for each artificial fish.

    D.Bacterial Foraging Optimization

    The BFO algorithm was first proposed by Passino in 2002[48].It models the foraging behavior of Escherichia coli in humans,is relatively new to the family of nature-inspired optimization algorithms.It optimizes through the competition and cooperation among bacterial populations and is utilized as a global random search algorithm.It has characteristics of simplicity and fast convergence,and there is no need to optimize the gradient information of the object during the optimization process [49].The BFO simulation of the bacterial population contains three steps: chemotaxis,reproduction,and elimination/dispersal,corresponding to the following three processes [50]:in chemotaxis,a bacterium moves by carrying on small steps while searching for nutrients,and a core concept of BFO is simulating the chemotactic movement of the virtual bacteria in the corresponding search space;then,in the reproduction step,the remaining bacteria multiply to keep the population size;finally,the local area where the individual bacteria live may suddenly change,which may lead to elimination/dispersal,i.e.,the collective death of the bacterial population staying in this local area or the collective migration to another local area[51].The new individual randomly created by the migration operation may be nearer to the global optimal solution,which is more valuable to jumping out of a local optimal solution.The main procedures in basic BFO are summarized as follows[48],[49].

    The size of the bacterial population isS,and the location of the bacteria indicates a candidate solution to this problem.The information of a bacteriumiis labelled with aD-dimensional vectorA chemotactic step is deemed as a tumble followed by a tumble,or a tumble followed by a run.In addition,θi(j+1,k,l) indicates the position of the bacteriumiafter thej-th tropism operation,k-th replication operation,andl-migration operation.Moreover,C(i)>0 represents the step size of forward swimming,and Φ(j) represents the unit direction vector (randomly selected after rotation).If the fitness at θi(j+1,k,l) is better than that at θi(j,k,l),Φ is kept unchanged,and the bacteria keep to swim in this direction until a position with the best fitness is discovered or the set number of tendencies is reached;otherwise,a new Φ is generated,and the next rotation is performed.

    When searching for food in their own way,each bacterium also receives attractive signals from other individuals in the population,and receives repulsive signals from nearby individuals for maintaining a safe distance.P(j,k,l) is the position of each bacterium in the populationSafter thej-th tropism operation,k-th replication operation,andl-migration operation.Jcc(θ,P(i,k,l)) indicates that the synthesized effects of attraction and repulsion among the bacteria are considered simultaneously.

    In the BFO algorithm,its population size does not change after a copy operation.The number of bacteria eliminated is usually set asSr=S/2.First,the bacteria are sorted according to their pros and cons;then,theSrbacteria that are ranked lower are eliminated.The remainingSrbacteria will selfreplicate.A newly generated individual has the same position as the original individual.The migration operation appears with a certain probability.If a bacterial individual meets the probability of migration,the bacterial individual would die,and a new individual would be randomly created at any position in the solution space.Destroyed individuals may own different positions.

    The BFO algorithms also have some extensions for subsequent research.Majhiet al.[52] proposed a highperformance forecasting model for the prediction of various stock indices,and they minimized the mean square error to optimize the connecting weights of the adaptive linear combiner.In [53],Niuet al.described a multi-objective BFO and introduced two different performance metrics,diversity and generational distance,to evaluate multi-objective optimization problems.Chenet al.[54] developed a multicolony BFO algorithm in which a single population bacterial foraging method was applied to the interacting multiple colony model (which links the chemotaxis behavior of individual bacterial cells to the intercellular communications of bacterial communities).Several improvements have been made by combining bacterial foraging optimization algorithms with other algorithms.Gollapudiet al.[55] developed a velocity-modulated BFO technique by combining the BFO algorithm and a PSO technique.Kimet al.[56] presents a hybrid technique involving genetic algorithms and BFO algorithms,focusing on mutations,crossovers,step changes,chemotaxis steps,and bacterial lifespans.

    E.Artificial Bee Colony

    A group of honey bees,called as a swarm,could complete different tasks through social collaboration.Karaboga proposed the ABC algorithm in 2005 to solve a multivariate function optimization problem,in view of the intelligent foraging behavior of a honey bee swarm [57].Bees form groups of insects,and although the behavior of a single insect is extremely simple,a group of single individuals represent very complex behavior.Real bee populations can collect nectar from food sources (flowers) with high efficiency in any environment; they can also simultaneously adapt to environmental changes [58].In the original ABC algorithm,it contains three kinds of individuals:employed bees,onlooker bees,and scout bees.Each employed bee corresponds to a certain food source (solution vector),and the field of the food source is iteratively searched [59].More formally,employed bees will search for a richer food source (dm) near the original food source (xm).They will evaluate the fitness and compare it with that of original food source.They would find a new food source by using a process corresponding to (8) [60],as follows:

    Each food sourcexmis a solution vector to a problem withncomponents(xmi,i=1,2,...,n);and ωmiis a random number within a certain range.The bees choose the food source based on the fitness value,using a greedy selection.The fitness can be calculated as follows [61]:

    In the above,fm(·) is the objective function for the optimization problem.

    In the onlooker bee phase,an onlooker bee selects a food resource,based on the probability values computed by the fitness values of the employed bees’ food resources.A roulette method based on the probability value is employed to observe and collect food with the onlooker bees.The probability value indicates how likely a food resource is to be selected by an onlooker bee,and it can be computed by using(10) [62],as follows:

    Here,Mindicates the population size of the food resources.During the scout bee phase,if the food source is not updated multiple times,the food source is abandoned,and the employed bee turns to a scout bee to randomly search for another food source [63].To prevent this method from falling into a local optimum,when the food source iteration limit is not improved,the food source is recorded in a taboo table.Simultaneously,the employed bee corresponding to the honey source is transformed into a scout bee,to randomly produce a new position for finding a new food resource [64].

    Based on observing the operations and structure,there are also drawbacks to the basic ABC algorithm.For example,the artificial bee could only move straight to one of the nectar sources,and the characteristic may decrease the set of explored zones.Hence,several novel ABC algorithms have been proposed to address some discovered disadvantages.The searching mechanism can be improved by using an efficient genetic selection method [65].This new discrete binary ABC algorithm can solve the dynamic clustering problem,and the number of clusters is determined automatically.In [66],a Gbest-guided ABC and global bee colony search are improved by taking the basis of the maximum fitness values instead of the maximum cycle numbers,and the results are superior.Panet al.[67] improved the method for generating food resources by using a self-adaptive strategy to ensure that the discrete ABC algorithm could be adopted to discrete spaces.The improvement in the local search approach is also shown to be effective in enhancing the local intensification capability.Its convergence speed can be enhanced by applying a control parameter,as discussed in [68].Experiments with this method have proven the solution quality and the convergence to global optimum.In [69],Karaboga and Gorkemli provided a combinatorial ABC algorithm that showed good performance in solving complex optimization problems regarding finding a Hamiltonian path with the lowest cost.Jiet al.[70] proposed a scale-free ABC algorithm that applied an innovative method for a hybrid ABC algorithm with a scale-free network to enhance the optimization performance without sacrificing efficiency.

    F.Other Swarm Intelligence Algorithms

    Other swarm intelligence algorithms are still mainly derived from the simulation of natural biomes,especially those based on insects and animals.They are optimized using artificial algorithm simulations of the processes of finding food or exchanging information in biomes,and most of them are directional iterative methods based on a probability search[71].With the development of swarm intelligence,the optimization algorithms are not limited to introducing biological population characteristics into the algorithm.Some studies have introduced human biological characteristics into swarm intelligence algorithms,such as the human immune mechanism [72].The new swarm intelligence algorithms and improvements to classic algorithms proposed in recent years mainly focus on reducing parameters,simplifying processes,increasing computation speeds,and improving search capabilities,and are aimed at high-dimensional and multiobjective optimization problems [73]–[76].The fields of application for swarm intelligence algorithms are becoming increasingly extensive,and the expansion of algorithm applications has also guided the development directions.

    In addition to the several classic swarm intelligence algorithms,there are many other extended algorithms that have been widely discussed and employed.In view of the obligate brood parasitic behavior of several cuckoo species in combination with the Lévy flight behavior of several birds and fruit flies,a meta-heuristic algorithm defined as the“cuckoo search”was formulated to deal with optimization problems[77].Pigeon-inspired optimization is a swarm intelligence optimizer utilized to resolve air-robot path planning problems.It presents a map and compass operator model on account of a magnetic field and the sun,whereas a landmark operator model is constructed on account of landmarks [78].The bat algorithm is an animal group/herd-based learning algorithm which utilizes the echolocation behavior of bats to generate solutions for single-and multi-objective domains within a continuous solution space [79].The grey wolf optimizer is a swarm intelligence algorithm inspired by grey wolves,and simulates their leadership hierarchy and hunting mechanisms in nature [80].An artificial immune system (ARTIS) contains several typical properties of natural immune systems,such as diversity,distributed computation,error tolerance,dynamic learning and adaptation,and self-monitoring.ARTIS acts as a complete framework for a distributed adaptive system and could,in principle,be utilized in different domains [81].The fruit fly optimization algorithm (FOA) is a method of seeking global optimization based on the foraging behavior of fruit flies.Fruit flies can rely on their outstanding sense of smell and vision to find better food sources and locations where their companions gather [82].Glow-worm swarm optimization is developed by incorporating the behavior of glow-worms into AI systems.It can be used for the simultaneous calculation of multiple optima of multimodal functions [83].Inspired by colonizing weeds,invasive weed optimization is a numerical stochastic optimization algorithm for mimicking the robustness,adaptation,and randomness of colonizing weeds,using an uncomplicated yet effective optimization algorithm [84].

    III.APPLICATIONS AND RESULTS

    The appearance of swarm intelligence algorithms has provided fast and reliable methods for obtaining solutions to different complex optimization problems.These swarm intelligence algorithms,including ACO,PSO,AFS,BFO,ABC,and other algorithms,reflect their ability in resolving nonlinear design problems in real-world applications thinking over nearly all areas of science,engineering,and industry.These applications include,e.g.,those concerning transport problems (such as those concerning unmanned aerial vehicles(UAVs)),network routing,route planning,robot systems for scheduling problems,power systems,fault diagnosis,parameter optimization,system identification,cluster analysis,data mining,image processing,layout optimization,and signal processing.The purpose of this section is to present some typical and relevant professional applications in which the use of swarm-based optimization algorithms has successfully been made in the literature (but not to present all possible applications).

    Fig.5.Various applications of swarm intelligence algorithms in different fields.

    Fig.5 depicts the various above-mentioned applications for swarm intelligence algorithms in different fields,based on the published studies explored in Google Scholar.If one paper satisfies applications of swarm intelligence algorithms,it is selected in this review for further introduction.Some highfrequency application domain keywords are used to refine literature classification.Based on this kind of search strategy and criteria,appropriate publications are recorded and reviewed.Among them,the direction of the arrow indicates the application of that type of swarm intelligence to resolve the corresponding problem;the line thickness is based on the number of studies found.The information shown in Fig.5 is basically consistent with that in Fig.1,e.g.,PSO and ACO are the two most proverbially used algorithms.In addition,we can find that scheduling problem,robot system,power system,parameter optimization,system identification,image processing,and signal processing utilize the swarm intelligence algorithms the most;these seven applications are explained in further detail below.To ensure the quality of review,it should be emphasized that the principle of adopting the corresponding swarm intelligence algorithm literatures in each type of application is that the citation rate is relatively high.In addition,the fan diagrams in Fig.6 are used to illustrate the proportion of representative applications of different swarm intelligence algorithms,and this format refers to [85].

    A.Scheduling Problem

    Scheduling acts as the procedure of arranging,controlling,and optimizing task,resource,workloads,etc.in a production procedure,manufacturing module,or mathematical computation [86].Most research on scheduling problems aims to maximize the efficiency of operations and to reduce costs.For example,in manufacturing,the objective of scheduling is to minimize the production time and costs by instructing a facility when to produce a product,with which staff,and on which equipment.

    Fig.6.Statistical diagrams to show the proportion of applications of swarm intelligence algorithms.

    The resource-constrained project scheduling problem(RCPSP) acts as one of the most investigated kinds of scheduling problems;its objective is to generate the resourceto-task assignments that make a finite project plan the cheapest or shortest.Myszkowskiet al.[87] studied the usability and robustness of the ACO and its hybrids with priority rules in resolving the RCPSP,through renewing pheromone values in view of both the best and worst solutions preserved by the ants.Based on PSO,García-Nietoet al.[88]determined successful cycle programs for traffic lights and applied the optimization technique for two extensive traffic networks placed in the metropolitan centers of Spanish cities.Pan [89] presented a new cooperative co-evolutionary ABC algorithm for addressing a novel steelmaking continuous casting (SCC) scheduling problem in the context of iron and steel production processing.Its effectiveness was confirmed using instances from real-world SCC program.Kumaret al.[90] used the AFS algorithm as a stochastic-based search algorithm for serviceably solving the optimal scheduling problem of energy generation among available renewable energy sources.Bamini and Enoch [91] applied a classification and regression tree as a workflow to increase an allocation of resources,and a modified BFO algorithm was implemented to efficiently schedule the resources among the submitted jobs and to minimize the overall time of execution.

    B.Robot System

    Robot research in swarm intelligence focuses on the robot design,the physical body group size,and the control behaviors.Generally,this research studies the coordination of multiple robots,as a system consisting of many simple physical robots [92].It is supposed that the constant feedback for the desired collective behaviors appears from the interactions between these robots,and their interactions with the environment.As compared to ordinary distributed robotic systems,swarm robotics emphasizes large numbers of robots and heightens scalability.

    Globally optimal path planning is a common problem in the navigation of autonomous mobile robots,and the navigation must be collision-free.Tanet al.[93] utilized the ACO algorithm to resolve this problem,and verified that it provided better performance in convergence speed,solution variation,dynamic convergence behavior,and calculative efficiency than the general path planning technique,on account of a realcoded genetic algorithm with an elitist model.Contreras-Cruzet al.[94] combined the ABC algorithm and an evolutionary programming algorithm to improve a feasible path generated by a series of local optimization procedures.Experiments demonstrated the statistical significance of the promotion obtained by their proposed method.Yaoet al.[95] structured a hybrid adaptive AFS algorithm based on adaptive enhanced prey behavior and the segmented adaptive strategy for an artificial fish’s view and steps,and the algorithm was verified based on path planning for a coal mine rescue robot.In addition,for other issues in robotics research,Pughet al.[96]explored the performance of PSO in noisy environments while paying attention to unsupervised robotic learning and suggested an augmentation of the original method for achieving maximization of performance.Aghajarianet al.[97] presented an innovative application of the BFO algorithm to generate a fuzzy controller for tracking control of a robot manipulator.It was used efficiently to constitute a rule base and membership functions.The simulation results represented the superiority of the BFO algorithm (as compared to the PSO algorithm) for this application.

    C.Power System

    Owing to the network connections of the power system and further innovations in the electricity market,the power system is in fact developing to be a large-scale nonlinear dynamic system,and several complex engineering calculations are required to be resolved [98].Note that swarm intelligence does not require any preconditions for centralized control and a global model,it is really appropriate for solving large-scale nonlinear optimization problems that are difficult to solve using traditional methods (and for which it is hard to establish effective formalized models).Recently,there have been several studies on the successful practical application of swarm intelligence in power systems.

    Economic dispatch is of great value for electrical power system operations.To solve this complex problem,Zhouet al.[99] described a multi-objective multi-population ACO algorithm for a continuous domain.It proved valuable in minimizing a total fuel cost subject to various physical and operational constraints,via appropriately allocating power demands to generator modules.Abido [100] employed a PSO algorithm to find the optimal settings of power-system stabilizer parameters.An eigenvalue analysis and the nonlinear simulation results demonstrated its effectiveness in damping out the local and inter-area modes of oscillations,and it performed availably over a wide range of loading conditions and system configurations.Gozdeet al.[101]provided the usage of the ABC algorithm as an AI-based optimization technique for optimizing an automatic generation control system,and provided a comprehensive analysis of its tuning performance and contributions to robustness.The market clearing price (MCP) acts as one of the main factors for load control and market monitoring in a power system.Li and Wang [102] applied the AFS algorithm to improve the ability to determine the MCP,and it provided advantages insofar as fast training,the overcoming of local extrema,swift access to a global optimum,and so on.Ali and Abd-Elazim[103] proposed BFO-based load frequency control for the suppression of oscillations in a power system.It was adopted to obtain optimal controller parameters by minimizing a time domain objective function.

    D.Parameter Optimization

    Parameter optimization considers a field of problems that seek optimal control parameters [104].Traditionally,the most elaborate systematic methods for parameter optimization are on account of generalization error estimations and gradient descents.In recent years,the development of parameter optimizations for various operating systems has been supported by AI techniques and evolutionary strategies.

    One of the most important research problems in an SVM is the generation of optimal parameters to develop an efficient SVM,so as to achieve the desired output with an acceptable level of accuracy.Zhanget al.[105] applied the ACO algorithm to establish a new SVM model to resolve this problem.It only explored a finite subset of the possible values,to obtain the parameters that minimized the generalization error.A flux-cored arc welding (FCAW)procedure is a fusion welding procedure in which the welding electrode is a tubular wire that is continuously fed to the weld area.The welding input parameters are truly critical in determining the quality of a weld joint.Pandaet al.[106]provided a model of the weld bead geometry in the FCAW course based on an ANN,and optimization of the process parameters utilizing the PSO algorithm.Focusing on the improvement of algorithm performance,Akay and Karaboga[107] introduced some improvements on the original ABC algorithm,and the performance of the improved ABC algorithm was investigated for real-parameter optimization in both basic and composite functions.Proportional integral derivative (PID) control,as a classical control strategy,has good robustness,and is widely adopted in process and motion control.Cheng and Hong [108] utilized the AFS algorithm for parameter optimization,which is the core of PID controller design.In addition,the decrease of product development cycle time is a core concern for industries intended to keep competitive in the marketplace;therefore,rapid fabrication techniques are increasingly important.Pandaet al.[109]applied the BFO algorithm to represent a theoretical combination of parameter settings for simultaneously achieving great strength in all procedure responses.

    E.System Identification

    System identification is the phrase adopted in the automatic control field to denote the optimal estimation of dynamic black-box models for systems,in view of measurements of the input and output signals of those systems [110].It aims to start from measurements of the system behavior and external influences,and attempts to confirm a mathematical relation between them,without investigating the details of what is actually occurring inside the system.

    Fuzzy logic has been used in the analysis of physical and engineering problems.In a fuzzy identification problem,several unknown parameters are required to be found,and the interval for each parameter is uncertain.Thus,Tsai and Chen[111] provided a novel ACO-based fuzzy identification algorithm for finding silted initial consequent parameters and an initial membership function matrix.In [112],Alfi and Modares proposed a method to find optimal system parameters and optimal control parameters utilizing a new adaptive PSO algorithm.Its main advantages were in achieving a faster convergence speed and better solution accuracy with the minimum increment in the computational burden.To distinguish unknown fractional-order chaotic systems from the view of optimization,Huet al.[113]presented an innovative parameter estimation scheme on account of a hybrid ABC algorithm;it was proven to be a feasible,effective,and promising method.Hanet al.[114]addressed parameter identification for a photovoltaic module in view of an enhanced AFS algorithm.Its feasibility was approved by identifying the parameters of one commercial photovoltaic submodule with a single diode model under various operating conditions.Majhi and Panda [115]developed a valuable identification scheme for nonlinear dynamic systems adopting combinatorial swarm intelligence algorithms,e.g.,BFO and PSO.Compared to functional link ANN-based methods,the combinatorial technique is faster,more accurate,and involves less computation.

    F.Image Processing

    Image processing is the use of a computer to process different images using algorithms [116].The main influencing factor in the development of image processing is the increasing demand for a comprehensive range of applications,e.g.,in medicine,the military,agriculture,and industry.Currently,increasing research has been conducted on natureinspired methods for image processing,such as swarm intelligence algorithms.

    An edge is one of the simplest and the most important features of an image,and it is widely adopted in image recognition,segmentation,enhancement,and compression.Tianet al.[117] successfully applied an ACO algorithm to image edge detection by combining gradients and the relative differences of statistical means.Based on PSO,Omranet al.[118] provided a novel dynamic clustering method for image segmentation.It could automatically determine the optimum cluster number and simultaneously cluster the data set,with minimal user interference.Image enhancement acts as the core phase in an image processing system and tries to improve both the visual and informational qualities of distorted images.Draa and Bouaziz [119] proposed an innovative ABC algorithm for image contrast enhancement,and its superiority was proven by comparing the obtained results with those of a genetic algorithm.For efficient color image quantization,El-Said [120] exploited the optimization capability of an improved AFS algorithm to overcome the deficiencies of a fuzzy C-means algorithm.The algorithm could also be used to obtain visually improved images after the quantization.In addition,to enhance color images using a fuzzy logic technique,Hanmandluet al.[121] used a bacterial foraging algorithm to generate the parameters for optimized entropy and the visual factors which are involved in the objective function.

    G.Signal Processing

    Signal processing plays an extremely key role in extracting valuable information from different detectors,and acts as an electrical engineering subfield which includes analyzing,modifying,and synthesizing signals such as sounds,images,and biological measurements [122].Relevant techniques are applied to promote the transmission,storage efficiency,and subjective quality of the signals,and also to emphasize or find components of interest in a measured signal.

    Signal processing in regard to brain activity has become a challenge for various researchers.Huanget al.[123] proposed an ACO algorithm to generate the distributions of the intensity of a targeted image,so as to obtain the spatiotemporal characteristics of brain activities.Shadmand and Mashoufi[124] used an ANN to classify electrocardiogram heartbeats of a patient into five typical types,and the network structure and weights were effectively optimized adopting a classic PSO algorithm.Koza and Karaboga [125] utilized an ABC algorithm for the construction of a quadrature mirror filter bank with a mitral valve Doppler signal,and used it in various subjects,such as image coding and signal processing.Jiang and Yuan [126] combined an AFS algorithm with a conventional wavelet threshold algorithm,and enhanced the effects of signal processing using optimization.Sahuet al.[127] resolved the problem of electric power systems being polluted by unwanted fluctuations in both the current and voltage signals,based on an improved BFO algorithm.The algorithm could also raise the classification accuracy of various power quality events.

    IV.DISCUSSIONS AND TRENDS

    Recently,swarm intelligence has been extensively discussed in the field of multi-objective optimization.Consequently,many related approaches have emerged.Along with the promising performances swarm intelligence has realized,the research literature has also indicated some important challenges,as well as inherent trends.These are summarized in Table I,and are described below.

    TABLEI ANALYSIS OF SWARM INTELLIGENCE ALGORITHMS

    A.Strengths

    Swarm intelligence algorithms present advantages as compared to traditional methods,and even as compared to other AI algorithms.Based on the relevant literature,some representative advantages are summarized as follows.

    1) The most basic feature of the swarm intelligence algorithms is the flexible number of individuals,facilitating higher scalability.Thus,their control mechanisms are not excessively based on swarm size,as long as it is not too small or too large.

    2) They present the capability to realise a relatively largescale search,and to refine the solutions in a single search process.In the phase of generating the optimal solution,they usually present excellent exploration and exploitation capabilities.

    3) Several simple individuals interacting locally among themselves can eventually lead to a sophisticated global behavior.Low-level interactions between simple individuals have been represented to make the percolation of information known by only a few individuals throughout an entire aggregate,and to promote the emergence of collective intelligence.This mechanism is particularly suitable for some frontiers in special fields,such as targeted medical therapies,bee colony drones,and geographic surveillance.

    4) Basically,the swarm operates without central control,and no specific individual is indispensable for the swarm to keep operation.This means that a malfunction in any part of this system would not increase the risk of a complete failure of the entire swarm system,providing good robustness.

    B.Limitations

    As concluded in the majority of the papers,swarm intelligence algorithms have proven to be valuable in dealing with various tasks and challenges.However,there are also several limitations worth mentioning,and the representative limitations are summarized as follows.

    1) Swarm intelligence algorithms are generally timeconsuming processes which are affected by factors such as the size of population,frequency of iterations,and pattern of iterations.Thus,these factors directly affect the efficiency of swarm intelligence algorithms relative to the size of relevant applications,and these algorithms may be worthless when the factors exceed certain sizes.

    2) Basic swarm intelligence algorithms can suffer from stagnation or a premature convergence to a local optimum,owing to the lack of central coordination.Therefore,they need to be improved to realise adaptive mechanisms for continuously exploiting and exploring the search space,and for balancing the searching speed.

    3) Blocking coordination mechanisms increase the calculation time.In particular,an agent has to stay for other agents to be evaluated before it could continue to another position and start exploring the search domain.For example,the basic PSO algorithm is computationally costly,owing to the waiting time between evaluation of the first particle and that of the rest of the swarm before the state could finally be updated.

    C.Future Research Directions

    The research in swarm intelligence algorithms has shown a positive evolution;in addition,it provides several future research directions that may contribute to further growth in the field of swarm intelligence,as discussed below.

    1) Novel extensions of swarm algorithms.The existing swarm intelligence algorithms form an older category of optimization techniques,i.e.,immense research has been introduced,and potential application areas have been expanded.These techniques are represented by an adequate number of references in the previous research.Deeper investigation into the normal,possibly universal,features of the collective behaviors of animals,bacteria,cells,molecular motors,and driven granular matter are necessary for extension of the swarm algorithm family.Appropriate and novel application areas are required to be sought out for other advanced and/or simplified,time-effective algorithms,as these are still unexplored and comparatively newer approaches to resolving complicated computation problems.

    2) Characteristic refinement of individual agent.Future research should also pay attention to the characteristic refinement of each individual agent,to promote the overall optimization performance of the swarm intelligence algorithms.For example,there are two proverbial kinds of interaction patterns:hierarchical,and egalitarian.Some research has indicated that an individual will switch interaction patterns in different situations.This switching mechanism is prospective for potential industrial applications in multi-robot system coordination,unmanned vehicle formation control,and other areas.Within these individual agent-based techniques,there is also immense potential for fine-tuning,i.e.,to further optimize outcomes based on context-specific requirements.Furthermore,the interplay of social and environmental information and their influences on population performance could be further investigated.

    3) Further improvement of application transformation.In nature,thousands (or even millions) of individuals could easily form a great range of patterns such as bacteria colonies,bird flocks,and insect swarms,purely based on local communications.In view of the designed methods and wireless communication,robotic systems could simulate several complex biological swarm structures.Using patterns and rules,these elements could dramatically amend the swarm shape based on the environment.However,creating a swarming robotic system with a large number of individual members at a micro-scale that embodies functional collective behaviors is also a challenge.Therefore,while expanding the algorithm,future research should also focus on hardware developments and algorithm transformations in the application area.

    4) Comprehensive consideration of negative factors.It is well-known that social interactions benefit swarm members in diverse ways,such as by increasing the foraging ability or reducing the predation risk.However,the social benefits should be weighed against possible adverse effects,such as de-valued individual information,and reduced sensitivity to changing environments.Individuals must weigh environmental information gathered from their senses along with social information or face a larger risk of isolation and predation.These possible negative factors from group members should be investigated,to make it easier to consider and simulate negative downside elements in an actual system when applying swarm intelligence algorithms.

    5) Normative systematization of refined datasets.Swarm intelligence is achieving popularity for resolving different optimization problems and is attracting the attention of an incremental number of researchers.Researchers should pay attention to comparing different swarm intelligence algorithms and their variants and providing in-depth empirical and analytic studies.All the innovative algorithms should be tested,and a normative systematization of refined datasets should be required and/or recommended.The refined datasets should be classified based on the different algorithms.In addition,a taxonomy of uniformly recognized evaluation methods (e.g.,for testing the performance,advantages,and drawbacks of specific algorithms) are necessary for improving comparisons of the results,and could greatly promote the application process.

    V.CONCLUSION

    Swarm intelligence algorithms mimic the swarm intelligence behavior of biological elements in nature,and have been gradually popular cross-discipline and/or research field in recent years.This review provides a comprehensive survey of swarm intelligence and represents a categorization scheme for analyzing the existing swarm intelligence studies with applications.It is evident from the study that the representative algorithms discussed here (ACO,PSO,AFS and BFO and ABC) are very popular and well-explored,with several literatures relating to each subject.Authors have also expressed details of lesser-known algorithms,among which the cuckoo search,pigeon-inspired optimization,bat,grey wolf optimizer,artificial immune system,and fruit fly optimization algorithms have formed an important and growing category.These algorithms have been utilized in numerous diverse domains for solving optimization problems,e.g.,scheduling problems,robots,power systems,parameter optimization,system identification,image processing,and signal processing.Despite the promising results reported so far,there is significant space for further advances.The strengths and limitations are discussed and future research trends are summarized,in an attempt to provide a standard for the evaluation and application of the algorithms and to display potential valuable directions for subsequent research.Swarm intelligence is still active,and its potential is far from being exhausted;research and applications will continue to increase in the future.

    ABBREVIATIONS

    In Table II,we summarize the abbreviations and notations used throughout the review in alphabetical order.

    TABLE IIABBREVIATIONS USED IN THIS OVERVIEW

    国产欧美日韩综合在线一区二区| 久久久久久久久久久久大奶| 精品一区二区三卡| 欧美午夜高清在线| 亚洲国产精品sss在线观看 | 中出人妻视频一区二区| 亚洲全国av大片| 老司机亚洲免费影院| 亚洲情色 制服丝袜| 久久这里只有精品19| 国产伦人伦偷精品视频| 在线天堂中文资源库| 天天影视国产精品| 女人高潮潮喷娇喘18禁视频| 天天躁夜夜躁狠狠躁躁| 美女福利国产在线| 免费观看精品视频网站| 最新美女视频免费是黄的| 精品久久久久久电影网| 午夜免费成人在线视频| 在线免费观看的www视频| 在线观看www视频免费| 首页视频小说图片口味搜索| 男女做爰动态图高潮gif福利片 | 欧美日韩黄片免| 亚洲自偷自拍图片 自拍| 成人18禁高潮啪啪吃奶动态图| 亚洲一区高清亚洲精品| 亚洲免费av在线视频| 免费看十八禁软件| 欧美成人性av电影在线观看| 日本一区二区免费在线视频| 中文字幕最新亚洲高清| 国产精品乱码一区二三区的特点 | 国产亚洲欧美98| 一区二区三区精品91| 又大又爽又粗| 日韩大码丰满熟妇| 最近最新免费中文字幕在线| 美女 人体艺术 gogo| 欧洲精品卡2卡3卡4卡5卡区| av电影中文网址| 一区二区三区激情视频| 国产又爽黄色视频| 美女高潮到喷水免费观看| 一级毛片高清免费大全| 午夜成年电影在线免费观看| 午夜福利影视在线免费观看| 欧美一区二区精品小视频在线| 精品一区二区三卡| 大陆偷拍与自拍| tocl精华| 国产片内射在线| 精品卡一卡二卡四卡免费| 又紧又爽又黄一区二区| 精品久久久久久,| 精品乱码久久久久久99久播| 黑丝袜美女国产一区| 五月开心婷婷网| 啪啪无遮挡十八禁网站| 日本一区二区免费在线视频| 亚洲欧美日韩高清在线视频| 欧美精品亚洲一区二区| avwww免费| 99在线人妻在线中文字幕| 一夜夜www| 后天国语完整版免费观看| 成人三级做爰电影| 国产极品粉嫩免费观看在线| 日本三级黄在线观看| 久久国产精品人妻蜜桃| 国产成人精品久久二区二区免费| 国产激情欧美一区二区| 亚洲av五月六月丁香网| 色在线成人网| 国产成人精品在线电影| tocl精华| 国产高清videossex| 韩国av一区二区三区四区| 一级a爱视频在线免费观看| 久久婷婷成人综合色麻豆| 精品久久久久久,| 黄色视频不卡| 色精品久久人妻99蜜桃| 精品久久久久久电影网| 免费不卡黄色视频| 亚洲第一欧美日韩一区二区三区| 高潮久久久久久久久久久不卡| 亚洲精品成人av观看孕妇| 国产男靠女视频免费网站| avwww免费| 99精国产麻豆久久婷婷| 国产成人av教育| 级片在线观看| 中文字幕精品免费在线观看视频| 亚洲全国av大片| 日韩欧美三级三区| 色综合站精品国产| 国产av又大| www.精华液| 亚洲欧美一区二区三区久久| 动漫黄色视频在线观看| 国产成人系列免费观看| 91av网站免费观看| 日本撒尿小便嘘嘘汇集6| 国产人伦9x9x在线观看| 亚洲va日本ⅴa欧美va伊人久久| 乱人伦中国视频| 成年人免费黄色播放视频| 电影成人av| 操美女的视频在线观看| 国产高清国产精品国产三级| 中文欧美无线码| 很黄的视频免费| 午夜福利免费观看在线| 亚洲,欧美精品.| 精品人妻在线不人妻| 老熟妇仑乱视频hdxx| 亚洲精品在线美女| 黑人巨大精品欧美一区二区蜜桃| 国产xxxxx性猛交| 国产又爽黄色视频| 欧美黑人欧美精品刺激| 亚洲狠狠婷婷综合久久图片| 啦啦啦 在线观看视频| 成人三级黄色视频| 999久久久精品免费观看国产| 精品久久久精品久久久| 欧美国产精品va在线观看不卡| 亚洲欧美精品综合久久99| 日本vs欧美在线观看视频| 可以免费在线观看a视频的电影网站| 老鸭窝网址在线观看| 99国产精品一区二区蜜桃av| 亚洲欧美一区二区三区黑人| 黄频高清免费视频| 免费看十八禁软件| 两性夫妻黄色片| 大香蕉久久成人网| 久久精品亚洲av国产电影网| 99久久国产精品久久久| 久久草成人影院| 久久午夜综合久久蜜桃| 日日摸夜夜添夜夜添小说| 老熟妇乱子伦视频在线观看| 国产xxxxx性猛交| 国产精品二区激情视频| 1024视频免费在线观看| 精品国产一区二区久久| 欧美最黄视频在线播放免费 | 亚洲va日本ⅴa欧美va伊人久久| 乱人伦中国视频| 国产一区二区激情短视频| 男人操女人黄网站| 国产精品日韩av在线免费观看 | 女性被躁到高潮视频| 99国产极品粉嫩在线观看| 国产欧美日韩精品亚洲av| 午夜福利在线观看吧| 国产精品久久电影中文字幕| 看黄色毛片网站| 12—13女人毛片做爰片一| 日韩欧美三级三区| 嫁个100分男人电影在线观看| 成人18禁在线播放| 在线观看www视频免费| 久久中文看片网| 免费搜索国产男女视频| 久久伊人香网站| av电影中文网址| 久久久久久免费高清国产稀缺| 欧美日韩黄片免| 亚洲精品一二三| 最近最新免费中文字幕在线| 黄片播放在线免费| 少妇粗大呻吟视频| 欧美午夜高清在线| 99热国产这里只有精品6| 女人被狂操c到高潮| 久久 成人 亚洲| 亚洲精品在线美女| 成人av一区二区三区在线看| 免费看a级黄色片| 久久中文看片网| 高清毛片免费观看视频网站 | 无限看片的www在线观看| 婷婷精品国产亚洲av在线| videosex国产| 99久久精品国产亚洲精品| 久久青草综合色| 妹子高潮喷水视频| 久久亚洲精品不卡| 久久中文字幕一级| 97超级碰碰碰精品色视频在线观看| 国产片内射在线| 日日摸夜夜添夜夜添小说| 亚洲欧美日韩另类电影网站| 亚洲男人天堂网一区| 两性夫妻黄色片| 国产麻豆69| av国产精品久久久久影院| av电影中文网址| 99久久人妻综合| 女性生殖器流出的白浆| 久久这里只有精品19| 欧洲精品卡2卡3卡4卡5卡区| 老司机在亚洲福利影院| 国产主播在线观看一区二区| 免费看十八禁软件| 欧美久久黑人一区二区| 神马国产精品三级电影在线观看 | 欧美中文综合在线视频| 精品少妇一区二区三区视频日本电影| 亚洲精品一区av在线观看| 国产精品九九99| 欧美日韩一级在线毛片| 午夜福利影视在线免费观看| 国产av一区二区精品久久| av超薄肉色丝袜交足视频| 搡老乐熟女国产| 97人妻天天添夜夜摸| 国产亚洲精品第一综合不卡| 男女下面插进去视频免费观看| www国产在线视频色| 亚洲五月色婷婷综合| 午夜福利,免费看| 成人黄色视频免费在线看| 国产成人av激情在线播放| 成人特级黄色片久久久久久久| 精品国产美女av久久久久小说| 在线观看午夜福利视频| 国产av一区在线观看免费| 欧美最黄视频在线播放免费 | 国产一区二区三区视频了| 亚洲视频免费观看视频| 男人的好看免费观看在线视频 | 国产日韩一区二区三区精品不卡| 9191精品国产免费久久| 激情在线观看视频在线高清| 国产高清视频在线播放一区| 国产亚洲欧美在线一区二区| 亚洲人成伊人成综合网2020| 亚洲av五月六月丁香网| 欧美激情 高清一区二区三区| 淫秽高清视频在线观看| 最新在线观看一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 黄片小视频在线播放| 99精品久久久久人妻精品| 精品卡一卡二卡四卡免费| 精品人妻1区二区| 一夜夜www| 男女高潮啪啪啪动态图| 精品福利观看| 久久国产乱子伦精品免费另类| 国产色视频综合| 成人黄色视频免费在线看| 亚洲欧美精品综合久久99| 淫妇啪啪啪对白视频| 免费日韩欧美在线观看| 三上悠亚av全集在线观看| 国产深夜福利视频在线观看| 午夜激情av网站| 午夜福利欧美成人| 国产有黄有色有爽视频| 夜夜夜夜夜久久久久| 精品久久久久久成人av| 国产成人精品久久二区二区免费| 国产成人免费无遮挡视频| 亚洲午夜精品一区,二区,三区| 乱人伦中国视频| 精品人妻在线不人妻| 香蕉丝袜av| 欧美日韩福利视频一区二区| netflix在线观看网站| 国产真人三级小视频在线观看| 又大又爽又粗| 一级毛片高清免费大全| 免费日韩欧美在线观看| 欧美一级毛片孕妇| 久9热在线精品视频| 夫妻午夜视频| 中出人妻视频一区二区| 韩国av一区二区三区四区| 伦理电影免费视频| 国产精品电影一区二区三区| 国产精品免费视频内射| 女人高潮潮喷娇喘18禁视频| 黄色视频不卡| 男人操女人黄网站| 国产精品av久久久久免费| 午夜精品在线福利| 美女扒开内裤让男人捅视频| 亚洲国产精品合色在线| www.自偷自拍.com| 天堂动漫精品| 日韩一卡2卡3卡4卡2021年| 曰老女人黄片| 19禁男女啪啪无遮挡网站| 十分钟在线观看高清视频www| 一二三四在线观看免费中文在| 中文字幕人妻丝袜一区二区| 久久这里只有精品19| 我的亚洲天堂| 亚洲少妇的诱惑av| 男女下面插进去视频免费观看| 亚洲欧美日韩高清在线视频| 欧美av亚洲av综合av国产av| 欧美精品一区二区免费开放| 国产精品乱码一区二三区的特点 | 欧美精品亚洲一区二区| 成年人免费黄色播放视频| 亚洲第一欧美日韩一区二区三区| 精品久久久久久电影网| 精品少妇一区二区三区视频日本电影| 国产免费现黄频在线看| 精品久久久久久成人av| 两性午夜刺激爽爽歪歪视频在线观看 | 婷婷丁香在线五月| 人妻丰满熟妇av一区二区三区| 亚洲欧美激情在线| 在线观看免费视频日本深夜| 99re在线观看精品视频| 99riav亚洲国产免费| 精品一区二区三区av网在线观看| 一级毛片高清免费大全| 12—13女人毛片做爰片一| 中文字幕人妻丝袜一区二区| 久久精品亚洲av国产电影网| 国产在线精品亚洲第一网站| 国产成+人综合+亚洲专区| 亚洲七黄色美女视频| 无人区码免费观看不卡| 国产免费男女视频| 天天影视国产精品| 国产av精品麻豆| 久久久国产成人免费| www.熟女人妻精品国产| 无限看片的www在线观看| 美女国产高潮福利片在线看| 搡老岳熟女国产| 叶爱在线成人免费视频播放| av福利片在线| 黄色毛片三级朝国网站| 男女午夜视频在线观看| 侵犯人妻中文字幕一二三四区| 亚洲精品在线观看二区| 男女下面进入的视频免费午夜 | 日本一区二区免费在线视频| 一区二区三区国产精品乱码| 动漫黄色视频在线观看| 99热只有精品国产| 一二三四社区在线视频社区8| 国产精品久久视频播放| 亚洲一码二码三码区别大吗| av在线天堂中文字幕 | 69av精品久久久久久| av天堂久久9| 777久久人妻少妇嫩草av网站| 国产成人系列免费观看| 久久久久九九精品影院| 韩国av一区二区三区四区| 午夜福利,免费看| 两人在一起打扑克的视频| 黄网站色视频无遮挡免费观看| 国产伦人伦偷精品视频| 91精品国产国语对白视频| 视频区欧美日本亚洲| 国产av一区在线观看免费| 午夜影院日韩av| 国产欧美日韩一区二区精品| 国产精品免费视频内射| 免费女性裸体啪啪无遮挡网站| 精品一区二区三区av网在线观看| 久久精品aⅴ一区二区三区四区| 日本一区二区免费在线视频| 中文字幕精品免费在线观看视频| 天堂动漫精品| 高清av免费在线| 国产欧美日韩一区二区三区在线| 1024视频免费在线观看| 高清欧美精品videossex| 欧美激情 高清一区二区三区| 亚洲欧美激情综合另类| 亚洲美女黄片视频| 国产精品美女特级片免费视频播放器 | 国产一区二区激情短视频| 精品第一国产精品| 亚洲人成77777在线视频| 亚洲男人的天堂狠狠| 亚洲成人精品中文字幕电影 | 国产亚洲精品久久久久久毛片| 高清毛片免费观看视频网站 | 日韩人妻精品一区2区三区| 黄色成人免费大全| 久久九九热精品免费| 亚洲精品成人av观看孕妇| 午夜免费激情av| 一区二区三区激情视频| 久99久视频精品免费| 欧美日韩黄片免| 国产精品久久视频播放| 久久人人爽av亚洲精品天堂| 欧美久久黑人一区二区| 国产真人三级小视频在线观看| 18禁黄网站禁片午夜丰满| 国产乱人伦免费视频| 啦啦啦 在线观看视频| 国产精品 欧美亚洲| ponron亚洲| 亚洲性夜色夜夜综合| 午夜免费鲁丝| 80岁老熟妇乱子伦牲交| 97超级碰碰碰精品色视频在线观看| 亚洲精品国产区一区二| 亚洲精品中文字幕一二三四区| 国产成人免费无遮挡视频| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| 亚洲全国av大片| 国产精品偷伦视频观看了| 欧美不卡视频在线免费观看 | 欧美在线一区亚洲| netflix在线观看网站| 长腿黑丝高跟| 国产人伦9x9x在线观看| 天堂动漫精品| 黑人操中国人逼视频| 免费观看人在逋| 日韩精品中文字幕看吧| 又大又爽又粗| 18禁观看日本| 中出人妻视频一区二区| 91麻豆精品激情在线观看国产 | 黑丝袜美女国产一区| 久久精品亚洲熟妇少妇任你| a级毛片在线看网站| 激情视频va一区二区三区| 操出白浆在线播放| 村上凉子中文字幕在线| 在线播放国产精品三级| 国产欧美日韩精品亚洲av| 在线国产一区二区在线| 男女下面插进去视频免费观看| 黑人巨大精品欧美一区二区mp4| 后天国语完整版免费观看| 97碰自拍视频| 国产亚洲欧美精品永久| 免费在线观看完整版高清| 欧美最黄视频在线播放免费 | 免费在线观看日本一区| 久久久国产精品麻豆| 成人18禁在线播放| 精品国产美女av久久久久小说| 亚洲精品久久午夜乱码| 国产有黄有色有爽视频| 女同久久另类99精品国产91| 操出白浆在线播放| 久久国产精品男人的天堂亚洲| 高清在线国产一区| 国产一区二区激情短视频| 成人永久免费在线观看视频| 视频在线观看一区二区三区| 999久久久国产精品视频| 女人被狂操c到高潮| 日韩精品免费视频一区二区三区| 午夜免费鲁丝| 中文字幕人妻丝袜制服| 亚洲成人免费电影在线观看| 精品久久久久久电影网| 男女做爰动态图高潮gif福利片 | 精品一区二区三区av网在线观看| 久久久国产精品麻豆| 欧美精品啪啪一区二区三区| 成年人黄色毛片网站| 久久欧美精品欧美久久欧美| 午夜久久久在线观看| 亚洲av成人av| 一本综合久久免费| 亚洲人成电影观看| 99久久人妻综合| 亚洲男人的天堂狠狠| 亚洲熟妇中文字幕五十中出 | 又大又爽又粗| 神马国产精品三级电影在线观看 | 久久久水蜜桃国产精品网| 亚洲一区二区三区不卡视频| 黑丝袜美女国产一区| 手机成人av网站| 成人手机av| 久久人人97超碰香蕉20202| 亚洲欧美精品综合久久99| 精品乱码久久久久久99久播| 咕卡用的链子| 国产成人影院久久av| 天天躁夜夜躁狠狠躁躁| 亚洲avbb在线观看| 日韩三级视频一区二区三区| 中文字幕精品免费在线观看视频| 亚洲五月婷婷丁香| 午夜影院日韩av| 天堂影院成人在线观看| 视频区欧美日本亚洲| 亚洲午夜精品一区,二区,三区| 国产av又大| 亚洲国产毛片av蜜桃av| 精品一区二区三区av网在线观看| 国产亚洲精品久久久久5区| 一级黄色大片毛片| x7x7x7水蜜桃| 久久天堂一区二区三区四区| 成人精品一区二区免费| 女生性感内裤真人,穿戴方法视频| 亚洲av片天天在线观看| 男女下面插进去视频免费观看| 亚洲精品中文字幕在线视频| 欧美黄色淫秽网站| 91麻豆精品激情在线观看国产 | av超薄肉色丝袜交足视频| 午夜精品国产一区二区电影| 精品日产1卡2卡| 欧美色视频一区免费| 亚洲精品在线观看二区| 亚洲欧洲精品一区二区精品久久久| 欧美日韩中文字幕国产精品一区二区三区 | 一进一出抽搐动态| 自线自在国产av| 中出人妻视频一区二区| 老汉色av国产亚洲站长工具| 免费在线观看视频国产中文字幕亚洲| 黄片大片在线免费观看| 涩涩av久久男人的天堂| 自拍欧美九色日韩亚洲蝌蚪91| av在线播放免费不卡| 日韩大尺度精品在线看网址 | 亚洲欧美精品综合久久99| 亚洲久久久国产精品| 日韩国内少妇激情av| 成人黄色视频免费在线看| 国产又爽黄色视频| 亚洲精品av麻豆狂野| 最近最新中文字幕大全免费视频| 少妇粗大呻吟视频| 亚洲国产中文字幕在线视频| 久久精品人人爽人人爽视色| 少妇 在线观看| 国产亚洲欧美98| 1024香蕉在线观看| 国产成+人综合+亚洲专区| 国产色视频综合| 国产又爽黄色视频| 黄色毛片三级朝国网站| 午夜福利在线观看吧| 亚洲国产欧美一区二区综合| 自线自在国产av| 国产高清videossex| 亚洲av电影在线进入| 亚洲人成网站在线播放欧美日韩| 日本vs欧美在线观看视频| 亚洲,欧美精品.| 精品人妻在线不人妻| 午夜日韩欧美国产| 亚洲一区二区三区欧美精品| 国产亚洲精品第一综合不卡| 法律面前人人平等表现在哪些方面| 宅男免费午夜| 丁香六月欧美| 超碰97精品在线观看| 婷婷六月久久综合丁香| 丰满人妻熟妇乱又伦精品不卡| 亚洲第一青青草原| 看片在线看免费视频| 91av网站免费观看| 国产精品99久久99久久久不卡| 后天国语完整版免费观看| 久久久久亚洲av毛片大全| 两人在一起打扑克的视频| 男女下面进入的视频免费午夜 | 国产亚洲精品久久久久久毛片| 国产精品偷伦视频观看了| 每晚都被弄得嗷嗷叫到高潮| 日本一区二区免费在线视频| 国产三级在线视频| 欧美日韩亚洲国产一区二区在线观看| 一进一出好大好爽视频| 免费在线观看日本一区| 大陆偷拍与自拍| 老司机在亚洲福利影院| 日韩高清综合在线| 一进一出抽搐gif免费好疼 | 99久久人妻综合| 老司机福利观看| 在线免费观看的www视频| 看片在线看免费视频| 黄频高清免费视频| 侵犯人妻中文字幕一二三四区| 久9热在线精品视频| 91成人精品电影| 久久伊人香网站| 免费在线观看完整版高清| 99国产精品99久久久久| 亚洲欧美一区二区三区黑人| 一边摸一边做爽爽视频免费| 欧美日韩福利视频一区二区| 免费看a级黄色片| 多毛熟女@视频| 天堂√8在线中文| 一本大道久久a久久精品| 后天国语完整版免费观看| 黄色女人牲交| 日韩国内少妇激情av| 国产成人一区二区三区免费视频网站| 88av欧美| av免费在线观看网站| 狂野欧美激情性xxxx|