Biogeography-Based Optimization of the Portfolio Optimization Problem with Second Order Stochastic Dominance Constraints

: The portfolio optimization problem is the central problem of modern economics and decision theory; there is the Mean-Variance Model and Stochastic Dominance Model for solving this problem. In this paper, based on the second order stochastic dominance constraints, we propose the improved biogeography-based optimization algorithm to optimize the portfolio, which we called ε BBO. In order to test the computing power of ε BBO, we carry out two numerical experiments in several kinds of constraints. In experiment 1, comparing the Stochastic Approximation (SA) method with the Level Function (LF) algorithm and Genetic Algorithm (GA), we get a similar optimal solution by ε BBO in [ 0, 0.6 ] and [ 0, 1 ] constraints with the return of 1.174% and 1.178%. In [ − 1, 2 ] constraint, we get the optimal return of 1.3043% by ε BBO, while the return of SA and LF is 1.23% and 1.26%. In experiment 2, we get the optimal return of 0.1325% and 0.3197% by ε BBO in [ 0, 0.1 ] and [ − 0.05, 0.15 ] constraints. As a comparison, the return of FTSE100 Index portfolio is 0.0937%. The results prove that ε BBO algorithm has great potential in the ﬁeld of ﬁnancial decision-making, it also shows that ε BBO algorithm has a better performance in optimization problem.


Introduction
In solving the problem of uncertainty, the Expected Utility Theory describes the rational people how to determine the optimal decisions and then take action when faced with the uncertainty of risk and return.The Expect Utility Theory assumes that the individual is risk aversion, that is, its utility function is concave.In modern economics, how to choose the optimal portfolio is one of its central problems.Based on the relationship between the risk of assets and the return on assets, Markowitz [1] solve the portfolio optimization problem through mathematic statistic method and propose the Mean-Variance model.The Mean-Variance model measures the expected revenue at the required rate of return, and measures the risk size with the variance of required rate of return.However, the model gets poor applications because of its harsh assumptions.
To supply a gap, Fishburn [2] propose the Stochastic Dominance theory and apply it into portfolio optimization problem.The Stochastic Dominance theory takes the risk appetite into consideration, which has been widely used in the financial field.Thereinto, second order stochastic dominance (SSD) is a branch of Stochastic Dominance theory.
Darinka Dentcheva and Andrzej Ruszczy ński [3,4] introduce the second order stochastic dominance into portfolio optimization model.They also take the portfolio selection theory as a constraint and apply it to the risk repugnance asset optimization problem.While the second order stochastic dominance constraint requires that the investable assets should be compared with each other.When the number of investable assets is large, the calculation requirements becomes very harsh.Rudolf and Ruszczy [5] propose a new duality theory for this problem and optimize the model by cutting plane methods.Homem-De-Mello and Mehrotra [6] and Fábián, Mitra and Roman [7] also use the cutting plane methods to solve the portfolio optimization problem.The cutting plane method has a great efficiency in portfolio optimization problem, while this method is too restrictive.Besides, the probability space is discrete, which greatly limits the application of this method.Dentcheva and Ruszczy ński [8] introduce a subsidiary variable into SSD model and transform it into a linear programming model.However, the number of variables in the model depends on the number of investable assets, which increases the pressure of calculation.Meskarian, Xu and Fliege [9] penalize SSD constraints to the objective under Slater's constraint qualification and apply the SA method and LF method to solve the penalized problem.Hu [10] apply the Sample Average Approximation method into SSD model, which introduces new variables into the model.Besides, Roman, Mitra and Zverovich [11] present an empirical study which analyses the effectiveness in the context of enhanced indexation.
Inspired by Biogeography [12], Dan [13] propose the Biogeography-based optimization (BBO) which is a sort of bionic optimization algorithm consists of three operators: migration operator, mutation operator and elimination operator.As a swarm intelligence optimization algorithm, the Biogeography-based optimization has advantages of good ability of robustness, easy to implement and the process of algorithm is relatively simple.Besides, the mutation operator helps to increase the diversity of solutions.Because the BBO algorithm does not depend on a particular problem, it is widely used in various fields, such as pathology test [14,15], distributed generation [16], intelligent identification [17][18][19] and other fields recently.While BBO algorithm also has some disadvantages, including easy to fall into the local optimization, poor exploration ability, strong randomness and lacking theoretical basis.In this paper, to increase the exploration ability of algorithm, we propose the εBBO algorithm based on the DE algorithm.We improve the mutation operator to prevent the algorithm falling into the local optimization.Besides, we introduce it into portfolio optimization problem.
The rest of the paper is organized as follows.In Section 2, we discuss the SSD constrains model.In Section 3, we discuss the εBBO algorithm.In Section 4, we apply the proposed methods to portfolio optimization problems and report some numerical test results.In Section 5, we summarize the performance of εBBO on SSD constrains model.

