Research Article  Open Access
Spacecraft MultipleImpulse Trajectory Optimization Using Differential Evolution Algorithm with Combined Mutation Strategies and BoundaryHandling Schemes
Abstract
Since most spacecraft multipleimpulse 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 boundaryhandling schemes are used to simultaneously address the same infeasible trial vector. Two typical spacecraft multipleimpulse trajectory optimization problems are studied and optimized using the proposed DE method. The experimental results demonstrate that the proposed DE method 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.
1. 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 multipleimpulse 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 multipleimpulse 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 gradientbased methods [4–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 lowthrust trajectory optimization is less attractive due to the large number of variables required to describe lowthrust 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 multipleimpulse trajectory optimization problems which are called samecircle 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 multipleimpulse 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 multipleimpulse 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 dynamicsize 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 multipleimpulse 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–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–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 multipleimpulse 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 realworld 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 boundaryhandling schemes may have discrepant peculiarities and significantly affect the optimization performance. Xu and RahmatSamii [28] analysed the effect of different boundaryhandling schemes in the PSO method, and multiple boundaryhandling 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 boundaryhandling 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 multipleimpulse trajectory optimization problems, we studied not only multiple mutation strategies but also multiple boundaryhandling schemes and extended the traditional DE method to a modified variant with combined mutation strategies and boundaryhandling 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 boundaryhandling schemes when solving spacecraft multipleimpulse trajectory optimization problems. Furthermore, we successfully obtain the global optima for the samecircle 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.
2. Model Description of the Optimization Problems
2.1. SameCircle Rendezvous
The samecircle rendezvous problem is a multipleimpulse multiplerevolution 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 ( is the period of the circular orbit), was first presented by Prussing and Chiu [30] with a fourimpulse optimum solution. Colasurdo and Pastrone [31] and Prussing [32] further obtained a better fourimpulse optimum solution using the method based on Lawden’s theory and gradientbased 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 arein which , , , and are the timeapplying impulses, and are the characteristic velocities of the first two impulses, and and determine the direction of the first two impulses. The moduli and direction of the last two impulses, and , are determined using Lambert algorithm.
An optimization will minimize the propellant cost, which is equivalent to minimizing the total characteristic velocity. The objective function is
2.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 (MGA1DSM). 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 MGA1DSM 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 MGA1DSM model is able to design a better trajectory compared with the MGA model. However, there are more variables in the MGA1DSM model and the problem becomes much more difficult considering the large number of local optima.
In the MGA1DSM 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 MGA1DSM model are as follows:where is the departure time. Here, , , and 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 as follows:where is the velocity of planet in heliocentric ecliptic inertial reference frame (HEIRF), and are the velocities of spacecraft before and after flying by the planet in HEIRF, and is the planet gravitational constant.
Optimization also minimizes the total characteristic velocity:where is the escape velocity from Earth, represent the middle impulses between each two gravity assist times, and is the final impulse used to rendezvous with the target planet.
3. Differential Evolution Algorithm with Combined Mutation Strategies and BoundaryHandling 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.
3.1.1. Population Initialization
Given a function with dimensional variables, the population is initialized with a size of. The th vector of the population 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 is a random number uniformly distributed between 0 and 1.
3.1.2. Mutation Operation
After initialization, each individual vector , 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/currenttorand/1:”
“DE/randtobest/1:”
“DE/currenttobest/1:”
The indices , , , , and are five different integers, not equal to randomly selected from . We use 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 .
3.1.3. Crossover Operation
Crossover is applied after the mutation operation. Crossover gives the new individual vector , which is called a trial vector, the genes from both the target vector , and the mutant vector . Binomial crossover is more frequently used in the crossover operation:
In (16), is an integer randomly selected from . This ensures that at least one component in comes from , which maintains the efficiency of each crossover operation. Here, 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 in the classical DE method.
3.1.4. Boundary Handling
In most problems, especially realworld problems, the vector is restricted by the boundary constraint. Hence, the trial vector should be in the range . If or , boundary handling must be applied in order to address the th component’s presence in the range. The common boundaryhandling 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 , then the vector enters the next generation. Otherwise, is retained in the next generation.
3.2. 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/currenttorand/1” strategies. The other class focuses on the local search. Examples of this include the “DE/best/1,” “DE/best/2,” and “DE/currenttobest/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/currenttorand/1,” and “DE/currenttopbest/1.” The “DE/currenttopbest/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 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/currenttobest/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/currenttopbest/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.
3.3. Combined BoundaryHandling Schemes
When addressing problems with boundary constraints, boundary handling is an important operation if a trial vector exceeds the range. Because boundary handling essentially repeats the initialization process in order to replace the infeasible component, the boundaryhandling scheme determines the distribution of the new trial vector 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 boundaryhandling 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
Scheme 2 (“current_rand”). Consider
Scheme 3 (“reflect_rand”). Consider
Scheme 4 (“cut_off”). Consider
In (17)–(20), is the beforehandling component and is the afterhandling component. Figure 3 illustrates these four kinds of boundaryhandling schemes. And because the handling schemes of a component exceeding the minimum or maximum are similar, the situation of is considered as an example in the following illustration.
(a) “Whole_rand”
(b) “Current_rand”
(c) “Reflect_rand”
(d) “Cut_off”
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 . Then a new component is randomly created between and . 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 boundaryhandling schemes than that of methods that use a single scheme.
3.4. SelfAdaption of Parameters
The parameters and 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 selfadaption 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 withwhere is a random number uniformly distributed between 0.8 and 1, is the assemblage of all successful mutation factors at the current generation, and represents the power mean of all successful mutation factors, which can be calculated by
In (23), is the number of successful mutation factors.
The selfadaption scheme used for is similar to the one used for . The difference is that is a random number taken from the Gaussian distribution
The fundamental idea of this selfadaption 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 , chosen from the Gaussian distribution, concentrates the value. If the value for or exceeds the range , it will be recreated until the value is within the range. If and 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 , selected from the uniform distribution , maintain updating robustness. The parameters and do not have fixed values. Numerous experiments show that and are the best choices for complex multimodal problems.
3.5. Algorithm Flow of DE_CMSHBS
Figure 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 iterations no longer produce better results. If attaining the prescribed fitness evaluation number is the stop criterion, fewer evolutionary generations are needed because more fitness evaluation numbers in each generation are contained in the DE_CMSHBS algorithm than those contained in other algorithms.
4. Simulation and Comparison
4.1. Simulation of SameCircle Rendezvous Case
4.1.1. Result Presentation
Based on the samecircle 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.27 m/s, surpassing the best solution of Luo et al. [10] using PSASM method with an objective function value of . 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 primermagnitude 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 samecircle rendezvous problem.
4.1.2. Algorithm Comparison and Analysis
Except for the DE_CMSBHS method, we also tried eight singlecombined DEs in solving the samecircle 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 boundaryhandling 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 boundaryhandling scheme “whole_rand.” “DE/rand/1” + PBHS denotes the combination of the mutation strategy “DE/rand/1” and the combined boundaryhandling schemes. The population size 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.(3)When FEs = 10^{5}, all of the results converge. In this case, DE_CMSBHS, DE_CMS + “whole_rand,” DE_CMS + “reflect_rand,” and “DE/currenttopbest/1” + DE_CBHS converge to the theoretical optimum point 1256.27 km/s.
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.30 m/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 boundaryhandling schemes, overcomes this challenge and seldom converges to this local solution.
4.2. Deep Space Gravity Assist Maneuvers Case
4.2.1. Result Presentation
Based on the MGA1DSM 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 4.

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 website is X = [−779.046754, 3.259114, 0.525976, 0.380865, 167.378952, 424.028254, 53.289741, 589.766955, 2200.000000, 0.769483, 0.513289, 0.027418, 0.263985, 0.599985, 1.348780, 1.050000, 1.307303, 69.809014, −1.593737, −1.959525, −1.554988, −1.513462]^{T} with an objective function value of .
The best solution obtained by DE_CMSBHS method is X = [−780.148378, 3.274998, 0.530680, 0.382069, 168.487299, 423.998725, 53.306666, 589.771686, 2200.000000, 0.774442, 0.533183, 0.391920, 0.057816, 0.888942, 1.360604, 1.050000, 1.306803, 69.812308, −1.594186, −1.959564, −1.554776, −1.513431]^{T}, which slightly surpasses the current optimum solution with an objective function value of .
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.

4.2.2. 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.(1)The comparison shows that the results of the DE_CMSBHS method not only exceed the modified DE variants but also the DE_CMS and DE_CBHS method. It verifies that the DE_CMSBHS method has a stronger global search ability for this problem.(2)The standard deviation (Std) of the DE_CMSBHS method in Table 6 is less than those of the SADE, JADE, and two classic DE methods while it is greater than the MDE_pBX method. This suggests that the optimization stability of the DE_CMSBHS method should be further improved.(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 boundaryhandling schemes combination method does.
When optimizing realworld 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 + “current_rand” methods converge 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 Figure 9. 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 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 boundaryhandling 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 boundaryhandling schemes.
5. 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 boundaryhandling schemes is proposed. By applying the DE_CMSBHS method in the samecircle rendezvous problem and deep space gravity assist maneuvers problem, we successfully obtain the global optimum for the 400 km height samecircle 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 singlecombined DE methods and some other popular EAs. The simulation results indicate that by simultaneously utilizing multiple mutation strategies and boundaryhandling 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 boundaryhandling 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.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the Natural Science Foundation of China (no. 11402295), the Hunan Provincial Natural Science Foundation of China (no. 2015JJ3020), and the TransCentury Training Programme Foundation for the Talents by the State Education Commission (NCET130159).
References
 B. A. Conway, Spacecraft Trajectory Optimization, Cambridge University Press, New York, NY, USA, 2010.
 S. P. Hughes, L. M. Mailhe, and J. J. Guzman, “A comparison of trajectory optimization methods for the impulsive minimum fuel rendezvous problem,” Advances in the Astronautical Sciences, vol. 113, pp. 85–104, 2003. View at: Google Scholar
 Y. Z. Luo, J. Zhang, and G. J. Tang, “Survey of orbital dynamics and control of space rendezvous,” Chinese Journal of Aeronautics, vol. 27, no. 1, pp. 1–11, 2014. View at: Publisher Site  Google Scholar
 P. Gage, R. Braun, and I. Kroo, “Interplanetary trajectory optimization using a genetic algorithm,” Journal of the Astronautical Sciences, vol. 43, no. 1, pp. 59–76, 1995. View at: Google Scholar
 V. CoverstoneCarroll, “Nearoptimal lowthrust trajectories via microgenetic algorithms,” Journal of Guidance, Control, and Dynamics, vol. 20, no. 1, pp. 196–198, 1997. View at: Publisher Site  Google Scholar
 J. W. Hartmann, V. L. CoverstoneCarroll, and S. N. Williams, “Optimal interplanetary spacecraft trajectories via a Pareto genetic algorithm,” Journal of the Astronautical Sciences, vol. 46, no. 3, pp. 267–282, 1998. View at: Google Scholar  MathSciNet
 B. Dachwald and B. Wie, “Solar sail kinetic energy impactor trajectory optimization for an asteroiddeflection mission,” Journal of Spacecraft and Rockets, vol. 44, no. 4, pp. 755–764, 2007. View at: Publisher Site  Google Scholar
 M. R. Sentinella and L. Casalino, “Cooperative evolutionary algorithm for space trajectory optimization,” Celestial Mechanics & Dynamical Astronomy, vol. 105, no. 1–3, pp. 211–227, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 Y.Z. Luo, G.J. Tang, Y.J. Lei, and H.Y. Li, “Optimization of multipleimpulse, multiplerevolution, rendezvousphasing maneuvers,” Journal of Guidance, Control, and Dynamics, vol. 30, no. 4, pp. 946–952, 2007. View at: Publisher Site  Google Scholar
 Y.Z. Luo, J. Zhang, H.Y. Li, and G.J. Tang, “Interactive optimization approach for optimal impulsive rendezvous using primer vector and evolutionary algorithms,” Acta Astronautica, vol. 67, no. 34, pp. 396–405, 2010. View at: Publisher Site  Google Scholar
 M. Pontani, P. Ghosh, and B. A. Conway, “Particle swarm optimization of multipleburn rendezvous trajectories,” Journal of Guidance, Control, and Dynamics, vol. 35, no. 4, pp. 1192–1207, 2012. View at: Publisher Site  Google Scholar
 T. Vinkó and D. Izzo, “Global optimization heuristics and test problems for preliminary spacecraft trajectory design,” ACT Technical Report, 2008. View at: Google Scholar
 M. Vasile, E. Minisci, and M. Locatelli, “An inflationary differential evolution algorithm for space trajectory optimization,” IEEE Transactions on Evolutionary Computation, vol. 15, no. 2, pp. 267–281, 2011. View at: Publisher Site  Google Scholar
 A. Gad and O. Abdelkhalik, “Hidden genes genetic algorithm for multigravityassist trajectories optimization,” Journal of Spacecraft and Rockets, vol. 48, no. 4, pp. 629–641, 2011. View at: Publisher Site  Google Scholar
 O. Abdelkhalik and A. Gad, “Dynamicsize multiple populations genetic algorithm for multigravityassist trajectories optimization,” Journal of Guidance, Control, and Dynamics, vol. 35, no. 2, pp. 520–529, 2012. View at: Publisher Site  Google Scholar
 S. M. Islam, S. Das, S. Ghosh, S. Roy, and P. N. Suganthan, “An adaptive differential evolution algorithm with novel mutation and crossover strategies for global numerical optimization,” IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, vol. 42, no. 2, pp. 482–500, 2012. View at: Publisher Site  Google Scholar
 R. Storn and K. V. Price, “Differential evolution—a simple and efficient adaptive scheme for global optimization over continuous spaces,” Tech. Rep. TR95012, ICSI, Berkeley, Calif, USA, 1995. View at: Google Scholar
 R. Storn and K. Price, “Differential evolution—a simple and efficient heuristic for global optimization over continuous spaces,” Journal of Global Optimization, vol. 11, no. 4, pp. 341–359, 1997. View at: Publisher Site  Google Scholar  MathSciNet
 J. Brest, S. Greiner, B. Bošković, M. Mernik, and V. Žumer, “Selfadapting control parameters in differential evolution: a comparative study on numerical benchmark problems,” IEEE Transactions on Evolutionary Computation, vol. 10, no. 6, pp. 646–657, 2006. View at: Publisher Site  Google Scholar
 A. K. Qin, V. L. Huang, and P. N. Suganthan, “Differential evolution algorithm with strategy adaptation for global numerical optimization,” IEEE Transactions on Evolutionary Computation, vol. 13, no. 2, pp. 398–417, 2009. View at: Publisher Site  Google Scholar
 J. Zhang and A. C. Sanderson, “JADE: adaptive differential evolution with optional external archive,” IEEE Transaction on Evolutionary Computation, vol. 13, no. 5, pp. 945–958, 2009. View at: Publisher Site  Google Scholar
 F. Simeoni and L. Casalino, “Evolutionary optimization of interplanetary trajectories: improvements from initial diversification,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 225, no. 11, pp. 1277–1288, 2011. View at: Publisher Site  Google Scholar
 D. Izzo, V. M. Becerra, D. R. Myatt, S. J. Nasuto, and J. M. Bishop, “Search space pruning and global optimisation of multiple gravity assist spacecraft trajectories,” Journal of Global Optimization, vol. 38, no. 2, pp. 283–296, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 M. Vasile, E. Minisci, and M. Locatelli, “A dynamical system perspective on evolutionary heuristics applied to space trajectory optimization problems,” in Proceedings of the IEEE Congress on Evolutionary Computation (CEC '09), pp. 2340–2347, May 2009. View at: Publisher Site  Google Scholar
 B. Addis, A. Cassioli, M. Locatelli, and F. Schoen, “A global optimization method for the design of space trajectories,” Computational Optimization and Applications, vol. 48, no. 3, pp. 635–652, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 M. F. Tasgetiren, P. N. Suganthan, P. QuanKe, R. Mallipeddi, and S. Sarman, “An ensemble of differential evolution algorithms for constrained function optimization,” in Proceedings of the IEEE Congress on Evolutionary Computation (CEC '10), pp. 1–8, Barcelona, Spain, July 2010. View at: Publisher Site  Google Scholar
 S. M. Elsayed, R. A. Sarker, and D. L. Essam, “An improved selfadaptive differential evolution algorithm for optimization problems,” IEEE Transactions on Industrial Informatics, vol. 9, no. 1, pp. 89–99, 2013. View at: Publisher Site  Google Scholar
 S. Xu and Y. RahmatSamii, “Boundary conditions in particle swarm optimization revisited,” IEEE Transactions on Antennas and Propagation, vol. 55, no. 3 I, pp. 760–765, 2007. View at: Publisher Site  Google Scholar
 T. Huang and A. S. Mohan, “A hybrid boundary condition for robust particle swarm optimization,” IEEE Antennas and Wireless Propagation Letters, vol. 4, no. 1, pp. 112–117, 2005. View at: Publisher Site  Google Scholar
 J. E. Prussing and J.H. Chiu, “Optimal multipleimpulse timefixed rendezvous between circular orbits,” Journal of Guidance, Control, and Dynamics, vol. 9, no. 1, pp. 17–22, 1986. View at: Publisher Site  Google Scholar
 G. Colasurdo and D. Pastrone, “Indirect optimization method for impulsive transfers,” in AIAA/AAS Astrodynamics Conference, Scottsdale, Ariz, USA, 1994. View at: Publisher Site  Google Scholar
 J. E. Prussing, “A class of optimal twoimpulse rendezvous using multiplerevolution Lambert solutions,” Journal of the Astronautical Sciences, vol. 48, no. 23, pp. 131–148, 2000. View at: Google Scholar
 S. Das and P. N. Suganthan, “Differential evolution: a survey of the stateoftheart,” IEEE Transaction on Evolutionary Computation, vol. 15, no. 1, pp. 4–31, 2011. View at: Publisher Site  Google Scholar
 H. Y. Fan and J. Lampinen, “A trigonometric mutation operation to differential evolution,” Journal of Global Optimization, vol. 27, no. 1, pp. 105–129, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 D. F. Lawden, Optimal Trajectories for Space Navigation, Butterworths, London, UK, 1963. View at: MathSciNet
Copyright
Copyright © 2015 Yuehe Zhu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.