A Path Optimization Algorithm for Multiple Unmanned Tractors in Peach Orchard Management

: In order to improve the management efﬁciency of peach orchards, this paper considers the cooperative operation scheme of multiple unmanned tractors. According to the actual situation, this paper constructs the path planning model of multiple unmanned tractors in a standard peach orchard, designs the objective function to optimize the total turning time and total operating time according to the tractor driving parameters, and solves it by improving the differential evolution algorithm. Aiming at the premature convergence problem, the permutation matrix is introduced to represent the driving paths of multiple unmanned tractors. Then, the dynamic parameters are adopted to make the parameters change with the number of iterations, and the elite selection strategy is used to eliminate the redundant feasible solutions. An Adaptive Elite Differential Evolution (AEDE) algorithm suitable for multi-tractor path optimization is proposed. The results show that, compared with the traditional Differential Evolution algorithm (Differential Evolution, DE), the total turning time and total operating time in the rectangular peach orchard optimized by AEDE are reduced by 3.34% and 0.87%, respectively. Compared with the block operation, the total turning time and total operating time of the AEDE-optimized rectangular peach orchard operation path were reduced by 37.37% and 9.47%, respectively. Experiments show that AEDE, which optimizes the operating path of multi-tractors in standard peach orchards, is able to improve the efﬁciency and reduce the operating time.


Introduction
The cultivation area of peaches in China is 55,317.23 (km 2 ), and the planting scale ranks first in the world [1]. Peach orchard management includes weeding, ditching, fertilizing, and soil cultivation, which requires a large amount of labor, contradicting the increasingly serious aging of the rural labor force. At present, China is vigorously promoting the integration of agricultural machinery and agronomy and standardized farmland transformation. The integrated path planning systems of agricultural unmanned tractors can perceive and act independently [2]. By connecting different agricultural tools to carry out different operations in standard peach orchards, the labor force can be reduced and the planting management efficiency can be improved. Optimizing the driving path is a key to improving unmanned production efficiency and operation quality [3].
Spekken et al. mentioned that during agricultural machinery operation, most of the time was used for non-productive operations, such as turning in non-working areas. Reasonable path planning can reduce the working time of agricultural machinery and improve operation efficiency [4]. At present, the B-Pattern model [5] proposed by Bochtis was adopted as the path planning model in most studies. Based on the B-Pattern model, the tractor operation scheduling is optimized by intelligent optimization methods, such as genetic algorithm [6][7][8], simulated annealing algorithm [9,10], ant colony algorithm [11,12], domain search [13], tabu search algorithm [14,15], particle swarm optimization algorithm [16], and Clark-Wright algorithm [17,18]. In addition, there is research on operation path optimization based on probability. Utama et al. introduced domain search and elite selection to improve the effect of distribution estimation algorithms in agricultural machinery operation path planning [19].
Fabre et al. [20] transformed the path planning problem of agricultural machinery into a Hamiltonian problem and proposed a path optimization algorithm based on a greedy strategy to minimize the repetitive operation area. Jensen et al. transformed the path planning problem of agricultural machinery into a Traveling Salesman Problem (TSP) problem and proposed an algorithm based on the state-space search to reduce the length of non-operational paths [21]. Muthukumaran et al. adopted the hybrid dragonfly cuckoo search algorithm to optimize the path coverage based on the TSP model [22], and the computational efficiency and effect were better than most intelligent optimization algorithms. Postmes transformed the operation coverage optimization of agricultural machinery robots into VRP (Vehicle Routing Problem) and tried to solve the problem with various approximation algorithms [23]. In order to reduce the complexity of agricultural machinery path planning, Oksanenand et al. used trapezoidal decomposition and an incremental method to simplify the path coverage problem [24].
At present, the studies on autonomous driving path planning in farmland mainly focus on food production, including global path planning and local path planning. The traditional B-Pattern planning model is mainly used in the field where the spacing of operation lines is equal (operation line refers to the route that the tractor starts from one end of the orchard and ends at the opposite end, such as transplanting and harvesting in the field). When the unmanned driving system is applied to the management of a peach or-chard (such as ditching, weeding, fertilization, soil tillage), the operation is usually one-sided, resulting in uneven spacing between operation lines. Therefore, the path autonomous decision-making in the peach orchards cannot be directly applied to the B-Pattern model [25].
The spacing of operation lines in the peach orchard is unequal, and peach trees are densely distributed. This paper studies a multi-machine cooperative operation path optimization method that is suitable for standard peach orchard management, which provides a reference for the continuous and stable operation of unmanned systems.

