Research Article  Open Access
Hybrid Differential Evolution Optimisation for Earth Observation Satellite Scheduling with TimeDependent EarlinessTardiness Penalties
Abstract
We study the order acceptance and scheduling (OAS) problem with timedependent earlinesstardiness penalties in a single agile earth observation satellite environment where orders are defined by their release dates, available processing time windows ranging from earliest start date to deadline, processing times, due dates, sequencedependent setup times, and revenues. The objective is to maximise total revenue, where the revenue from an order is a piecewise linear function of its earliness and tardiness with reference to its due date. We formulate this problem as a mixed integer linear programming model and develop a novel hybrid differential evolution (DE) algorithm under selfadaptation framework to solve this problem. Compared with classical DE, hybrid DE employs two mutation operations, scaling factor adaptation and crossover probability adaptation. Computational tests indicate that the proposed algorithm outperforms classical DE in addition to two other variants of DE.
1. Introduction
Over the last decade, the order acceptance and scheduling (OAS) problem has been widely applied in maketoorder systems. In the existing literature, an order completed past the due date incurs a penalty. There is no reward or penalty for order completion before the promised date, for example, the printing and laminating company described by OĒ§uz et al. [1].
The justintime philosophy, which is attracting more attention currently, refers to the opinion that earliness, in addition to tardiness, should be discouraged. It has promoted the study of scheduling problems in which orders are preferred to be completed just at their respective due dates, and both early and tardy products are penalised. We indeed need to consider the penalty for the early completion of orders in an OAS system for some cases, which is the motivation in Section 2.
A variety of OAS problems under different settings and with various objectives have been reviewed by Slotnick [2]. The problem is classified into categories in terms of problem characteristics, objectives, and solution methods. There are also some papers to study a single machine scheduling problem with the objective of minimising the total tardiness [3, 4]. In the existing literature, we focus on reviewing studies on the OAS problem with timerelated penalties and sequencedependent setup times.
For the single machine total tardiness problem with sequencedependent setup times, Anghinolfi and Paolucci [5] proposed a novel discrete particle swarm optimisation (DPSO) approach to solve this nondeterministic polynomial time (NP) hard problem, and Tasgetiren et al. [6] provided a discrete differential evolution (DE) algorithm for the same problem. OĒ§uz et al. [1] studied the OAS problem, where orders were defined by their release dates, processing times, due dates, deadlines, sequencedependent setup times, revenues, and tardiness penalty weights. The authors assumed that the revenue of each accepted order is a function of its tardiness and deadline. They provided a mixed integer linear programming (MILP) formulation of the problem and developed three heuristic algorithms to solve this problem. Cesaret et al. [7] studied the OAS problem with the same setting as OĒ§uz et al. [1] and developed a tabu search algorithm. The proposed tabu search outperforms the three heuristics proposed by OĒ§uz et al. [1] in terms of solution quality. For the same problem setting, Lin and Ying [8] and Cheng et al. [9] developed an artificial bee colony algorithm and improved genetic algorithm. The deviation from the upper bound of the benchmark instances was further improved by the two algorithms.
There are some studies on the scheduling problem with earliness/tardiness penalties in a single machine environment. Hassin and Shani [10] studied several scheduling problems. They considered earliness and tardiness (E/T) penalties and developed polynomial time algorithms for special cases of these problems with nonexecution penalties. Chen et al. [11] formulated a mathematical model for production scheduling subject to multiple constraints, including sequencedependent setup costs, nonexecution penalties, and E/T penalties. The authors considered the E/T penalty weight as a special value (equal to one) to calculate. Baker [12] assumed that processing times follow normal distributions and due dates are decisions. The author developed a branch and bound algorithm to determine optimal solutions to minimise the total expected E/T costs and tested some heuristic procedures. Shabtay [13] analysed a model that integrates due date assignment and scheduling decisions. The objective was to minimise the total weighted earliness, tardiness, and due date assignment penalties.
To the best of our knowledge, the tardiness penalty and earliness penalty are seldom jointly considered for OAS problems in current research. In this study, we simultaneously consider the tardiness penalty and earliness penalty for the OAS problem with a sequencedependent setup time on a single agile earth observation satellite. Inspired by agile earth observation satellite scheduling, we formulated a general MILP model for this problem. The formulated general MILP model can be also applied in manufacturing system, logistics management, and maketoorder systems, as oversubscribed scheduling problem also exists in these fields.
The remainder of the paper is organised as follows. We illustrate the motivation for this paper in Section 2. In Section 3, we formally define the OAS problem with timedependent earlinesstardiness penalties. In Section 4, we present the mathematical model for this problem. In Section 5, we provide details for the proposed algorithm for solving the problem. In Section 6, we present the experimental studies. In Section 7, we provide the conclusion.
2. Motivation
With the expansion of earth observation satellites (EOSs) application, such as earth resource exploration, forest fire warning, volcano eruption discovery, and earthquake surveillance, EOSs become more limited and expensive relative to numerous orders submitted by users from different departments. The satellite operation firms that can operate on a maketoorder basis realize benefits in terms of timely image availability and decreased risk to imaging randomness. On the other hand, limitations on satellite resource force these firms to be selective on user orders. Users giving orders typically quote a due image quality and may penalise the firm for notgood quality image. This translates into reduced revenue and even loss of users.
Our motivation originates from the agile earth observation satellite (AEOS) scheduling problem. Compared with the traditional earth observation satellite, the AEOS is mobile along three axes (roll, pitch, and yaw). This mobility considerably increases the length of the time windows for imaging targets. The pitch angle required for the satellite to image the target changes in the time window. When the satellite flies over the target, the pitch angle required is zero and the best quality image is obtained. If the target is observed earlier or later than this time, a relatively larger pitch angle is required and a deterioration of image quality is incurred, as illustrated in Figure 1, which is partially taken from LemaĆ®tre et al. [14].
Consequently, the order processing time window for this target is from the earliest start time to the latest end time. Earth observation satellites are scarce and oversubscribed resources. A satellite can at most image one target at any time moment and cannot image all the userrequested targets because of its limited capacity and ability. When imaging targets in reality, the pitch angles required always differ from one target to another; therefore, the setup times depend on the imaging sequence of the targets. Thus, simultaneous decisions have to be made to choose which targets are to be imaged and fix the appropriate start time of imaging. We formulate this problem as an OAS problem, considering the timedependent earliness and tardiness penalties with the objective of maximising total revenue. In the remainder of this paper, we refer to this problem as the OAS with timedependent earlinesstardiness penalties (OASET) problem.
3. Problem Definition
In a single AEOS environment, a set of jobs are initialised at the beginning of the planning horizon. The satellite can process one job at most at any time moment without preemption. There are no precedence constraints among the jobs, and the production system does not break down. The description of the job is shown in Figure 2.
For convenience, the notations used in the remainder of the paper are summarised in the Notations.
Note that the setup time depends on the processing sequence of two consecutive jobs. Thus is not necessarily equal to .
The revenue of a job is defined as a piecewise linear function that depends on the earliest start date, processing time, due date, and deadline as illustrated in Figure 3. The job is not processed yet in region A and the job is not completed yet in region B. In both regions, we obtain no revenue. In regions C and D, the revenue gained decreases linearly with its earliness or tardiness, respectively. In region E, the manufacturer does not accept the job, which is completed beyond its deadline, and no revenue is gained. If the job is completed exactly at the due date, the maximum revenue is obtained.
Given a sequence of the selected job set , we can calculate the completion time of each job . Let be the tardiness of job in sequence . We can calculate the tardiness of job using the formula . The earliness of the job is . Thus, the revenue of job is calculated as . Consequently, the total revenue gained from the sequence is . The OASET problem is to determine the set and sequence when is maximum.
This problem reduces to the OAS problem studied by OĒ§uz et al. [1] if the unit earliness penalty cost is set to zero and the release date equals the earliest start date for all orders. The OASET problem generalises the OAS problem, which has been proved to be NPhard; thus it is still of NPhard complexity.
4. The Mathematical Model
In this section, we formulate the OASET problem as an MILP model. We apply two binary variables, and , to denote order acceptance and sequencing decisions, respectively:
MILP is
In this model, we add the definition of earliness for each job with constraints set (12), (13), and (14) and then update constraint (15) to calculate the revenue the firm can gain when job is accepted and incurs tardiness or earliness . Additionally, the procession time window is characterised by the earliest start date and deadline; therefore, we use the earliest start date in constraint (6).
5. Hybrid Differential Evolutionary Algorithm
Differential evolution (DE), which was first proposed by Storn and Price [15], has become arguably one of the most efficient stochastic realparameter optimisation algorithms. DE is very simple and easy to implement. Despite this, DE has exhibited better performance than several other algorithms [16]. Based on classical DE, many variants have greatly improved performance for a variety of applications. We also propose a hybrid DE (HDE) to solve the problem in Section 5.5.
5.1. Parameter Vector
For an OASET problem with all jobs, an dimensional real vector represents an individual, where . The value of for job represents the proportion of the part time before job is completed in the completion interval . The completion time of job is calculated by . In such a representation, the completion time of each job is located in the completion time window .
5.2. Graph Representation
Given a parameter vector, the completion time of each job is fixed. Therefore, the revenue that might be gained from each job is also fixed. To maximise the total revenue, we need to select a subset of jobs that satisfy the sequencedependent setup time constraints. First, we introduce a directed acyclic graph of the problem under a given parameter vector. The procedure is described as follows.
Step 1. Compute the procession start time and completion time of each job under a given parameter vector.
Step 2. Sort the jobs in ascending order by their procession start times, suppose the sorted sequence of all the jobs is , and calculate the penalised revenue for each job .
Step 3. First, add nodes to graph , where each node represents job , then as illustrated in the mathematical model, add two dummy nodes including one denoted by ā0ā that represents the source node and another one denoted by āā that represents the sink node to .
Step 4. For each job in and for each job after job in , if the sequencedependent setup time constraint is satisfied, that is, , then add a directed weighted to .
Step 5. For each job in , add a directed weighted and a directed weighted zero to .
Consider the following example: suppose there are four incoming jobs and the data for the problem are as follows: , , , , , , and . Suppose the parameter vector , and we can obtain the start and completion times as and , respectively; then the penalised revenues of the four jobs are calculated as , , , . Figure 4 illustrates the jobs sorted in ascending order by their start times. We can observe that Jobs 1 and 2 are penalised by earliness, whereas Jobs 3 and 4 are penalised by tardiness. The constructed directed graph from the jobs is illustrated in Figure 5.
5.3. GraphBased Fitness Evaluation
Based on the directed acyclic graph constructed under a given parameter vector, selecting a subset of jobs maximising the total net revenue is equivalent to determining the longest path from the source node ā0ā to sink node āā [17]. The polynomial time algorithm exists to solve the longest path problem in a directed acyclic graph. We modify the BellmanFord shortest path algorithm to determine the longest path in the graph. The nodes, except the source node and sink node in the path, represent the selected jobs in the solution. In Figure 5, the longest path, which is illustrated by bold arcs, is 0245 with total weight .
5.4. Classical DE
The classical DE algorithm proposed by Storn and Price [15] can be summarised as follows. The initial population comprising individuals , , is randomly generated, where and is the dimension of the problem. DE generates a new population through mutation, crossover, and selection operators.
(1) Mutation. This operation generates new target vectors according towhere , , and are different integers randomly chosen from and different from the vector index . is a constant factor and is located in , according to Storn and Price [15].
(2) Crossover. This operation forms the trial vectors according towhere is the constant crossover rate, which is often set to 0.9, and is a randomly chosen index from 1 to to ensure that the trial vector does not duplicate .
(3) Selection. This operation determines which parent vector is selected for the next generation from vector and trial vector according to
5.5. Hybrid DE Method
Our proposed HDE combines several components of other algorithms under a selfadaptation framework. The DE/currentto with archive in Adaptive Differential Evolution with Optional External Archive (JADE) [18] and the widely used DE/rand are combined under selfadaptation strategies. Considering the plateaus in the fitness landscape, the neighbourhood search in Neighbourhood Search Differential Evolution (NSDE) [19], which adopts two random distribution generators for tuning the scaling factor, is applied in our algorithm to help the algorithm escape from plateaus.
5.5.1. Mutation
In addition to (20), there are several other commonly used mutation operations [20]:
We apply selfadaptation strategies [21], and two mutation operations are selected as candidates. One is DE/ran/1, as shown in (20), and the other is DE/currentto proposed by Zhang and Sanderson [18], given aswhere is randomly chosen as one of the top individuals in the current population and is randomly chosen from the union of the current population and archive. The archive is initiated to be empty. After each generation, the parent solutions that fail in the selection of (22) are added to the archive. We randomly remove some individuals in the archive to keep the size within if the size exceeds [19].
The mutation process applied here is the same as Selfadaptive Differential Evolution (SaDE) [21], except that (25) used in SaDE is replaced by (26). The trial vector is produced aswhere is initialised as 0.5. After each generation, the number of trial vectors successfully entering the next generation generated by (20) and (26) is recorded as and . The number of trial vectors that fail entering the next generation generated by (20) and (26) is recorded as and . The four numbers are aggregated within a specific number of generations and the probability is updated as
5.5.2. Scaling Factor Adaptation
For our problem described in Section 3, we can determine that there are possibly many plateaus in the fitness landscape; that is, there are many vectors with the same objective value. As the population converges, the difference between two vectors tends to be small. To escape from a plateau, various scaling factors are required. We apply the selfadaptive neighbourhood search strategy [19] to set the value of scaling factor:where is a randomly generated number according to a normal distribution of mean 0.5 and standard deviation 0.3, denotes a Cauchy random number generator with scale parameter , and is initialised as 0.5 and updated using (28), as does SaDE.
5.5.3. Crossover Probability Adaptation
The crossover probability for each individual at generation is generated using the procedure in JADE. In JADE, is generated according to a normal distribution of mean and standard deviation 0.1,and truncated to , where is initialised as 0.5 and at each generation updated according towhere is the set of all successful ās at generation and is the usual arithmetic mean of . We set as a constant value 0.1.
6. Experimental Studies
6.1. Data Generation
Five datasets are generated with 20, 40, 60, 80, and 100 jobs. In each dataset, 10 instances are generated. Each instance is denoted as , where is the number of jobs and is the index of the instance in the dataset with jobs. BĆ¼lbĆ¼l et al. [22] proposed a method for generating test instances for the single machine earliness/tardiness scheduling problem. We apply this method and add the deadline and setup time generation methods to generate our instances. On each instance, the processing time is a random number generated uniformly in ; the due date is generated in , where is the total processing time of all jobs; is the tardiness factor, which is a coarse measure of the proportion of jobs that might be tardy; determines the interval length of the due date; and the term is added to to prevent the possible initial infeasibility of . Earliest start date is generated uniformly in , setup time is generated in , deadline is set as , and revenue is a randomly generated number in . These parameters are listed in Table 1.

