p -Optimality-Based Multiobjective Root System Growth Algorithms for Multiobjective Applications

. In this paper, a multiobjective root system growth algorithm-based p -optimality ( p -MORSGA) is proposed. The proposed p -MORSGA extended original root system growth algorithm with multiobjective nondomination strategy. To enhance its eﬀect of convergence of solution groups, the p -optimality criterion is employed to determine the solutions of last nondominated front into the next generation group. In the evolution process, global general (GG), concerning the margin information and population density, is selected as the suitable optimality criterion of evaluating the performance of solutions. Application of the new p -MORSGA on several multiobjective benchmark functions shows a marked improvement in performance over the modiﬁed classical MOEAs with such criterion. Finally, the proposed p -MORSGA is applied to solve two real-world problems, multi-objective portfolio optimization problems (MOPOPs) and multiobjective optimal power ﬂow (OPF) problems. The experimental results demonstrate that p -MORSGA is extremely eﬀective for real-world application problems.


Introduction
For multiobjective optimization problems (MOPs), many objectives should be optimized simultaneously under some constraints.Owing to the problems' particularity, researchers prefer to find an approximation of the Pareto optimal set to support the effectiveness of their proposed approaches.At present, MOPs have made some achievements in academia and industry [1,2], but the complexity is still a noticeable puzzle.Since the traditional linear nonlinear programming method cannot effectively solve MOPs, MOEA has become the main method to obtain approximate Pareto frontier (PF).erefore, it has been widely used in power dispatch [3], logistics system optimization [4], motor design [5], artificial neural network optimization [6], and other applications.
In the existing literature [10], root system growth algorithm (RSGA) is presented to solve single-objective problems.e root system growth algorithm is inspired by natural plant root system.Plant root system is characterized by the ability to adapt to the changing diversity of the environment [11][12][13].e goal of this paper is to use RSGA to solve multiobjective problems such as MOEAs.As conventional evolutionary algorithms for MOPs, RSGA should design multiobjective mechanism to balance its diversity and convergence.
e constructed and updated Pareto nondominant set can be approximately represented by the real optimal frontier.is method will obtain an acceptable result when the advantage relationship in the target space is obvious, while it will show poor performance when the advantage relationship is weak or when the multimodal and deceptive problems are optimal in isolation [14].
In view of this situation, Pareto dominance processing technology is divided into the following types: (i) Randomized schemes, such as roulette selection [15], result in inefficiencies and poor solution caliber.(ii) Pareto-based relaxed domination methods, such as ε-domination [16], r-domination [17], grid-based method [18], etc. (iii) Reference-point-based dominance methods [19,20] adjusting solutions according to reference points.(iv) In the k-optimality used in literature [21], the dominant solution may not perform well on some specific targets, which is acceptable to decision makers.e number of targets was previously specified by adjusting the k value.(v) Optimality criterion in [22] is adopted to search the optimal spatial position meeting certain predetermined criteria, and this hypothesis is assumed to be the corresponding optimal result of the problem.(vi) Average ranking (AR) [23], preference order ranking (PO) [24], global general (GG) [25], and so on.
In this paper, p-optimality criteria combining RSGA are proposed to solve MOP.p-Optimality criterion is applied to selection operators in the evolutionary process.Such criteria are helpful to distinguish the viable solutions among solutions at the same nondominated level.erefore, it is possible for them to converge better on feasible solutions at a later stage of evolution.Meanwhile, global general (GG) criterion is selected as the suitable optimality criterion of evaluating the performance of solutions.Global general, a ranking-based dominant mechanism, uses the margin information of solutions with concerning population density (Euclidean distance among the solutions).
With the test against several multiobjective benchmark functions, the results demonstrate that the proposed p-optimality-based multiobjective root system growth algorithm (p-MORSGA) outperforms NSGA-II, MOEA/D, MOPSO, rNSGA-II, and RPDNSGA-II (classical multiobjective algorithms with p-optimality criteria) on convergence metric and diversity metric.In order to test the effect of real-world multiobjective applications with p-MORSGA, multiobjective portfolio optimization problems (MOPOPs) and multiobjective optimal power flow (OPF) are employed.
In [26], the proposed multiobjective portfolio selection model has been transformed into a single-objective programming model using fuzzy normalization and uniform design method.
e model is then applied to monthly data of top 50 companies of Tehran Stock Exchange Market in 2013.
e experiment concluded that IWO algorithm has better solving time than PSO algorithm.Mehlawat uses the return, risk, and liquidity model in [27] to better solve the problem of fuzzy combination parameters in financial markets.Jiang proposed a general multiobjective CVaR model in time interval for loan portfolios in [28].e results prove that this model can better deal with the actual financial risk problem.
Yuan et al. [29] proposed an improved strength Pareto evolutionary algorithm and applied it to the OPF problem.
e experiment results show that the improved algorithm has obvious performance.In [30], the paper presents an enhanced self-adaptive differential evolution with mixed crossover (ESDE-MC) algorithm to solve the multiobjective OPF problems with conflicting objectives that reflect the minimization of total production cost, emission pollution, L-index, and active power loss.Teeparthi and Vinod Kumar [31] proposed a new low level with teamwork heterogeneous hybrid particle swarm optimization and artificial physics optimization (HPSO-APO) algorithm.e algorithm balances global and local search capabilities.Hence, the proposed hybrid method can be used for the large interconnected power system to solve MO-SCOPF problem with integration of wind and thermal generators.
Portfolio management refers to the investment manager's diversification management of assets according to asset selection theory and portfolio theory to achieve the purpose of diversifying risks and improving efficiency.e application of EAs in solving MOPOP has attracted wide attention, yet their method is usually to convert MOPOP into single-objective POP by weighted sum.erefore, OPF problem is to schedule the promised generator set to meet system load requirements with minimal costs under some constraints.OPF problem is a multiobjective nonlinear constrained optimization problem of competition.e experimental results show that p-MORSGA outperforms MOEAs in terms of accuracy and computation.e remainder of this paper is organized as follows.Section 2 introduces related works of this paper, p-optimality criterion and global general criterion.Section 3 describes original root system growth algorithm.Section 4 describes a novel optimization model called p-MORSGA.In Section 5 provides the extensive experimental results and discussion.In Section 6, p-MORSGA for MOPOP and OPF is presented.Finally, Section 6 gives the conclusion.[32,33] are defined as follows:

p-Optimality Criterion. p-Optimality criteria
where S i (v) is the rank according to the objective function f i by means of the quick-sort algorithm and |K| is the size of feasible solutions.erefore, the goal is to search for a vector Consider the following functions: where p > 0, and p i (a) is calculated as equation (1).Such criterion promises the most feasible solution in the same nondominant rank.

Global General Criterion
Definition 1 (global margin, GM).GM is defined as the sum of all the individual target values of the difference: where M is set to the number of objectives and X i and X j are two different solutions.
According to Pareto dominance, the smaller GM(X i ), the more X i will dominate.
e framework combines each goal with the information of all individuals to obtain ranking values.By comparing with the solution pairs, the sum of the good or bad parts of the solution pairs is calculated, as shown in Figure 1.

Definition 2 (global density, GD).
where GM(X i ) denotes the global density of the particle X i and d i,j being the Euclidean distance.
where GG(X i ) represents the overall ranking of X i .e smaller GG(X i ) indicates that the decent GD of X i has a better distribution.

Multiobjective Root System
Growth Algorithm is section describes the p-optimality-based multiobjective root system growth algorithm (p-MORSGA).
e root system is defined as a group of tips, which are expressed as where θ t i � <x t i , f t i , n t i , α t i > denotes a single root tip; P t denotes the root tips' number; T is the last time; and θ t i has its own position x t i , fitness f t i , nutrient n t i , and the auxin α t i .

Auxin Concentration.
In initialization, P 0 are generated in D-dimensional space.e root tips i will forage for nutrients, which can be updated in the following ways: where GG is the global general metric introduced in Section 2.
In the initial stage, the nutritional value of each root tip is set as 0. During root growth, if its new position is better than the last one, the tip will receive nutrients from the environment.Inversely, the tips will lose nutrients and the nutrient content will be decreased by one.
en, the auxin concentration of auxin α t i is manipulated according to the following equation, which combined with the health status: where J worst and J best are the worst and most suitable place and GG is the global general metric.All root taps are classified according to the auxin concentration values.e stronger tip has a higher probability of becoming the main root.e size of main roots is limited as follows: where S t m is the number of the main root group, P t is the total number of root tips, and Q is the selection factor.e other S t l � P t − S t m root tips are considered as lateral roots.