Path Optimization Model of Unmanned Tractor
The distribution of peach trees in the standard peach orchard is shown in Figure 1. The peach orchard is divided into Upper headland (UH), working area, and Lower headland (LH). The working area is the area where the tractors manage the peach trees, and the headland is the area where the tractor turns. Suppose there are n rows of peach trees in the peach orchard-vectors x k and y k represent the x-coordinate and y-coordinates of each peach tree in row k, respectively. Additionally, x m k represents the x-coordinate of the peach tree m in row k (counted from bottom to top, consistent with the positive direction of Y-axis), where 1 ≤ k ≤ n, 1 ≤ m ≤ n k, n k represents the number of peach trees in row k. The coordinates of peach trees at both ends of each row are represented by n × 2-dimensional matrices U and L, and U k,1 and U k,2 respectively represent the horizontal and vertical coordinates of peach trees at the junction of headland UH and working area in row k, i.e., x k n k and y k n k . Similarly, L k,1 and L k,2 respectively represent the horizontal and vertical coordinates of the peach trees in row k at the junction of the headland LH and working area, i.e., x k 1 and x k 1 . The boundary of the headland and the working area in the standard peach orchard can be determined according to the matrices U and L. does not work in sequence on both sides of the peach trees. Therefore, when the tractor works in the standard peach orchard, the driving path is not the equidistant operation line in the traditional B-Pattern model. If the peach orchard has n rows of peach trees, the tractor needs to make at least 2n to complete the management of all peach trees. The dashed line in Figure 1 represents the track of the tractor when it is working in the standard peach orchard. Although it is also composed of several parallel lines, they are not equidistant from each other. Unmanned tractors need to make turns in the field before working on the next row of peach trees. There are three main types of turns, including the Π-Turn, Ω-turn, and Tturn. As the T-turn requires the tractor to reverse and stop, it will not only consume time, but also aggravate the tractor wear. This paper considers only the Π-Turn and Ω-turn, as shown in Figure 2a,b. In Figure 2, w and d represent the horizontal and vertical distance of two adjacent rows of peach trees, respectively, and r represents the turning radius of the tractor. Based on the above path planning model, this paper designs a path optimization scheme for multiple unmanned tractors. The design of a multi-tractor driving route can avoid incomplete operation or repeated operation. The simplest scheme is to divide the peach orchard into several sub areas, and then assign different tractors to operate in the divided areas. The schematic diagram of block operation is shown in Figure 3.
Block operation is the operation mode of an S-type operation for multiple tractors, which requires multiple turns and has low operation efficiency. Reasonably arranging the In order to ensure that peach trees have enough space for growth and management, the spacing of peach trees is much larger than the width of tractors. As shown in Figure 1, the tractor needs to be driven at least twice to manage a row of peach trees, and probably does not work in sequence on both sides of the peach trees. Therefore, when the tractor works in the standard peach orchard, the driving path is not the equidistant operation line in the traditional B-Pattern model. If the peach orchard has n rows of peach trees, the tractor needs to make at least 2n to complete the management of all peach trees. The dashed line in Figure 1 represents the track of the tractor when it is working in the standard peach orchard. Although it is also composed of several parallel lines, they are not equidistant from each other.
Unmanned tractors need to make turns in the field before working on the next row of peach trees. There are three main types of turns, including the Π-Turn, Ω-turn, and T-turn. As the T-turn requires the tractor to reverse and stop, it will not only consume time, but also aggravate the tractor wear. This paper considers only the Π-Turn and Ω-turn, as shown in Figure 2a,b. In Figure 2, w and d represent the horizontal and vertical distance of two adjacent rows of peach trees, respectively, and r represents the turning radius of the tractor.
In order to ensure that peach trees have enough space for growth and manag the spacing of peach trees is much larger than the width of tractors. As shown in 1, the tractor needs to be driven at least twice to manage a row of peach trees, and pr does not work in sequence on both sides of the peach trees. Therefore, when the works in the standard peach orchard, the driving path is not the equidistant operat in the traditional B-Pattern model. If the peach orchard has n rows of peach tre tractor needs to make at least 2n to complete the management of all peach tre dashed line in Figure 1 represents the track of the tractor when it is working standard peach orchard. Although it is also composed of several parallel lines, t not equidistant from each other. Unmanned tractors need to make turns in the field before working on the ne of peach trees. There are three main types of turns, including the Π-Turn, Ω-turn, turn. As the T-turn requires the tractor to reverse and stop, it will not only consum but also aggravate the tractor wear. This paper considers only the Π-Turn and Ωshown in Figure 2a,b. In Figure 2, w and d represent the horizontal and vertical d of two adjacent rows of peach trees, respectively, and r represents the turning ra the tractor.  Based on the above path planning model, this paper designs a path optim scheme for multiple unmanned tractors. The design of a multi-tractor driving ro avoid incomplete operation or repeated operation. The simplest scheme is to div peach orchard into several sub areas, and then assign different tractors to operat divided areas. The schematic diagram of block operation is shown in Figure 3.
Block operation is the operation mode of an S-type operation for multiple t which requires multiple turns and has low operation efficiency. Reasonably arrang Based on the above path planning model, this paper designs a path optimization scheme for multiple unmanned tractors. The design of a multi-tractor driving route can avoid incomplete operation or repeated operation. The simplest scheme is to divide the peach orchard into several sub areas, and then assign different tractors to operate in the divided areas. The schematic diagram of block operation is shown in Figure 3.