6.2. Experimental Results
Our proposed HDE was implemented in and run on a Windows 7 laptop with Intel i5 2.4āGHz CPU, 4āGB RAM. The initial population size of all the problems was set to 40. The algorithm stopped if it reached a maximum of 2,000 generations. On each instance, the algorithm ran 10 times () independently. The maximum, minimum, and average fitness in the 10 runs were recorded. The results of DE [15], NSDE [19], and JADE [18] were also recorded. The results of datasets with different jobs are given in Tables 2ā5. The average relative percentage deviation () was calculated as follows:where and are the fitness function values generated by the HDE algorithm and the reference algorithms (including DE, NSDE, and JADE algorithms), respectively, for each run and is the number of runs. Similarly, and denote the minimum and maximum of the relative percentage deviations, respectively, from the reference fitness function values taken.




6.3. Performance Comparisons
For small problems with 20 jobs, Tables 2 and 3 show that the HDE algorithm performed better than DE and NSDE and slightly better than JADE on all instances. Compared with DE and NSDE, the average relative percentage improvements for HDE are approximately 3 generally, of which the best is 5.08 and the worst is 0.80. Compared with JADE, the average relative percentage improvements for HDE are all less than 1, of which the best is 0.88 and the worst is 0.03. There is no significant difference between the performance of DE and NSDE, but for small problems with 40 jobs, the maximum and average fitness function values of NSDE and JADE are much better than those for DE, as shown in Table 2. HDE achieves the maximum fitness among the four algorithms. In Table 3, the average relative percentage improvement of HDE compared with DE is appropriately 10 times better than that compared with NSDE and JADE, and JADE and NSDE achieve similar performance. Compared with DE, the average relative percentage deviations for HDE are almost 14 generally, and the minimum relative percentage deviations are nearly 18. Compared with NSDE and JADE, the average relative percentage deviations for HDE are approximately 1.5 generally, and the minimum relative percentage deviations are approximately 1.6.
For medium sized problems with 60 and 80 jobs, compared with NSDE and JADE, our proposed HDE provides significantly better results on instances 80_1, 80_3, 80_4, 80_6, 80_7, 80_8, 80_9, and 80_10, as shown in Tables 1 and 4. On those instances, even the worst fitness function value determined by HDE is better than the best fitness function values determined by NSDE and JADE. On other instances, the maximum and average fitness function values of HDE are better than those of NSDE and JADE. Both NSDE and JADE outperform classical DE; in particular, the average relative percentage improvements by HDE compared with DE are almost five times better than the ones compared to NSDE and JADE. NSDE performs slightly better than JADE in terms of all three indicators.
For large problems with 100 jobs, the results are shown in Tables 1 and 5. Our proposed HDE significantly outperforms DE, NSDE, and JADE. In 10 independent runs, the worst fitness value of HDE is even better than the best fitness value of DE, NSDE, and JADE on instances 100_1, 100_2, 100_6, 100_7, 100_8, 100_9, and 100_10. On the other three instances 100_3, 100_4, and 100_5, the maximum and average fitness of 10 runs are better than those of DE, NSDE, and JADE (the details of result are in http://pan.baidu.com/s/1eSf8jzs).
7. Conclusions
To the best of our knowledge, this is the first attempt to study the OAS problem with timedependent earlinesstardiness penalties and sequencedependent setup times, which is inspired by AEOS scheduling. We proposed hybrid DE algorithm under selfadaptation framework, including selfadaptation strategy for mutation operator, the selfadaptive neighbourhood search strategy for scaling factor, and crossover probability adaptation. The significance of selfadaptation framework is proved by the experimental results. Furthermore, we presented a novel method to evaluate the fitness for this problem. Firstly a directed acyclic graph (DAG) is constructed under a given parameter vector, and then the fitness evaluation phase is transformed to finding the longest path from the source node to sink node in the DAG. The graphbased fitness evaluation provides an efficient method for the similar problems. Computational results show that our proposed algorithm outperforms classical DE in addition to two variants: NSDE and JADE. The proposed method not only solves the OASET problem but also suggests an effective way to apply a continuous algorithm to combinatorial optimisation problems. In regard to future research, the dynamic OAS problem catches our attention for considering the orders come stochastically during the planning horizon.
Notations
Indices:  Job index, 
:  Set of all jobs 
:  Set of selected jobs, 
:  Sorted sequence of jobs in 
:  Sequence of jobs in 
:  Release date of job 
:  Earliest start date of job 
:  Processing time of job 
:  Due date of job , ; is not necessarily in the middle point of and 
:  Deadline of job , 
:  Maximum revenue of job 
:  Unit tardiness penalty cost of job ,āā 
:  Unit earliness penalty cost of job ,āā 
:  Sequencedependent setup time for job being processed immediately after job () 
:  Procession time window of job , 
:  Completion time window of job , 
:  Procession start time of job 
:  Completion time of job . 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research is supported by the National Natural Science Foundation of China under Grant nos. 61603400, 61473301, 61203180, 71501179, and 71501180.
References
 C. Oǧuz, F. Sibel Salman, and Z. Bilgintürk Yalçin, āOrder acceptance and scheduling decisions in maketoorder systems,ā International Journal of Production Economics, vol. 125, no. 1, pp. 200ā211, 2010. View at: Publisher Site  Google Scholar
 S. A. Slotnick, āOrder acceptance and scheduling: a taxonomy and review,ā European Journal of Operational Research, vol. 212, no. 1, pp. 1ā11, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 P. Guo, W. Cheng, and Y. Wang, āA general variable neighborhood search for singlemachine total tardiness scheduling problem with stepdeteriorating jobs,ā Journal of Industrial and Management Optimization, vol. 10, no. 4, pp. 1071ā1090, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Y. K. Lin and C. S. Chong, āA tabu search algorithm to minimize total weighted tardiness for the job shop scheduling problem,ā Journal of Industrial and Management Optimization, vol. 12, no. 2, pp. 703ā717, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 D. Anghinolfi and M. Paolucci, āA new discrete particle swarm optimization approach for the singlemachine total weighted tardiness scheduling problem with sequencedependent setup times,ā European Journal of Operational Research, vol. 193, no. 1, pp. 73ā85, 2009. View at: Publisher Site  Google Scholar
 M. F. Tasgetiren, Q.K. Pan, and Y.C. Liang, āA discrete differential evolution algorithm for the single machine total weighted tardiness problem with sequence dependent setup times,ā Computers and Operations Research, vol. 36, no. 6, pp. 1900ā1915, 2009. View at: Publisher Site  Google Scholar
 B. Cesaret, C. Oǧuz, and F. Sibel Salman, āA tabu search algorithm for order acceptance and scheduling,ā Computers and Operations Research, vol. 39, no. 6, pp. 1197ā1205, 2012. View at: Publisher Site  Google Scholar
 S.W. Lin and K.C. Ying, āIncreasing the total net revenue for single machine order acceptance and scheduling problems using an artificial bee colony algorithm,ā Journal of the Operational Research Society, vol. 64, no. 2, pp. 293ā311, 2013. View at: Publisher Site  Google Scholar
 C. Cheng, Z. Yang, L. Xing, and Y. Tan, āAn improved genetic algorithm with local search for order acceptance and scheduling problems,ā in Proceedings of the 2013 IEEE Symposium on Computational Intelligence in Production and Logistics Systems, CIPLS 2013  2013 IEEE Symposium Series on Computational Intelligence, SSCI 2013, pp. 115ā122, sgp, April 2013. View at: Publisher Site  Google Scholar
 R. Hassin and M. Shani, āMachine scheduling with earliness, tardiness and nonexecution penalties,ā Computers & Operations Research, vol. 32, no. 3, pp. 683ā705, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 Y.W. Chen, Y.Z. Lu, and G.K. Yang, āHybrid evolutionary algorithm with marriage of genetic algorithm and extremal optimization for production scheduling,ā International Journal of Advanced Manufacturing Technology, vol. 36, no. 910, pp. 959ā968, 2008. View at: Publisher Site  Google Scholar
 K. R. Baker, āMinimizing earliness and tardiness costs in stochastic scheduling,ā European Journal of Operational Research, vol. 236, no. 2, pp. 445ā452, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 D. Shabtay, āOptimal restricted due date assignment in scheduling,ā European Journal of Operational Research, vol. 252, no. 1, pp. 79ā89, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 M. Lemaître, G. Verfaillie, F. Jouhaud, J.M. Lachiver, and N. Bataille, āSelecting and scheduling observations of agile satellites,ā Aerospace Science and Technology, vol. 6, no. 5, pp. 367ā381, 2002. View at: Publisher Site  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
 S. Das and P. N. Suganthan, āDifferential evolution: A survey of the stateoftheart,ā IEEE Transactions on Evolutionary Computation, vol. 15, no. 1, pp. 4ā31, 2011. View at: Publisher Site  Google Scholar
 E. M. Arkin and E. B. Silverberg, āScheduling jobs with fixed start and end times,ā Discrete Applied Mathematics. The Journal of Combinatorial Algorithms, Informatics and Computational Sciences, vol. 18, no. 1, pp. 1ā8, 1987. View at: Publisher Site  Google Scholar  MathSciNet
 J. Zhang and A. C. Sanderson, āJADE: Adaptive differential evolution with optional external archive,ā IEEE Transactions on Evolutionary Computation, vol. 13, no. 5, pp. 945ā958, 2009. View at: Publisher Site  Google Scholar
 Z. Yang, K. Tang, and X. Yao, āSelfadaptive differential evolution with neighborhood search,ā in Proceedings of the 2008 IEEE Congress on Evolutionary Computation, CEC 2008, pp. 1110ā1116, chn, June 2008. View at: Publisher Site  Google Scholar
 K. V. Price, R. M. Storn, and J. A. Lampinen, Differential Evolution: A Practical Approach to Global Optimization, Springer, Berlin, Germany, 2005. View at: MathSciNet
 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
 K. Bülbül, P. Kaminsky, and C. Yano, āPreemption in single machine earliness/tardiness scheduling,ā Journal of Scheduling, vol. 10, no. 45, pp. 271ā292, 2007. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2017 Guoliang Li 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.