[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN109144102B - Unmanned aerial vehicle route planning method based on improved bat algorithm - Google Patents

Unmanned aerial vehicle route planning method based on improved bat algorithm Download PDF

Info

Publication number
CN109144102B
CN109144102B CN201811091770.6A CN201811091770A CN109144102B CN 109144102 B CN109144102 B CN 109144102B CN 201811091770 A CN201811091770 A CN 201811091770A CN 109144102 B CN109144102 B CN 109144102B
Authority
CN
China
Prior art keywords
bat
aerial vehicle
unmanned aerial
algorithm
route
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811091770.6A
Other languages
Chinese (zh)
Other versions
CN109144102A (en
Inventor
林娜
唐嘉诚
赵亮
拱长青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang Aerospace University
Original Assignee
Shenyang Aerospace University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenyang Aerospace University filed Critical Shenyang Aerospace University
Priority to CN201811091770.6A priority Critical patent/CN109144102B/en
Publication of CN109144102A publication Critical patent/CN109144102A/en
Application granted granted Critical
Publication of CN109144102B publication Critical patent/CN109144102B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

An unmanned aerial vehicle route planning method based on an improved bat algorithm is based on a traditional bat algorithm, an optimization success rate is introduced to change a bat individual speed updating mode, meanwhile, a chaos method is adopted to initialize distribution of bat individuals in a search space, a gravitational field of a terminal point, a starting point and a repulsive force field of an obstacle are simulated by utilizing the concept of an artificial potential field, the speed of the bat individuals moving to an optimal solution is accelerated, and finally, the improved bat algorithm based on the chaos artificial potential field is provided. In the unmanned aerial vehicle route planning task, compared with the traditional bat algorithm, the method has the advantages that the route length is 36.56% shorter, the planning time is 56.05% shorter, and the obstacle avoidance effect adaptability value is 49.53% lower; compared with a differential evolution bat algorithm, the flight path length is 27.16 percent, the planning time is 27.30 percent, and the obstacle avoidance effect adaptability value is 42.46 percent. The method is an air route planning algorithm with practical application significance.

Description

Unmanned aerial vehicle route planning method based on improved bat algorithm
Technical Field
The invention belongs to the technical field of intelligent control, and relates to an unmanned aerial vehicle route planning method, in particular to an unmanned aerial vehicle route planning method based on an improved bat algorithm.
Background
The unmanned aerial vehicle route planning generally refers to a process of finding a flyable route which is from a starting point to a target point and meets the performance index of the unmanned aerial vehicle under a specific constraint condition. The algorithm adopted for the unmanned aerial vehicle route planning problem directly influences the success rate and efficiency of route planning. The swarm intelligence algorithm generally has the advantages of high convergence speed, strong robustness, potential parallelism and the like, and is widely applied to unmanned aerial vehicle route planning. However, when the group intelligent algorithms are used for solving the route planning problem, the defects that the solving precision is not high enough, the planned track is not smooth and the local optimization is easy to happen generally exist. Therefore, an intelligent algorithm with excellent optimization performance is selected and pertinently improved, so that the unmanned aerial vehicle routing problem is better.
The swarm intelligence algorithm is widely applied to solving the problem of unmanned aerial vehicle route planning by researchers at home and abroad, and a typical algorithm is selected and introduced as follows. Shibo Li and the like firstly use a particle swarm algorithm and combine fuzzy logic to solve the problem of two-dimensional unmanned aerial vehicle route planning, and carry out comparative analysis on different planning effects after parameter change in the particle swarm algorithm, but the defect that the two-dimensional experimental environment is too simple to set exists, and meanwhile, a three-dimensional simulation experiment does not exist; hao Meng and the like use a genetic algorithm to combine a digital elevation map and an RBF neural network to solve three-dimensional unmanned aerial vehicle route planning, provide a strict mathematical demonstration process, and obtain an intuitive three-dimensional route planning effect, but the method is limited by the genetic algorithm, and the convergence rate of the three-dimensional route planning and the route planning effect have a great space; the differential evolution algorithm is applied to the problem of multi-unmanned aerial vehicle collaborative route planning by Ioannis K.Nikolos and the like, a better solution is provided for the unmanned aerial vehicle formation to jointly complete tasks, but experimental results show that the route which is only planned by the differential evolution algorithm has the characteristic of Levy flight trajectory and cannot be directly used as the unmanned aerial vehicle flight route, in addition, the solving speed of the algorithm is not fast enough, and the execution stability of the algorithm is not high enough; GaigeWang et al apply bat algorithm to optimize flight cost function optimization in the unmanned aerial vehicle air route planning problem, carry on the contrastive analysis with other optimization effects of traditional group intelligence algorithm, prove that bat algorithm have solving fast, convergence precision high, algorithm stability advantage such as being good in optimizing the flight cost function, but lack the simulation experiment of the intuitive two-dimentional and three-dimensional air route planning effect; in subsequent researches, Gai-Ge Wang applies a bat algorithm based on a differential evolution algorithm to unmanned aerial vehicle route planning and provides two-dimensional and three-dimensional simulation experiment results, and compared and demonstrated the route planning effect of the bat algorithm combined with the differential evolution algorithm and the route planning effect of only using the bat algorithm, the superiority of the proposed algorithm is verified, but the experimental environment setting has the defect of over simplicity, which is not beneficial to the real optimizing and obstacle avoiding performance of the detection algorithm, and in addition, the planned route also has the characteristic of a Levy flight track, and route smoothing processing is not carried out.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides an unmanned aerial vehicle route planning method based on an improved bat algorithm, and the flow chart of the algorithm is shown in figure 1. The execution process of the route planning method comprises five main steps of improving a bat algorithm to carry out route planning, accelerating an iteration process by an artificial potential field method, homogenizing population distribution of a Logistic function of a chaotic strategy, establishing a route planning constraint condition and an objective function, and smoothing a route track. As will be explained in detail below.
Firstly, improving bat algorithm
The bat algorithm is a new meta-heuristic intelligent algorithm published by Yang.X.S in 2010, and the inspiration source of the bat algorithm is the behavior of bat in nature for capturing prey and avoiding obstacles and natural enemies through echo positioning; the bat algorithm follows the following three rules:
1. each bat individual in the search space is positioned only by relying on ultrasonic echo without involving other sensory factors such as vision and the like; in addition, the bat can judge the difference between food, game and obstacles according to the echo difference.
2. Flight optimization parameter setting of the bat individuals: the bat i flies at a speed vi, is located at a real-time position denoted by xi, and has an initial sound frequency fmin. The bat can dynamically adjust the frequency f and the volume A of the ultrasonic wave emitted by the bat according to the degree of approaching a target in the searching process0
3. Ideal restriction is made on the sound volume of bat sound production, namely the sound volume of bat sound production is assumed to be from a great positive number A0To the minimum value Amin without going out of bounds.
For the virtual bat in a d-dimensional search space, the bat sound production frequency f at the time tiVelocity viAnd position xiThe update formula of (2) is:
fi=fmin+(fmax-fmin)×β (1)
vi t=vi t-1+(xi t-1-x*)×f (2)
xi t=xi t-1+vi t (3)
for the update of the voicing frequency of equation (1), β is subject to a uniformly distributed random variable and satisfiesβ∈[0,1];fmaxAnd fminThe maximum value and the minimum value of the initial set sounding frequency range. X in formula (2)*Is the position of the current global optimal solution, and the value is the optimal value obtained by comparing fitness values of all the individuals in the bat population. The bat individual is positioned according to the position of the bat individual at the last moment
Figure GDA0003041186690000031
And global optimal solution x*Measure the acceleration of its movement to the optimal solution (urgency); the moving speed at the next moment also takes the inertia of the bat body into consideration, namely the moving speed at the next moment, besides the acceleration generated by referring to the near-global optimal solution
Figure GDA0003041186690000034
And also receives the speed of the last moment
Figure GDA0003041186690000032
The influence of (c). Formula (3) describes the process that the bat body is located and continuously shifts along with the movement.
The above is an iterative mechanism followed by the bat population when global search is performed in a solution space, and bats near the global optimal solution generate a local new solution by adopting a random walk rule:
xnew=xold+εAt (4)
in the formula (4), ε ∈ [ -1,1]Is a random number;
Figure GDA0003041186690000033
is the average volume of all bats at time t.
The update formula of the sound emission frequency and loudness at this time is described as follows:
Figure GDA0003041186690000041
ri t+1=ri 0[1-exp(-γt)] (6)
both α and γ in formula (5) and formula (6) are constants, and typically, α ═ γ ═ 0.9 is taken. From the expression, it can be analyzed that as the bat infinitely approaches the optimal solution, the loudness of sounding is continuously reduced and sounding is suspended when the optimal solution is reached; the sounding frequency is continuously close to the initial pulse rate r with the timei 0
Similar to exploration operation and mining operation in heuristic search algorithm, the swarm intelligence algorithm has global search and local search processes in the optimizing process. Global search is carried out to determine the approximate range of the optimal solution; while local searches are refined within these ranges.
The trade-off between global search and local search directly affects the search efficiency and the optimization accuracy. In order to better control the speed of global search of bat individuals, one of the innovation points of the invention is to introduce an adaptive inertia weight w to rewrite a bat speed updating formula (2), as shown in a formula (7):
vi t=wvi t-1+(xi t-1-x*)×f (7)
and the concept of the optimizing success rate is introduced, so that the inertia weight is adaptively adjusted along with the optimizing success rate of the bat group. The adaptive inertial weight based on the optimizing success rate is defined as follows:
Figure GDA0003041186690000042
wherein,
Figure GDA0003041186690000043
the success rate of the bat group optimization is achieved; n represents the bat population number;
Figure GDA0003041186690000044
shows the optimizing result of the bat individual i in the process of the t-th iteration,
Figure GDA0003041186690000045
the adaptive value of t generation is superior to t-1 generation, and a better solution is searched, if not, the optimal solution is obtained
Figure GDA0003041186690000051
Second, artificial potential field method
The artificial potential field method is originally published by Khatib in 1994, is applied to obstacle avoidance motion planning of a robot operating arm, and is then largely applied to path planning of a mobile robot; at present, no precedent exists that an artificial potential field method and an improved bat algorithm are fused and applied to the unmanned aerial vehicle route planning problem. The inspiration of the artificial potential field method is derived from the principle that electrostatic field heterogeneous charges generate attraction and homogeneous charges generate repulsion, and the acting force between an obstacle and an aircraft in a search space is defined as repulsion and the acting force between a target point and the aircraft. The stress diagram of the aircraft in the artificial potential field is shown in FIG. 3:
in the context of figure 3, it is shown,
Figure GDA0003041186690000052
which is indicative of the repulsive force,
Figure GDA0003041186690000053
which is indicative of the force of attraction,
Figure GDA0003041186690000054
representing the resultant force to which the aircraft is subjected, directly affecting the motion of the aircraft. According to the gradient descent method of the artificial potential field, the repulsive force and the attractive force can be expressed as formula (9) and formula (10):
Figure GDA0003041186690000055
Figure GDA0003041186690000056
let n be any point in the search space that is subject to an attractive force Fa(x) And repulsive force Fr(x) Can be expressed as:
Fa(x)=-grad(Ua(x))=kρG(q)(11)
Figure GDA0003041186690000057
Third, chaos strategy and Logistic mapping
The basis of the chaotic algorithm is a logistic mapping:
xn+1=μxn(1-xn-1)n=1,2,..... (13)
when μ is 4, the mapping is a perfect mapping over the interval [0, 1], the iteratively generated value is in a pseudo-randomly distributed state, and when other values are taken, the generated value converges to a specific value after a certain number of iterations. The bat populations are randomly distributed in the normalized search space, so that the overall search efficiency of the populations is greatly improved; meanwhile, the bat individuals are extremely difficult to have position transition in a large range in the respective local search process, and the uniform and random distribution of the initial population effectively solves the problem that the bat population is easy to fall into a local optimal solution in the later period of search.
Therefore, when an artificial potential field is combined, the method provides a second innovation point, namely, Logistic mapping based on a chaos strategy is adopted to initialize distribution of the bat populations in a solution space, so that the overall convergence speed of the algorithm is accelerated, and the possibility of falling into a local optimal solution is greatly reduced.
Constraint condition and objective function of route planning task
The definition of unmanned aerial vehicle route planning is a process of finding a route which has the minimum cost from a flight starting point to an end point and meets the performance index of the unmanned aerial vehicle under certain specific constraint conditions. The essence of the method is that under the condition of multiple constraints, the extremum is solved by a multi-objective function. The following analysis was made for each of the routings problems.
The cost function in the flight process of the unmanned aerial vehicle can be divided into three parts, namely, the route length cost, the threat cost and the fuel consumption cost. The total cost function is denoted by J, and the minimization problem with J is defined as follows:
min J=k1JL+k2JT+(1-k1-k2)JF (14)
wherein, JLIs the cost of the length of the flight path, JTIs the cost of the threat, JFRefers to the fuel consumption penalty incurred in flight. k is a radical of1,k2Are all normal numbers and satisfy k being more than or equal to 01≤1,0≤k2≤1。
For the length of flight path cost JLThe definition is as follows:
Figure GDA0003041186690000061
wherein L is the length of the total flight path, LijThe length of the track segment is used for calculating the length of the air route which cannot be calculated by integration before the air route is smoothly carried out.
Cost J for threatTThe definition is as follows:
Figure GDA0003041186690000071
wherein, tkThe threat factor is a measure of the threat degree of the threat source to the unmanned aerial vehicle; n is a radical oftRepresenting the total number of threat sources; the current coordinate of the unmanned aerial vehicle is (x, y); the center coordinate of the threat source is (x)k,yk)。
For fuel consumption cost JFThe definition is as follows:
Figure GDA0003041186690000072
wherein,
Figure GDA0003041186690000073
k in (1) indicates that the unmanned aerial vehicle travels along a unit length pathA fuel cost of consumption;
Figure GDA0003041186690000074
in H is the height of the unmanned aerial vehicle flight safety circle, the flight height of the unmanned aerial vehicle should not exceed the value, w0The fuel cost that the unmanned aerial vehicle needs to consume when maintaining the height H is shown, and H shows the current height of the unmanned aerial vehicle.
The unmanned aerial vehicle needs to accord with the dynamic characteristics of the unmanned aerial vehicle in the flying process, so that the preplanned air route needs to meet the limits of some constraint conditions. In order to more intuitively illustrate the constraints, a part of the waypoints in the flight path is selected as a model for illustration. Let A (x)i-1,yi-1,zi-1) As the last waypoint, B (x)i,yi,zi) As the current waypoint, C (x)i+1,yi+1,zi+1) For the next waypoint and records
Figure GDA0003041186690000075
And (4) a track point migration vector.
When the unmanned aerial vehicle wants to change the self height in the flying process, climbing or diving operation is needed, the invention assumes that the maximum climbing and diving angle is theta, and then the constraint on the climbing/diving angle of the unmanned aerial vehicle is as follows:
Figure GDA0003041186690000076
when the unmanned aerial vehicle avoids the obstacle, the situation that the turning radius corresponding to the planned track inscribed circle is too small may occur, and the flight characteristics of the unmanned aerial vehicle are violated. Therefore, a limit is made on the minimum turning radius of the unmanned aerial vehicle in flight, and the turning radius limit is schematically shown in fig. 4:
route angle
Figure GDA0003041186690000077
The minimum turning angle may be determined by:
Figure GDA0003041186690000081
wherein r isminIs the minimum turning radius.
The fly height limit depends on the specific characteristics of the mission being performed by the drone; in order to reduce consumption and ensure that the air route is hidden as much as possible, the absolute height H of the ground is limited to be the maximum value H which is not more than H, wherein H is the height of the plane where the flight safety circle is located; and in order to react to terrain changes in time, the relative height h of the unmanned aerial vehicle from the terrain surface is requiredi≥hmin
The unmanned aerial vehicle detects static obstacles and dynamic obstacles by utilizing a sensor of the unmanned aerial vehicle, and an avoidance mechanism of the obstacles after the obstacles are identified is a key factor influencing the route planning effect. As shown in fig. 5:
in order to meet the characteristics of the flight motion direction and the yaw angle of the unmanned aerial vehicle, the avoidance mechanism of the unmanned aerial vehicle on the obstacle is explained by using the motion trail under the polar coordinate system. In the figure, O is the track start point and G is the track end point, L1To LkThe track segment between every two adjacent segments is generated after equal segmentation based on the polar axis, and the starting point of each segmented track segment is identified by a square. The mechanism for measuring and avoiding threats between every two track points is visually shown in figure 6:
segmenting the sub-track section again according to the number of the threat sources to calculate the threat cost
Figure GDA0003041186690000082
In general, this threat cost may be expressed as:
Figure GDA0003041186690000083
wherein N isTRepresenting the number of threat sources, tjThreat factor representing a source of threat, dk(i, j) represents the linear distance between the starting points i and j of the kth track segment.
The obstacle avoidance process is the process of minimizing the threat cost under the condition of meeting the constraint condition, the advantage of the optimization problem is solved by using an intelligent algorithm, and the pre-planned optimal route can effectively avoid the static obstacle threat in the space. For avoiding the dynamic obstacle, the unmanned aerial vehicle is required to return to the pre-planned optimal route again after recognizing the dynamic obstacle and emergently avoiding the dynamic obstacle, and the process is also called route following.
And (4) carrying out route flyability detection on a result of the pre-planned route, namely checking whether the planned route meets the flight constraint condition of the unmanned aerial vehicle. The planned route meeting the constraint conditions meets the Levy flight trajectory, but the phenomenon of unsmooth route may exist, so an effective route smoothing method is needed.
Fifth, route smoothing method
This patent adopts two points cubic Hermite interpolation method to realize the route smoothing, assumes that the curve to be smoothed is a broken line segment from A to B to C under the rectangular coordinate system, as shown in FIG. 7:
known interpolation node x0,x1The function value and the derivative value are respectively as follows:
yi=f(xi),mi=f′(xi),(i=0,1) (21)
a polynomial H with a degree not exceeding 3 is required3(x) So as to satisfy the following conditions:
H3(xi)=yi,H3′(xi)=mi,(i=0,1) (22)
introducing an interpolation basis function alpha0(x),α1(x),β0(x),β1(x) To obtain H3(x) In that respect Can be expressed as:
H3(x)=y0α0(x)+y1α1(x)+m0β0(x)+m1β1(x) (23)
to satisfy the interpolation condition, the difference basis function has the following limitations:
Figure GDA0003041186690000091
and solving the expressions of the interpolation basis functions by using an undetermined coefficient method:
Figure GDA0003041186690000092
to-be-planned route broken line segment
Figure GDA0003041186690000093
After two-point three times Hermite interpolation, the arc segment is used
Figure GDA0003041186690000094
The method replaces the prior broken line flight path, and realizes the pre-planning route smoothing problem.
Has the advantages that: in previous research work, an artificial potential field method and a chaos strategy are applied to the unmanned aerial vehicle route planning problem, but a heuristic A-star algorithm is adopted as a main method for searching and expanding the route, and compared with a swarm intelligence algorithm, the method is easy to fall into local optimization, and has the problems of low convergence speed and low convergence precision. Therefore, on the basis of the previous research, the method introduces an optimization success rate to the traditional bat algorithm to change the speed updating mode of the bat individuals, simultaneously adopts a chaos method to initialize the distribution of the bat individuals in a search space, utilizes the concept of an artificial potential field to simulate a gravitational field of a terminal point, a starting point and a repulsive field of an obstacle, accelerates the bat individuals to move to an optimal solution speed, and finally provides the improved bat algorithm based on the chaos artificial potential field to solve the problem of unmanned aerial vehicle route planning under the constraint condition.
In the proposed improved algorithm, the introduction of the self-adaptive inertia weight balances the global search and the local search, and avoids the algorithm from falling into the local optimal solution; the chaos strategy is combined to randomize and homogenize the initial distribution of the bat population, so that the solution space is searched and traversed more comprehensively; the artificial potential field is introduced by referring to the characteristics of the unmanned aerial vehicle route planning problem, the information of the terrain, the starting point and the ending point is utilized to a greater extent, and the algorithm execution efficiency is improved.
Drawings
FIG. 1 is a flow chart of a route planning task using the improved algorithm set forth in this patent;
FIG. 2 is a pseudo-code diagram of an improved bat algorithm based on an artificial potential field and a chaos strategy;
FIG. 3 is a diagram of an aircraft in an artificial potential field;
FIG. 4 is a schematic view of a turn radius limitation;
fig. 5 is a schematic diagram of an unmanned aerial vehicle avoiding an obstacle;
FIG. 6 is a schematic diagram of threat cost measurement between split track points;
FIG. 7 is a schematic diagram of Hermite interpolation for implementing airway smoothing;
FIG. 8 shows comparison of the routing effects of BA, DEBA and CPFIBA in a two-dimensional terrain environment;
FIG. 9 is a comparison of convergence curves of the BA, DEBA, CPFIBA route cost functions in the two-dimensional terrain environment;
FIG. 10 illustrates a perspective view and a side view of a three-dimensional terrain environment;
FIG. 11 is a comparison of the effects of the BA, DEBA, CPFIBA route planning in the three-dimensional terrain environment;
FIG. 12 is a comparison of the convergence curves of the fitness function of the BA, DEBA and CPFIBA route planning cost functions in the three-dimensional terrain environment.
Detailed Description
Implementation hardware conditions:
the implementation environment is a PC with 64-bit windows 7 operating system, performance parameters Intel (R) core (TM) i5-3470 CPU 3.20GHz, and ROM size 8 GB. The simulation experiments under the two-dimensional terrain and the three-dimensional terrain are realized by programming under Matlab R2016a software.
Setting related algorithm parameters:
representing an original bat algorithm by BA; DEBA stands for an improved bat algorithm based on a differential evolution algorithm; CPFIBA represents the improved bat algorithm based on the chaos strategy and artificial potential field proposed by this patent.
Basic bat algorithm iterative formula parameter setting: r isi 0=0.6,Ai 00.95, β 0.9; the population size N is 90; the iteration number NC under the two-dimensional environment is 30; the iteration number NC in the three-dimensional environment is 30.
The differential evolution algorithm used for comparison iterates formula parameter settings: NP is 30, a is 0.95, Q is r is F is 0.5; the population size N is 90; the iteration number NC under the two-dimensional environment is 30; the iteration number NC in the three-dimensional environment is 30.
Parameter setting in the improved bat algorithm based on artificial potential field and chaos strategy is mu-4, k-1, m-1, rho00.5, 90 for population size N; the iteration number NC under the two-dimensional environment is 30; the iteration number NC in the three-dimensional environment is 30.
Example 1:
in the two-dimensional environment simulation experiment, the two-dimensional terrain environment overall appearance can be reflected from the air route planning effect diagram, and the terrain diagram is not presented independently. The two-dimensional rectangular coordinates of the starting point in the two-dimensional terrain environment are (0, 0), the coordinates of the target point are (110, 100), and the respective obstacle settings are shown in table 1:
TABLE 1 two-dimensional environmental simulation experiment threat source settings
Figure GDA0003041186690000111
Figure GDA0003041186690000121
The optimal airway length, the airway cost fitness function value and the algorithm execution time comparison (30 times of single experiment iteration, and ten times of independent experiment average) of BA, DEBA and CPFIBA in the two-dimensional environment are shown in table 2:
TABLE 2 comparison of Performance indicators of improved Algorithm in two-dimensional Environment
Figure GDA0003041186690000122
The three algorithms of BA, DEBA and CPFIBA are applied to the unmanned aerial vehicle route planning problem under the two-dimensional complex environment, the two-dimensional route planning simulation effect is obtained and is shown in figure 8, and the route cost fitness function convergence curve is shown in figure 9.
The following analysis can be made for experimental data, simulation result graphs and convergence curve analysis in a two-dimensional environment: under the same terrain condition, the length of the air route planned by CPFIBA is 197.35km, which is shortened by 7.30% compared with 212.89km of the air route planned by DEBA and 25.73% compared with 265.73km of the air route planned by BA, namely the air route planned by CPFIBA is shorter; under the same constraint condition, the flight cost function under the same iteration number is converged. The fitness function for the route cost of CPFIBA is 273.83, 9.26% lower than the fitness value 301.76 of DEBA and 17.06% lower than the fitness value 330.14 of BA. The final fitness convergence value reflects the optimizing precision of each algorithm, and the optimizing precision of CPFIBA obtained from the airway cost fitness function value is higher than that of DEBA and BA. The slope of the curve during initial iteration can be used for obtaining that CPFIBA has higher stability compared with DEBA and BA; from the algorithm execution time perspective analysis, the CPFIBA execution time of 423.83ms saves 13.35% of the execution time compared to the DEBA execution time of 543.92ms, and 59.84% of the execution time compared to the BA execution time of 617.62 ms. From the iteration times when the objective function tends to be stable, the CPFIBA convergence rate is faster than that of DEBA and BA; the comparison of the convergence of the cost function to the final value can result in that CPFIBA is more excellent in the optimizing accuracy. In conclusion, the performance of the CPFIBA in the problem of the two-dimensional air route planning of the unmanned aerial vehicle is better than that of the DEBA and the BA, and the CPFIBA has better applicability.
Example 2:
in the three-dimensional environment simulation experiment, the starting point coordinates are (0, 0, 100), and the target point coordinates are (100, 100, 100). The concept of a safety circle of flight is introduced, i.e. the flying height of the drone cannot be higher than the level of the safety circle of flight, where the safety circle is set parallel to the plane and z is 600 m. The panoramic view and the side view of the three-dimensional unmanned aerial vehicle flight simulation environment terrain model set by using the digital elevation map are shown in fig. 10.
The parameter settings for BA, DEBA and CPFIBA were consistent with the two-dimensional experiments. The three algorithms are respectively utilized to carry out the route planning in the three-dimensional terrain space, the obtained route planning simulation result is shown in fig. 11, and the route cost fitness function convergence curve is shown in fig. 12.
Optimal path lengths, airway cost fitness function values and algorithm execution time comparison of the BA, the DEBA and the CPFIBA in the three-dimensional environment (30 times of single experiment iteration, and ten times of independent experiment average values) are shown in Table 3:
TABLE 3 comparison of Performance indicators of improved Algorithm in three-dimensional Environment
Figure GDA0003041186690000131
The three-dimensional simulation experiment result also comprises experiment data, a simulation effect graph and a flight cost convergence curve. Under the same three-dimensional terrain setting, the length of the air route planned by CPFIBA is 567.37km, which is 27.16% shorter than the length 778.91km of the air route planned by DEBA, and is 36.56% shorter than the length 894.38km planned by BA, and the simulation result shows that the length of the air route planned by CPFIBA is shorter than that of DEBA and BA; from the analysis of the angle of the airway cost adaptability value, the airway cost adaptability function value of the CPFIBA is 110.35, is 42.46 percent less than the DEBA adaptability value 153.47 and 49.53 percent less than the BA adaptability value 218.64, and the final convergence value of the CPFIBA is lower than that of the DEBA and the BA, namely the convergence precision of the CPFIBA is higher than that of the DEBA and the BA; from the analysis of algorithm implementation time, the CPFIBA execution time is 521.34ms, which saves 27.3% than the DEBA execution time 639.05ms, and 56.05% than the BA execution time 813.54ms, and the convergence time of the CPFIBA is shorter than that of the DEBA and the BA, namely the convergence speed is faster than that of the DEBA and the BA. As can be seen from the effect of route planning, the CPFIBA can avoid falling into local optimum, and meanwhile, the convergence precision is greatly superior to DEBA and BA; the convergence curve of the airway cost fitness function verifies the superiority of the convergence precision, the convergence speed and the convergence stability of the CPFIBA relative to the DEBA and the BA. Simulation experiment results show that the CPFIBA algorithm has better applicability than DEBA and BA in three-dimensional environment unmanned aerial vehicle air route planning, and is an air route planning algorithm with practical application significance.