The Model of Second Order Stochastic Dominance
In this section, we discuss the SSD constraints model for portfolio optimization problem.We start by discussing preliminaries needed.

The Expected Utility Theory and Mean-Variance Model
The expected utility is proposed by Bernoulli in 1738.Based on the explanation to St. Petersburg Paradox, Neumann and Morgenstern [20] discuss the strict axiomatic assumptions and propose the Expected Utility Theory.Let U denote the expected utility of a certain decision: where u is utility, p i and w i is the probability and return of various circumstances.However, the Expected Utility Theory can not explain some investors' behavior, such as in buying lottery tickets and insurance at the same time.The emergence of Allais Paradox and Ellsberg Paradox proves that the Expected Utility Theory has a big loophole in derivation of theoretical model.Laciana and Weber [21] introduce the Regret and Disappointment Theory into expected utility function, which explain the Allais Paradox and Ellsberg Paradox.Generally speaking, the main problem of Expected Utility Theory is lacking practical application.
It is generally believed that modern investment theory begins with the portfolio selection theory proposed by Markowitz.Markowitz mainly solves the two problems: why we need to select portfolio and how to choose optimal portfolio.Therefore, Markowitz propose the following portfolio optimization model: The mean-variance model is based on several assumptions, such as effective securities market, risk-aversion investors and so on, which means the mean-variance model has poor applications.Therefore, in this paper, we mainly discuss the portfolio optimization problem based on the stochastic dominance theory.

The Stochastic Dominance Theory
The stochastic dominance theory originates from discrete stochastic variables optimization theory, which is later developed into generalized stochastic variables optimization theory and widely used in economics and finance.In the stochastic dominance theory, the comparison among the random variables is carried out the point-wise comparison of their k-order distribution function F (k) .Assuming that there are only two portfolios X and Y, and the utility function of all investors is monotonically increasing.If all investors prefer portfolio X or believe that only part of the portfolio X and Y is no difference, then we can say the portfolio X stochastically dominates portfolio Y in the first order.
From the mathematical point of view, let x and y be the decision vectors and ξ be the random variable.It is said that g(x, ξ) stochastically dominates g(y, ξ) in the first order, denoted by g(x, ξ) 1 g(y, ξ), if F(g(x, ξ); η) ≤ F(g(y, ξ); η), ∀η ∈ R where g(x, ξ) is the concave continuous function both in x and ξ, F(g(x, ξ); η) and F(g(y, ξ); η) are the cumulative distribution function of g(x, ξ) and g(y, ξ).For a random variable X ∈ R, the first order distribution function of X is its right-continuous cumulative distribution function: The first order of stochastic dominance is applicable to investors with arbitrary risk appetite, while the second order stochastic dominance is only for the risk averse decision makers whose marginal utility of expected utility function decreases.Based on the definition of first order stochastic dominance, similarly, g(x, ξ) stochastically dominates g(y, ξ) in the second order, denoted by g(x, ξ) The corresponding strict dominance relation (k) is defined in the usual way: In this paper, we mainly discuss the second order stochastic dominance.Assuming that the rate of return on g(x, ξ) and g(y, ξ) in portfolio x and y having a finite expectation.In the portfolio optimization problem, it is said that portfolio x stochastically dominates portfolio y in the second order, denoted by g(x, ξ) 2 g(y, ξ), its requirement is the same as model (5).