Root Branching.
e threshold B G is to determine whether main root can be regarded as branch.If α t i is greater than B G , the root is considered as the branch root.e number of branching w i of θ t i is decided by the following formula: Mathematical Problems in Engineering where S max and S min are the maximum and minimum limits.e new branching tips foraging in the new regions will grow.e fitness f is the result of applying the following functions to each successive group of three components: where l max is the maximum of root elongation length.

Tropisms.
Trajectory of the root is influenced by different tropisms.In the RSGA model, two typical orientations are realized, i.e., hydrotropism and gravitropism.
Half of the main roots will grow to the optimum position, and the water content in the roots will be the largest, which is given by the following formula: where i ⊆ [1, (S t m /2)] and x t best is the best position in the root tip group.

Root Tip Death.
Lower auxin concentration indicates that the root tip does not acquire as much nutrients as possible during its foraging and therefore is not active and continuous growing is not possible.e algorithm can be divided into six parts.e first part is to initialize the population and calculate the fitness and GG values.e second part is to calculate the auxin concentration values according to the content of the chapter II and then sort them to divide the main root group and the lateral root group.e third and fourth parts are the growth of the main root and the lateral root.Elimination of dead roots is the fifth part.We remove the dead root tips (auxin concentration 0) from the root tip group.e sixth part is to determine the solution of the next iteration: first combine parent group and offspring group Q t � P t ∪ Q t , secondly, perform nondominated sorting, and get {F 1 , F 2 , . ..}, in the critical layer F i , where en, the p-optimality will be used to sort the individuals in the critical layer, where the smaller the p-function is, the better the individual is.Finally, select the smallest to make the number of P t+1 equal to the number of P 0 .