Optimization Model of Minimum Operating Time for Unmanned Tractor
The term "operation line" refers to the route that the tractor starts from one end of the orchard and ends at the opposite end. E represents the set of all operation lines, V is the set of end points of all operation lines, and the driving path is represented by a directed graph G = {V,E}. Graph G is mainly composed of vertices and connecting edges of parallel lines, so it is connected. To improve the management efficiency of the standard peach orchard, it is considered to formulate a cooperative scheme for multiple tractors to ensure the shortest working time. As shown in Equation (1), this paper establishes an objective function to minimize the total turning time and working time of multiple tractors, which ignores the time of tractors entering and leaving the orchard.
where, tturn-total turning time of multiple unmanned tractors twork-total operating time of unmanned tractors M-the number of unmanned tractors tall-weighted average of total turning time and total operating time z-the weight variable represents farmers' attention to the total turning time and the total operating time tturn can be calculated according to Equation (2): where,  Block operation is the operation mode of an S-type operation for multiple tractors, which requires multiple turns and has low operation efficiency. Reasonably arranging the operation sequence of the tractor and reducing the Ω-turn is conducive to improving the overall operation efficiency.

Optimization Model of Minimum Operating Time for Unmanned Tractor
The term "operation line" refers to the route that the tractor starts from one end of the orchard and ends at the opposite end. E represents the set of all operation lines, V is the set of end points of all operation lines, and the driving path is represented by a directed graph G = {V,E}. Graph G is mainly composed of vertices and connecting edges of parallel lines, so it is connected.
R represents the set of tractor paths, and R k represents the driving path of the k-th tractor. R k is composed of a vertex set <R 1 k , R 2 k , . . . , R k m k >, where m k represents the number of vertices of path R k .
To improve the management efficiency of the standard peach orchard, it is considered to formulate a cooperative scheme for multiple tractors to ensure the shortest working time. As shown in Equation (1), this paper establishes an objective function to minimize the total turning time and working time of multiple tractors, which ignores the time of tractors entering and leaving the orchard. t work can be calculated according to Equation (3): -operating time of tractor in the i-th operation line, which is calculated by the operation line length and tractor speed. Assuming that only one tractor is allowed to work in each operation line, each operation line is traversed once, and constraints (4) and (5) can be obtained. Under the constraint of this model, the number of operation lines is 2n, and the number of turning paths is 2n-M.
When an unmanned tractor enters the next operation line from one operation line, C i,j describes the time cost of the tractor from line i to line j, and its definition is shown in Equations (6) and (7).
When i = j, C i,j represents the work time in the i-th operation line. When i = j, C i,j means turning from the i-th operation line to the j-th operation line. At this time, there are two possibilities: 1 the tractor changes the operation line at the headland UH; 2 the tractor converts the operation line in the headland LH. Considering that peach trees may be asymmetrically distributed on both sides of the headland, the time cost of changing the operation line at headland UH may be different from that at headland LH; that is, c i,j UH may not be equal to c i,j LH . Based on the working time of the tractor in the operation line and the turning time of the tractor in the headland, the matrix C representing the time cost of the tractor changing the operation line is constructed in this paper: The calculation of c i,j UH and c i,j LH is related to the turning model used in path planning. In order to reduce turning time and improve operation efficiency, this paper adopts the Π-Turn model and Ω-Turn model. Assuming that c i,j D , D ∈ {UH,LH} represents the time cost of changing from line i to line j at the headland, and the time cost of the Π-Turn model and Ω-Turn model can be calculated by Equations (9) and (10), respectively.
where, d i,j D -vertical coordinate difference between the i-th operation line and the j-th operation line at the endpoint near the headland D.
w i,j D -abscissa coordinate difference between the i-th operation line and the j-th operation line at the endpoint near the headland D.
v t -the turning speed of the tractor at the end of the headland. v s -the straight speed of the tractor at the end of the headland.
x-coordinate of the j-th operation line. After calculating all U k,1 ± W (1 ≤ k ≤ n), X i U and X j U are sorted, and similarly, X i L is obtained, where W is the working width of the tractor.
According to W and r, c i,j D can be obtained: Based on the number |E| of operation lines, the number M of unmanned tractors, and the cost matrix C , the objective function is obtained as shown in Equation (1), where the set of traveled operation lines by each tractor is used as the variable.
A reasonable path optimization scheme for multi-unmanned tractors can reduce the total operating time and improve work efficiency [26]. However, this problem belongs to the TSP [27,28], and it is difficult to calculate the optimal solution in a short time. Differential evolution algorithm is an intelligent optimization algorithm based on population evolution proposed by Storn [29,30]. It has the characteristics of fewer operators and strong search ability. The solution can be obtained through mutation, crossover, and selection to update the population [31][32][33]. This paper designs a differential evolution algorithm to obtain multi-tractor paths.

Differential Evolution Algorithm
Differential evolution algorithm needs to balance the global search and local search in the iterative process [34]. A global search is to search different spaces that have not appeared. A local search means the ability to improve the solution quality according to the information of the current solution. Differential evolution algorithms can obtain effective solutions through mutation and crossover operations. However, when differential evolution algorithm pays attention to local search, it is easy to fall into local optimization, resulting in premature convergence. On the other hand, if the differential evolution algorithm focuses on global search, the convergence will be very slow [35].
The mutation is also called mutation operator, which calculates the position of the new individual according to the individual position of the last iteration and the optimal position of all individuals. Differential evolution algorithm mainly has six mutation operators, and Equation (12) gives a common mutation operator [36].
where, t-iteration number. i-individual ID in algorithm population. κ i t -the position of the i-th individual in the t-th iteration, |E|-dimensional column vector.
κ best t -the individual position with the best fitness in the t-th iteration, |E|-dimensional column vector.
F-a scaling factor that controls the scale of the difference. The crossover of differential evolution algorithm is shown in Equation (13), where there is a random variable subject to the uniform distribution. When the random variable is less than the crossover probability, accept the mutation operation, otherwise, keep the original value.
where, CR-crossover probability. rand-random number obeying uniform distribution [0,1). The selection can be represented by Equation (14). If the fitness of the new solution is less than that of the original solution, the new solution is accepted, otherwise the original solution is retained.
Although the differential evolution algorithm can obtain the minimum driving cost solution of the unmanned tractors in a limited time, it has the problem of premature convergence and insufficient exploration of solution space. In addition, the mutation operator generates a sequence of real numbers, not an integer sequence in a given range. Multi-tractor cooperation requires a finite number of integers to indicate the number of unmanned tractors and the current scheduling state. The traditional differential evolution algorithm cannot be directly used to solve the multi-tractor path optimization problem. When solving this kind of problem, they usually optimize the continuous variables first, and then indirectly obtain the discrete operation line travelling scheme through the round rather than directly optimizing the discrete operation line travelling scheme. Therefore, this paper provides a better balance between global search and local exploration by improving the mutation operator and crossover operator to optimize the differential evolution algorithm so as to obtain an efficient multi-tractor cooperative path optimization scheme.

Mutation Operator
A permutation matrix is not only an invertible matrix but also an orthogonal matrix. It can greatly reduce the convergence time by transposing when finding the inverse. In the cooperative path optimization problem for multi-tractors, the tractor sequence κ i is one of the permutations composed of 1~|E|, and κ i is a column vector. According to the definition of the permutation matrix, the column vector x i obtained by left multiplying a permutation matrix is still one of the permutations composed of 1~|E|; that is, it is still a legal operation sequence after the transformation of the permutation matrix. Therefore, this paper introduces a permutation matrix to optimize the mutation operator.
Let the position κ of the i-th individual in the t iteration be: where, P i t -the permutation matrix of the individual i in the t-th iteration, with the size of |E|×|E|.
α-|E|-dimensional column vector, the expression is [1 2 . . . |E|] T . Substitute Equation (15) into mutation operator (Equation (12)). Since α is a non-zero vector, the simplified operator (16) is obtained: However, when the permutation matrices are added or subtracted from each other, it is difficult to ensure the orthogonality of the matrix, so it can not be calculated directly by Equation (16).
Two permutation matrices are still permutation matrices after point multiplication, and the inverse cost is very low. Point multiplication can be used instead of plus; dot inverse matrix can be used instead of subtraction; let F = 1; Equation (16) can be rewritten as: However, in Equation (17), the time complexity of calculating the inverse matrix is at least the cubic, and the calculation cost is large, which is not conducive to the fast acquisition of planning paths. The inverse operation of the orthogonal matrix is equivalent to the transpose operation, and the transpose operation is simpler. Therefore, this paper rewrites Equation (17) into Equation (18) with clearer mathematical meaning and less calculation: This paper takes Equation (18) as the mutation operator of the improved differential evolution algorithm so as to optimize the operation path of the unmanned tractors.

Dynamic Parameters
In the traditional differential evolution algorithm, the crossover probability CR and scaling factor F are fixed and do not change with iteration. CR affects the probability of mutation. The smaller its value, the smaller the probability of mutation. Currently, the local search ability of the algorithm is strong. On the contrary, the greater the mutation probability, the stronger the global search ability of the algorithm. If F is fixed, the global search and local search cannot change adaptively, resulting in premature convergence and low search accuracy. As shown in Equation (19), the dynamic parameter scheme is adopted in this paper to make CR changes with the number of iterations. n is the number of operation lines, t is the current number of iterations, and T is the maximum number of iterations.
When t increases from 1 to T, the crossover probability CR gradually decreases to 1/n. The global search of the algorithm changes from strong to weak, and the local search changes from weak to strong. This can avoid premature convergence and improve search accuracy.
Similarly, too small a scaling factor F will also affect the search accuracy and convergence speed. This paper designs an adaptive function as shown in Equation (20). When t increases from 1 to T, the scaling factor F gradually increases to 1. The global search of the mutation operator changes from strong to weak, and the local search changes from weak to strong, which promotes the convergence of the algorithm.

Elite Selection
Elite selection can improve the convergence speed and local search performance. It is assumed that κ n−2 t , κ n−1 t , κ n t are the first three individuals with the worst fitness in the t-th iteration and κ 1 t−1 , κ 2 t−1 , κ 3 t−1 are the first three individuals with the best fitness in the (t − 1)-th iteration. κ n−2 t , κ n−1 t , κ n t can be replaced by κ 1 t−1 , κ 2 t−1 , κ 3 t−1 to realize the elite selection.
This paper introduces the permutation matrix, adopts dynamic parameters and elite selection strategy to improve the differential evolution algorithm, and proposes an Adap- information, tractor parameters, and the number of tractors. Then, the crossover probability CR and scaling factor F of the current iteration are calculated, and the mutation operator and crossover operator designed in this paper are used to generate a new multitractor cooperative operation path optimization. Then, the solution with the worst fitness is eliminated through the elite selection, and the three schemes with the best fitness in the last iteration are introduced. The scheme with better fitness is selected through the selection operation to replace the current scheme until the termination condition is reached. This scheme with the best current fitness is output as the optimal multi machine cooperative path planning scheme solved by the algorithm.

Experiment and Analysis
This paper took the John Deere 6B-1204 tractor commonly used in planting as the research object, and the total number of tractors was six. The working width of each tractor was W = 3.5 m, the minimum turning radius was rmin = 5.5 m, the working speed was vw = 1 m·s −1 in the working area, the straight speed was vs = 5 m·s −1 in the non-working area, and the turning speed was vt = 2 m·s −1 . The peach orchard was rectangular, and the row spacing of peach trees was 4 × 5 m. As the peach tree distribution will affect the algorithm performance, this paper tests the path planning effect of the proposed model and the improved algorithm with different peach tree distributions. Figure 5 shows the distribution of peach trees in three rectangular standard peach orchards. Figure 5a shows the peach orchards with the rectangular distribution of peach trees. There were 20 rows of peach trees in the peach orchard, with 15 trees in each row. Figure 5b shows peach orchards with the approximately trapezoidal distribution of the peach trees, and Figure  5c shows peach orchards with an irregular distribution of peach trees. The experiment assumes that the tractor works in a peach orchard without obstacles and the driving process is continuous. The multi-tractor cooperative path optimization method studied in this paper firstly generates several feasible driving paths randomly according to the input orchard information, tractor parameters, and the number of tractors. Then, the crossover probability CR and scaling factor F of the current iteration are calculated, and the mutation operator and crossover operator designed in this paper are used to generate a new multi-tractor cooperative operation path optimization. Then, the solution with the worst fitness is eliminated through the elite selection, and the three schemes with the best fitness in the last iteration are introduced. The scheme with better fitness is selected through the selection operation to replace the current scheme until the termination condition is reached. This scheme with the best current fitness is output as the optimal multi machine cooperative path planning scheme solved by the algorithm.

Experiment and Analysis
This paper took the John Deere 6B-1204 tractor commonly used in planting as the research object, and the total number of tractors was six. The working width of each tractor was W = 3.5 m, the minimum turning radius was r min = 5.5 m, the working speed was v w = 1 m·s −1 in the working area, the straight speed was v s = 5 m·s −1 in the non-working area, and the turning speed was v t = 2 m·s −1 . The peach orchard was rectangular, and the row spacing of peach trees was 4 × 5 m. As the peach tree distribution will affect the algorithm performance, this paper tests the path planning effect of the proposed model and the improved algorithm with different peach tree distributions. Figure 5 shows the distribution of peach trees in three rectangular standard peach orchards. Figure 5a shows the peach orchards with the rectangular distribution of peach trees. There were 20 rows of peach trees in the peach orchard, with 15 trees in each row. Figure 5b shows peach orchards with the approximately trapezoidal distribution of the peach trees, and Figure 5c shows peach orchards with an irregular distribution of peach trees. The experiment assumes that the tractor works in a peach orchard without obstacles and the driving process is continuous. The parameters of the original differential evolution algorithm DE are as follows: population size N = 100, crossover probability CR = 0.1, scaling coefficient F = 0.5, and the maximum number of iterations is 500. Equation (1) is used as the fitness function. The experimental environment is Windows 10 and Python 3.7.

Analysis of Experimental Results
The evaluation indexes of the model adopt the total turning time decrease rate (TTTDR) and total operating time decrease rate (TOTDR). The higher the value, the better the optimization effect of the algorithm. TTTDRn represents the reduction rate of total turning time between AEDE optimized and control group when n tractors work together. Similarly, TOTDRn represents the reduction rate of operation time after AEDE optimization and comparison algorithms when n unmanned tractors work together. The calculation formulas of TTTDRn and TOTDRn are shown in Equations (21) and (22).
In addition, we also use the following indexes to evaluate the algorithm performance. TTTOn is the total turning time of the n combined after control group optimization. TTTIn is the total turning time of the n combined after AEDE optimization. TOTOn is the total operating time of the n combined after control group optimization. TOTIn is the total operating time of the n combined after AEDE optimization.

Performance Comparison of Elite Selection Algorithm
Using traditional DE as the control group, we compared it with the AEDE proposed in this paper in terms of total turning time and total operating time, respectively. Both DE and AEDE are intelligent evolutionary algorithms; different peach tree distributions do not affect the solution process. Therefore, this comparison is based on the rectangular distribution of peach trees. The results are shown in Table 1. According to Table 1, the total turning time of AEDE method proposed in this paper was reduced by 3.34%, and the total operating time was reduced by 0.87%. It shows that after improving DE, the total turning time and total operating time are optimized to a certain extent. The parameters of the original differential evolution algorithm DE are as follows: population size N = 100, crossover probability CR = 0.1, scaling coefficient F = 0.5, and the maximum number of iterations is 500. Equation (1) is used as the fitness function. The experimental environment is Windows 10 and Python 3.7.

Analysis of Experimental Results
The evaluation indexes of the model adopt the total turning time decrease rate (TTTDR) and total operating time decrease rate (TOTDR). The higher the value, the better the optimization effect of the algorithm. TTTDR n represents the reduction rate of total turning time between AEDE optimized and control group when n tractors work together. Similarly, TOTDR n represents the reduction rate of operation time after AEDE optimization and comparison algorithms when n unmanned tractors work together. The calculation formulas of TTTDR n and TOTDR n are shown in Equations (21) and (22).
TOTDR n = TOTO n − TOTI n TOTO n × 100% In addition, we also use the following indexes to evaluate the algorithm performance. TTTO n is the total turning time of the n combined after control group optimization. TTTI n is the total turning time of the n combined after AEDE optimization. TOTO n is the total operating time of the n combined after control group optimization. TOTI n is the total operating time of the n combined after AEDE optimization.

Performance Comparison of Elite Selection Algorithm
Using traditional DE as the control group, we compared it with the AEDE proposed in this paper in terms of total turning time and total operating time, respectively. Both DE and AEDE are intelligent evolutionary algorithms; different peach tree distributions do not affect the solution process. Therefore, this comparison is based on the rectangular distribution of peach trees. The results are shown in Table 1. According to Table 1, the total turning time of AEDE method proposed in this paper was reduced by 3.34%, and the total operating time was reduced by 0.87%. It shows that after improving DE, the total turning time and total operating time are optimized to a certain extent. Block is a common multi-tractor scheduling scheme, which divides orchards into different areas and is operated by a single tractor. Using block operation as the control group, this paper compares it with AEDE in terms of total turning time and total operating time, respectively. The results are shown in Table 2. From Table 2, the total turning time of AEDE after optimization is reduced by more than 30%, and the total operating time is reduced by about 10%. This is because when working in block mode, the tractor needs to make frequent Ω-turns on the headland, which increases a lot of non-working time compared with the Π-Turn. AEDE avoids a lot of turning by optimizing the driving path of multiple tractors so as to reduce the total turning time and total operating time. In addition, AEDE can get good results whether the shape of the peach orchard is regular or irregular, which shows that AEDE is not sensitive to peach tree distribution and can be used to optimize peach orchard with complex peach tree distribution.

Comparison of Effective Operation Ability
Effective operation and average operation are important indexes to evaluate the cooperative operation efficiency of multiple tractors [37]. In a standard peach orchard, the effective operation refers to the ratio of the peach orchard area to the total operating time, and the average operation refers to the ratio of the effective operating capacity to the number of tractors. Figure 6 shows the changes of effective operating capacity and average effective operating capacity after block and AEDE optimization when the number of tractors working in rectangular orchards increases from two to six.

Comparison of Effective Operation Ability
Effective operation and average operation are important indexes to evalua cooperative operation efficiency of multiple tractors [37]. In a standard peach orcha effective operation refers to the ratio of the peach orchard area to the total operatin and the average operation refers to the ratio of the effective operating capacity number of tractors. Figure 6 shows the changes of effective operating capacity and a effective operating capacity after block and AEDE optimization when the num tractors working in rectangular orchards increases from two to six. In Figure 7, the effective operating capacity shows an increasing trend b increasing the number of tractors reduces the total operating time in the orchard. reduces the total turning time by optimizing the turning path of the tractor, which reduces the total operating time. When the number of tractor operations is tw effective operation capacity and average effective operation capacity calculated by operation are 0.0124 km 2 ·h −1 and 0.0062 km 2 ·h −1 , respectively. The effective ope capacity and average effective operation capacity calculated by AEDE are 0.0138 k and 0.0069 km·h −1 , respectively, which are 11.5% and 11.5% higher than respectively. When the number of tractors is six, the effective operation capaci average effective operation capacity calculated by AEDE are increased by 8.6% and respectively, compared with block. In Figure 7, the effective operating capacity shows an increasing trend because increasing the number of tractors reduces the total operating time in the orchard. AEDE reduces the total turning time by optimizing the turning path of the tractor, which in turn reduces the total operating time. When the number of tractor operations is two, the effective operation capacity and average effective operation capacity calculated by block operation are 0.0124 km 2 ·h −1 and 0.0062 km 2 ·h −1 , respectively. The effective operation capacity and average effective operation capacity calculated by AEDE are 0.0138 km 2 ·h −1 and 0.0069 km·h −1 , respectively, which are 11.5% and 11.5% higher than block, respectively. When the number of tractors is six, the effective operation capacity and average effective operation capacity calculated by AEDE are increased by 8.6% and 8.6%, respectively, compared with block. Effective operation and average operation are important indexes to evaluate th cooperative operation efficiency of multiple tractors [37]. In a standard peach orchard, th effective operation refers to the ratio of the peach orchard area to the total operating time and the average operation refers to the ratio of the effective operating capacity to th number of tractors. Figure 6 shows the changes of effective operating capacity and averag effective operating capacity after block and AEDE optimization when the number o tractors working in rectangular orchards increases from two to six. In Figure 7, the effective operating capacity shows an increasing trend becaus increasing the number of tractors reduces the total operating time in the orchard. AEDE reduces the total turning time by optimizing the turning path of the tractor, which in turn reduces the total operating time. When the number of tractor operations is two, th effective operation capacity and average effective operation capacity calculated by block operation are 0.0124 km 2 ·h −1 and 0.0062 km 2 ·h −1 , respectively. The effective operation capacity and average effective operation capacity calculated by AEDE are 0.0138 km 2 ·h − and 0.0069 km·h −1 , respectively, which are 11.5% and 11.5% higher than block respectively. When the number of tractors is six, the effective operation capacity and average effective operation capacity calculated by AEDE are increased by 8.6% and 8.6% respectively, compared with block.

Conclusions
In order to improve the operating efficiency of tractors in standard peach orchards, this paper studies a path optimization algorithm based on differential evolution for multiple tractors.
(1) In the rectangular peach orchard, the total turning time and total operating time in the process of management operation are reduced by 3.34% and 0.87%, respectively. Compared with the non-optimized differential evolution algorithm, it shows that the performance of the algorithm studied in this paper is better.
(2) Compared with block operation, the total turning time and total operating time in the process of management operation obtained by this algorithm are reduced by 37.37% and 9.47%, and the effect is better. It not only improves the operation efficiency of the tractors in the standard peach orchard, but also reduces the operating loss of the tractors.
However, when the tractors meet in the turning area, they will automatically avoid obstacles by adjusting the path. As the adjusted path is shorter than the length of the turning path and working path, the mutual interference of tractors in the turning area is ignored in this paper. The time consumption of tractors entering and leaving the peach orchard is also not considered in the model. In further work, we plan to use the time window to construct the turning constraint so as to avoid the situation that tractors meet in the turning area. At the same time, a more accurate path optimization method for multi-tractor cooperative work is realized by setting fixed starting points and ending points of tractors.