Claims (1)

1. An unmanned aerial vehicle route planning method based on an improved bat algorithm, the improved bat algorithm: firstly, based on the traditional bat algorithm, the optimization success rate is introduced to change the speed updating mode of the bat; secondly, a chaos method is adopted to initialize the distribution of the bat individuals in a search space, so that the search efficiency is improved; thirdly, simulating a gravitational field at a terminal point, a repulsive force field at a starting point and an obstacle by utilizing the concept of an artificial potential field, and accelerating the speed of the bat individual moving to an optimal solution; an improved bat algorithm based on the chaotic artificial potential field is provided by combining the three points; the method is characterized in that:
the method is applied to unmanned aerial vehicle route planning, and specifically comprises the following steps:
step 1: setting necessary parameters in the execution process of the bat algorithm, wherein the necessary parameters comprise iteration times, initial speed, initial position, initial sounding frequency and initial loudness; setting chaotic strategy parameters, and uniformly distributing; setting an artificial potential field, and defining a function of attractive force and repulsive force of a starting point, a target point and an obstacle;
step 2: establishing an unmanned aerial vehicle flight model according to the unmanned aerial vehicle flight cost function and the constraint condition, and determining a performance evaluation mechanism of an unmanned aerial vehicle route planning task;
and step 3: optimizing an unmanned aerial vehicle route planning cost function by using an improved bat algorithm based on an artificial potential field and a chaos strategy to finally obtain a flight track of the unmanned aerial vehicle route planning;
the specific steps of the step 1 are as follows:
step 1.1, under d-dimensional search space, the bat sound frequency f at the time tiVelocity viAnd position xiThe update formula of (2) is:
fi=fmin+(fmax-fmin)×β
vi t=vi t-1+(xi t-1-x*)×f
xi t=xi t-1+vi t
for the update of the sounding frequency, beta is a random subject to uniform distributionA variable, and satisfies beta e [0, 1]];fmaxAnd fminThe maximum value and the minimum value of the initial set sounding frequency range; x in the formula*The position of the current global optimal solution is the optimal value obtained by comparing fitness values of all individuals in the bat population; the bat individual is positioned according to the position of the bat individual at the last moment
Figure FDA0003041186680000021
And global optimal solution x*Measure the acceleration of its movement to the optimal solution; speed of movement at the next moment
Figure FDA0003041186680000022
And also receives the speed of the last moment
Figure FDA0003041186680000023
The influence of (a);
the above is an iterative mechanism followed by the bat population when global search is performed in a solution space, and bats near the global optimal solution generate a local new solution by adopting a random walk rule:
xnew=xold+εAt
in the formula, epsilon [ -1,1 ∈ [ ]]Is a random number;
Figure FDA0003041186680000024
is the average loudness of all bats at time t; r ist=<ri tIs the average sound frequency of all bats at the time t;
the updated formula of the bat individual sound frequency and loudness at the time t is described as follows:
Figure FDA0003041186680000025
ri t+1=ri 0[1-exp(-γt)]
both α and γ are constants, typically taken as 0.9;
step 1.2: the inspiration of the artificial potential field method is derived from the principle that electrostatic field heterogeneous charges generate attraction and homogeneous charges generate repulsion, the acting force between an obstacle and an aircraft in a search space is defined as repulsion, and the acting force between a target point and the aircraft is defined as attraction; applying an artificial potential field method to the problem of route planning, setting barriers and threat sources as repulsive fields, and setting target points as gravitational fields, so that the convergence speed of the whole route planning algorithm is accelerated by the strategy;
Figure FDA0003041186680000026
which is indicative of the repulsive force,
Figure FDA0003041186680000027
which is indicative of the force of attraction,
Figure FDA0003041186680000028
representing the resultant force to which the aircraft is subjected, directly affecting the motion of the aircraft; according to the gradient descent method of the artificial potential field, the repulsive force and the attractive force can be expressed as:
Figure FDA0003041186680000029
Figure FDA00030411866800000210
let n be any point in the search space that is subject to an attractive force Fa(x) And repulsive force Fr(x) Can be expressed as:
Figure 1
step 1.3: the basis of the chaotic algorithm is a logistic mapping:
xn+1=μxn(1-xn-1),n=1,2,.....
when mu is 4, the mapping is a full mapping on the [0, 1] interval, the iteratively generated value is in a pseudo-random distribution state, and when other values are taken, after a certain number of iterations, the generated value converges to a specific value, the bat populations are randomly distributed in the normalized search space, and the global search efficiency of the populations is greatly improved; meanwhile, the bat individuals are extremely difficult to have position transition in a large range in the respective local search process, and the uniform and random distribution of the initial population effectively solves the problem that the bat population is easy to fall into a local optimal solution in the later period of search;
therefore, when an artificial potential field is combined, Logistic mapping based on a chaos strategy is adopted to initialize distribution of the bat population in a solution space, so that the overall convergence speed of the algorithm is accelerated, and the possibility of falling into a local optimal solution is greatly reduced;
the specific steps of the step 2 are as follows:
step 2.1: the group intelligent algorithm has the processes of global search and local search in the optimizing process, and the global search is used for determining the approximate range of the optimal solution; and the local search is carried out to precisely optimize in the ranges;
the balance between the global search and the local search directly influences the search efficiency and the optimization precision, and the bat speed updating formula is rewritten by using the adaptive inertia weight w as follows:
vi t=wvi t-1+(xi t-1-x*)×f
and introducing the concept of the optimizing success rate, so that the inertia weight is adaptively adjusted along with the optimizing success rate of the bat group, and the adaptive inertia weight based on the optimizing success rate is defined as follows:
Figure FDA0003041186680000041
wherein,
Figure FDA0003041186680000042
the success rate of the bat group optimization is achieved; n represents the bat population number;
Figure FDA0003041186680000043
shows the optimizing result of the bat individual i in the process of the t-th iteration,
Figure FDA0003041186680000044
the adaptive value of t generation is superior to t-1 generation, and a better solution is searched, if not, the optimal solution is obtained
Figure FDA0003041186680000045
The self-adaptive inertia weight of the optimizing success rate is introduced, and the method is different from the traditional linear inertia weight, and better balances the proportion of global search and local search. The introduction of the optimizing success rate is a parameter adjusting method with a feedback mechanism, the optimizing result of the global search influences the value of the self-adaptive inertia weight, so that the speed of the bat individual and the time for the bat individual to enter the local search are changed;
step 2.2: the cost function of the unmanned aerial vehicle in the flying process can be divided into three parts, namely, the route length cost, the threat cost and the fuel consumption cost; the total cost function is denoted by J, and the minimization problem with respect to the total cost function J is defined as follows:
min J=k1JL+k2JT+(1-k1-k2)JF
wherein, JLIs the cost of the length of the flight path, JTIs the cost of the threat, JFRefers to the fuel consumption penalty incurred in flight; k is a radical of1,k2Are all normal numbers and satisfy k being more than or equal to 01≤1,0≤k2≤1;
For the length of flight path cost JLThe definition is as follows:
Figure FDA0003041186680000046
wherein L is the length of the total flight path, LijThe length of the track segment is used for calculating the length of the air route which cannot be calculated by integral calculation before the air route is smoothly carried out;
cost J for threatTThe definition is as follows:
Figure FDA0003041186680000051
wherein, tkThe threat factor is a measure of the threat degree of the threat source to the unmanned aerial vehicle; n is a radical oftRepresenting the total number of threat sources; the current coordinate of the unmanned aerial vehicle is (x, y); the center coordinate of the threat source is (x)k,yk);
For fuel consumption cost JFThe definition is as follows:
Figure FDA0003041186680000052
wherein,
Figure FDA0003041186680000053
k in (1) represents the fuel cost consumed by the unmanned aerial vehicle for driving a path with a unit length;
Figure FDA0003041186680000054
in H is the height of the unmanned aerial vehicle flight safety circle, the flight height of the unmanned aerial vehicle does not exceed the value, w0The fuel cost required to be consumed when the unmanned aerial vehicle maintains the height H is represented, and H represents the current height of the unmanned aerial vehicle;
the specific steps of the step 3 are as follows:
step 3.1: the unmanned aerial vehicle needs to accord with the self dynamic characteristics in the flying process, a part of route points in the flight path are selected as a model for explanation, and A (x) is seti-1,yi-1,zi-1) As the last waypoint, B (x)i,yi,zi) As the current waypoint, C (x)i+1,yi+1,zi+1) For the next waypoint and records
Figure FDA0003041186680000055
A track point migration vector is obtained;
when wanting to change self height among the unmanned aerial vehicle flight process, need do and climb or dive operation, assume that the biggest climbing is theta with dive angle, then to unmanned aerial vehicle climb/dive angle's restraint do:
Figure FDA0003041186680000056
when the unmanned aerial vehicle avoids the obstacle, the situation that the turning radius corresponding to the planned track inscribed circle is too small may occur, and the flight characteristics of the unmanned aerial vehicle are violated; therefore, the minimum turning radius of the unmanned aerial vehicle in flight is limited, and the included angle of the air route
Figure FDA0003041186680000057
The minimum turning angle may be determined by:
Figure FDA0003041186680000061
wherein r isminIs the minimum turning radius;
the absolute height H of the ground has the maximum value that H is limited to be less than or equal to H, wherein H is the height of the plane where the flight safety circle is located; and in order to react to terrain changes in time, the relative height h of the unmanned aerial vehicle from the terrain surface is requiredi≥hmin
Step 3.2: the unmanned aerial vehicle detects static obstacles and dynamic obstacles by utilizing a sensor thereof, and an avoidance mechanism of the obstacles after the obstacles are identified is a key factor influencing the route planning effect;
segmenting the sub-track segment again according to the number of the threat sources to calculate the threatCost of
Figure FDA0003041186680000062
In general, this threat cost may be expressed as:
Figure FDA0003041186680000063
wherein N isTRepresenting the number of threat sources, tjThreat factor representing a source of threat, dk(i, j) represents a straight-line distance between the starting points i and j of the kth branch track segment; liRepresenting the number of track segments;
the obstacle avoidance process is the process of minimizing the threat cost under the condition of meeting the constraint condition, the advantage of the optimization problem is solved by using an intelligent algorithm, and the pre-planned optimal route can effectively avoid the static obstacle threat in the space; for avoiding the dynamic obstacle, the unmanned aerial vehicle is required to return to the preplanned optimal route again after recognizing the dynamic obstacle and emergently avoiding the dynamic obstacle, and the process is also called route following;
step 3.3: a two-point three-time Hermite interpolation method is adopted to realize airway smoothing;
known interpolation node x0,x1The function value and the derivative value are respectively as follows:
yi=f(xi),mi=f′(xi), i=0,1
a polynomial H with a degree not exceeding 3 is required3(x) So as to satisfy the following conditions:
H3(xi)=yi,H3′(xi)=mi, i=0,1
introducing an interpolation basis function alpha0(x),α1(x),β0(x),β1(x) To obtain H3(x) (ii) a Can be expressed as:
H3(x)=y0α0(x)+y1α1(x)+m0β0(x)+m1β1(x)
to satisfy the interpolation condition, the difference basis function has the following limitations:
Figure FDA0003041186680000071
βi(xj)=0,βi′(xj)=δij, i,j=0,1
and solving the expressions of the interpolation basis functions by using an undetermined coefficient method:
Figure FDA0003041186680000072
Figure FDA0003041186680000073
to-be-planned route broken line segment
Figure FDA0003041186680000074
After two-point three times Hermite interpolation, the arc segment is used
Figure FDA0003041186680000075
The method replaces the prior broken line flight path, and realizes the pre-planning route smoothing problem.
CN201811091770.6A 2018-09-19 2018-09-19 Unmanned aerial vehicle route planning method based on improved bat algorithm Active CN109144102B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811091770.6A CN109144102B (en) 2018-09-19 2018-09-19 Unmanned aerial vehicle route planning method based on improved bat algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811091770.6A CN109144102B (en) 2018-09-19 2018-09-19 Unmanned aerial vehicle route planning method based on improved bat algorithm