The Second Order Stochastic Dominance Model
Dentcheva and Ruszczy ński [4] propose the following fundamental model of second order stochastic dominance: Here f : X → R is a concave continuous function.In particular, we may use where g : R n × R k → R are concave continuous functions both in x and ξ, x ∈ X is a decision vector with X being a nonempty convex subset of R n , and ξ is a random vector defined on probability space (Ω, F, P) with support Ξ, y ∈ X is a predefined vector, E(•) denotes the expected value w.r.t. the probability distribution of ξ [9].The constraints of model (7) are based on the definition of second order stochastic dominance, Fábián [7] and Ogryczak [22] prove that the constraint g(x, ξ) 2 g(y, ξ) is equivalent to the following two inequalities: E(U(R x )) ≥ E(U(R y )), for any increasing and concave utility function U (10) Combined with the constraint (9), the model ( 7) can be formulated as a stochastic semi-infinite programming problem [3,9] : For the random variable η ∈ R, it is clear that there is no x ∈ X seriously satisfies the constraint condition when η takes a sufficiently small value.So the inequality does not hold and the model (11) does not satisfy the Slater condition.In the actual research, in order to overcome the serious technical difficulties and seek the optimal solution of algorithm, we get the following relaxed model: where [a, b] is a closed interval in R. Darinka Dentcheva et al. [3] propose that model ( 12) is equivalent to model (11) for some appropriate interval [a, b] if ξ has uniformly bounded distribution.As a relation of model (11), model (12) has a larger set of feasible solutions and subsequently its optimal value gives a lower bound for model.Besides, the relaxed model is more likely to satisfy the Slater condition which is closely related to numerical stability [9].
Although the constraints are relaxed, the model ( 12) is a stochastic semi-infinite and non-smooth programming problem in nature.So the difficulty of solving is still large.Dentcheva and Ruszczynski [4] change the model into a liner programming problem which solves the problem of non-smooth of SSD constraints.Besides, this method also applies to non-linear cases.However, facing the multi-dimensional problems, due to the increase in variables, the efficiency of the method will be greatly affected.
In this paper, we mainly discuss the discrete condition, assuming the random variable ξ satisfies the discrete distribution, that is P(ξ = ξ i ) = p i , i = 1, . . ., m. Fábián [7] proposes the following model: where η j = g(y, ξ j ), j = 1, . . ., m .In particular, p i = 1 m .In this case, model ( 13) is the sample average approximation of model (12).However, the model ( 13) is still non-smooth and does not satisfy the Slater condition.

Portfolio Optimization Model
In the portfolio optimization problem, our aim is to invest our capital in some assets in order to obtain some desirable characteristics of the total return on investment.
Let n denote the number of assets available for investment at the beginning of a fixes time period and assume that we have a fixed capital to be invested in n assets.To simplify the discussion, we use x = (x 1 , . . ., x n ) T to denote the fractions of initial capital invested in different assets.i = w i w (i = 1, . . ., n) where w i is the capital invested in asset i and w is the total amount of capital to be invested.Let X denote a set of feasible portfolios and it is clear that X ∈ R n is a bound convex polyhedron.Let R i (ξ), i = 1, . . ., n denote the return of asset i in the case of discrete distribution, we assume that E[| R j |] < ∞ for all j = 1, . . ., n.Therefore, if we have a fixed capital, the return of portfolio can then be formulated as: In economics and related disciplines, the transaction cost is a cost in making any economic trades when participating in the market.If we take the transaction cost into account, let c(x) denote the V-type function of transaction cost and the return of portfolio can be formulated as: where x is the investment proportion vector and ξ is a vector represents the return of assets.Clearly, the set of possible asset allocations can be defined as follows: In order to ensure the diversity of portfolio, the bounds are always set to 0.6 and 0: The two constraints proposed above based on the short-selling is prohibited.Short-selling is the sale of a security that is not owned by the seller, or the seller has borrowed.If the short-selling is allowed and investors do not want to invest very small amounts in an asset, the upper and lower bounds on the fraction of capital invested in each assets are set to 2 and −1 [9]: Assumed that a reference random return Y having a finite expected value which may be an index such as FTSE 100 index, or a current portfolio.Let y denote a benchmark investment, normally y i = 1 m , for i = 1, . . ., n.Based on the model ( 13) and the discussion above, we propose the following portfolio optimization model with SSD constraints: Firstly we consider the simplest case: Then we consider a sightly more complicated condition.In view of the impact of transaction cost and the small amounts investment in portfolio, we consider the following performance function [9]:

The εBBO Algorithm for SSD Model
Aiming at the portfolio optimization problem based on SSD constraints, drawn on the experience of mutation operator of the Differential Evolution (DE) algorithm, we propose the εBBO algorithm.In this section, we discuss the fitness function, migration operator and mutation operator for SSD model.Besides, we introduce the main procedure of εBBO algorithm for SSD model.

The Fitness Function for SSD Model
BBO is a population-based algorithm in which a population of candidate individuals is used for solving the global optimization problem.In the BBO algorithm, each habitat is considered to be an individual and has its habitat suitability index (HSI) to show the degree of its goodness.Habitat with a high level of HSI represents a good solution to the problem, while habitat with a low level of HSI represents a poor solution to the problem.The HSI is determined by several vectors, such as rainfall and the diversity of vegetation, which is named as Suitability Index Variable (SIV).
In this paper, we mainly discuss the portfolio optimization problem based on the SSD constraints.Therefore, the SIV of the portfolio optimization problem is the return of assets and transaction cost.If we do not consider the influence of transaction cost, the HSI can be express as: In the actual investment market, in order to avoid petty investments, we take the transaction cost into account and the HSI can be express as: where λ is a constant which is determined by physical truth.

Migration Operator for SSD Model
Modified migration operator is a generalization of standard BBO operator.The idea inspired by blended migration operator [23], the coefficient of solution H i is a constant.In our algorithm, along the number of iterations increases, solution H i is much fitter than solution H j , and solution H i is more affected by itself.The migration operator based on the number of iterations is designed to accelerate the speed of convergence to global optimal solution.In εBBO, the Modified migration is defined as: )H e (j) (24) where H i is a immigrated island, H e is a emigrated island, H i (j) is the jth dimension of the ith solution, and t is the number of iterations, t max is the maximum number of iterations.Equation (24) means that the features of solution are changed by solution H i and solution H e .At first, the new solution is more affected by another better solution which accelerates the convergence characteristics of the algorithm.
Along with the growth of generation, the new solution is more affected by itself.
After the migration operation, the portfolio in the immigrated island (a bad solution) accepts the sharing information from the emigrated island (a better solution).However, another portfolio in immigrated island may have the same value as immigrated portfolio.To keep the new solution feasible, therefore, we need to adjust the portfolio which has identical vector as the sharing information.That is to say, if the updated portfolio has the same value as another portfolio in the immigrated island, the portfolio swaps places with the immigrated portfolio.

Mutation Operator for SSD Model
In BBO, if a solution is selected for mutation, it will be replaced by another randomly generated new solution set.The hybridization between the BBO and Differential Evolution Algorithm (DE) has achieved many great results [24][25][26].However, they mostly incorporate DE into the migration procedure.Our algorithm incorporates the DE into the mutation procedure.The algorithm adjusts near the optimal solutions, so that it can find the global optimal solution.In εBBO, A mutated individual (H i (j)) is generated according to the following equation: where '| • |' means absolute value, H i (j) is selected for mutation, c 1 is the mutation scaling factor and its value is usually set as 0.5.H r1 (j) and H r2 (j) are the two solutions randomly selected, H best (j) is the best solution in this generation.In εBBO, this mutation scheme tends to increase the diversity among the population.It acts as a fine tuning unit and helps to achieve global optimal solution.

