Research Article  Open Access
Q. H. Zhai, T. Ye, M. X. Huang, S. L. Feng, H. Li, "Whale Optimization Algorithm for Multiconstraint SecondOrder Stochastic Dominance Portfolio Optimization", Computational Intelligence and Neuroscience, vol. 2020, Article ID 8834162, 19 pages, 2020. https://doi.org/10.1155/2020/8834162
Whale Optimization Algorithm for Multiconstraint SecondOrder Stochastic Dominance Portfolio Optimization
Abstract
In the field of asset allocation, how to balance the returns of an investment portfolio and its fluctuations is the core issue. Capital asset pricing model, arbitrage pricing theory, and Fama–French threefactor model were used to quantify the price of individual stocks and portfolios. Based on the secondorder stochastic dominance rule, the higher moments of return series, the Shannon entropy, and some other actual investment constraints, we construct a multiconstraint portfolio optimization model, aiming at comprehensively weighting the returns and risk of portfolios rather than blindly maximizing its returns. Furthermore, the whale optimization algorithm based on FTSE100 index data is used to optimize the above multiconstraint portfolio optimization model, which significantly improves the rate of return of the simple diversified buyandhold strategy or the FTSE100 index. Furthermore, extensive experiments validate the superiority of the whale optimization algorithm over the other four swarm intelligence optimization algorithms (gray wolf optimizer, fruit fly optimization algorithm, particle swarm optimization, and firefly algorithm) through various indicators of the results, especially under harsh constraints.
1. Introduction
In the field of finance, the portfolio optimization problem has drawn a lot of attention since the meanvariance (MV) model was proposed by Markowitz [1]. Based on the probability theory, the MV model takes the first two moments of the return rate distribution, including expected returns and variance, into account. However, subsequent research studies find that the MV model is not suitable for the practical financial environment. For example, Chunhachinda et al. [2] pointed out that the returns to the major international stock market are not normally distributed, while the MV model assumes that the expected returns of the portfolio have a symmetric normal distribution. Besides, the MV model takes the variance as a risk measure, which counts both upward and downward deviations, which is contrary to the definition of the investment risk [3]. Hence, Markowitz [4] replaced the risk measure with the semivariance, which is more suitable in the case of asymmetric distribution. In addition, there are many other alternative risk measures, such as the mean absolute deviation (MAD), value at risk (VaR), and conditional value at risk (CVaR).
Other than the traditional MV model, the portfolio optimization model based on the stochastic dominance (SD) relation is highly applied. To be specific, the relationship between the firstorder stochastic dominance (FSD) and utility theory was first discussed by Quirk and Saposnik [5]. Hadar and Russell [6] extended the above relationship to the secondorder stochastic dominance (SSD). Thereinto, differentorder SD relations relate to different kinds of utility functions. In the meantime, SD relation is of interest because it takes the risk appetite into consideration and is more suitable to investors in realistic financial environments. Specifically speaking, FSD reflects the behavior of rational investors, while SSD reflects the behavior of rational riskaverse investors. Dentcheva and Ruszczynski [7] first proposed the portfolio optimization model with SD constraints. Then, Roman and Mitra [8] pointed out that the portfolio optimization model with the SSD rule improves the performance of the traditional MV model. Besides, Leshno and Levy [9] established an almost stochastic dominance (ASD) relation, which reveals the preference for most investors rather than all of them. Moreover, Fabian et al. [10] proposed the SSD portfolio optimization model, which is called TSSD, using tail risk measures at different confidence levels. Recently, Javanmardi and Lawryshyn [11] proposed the SSDDP model, which does not require a benchmark portfolio.
However, only considering the SSD constraint is not enough to get an effective portfolio in the complex investment environment. Therefore, additional reallife constraints, such as transaction costs, higher moments, diversification, and boundary constraints, should be taken into consideration in portfolio optimization. In any financial market, investors are obliged to pay transaction costs when buying or selling the securities. As one of the key factors affecting the net return, transaction costs are of great importance to investors, and ignoring them will result in an ineffective portfolio [12, 13]. Besides, higher moments of return distribution, including skewness and kurtosis, have drawn investors’ attention, which takes the extreme value into account rather than focusing solely on the average [14]. Leung et al. [15] pointed out that neglecting skewness will lead to an inefficient portfolio. Yu and Lee [16] showed that the model with higher moments performs better. As for diversification, it is a strategy that integrates a wide variety of investments into the portfolio. Models considering diversity are propitious to eliminate the unsystematic risk in the portfolio and make asset allocation more feasible [17, 18]. Normally, Shannon’s entropy and Yager’s entropy are used to measure the diversity of a portfolio. Furthermore, given that small proportions held in the portfolio have little impact on its performance and high proportions go against the flexibility, the boundary constraint, also known as the buyin threshold constraint, is included in the portfolio optimization model to limit the upper and lower bounds of the investment ratio in individual assets [19, 20].
After introducing several constraints into the model, the portfolio optimization problem converts to an NPcomplete problem, which limits the computational efficiency as the problem size increases. Besides, SSD constraints require the comparison of any two investable assets. As the number of investable assets increases, the computational requirements become more demanding. Therefore, heuristic algorithms are widely used to optimize the multiconstraint portfolio optimization problem, such as the particle swarm optimization (PSO) [21], firefly algorithm (FA) [22], biogeographybased optimization (BBO), and artificial bee colony (ABC) [23]. Recently, Babazadeh and Esfahanipour [24] developed a nondominated sorting genetic algorithm and applied it to a multiperiod meanvar portfolio optimization model under the cardinality, budget, floor, and ceiling constraints. Chen et al. [25] combined the FA and the genetic algorithm (GA) and then applied it to a meanvarianceskewness portfolio selection model under the transaction costs, bounds on holdings, cardinality, and transaction lots constraints.
Whale optimization algorithm (WOA) is a metaheuristic optimization algorithm introduced by Mirjalili and Lewis [26], which is inspired by the bubblenet hunting strategy of humpback whales. Besides, WOA was tested for 29 mathematical optimization problems and 6 structural design problems, and the results prove that WOA is very competitive with the existing metaheuristic algorithms as well as conventional methods. The most important is that WOA mainly mimics the hunting behavior of humpback whales in searching for and attacking preys called the bubblenet feeding behavior [27], which improved the performance of candidate solutions in each step. Because of its characteristics of simple theory, easier operation, less parameter setting, and no special requirements for optimized function, WOA is widely used in image segmentation [28], parameter estimation [29], sizing optimization [30], and global optimization [31]. Recently, AbdelBasset et al. [32] integrated WOA with a local search strategy for solving the problem of permutation flow shop scheduling. Reddy et al. [33] applied WOA to profitbased unit commitment problems in competitive electricity markets with good performance.
However, in these papers, WOA has still not been applied to solving the multiconstraint SSD portfolio optimization model. The motivation for proposing WOA for the multiconstraint SSD portfolio optimization model in this research is twofold. First, based on the above literature review, there have been several successful applications based on WOA. Second, in the preliminary research work, we have achieved decent performance on the SSD portfolio optimization model by GWO. In this paper, a new improved multiconstraint SSD portfolio optimization model is proposed and optimized by WOA.
For the multiconstraint SSD portfolio optimization model, the area of the feasible region is very small, which makes it difficult to produce feasible solutions. At the same time, considering harsh constraints makes the traditional optimization methods tend to converge to the local optimum, which results in the loss of population diversity and poor optimization ability. Compared with the traditional evolutionary algorithm, such as GA, on the one hand, WOA has a stronger ability to jump out of the local optimum and explore the globally optimal solution; on the other hand, the bubblenet attack process ensures that WOA has higher accuracy than other swarm intelligence algorithms, such as the gray wolf optimizer (GWO) and the fruit fly optimization algorithm (FOA). Therefore, WOA is used to optimize the above multiconstraint SSD portfolio optimization model. Furthermore, we perform numerical experiments based on FTSE100 index data and compare the performance of the optimal portfolio obtained by WOA with that obtained by GWO, FOA, PSO, and FA.
The novel contributions of this paper are as follows.
Firstly, the multiconstraint SSD portfolio optimization model is proposed. We incorporate several realistic constraints into the SSD portfolio optimization model, including the transaction cost, skewness, kurtosis, diversification, and boundary constraints. Secondly, WOA is successfully used to optimize the above multiconstraint SSD portfolio optimization model. Lastly, an optimal portfolio strategy obtained by WOA is found to be greater than that obtained by other four different algorithms. The rest of this paper is organized as follows. In Section 2, we discuss the portfolio optimization model with SSD and other reallife constraints. In Section 3, we discuss the WOA for solving the proposed model. In Section 4, we present a numerical experiment and analyse the performance of WOA. Finally, a brief conclusion is illustrated in Section 5.
2. The Multiconstraint Portfolio Optimization Model
In this section, a detailed discussion of the SSD portfolio optimization model under several realistic or higher moment constraints is presented. Specifically speaking, it is well known that only considering returns and risk is not enough to get an effective portfolio that meets the complicated investment environment. Besides, in the realworld financial market, the return distribution of the portfolio is not a normal distribution, which means that higher moments of return distribution should not be neglected. Therefore, several realistic or higher moment constraints, including the transaction cost, skewness, kurtosis, diversification, and boundary constraints, are included in the framework of SSD portfolio optimization.
2.1. The Returns, Net Return, and Excess Mean Return
As one of the most basic elements of portfolio optimization problems, returns play an important part in indicating the performance of a portfolio. In the financial market, the returns of investment can be defined as the ratio of the net gain or loss to the initial cost of the investment, separately known as positive returns and negative returns, which can be formulated as follows:where and are the price of asset at period and , respectively, and is the dividend of asset and period . Let denote the number of assets, which are available for investment at the beginning of a fixed period and denote the fraction of the initial capital invested in . Then, we use to denote a set of feasible portfolios, which is a bounded convex polyhedron clearly. Besides, let denote the returns of asset under discrete distribution, where a random vector is on the probability space . Assuming that the capital available for investment is fixed and , the returns of portfolio can be expressed as follows [3]:
However, we always need to pay fees when trading assets, such as brokerage fees, bidask spreads, taxes, and fund load, which are known as transaction costs. Arnott and Wanger [12] pointed out that ignoring the transaction costs would lead to an ineffective portfolio. In this paper, we employ a Vshaped transaction cost function to express the transaction costs between the new portfolio and the existing portfolio . Let denote the total transaction costs of assets, which can be formulated as follows:where is the unit transaction cost of asset . It is set to 0.9%, among which the commission is 0.4% and the stamp duty is 0.5%. Thus, the net return of portfolio can be displayed as follows [34]:
Furthermore, the excess mean return (EMR) is always used to describe performance over the existing portfolio, which is an average of the difference between the returns of the new portfolio and the existing portfolio , and can be formulated as follows:where and denote the returns of the portfolio and in the period . In a word, all of the above three indicators can be used to evaluate the performance of a new portfolio. Specifically speaking, the returns of portfolio are used to calculate constraints, such as skewness. The net return of the portfolio is used as an objective function, and the EMR is only used as an evaluation indicator.
Fama and French [35] put forward a threefactor pricing model and recently put forward a fivefactor pricing model [36]. It builds its portfolio at the end of June and holds it for a year, based on financial data from the previous year. In order to simplify the rebalancing problem, this paper only considers the construction of a singleperiod portfolio and assumes that the transaction cost required to construct each portfolio is fixed. Therefore, the returns of portfolio are used as the objective function of the optimization model in this paper.
2.2. Risk Measures and Stochastic Dominance Constraint
As another basic element of portfolio optimization problems, risk presents the uncertainty of returns. The risk measures have many forms, including variance, semivariance (SV), MAD, VaR, and CVaR. VaR is one of the most wellknown downside risk measures, which measures the worst returns a portfolio may potentially suffer. With a fixed confidence level , is defined as the quantile of the cumulative distribution function, which can be defined as follows:where represents the returns of the portfolio in each period. Meanwhile, as a coherent risk measure, CVaR measures the conditional expectation of losses beyond VaR, which can be displayed as follows:
Unlike the above risk measures, the SD constraint is not an indicator, which establishes relative risk advantages by comparison between the korder distribution function of any two portfolios. Assuming that the utility function of all investors is monotonically increasing, the portfolio stochastically dominates the portfolio in the first order when all the investors prefer portfolio to portfolio or that there is no difference between a part of them [37]. Let and be the decision vectors and be a random variable; then is preferred to weakly in , firstorder stochastic dominance, denoted by , if and only ifwhere is the returns of the portfolio , which is a concave continuous function both in and , and is the cumulative distribution function of . Similarly, is preferred to weakly in the SSD, denoted by , if and only if
Therefore, the strict dominance relation succ_{k} is defined as follows:
There are several equivalent characterizations of the SSD constraint. Hadar and Russell [6] pointed out that for any nondecreasing and concave utility function , if and only ifwhere and are the two random variables, which generally represent the returns of portfolios and , and is the expected value with respect to the probability distribution of [10]. Ogryczak and Ruszczyński [38] pointed out that if and only ifwhere . Ogryczak and Ruszczyński pointed [39] out that is equivalent to the continuum of CVaR constraints for all confidence levels :
In this paper, we use the last equivalent SSD relation, VaR, and CVaR as evaluation indicators.
2.3. Higher Moment Constraints
The traditional MV model formulated by Markowitz [1] only takes the first two moments of return distribution into account for portfolio optimization. Arditti and Levy [40] and Rubinstein [41] have argued that the higher moments of return distribution should not be neglected because the distribution of returns in the financial market is not a normal distribution. From the theoretical and empirical point of view, Arditti [42, 43] proved that investors demand higher (lower) returns for investments with negative (positive) skewness of income distribution. Further, Scott and Horvath [44] extended this analysis to the higherorder moments of the return distribution and proved that the positive values of even (odd) order moments bring a positive (negative) risk premium, and vice versa. Referring to the recent studies by Chen and Yue et al. [25, 45], the skewness and kurtosis constraints are incorporated into the optimization model in this paper.
Let be an uncertain variable with a finite expected value , and the skewness and kurtosis of are, respectively, defined by
2.4. Diversification Constraint
In the financial market, portfolio diversification is the process of allocating the capital in a way that reduces the exposure to any one particular asset or risk, which implies that the idiosyncratic risk of the portfolio can be reduced as the assets included in the investment increase. Considering the low diversity of the portfolio which may lead to losses, we include the diversification constraint into the portfolio optimization model. SE and Yager’s entropy are widely used to measure the diversity of a portfolio [17, 46]. Let denote the weight of assets , and SE can be formulated as follows:and Yager’s entropy can be calculated as follows:where is a constant, and in this paper, it is set to 2. Obviously, the higher the value of SE or Yager’s entropy, the better the diversity of the portfolio. In particular, both SE and Yager’s entropy get the maximum value when .
2.5. Boundary Constraint
The boundary constraint, also known as the buyin threshold constraint, means that each asset should be invested in a specific range, while the lower bounds are used to reduce the brokerage costs and monitoring costs and upper bounds are used to increase the flexibility [19]. Let denote the lower bounds and denote the upper bounds; then the boundary constraint of the portfolio can be described as follows:where is the weight of assets and ; if , ; otherwise, .
2.6. Budget Constraint
Budget constraint means that the amount of capital to be invested is fixed and all capital should be invested, which can be expressed as follows:
2.7. No Short Selling Constraint
The short selling refers to the sale of securities borrowed by the seller, which is an act of speculation with high risk in the financial market. Therefore, the short selling is not considered in this paper, and the no short selling constraint is represented by
2.8. Multiconstraint SecondOrder Stochastic Dominance Portfolio Optimization Model
Given all that, the multiconstraint SSD portfolio optimization model is provided below:
3. The Whale Optimization Algorithm (WOA) for Solving the Proposed Model
3.1. Background
WOA is a swarm intelligence optimization algorithm proposed by Australian scholars Mirjalili and Lewis [26] in 2016. The algorithm simulates the process of humpback whales searching and capturing food by establishing mathematical models. WOA was inspired by the bubblenet attack strategy of humpback whales. Whales surround their prey by spiraling up and bubbling up as they dive about 12 meters below the surface. WOA has the advantages of a simple principle, few parameters, and strong searching ability. Since its inception, it has been widely used in engineering optimization, parameter extraction, feature selection, and other such aspects. Figure 1 shows the bubblenet feeding behavior of humpback whales.
3.2. Mathematical Model
In WOA, the position of the ith whale (search agent) is described as , where is required to be given in advance, and the size of is equal to dim (the dimension of the problem). The steps of the WOA are as follows.
First, initial search agents with positions are generated. Then, the position and the score of the optimal search agent with position are initialized. For solving the minimum problem, the bestscore (the initial score of the optimal search agent) is set to be instead of .
The algorithm loop is entered. For each search agent, each variable is checked whether or not it is across the border. For a variable that is out of bounds, its value is returned to the boundary.
Real numbers and are randomly generated and and are used to calculate and . The calculation formulas of and are as follows:where is a real number that goes linearly from 2 to 0 in the iteration.
3.2.1. Random Prey
If and , the stage of random search for the prey is entered. In this process, the search agent is forced to move away from its current location and wander randomly through the space in search of the prey. The mathematical model and position transformation formula are as follows:where is a randomly generated position vector within the boundary range, is the generation of the ith search agent’s position vector, and is the ith search agent’s position vector generation .
3.2.2. Encircling Prey
If and , the stage of surrounding the prey is entered. In this process, the location of the prey is identified, and the prey is surrounded. The search agent moves closer to the location of the optimal search agent. The mathematical expression of its position update iswhere is the generation t of optimal search agent’s position vector.
3.2.3. BubbleNet Attacking
If , the stage of spiral contraction encircling and bubblenet attack is entered. Firstly, the distance between the search agent and the optimal search agent is calculated, and then a spiral mathematical model is established to update the search agent’s position. The formula is as follows:where , is the defined logarithmic helix shape constant, and is a random real number in the range .
3.3. The WOA for Solving the Proposed Model
When it comes to WOA for solving this model, we need to add judgment constraints into the WOA iteration process. That is, the search agent that meets the constraint and has better fitness is allowed to be updated to the leading search agent. There are several constraints in the model. In order to ensure the quality of the final solution, we added the initialization coefficient in the initialization process, which makes the quality of the initial solution reach our desired effect. Figure 2 shows the procedure of the WOA algorithm. Meanwhile, the pseudocode of the main procedure of WOA is shown in Algorithm 1, and the pseudocode of WOA for solving the proposed model is shown in Algorithm 2.