Publications (2)

Publication Number Publication Date
CN109144102A CN109144102A (en) 2019-01-04
CN109144102B true CN109144102B (en) 2021-08-20

Family

ID=64814948

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811091770.6A Active CN109144102B (en) 2018-09-19 2018-09-19 Unmanned aerial vehicle route planning method based on improved bat algorithm

Country Status (1)

Country Link
CN (1) CN109144102B (en)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109870906B (en) * 2019-02-25 2020-06-12 北京航空航天大学 BBO (broadband barrier optimization) based high-speed rotorcraft path planning method for optimizing artificial potential field
CN111832725B (en) * 2019-04-15 2023-05-26 电子科技大学 Multi-robot multi-task allocation method and device based on improved genetic algorithm
CN110274607A (en) * 2019-05-21 2019-09-24 深圳壹账通智能科技有限公司 Intelligent paths planning method, device and computer readable storage medium
CN110274608A (en) * 2019-05-21 2019-09-24 深圳壹账通智能科技有限公司 Intelligent paths planning method, device and computer readable storage medium
CN110162077B (en) * 2019-06-18 2022-04-05 哈尔滨工程大学 Unmanned aerial vehicle flight path planning method based on flying fish algorithm
CN110427046B (en) * 2019-07-26 2022-09-30 沈阳航空航天大学 Three-dimensional smooth random-walking unmanned aerial vehicle cluster moving model
CN112444248B (en) * 2019-08-27 2022-12-27 广州极飞科技股份有限公司 Route generation method, device, equipment and storage medium
CN110632922B (en) * 2019-09-10 2022-06-17 青岛理工大学 Path planning method based on bat algorithm and reinforcement learning
CN110989656A (en) * 2019-11-13 2020-04-10 中国电子科技集团公司第二十研究所 Conflict resolution method based on improved artificial potential field method
CN110967016B (en) * 2019-11-22 2021-10-22 中国人民解放军63629部队 Off-line planning method and device for aircraft route and computer equipment
CN111123976A (en) * 2019-12-24 2020-05-08 一飞智控(天津)科技有限公司 Unmanned aerial vehicle cluster path planning processing method based on artificial potential field and unmanned aerial vehicle
CN111189455B (en) * 2020-01-14 2022-02-01 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) Unmanned aerial vehicle route planning method, system and storage medium
DE112021001104T5 (en) * 2020-02-19 2023-01-12 Fanuc Corporation MOVEMENT PLANNING METHOD WITH COLLISION AVOIDANCE FOR INDUSTRIAL ROBOTS
CN111813144B (en) * 2020-06-11 2022-02-18 南京航空航天大学 Multi-unmanned aerial vehicle collaborative route planning method based on improved flocks of sheep algorithm
CN111880561B (en) * 2020-07-16 2023-03-28 河南大学 Unmanned aerial vehicle three-dimensional path planning method based on improved whale algorithm in urban environment
CN111930121B (en) * 2020-08-10 2022-10-25 哈尔滨工程大学 Mixed path planning method for indoor mobile robot
CN111969488A (en) * 2020-08-21 2020-11-20 福州大学 Coordinate transformation-based unmanned aerial vehicle inspection obstacle avoidance method for high-voltage line overhead area
CN112462803B (en) * 2020-11-27 2022-06-17 北京工商大学 Unmanned aerial vehicle path planning method based on improved NSGA-II
CN112880688B (en) * 2021-01-27 2023-05-23 广州大学 Unmanned aerial vehicle three-dimensional track planning method based on chaotic self-adaptive sparrow search algorithm
CN113190037A (en) * 2021-04-08 2021-07-30 上海吞山智能科技有限公司 Unmanned aerial vehicle optimal path searching method based on improved fluid disturbance and sparrow algorithm
CN113759958B (en) * 2021-07-07 2023-08-01 哈尔滨工程大学 Unmanned aerial vehicle formation track planning method based on positioning precision
CN113359773B (en) * 2021-07-07 2024-08-09 大连海事大学 Unmanned ship navigation path decision method and system
CN114527746A (en) * 2022-01-12 2022-05-24 燕山大学 Intelligent automobile path planning and obstacle avoidance tracking method
CN114895707B (en) * 2022-05-13 2023-06-30 华南农业大学 Agricultural unmanned aerial vehicle path planning method and system based on variable frequency bat algorithm
CN114779821B (en) * 2022-05-25 2023-06-27 四川大学 Unmanned aerial vehicle self-adaptive repulsive force coefficient path planning method based on deep learning
CN116400737B (en) * 2023-06-02 2023-08-25 中国传媒大学 Safety path planning system based on ant colony algorithm
CN117556979B (en) * 2024-01-11 2024-03-08 中国科学院工程热物理研究所 Unmanned plane platform and load integrated design method based on group intelligent search
CN117647997B (en) * 2024-01-29 2024-04-16 中国人民解放军战略支援部队航天工程大学 Knowledge bidirectional migration unmanned aerial vehicle collaborative track local re-planning method and system
CN118113068A (en) * 2024-03-01 2024-05-31 中国人民解放军海军工程大学 Unmanned cluster route planning method and system based on improved particle swarm algorithm
CN117913962B (en) * 2024-03-19 2024-05-14 哈尔滨工业大学 MPPT power supply control method for satellite based on improved BA-P & O hybrid algorithm

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104406593A (en) * 2014-12-03 2015-03-11 广西民族大学 Method for determining optimal route of airway of unmanned aerial vehicle
CN105427348A (en) * 2015-12-03 2016-03-23 山东理工大学 Video object tracking method based on bat algorithm
CN107103397A (en) * 2017-06-26 2017-08-29 广东工业大学 A kind of traffic flow forecasting method based on bat algorithm, apparatus and system
CN107253194A (en) * 2017-07-31 2017-10-17 中南大学 A kind of carrying machine human arm manipulation multiple spot mapping intelligent control method and system
CN107272705A (en) * 2017-07-31 2017-10-20 中南大学 A kind of multiple neural network controlling planning method of robot path under intelligent environment
CN107392317A (en) * 2017-07-28 2017-11-24 中南大学 A kind of neutral net colony mixing computational methods of intelligent environment carrying robot identification floor

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104406593A (en) * 2014-12-03 2015-03-11 广西民族大学 Method for determining optimal route of airway of unmanned aerial vehicle
CN105427348A (en) * 2015-12-03 2016-03-23 山东理工大学 Video object tracking method based on bat algorithm
CN107103397A (en) * 2017-06-26 2017-08-29 广东工业大学 A kind of traffic flow forecasting method based on bat algorithm, apparatus and system
CN107392317A (en) * 2017-07-28 2017-11-24 中南大学 A kind of neutral net colony mixing computational methods of intelligent environment carrying robot identification floor
CN107253194A (en) * 2017-07-31 2017-10-17 中南大学 A kind of carrying machine human arm manipulation multiple spot mapping intelligent control method and system
CN107272705A (en) * 2017-07-31 2017-10-20 中南大学 A kind of multiple neural network controlling planning method of robot path under intelligent environment

