Abstract
In dynamic multiobjective optimization problems, the environmental parameters change over time, which makes the true pareto fronts shifted. So far, most works of research on dynamic multiobjective optimization methods have concentrated on detecting the changed environment and triggering the population based optimization methods so as to track the moving pareto fronts over time. Yet, in many realworld applications, it is not necessary to find the optimal nondominant solutions in each dynamic environment. To solve this weakness, a novel method called robust paretooptimal solution over time is proposed. It is in fact to replace the optimal pareto front at each timevarying moment with the series of robust paretooptimal solutions. This means that each robust solution can fit for more than one timevarying moment. Two metrics, including the average survival time and average robust generational distance, are present to measure the robustness of the robust pareto solution set. Another contribution is to construct the algorithm framework searching for robust paretooptimal solutions over time based on the survival time. Experimental results indicate that this definition is a more practical and timesaving method of addressing dynamic multiobjective optimization problems changing over time.
1. Introduction
In many practical fields, such as engineering design, scientific computing, social economy, and network communication, there exist a large number of complex optimization problems. Particularly, many optimization problems contain multiple objective functions and dynamic parameters that make the objective functions change over time. Moreover, the number of objective functions and constraints may also vary from time to time. We call multiobjective optimization problems with above uncertain factors as dynamic multiobjective optimization problems (DMOPs). In this paper, we focus on the DMOPs with continuously changed dynamic parameters. Suppose is the objective vector in DMOPs and is the dynamic parameters depending on . The aim of DMOPs is to find the pareto front () approximating to the true pareto front of as soon as possible for all of the dynamic environment.
Aiming at tracking the moving true pareto fronts over time, dynamic multiobjective evolutionary optimization algorithms (DMOEA) were proposed [1–12]. The universal framework of DMOEA is presented in Algorithm 1. First of all, we need to accurately judge whether the environment has changed. It is the basic premise of using evolutionary optimization methods to respond to the environmental changes. The most common changedetection approach is to reevaluate the detectors. The detectors can be the current best solutions, a memorybased subpopulation, or a feasible subpopulation [5].

