Abstract

With the development of economy, the requirement of financial planning for individuals or families is emerging. In the era of the Internet, individual investors can conveniently enter the market and purchase financial products. Traditional portfolio management models focus on risky markets such as stock markets. However, risk-averse investors, such as normal families, may concern appropriate long-term financial planning. This paper considers the problem of portfolio management of bank financial products with a long-term planning horizon. By taking into account the final return and the flexibility, a multiobjective model of long-term portfolio is proposed. A multiobjective evolutionary approach is employed for the handling of conflicting objectives. Test instances are generated to illustrate the problem. Experiment results show that the presented algorithm can efficiently find trade-off solutions. Our experimental results also show that crossover probabilities should be separately implemented for long-term portfolio problems with hybrid encoding. Performance comparison of different crossover operators suggest that, for a real-valued encoding part, the simulated binary crossover (SBX) has a better performance than BLX- operator. While for a binary encoding part, a uniform crossover operator might be appropriate for large-scale instances. The proposed multiobjective model in this paper provides risk-averse investors with an appropriate decision support model for the long-term financial planning and management.

1. Introduction

Research studies on portfolio problems have gained much attention since the 1950s. Usually, portfolio return and risk are simultaneously considered as objectives to be optimized. By formulating risk measures from different perspectives, various models have been proposed in the literature, such as Mean-Variance (M-V) model [1], Mean-Semivariance (M-SV) model [2], Mean-Value at Risk (M-VaR) model [3], and Mean-Mean Absolute Deviation (M-MAD) model [4]. These models are widely employed for portfolio management in high-risk markets such as stock markets [57].

In the era of the Internet, the financial market becomes more diversified. Specifically, individual investors can conveniently enter the market. Moreover, various financial products are emerging with the development of the Internet and information technologies [8, 9]. For instance, in China, almost each bank releases financial products via webs or mobile applications. Customers can purchase financial products anywhere and anytime. Compared to stock or fund products, financial products released by banks have lower risks. Many risk-averse investors prefer to choose such types of financial products. For example, most Chinese families tend to save their idle capital in the bank. In such a case, risk might not be the major concern of these risk-averse customers. However, when customers make their long-term financial planning, they still need to decide the type and the amount of products to be purchased, as well as the amount of capital should be reserved. In contrary, investors should manage the balance between the final return and the flexibility in the long-term financial planning.

Compared to traditional portfolio problems, models for a low-risk market with long-term planning horizon are rather scant. In this paper, we address the problem of long-term portfolio management with bank financial products. Based on real data of Chinese banks, this paper proposes a new portfolio model in which the final return and the flexibility are simultaneously optimized. In the proposed model, consideration of multiperiod investment decisions and realistic constraints such as the buy-in threshold increases the complexity of the problem. A multiobjective evolutionary algorithm (MOEA) is presented to solve the problem.

The contribution of the paper is twofold. Firstly, a new multiobjective model of portfolio is proposed. With consideration of realistic requirement of risk-averse investors, flexibility is designed as one of the objectives. The proposed model expands the scope of classic portfolio problems and can provide potential customers in financial markets with effective decision support model. Secondly, an effective algorithm is presented to solve the problem. The performance of different crossovers is investigated. A set of test instances are generated to illustrate the addressed problem and validate the effectiveness of the algorithm.

The remainder of this paper is organized as follows. Section 2 gives a brief review of the related work. A problem description and formulation are presented in Section 3. The solution technique is described in Section 4. Experimental results on generated test instances are reported in Section 5. Finally, summary and some conclusions are provided in Section 6.