Adaptive BBO for SSD Model
There are modification probability and mutation probability factors in the BBO algorithm.Modification probability is denoted as P modi f and mutation probability is denoted as P muta .The two factors which range from 0 to 1 are defined by users.When generated random number is less than modification probability, the program executes migration operator; when generated random number is less than mutation probability, the program executes mutation operator.In original BBO algorithm, the modification probability and mutation probability are user-defined parameters.The setting of the parameters are very relevant to the experience of the user, and they are unfavorable for the selection of migration individuals.In order to choose better migration individuals, these parameters are related to the fitness function.When the value of fitness function is better than the average, we will reduce the probability of modification and mutation operations, or else we will increase the probability of modification and mutation operations.In other terms, in migration operations, when the value of fitness function is larger than the average, we set constant factor k 1 to smaller value, and the probability of modification operations is decided by Formula (26).When the value of fitness function is less than the average, we set constant factor k 2 to larger value, so as to search extensively in the solution space.Similarly, the selection of mutation probability is similar to the modification probability.When the value of the fitness function is better than the average, we set constant factor k 3 to smaller value.Otherwise, we set constant factor k 4 to larger value.As the mutation operation is adjusted near the optimal solution, we set the maximum mutation probability k 4 to 0.25.So modification probability and mutation probability are adjusted adaptively by the fitness function.
In the adaptive Biogeography-based optimization algorithm, modification probability and mutation probability are modified according to the relation between the cost of fitness function of randomly selected habitat and average cost of fitness function of all habitats in last generation.In other words, if the cost of fitness function is equal or greater than average cost, modification probability and mutation probability are modified by Equations ( 26) and ( 27) respectively.Otherwise, modification probability and mutation probability are taken as constant k 2 and k 4 respectively.In Equations ( 26) and ( 27), constant factors k 1 , k 2 , k 3 and k 4 which range from 0 to 1 are defined by users.In another literature [27], we used the values of adaptive factors and achieved good results.Therefore, in this paper, we adapt the same adaptive factors.We set k 1 = 0.4, k 2 = 0.95, k 3 = 0.1, k 4 = 0.25.

Main Procedure of εBBO for SSD Model
By incorporating the above-mentioned migration operator and mutation operator into Biogeography-based optimization, the inverse operation is performed after the mutation operator; in addition, modification probability and mutation probability are modified in term of Equations ( 26) and ( 27) respectively.Meanwhile, immigration rate and emigration rate based on cosine curve are modified.The εBBO algorithm is able to explore the new search space with the mutation operator of DE algorithm and to exploit the population information with the migration operator of BBO algorithm.This feature overcomes the lack of exploration of the original BBO algorithm.
The procedure chat of the εBBO algorithm is shown in Figure 1.The εBBO approach is described in Algorithm 1 and SortPopulation (Population P) is described in Algorithm 2.
Algorithm 1 the main procedure of εBBO for the SSD.P ← replaceWorst(P , Elite)

Numerical Experiments
We have carried out a number of numerical tests by using εBBO algorithms in MATLAB2014b install of a Lenovo PC with Windows 8.1 operating system and 8 GB of RAM.In this section, we report the test results.