4. Numerical Experiments
In this section, we carry out several numerical experiments based on FTSE100 index stock historical data. Specifically speaking, the historical return rate of FTSE 100 index assets prior to December 2018 is collected to construct the portfolio strategy. Besides, a series of indicators are introduced to evaluate the performance of portfolios by Goel and Sharma [47], such as EMR, downside deviation (DD), and Sharpe ratio. Moreover, the algorithms are coded in MATLAB 2016a, and all tests are performed on a PC with a Windows 10 operating system and 8 GB of RAM.
4.1. Performance Measures
We mainly evaluate a portfolio from two aspects: absolute indicator and relative indicators, which are described briefly as follows.
As a measure of the downside risk, DD focuses on the returns below the minimum threshold or the minimum acceptable returns. Let and denote the returns of portfolios and in the period ; then DD can be formulated as follows:where is the FTSE100 index in this paper, and the lower value of DD means a better performance of portfolio . Besides, DD is also used to calculate the Sortino ratio. Sortino ratio is a variation of the Sharpe ratio, which uses DD rather than the standard deviation as a risk measure. The Sortino ratio and the Sharpe ratio can be calculated as follows:where denotes the net return of portfolio , denotes the riskfree rate, and denotes the standard deviation. Obviously, a higher value of the Sortino ratio or the Sharpe ratio is desirable. Moreover, the STARR ratio (STARR) and the information ratio (IR) are alternatives to the Sharpe ratio. STARR also considers the major drawback of the standard deviation as a risk measure and employs the CVaR for the performance adjustment, which is defined aswhere is the confidence level. Considering that CVaR is usually negative, therefore, the lower value of STARR is preferable. Furthermore, IR is a measure of the returns on a portfolio beyond a benchmark, and the benchmark is typically an index, such as the FTSE100 index, which is given by
Similarly, the higher value of IR is desirable. In this paper, the confidence level is set to 5%, the riskfree rate is set to 4%, and the FTSE100 index is used as a benchmark.
4.2. Numerical Results
In this section, we report the detailed experimental results. Specifically, the return rate of FTSE100 index stocks during 2018 is collected as historical data. Besides, in order to ensure the diversity of the portfolio, the upper bound δ of the portfolio is set to 5%, 7%, or 10%. Above all, we propose the following multiconstraint portfolio optimization model under SSD:
To evaluate the performance of the WOA, we compare it with the standard GWO, FOA, PSO, and FA. In order to make the running time roughly equal, except for FA, the number of initial solutions for the other four algorithms is 3000 (in FA, it is set to 1500). And, the initialization coefficient . In WOA, GWO, and PSO, the number of iterations is 30, while it is 20 in FOA and FA. The experiments are repeated 30 times.
Comparing with the FTSE100 index, we get the experiment results as shown in Tables 1–3 under δ = 0.05, 0.07, and 0.10, respectively, where the “Mean” column describes the average indicator value of optimal portfolios during repeated experiments, and the “Optimal” column describes the indicator value of optimal portfolios among repeated experiments. For the simplicity of data representation, the order of magnitude of EMR in the table is 10^{−4}. The order of magnitude of DD is 10^{−2}.