Tests
4.1.Experimental Setting.In the two-objective test experiment, population size is set as 200, and the number of function evaluations (FEs) is set as 40000.In the threeobjective test experiment, population size is set as 300, and (1) 1: Initialization.
(2) Iteration � 0; (3) Initialize P 0 , the positions of root tip group; (4) Evaluate the fitness of the group; (5) Evaluate the global general (GG) of each solution by equations ( 4)-( 6) (6) 2: Do while (unsatisfactory conditions) (7) 3: Auxin concentration calculation (8) Compute the auxin concentration values of the root tip group by equation ( 9); (9) Sort the root tip group in order of descending accumulated auxin concentration; (10) Divide the root tip group into two subgroups, i.e. the main root group and the lateral root group of equation (10).( 11) 4: Root branching (12) Determine the branch number of each branching root tip by equation ( 11); (13) 5: Main roots growing (14) Sort the main root subgroup in order of increasing cumulative auxin concentration; (15) e first half main roots are hydrotropic according to equation ( 12); ( 16) 6: Lateral roots growing (17) All the lateral root tips perform random walk operation according to equation ( 13 1, four metrics are used to compare the performance of the algorithms.In order to facilitate reading, the best test results are marked in bold.In general, p-MORSGA performs best, followed by the MOEA/D.Figures 2-6 show that MOEAs produce poor results on these test functions and they are almost impossible to achieve true PF, while p-MORSGA has great potential to approach true PF.Obviously, for the uniformity of the final solution, p-MORSGA has better results for the ZDT series test function.At the same time, MOEA/D has the secondbest effect, and rNSGA and RPDNSGA-II perform poorly.In general, p-MORSGA is superior in terms of uniformity and coverage, followed by MOEA/D; NAGA-II is not uniform, rNSGA-II and RPDNSGA-II are not ideal, and MOPSO algorithm has the worst performance except ZDT3.e true Pareto frontier is shown in Figure 7 for three-objective functions.e results    2 shows the test result on three-objective experience.

ree-Objective Functions.
It can be seen from Figure 8 that p-MORSGA and NSGA-II have relatively complete edge shapes, the MOEA/D algorithm's distribution is relatively uniform, followed by rNSGA-II , and MOPSO and RPDNSGA-II do not get a good frontier.
From Figure 9, p-MORSGA and MOEA/D have obtained a good frontier, and the MOEA/D distribution is uniform.As can be seen from Table 2, p-MORSGA has better convergence and diversity.As can be seen from Figure 10, p-MORSGA can get a better solution, followed by MOEA/D.In addition, MOEA's performance is slightly worse than that of p-MORSGA.However, p-MORSGA is superior to MOEAs for the problem.

Time Complexity Analysis.
In order to demonstrate the difference in the time complexity of these algorithms, Figures 11 and 12 plot the average CPU time over 10 runs. Figure 11 shows the CPU time for six algorithms and gives the results on five benchmark functions (ZDT1, ZDT2, ZDT3, ZDT4, and ZDT6).Overall, rNSGA-II performs poorly, and the running time is 3.02-3.56times that of other algorithms.
e time complexity of p-MORSGA is better than rNSGA-II and RPDNSGA-II and worse than NSGA-II.NSGA-II has the best performance with time complexity, and p-MORSGA is about 1.04-1.49times that of NSGA-II.It is worth mentioning that although p-MORSGA consumes more running time than NSGA-II, the performance has been significantly improved.
Figure 12 shows the CPU time for six algorithms and gives the results on three benchmark functions (DTLZ1, DTLZ2, and DTLZ3).
ere is some difference as the number of objectives increases.
e time complexity of rNSGA-II and other algorithms is reduced.rNSGA-II is about 1.01 to 2.67 times that of other algorithms.In each test function, the best performance of time complexity is still NSGA-II.Secondly, the best performance of time complexity is p-MORSGA, and p-MORSGA is about 1.02-1.26times that of NSGA-II.
Compared to Figure 11, the running time of the algorithm increases overall.It can be seen that the threeobjective benchmark functions are more difficult to calculate than the two-objective benchmark functions.[38] proposes an MV portfolio model, which is a formal twoobjective portfolio model.e model uses the covariance and average return of assets to describe the risk and return of investment, respectively.Two conflicting aspects are considered by this model: maximizing the portfolio's expected return while minimizing its risk [39].

Two-Objective Portfolio Model. Markowitz
e MV model is constructed as follows [38]: where σ 2 are portfolios, w i and w j are the weights, σ ij expresses covariance between two assets, N denotes quantity, and r i indicates expected return of the ith asset.e assumptions of several portfolio models introduced are too harsh to meet the actual needs of the securities market, which leads to the deviation between the results of the model operation and the actual situation.On the basis of revising the hypothesis of the classical portfolio model, semivariance is introduced to replace variance, which makes the model more reasonable, more in line with the investment and financing environment of China's financial market, and provides more effective auxiliary tools for investors.e MV model has the function as follows:

ree-Objective Portfolio
e weight of each asset portfolio should be the sum of 1.

Mathematical Problems in Engineering
In [39], semivariance is described as follows: where B represents the comparative return and R it represents the earnings of asset I at period T. e indirect method is utilized to represent transaction cost: erefore, the return-risk-cost portfolio model is constructed as follows: , minimizing expected cost. ( Budget constraints are given as follows: It should be noted that when appropriate indicators are available, the number of model objectives can also increase.

Applications for Portfolio Problem.
e experiment uses daily historical data of 12 kinds of assets from Shanghai Stock Exchange, which are collected at the monthly rate of each stock from January 2010 to December 2016.e frontier of p-MORSGA is searched and shown in Figures 13 and 14.In addition, since the true PF of the portfolio problem is not yet clear about p-MSMOEAs, all acquired PS are considered to be true PF [40].e comparison results are listed in Tables 3 and 4, where hypervolume is used to test their performance.
Figure 13 shows the nondominant solution of the MV model generated by p-MORSGA.Investors can choose portfolio selection methods according to their risk preferences in order to obtain corresponding returns.In Figure 13, the effective curves of NSGA-II and the proposed algorithm obviously are continuous, but MODE locates the opposite.As can be seen from Table 3, the performance of p-MORSGA, NSGA-II, and MOEA/D in hypervolume index decreases in turn.
From Figure 14, when considering, risk and expected return are nearly contrary.e case of costs and risks, costs and expected benefits, regardless of another variable, is nearly the same.For facing all three variables, investors can choose portfolio selection methods according to their preferences.
As is visible from Figure 14, nondominant solutions' distribution in Figures 14(a) and 14(c) is more uniform and diversified.Obviously, as can be seen from Table 4, p-MORSGA outperforms NSGA-II and MODE in terms of hypervolume indicators.

Optimal Power Flow Problem Formulation. OPF's main
objective is to optimize restriction variables' setting while fulfilling several inequality constraints and equality.Generally speaking, OPF issues can be expressed mathematically as follows: where F represents the objective function, g(x, u) is the equality constraint, and h(x, u) is the system operation constraint.Here, u is a vector of independent control variables, including where P G1 is the active power output of the generator; V L is the load (PQ) bus voltage; Q G is the reactive power output of the generator; S L is the transmission line load; N pq denotes the number of PQ buses; and N l denotes the total number of transmission lines.

Objective Function.
In this paper, there are three competing objective functions regarding the OPF problem, namely, total fuel cost, total power loss, and total emission cost, while satisfying several inequality constraints and equality.e problem is usually expressed as follows.
(1) Minimize Total Fuel Costs.e curve is expressed as follows: where a i , b i , and c i are equivalent to the fuel cost coefficient of the ith generator and P Gi is the actual power output of the ith generator.(2) Minimization of Total Power Losses.Minimize total power loss as follows: where N l represents the number of transmission lines, V expresses the voltage magnitudes, and δ represents the voltage angles.(3) Minimization of Total Emission Cost.e total cost of emission is defined as [19] f emission � where f emission represents total emission cost.

Equality Constraints.
From (20), the equality constraint g(x, U) is defined as follows: where P and Q are the active and reactive power loads, respectively, and G and B are the real and imaginary parts.

Inequality Constraints.
e inequality constraints are the power system handling.
(1) Generator Constraints.Generator active power PG, generator reactive power QG, and generator voltage amplitude QG: (2) Transformer Constraints.Transformer taps have restrictions as follows: (3) Switchable VAR Sources.e switchable VAR sources have minimum and maximum setting limits: (4) Security Constraints.e limitations of load bus voltage amplitude and transmission line flow limitation are as follows:

Multiobjective Optimal Power Flow Based on p-MORSGA.
To verify the numerical correctness and the effectiveness of p-MORSGA, the multiobjective OPF problem will be used to implement simulation experiments.is algorithm will be applied to two objectives and three objectives of the OPF problem, and this problem is solved.

Steps for p-MORSGA for OPF Problems
Step 1. Input the system parameters and the minimum and maximum limits of control variables.
Step 2. Input parameters of p-MORSGA and the lower and upper limits of each variable.Mathematical Problems in Engineering Step 3. Input the minimum and maximum limits of variables.
Step 4. Produce the initial population.
Step 5. Cluster the population by K-means and compute the fitness.
Step 6. Store nondominated solutions in the memory table.
Step 8. Handle the constraints.
Step 9. Compute the fitness of the solutions and store nondominated solutions into EA.
Step 10.Update each EA of each hive.
Step 11.Judge the condition of reclustering.If it meets the condition, the number of clusters can be reset and go to Step 4; if not, go to Step 7.

Simulation Results.
e standard IEEE 30 bus system has been applied as a test system [23].e system contains six generators, 41 transmission lines, and 4 transformers.
e system data are shown in Table 5 [41].
In this section, we consider the following two test system scenarios: (1) Two objectives: the emission cost, the emission loss, and the loss cost are considered separately.(2) ree objectives: emissions, costs, and losses.
e new algorithm is compared with three MOEAs in three aspects: minimum cost, minimum emission, and minimum loss.Figures 15-17 show the performance of four algorithms, respectively.show the two-objection compromise solution data.Finally, Table 7 shows the best solutions for minimum cost and emission.Table 8 provides data on emission losses and loss costs.Table 9 shows the best solutions for loss and emission.
First, two competing objectives need to be considered: fuel costs and emissions.From Figure 15, we can make out that the p-MORSGA method can provide a uniformly distributed solution.From Table 7, it is clear that the minimum fuel cost that is obtained by p-MORSGA is 605.62 $/h, which is more wonderful than the results that are e objective function of cost loss is shown in Figure 16.We can see from Figure 16 that the proposed algorithm obtains better results than the other three algorithms.
From Table 9, it can be concluded that the p-MORSGA method can provide lower economic fuel cost when considering the minimum loss, while the p-MORSGA method can provide smaller loss when considering the optimal fuel cost effect.
For the emission loss objective function, we can come to know that the algorithm achieves performance ranking similar to that of the fuel cost emission objective function.Mathematical Problems in Engineering (2) Case II: Tree-Objective OPF Optimization.In this instance, each algorithm needs to optimize three objectives.Figure 18 is the convergence curve and Table 13 is the list of data results.Compared with the four convergence effects of Figure 18, the Pareto optimal front solutions obtained by p-MORSGA are uniformly distributed and widely distributed, while the results obtained by MOABC in the other three algorithms are better.Judging from the data in Table 13, for p-MORSGA, when F1 (cost) is 610.99 $/h, F2 (emission) is 0.2503 ton/h, and F3 (loss) is 2.0021 MW, F1 is greater than 610.99 $/h for the other three algorithms and the other two objectives are relatively poor when F1 is the smallest.It can be concluded that p-MORSGA can obtain better solutions than other algorithms on this three-objective problem.
Referring to the data in Table 6, the proposed p-MORSGA obtains the best pollution emission value and generation cost.Comparatively speaking, p-MORSGA has the best optimization effect.

Conclusion
In this paper, a multiobjective optimizer based on p-optimality and global general (GG) is proposed, which is called p-MORSGA.
rough eight mathematical benchmark functions, p-MORSGA has an outstanding performance on accuracy and convergence.It can be inferred that p-MORSGA is simple in structure, easy to implement, and has great potential to solve complex MOP.In order to prove its effectiveness in solving multiobjective realistic problems, the performance of p-MORSGA is tested by using multiobjective combinatorial optimization problem (MOPOP) and multiobjective optimal power flow (OPF) problems.e MOPOP includes three objectives: expected return, risk, and expected cost.
e experimental results show that the proposed p-MORSGA is considered to be able to obtain the distribution of multiple solutions.For OPF problem, a 30 bus IEEE test system is used to test the proposed algorithm.Compared with three other OPF problems with different objectives, the proposed method can provide uniformly distributed Pareto optimal solutions.

Figure 1 :
Figure 1: A sketch map on global margin ranking.

3. 5 .
Pseudocode of p-MORSGA.p-MORSGA algorithm code and flow chart are listed in Algorithm 1.
Model.Different from the 2objective model, the 3-objective portfolio model includes other objective: minimizing expected cost.

Figure 11 :Figure 12 :
Figure 11: CPU time of MOEAs on two-objective benchmark functions.

Figure 13 :Figure 14 :
Figure 13: PF obtained by p-MORSGA and p-MOEAs on MV portfolio model: (a) PF obtained by p-MORSGA on MV portfolio model; (b) PF obtained by NSGA-II on MV portfolio model; (c) PF obtained by MOEA/D on MV portfolio model.

: Determine the solutions of the next iteration
Remove the dead root tips (auxin concentration 0) from the root tip group; (20) 8: Generate offspring group (21) Gather the new generated solution as offspring group O t ; (22) 9(23) Combine parent group and offspring group Q t � P t ∪ Q t (24) Non-dominated sort on Pareto dominance {F 1 , F 2 , . ..} (25) Add nondominated front F i into the next group P t+1 one by one until |P t+1 ∪ F i | ≥ |P 0 | (26)Add the solutions of last nondominated front F i with the largest p-optimality criterion one by one until|P t+1 | � |P 0 |.(27) 8: Memorize the best solution to date (28) Iteration � iteration + 1; (29) End while (30) 9: Output the best solution achieved.ALGORITHM 1: Pseudocode of p-MORSGA.4 Mathematical Problems in Engineering the number of function evaluations (FEs) is set as 90000.e number of independent runs of the experiment is 10 times.For p-MORSGA, NSGA-II, MOEA/D, MOPSO, rNSGA-II, and RPDNSGA-II, parameter settings are the same as the ones in [8, 17, 34-37], respectively.For p-MORSGA, the whole swarm has 200 individuals, and parameter settings are the same as the ones in [32].4.2.Results and Discussion 4.2.1.Two-Objective Functions.Table 1 and Figures 2-6 demonstrate the results of p-MORSGA and MOEAs against ZDT series.In Table

Table 1 :
Test results of two objective functions.

Table 2 :
Test results on three-objective functions.
Min σ 2 �  N i�1  N j�1 w i w j σ ij , minimizing risk, Max r �  N i�1 w i r i , maximizing expected return, subject to  N i�1

Table 3 :
Comparison of performance on MV portfolio model.

Table 4 :
Comparison of performance on return-risk-cost portfolio model.

Table 5 :
Characteristics of the generation units.

Table 6 :
e best compromise solutions for cost, emission, and loss.

Table 7 :
e best solutions for cost and emission.

Table 8 :
e best solutions for cost and loss.

Table 9 :
e best solutions for emission and loss.

Table 10 :
e best compromise solutions for emission and cost.

Table 11 :
e best compromise solutions for loss and cost.

Table 12 :
e best compromise solutions for loss and emission.otherthreealgorithms.And when the best fuel cost effect case is considered, the p-MORSGA method can provide smaller loss.e best compromise solutions are 614.8302$/hand0.1998 t/h (p-MORSGA), as shown in Table10.Once again, testify that the method can get a better solution.

Table 13 :
e best solutions for cost, emission, and loss.