Example 1
At first, we consider the portfolio optimization problem in low dimensional.We consider a history of percentage returns, for m = 10 time periods and a group of n = 5 assets in Table 1 [9].Our aim is to find an optimal investment strategy for a fixed capital in n assets which maximizes the expected return.If we do not consider the impacts of transaction cost, based on Table 1 and model (19), we propose the following SSD constraints model: where X has three conditions which we have discussed above.We apply the εBBO algorithm to the model ( 28), and we let y = (0.The Figures 2-4 show the performance of the εBBO and the optimal portfolio in three kinds of X constraints of model (28).
In model (28) we do not consider the transaction cost, if we take it into account, in order to weaken the contribution of small transactions we propose the following model: The same as the model ( 28), we get the test results of model (29) shown in  Return % 0.9 optimal portfolio reference portfolio SA/LF portfolio (b) The return of optimal portfolio by εBBO in each period with benchmark portfolio and the optimal portfolio by SA and LF algorithm Optimal return % 1.17 Return % 0.9 optimal portfolio reference portfolio GA100 portfolio GA1000 portfolio (b) The return of optimal portfolio by εBBO in each period with benchmark portfolio and the optimal portfolio by GA algorithm Figure 3.The performance of the εBBO and the optimal portfolio in X ∈ [0, 1] of model (28).Return % 0.8 0.9 optimal portfolio reference portfolio (b) The return of optimal portfolio by εBBO in each period with benchmark portfolio  The Figures 5-7 show the performance of the εBBO and the optimal portfolio in the three kinds of X constraints of model (29).Return % 0.9 optimal portfolio reference portfolio SA/LF portfolio (b) The return of optimal portfolio by εBBO in each period with benchmark portfolio and the optimal portfolio by SA and LF algorithm

Example 2
In the practical application of portfolio optimization, often faced with a large number of assets for investors to choose, which requires the algorithm maintain a good performance to large-scale operations.In order to test the computational power of εBBO algorithm on large-scale SSD constraints, we collect 249 daily historical returns of 101 FTSE100 assets prior to December 2016 to construct the portfolio strategy.Therefore, we get a historical of percentage returns for m = 249 time periods and a group of n = 101 assets.Besides, we use y = ( 1 101 , . . ., 1 101 ) as the benchmark and we discuss the following SSD constraints model: (g(y, ξ j ) − g(y, ξ i )) + , j = 1, . . ., 249 x ∈ X Due to the fact that the stock certificate can not always guarantee positive returns, we discuss the suitability function as 2 − g(x, ξ) to ensure its value is positive.In order to ensure the diversity of the portfolio, we discuss the X constraint as [0, 0.1] and [−0.05, 0.15] which stand for two conditions: short-selling is allowed and short-selling is prohibited.Comparing with the FTSE100 Index we get the test results shown in Table 4.The specific index weight data of FTSE100 Index is shown in Appendix A. The Figures 8 and 9 show the specific asset structure of optimal portfolio in two kinds of X constraints of model (30).The Figures 10-13 show the performance of the εBBO and the optimal portfolio in the two kinds of X constraints of model (30).

Numerical Analysis
In the numerical experiment, we try the number of initial populations in different orders of magnitude.Taking the efficiency of algorithm into account, we finally set the number of initial population to 100 in the experiment.Besides, we hold the best 10 populations in each iteration while the others carry on migrating and mutating.The model (28) and model (29) discusses the level of optimization of εBBO algorithm under the SSD in low dimension.Form the result of model (28), under the constraint of [0, 0.6], the return of optimal portfolio by εBBO is 1.1736%, while the returns of portfolio by SA and LF are 1.147% and 1.148%.Besides, the εBBO is faster in running, but the number of iteration is much more than LF.Under the constraint of [0, 1], the return of optimal portfolio by εBBO 1.1779%, while the returns of portfolio by GA are 1.165% and 1.177%.The GA(100) has the same initial population number as εBBO, while the efficiency of εBBO is better than the GA.Under the constraint of [−1, 2], because the short-selling is allowed, the return of optimal portfolio gets greatly improved.Model ( 29) is an improvement to model (28) essentially, whose enhanced suitability function can weak the small transactions.Under the constraint of [0, 0.6], three kinds of algorithm get the similar portfolio, while the εBBO is faster in running than SA and LF.Under the constraint of [0, 1], the return of optimal portfolio by εBBO is 1.1779%, which is better than the model (28).Under the constraint of [−1, 2], the return of optimal portfolio by εBBO is 1.3043%, while the returns of portfolio by SA and LF are 1.23% and 1.26%.The εBBO is also fast in running.In the comparison of model (28) and model (29), model (29) can approximate the optimal portfolio in less number of iteration than model (28).Because the model ( 29) is more complex, the running speed of model ( 29) is slower.However, the portfolio performance of model ( 29) is much better.Compared with the SA, the εBBO has great advantages in all respects.Compared with the LF, the εBBO requires large number of iteration to approximate the optimal portfolio due to its randomness.Therefore, the LF has a advantage of iterative procedure.Besides, the performance of GA is similar to LF.If having the same initial population, the efficiency of εBBO is higher than GA. the εBBO also has a strong ability to approximate the optimal portfolio.
The model (30) discusses the effectiveness of εBBO algorithm in high dimension.Under the constraint of [0, 0.1], the return of optimal portfolio by εBBO is 0.1325%, the number of asset whose investment proportion more than 2% is 12, the number of asset whose investment proportion more than 1% is 34, the number of asset whose investment proportion less than 0.5% is 39, and the maximum investment proportion of optimal portfolio is 5.17%.Under the constraint of [−0.05, 0.15], the return of optimal portfolio by εBBO is 0.3197%, the number of asset whose investment proportion more than 4% is 12, the number of asset whose investment proportion more than 2% is 35, the number of short-selling asset is 29, the maximum investment proportion of optimal portfolio is 7.97% and the minimum investment proportion of optimal portfolio is −4.79%.As a comparison, the return of FTSE100 Index portfolio is 0.0937%, the number of asset whose investment proportion more than 2% is 16, and the maximum investment proportion of portfolio is 7.3%.Compared with the FTSE100 Index, the optimal portfolio of εBBO is inclined to the concentrated investment and avoid the small amount of investment.Besides, once the short-selling is allowed, the return of the optimal portfolio has been greatly improved, and the capital is tending to focus on some assets with excellent performance.It can be seen from the experimental results that the εBBO algorithm is less efficient than its performance in the low dimension.
Generally speaking, the εBBO algorithm gets an excellent performance.However, the εBBO algorithm has a large chance of searching the optimal solution.The multi-directionality of the mutation and migration leads to the inability to guarantee the positive optimization of the solution.Besides, we can increase the number of initial populations to reduce the number of iterations but it will increase the performance period of the algorithm.From the performance of the optimal portfolio, the optimal portfolio increase the risk of portfolio which can not guarantee the stability of the benefits in each period.

Conclusions and Future Research
The second order stochastic dominance poses a great challenge to the numerical processing ability of the algorithm because of its unique judgment method.Based on the BBO algorithm, in this paper we discuss the portfolio optimization problem with the second order stochastic dominance constraints by using the improved εBBO algorithm, and give a detailed numerical experiment.The experimental results show that the εBBO algorithm is effective on portfolio optimization problem.In this paper, we prove the feasibility in the field of economic decision-making and explore the practice scope of BBO algorithm.
In this paper, we mainly discuss the simple target portfolio optimization problem.In the future, we will try the multiple target portfolio optimization problem.The target is not only the return of portfolio, but also can be the investment cycle, conditional value at risk (CVaR) and so on.Besides, in terms of transaction costs, we will discuss different types of transaction cost functions and introduce them into the second order stochastic dominance model.Furthermore, we also try to combine the εBBO algorithm with some prediction algorithms, such as Neural network algorithm, to predict the performance of the portfolio in the future.

Figure 1 .
Figure 1.The procedure chat of the εBBO algorithm.

Figure 2 .
Figure 2.The performance of the εBBO and the optimal portfolio in X ∈ [0, 0.6] of model(28).

1 . 18 (
a) The change of the the return of optimal portfolio between different iterative

Figure 5 .
Figure 5.The performance of the εBBO and the optimal portfolio in X ∈ [0, 0.6] of model (29).

Figure 6 .
Figure 6.The performance of the εBBO and the optimal portfolio in X ∈ [0, 1] of model (29).

Figure 11 .Figure 12 .Figure 13 .
Figure 11.The return of optimal portfolio by εBBO in each period with benchmark portfolio and the FTSE100 Index portfolio in X ∈ [0, 0.1] of model (30).

Table 1 .
Rates of return on five assets over ten period.

Table 2 .
(28)model(28)using date in Table1.Time is in minute and the expected return of the benchmark portfolio E[g(y, ξ)] = 1.073.

Table 3 .
The model (29) using date in Table1.Time is in minute and the expected return of the benchmark portfolio E[g(y, ξ)] = 1.073.

Table 4 .
The model (30) using date of 249 daily historical returns of 101 FTSE100 assets prior to December 2016.Time is in minute and the expected return of the benchmark portfolio E[g(y, ξ)] = 0.0595.