In Figures 3(a)–3(c), the histograms, respectively, show the specific asset structure of the optimal portfolio of 101 kinds of assets when the model upper bound = 0.05, 0.07, and 0.10 by the WOA algorithm. The histogram in Figure 3(d) shows the specific asset structure of FTSE100. From Figures 4–21, it can be seen that the line charts show the daily rate of return of the optimal portfolio of WOA, GWO, FOA, PSO, and FA when the model upper bound = 0.05, 0.07, and 1.10 and compared with the daily rate of return of FTSE100.
(a)
(b)
(c)
(d)
Figures 22–24 show the box plots of the net return, skewness, and Shannon’s entropy obtained by five different algorithms when the model upper bound = 0.05, 0.07, and 0.10, respectively.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
4.3. Numerical Analysis
We perform the numerical experiments of the multiconstraint portfolio optimization model under SSD by WOA, GWO, FOA, PSO, and FA. In this part, the performance of the portfolio obtained by each algorithm is evaluated through returns, diversity, risk, relative risk value, and other indicators.
First of all, it can be seen from Table 1 that the mean net return of portfolios obtained by WOA is 0.06737, while the portfolios obtained by GWO, FOA, PSO, and FA are 0.05907, 0.04590, 0.05327, and 0.05043, respectively, under δ = 0.05. Besides, the mean net return and EMR of portfolios obtained by WOA under different upper bounds are also higher than portfolios obtained by the other four algorithms. If only considering the optimization results, it is obvious that WOA has a better optimization capability than the other four algorithms under strong constraints. Although the mean SE of the portfolio obtained by WOA is slightly lower than the value obtained by FOA, the mean Sharpe ratio of FOA is much lower than the performance obtained by WOA. For investors, they tend to choose the portfolio with a higher Sharpe ratio when there is a little difference in SE respect. Meanwhile, the primary goal of investors is to obtain the maximum returns, while the returns of the portfolio obtained by FOA are far lower than that obtained by WOA. Taken together, FOA’s results are unsatisfactory in that they give us a slight advantage in diversity at the expense of returns.
In Table 2, when δ = 0.07, the mean and optimal net returns of the portfolio obtained by WOA are higher than those obtained by the other four algorithms. In Table 3, the performance obtained by WOA is superior to that obtained by the other four algorithms under δ = 0.10. The numerical analysis results in Table 2 and Table 3 are similar to those in Table 1.
Secondly, it can be seen from Figures 3(a)–3(c) the 101 asset structure of the optimal portfolio obtained by the WOA algorithm under δ = 0.05, 0.07, and 0.10. Figure 3(b) shows the asset portfolio obtained by FTSE100. The portfolio optimized by WOA has a more even distribution of assets, and its upper bound is lower than the upper bound of FTSE100, which proves that its diversity is better.
In Figures 4–21, the horizontal axis represents 253 trading days and the vertical axis represents the return rate. These figures show the daily rate of return of the optimal portfolio of WOA, GWO, FOA, PSO, and FA when δ = 0.05, 0.07, and 0.10 and compare it with the daily rate of return of FTSE100. It can be seen that the return rate of the portfolio optimized by these five different algorithms is mostly higher than the return rate of FTSE100 under different values of δ. From Figures 4–9, it is seen that when δ = 0.05, the return rate of the portfolio obtained by the five different algorithms is mostly better than the return rate of FTSE100 in 253 trading days. From Figures 10–15, it is seen that the portfolio return rate obtained by the five different algorithms is mostly better than the return rate of FTSE100 under δ = 0.07. From Figures 16–21, it is seen that the optimized portfolio return rate is mostly better than the return rate of FTSE100 under δ = 0.10. From Figures 4, 10, and 16, for approximately 80% of the trading days, the return rate line of the portfolio obtained by WOA (red line) is higher than the FTSE100’s return rate line (blue dotted line). Compared with the return rate line of the portfolio obtained by other four kinds of algorithm, it can be proved that the returns of the portfolio optimized by WOA are higher than that of the other four algorithms.
Further, from Figures 22–24, we can see the box plots of net return, Shannon’s entropy, and skewness of the portfolios optimized by the five different algorithms. It can be seen that when δ = 0.05, 0.07, and 0.10, respectively, the mean returns of the portfolio optimized by WOA are higher than that of the other four algorithms. For investors, they tend to choose the portfolio with a higher net return, so the portfolio optimized by WOA is a useful guideline to investors.
Although the mean Shannon’s entropy of the portfolio obtained by WOA is slightly lower than the value obtained by the other four algorithms, the mean net return of WOA is the highest of these algorithms. It means that in order to achieve a higher net return, a part of the diversity needs to be sacrificed.
When δ = 0.05, the mean skewness of the portfolio optimized by WOA is only 0.0124 lower than the maximum of that obtained by the other four algorithms. When δ = 0.07 and 0.10, the portfolios optimized by WOA have a greater skewness than the other algorithms. At the same time, the mean net return of WOA is the highest of these algorithms. For investors, they prefer the portfolio with larger skewness and higher returns. Therefore, the portfolio obtained by WOA is more suitable for investors to choose.
In addition, besides using FTSE100 as the comparison criterion, equally weighted indicators are also used as the comparison object. When investors build portfolios, more reference schemes are provided. The indicators of FTSE100 and equal weight are also given in Table 4.
From Table 4, it can be seen that the net return under FTSE100 and equal weight is −0.1038 and −0.1057, which is lower than that optimized by these algorithms. Meantime, their skewness is all lower than that optimized by these algorithms. Their kurtosis is greater than 3, indicating that the solutions present a spike pulse. Their kurtosis is all higher than that optimized by these algorithms. By comparing with the performance indicators in the case of FTSE100 and equal weight, it indicates that the performance of the optimized portfolio by WOA has a greater reference value for investors.