Also Published As

Publication number Publication date
CN109144102A (en) 2019-01-04

Similar Documents

Publication Publication Date Title
CN109144102B (en) Unmanned aerial vehicle route planning method based on improved bat algorithm
US11727812B2 (en) Airplane flight path planning method and device based on the pigeon-inspired optimization
Lin et al. A Novel Improved Bat Algorithm in UAV Path Planning.
CN107504972B (en) A kind of aircraft's flight track method and device for planning based on dove group's algorithm
CN111381600B (en) UUV path planning method based on particle swarm optimization
Wen et al. UAV online path planning algorithm in a low altitude dangerous environment
CN110134140B (en) Unmanned aerial vehicle path planning method based on potential function reward DQN under continuous state of unknown environmental information
JP2023546810A (en) Vehicle trajectory planning method, vehicle trajectory planning device, electronic device, and computer program
CN109871031B (en) Trajectory planning method for fixed-wing unmanned aerial vehicle
CN111930121B (en) Mixed path planning method for indoor mobile robot
Tordesillas et al. Real-time planning with multi-fidelity models for agile flights in unknown environments
CN109597425A (en) Navigation of Pilotless Aircraft and barrier-avoiding method based on intensified learning
CN108919818B (en) Spacecraft attitude orbit collaborative planning method based on chaotic population variation PIO
CN113848919A (en) Ant colony algorithm-based indoor AGV path planning method
CN113268074B (en) Unmanned aerial vehicle flight path planning method based on joint optimization
CN116804879B (en) Robot path planning framework method for improving dung beetle algorithm and fusing DWA algorithm
CN114518770B (en) Unmanned aerial vehicle path planning method integrating potential field and deep reinforcement learning
CN111813144B (en) Multi-unmanned aerial vehicle collaborative route planning method based on improved flocks of sheep algorithm
Li et al. A behavior-based mobile robot navigation method with deep reinforcement learning
CN109870906B (en) BBO (broadband barrier optimization) based high-speed rotorcraft path planning method for optimizing artificial potential field
CN113391633A (en) Urban environment-oriented mobile robot fusion path planning method
CN115903888A (en) Rotor unmanned aerial vehicle autonomous path planning method based on longicorn swarm algorithm
Cui et al. UAV path planning method for data collection of fixed-point equipment in complex forest environment
Chronis et al. Dynamic navigation in unconstrained environments using reinforcement learning algorithms
CN116698037B (en) Unmanned aerial vehicle track planning method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant