Spacecraft Multiple-Impulse Trajectory Optimization Using Differential Evolution Algorithm with Combined Mutation Strategies and Boundary-Handling Schemes

Since most spacecraft multiple-impulse trajectory optimization problems are complex multimodal problems with boundary constraint, finding the global optimal solution based on the traditional differential evolution (DE) algorithms becomes so difficult due to the deception of many local optima and the probable existence of a bias towards suboptimal solution. In order to overcome this issue and enhance the global searching ability, an improved DE algorithm with combined mutation strategies and boundaryhandling schemes is proposed. In the first stage, multiple mutation strategies are utilized, and each strategy creates a mutant vector. In the second stage, multiple boundary-handling schemes are used to simultaneously address the same infeasible trial vector. Two typical spacecraft multiple-impulse trajectory optimization problems are studied and optimized using the proposed DE method. The experimental results demonstrate that the proposed DEmethod efficiently overcomes the problem created by the convergence to a local optimum and obtains the global optimum with a higher reliability and convergence rate compared with some other popular evolutionary methods.


Introduction
The subject of spacecraft trajectory optimization has a long history [1].In its early stage, spacecraft trajectory was primarily optimized using analytical theory and gradientbased optimization algorithms.Hughes et al. [2] concluded that these optimization approaches were efficient in some simple cases but did not work well on more complex ones, such as spacecraft multiple-impulse trajectory optimization problems, because most of these cases present complex multimodal peculiarities, which easily converge to the local optimum and make the search of the global optimum difficult.
Since the 1990s, evolutionary algorithms (EAs) have rapidly emerged.Varieties of EAs have been used for spacecraft multiple-impulse trajectory optimization problems.These include genetic algorithm (GA), simulated annealing (SA), differential evolution (DE), and particle swarm optimization (PSO) [3].Initial applications of EAs to space trajectory optimization mainly employed GAs in conjunction with gradient-based methods [4][5][6].A combination of artificial neural networks with EAs has also been applied to solar sail trajectories [7].Sentinella and Casalino [8] mentioned that the use of EAs for low-thrust trajectory optimization is less attractive due to the large number of variables required to describe low-thrust trajectories with sufficient accuracy, while EAs are better suited to the optimization of impulse trajectories which can be described by a limited number of variables and the number of function evaluations that are required to obtain the optimal solution is usually acceptable.
In this paper, we focus our research on two typical multiple-impulse trajectory optimization problems which are called same-circle rendezvous problem and deep space gravity assist maneuvers problem, respectively.Some researchers have studied these problems based on EAs recently.Luo et al. [9,10] applied an improved SA method to optimize the multiple-impulse rendezvous problem.Pontani et al. [11] used a PSO method to solve similar cases.Vinkó and Izzo [12] formulated a model of gravity assist using deep space maneuvers and applied several kinds of EAs to optimize the multiple-impulse gravity assist problem.Vasile et al. [13] furthered the work of Vinkó and Izzo [12] by applying a hybrid algorithm with the DE method and monotonic basin hopping to the original model.Gad and Abdelkhalik [14,15] presented a software tool based on their developed hidden genes GA and dynamic-size multiple population GA method to solve the deep space gravity assist maneuvers problems.These employed EAs tend to have a stronger searching ability than traditional approaches but cannot yet converge to the global optimum when optimizing multimodal spacecraft multiple-impulse trajectory problems.Therefore, a more powerful and efficient EA is necessary in order to search the global optimum when solving these problems.
The DE algorithm, as a simple yet efficient optimizer with fewer parameters, is currently one of the most widely used EAs.First proposed by Storn and Price [17,18], it has received a substantial amount of attention.Many researchers have studied DE and proposed numerous notable variants which have been verified on a series of numerical problems [16,[19][20][21].Thanks to its high efficiency, DE has been applied to solve optimization problems in many fields, including the spacecraft trajectory optimization [8,12,13,[22][23][24][25].However, because only one search operator is used in most modified DE variants, they may have excellent optimization performance on unimodal problems but their capability to solve the multimodal problems is not so outstanding, especially for the complex multiple-impulse trajectory optimization problem.Hence, some researchers applied multiple search operators in the algorithm.Tasgetiren et al. [26] developed an ensemble DE in such a way that each individual was assigned to one of two distinct mutation strategies or a variable parameter search.Elsayed et al. [27] proposed an improved differential evolution algorithm that uses a mix of different mutation operators to solve the benchmark and some real-world optimization problems.
Moreover, for the problems with boundary constraint, the boundary handling also plays an important role in the evolutionary process except the mutation operation.Different boundary-handling schemes may have discrepant peculiarities and significantly affect the optimization performance.Xu and Rahmat-Samii [28] analysed the effect of different boundary-handling schemes in the PSO method, and multiple boundary-handling schemes were applied by Huang and Mohan [29] in order to improve the algorithm's robustness.However, little literature has reported the usage of different boundary-handling schemes on DE's multiple search operators.
In our research, in order to further enhance the optimization performance of DE algorithm and then to obtain global optimal solutions of the spacecraft multiple-impulse trajectory optimization problems, we studied not only multiple mutation strategies but also multiple boundary-handling schemes and extended the traditional DE method to a modified variant with combined mutation strategies and boundary-handling schemes.By comparing our modified DE method with some other popular evolutionary methods in amount of simulations, we found that the global searching ability of DE algorithm could be efficiently improved by simultaneously combining multiple mutation strategies and boundary-handling schemes when solving spacecraft multiple-impulse trajectory optimization problems.Furthermore, we successfully obtain the global optima for the same-circle rendezvous case and deep space gravity assist maneuvers case using the proposed DE method.The reminder of the paper is organized as follows.In Section 2, the model Description of the Optimization problems is given out.The modification mechanism of the improved DE method is detailed in Section 3. In Section 4, the simulation results are presented first and the comparison and analysis of the employed algorithms followed then.Finally, conclusions are presented in Section 5.

Model Description of the Optimization Problems
2.1.Same-Circle Rendezvous.The same-circle rendezvous problem is a multiple-impulse multiple-revolution rendezvous problem.In this, a chaser spacecraft and a target spacecraft have the same circular orbit.The chaser has an initial separation angle of 180 ∘ behind the target, which is shown in Figure 1.The chaser must rendezvous with the target at a prescribed time using impulse maneuvers.This problem, using a transfer time of 2.3 *  ( is the period of the circular orbit), was first presented by Prussing and Chiu [30] with a four-impulse optimum solution.Colasurdo and Pastrone [31] and Prussing [32] further obtained a better four-impulse optimum solution using the method based on Lawden's theory and gradient-based optimization algorithms.Recently, Luo et al. [10] and Pontani et al. [11] produced similar results using the global convergence ability of the evolutionary algorithms, in which a parallel simulated annealing that utilizes a simplex method was presented in Luo et al. [10] and a PSO method was applied in Pontani et al. [11].
The variables of this problem are  An optimization will minimize the propellant cost, which is equivalent to minimizing the total characteristic velocity.The objective function is (2)

Deep Space Gravity Assist Maneuvers.
Gravity assist maneuvers can be divided into two categories, which are called multiple gravity assist (MGA) maneuvers and multiple gravity assist using deep space maneuvers (MGA-1DSM).
Their corresponding optimization models have been well formulated by Vinkó and Izzo [12].The MGA manoeuver represents the interplanetary trajectory of a spacecraft equipped with chemical propulsion that can only thrust during its planetocentric phases.This type of problem is easier to optimize because it has fewer design variables.The MGA-1DSM refers to the model in which a spacecraft is able to thrust its engine once at any time between each trajectory leg.Generally, for the same gravity assist sequence, using the MGA-1DSM model is able to design a better trajectory compared with the MGA model.However, there are more variables in the MGA-1DSM model and the problem becomes much more difficult considering the large number of local optima.
In the MGA-1DSM model, the trajectory between two gravity assist maneuvers is divided into two parts, as can be seen in Figure 2. Kepler algorithm is used to propagate the first part, and the second part should be solved by Lambert algorithm.The variables of the MGA-1DSM model are as follows: where  0 is the departure time.Here,  ∞ , , and V define the heliocentric direction of the departure hyperbolic velocity,   is the time rendezvous with each planet,   is a coefficient between 0 and 1 that determines the flight time of the Kepler part, and   and   are used to calculate each velocity increment ΔV GA as follows: ΔV out     ΔV out     = cos i  + sin   sin i  + cos   sin i  , where V pla is the velocity of planet in heliocentric ecliptic inertial reference frame (HEIRF), V in and V out are the velocities of spacecraft before and after flying by the planet in HEIRF, and  pla is the planet gravitational constant.Optimization also minimizes the total characteristic velocity: where  ∞ is the escape velocity from Earth, Δ DSMi represent the middle impulses between each two gravity assist times, and Δ final is the final impulse used to rendezvous with the target planet.

Differential Evolution Algorithm with Combined Mutation Strategies and Boundary-Handling Schemes (DE_CMSBHS)
3.1.Classic DE Algorithm.The DE algorithm is a simple real parameter optimization algorithm [33].It proceeds through a simple cycle of stages, which include mutation, crossover, boundary handling, and selection.

Population Initialization.
Given a function (X) with -dimensional variables, the population is initialized with a size of.The th vector of the population X  for each generation can be expressed as The search space is constrained by the prescribed minimum and maximum bounds: The th component of the th vector in ( 6) is initialized as follows: where rand  [0, 1] is a random number uniformly distributed between 0 and 1.

Mutation Operation.
After initialization, each individual vector X  , called a target vector, employs the mutation operation in order to create a corresponding mutant vector according to the prescribed mutation strategy.The following mutation strategies are frequently employed in the classic DE. "DE/rand/1:" "DE/rand/2:" "DE/best/1:" "DE/best/2:" "DE/current-to-rand/1:" "DE/rand-to-best/1:" "DE/current-to-best/1:" The indices 1, 2, 3, 4, and 5 are five different integers, not equal to  randomly selected from [1, 𝑁𝑝].We use X best to denote the best individual vector of the current generation, and   is the mutation factor of the th vector.The mutation factor is an important parameter in the algorithm.In the classical DE method,   is a fixed value within the range [0, 1].

Crossover Operation.
Crossover is applied after the mutation operation.Crossover gives the new individual vector U  = [ 1 ,  2 ,  3 , . . .,   ]  , which is called a trial vector, the genes from both the target vector X  , and the mutant vector V  .Binomial crossover is more frequently used in the crossover operation: In ( 16),  rand is an integer randomly selected from [1, 𝐷].This ensures that at least one component in U  comes from V  , which maintains the efficiency of each crossover operation.Here, CR  is the crossover probability of the th vector.This crossover probability is another important parameter in the algorithm and has a fixed value within the range [0, 1] in the classical DE method.

Boundary Handling.
In most problems, especially realworld problems, the vector X  is restricted by the boundary constraint.Hence, the trial vector U  should be in the range [X min , X max ].If   <  min  or   >  max  , boundary handling must be applied in order to address the th component's presence in the range.The common boundary-handling scheme recreates a new component, as in (4), to replace the original one.
3.1.5.Selection.After the previous steps have been completed, each trial vector must be evaluated.A "greedy" strategy is applied in the selection operation.This means that as long as (U  ) ≤ (X  ), then the vector U  enters the next generation.Otherwise, X  is retained in the next generation.

Combined Mutation Strategies.
Mutation is the most important step in a DE algorithm.Equations ( 9)-( 15) are some efficient mutation strategies frequently used in various DE methods.The effectiveness of some other strategies, such as the triangular mutation strategy proposed by Fan and Lampinen [34], have also been verified.In general, all mutation strategies can be divided into two classes.One group emphasizes the global search, similar to the "DE/rand/1, " "DE/rand/2, " and "DE/current-to-rand/1" strategies.The other class focuses on the local search.Examples of this include the "DE/best/1, " "DE/best/2, " and "DE/current-tobest/1" strategies.
For complex optimization problems, especially multimodal problems, using one of these strategies alone cannot ensure that the algorithm's global and local searches are effective and efficient.Qin et al. [20] employed multiple mutation strategies in the SADE method.A strategy is chosen based on its probability, which is updated according to the success ratios of the past generations.This approach is selfadaptive to a certain extent.However, only one of the available strategies is used in each mutation operation.These strategies are not used to their full potential, and choosing a strategy based on its probability is an unreliable approach.For a given target vector in a current generation, the strategy with a higher success probability may not create a better trial vector than one with a lower success probability.
In this case, in order to fully utilize the advantages of multiple mutation strategies, combined mutation strategies are applied in DE CMSBHS to simultaneously improve the global and local searches.Testing demonstrates that the combination of the following four strategies produces a better performance: "DE/rand/1, " "DE/rand/2, " "DE/currentto-rand/1, " and "DE/current-to-pbest/1. " The "DE/current-topbest/1" strategy is used in the JADE algorithm, which was proposed by Zhang and Sanderson [21].The pbest vector is randomly chosen from the top p% of individuals in the population.According to the conclusions in Zhang and Sanderson [21], the parameter p in the range [5,20] performs better.We used p = 20 in the TP SDE algorithm in order to enhance its global search ability.
The first three global searching strategies can help maintain the population's diversity to some degree; the fourth local searching strategy can improve the algorithm's search speed and accuracy.Even when the fourth strategy converges to a local optimum, the first three strategies increase the possibility of escaping the local optimum in order to evolve to the better solution.Their relationship can be regarded as complementation.However, if the fourth strategy is replaced with "DE/best/1" or "DE/current-to-best/1, " the algorithm's optimization performance significantly worsens.Because these two "best" strategies are too greedy, the first three strategies will be overwhelmed by the "best" strategy and become useless.Once the algorithm converges to the local optimum, it cannot progress further.Thus, we apply the less greedy strategy, "DE/current-to-pbest/1. " On one hand, this strategy can decrease the probability of premature convergence.On the other hand, the first three strategies will be used more frequently in the evolution procedure in order to maintain the population's diversity.The comparisons made in Section 4 demonstrate that the optimization performance of combined mutation strategies is better than methods that only use a single strategy.

Combined Boundary-Handling Schemes.
When addressing problems with boundary constraints, boundary handling is an important operation if a trial vector U  exceeds the range.Because boundary handling essentially repeats the initialization process in order to replace the infeasible component, the boundary-handling scheme determines the distribution of the new trial vector U  in the design space.Different schemes can significantly affect the final result.
Currently, only one scheme, which is the same as the one in (8), is used to operate boundary handling in most evolutionary algorithms.Randomly creating a new component with this scheme does not appear to be sufficiently efficient.
Although an individual with a component that exceeds the range is not a feasible solution, it is still valuable because it contains information about the mutation.In order to fully utilize this information about the current and infeasible individuals, multiple boundary-handling schemes are applied in the DE CMSBHS process.Testing demonstrates that the combination of the following four schemes performs better.
Scheme 1 ("whole rand").Consider if (  <  min  or   >  max  ) , Scheme 4 ("cut off ").Consider In ( 17)-( 20),   is the before-handling component and   is the after-handling component.Figure 3 illustrates these four kinds of boundary-handling schemes.And because the handling schemes of a component exceeding the minimum or maximum are similar, the situation of   >  max  is considered as an example in the following illustration.
As is shown in Figure 3, "whole rand" is the retaining scheme used to ensure the potential selection of any value within the given range.The "current rand" scheme randomly creates a new component between the current component and the maximum, which efficiently utilizes the information of the current individual.The "reflect rand" scheme first creates an    .   , in which the   are symmetric with respect to  max  .Then a new component is randomly created between    and  max  .This scheme efficiently uses the information of infeasible individuals.The "cut off " scheme directly replaces the component with the maximum, which is extremely efficient if some components of the optimum solution are on the boundary.The comparisons made in Section 4 also demonstrate a better optimization performance of combined boundary-handling schemes than that of methods that use a single scheme.
u ij and u  ij are symmetric with respect to x maxj (c) "Reflect rand"

Self-Adaption of Parameters.
The parameters  and CR in the classic DE algorithm are usually determined by testing or experience.They are always fixed values in any evolution procedure.However, the efficiency and reliability of this situation seems insufficient.For this reason, many studies on the selfadaption of parameters exist.The self-adaption of parameters scheme applied in this paper refers to the one proposed by Islam et al. [16].The parameter  adapts as follows: In (21),   is a random number selected from a Cauchy distribution with mean and standard deviation 0.1.Here,   is updated in every generation in accordance with where   is a random number uniformly distributed between 0.8 and 1,   is the assemblage of all successful mutation factors at the current generation, and mean(  ) represents the for (j = 0; j < 4; j++) Create U j by crossover V 1 using "DE/rand/1" V 2 using "DE/rand/2" using "DE/current-to-rand/1" using "DE/current-to-pbest/1" Create mutant vector V 3 Create mutant vector V 4 power mean of all successful mutation factors, which can be calculated by mean In (23), |  | is the number of successful mutation factors.The self-adaption scheme used for CR is similar to the one used for .The difference is that CR  is a random number taken from the Gaussian distribution CR  = Gaussian (CR  , 0.1) , The fundamental idea of this self-adaption scheme is that the parameters of the next generation are guided according to the parameters of all successful scale factors in the current generation.The parameter  selected from the Cauchy distribution maintains the population's diversity, and CR, chosen from the Gaussian distribution, concentrates the value.If the value for  or CR exceeds the range (0, 1], it will be recreated until the value is within the range.If   and CR  are 0.5 and 0.75, respectively, at the first generation, then they are adapted according to the situation of the following generations.The values for   and  CR , selected from the uniform distribution [0.8, 1], maintain updating robustness.The parameters  and  do not have fixed values.Numerous experiments show that  = 2 and  = 1.5 are the best choices for complex multimodal problems.4 presents the DE CMSHBS algorithm flow.The stop criterion can be set according to the problem's need.In general, the algorithm stops when the prescribed fitness evaluation number approaches or converges to a plateau in which successive

Result Presentation.
Based on the same-circle rendezvous model that is formulated in Section 2, the circular orbit is set to 400 km in height, which is the same as the simulation condition in Luo et al. [10].The boundaries of each variable are established in Table 1.
The best solution obtained by the DE CMSBHS method is X = [0.000000,0.207320, 0.792680, 1.000000, 392.546673, 0.500000, 0.547843, 235.587940, 0.500000, 0.142865] T with an objective function value of 1256.27m/s, surpassing the best solution of Luo et al. [10] using PSASM method with an objective function value of  = 1256.32m/s.Each impulse time and characteristic velocity optimized by DE CMSBHS and PSASM method are listed in Table 2.The chaser's rendezvous trajectory optimized by the DE CMSBHS method is given out in Figure 5.
Figure 6 illustrated the primer-magnitude time history of the obtained best solution.Apparently it satisfies Lawden's conditions, which is a necessary condition for an optimal transfer orbit [35].It indicates that DE CMSBHS has successfully obtained the optimum solution of the same-circle rendezvous problem.

Algorithm Comparison and Analysis.
Except for the DE CMSBHS method, we also tried eight single-combined DEs in solving the same-circle rendezvous problem to compare with the PSASM method that is used by Luo et al. [10], including four kinds of DEs with combined mutation strategies and four kinds of DEs with combined boundary-handling schemes.For convenience, the singlecombined DEs are abbreviated as follows.If we use DE CMS + "whole rand" and "DE/rand/1" + DE CBHS, for instance, DE CMS + "whole rand" represents the combination of the combined mutation strategies with the boundary-handling scheme "whole rand." "DE/rand/1" + PBHS denotes the combination of the mutation strategy "DE/rand/1" and the combined boundary-handling schemes.The population size NP is uniformly selected to be 100.We carry out 30 independent trials for each method and the results are listed in Table 3.The best solutions of each fitness evaluation numbers (FEs) are in boldface.
Comparing the results in Table 3, several conclusions can be drawn.
(1) In general, the results of all the combined DE methods are better than the PSASM method [10].Furthermore, the DE CMSBHS surpasses all other combined DEs.
(2) When the FEs are set to be 3 × 10 4 and 5 × 10 4 , the mean and standard deviation of the DE CMSBHS are less than the corresponding values of the combined DEs and PSASM.This indicates that the DE CMSBHS has a higher optimization speed and is more stable for this problem.
It is worth mentioning that there is a local optimum that also satisfies Lawden's conditions.Specifically, the solution is X = [0.000000,0.345262, 0.654738, 1.000000, 533.945461, 0.500000, 0.015312, 191.202864, 0.500000, 0.682373] T , and the objective function value is 1450.30m/s.Our experiment shows that results with larger standard deviations in Table 3 converge to this local optimum several times during the 30 trials.However, the DE CMSBHS, which combines both multiple mutation strategies and boundary-handling schemes, overcomes this challenge and seldom converges to this local solution.4.

Deep Space Gravity Assist Maneuvers
This case has been solved by Vasile et al. [13] and Gad and Abdelkhalik [14].The objective function value of their optimum solution was 8.3889 km/s and 8.3850 km/s, respectively.Solutions optimized by some other researchers are listed on the ESA's official website: http://www.esa.int/gsp/ACT/inf/op/globopt.Currently, the optimum solution on the    Figure 7 presents the optimized trajectory from departure to rendezvousing with Saturn.Each velocity increment from impulse or gravity assist is shown in Table 5.From Table 5 we can find that most characteristic velocity is donated by the initial departure and the final rendezvous.With the help of gravity assist, few velocity increments of the impulses are needed in the middle transfer.

Algorithm Comparison and Analysis.
Islam et al. [16] optimized the problem using several DE variants and their results are listed in Table 6.As a comparison for our results that optimized by DE CMSBHS method, the population size  was uniformly selected to be 100.Fifty independent trials were performed.Moreover, a comparison between the DE CMSBHS, four kinds of DE CMS, and four kinds of DE CBHS methods are displayed in Table 7.
From Tables 6 and 7, we can make the following observations.
(  (3) In general, the results of all the DE CMS exceed those of the DE CBHS method.This can be seen in Table 7.It could be argued that the mutation strategies combination more efficiently improves the DE algorithm than the boundary-handling schemes combination method does.
When optimizing real-world problems, we prefer to obtain an optimum solution.Hence, a further comparison is carried out and the algorithms do not stop until it converges to a plateau, in which successive iterations no longer produce better results.The best results for the compared algorithms in fifty independent trials are listed in Table 8.
It is apparent that the best solutions optimized by MDE pBX, DE CMS + "current rand, " and DE CMSBHS methods surpass the published optimum solution and the DE CMSBHS's solution is slightly better than other two methods.Figure 8 displays the convergence numbers used to determine the optimum solution of the MDE pBX, DE CMS + "current rand, " and DE CMSBHS methods in the 50 trials.The DE CMSBHS method converges to the optimum solution in 14 of the trials whereas the MDE pBX and DE CMS +  8 repeatedly converges to this local optimum.Further analysis of this local optimum demonstrates that the second component,  ∞ , was selected to be the minimum, which was 3 km/s.Because this component is directly added to the objective function, it seems to be more sensitive to the fitness value than the other components do.Also, because the DE CMS + "cut -off " method always replaces infeasible components with the lower bound (LB) or upper bound (UB), it clearly increases the possibility of obtaining this local optimum.The DE CMSBHS method, which simultaneously uses multiple boundary-handling schemes, seldom converges to this local optimum.This indicates that the global search ability of the DE CMSBHS method is efficiently improved by combining multiple boundary-handling schemes.

Conclusion
In order to improve the optimization performance of DE algorithms when solving multimodal spacecraft multipleimpulse trajectory optimization problems, a modified DE variant called DE CMSBHS which combines four mutation strategies and four boundary-handling schemes is proposed.By applying the DE CMSBHS method in the same-circle rendezvous problem and deep space gravity assist maneuvers problem, we successfully obtain the global optimum for the 400 km height same-circle rendezvous case and find an optimal solution for the Cassini mission which is currently the best solution as far as I know.Furthermore, we compared the DE CMSBHS method with eight kinds of single-combined DE methods and some other popular EAs.The simulation results indicate that by simultaneously utilizing multiple mutation strategies and boundary-handling schemes DE CMSBHS efficiently avoid the issue of converging to the local optimum and tend to the global optimum with a higher optimization speed and success rate.
However, the four mutation strategies and four boundary-handling schemes applied in the DE CMSBHS method may not produce the most efficient combination for other problems.Future research will focus on a method for determining the most efficient combination according to different problems and will extend the DE CMSBHS to multiobjective spacecraft trajectory optimization problems.

Figure 3 :
Figure 3: Illustration of the four kinds of boundary-handling schemes.

Figure 5 :Figure 6 :
Figure 5: The trajectory of same-circle rendezvous optimized by DE CMSBHS. Case

4. 2 . 1 .
Result Presentation.Based on the MGA-1DSM model that is formulated in Section 2, we study and optimize the classic Cassini mission.The flight sequence of Cassini mission is Earth-Venus-Venus-Earth-Jupiter-Saturn with 22 variables.The boundaries are established in Table

Figure 7 :
Figure 7: The trajectory of Cassini mission optimized by DE CMSBHS.

Figure 8 :Figure 9 :
Figure 8: Comparison of convergence numbers between the DE CMSBHS and other two methods.
in which  1 ,  2 ,  3 , and  4 are the time-applying impulses, Δ 1 and Δ 2 are the characteristic velocities of the first two impulses, and  1 ,  1 and  2 ,  2 determine the direction of the first two impulses.The moduli and direction of the last two impulses, Δ 3 and Δ 4 , are determined using Lambert algorithm.

Table 1 :
Design boundaries of same-circle rendezvous model.

Table 2 :
Impulse time and characteristic velocity of the best solution.

Table 3 :
Results of same-circle rendezvous model.

Table 4 :
Design boundaries of Cassini mission.

Table 5 :
The best result of each velocity increment from impulse/gravity assist.T , which slightly surpasses the current optimum solution with an objective function value of  = 8.3829 km/s.

Table 6 :
[16]lts of Cassini mission solved by some DE variants[16]and optimized by DE CMSBHS.

Table 7 :
Results of Cassini mission optimized by the DE CMS, DE CBHS, and DE CMSBHS.

Table 8 :
Best results of Cassini mission optimized by the compared algorithms.currentrand"methodsconverge only two and three times, respectively.This indicates that the DE CMSBHS method has a higher global optimization success rate.The convergence history of the three best results is shown in Figure9.The graph lets us see that the DE CMSBHS method converges faster than other two methods.In addition, there is a local optimum causing strong interference.The components are X = [−805.733015,3.000000, 0.616195, 0.384658, 195.117030, 422.971282, 53.293530, 589.769333, 2200.000000,0.113728, 0.514959, 0.047143, 0.013736, 0.026443, 1.270316, 1.050000, 1.307191, 69.809127, −1.616293, −1.959523, −1.554919, −1.513431] T , and the objective function value is 8.61 km/s.The DE CMS + "cutoff " method in Table "