Various practical models of portfolio have been reported in the existing literature, such as uncertain portfolio selection problems [1012]. The methodology of uncertain portfolio can be used for other practical problems such as multiobjective multi-item fixed charge solid transportation problem [13]. Practical portfolio problems are usually characterized by nonconvex and discrete objectives as well as constraints. Thus, metaheuristic algorithm such as evolutionary computation (EC) has become a popular solution technique for realistic portfolio problems. Among the first attempts, Arnone et al. [14] suggested the possibility of applying genetic approach to solve portfolio problems and presented some ideas. Maringer and Kellerer [15] developed a hybrid local search algorithm for cardinality constrained portfolios. Streichert et al. [16, 17] compared different crossovers and chromosome representations for constrained portfolio selection problems. Cura [18] presented a particle swarm optimization approach to solve portfolio optimization problems. Pai and Michel [19] studied portfolio problems with different types of constraints. The problems were converted to a single-objective optimization problem with an aggregation parameter. The proposed approach employed k-means cluster analysis to handle the cardinality constraint and used an evolution strategy to serve as an optimizer. Ruiz-Torrubiano and Suárez [20] proposed a memetic algorithm combining a genetic algorithm and quadratic programming for solving portfolio problems with cardinality constraints and transaction costs. Guo et al. [21] addressed a fuzzy multiperiod portfolio selection problem in which assets have different investment horizons. The authors designed a simulation-based genetic algorithm to solve the problem with different objectives, i.e., maximizing mean total return and minimizing variance.

Portfolio problems inherently have the trait of multiobjective optimization. The last two decades have witnessed the growth of taking advantage of MOEAs for solving portfolio problems. Doerner et al. [22] developed a Pareto ant colony optimization approach to solve project portfolio problems with objectives of benefit and remaining resources. Branke et al. [23] proposed an envelope-based multiobjective evolutionary algorithm to solve portfolio selection problems with cardinality constraints and minimum buy-in thresholds. With the rapid emerging of related work, some surveys focused on solving portfolio problems by means of MOEAs [24, 25]. Mishra et al. [26] addressed a realistic portfolio asset selection problem with budget, floor, ceiling, and cardinality constraints. The authors proposed a nondominated sorting multiobjective particle swarm optimization to solve the problem and conducted comparative performance assessment of a set of MOEAs. Kar et al. [27] proposed a fuzzy portfolio selection model for which Sharp ratio (SR) and Value at Risk ratio (VR) were simultaneously optimized by means of a multiobjective cellular genetic algorithm. Liagkouras [28] presented a three-dimensional encoding in an MOEA for portfolio problems. The proposed encoding technique can keep the processing time invariant to the problem size.

Although portfolio optimization is a classic problem in the area of finance and operations research, studies focusing on long-term portfolio planning are rather scant. Multiperiod portfolio can be classified as a type of long-term portfolio. Woodside-Oriakhi et al. [29] considered the problem of rebalancing an existing financial portfolio with transaction costs. Meghwani and Thakur [30] proposed a triobjective portfolio optimization model with risk, return, and transaction costs. The authors adapted classic MOEAs for solving the problem. However, either in single-period or multiperiod portfolio problems, objectives are usually portfolio return and related risks. These models are appropriate for the stock market. To the best of our knowledge, there is no model to consider portfolio flexibility in long-term financial planning in a low-risk market.

3. Problem Statement

3.1. Real Data

In this section, we first present some real data of financial products released by some banks in China. Usually, according to whether a product period exists, bank financial products can be classified as two types: open-ended financial products and close-ended financial products. For the latter type of products, the capital allocated to a product cannot be redeemed until the period is ended. According to the data released by Agricultural Bank of China (ABC) on 07 May 2019 (http://ewealth.abchina.com/fs/filter/), there are 93 financial products available in the market. Most financial products provided by ABS belong to the close-ended type. We select 20 close-ended products and report the detail information in Table 1. We also report similar data from Bank of Communications (BC)(http://www.bankcomm.com/BankCommSite/default.shtml). Products released by BC include close-ended and open-ended types. We select 10 close-ended products and 10 open-ended products and report the information in Table 2.

One can see from the data that products released by banks have three main attributes: minimum buy-in amount, product period, and interest rate. According to the data of ABC and BC, the product period can be modeled by a finite set with unit of month. It can also be seen from the abovementioned tables that interest rates of financial products can be fixed or floating. However, for products with floating interest rates, the floating ranges are rather small. Thus, all financial products can be considered as low risky or medium risky. When a risk-averse investor plans to purchase these bank financial products, she might not concern an immediate higher return of a product with a floating interest rate. What she most concerns is how to appropriately allocate her capital in the long run and obtain a high final return.

3.2. Basic Definitions

Based on the real data provided in Section 3.1, the formulation of the problem studied in this paper is refined. We consider an investment planning in a time horizon with periods. For the sake of simplicity, similar to Gutjahr et al. [31], we assume that each period indicates one month. There are financial products released by commercial banks. Each product has the following basic attributes:(i): buy-in threshold(ii): product period(iii): return rate per time unit

In this paper, we assume that return rate is deterministic and pregiven for each product. The new arrived capital of an investor at each time point is . The total available capital at each time point is . At each time point, the investor needs to allocate her capital to different products. Thus, the problem can be categorized as a multiperiod portfolio problem.

3.3. Example

In multiperiod portfolio problems, an investor reallocates her wealth among available products at the beginning of each period [32]. In the presence of the product period, we have to optimize the investment plan in the long run. We use a simple example to illustrate the problem. We suppose that there are three available products:(i)Product is with a product period of 12 months. The buy-in threshold is 5 units and the return rate is 0.3% per period (one month).(ii)Product is with a product period of 6 months. The buy-in threshold is 10 units and the return rate is 0.4% per period (one month).(iii)Product is with a product period of 1 month. The buy-in threshold is 2 units and the return rate is 0.2% per period (one month).

The investment horizon is 1 year, i.e., 12 periods. We suppose at the beginning of the first and the seventh period, there is available capital with an amount of 8. If there is no period constraint for each product, the optimal investment plan can be shown, as in Figure 1. At each decision point, the investor can choose products with the highest interest rate if the buy-in threshold constraint is satisfied. In this plan, the final return with a compound rate at the end of the planning horizon is .

Then, we consider the situation that the product period constraint is present. At the first decision point, if the investor makes the same decision as in Figure 1, i.e., investing all available budget on the product with the highest interest rate, the optimal investment plan given the first-period decision can be shown, as in Figure 2. At the seventh period, new budget is available. However, at this decision point, the cash invested on product cannot be redeemed because of the product period constraint. Thus, the best choice of the investor at the second decision point is product . With this investment plan, the final return is  = 16.3893.

Clearly the abovementioned plan is not optimal. The optimal plan can be shown, as in Figure 3. The final return of the optimal plan is . For the plan shown in Figure 2, all available capital is allocated to products with the highest return rates. However, this will cause inflexibility in the long run and result in possible loss of final return. Thus, in the long-term financial planning, it is worth considering the flexibility of a plan and making decisions with a global optimal perspective.

3.4. Multiobjective Model

In this section, we formally present the multiobjective model of the addressed long-term portfolio problems.

3.4.1. Decision Variables

At each time point , the investor needs to decide:(i): the capital amount allocated to product ,

3.4.2. Constraints

(i)Buy-in threshold constraint: this type of constraint is similar to a floor-ceiling constraint which imposes lower and/or upper bounds on the capital proportion allocated to each asset. In our problem, this constraint demands that if a product is selected, at least some amount should be bought, i.e., if is not zero.(ii)Total available capital constraint: this constraint ensures that the sum of capital allocated to each product does not exceed the total available amount at each time point . Then, the constraint can be formulated as .

3.4.3. Objectives

For portfolio problems, the “maximizing return-minimizing risk” strategy is generally used as the definition of optimization objective(s) [25]. However, in this paper, we address the problem where the return rate of each product is deterministic. Then, the objective of minimizing risk is not the major task of the problem. Investors might pay much attention on the trade-off between the return and the flexibility of an investment plan. Thus, in our problem, maximizing return and minimizing inflexibility are the two objectives taken into account.

The portfolio return is calculated at the end of planning horizon. Recall that the planning horizon consists of periods. If an investor allocates the cash of to the th product at the beginning of the plan and it remains unchanged during the whole planning horizon, by the end of the period the product with the fixed return rate will be worth . However, investors may rebalance their investment at subsequent decision points according to the newly available capital. Then, the analytical formation of the portfolio return is difficult to present and the return can only be revealed at the end of the whole planning horizon.

Research studies on flexibility have gained much attention in many areas such as economics, organization, and manufacturing. Despite different definitions of flexibility with various backgrounds and perspectives, flexibility is usually used to measure the ability of dealing with future uncertainties. Some researchers measure the flexibility of capital investments by explicitly taking into account the cost of adopting these investments to possible changes in the future environment [33]. For the addressed problem, if an investor allocates all available capital on close-ended products with long periods, she can obtain a high final return at the end point of the planning horizon. However, there is no reserved capital until invested products expire. During this period, the investor is of weak ability to deal with uncertain events, such as configuring other household expenses. We consider such portfolio plans as inflexible. In other words, a higher final return of a portfolio will sacrifice the flexibility of the plan. Thus, the final return and the flexibility are in conflict with each other to certain extend. Then, an investor needs to balance between the final portfolio return and the reserved capital she holds during the planning horizon. Specifically, an investor tends to maximize the final return and minimize the inflexibility of a portfolio.

Before mathematically formulating the inflexibility measure, we first introduce a preference-based parameter which indicates an investor’s acceptable minimum level of holding the reserved capital. Let be the reserved capital at decision point , i.e., capital which is not allocated to any financial products. If , it is considered as a risky situation. Then, an investment plan is considered as flexible if for all . In other words, the optimization task is to minimize a specific measure relating to risky points. Let be the set of risky reserved capital, i.e., . The gap set is represented as . The number of elements in is denoted as . The mean and standard deviation of is, respectively, denoted as and . Then, we formulate the inflexibility measure (IFM) of a portfolio by a modified model as follows:where is a coefficient and is set to be 0.8 in this paper. From the abovementioned equation, one can see that the inflexibility measure corresponds to three parts: average gap to an investor’s acceptable minimum level of holding reserved capital, variation of the gap, and the number of risky points.

It is clear that the inflexibility measure is to be minimized. For the sake of simplicity, we convert the objective of maximization of to minimization of during the solving process. Then, objectives of the addressed problem can be presented as follows:

4. Methodology

The addressed problem is combinatorial and the trade-off surface in the objective space might be nonconvex and disconnected. For such a type of problem, multiobjective evolutionary algorithms are preferable as solution techniques [3437]. In this paper, we develop a solution technique based on a classic multiobjective genetic algorithm named NSGA-II [38], which is widely used to solve other financial problems such as feature selection in credit scoring [39]. In this paper, the chromosome and genetic operators are designed for our problem.

4.1. Chromosome Representation

For portfolio problems, genetic representation should include two types of information: selected product types and capital amount allocated to the selected products. Streichert et al. [16, 17] suggest to use a hybrid binary/real-valued encoding in which the binary part indicates whether an asset is included in the portfolio and the real-valued part is used to represent the weights allocated to individual assets. This type of representation is also employed by other approaches [23]. If there are products, the binary part and the real-valued part can be denoted as and , respectively. Note that, for our problem, there are decision points. Thus, the dimension of the chromosome expands to . Then, we represent the chromosome by hybrid binary/real-valued matrices, as shown in Figure 4.

4.2. Genetic Operators

Crossover and mutation are two main operators for genetic algorithms. For the addressed problem, genetic operators are independently carried out on real-valued and binary matrix chromosomes.

Many types of crossover operators are designed for hybrid encoding of portfolio problems. However, comparative results of the effectiveness of different operators are rather scant. Two popular crossovers for real-valued chromosome are BLX- and simulated binary crossover (SBX). Experimental results in Streichert et al. [16, 17] suggest that BLX-0.5 crossover performed best amongst three investigated recombination operators. In a recent work, SBX is used for crossing real-valued vectors [28]. In our experiments, we will analyze performances of these two different crossovers. The procedure of crossover for real-valued matrix is presented in Algorithm 1.

(1)for t = 1 to T do
(2)for n = 1 to N do
(3)  Generate from [0, 1]
(4)  if then
(5)   Apply BLX-/SBX crossover operator on element
(6)   end if
(7)  end for
(8)end for

In classical portfolio problems, single-point crossover is usually employed to crossing binary vector [28, 40]. In our problem, crossover is carried out at each decision point, i.e., on each row of the binary matrix. Then, the procedure of single-point crossover is presented as in Algorithm 2.

(1)for t = 1 to T do
(2) Generate from [0, 1]
(3)if then
(4)  Randomly generate a crossover point .
(5)  Apply single-point crossover operator on the th row with .
(6)end if
(7)end for

Note that, in classic single-point crossover, two offspring are generated by swapping bits in the right part of the crossover point between two parents. For binary strings, swapping the right parts or the left parts of the parents does not affect the generated offspring. However, in our problem, the chromosome is a binary matrix. Swapping the fixed parts of parents might generate biased children. Thus, a uniform crossover might be more appropriate for our problem. With a uniform crossover, offspring is obtained by randomly inheriting chromosome parts from two parents. In our experiments, we will compare the performances of single-point crossover and uniform crossover with a single crossover point, as well as two-point crossover and uniform crossover with two crossover points.

For the real-valued matrix, a widely used operator named polynomial mutation is employed for mutating the solution. The mutation operator for binary matrix is shown in Algorithm 3.

(1)for t = 1 to T do
(2) Generate from [0, 1]
(3)if then
(4)  Randomly generate a position index from [1, N].
(5)  if then
(6)   
(7)  else
(8)   
(9)  end if
(10)end if
(11)end for

5. Experimental Results

5.1. Test Instances Generation

Compared to classic portfolios, the addressed problem has two different characteristics: (1) product might have different periods and (2) budget is available at each decision point. Thus, existing benchmark of portfolio cannot be used in our problem. Then, we randomly generated test instances to illustrate the problem and investigate the performance of the algorithm. In real cases, return rates usually depend on two parameters: product period and buy-in threshold . From the real data reported in Section 3.1, for most cases, products with higher values of and usually have higher values of . Then, in the test generation, we explicitly model the relation of parameters of , , and as follows:where is baseline return rate and is the upper bound of the buy-in threshold for all products. We are aware of that this formulation might be different with real cases. However, in test instance generation, this will not affect the analysis process and related conclusive results. In future research, the relation between different parameters will be refined with more real data.

We generated two sets of instances, denoted as and . The includes two instances, denoted as and , with four products and a planning horizon of 24 periods. In both instances, in , parameters for the 4 products are the same, as shown in Table 3. The difference of the two instances lies in the available capital set. For , the set is , while the set is for . The frequency of capital arrival in is lower than that in . Intuitively, has a more complex structure compared to .

The consists of six instances. Table 4 shows parameters in . The product number varies from 10 to 35. The planning period is 60 units. The other parameters are randomly drawn from intervals given in the table. For example, the parameter of instance 1 is randomly generated between the lower bound 1 and the upper bound 50 for each time unit.

5.2. Parameter Settings

For both instance sets, the population size was set as 200. The algorithm was terminated after 100 generations for . Since has a larger size, the generation number was set as 200. As suggested in other studies [16, 17], crossover rates and were set as 0.9 and the mutation rate was 0.1. In our experiments, sensitivity of the performance of the algorithm on the crossover rates was analyzed on with a constant mutation rate 0.1. In the calculation of hypervolume and IGD, objective values were normalized. For , maximum and minimum values for horizontal axis were −200 and −800 and the values for vertical axis were 50 and 0. These values for were set as −1000 and −3500 for horizontal axis and 60 and 0 for vertical axis, respectively. These values were determined with analysis of results obtained by a single run of the algorithm on and . All algorithms were repeated for 30 independent runs.

5.3. Four-Product Examples
5.3.1. Results of

We first plot the Pareto front of with different , as shown in Figure 5. We plot portfolio return (PR) on the horizon axis and inflexibility measure (IFM) on the vertical axis. For each , the two endpoints of the Pareto front correspond to the portfolio plans with the maximum and the minimum averaged reserve capital. We list the two endpoints as well as the middle point of each Pareto front in Table 5.

From Figure 5 and Table 5, one can see that the values of the maximum portfolio return are comparable despite different levels of . However, the value of IFM varies from 2.47 to 12.45 . This is because the calculation of IFM depends on the predefined value of while the return is only determined by investing structures of portfolios. Figure 6 shows investing structures of endpoints with the maximum return under different levels of for . It can be seen that the structures of different endpoints with the maximum return are very similar.

Structures of endpoints with the minimum return under different levels of for are shown in Figure 7. The main difference of portfolio structures in Figure 7 lies in the level of reserved capital at each decision point. Endpoints with the minimum return correspond to the portfolios in which the average reserved capital approaches to the pregiven level of . If is set at a higher level, less capital is available for investing on products. Thus, the return of the corresponding portfolio will be lower compared to other portfolios with lower levels of . This can be seen in Figure 5 and Table 5, and returns of the endpoint with is 224.17 while the return is only 203.35 for . From the perspective of decision making, an investor with a high level of can be considered as a risk-averse decision maker. While the investor with a low level of is a risk-taking decision maker. It is reasonable that, for a risk-taking investor she might not accept portfolios with low return. While for a risk-averse investor, portfolios with low return but high level of flexibility will be acceptable. As shown in Figure 5, given limited computational resource, obtained nondominated solutions focus on the area of high returns for .

5.3.2. Results of

We then plot the obtained approximate Pareto front of with different in Figure 8. One can see from Figure 8 that the shapes of the obtained approximate efficient fronts are closer to be linear, compared to the efficient fronts shown in Figure 5. The difference arises from the structure of the problem instances. For , a lower level of capital arrival frequency decreases the complexity of the problem and makes the investment policy simpler and clearer.

For , investing structures of obtained endpoints with maximum return under different levels of are shown, as in Figure 9. We recall that, for , there is sufficient available capital at each decision point. Thus, an investor can allocate all capital on one or more products with high return rates. Theoretically, endpoints with the maximum return correspond to portfolios in which the reserved capital at each point remains the level of 0. This can be seen from Figure 9 that the level of reserved capital is rather flat and approaches to 0.

5.3.3. Comparison of BLX and SBX Crossovers

In this section, we compare BLX- and SBX operators on . As suggested in Streichert et al. [16, 17], the value of was set as 0.5 for BLX- crossover. The value of was set as 10. Obtained nondominated solutions for and with different crossover operators are plotted in Figure 10. It is clear that SBX performed better than BLX-0.5 for both instances. Note that, for , the gap between nondominated solutions obtained with BLX-0.5 and SBX is larger than that of . As aforementioned, is less complex than and thus is easy to be solved. Since SBX outperforms BLX-0.5, we employed SBX as the crossover operator in the following experiments.

5.4. Results on

In this section, we first investigated performances with different crossover probabilities. Then, performances with different crossover mechanism for binary matrix were compared. In the performance comparison, two most often used indicators named hypervolume [41] and Inverted Generational Distance (IGD) [42] were employed as performance measures. A larger value of hypervolume indicates a better performance and it is in contrast for IGD value. In the calculation of IGD, an approximating Pareto optimal set was obtained by merging all nondominating solutions in each run. Note that, for , the value of was 30 for all instances.

5.4.1. Sensitivity Analysis of Crossover Probability

In this section, we provide sensitivity analysis of the crossover probability. As aforementioned, SBX crossover is employed for real-valued matrix. We first chose the single-point crossover for binary matrix in the analysis. The chromosome consists of real-valued and binary parts. Thus, it is worth investigating the sensitivity of the crossover probability for different parts. In previous research studies such as Liagkouras [28], the crossover probability is set identical for different chromosome parts. Firstly, we analyze the sensitivity of the crossover probability for the real-valued chromosome part. Figures 11 and 12 show the results of hypervolume and IGD with different levels of the crossover probability for the real-valued matrix. In the experiments, the crossover probability for the binary matrix and the mutation probabilities were constant and set as , , and . One can see from the two figures that the algorithm performed worst with . What should be noted is that a high level of the real-valued matrix crossover probability, , as suggested in other studies, could not achieve the best performance. While with an intermediate level of or , the algorithm performed better. If we look at the average values of hypervolume and IGD in 30 runs, compared with , the performance with shows some advantages for the first 4 instances. However, for and , the algorithm with performed better than . This difference may be derived from the scale of the instance. For large scale instances, a relatively lower level of might be more effective. Then, in the following experiments, we fixed for to and for and .We then analyzed the sensitivity of with constant . Box charts of hypervolume and IGD values are reported in Figures 13 and 14, respectively. One can see that, for most instances, the algorithm with performed best. Only for , the algorithm performed best with and the performance was followed by the algorithm with . Then, we fixed in the following experiments. It is noteworthy that effective crossover probabilities for real-valued and binary parts of chromosome are not identical. However, in the literature, the crossover probabilities for both parts were usually set to be the same. Based on the abovementioned results, we suggest to independently investigate the crossover mechanism for portfolio problems with a hybrid encoding.

5.4.2. Performance Comparison of Crossover Operators

Finally, we investigated the effectiveness of four types of crossovers for the binary matrix. Metric values of hypervolume and IGD for single-point crossover, uniform single-point crossover, two-point crossover, and uniform two-point crossover are presented in Figures 15 and 16, respectively. It can be seen from the figures that it is difficult to determine which crossover performed the best for the test instances. In contrary, it can be considered that these four crossovers have no statistically significant differences. However, decision makers concern the obtained nondominated solutions to a great extent. Thus, we plotted the obtained nondominated solutions in the objectives for all six test instances in Figure 17. One can see from the figure that, for and , the algorithm with the single-point crossover obtained a better approximated Pareto front. While uniform crossovers performed better for , , and . The four crossovers had no remarkable difference for . From the abovementioned results, we may cautiously suggest the usage of uniform crossover for large-scale problems. It can be seen that the mechanism of the crossover for binary matrix affects the quality of obtained nondominated solutions to a large degree. Thus, in future research studies, it is worthwhile to investigate the designing of more effective crossover mechanisms for the addressed problem.

6. Conclusions

In this paper, we address the problem of portfolio management of bank financial products with a long-term planning horizon. Instead of considering portfolio risk, we concentrate on the flexibility of a portfolio. A multiobjective optimization model with objectives of portfolio return and flexibility is presented. The proposed model is also characterized by multiperiod investment where capital is available step by step and the decision maker needs to decide her investment behavior at each decision point. Due to the nonconvex and discrete characteristics of the problem, a multiobjective evolutionary algorithm is presented as the solution technique. As suggested in the literature, a hybrid binary/real-valued chromosome is employed for encoding. In the experimental analysis, two sets of test instances are randomly generated for illustrating the addressed problem and testing the proposed approach. Experimental results of parameter sensitivity suggest that crossovers on different parts of the chromosome should be separately considered. Performance comparison of different crossovers show that the effective crossover mechanism for the problem is still a quite open issue. Additionally, phenotype analysis is presented based on two small size instances.

The presented model and solution technique proposed in this paper provides risk-averse investors with an appropriate decision support tool. When the investors make a long-term financial planning, they can trade off the final return and the flexibility of a portfolio. However, we are aware of that some issues need to be considered in the model, such as unexpected changes of financial products, uncertainty of return rates, and dynamic or incremental decision in the long run. We will focus on these topics in future research studies.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant nos. 71501181, 71871185, 71910107002, and U1811462.