Once the various environment parameters happened, the new evolutionary optimization process was triggered. Many evolutionary algorithms with good performances on the static multiobjective optimization problems have been introduced into the DMOPs one after another, such as genetic algorithm [6], particle swarm algorithm [7], differential evolution algorithm [8], quantum immune clonal coevolutionary algorithm [9], and memetic algorithm [10]. As we know, the goal of the static optimization problems is to make the population gradually converge to optimal nondominant solution. The diversity of the population becomes weakened during the optimization methods. How to improve and maintain the evolutionary algorithm’s ability to adapt to the various environment is a major challenge faced in the dynamic evolutionary computation. In recent years, various methods have been used to improve the diversity of the population. Chen et al. [11] extended additional objectives to deal with the DMOPs. Individual diversity is used as an additional objective to provide the historical information. Zhou et al. [12] proposed a population prediction strategy to improve the DMOEA’s performance when the new environment is detected. The reinitial population was formed by a center point and a manifold. The track of center points was preserved to train autoregressive model so as to predict the center point in the new environment. A new coevolutionary paradigm [2] combining competition with cooperation was proposed to track the true pareto front in dynamic environment.
The conventional methods are mostly to trigger the multiobjective optimization process after detecting the change and then finish the evolution process as the following new environment occurs. Yet, this is impractical in many realworld optimization problems due to the following reasons. These methods are not suitable for rapidly changing environments, in which the environmental parameters vary quickly or frequently. Moreover, it is difficult to find the satisfied pareto fronts before detecting the change of fitness landscape because the methods are extremely timeconsuming.
To address the above concern, several methods to find robust pareto front of multiobjective optimization problems with noise were presented by [13–17]. The emphasis is to seek an insensitive robust pareto front instead of the global optimal pareto front. The detailed definition about the robustness is illustrated in Section 2. For dynamic scalar optimization problems, robust optimization over time (ROOT) has been defined clearly by Yu et al. [18]. The task for ROOT is to find a sequence of robust solutions over time intervals. They have acceptable qualities and are relatively insensitive to the dynamic environment. By ROOT, the uncertainties in the parameter space and their cumulative effect on objective space are considered simultaneously. Furthermore, Jin et al. [19] gave a framework of ROOT which consists of a populationbased optimization algorithm, the database, a fitness approximator, and a fitness predictor. A solution’s robustness over time is estimated by both its past and its future performances. Subsequently, Fu et al. [20] provided a feasible algorithm in order to find robust optimal solutions over time and gave a detailed definition of the survival time and the average fitness. The robust solutions are expected to have longer survival time or larger average fitness. Though ROOT is easy to be realized and computed, it only fits for dynamic scalar optimization problems.
For DMOPs, the detailed robust definition over time is still an open issue. There were seldom relevant literatures about this area. We will introduce this idea about ROOT into DMOPs. The main contribution in this paper is to propose a novel concept on robust paretooptima over time (RPOOT) to DMOPs that search for robust paretooptimal solution set for all dynamic environments. Subsequently, the populationbased multiobjective evolutionary algorithm is introduced to find RPOOT in terms of the nondomination solutions’ robust performance. The robustness is measured by the survival time derived from the robust index given by Deb and Gupta [14]. We believe it is a more practical way of addressing continuously changed DMOPs.
The remainder of the paper is structured as follows. Section 2 presents a brief introduction to research on existing robust optimization methods and analyzes the existing problems in detail. Section 3 presents a class of DMOPs and describes formally the definition of RPOOT. In Section 4, the robustness and the performance metrics are defined for RPOOT. Furthermore, a populationbased evolutionary algorithm to find robust paretooptimal solutions over time is presented. Section 5 provides a brief overview of the existing benchmark functions for DMOPs and the experimental result for RPOOT. Conclusions and future work are presented in Section 6.
2. Related Works
In order to solve multiobjective optimization problems with uncontrollable variations, Li et al. [13] presented a robust multiobjective genetic algorithm by considering two objective functions: the fitness value , which measures a solution’s performance by a combined objective, and the robustness index . They investigated the tradeoff between the convergence and robustness of the nondominant solutions. Furthermore, Li proposed the OuterInner optimization structure. The outer subproblem was to simultaneously minimize the fitness value and to maximize the robustness index. The inner subproblem calculates the radius , which represents a solution’s robustness.
Deb and Gupta [14] extended an existing approach that finds robust solutions for singleobjective optimization problems to MOPs with dynamic parameters. They defined the mean effective objective functions instead of the original objective functions.
Consider a multiobjective optimization problem as follows:
In order to avoid obtaining the global optimal solutions which are quite sensitive to such variable perturbation in their vicinity, the following two approaches are defined for robust optimization by Deb and Gupta [14].
Definition 1 (multiobjective robust solution of type I (MORS1)). A solution is called a multiobjective robust solution of type I if it is the global feasible paretooptimal solution to the following multiobjective minimization problem (suppose a neighborhood of a solution x is ):where is defined as follows:
Definition 2 (multiobjective robust solution of type II (MORS2)). A solution is called a multiobjective robust solution of type II if it is the global feasible paretooptimal solution to the following multiobjective minimization problem:
MORS1 replaces the original objective function with the effective objective function . Where is the mean of th objective function values in the vicinity of . In Definition 2, the original objective functions need to be optimized. At the same time, the feasible solutions must satisfy the constraint; that is, the objective function values among neighboring solutions are limited to a userdefined threshold . Subsequently, the solution’s robustness is judged. Two kinds of performances can not be analyzed at the same time. Both definitions for robustness took the variable perturbation into account. However, this is not DMOPs in fact.
Furthermore, Barrico and Antunes [15] defined the degree of robustness based on the solutions’ behavior in their neighborhood in the decision space. The degree of robustness was also used to evaluate the solutions’ behavior in the neighborhood of the reference scenario in the space of the objective functions’ coefficients [16]. The weakness of above definition about the degree of robust is that the paretooptimal set has to be known in advance. To solve this problem, Cromvik et al. [17] put forward a new definition for robustness index and introduced the utility function to convert multiobjective optimization problem into an approximation for a single decision maker. In this method, the robust index cannot be used as an objective during the optimization.
In a word, the above methods can solve the multiobjective optimization problems with perturbation in the decision space or the space of the objective function coefficients. However, only a robust pareto front meeting all multiobjective optimization problems with the disturbance is found. In this paper, we will discuss the robust solution set to a class of DMOPs with changing parameters.
3. The Definition of RPOOT
In this paper, we focus on the DMOPs, in which the environmental parameters continuously change over time and keep stationary between two timevarying moments. Obviously, the true pareto fronts shift from time to time. In other words, the objective functions depending on are deterministic during each changing stage. Hence, this kind of DMOPs with continuous parameters can be discretized into a series of multiobjective optimization problems (MOPs) at each timevarying moment. Namely, the pareto fronts at each timevarying moment are regarded as the basis of optimization.
Without loss of generality, dynamic multiobjective optimization problem is defined as follows:where stands for the decision variable vector, is the decision space, and is the objective space. consists of objective functions . is the timedepending parameter vector. represents the time. Suppose is the frequency for dynamic environments. There may be different time steps. DMOPs are divided into MOPs, denoted by .
Based on the definition of robust solutions’ index given by Deb and Gupta [14] and the description of ROOT presented by Yu et al. [18], we propose a novel concept, called robust paretooptimal over time (RPOOT) for DMOPs. It consists of a sequence of robust paretooptimal sets, denoted by . Each robust paretooptimal solution , may not approximate to the true pareto front at each period between changes but is relatively closer to all true pareto fronts during at least two consecutive time intervals.
Two hypotheses for RPOOT are that the environmental parameters change over time with stationary periods between changes and that the robust solution fits for the consecutive changes during . With and , we obtain the problemwhere . The dynamics parameter is a random variable obeying a certain distribution such as a Gaussian distribution or a Uniform distribution. Let be the probability density function of . The nondominant solutions’ performances are measured by the following two definitions.
Definition 3 (the average fitness). Consider
Definition 4 (the variance of fitness). Consider
In essence, measures the average performance of each objective function within the time interval . is the degree of the performance influenced by changing the timevarying environment.
Definition 5 (robustness). A solution is called a robust pareto solution if it is nondominant individuals satisfied to the following DMOPs:
In fact, all solutions composed of robust paretooptimal solution sets must be compromise. Subsequently, the nondominant relationship is defined based on the expected fitness vector over time. Define if , and , is satisfied. and are compromise, if , and , . Consequently, we define robust paretooptimal solutions as . Obviously, RPOOT takes both the approximation and the robustness into account.
In Figure 1, the parameters vary in different periods, which map the uncertain area in the objective space. Obviously, and . and are nondominated. However, is not robust due to its large variance. As shown in Figure 2, the true paretooptimal sets for three consecutive timevarying periods are , , plotted by the dash lines. plotted by real line in Figure 2 is the robust paretooptimal solution. Obviously, is not the best paretooptimal set during this period but the satisfied pareto front over time.
4. The Robust Solutions’ Performance
In this paper, a robust paretooptimal solution set needs not only to approximate to the true pareto front during each stage to the largest extent but also to satisfy the requirements of robustness. We would like to make a clear distinction between the definition of robustness for solutions and the approximation of each solution.
4.1. The Robustness
For DMOPs, the task of the traditional optimization methods is to find the paretooptimal solutions after detecting the new environment. This is timeconsuming. If the paretooptimal front can meet the requirements of more than one kind of environment in a certain accuracy, the cost for search will be less. Consequently, the robustness has two means. One is the insensitivity to the fluctuation parameters. Deb and Gupta [14] has presented a general definition based on neighborhood perturbation. The other is the survival time, which reflects how many consecutive changed environments this solution can fit for. Based on the above two aspects, corresponding metric called the survival time is proposed to measure the robustness in DMOPs.
Suppose is a nondomination solution at time . The robustness of is defined by maximum survival time starting from time when all fitness values of from to belong to neighborhood of .
Define the survival time of solution’s robustness as ;
stands for the survival time of any optimal solution in robust paretooptimal set . is the userdefined threshold, which is the key parameter having a direct impact on . The larger means the tolerance to the varieties of the fitness is better. It makes larger. operator is used to measure the distance between the current fitness values and the future predictive fitness values. Here, the Euclidean norm is used. But any other suitable norm can also be adopted. is the approximated fitness instead of the real fitness value by a predictor [18]. Because the timedepending parameters vary randomly, each nondominated solution’s real fitness values for each timevarying environment cannot be all evaluated accurately. So we need a prediction method to approximate the fitness in the future dynamic environments. As shown in Figure 3, if the robust paretooptimal solutions at time are satisfied during the consecutive timevarying moments from to , the fitnesses are all constricted to the neighborhood of .
4.2. The Average Fitness Function
Generally, the integral part of (7) is not easy to be calculated since there is little knowledge to get the accurate probability density function. So we take the following average performance as the robustness of solution over the considered time interval.
We define as a robust paretooptimal solution over time if it is the global feasible paretooptimal solution to the following multiobjective optimization problem:where stands for the average fitness value of th objective during the consecutive timevarying moments from to . Subsequently, measures the average performance of in each objective during the time interval . The smaller average fitness value means that is more approximate to the true pareto fronts of corresponding timevarying environments.
In a word, the average survival time reflects the robustness of on the time space; that is, how many timevarying environments it fits for. The average fitness value measures how the paretooptimal fronts approximate to the true pareto fronts in the objective space during these dynamic stages. Consequently, a robust paretooptimal solution must be the one with the minimum average fitness values.
4.3. The PopulationBased Optimization Method
Taking above metrics as the objectives, a novel framework to solve RPOOT problems is constructed for DMOPs. In the timevarying moments, many populationbased evolutionary algorithms may be adopted to find a set of robust paretooptimal fronts over time, denoted by .
It is worth noting that we must take into account the future performance of all objectives during the calculation of robust survival time and robustness performance. We assume that future performance can be estimated by a database and a predictor. The database is used to store historical data, and the task of the predictor is to estimate a solution’s future performances [19]. A framework of populationbased optimization algorithm for RPOOT is presented as following Algorithm 2. By this method, the robust pareto solution could not only approximate to the true pareto front as close as possible but also fit for more than one dynamic environment.

4.4. The Measurement of Algorithm Performance
It is important to measure the performance of the RPOOT algorithm. We should not only consider the robustness of the robust paretooptimal solutions but also consider the accuracy of the solutions to the true pareto front. On the one hand, the average survival time measures the solutions’ robustness on the time scale. On the other hand, the average robust metric measures the approximation to the true pareto front in the objective space.
4.4.1. The Average Survival Time
In DMOPs, the robustness of the robust optimal pareto solutions is measured by the average survival time. Moreover, the algorithms will be compared across the whole time period . Based on the robustness defined in (11), the robustness of the algorithm performance is defined as follows:
is the robustness survival time of the nondomination solutions obtained from the algorithm during the timeline. Obviously, the longer the average survival time is, the better the robustness of solutions is. The robustness of optimal parato front on the time scale can be reflected by the average survival time. Moreover, depends on the threshold . Therefore, a more exhaustive analysis is necessary for the robust optimal pareto solutions under different .
4.4.2. The Average Robust Generational Distance
This metric reflects the quality of the robust nondominant solution sets. The general distance (GD) [21] indicates how close the obtained PFs are to the true pareto front in multiobjective optimization problems. Furthermore, the inverted generational distance (IGD) [21] is used to assess the approximation performance of the algorithms. IGD measures both the diversity and the convergence of the population. In our experimental studies, robust generational distance (RGD) and robust inverted generational distance (RIGD) indicate the average distance between each robust optimal pareto front and the true pareto fronts within its survival time. Suppose is the size of robust nondominant pareto solution set during the whole time interval. The RGD and RIGD metrics are defined as follows:where stands for the time step of robust nondominant solution set. is their survival time. and . is a set of uniformly distributed optimal solutions in the true at ; is the solutions obtained at . is the distance between and . is the distance between and . and are the cardinalities of and . In our experiments, we select 100 evenly distributed solutions in .
5. Analysis of the Experimental Results
In this section, eight dynamic multiobjective benchmark functions are adopted in the experiments. Simulation results and further analysis on solutions’ performance are conducted in Section 5.2.
5.1. Benchmark Functions
Eight dynamic multiobjective benchmark functions are adopted to test whether or not the algorithm can find robust paretooptimal solutions set. The first five functions are FDA1FDA5 presented by Farina et al. [3]. The other three functions are DMOP1, DMOP2, and DMOP3 [22]. FDA4 and FDA5 contain three objectives, and the others include two objectives. is the generation counter. is the number of iterations under the time window . is the number of distinct steps under , which controls the distance between two consecutive PSs.
The first type of benchmark are FDA1, FDA4, and DMOP3. For Type I problem, only the pareto sets (PSs) in the decision space dynamically change over time. However, the corresponding pareto fronts (PFs) in the objective space do not change with time. At any moment, the optimal pareto fronts are respectively , , and . FDA2 and DMOP1 belong to Type III problems, in which only the PFs in the objective space change while the PSs in the decision space remain the same. The optimal pareto front of FDA2 is that changes from a convex to a nonconvex shape. DMOP1 has a convex optimal pareto front . FDA3, FDA5, and DMOP2 belonging to Type II problems have changing PSs and PFs: , , . The definitions of these dynamic multiobjective benchmark functions are summarized in Table 1. Their true PFs when are shown in Figure 4.
(a) FDA1
(b) FDA2
(c) FDA3
(d) FDA4
(e) FDA5
(f) DMOP1
(g) DMOP2
(h) DMOP3
In this paper, we adopt a multiobjective evolutionary algorithm based on decomposition (MOEA/D) [21] to track the moving optimal pareto front over time. In MOEA/D, the penaltybased boundary intersection (PBI) approach is used as the surrogate model. For the benchmark functions, the population size is 100. In all experiments, timevarying moment alters with the evolutionary generation and is associated with the parameter . The larger means that the environmental parameters change more infrequently and multiobjective optimization algorithms can spend more sufficient time tracking the new pareto front. Otherwise, the less makes multiobjective optimization algorithm triggered after detecting the new environment, hardly finding the satisfied pareto solutions closed to the true pareto front during the limited iterations. Now we discuss the algorithm performances under different frequencies. As shown in Figure 5, when the environment changes every 20 generations, the algorithm can track dynamic paretooptimal fronts better. But within 5 generations, the satisfied PFs for the new environment are difficultly obtained. As a result, for all experiments, . The corresponding optimal pareto fronts of each benchmark function gotten by MOEA/D are shown below in Figure 6.
(a) FDA1
(b) FDA1
(c) FDA2
(d) FDA2
(e) FDA3
(f) FDA3
(g) FDA4
(h) FDA4
(i) FDA5
(j) FDA5
(k) DMOP1
(l) DMOP1
(m) DMOP2
(n) DMOP2
(o) DMOP3
(p) DMOP3
(a) FDA1
(b) FDA2
(c) FDA3
(d) FDA4
(e) FDA5
(f) DMOP1
(g) DMOP2
(h) DMOP3
5.2. Simulation Results and Analysis
In this section, two groups of experiments have been done. In the first group, the optimal nondomination solutions starting from the time are obtained. If the minimum survival time of this moment is , robust pareto solution fits for the environment from to . The process is repeated until the last environment occurs. The purpose of the second group of experiments is to obtain each robust optimal pareto front, respectively, for RPOOT at 100 timevarying moments.
5.2.1. The Effect of Neighborhood Size
The neighborhood size directly influences the evaluation criterion for the robust optimal pareto fronts. In the first group of experiments, we compare and analyze the performances of the robust optimal pareto fronts under different thresholds . The numbers of robust pareto fronts (NRPFs) on 100 timevarying moments are listed in Table 2. The statistical average results of RGD and RIGD on eight benchmarks over 15 runs also can be found in Table 3 under different neighborhood sizes .
As shown in Tables 2 and 3, with the increasing of the neighborhood size , less robust pareto fronts are contained in the whole period. It means that the robustness of RPOOT is better. At the same time, their average inverted generational distances become larger, which means that the convergence of the RPOOT is worse. For each benchmark function, the standard errors of NRPF, RGD, and RIGD over 15 runs are given in Tables 2 and 3. The means and standard errors shown in these tables indicate that the stability of RPOOTs is good enough. We choose in the following experiments. The robust optimal pareto fronts obtained from MOEA/D are shown in Figure 7.
(a) 37 robust PFs of FDA1
(b) 7 robust PFs of FDA2
(c) 51 robust PFs of FDA3
(d) 43 robust PFs of FDA4
(e) 59 robust PFs of FDA5
(f) 1 robust PFs of DMOP1
(g) 40 robust PFs of DMOP2
(h) 80 robust PFs of DMOP3
From Figures 7(a)–7(h) and Table 4, we find that the GD and IGD of tracking multiobjective (TMO) algorithm are less than RGD and RIGD of RPOOT algorithm. But the number of robust pareto solution sets of RPOOT is far less than 100.
5.2.2. The Average Survival Time
The second group of experiments record each robust pareto front in RPOOT at each timevarying moment. Independent 15time run is done for each benchmark. The results of the second group of experiments are plotted in Figure 8. It can be seen from Figure 8 that the average survival time of the robust pareto front at each timevarying moment is far larger than 1. The results of the average survival time listed in Table 5 are all far longer than changing time 100.
(a) Average survival time of FDA1
(b) Average survival time of FDA2
(c) Average survival time of FDA3
(d) Average survival time of FDA4
(e) Average survival time of FDA5
(f) Average survival time of DMOP1
(g) Average survival time of DMOP2
(h) Average survival time of DMOP3
6. Conclusions
Dynamic multiobjective optimization problems with changing parameters widely exist in real life. The aim of the traditional optimization algorithms is to track the optimal pareto solution set after detecting the environment change efficiently. These algorithms may not obtain the satisfied nondominant solutions between two timevarying moments. In this paper, we proposed a new perspective for solving DMOPs with consecutive timevarying periods. Its goal is to find the robust pareto solution set over time. Three contributions are contained in RPOOT. At first, the detailed concept of the robust paretooptimal over time is pointed out by the robustness definition of nondomination solution in the time scale. Secondly, we developed the new definition survival time which means how many timevarying environments it fits for. Thirdly, a framework for finding robust pareto fronts is proposed, and a MOEA/D is employed as an optimizer. Lastly, eight dynamic multiobjective benchmark functions are adopted to indicate the feasibility of the algorithm. From the simulation results, we receive that the robustness of RPOOT depends on the parameter . Less robust pareto fronts are contained in the whole period with the increasing of the neighborhood size . It means that the robustness is better. At the same time, the convergence of the RPOOT is worse. Moreover, the survival time of the robust pareto front at each timevarying moment is far larger than 1. That is, each robust pareto front of the robust paretooptimal solution set can fit for more than 1 consecutive changed environments.
Due to the fact that the dynamic system is performed online, we do not obtain the future landscape on the current environment. The future fitness values of the solutions considered in RPOOT should be predicted through the past fitness values. So, the estimation and prediction task are inevitable in the future work.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This research is jointly supported by the Fundamental Research Funds for the Central Universities (no. 2013XK09) and by National Basic Research Program of China (2014CB046300).