The specific index weight data of the FTSE100 index are shown in Table 5. In this paper, we use the Sortino ratio, STARR ratio, Sharpe ratio, and IR to measure the relative risk value, which is used to quantify the net return or the excess mean return from taking on the unit risk. When the upper bound is set to 0.05, it can be seen from Table 1 that the mean Sortino ratio of portfolios obtained by WOA is 0.3004, while the portfolios obtained by GWO, FOA, PSO, and FA are 0.2930, 0.2732, 0.2817, and 0.2746, respectively. Besides, the mean STARR ratio, Sharpe ratio, and IR of portfolios obtained by WOA are also better than the other four algorithms, which means that the portfolio obtained by WOA has a better relative risk value. Above all, under several constraints, WOA has an efficient search ability to find the optimal portfolio that meets the constraints. In addition, in terms of the algorithm structure, WOA also has the advantages of a simple structure and few parameter settings. From the upper bound of results, the upper bound of WOA is higher than other algorithms, which indicates that WOA has a good convergence effect and strong local search ability.

Furthermore, the SSD optimization model proposed in this paper has a strong constraint on the solution. In other words, the existence of constraints will lead to infeasible regions in the search space of decision variables. In fact, in the field of finance, realistic constraints are pretty important. The existence of constraints will lead to infeasible regions in the search space of decision variables. The algorithm needs to balance constraints and optimization. Under strong constraints, the area of the feasible region is very small, and a large number of solutions are considered to be infeasible because they do not meet the constraints, which makes it a very difficult problem to produce feasible solutions. At the same time, the existence of constraints may also make the original problem produce many new local optima [48]. In this case, the optimization algorithm is easy to converge too quickly, resulting in the loss of population diversity in the early stage, and the global optimal solution cannot be found. Therefore, the strong constraint optimization problem has a high demand on the ability of the algorithm to jump out of the local optimum and the ability of global search [49]. According to the experimental results in this section, under different upper bounds, the value of the objective function obtained by WOA is better in both the optimal conditions [50, 51]. This shows that the performance of WOA is better than that of the other four algorithms in strongly constrained optimization problems, which has an excellent application prospect.
5. Conclusion
Since Fama and French [35] put forward the threefactor pricing model, there has been a boom in searching for factors in academic circles recently. Factor model is essentially an asset pricing model and can predict the future returns on assets to a certain extent. However, no matter how many factors there are, we should strive to achieve higher returns on the premise of minimizing portfolio risks, namely, asset allocation. With the development of machine learning, more and more advanced algorithms are used to explore the nonlinear relationship between factors and returns on assets. Similarly, machine learning can be applied to asset allocation.
In this paper, we propose a portfolio optimization model under SSD and several realistic constraints, and it is optimized by WOA. As a matter of fact, the real financial environment is much more complex. Therefore, how to broaden the practicability of the intelligent algorithm in a complex financial environment is still a long way to go. Based on FTSE100 index stock data, the experimental results prove the outstanding performance of WOA during harsh constraints.
However, the proposed portfolio optimization model only considers a singleperiod problem. In fact, portfolio construction is a continuous process. In the future work, establishing a multiperiod portfolio model and verifying the performance of WOA or other intelligent algorithms are feasible. Moreover, WOA has a comparatively huge rise space in global searching, and we find that the WOA has the problem of premature convergence. We will pay attention to this topic in the future research.
Abbreviations
WOA:  Whale optimization algorithm 
GWO:  Gray wolf optimizer 
FOA:  Fruit fly optimization algorithm 
PSO:  Particle swarm optimization 
FA:  Firefly algorithm 
EMR:  Excess mean return 
DD:  Downside deviation 
STARR:  STARR ratio 
IR:  Information ratio 
SE:  Shannon’s entropy 
YE:  Yager’s entropy 
VaR:  Value at risk 
CVaR:  Conditional value at risk 
SD:  Stochastic dominance 
FSD:  Firstorder stochastic dominance 
SSD:  Secondorder stochastic dominance 
ASD:  Almost stochastic dominance 
MV:  Meanvariance. 
Data Availability
The composition of FESE100 is obtained from the London stock exchange (http://www.londonstockexchange.com). The daily closing price of the constituent shares and dividends is obtained from Google Finance (http://www.google.com/finance).
Conflicts of Interest
The authors declare no conflicts of interest.
Authors’ Contributions
Qianhao Zhai and Tao Ye contributed equally to this work.
Acknowledgments
The authors gratefully acknowledge the support of the National Key Research and Development Plan of China (no. 2018YFB1404400) in carrying out this research.
References
 H. Markowitz, “Portfolio selection,” The Journal of Finance, vol. 7, no. 1, pp. 77–91, 1952. View at: Publisher Site  Google Scholar
 P. Chunhachinda, K. Dandapani, S. Hamid, and A. J. Prakash, “Portfolio selection and skewness: evidence from international stock markets,” Journal of Banking & Finance, vol. 21, no. 2, pp. 143–167, 1997. View at: Publisher Site  Google Scholar
 Y. Ren, T. Ye, M. Huang, and S. Feng, “Gray wolf optimization algorithm for multiconstraints secondorder stochastic dominance portfolio optimization,” Algorithms, vol. 11, no. 5, p. 72, 2018. View at: Publisher Site  Google Scholar
 H. M. Markowitz, Portfolio Selection: Efficient Diversification of Investments, John Wiley & Sons, New York, NY, USA, 1959.
 J. P. Quirk and R. Saposnik, “Admissibility and measurable utility functions,” The Review of Economic Studies, vol. 29, no. 2, pp. 140–146, 1962. View at: Publisher Site  Google Scholar
 J. Hadar and W. R. Russell, “Rules for ordering uncertain prospects,” American Economic Review, vol. 59, no. 1, pp. 25–34, 1969. View at: Google Scholar
 D. Dentcheva and A. Ruszczynski, “Portfolio optimization with stochastic dominance constraints,” Journal of Banking & Finance, vol. 30, no. 2, pp. 433–451, 2006. View at: Publisher Site  Google Scholar
 D. Roman and G. Mitra, “Portfolio selection models: a review and new directions,” Wilmott Journal, vol. 1, no. 2, pp. 69–85, 2009. View at: Publisher Site  Google Scholar
 M. Leshno and H. Levy, “Preferred by “all” and preferred by “most” decision makers: almost stochastic dominance,” Management Science, vol. 48, no. 8, pp. 1074–1085, 2002. View at: Publisher Site  Google Scholar
 C. I. Fabian, G. Mitra, and D. Roman, “Processing secondorder stochastic dominance models using cuttingplane representations,” Mathematical Programming, vol. 130, no. 1, pp. 35–57, 2011. View at: Publisher Site  Google Scholar
 L. Javanmardi and Y. Lawryshyn, “A new rank dependent utility approach to model risk averse preferences in portfolio optimization,” Annals of Operations Research, vol. 237, no. 12, pp. 161–176, 2016. View at: Publisher Site  Google Scholar
 R. D. Arnott and W. H. Wagner, “The measurement and control of trading costs,” Financial Analysts Journal, vol. 46, no. 6, pp. 73–80, 1990. View at: Publisher Site  Google Scholar
 A. Yoshimoto, “The meanvariance approach to portfolio optimization subject to transaction costs,” Journal of the Operations Research Society of Japan, vol. 39, no. 1, pp. 99–117, 1996. View at: Publisher Site  Google Scholar
 M. K. Mehlawat, A. Kumar, S. Yadav, and W. Chen, “Data envelopment analysis based fuzzy multiobjective portfolio selection model involving higher moments,” Information Sciences, vol. 460–461, pp. 128–150, 2018. View at: Publisher Site  Google Scholar
 M. T. Leung, H. Daouk, and A.S. Chen, “Using investment portfolio return to combine forecasts: a multiobjective approach,” European Journal of Operational Research, vol. 134, no. 1, pp. 84–102, 2001. View at: Publisher Site  Google Scholar
 J.R. Yu and W.Y. Lee, “Portfolio rebalancing model using multiple criteria,” European Journal of Operational Research, vol. 209, no. 2, pp. 166–175, 2011. View at: Publisher Site  Google Scholar
 J.R. Yu, W.Y. Lee, and W.J. P. Chiou, “Diversified portfolios with different entropy measures,” Applied Mathematics and Computation, vol. 241, pp. 47–63, 2014. View at: Publisher Site  Google Scholar
 M. K. Mehlawat, P. Gupta, A. Kumar, S. Yadav, and A. Aggarwal, “Multiobjective fuzzy portfolio performance evaluation using data envelopment analysis under credibilistic framework,” IEEE Transactions on Fuzzy Systems, In press. View at: Publisher Site  Google Scholar
 O. Ertenlice and C. B. Kalayci, “A survey of swarm intelligence for portfolio optimization: algorithms and applications,” Swarm and Evolutionary Computation, vol. 39, pp. 36–52, 2018. View at: Publisher Site  Google Scholar
 P. Gupta, M. K. Mehlawat, S. Yadav, and A. Kumar, “A polynomial goal programming approach for intuitionistic fuzzy portfolio optimization using entropy and higher moments,” Applied Soft Computing, vol. 85, Article ID 105781, 2019. View at: Publisher Site  Google Scholar
 M. Corazza, G. Fasano, and R. Gusso, “Particle swarm optimization with nonsmooth penalty reformulation, for a complex portfolio selection problem,” Applied Mathematics and Computation, vol. 224, pp. 611–624, 2013. View at: Publisher Site  Google Scholar
 N. Bacanin and M. Tuba, “Firefly algorithm for cardinality constrained meanvariance portfolio optimization problem with entropy diversity constraint,” The Scientific World Journal, vol. 2014, Article ID 721521, 16 pages, 2014. View at: Publisher Site  Google Scholar
 D. Kumar and K. K. Mishra, “Portfolio optimization using novel covariance guided Artificial Bee Colony algorithm,” Swarm and Evolutionary Computation, vol. 33, pp. 119–130, 2017. View at: Publisher Site  Google Scholar
 H. Babazadeh and A. Esfahanipour, “A novel multi period meanVaR portfolio optimization model considering practical constraints and transaction cost,” Journal of Computational and Applied Mathematics, vol. 361, pp. 342–313, 2019. View at: Publisher Site  Google Scholar
 W. Chen, Y. Wang, P. Gupta, and M. K. Mehlawat, “A novel hybrid heuristic algorithm for a new uncertain meanvarianceskewness portfolio selection model with real constraints,” Applied Intelligence, vol. 48, no. 9, pp. 2996–3018, 2018. View at: Publisher Site  Google Scholar
 S. Mirjalili and A. Lewis, “The whale optimization algorithm,” Advances in Engineering Software, vol. 95, pp. 51–67, 2016. View at: Publisher Site  Google Scholar
 I. Aljarah, H. Faris, and S. Mirjalili, “Optimizing connection weights in neural networks using the whale optimization algorithm,” Soft Computing, vol. 22, no. 1, pp. 1–15, 2018. View at: Publisher Site  Google Scholar
 M. A. E. Aziz, A. A. Ewees, and A. E. Hassanien, “Whale optimization algorithm and mothflame optimization for multilevel thresholding image segmentation,” Expert Systems with Applications, vol. 83, pp. 242–256, 2017. View at: Publisher Site  Google Scholar
 D. Oliva, M. Abd El Aziz, and A. Ella Hassanien, “Parameter estimation of photovoltaic cells using an improved chaotic whale optimization algorithm,” Applied Energy, vol. 200, pp. 141–154, 2017. View at: Publisher Site  Google Scholar
 A. Kaveh and M. I. Ghazaan, “Enhanced whale optimization algorithm for sizing optimization of skeletal structures,” Mechanics Based Design of Structures and Machines, vol. 45, no. 3, pp. 345–362, 2017. View at: Publisher Site  Google Scholar
 I. N. Trivedi, J. Pradeep, J. Narottam, K. Arvind, and L. Dilip, “Novel adaptive whale optimization algorithm for global optimization,” Indian Journal of Science and Technology, vol. 9, no. 38, pp. 319–326, 2016. View at: Publisher Site  Google Scholar
 M. AbdelBasset, G. Manogaran, D. ElShahat, and S. Mirjalili, “A hybrid whale optimization algorithm based on local search strategy for the permutation flow shop scheduling problem,” Future Generation Computer Systems, vol. 85, pp. 129–145, 2018. View at: Publisher Site  Google Scholar
 K. S. Reddy, L. Panwar, B. K. Panigrahi, and R. Kumar, “Binary whale optimization algorithm: a new metaheuristic approach for profitbased unit commitment problems in competitive electricity markets,” Engineering Optimization, vol. 51, no. 3, pp. 369–389, 2019. View at: Publisher Site  Google Scholar
 W. Chen and W. Xu, “A hybrid multiobjective bat algorithm for fuzzy portfolio optimization with realworld constraints,” International Journal of Fuzzy Systems, vol. 21, no. 1, pp. 291–307, 2019. View at: Publisher Site  Google Scholar
 E. F. Fama and K. R. French, “Common risk factors in the returns on stocks and bonds,” Journal of Financial Economics, vol. 33, no. 1, pp. 3–56, 1993. View at: Publisher Site  Google Scholar
 E. F. Fama and K. R. French, “Dissecting anomalies with a fivefactor model,” Review of Financial Studies, vol. 29, no. 1, pp. 69–103, 2016. View at: Publisher Site  Google Scholar
 T. Ye, Z. Yang, and S. Feng, “Biogeographybased optimization of the portfolio optimization problem with second order stochastic dominance constraints,” Algorithms, vol. 10, no. 3, p. 100, 2017. View at: Publisher Site  Google Scholar
 W. Ogryczak and A. Ruszczyński, “From stochastic dominance to meanrisk models: semideviations as risk measures,” European Journal of Operational Research, vol. 116, no. 1, pp. 33–50, 1999. View at: Publisher Site  Google Scholar
 W. Ogryczak and A. Ruszczynski, “Dual stochastic dominance and related meanrisk models,” SIAM Journal on Optimization, vol. 13, no. 1, pp. 63–78, 2002. View at: Publisher Site  Google Scholar
 F. D. Arditti and H. Levy, “Portfolio efficiency analysis in three moments: the multiperiod case,” Financial Dec Making under Uncertainty, Academic Press, New York, NY, USA, 1977. View at: Publisher Site  Google Scholar
 M. E. Rubinstein, “A comparative statics analysis of risk premiums,” The Journal of Business, vol. 46, no. 4, pp. 605–615, 1973. View at: Publisher Site  Google Scholar
 F. D. Arditti, “Risk and the required return on equity,” The Journal of Finance, vol. 22, no. 1, pp. 19–36, 1967. View at: Publisher Site  Google Scholar
 F. D. Arditti, “Another look at mutual fund performance,” The Journal of Financial and Quantitative Analysis, vol. 6, no. 3, pp. 909–912, 1971. View at: Publisher Site  Google Scholar
 R. C. Scott and P. A. Horvath, “On the direction of preference for moments of higher order than the variance,” The Journal of Finance, vol. 35, no. 4, pp. 915–919, 1980. View at: Publisher Site  Google Scholar
 W. Chen and Y. Wang, “A new fuzzy multiobjective higher order moment portfolio selection model for diversified portfolios,” Physica A, vol. 465, pp. 124–140, 2017. View at: Google Scholar
 W. Yue and Y. Wang, “A new fuzzy multiobjective higher order moment portfolio selection model for diversified portfolios,” Physica A: Statistical Mechanics and Its Applications, vol. 465, pp. 124–140, 2017. View at: Publisher Site  Google Scholar
 A. Goel and A. Sharma, “Deviation measure in secondorder stochastic dominance with an application to enhanced indexing,” International Transactions in Operational Research, In press. View at: Publisher Site  Google Scholar
 Y. Chen, Z. Tian, H. Zhang, J. Wang, and D. Zhang, “Strongly constrained discrete hashing,” IEEE Transactions on Image Processing, vol. 29, pp. 3596–3611, 2020. View at: Publisher Site  Google Scholar
 L. Ma, H. Wang, B. Lu, and C. Qi, “Application of strongly constrained space particle swarm optimization to optimal operation of a reservoir system,” Sustainability, vol. 10, no. 12, p. 4445, 2018. View at: Publisher Site  Google Scholar
 B. Zeng, X. Li, L. Gao, Y. Zhang, and H. Dong, “Whale swarm algorithm with the mechanism of identifying and escaping from extreme points for multimodal function optimization,” Neural Computing and Applications, vol. 32, no. 9, pp. 5071–5091, 2020. View at: Publisher Site  Google Scholar
 A. A. Heidari, I. Aljarah, H. Faris, H. Chen, J. Luo, and S. Mirjalili, “An enhanced associative learningbased exploratory whale optimizer for global optimization,” Neural Computing and Applications, vol. 32, no. 9, pp. 5185–5211, 2020. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Q. H. Zhai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.