# Modeling Portfolio Optimization Problem by Probability-Credibility Equilibrium Risk Criterion

Authors: Ye Wang; Yanju Chen; YanKui Liu
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

This paper studies the portfolio selection problem in hybrid uncertain decision systems. Firstly the return rates are characterized by random fuzzy variables. The objective is to maximize the total expected return rate. For a random fuzzy variable, this paper defines a new equilibrium risk value (ERV) with credibility level beta and probability level alpha. As a result, our portfolio problem is built as a new random fuzzy expected value (EV) model subject to ERV constraint, which is referred to as EV-ERV model. Under mild assumptions, the proposed EV-ERV model is a convex programming problem. Furthermore, when the possibility distributions are triangular, trapezoidal, and normal, the EV-ERV model can be transformed into its equivalent deterministic convex programming models, which can be solved by general purpose optimization software. To demonstrate the effectiveness of the proposed equilibrium optimization method, some numerical experiments are conducted. The computational results and comparison study demonstrate that the developed equilibrium optimization method is effective to model portfolio selection optimization problem with twofold uncertain return rates.


---

## Body

## 1. Introduction

Based on mean-variance criterion, Markowitz [1] first established portfolio theory. In Markowitz’s mean-variance model, the returns of individual securities are taken as random variables and the expected value and variance of the random return are taken as the investment return and risk, respectively. The mean-variance theory has been accepted widely as a tool for building portfolio optimization problems, and many interesting works have further developed the mean-variance method. Considering that the situation of the expected return is lower than the average value, Markowitz [2] replaced variance with semivariance as the risk measure to describe the risk. Semivariance is an improvement of variance, which allows us to obtain the higher expected return rate under the same level of risk. Konno and Yamazaki [3] used mean absolute deviation model in Tokyo stock market. Simaan [4] showed that estimation error is large in both mean-variance model and mean absolute deviation model and severe in small samples and for investors with high risk tolerance. In order to describe the normal fluctuations in the market, the concept of value-at-risk (VaR) was taken as the risk measure. VaR is the maximum possible loss of a financial asset or portfolio of securities. In recent years, VaR gradually progresses as a positive risk management tool and plays an increasingly important role in the portfolio optimization problems. Jorion [5] studied VaR as a risk measure and applied mean-VaR model in finance industry. Mean-variance-skewness models [6, 7] were also proposed to study portfolio selection problem. In addition, there are other measures used to characterize the risk [8]. When the restriction that an asset cannot be traded for intervals of uncertain duration results in risk, a model of optimal allocation to liquid and illiquid assets was presented by Ang et al. [9]. Shen et al. [10] discussed a mean-variance portfolio selection problem under a constant elasticity of variance model and expressed explicit expressions of the optimal portfolio strategy, the value function, and the efficient frontier of the mean-variance problem.The conventional portfolio methods assume that the security returns are random variables. Probability distributions of random variables are usually derived from historical data. However, in the real investment environment, the security returns often present vagueness and ambiguity. Since the seminal works of Zadeh [11, 12], many achievements have been made based on fuzzy theory. Watada [13] discussed portfolio selection by using fuzzy decision theory. Tanaka and Guo [14] used possibility distributions to model the uncertain returns. Inuiguchi and Ramík [15] exemplified the advantages and disadvantages of fuzzy mathematical programming approaches in the setting of an optimal portfolio selection problem. When the returns of assets are taken as trapezoidal fuzzy numbers, a possibilistic approach for selecting portfolios with the highest utility value was introduced by Carlsson et al. [16]. Considering transaction costs, Fang et al. [17] proposed a portfolio rebalancing model based on fuzzy decision theory. Mean-variance models were constructed for fuzzy portfolio selection problems by Chen et al. [18] and Zhang et al. [19]. Qin et al. [20] discussed the Kapur cross-entropy minimization model for portfolio selection problem under fuzzy environment. Dastkhan et al. [21] studied a linguistic-based portfolio selection model by weighted max-min operator. Wu and Liu [22] developed a robust method to describe fuzzy returns by employing parametric possibility distributions. Chen et al. [23] applied the absolute semideviation as a new risk measure and developed three classes of fuzzy portfolio optimization models. Kamdem et al. [24] introduced notions of moments and semimoments for credibilistic portfolio problems to measure the leptokurtocity of credibilistic portfolio return and presented a mean-variance-skewness-semikurtosis model. Mehlawat and Gupta [25] addressed fuzzy portfolio selection problem from a perspective of chance constrained multiobjective programming. Deng and Li [26] proposed a biobjective nonlinear portfolio selection model, which aimed to maximize the future expected return and minimize the future expected risk. They also proposed a gradually tolerant constraint method to solve this model. Li et al. [27] developed a fuzzy portfolio selection model with background risk. Vercher and Bermúdez [28] introduced a cardinality constrained multiobjective optimization problem for generating efficient portfolios within a fuzzy mean absolute deviation framework. Considering transaction cost, Chen and Wang [29] proposed a two-stage fuzzy model for portfolio selection problem. For fuzzy portfolio selection problems, there are many other different optimization methods.In modern financial markets, there is not only randomness but also fuzziness that affects the total decision-making process. Based on the above works, this paper models the portfolio selection problem where randomness and fuzziness are considered simultaneously. Huang [30] studied portfolio selection in a random fuzzy environment and proposed two portfolio selection models with random fuzzy returns. We adopt the equilibrium optimization method [31] to model portfolio selection problems. In order to integrate uncertainties resulting from both subjective consciousness and objective factors, random fuzzy variables [32] are used to characterize the return rates. In this situation, the security return rates are assumed to be random variables with fuzzy information. This paper employs the expected value (EV) [33] to represent the investment return and introduces a new index called ERV for random fuzzy return rate to measure the investment risk. Then a random fuzzy EV-ERV model is proposed for portfolio selection problem. Since optimization models with twofold uncertainty are difficult to solve, the heuristic algorithms are often designed to solve the proposed optimization models. An advantage of the proposed equilibrium optimization framework is to turn the original equilibrium portfolio optimization problem to a computationally tractable problem. Firstly, the proposed model is translated into an equivalent credibilistic optimization model. Then the credibilistic model is reduced to its equivalent deterministic convex programming models under mild assumptions on fuzzy parameters.Compared with the existing literature, the main contributions of this paper consist of the following three aspects. Firstly, this paper defines a new ERV with credibility level beta and probability level alpha to measure the investment risk. The adopted risk measure quantifies the uncertainties of randomness and fuzziness simultaneously. This method shows the qualitative and quantitative analysis about the uncertainty of return rate. The proposed EV-ERV model is a useful optimization method for a practical investor. Secondly, the global optimal solution of the proposed EV-ERV model is obtained. In the case that the randomness of uncertain return rates follows normal distributions with deterministic covariance matrix, and the fuzziness is characterized by trapezoidal fuzzy variables, triangular fuzzy variables, or normal fuzzy variables, the proposed EV-ERV model is transformed into its deterministic convex programming models, which can be solved by general purpose optimization software. Thirdly, when the randomness of uncertain return rate vector follows multivariate normal distribution, the covariance matrix can reflect the interactions and correlation degrees among securities.This paper is organized as follows. In Section2, after introducing a new index called ERV for random fuzzy variable, the random fuzzy EV-ERV model is proposed for portfolio selection problems. In Section 3, under mild assumptions on possibility and probability distributions, the EV-ERV equilibrium optimization model is proved to be a convex programming. Section 4 deals with the equivalent deterministic convex programming models under three situations. In Section 5, we perform some comparison study via numerical experiments, give the sensitivity analysis about model parameters, and show the relationship between ERV and EV. The last section gives the conclusions of the paper.


## 2. Equilibrium Portfolio Selection Problem

### 2.1. Random Fuzzy Variable

LetΓ be an abstract space of generic elements γ and P(Γ) the power set of Γ. If Pos is a possibility measure on P(Γ), then the credibility measure Cr is defined as follows [34]: (1)CrA=121+PosA-PosAc,A∈PΓ,where Ac=Γ∖A. The triplet (Γ,P(Γ),Cr) is called a credibility space.Definition 1 (see [32]).
Let(Γ,P(Γ),Cr) be a credibility space and (Ω,Σ,Pr) a probability space. If ξ is a map defined on Γ×Ω such that, for each fixed γ∈Γ, ξ(γ,ω) as a function of ω is a random variable defined on the probability space (Ω,Σ,Pr), then one calls ξ a random fuzzy variable.Ifξ1,ξ2,…,ξn are random fuzzy variables defined on the space Γ, then ξ=(ξ1,ξ2,…,ξn) is a random fuzzy vector defined on the space Γ.A multivariate normal distributionN(μ,Σ) with uncertain mean μ and covariance Σ is an example of random fuzzy vector. In this case, the components of μ and Σ are represented as fuzzy variables with known possibility distributions.Definition 2 (see [33]).
Letξ be a random fuzzy variable. The expected value E[ξ] of ξ is defined as the following fuzzy integral: (2)Eξ=∫0+∞Crγ∈Γ∣Eξγ≥rdr-∫-∞0Crγ∈Γ∣Eξγ≤rdr,where E[ξγ] is the expected value of random variable ξγ for any given γ∈Γ.Definition 3 (see [31]).
Letξ=(ξ1,ξ2,…,ξn) be a random fuzzy vector, and let fj:Rn→R be continuous function for j=1,2,…,m. Then the equilibrium chance Ch of a random fuzzy event is defined as (3)Chfjξ≤0,j=1,2,…,m=supα,β∈0,12α∧β∣Crγ∈Γ∣Prfjξγ≤0,j=1,2,…,m≥α≥β.Equation (3) can be rewritten as the following equivalent form: (4)Chfjξ≤0,j=1,2,…,m=supα∈0,1α∧Crγ∈Γ∣Prfjξγ≤0,j=1,2,…,m≥α.According to the property of equilibrium chance, in the case ofα=β, one has (5)Crγ∣Prξγ≥z≥α≥α⟺Chξ≥z≥α.We next introduce the equilibrium risk value for a random fuzzy variable.Definition 4.
Letξ be a random fuzzy variable. The equilibrium risk value of ξ is defined as (6)VaRα,βξ=sup⁡z∣Crγ∣Pr⁡ω∣ξγω≥z≥α≥β,where α,β are prescribed confidence levels and take their values in the unit interval [0,1].In (6), the meanings of parameters α and β are different: the parameter α represents probability level, while the parameter β represents credibility level.If the random fuzzy variableξ reduces to a random variable, then the equilibrium risk value defined by (6) reduces to (7)ξsupα=sup⁡z∣Pr⁡ω∣ξω≥z≥α,which is the α-optimistic value of random variable ξ. On the other hand, if the random fuzzy variable ξ reduces to a fuzzy variable, then the equilibrium risk value defined by (6) reduces to (8)ξsupβ=sup⁡z∣Crγ∣ξγ≥z≥β,which is the β-optimistic value of fuzzy variable ξ. Therefore, the equilibrium risk value of random fuzzy variable is a natural extension of both the optimistic value of random variable and the optimistic value of fuzzy variable.


### 2.2. Formulation of Equilibrium Optimization Model

With the rapid development of the economic and society, more and more investors realize the importance of portfolio optimization problem under uncertainty. They hope to use the limited funds to get the maximum benefit and, at the same time, to minimize the investment risk as much as possible. However, under the market economy, there are various uncertain factors to affect investment market. Thus, the uncertainties should be taken into account during the modeling process.A rational investor should pursuit the maximum profit with the minimum risk. However, in real investment process, increasing the return, the investors have to tolerate greater risk; the lower risk corresponds to the less return. That is, investors need to make a trade-off between return and risk. Assume that the investor has a collection ofn optional risky assets. The return rate of investment refers to the result of the net income divided by the initial funds. We use random fuzzy variable ηi to denote the return rate of risky asset i, i=1,2,…,n. Then η=(η1,η2,…,ηn) is the return rate vector. When the return rates of assets are characterized by joint normal distribution N(μ,Σ), the fuzzy parameters μ and Σ stand for the average return rate vector and the correlations among the return rates of assets, respectively. Since the parameters μ and Σ often cannot be determined exactly due to the lack of the historical data, it is appropriate to represent the components of μ and Σ as fuzzy variables with known possibility distributions.Letxi be the investment ratio for risky asset i and let x=(x1,x2,…,xn) be the investment ratio vector such that ∑i=1nxi=1, where xi≥0, i=1,2,…,n. The vector x=(x1,x2,…,xn) is called a portfolio. So ηTx is the return rate during the holding period. We take E[ηTx] as objective function. Then, based on Definition 2, the objective function is computed by(9)EηTx=∫0+∞Crγ∈Γ∣EηγTx≥rdr-∫-∞0Crγ∈Γ∣EηγTx≤rdr.We take ERV as the risk index to measure the investment risk. Based on Definition4, the risk of portfolio x is (10)VaRα,βηTx=sup⁡z∣Crγ∣Pr⁡ω∣ηγTωx≥z≥α≥β,where α,β∈[0,1] are prescribed confidence levels.Based on the notations above, if the investor wants to maximize the expected return rate, the equilibrium portfolio optimization problem in hybrid uncertain decision systems is built as the following EV-ERV model:(11)maxEηTxs.t.VaRα,βηTx≥κ∑i=1nxi=1xi≥0,i=1,2,…,n,where parameter κ is the prescribed level of ERV.Introducing an additional variablez, the original EV-ERV model (11) can be equivalently represented as (12)maxEηTxs.t.Crγ∣Pr⁡ηγTx≥z≥α≥βz≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.


## 2.1. Random Fuzzy Variable

LetΓ be an abstract space of generic elements γ and P(Γ) the power set of Γ. If Pos is a possibility measure on P(Γ), then the credibility measure Cr is defined as follows [34]: (1)CrA=121+PosA-PosAc,A∈PΓ,where Ac=Γ∖A. The triplet (Γ,P(Γ),Cr) is called a credibility space.Definition 1 (see [32]).
Let(Γ,P(Γ),Cr) be a credibility space and (Ω,Σ,Pr) a probability space. If ξ is a map defined on Γ×Ω such that, for each fixed γ∈Γ, ξ(γ,ω) as a function of ω is a random variable defined on the probability space (Ω,Σ,Pr), then one calls ξ a random fuzzy variable.Ifξ1,ξ2,…,ξn are random fuzzy variables defined on the space Γ, then ξ=(ξ1,ξ2,…,ξn) is a random fuzzy vector defined on the space Γ.A multivariate normal distributionN(μ,Σ) with uncertain mean μ and covariance Σ is an example of random fuzzy vector. In this case, the components of μ and Σ are represented as fuzzy variables with known possibility distributions.Definition 2 (see [33]).
Letξ be a random fuzzy variable. The expected value E[ξ] of ξ is defined as the following fuzzy integral: (2)Eξ=∫0+∞Crγ∈Γ∣Eξγ≥rdr-∫-∞0Crγ∈Γ∣Eξγ≤rdr,where E[ξγ] is the expected value of random variable ξγ for any given γ∈Γ.Definition 3 (see [31]).
Letξ=(ξ1,ξ2,…,ξn) be a random fuzzy vector, and let fj:Rn→R be continuous function for j=1,2,…,m. Then the equilibrium chance Ch of a random fuzzy event is defined as (3)Chfjξ≤0,j=1,2,…,m=supα,β∈0,12α∧β∣Crγ∈Γ∣Prfjξγ≤0,j=1,2,…,m≥α≥β.Equation (3) can be rewritten as the following equivalent form: (4)Chfjξ≤0,j=1,2,…,m=supα∈0,1α∧Crγ∈Γ∣Prfjξγ≤0,j=1,2,…,m≥α.According to the property of equilibrium chance, in the case ofα=β, one has (5)Crγ∣Prξγ≥z≥α≥α⟺Chξ≥z≥α.We next introduce the equilibrium risk value for a random fuzzy variable.Definition 4.
Letξ be a random fuzzy variable. The equilibrium risk value of ξ is defined as (6)VaRα,βξ=sup⁡z∣Crγ∣Pr⁡ω∣ξγω≥z≥α≥β,where α,β are prescribed confidence levels and take their values in the unit interval [0,1].In (6), the meanings of parameters α and β are different: the parameter α represents probability level, while the parameter β represents credibility level.If the random fuzzy variableξ reduces to a random variable, then the equilibrium risk value defined by (6) reduces to (7)ξsupα=sup⁡z∣Pr⁡ω∣ξω≥z≥α,which is the α-optimistic value of random variable ξ. On the other hand, if the random fuzzy variable ξ reduces to a fuzzy variable, then the equilibrium risk value defined by (6) reduces to (8)ξsupβ=sup⁡z∣Crγ∣ξγ≥z≥β,which is the β-optimistic value of fuzzy variable ξ. Therefore, the equilibrium risk value of random fuzzy variable is a natural extension of both the optimistic value of random variable and the optimistic value of fuzzy variable.


## 2.2. Formulation of Equilibrium Optimization Model

With the rapid development of the economic and society, more and more investors realize the importance of portfolio optimization problem under uncertainty. They hope to use the limited funds to get the maximum benefit and, at the same time, to minimize the investment risk as much as possible. However, under the market economy, there are various uncertain factors to affect investment market. Thus, the uncertainties should be taken into account during the modeling process.A rational investor should pursuit the maximum profit with the minimum risk. However, in real investment process, increasing the return, the investors have to tolerate greater risk; the lower risk corresponds to the less return. That is, investors need to make a trade-off between return and risk. Assume that the investor has a collection ofn optional risky assets. The return rate of investment refers to the result of the net income divided by the initial funds. We use random fuzzy variable ηi to denote the return rate of risky asset i, i=1,2,…,n. Then η=(η1,η2,…,ηn) is the return rate vector. When the return rates of assets are characterized by joint normal distribution N(μ,Σ), the fuzzy parameters μ and Σ stand for the average return rate vector and the correlations among the return rates of assets, respectively. Since the parameters μ and Σ often cannot be determined exactly due to the lack of the historical data, it is appropriate to represent the components of μ and Σ as fuzzy variables with known possibility distributions.Letxi be the investment ratio for risky asset i and let x=(x1,x2,…,xn) be the investment ratio vector such that ∑i=1nxi=1, where xi≥0, i=1,2,…,n. The vector x=(x1,x2,…,xn) is called a portfolio. So ηTx is the return rate during the holding period. We take E[ηTx] as objective function. Then, based on Definition 2, the objective function is computed by(9)EηTx=∫0+∞Crγ∈Γ∣EηγTx≥rdr-∫-∞0Crγ∈Γ∣EηγTx≤rdr.We take ERV as the risk index to measure the investment risk. Based on Definition4, the risk of portfolio x is (10)VaRα,βηTx=sup⁡z∣Crγ∣Pr⁡ω∣ηγTωx≥z≥α≥β,where α,β∈[0,1] are prescribed confidence levels.Based on the notations above, if the investor wants to maximize the expected return rate, the equilibrium portfolio optimization problem in hybrid uncertain decision systems is built as the following EV-ERV model:(11)maxEηTxs.t.VaRα,βηTx≥κ∑i=1nxi=1xi≥0,i=1,2,…,n,where parameter κ is the prescribed level of ERV.Introducing an additional variablez, the original EV-ERV model (11) can be equivalently represented as (12)maxEηTxs.t.Crγ∣Pr⁡ηγTx≥z≥α≥βz≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.


## 3. Analysis of the Equilibrium Portfolio Model

In this section, we discuss the properties of the objective function and constrains in equilibrium portfolio optimization model (12), which facilities finding the appropriate algorithm to solve the model.


### 3.1. Computing Expected Return Rate

Since the return rates are denoted by the random fuzzy vectorη, for any given γ∈Γ, ηγ is the random vector (η1,γ,η2,γ,…,ηn,γ). Thus, the expected return rate is calculated by the following formula: (13)EηTx=EγEωηγTx=Eγ∑i=1nxiEωηi,γ.The first equality in (13) holds true due to the definition of expected value operator of random fuzzy variable [33], and the second equality in (13) is based on the linearity of expected value operator of random variable. Note that ∑i=1nxiE[ηi,γ] is the linear combination of fuzzy variables E[ηi,γ], i=1,2,…,n. To further compute Eγ[∑i=1nxiEω[ηi,γ]], the independence condition or comonotonicity condition for fuzzy variables Eω[ηi,γ], i=1,2,…,n is required. We next discuss the two situations, respectively.Case 1 (independence condition [35]).
According to [35], the expected value operator of fuzzy variable has the property of independence linearity. Based on this property, if E[η1,γ],E[η2,γ],…,E[ηn,γ] are mutual independent, one has(14)Eγ∑i=1nxiEωηi,γ=∑i=1nxiEγEωηi,γ=∑i=1nxiEξi.Case 2 (comonotonicity condition [34]).
According to [34, Theorem 1], the expected value operator of fuzzy variable has the property of comonotonic linearity. Based on this property, if E[η1,γ],E[η2,γ],…,E[ηn,γ] are comonotonic, one has(15)Eγ∑i=1nxiEωηi,γ=∑i=1nxiEγEωηi,γ=∑i=1nxiEξi.In the above two cases, the expected value operator of fuzzy variable has linear property,E[ηTx]=∑i=1nxiE[ηi]. In general case, we may compute Eγ[∑i=1nxiEω[ηi,γ] by using fuzzy simulation or approximation method [33, 36] based on the possibility distributions of fuzzy variables Eω[ηi,γ], i=1,2,…,n.


### 3.2. Processing Equilibrium Risk Value

In this subsection, we will handle the following probability constraint in model (12):(16)Pr⁡ηγTx≥z≥α,where the fuzzy parameter γ∈Γ is given in advance.It is assumed that random vectorηγ follows a multivariate normal distribution N(μγ,Σγ). The covariance matrix Σγ is a symmetric positive semidefinite matrix. So there exists a lower triangular and nonsingular matrix Dγ, μγ∈Rn, and a random vector ξ with components being independent standard normal random variables such that ηγ=Dγξ+μγ. Thus, one has E[ηγ]=μγ and Σγ=DγDγT. If Σγ is positive definite, then the multivariate normal distribution is nondegenerate. This case appears if and only if Dγ has full row rank. Otherwise the distribution is degenerate or singular. The multivariate normal distribution is uniquely determined by the expected value vector μγ and the covariance matrix Σγ.If the multivariate normal distribution is nondegenerate, then it is absolutely continuous. Thus, one hasηγTx~N(μγTx,xTDγDγTx).LetG(x)=Pr⁡{ηγTx≥z}. If xTDγDγTx>0, by standardization, G(x) has the following equivalent representation: (17)Pr⁡ηγTx≥z=Pr⁡ηγTx-μγTxxTDγDγTx≥z-μγTxxTDγDγTx=ΦμγTx-zxTDγDγTx.Taking all the cases ofG(x) into account, the following analytical expression is obtained: (18)Gx=1,ifxTDγDγTx=0,μγTx≥z,0,ifxTDγDγTx=0,μγTx<z,ΦμγTx-zxTDγDγTx,ifxTDγDγTx≠0.For any givenγ∈Γ, if xTDγDγTx≠0, then one has (19)Gx≥α⟺ΦμγTx-zxTDγDγTx≥α⟺μγTx-zxTDγDγTx≥Φ-1α⟺μγTx≥Φ-1αxTDγDγTx+z.From the above analysis, ifE[η1,γ],E[η2,γ],…,E[ηn,γ] are comonotonic or mutual independent, then model (12) can be turned into the following credibilistic programming problem: (20)max∑i=1nxiEμis.t.Crγ∣μγTx≥Φ-1αxTDγDγTx+z≥βz≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.


### 3.3. The Convexity of Feasible Region

Section3.2 has dealt with the probability risk. In this subsection, we will address the equivalent representation of credibilistic constraint and discuss the convexity of feasible region.LetC={x∣Cr{μTx≥Φ-1(α)xTDDTx+z}≥β}. In real-life portfolio selection problems, small confidence levels are meaningless. Thus, in this paper, we consider the case where α,β are in the interval [0.5,1].Theorem 5.
Assume thatη~N(μ,Σ), where Σ is deterministic matrix, μi, i=1,2,…,n, are mutual independent fuzzy variables, and μTx is continuous. In the case of α,β≥0.5, the following results hold.(i)
Cr{μTx≥Φ-1(α)xTDDTx+z}≥β⇔Φ-1(α)xTDDTx+z-(μTx)sup(β)≤0.(ii)
C={x∣Cr{μTx≥Φ-1(α)xTDDTx+z}≥β} is a convex set.Proof.
We first prove the necessity of assertion (i). By the definition of the optimistic value of fuzzy variableμTx, one has (21)μTxsupβ≥Φ-1αxTDDTx+z,which implies (22)Φ-1αxTDDTx+z-μTxsupβ≤0.The necessity of assertion (i) holds true.
We next prove the sufficiency of assertion (i). According to the definition of optimistic value, we know the credibilistic constraintCr{μTx≥(μTx)sup(β)}≥β holds. Since (23)μTxsupβ≥Φ-1αxTDDTx+zis equivalent to (24)Φ-1αxTDDTx+z-μTxsupβ≤0,we have (25)CrμTx≥Φ-1αxTDDTx+z≥βdue to the continuity of μTx. The sufficiency is proved.
In what follows, we prove assertion (ii).
According to the property of optimistic value, the equality(μTx)sup(β)=∑i=1nxiμi,sup(β) holds. According to assertion (i), the set C is equivalent to(26)x∣Φ-1αxTDDTx+z-∑i=1nxiμi,supβ≤0.
It is easy to know thatΦ-1(α)≥0 for α≥0.5, xTDDTx is a convex function with respect to x, and ∑i=1nxiμi,sup(β) is a linear function with respect to xi. According to the properties of convex function, {x∣Φ-1(α)xTDDTx+z-∑i=1nxiμi,sup(β)≤0} is a convex set. The proof of assertion (ii) is complete.Based on the above theoretical analysis, the following theorem gives the deterministic equivalent model of model (20).Theorem 6.
Assume thatη~N(μ,Σ), where Σ is a deterministic matrix, μi, i=1,2,…,n, are mutual independent fuzzy variables, and μTx is continuous. Then, in the case of α,β≥0.5, EV-ERV model (20) is equivalent to the following programming model: (27)max∑i=1nxiEμis.t.Φ-1αxTDDTx+z-∑i=1nxiμi,supβ≤0z≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.On the basis of Theorems5 and 6, we arrive at the following result.Theorem 7.
Assume thatη~N(μ,Σ), where Σ is a deterministic matrix, μi, i=1,2,…,n, are mutual independent fuzzy variables, and μTx is continuous. Then, in the case of α,β≥0.5, model (27) is a convex programming model.Proof.
In the case ofα,β≥0.5, by Theorem 5, C={x∣Cr{μTx≥Φ-1(α)xTDDTx+z}≥β} is a convex set. Then the feasible region is convex. Due to the linearity of the objective function and the convexity of the feasible region, model (27) is a convex programming model.


## 3.1. Computing Expected Return Rate

Since the return rates are denoted by the random fuzzy vectorη, for any given γ∈Γ, ηγ is the random vector (η1,γ,η2,γ,…,ηn,γ). Thus, the expected return rate is calculated by the following formula: (13)EηTx=EγEωηγTx=Eγ∑i=1nxiEωηi,γ.The first equality in (13) holds true due to the definition of expected value operator of random fuzzy variable [33], and the second equality in (13) is based on the linearity of expected value operator of random variable. Note that ∑i=1nxiE[ηi,γ] is the linear combination of fuzzy variables E[ηi,γ], i=1,2,…,n. To further compute Eγ[∑i=1nxiEω[ηi,γ]], the independence condition or comonotonicity condition for fuzzy variables Eω[ηi,γ], i=1,2,…,n is required. We next discuss the two situations, respectively.Case 1 (independence condition [35]).
According to [35], the expected value operator of fuzzy variable has the property of independence linearity. Based on this property, if E[η1,γ],E[η2,γ],…,E[ηn,γ] are mutual independent, one has(14)Eγ∑i=1nxiEωηi,γ=∑i=1nxiEγEωηi,γ=∑i=1nxiEξi.Case 2 (comonotonicity condition [34]).
According to [34, Theorem 1], the expected value operator of fuzzy variable has the property of comonotonic linearity. Based on this property, if E[η1,γ],E[η2,γ],…,E[ηn,γ] are comonotonic, one has(15)Eγ∑i=1nxiEωηi,γ=∑i=1nxiEγEωηi,γ=∑i=1nxiEξi.In the above two cases, the expected value operator of fuzzy variable has linear property,E[ηTx]=∑i=1nxiE[ηi]. In general case, we may compute Eγ[∑i=1nxiEω[ηi,γ] by using fuzzy simulation or approximation method [33, 36] based on the possibility distributions of fuzzy variables Eω[ηi,γ], i=1,2,…,n.


## 3.2. Processing Equilibrium Risk Value

In this subsection, we will handle the following probability constraint in model (12):(16)Pr⁡ηγTx≥z≥α,where the fuzzy parameter γ∈Γ is given in advance.It is assumed that random vectorηγ follows a multivariate normal distribution N(μγ,Σγ). The covariance matrix Σγ is a symmetric positive semidefinite matrix. So there exists a lower triangular and nonsingular matrix Dγ, μγ∈Rn, and a random vector ξ with components being independent standard normal random variables such that ηγ=Dγξ+μγ. Thus, one has E[ηγ]=μγ and Σγ=DγDγT. If Σγ is positive definite, then the multivariate normal distribution is nondegenerate. This case appears if and only if Dγ has full row rank. Otherwise the distribution is degenerate or singular. The multivariate normal distribution is uniquely determined by the expected value vector μγ and the covariance matrix Σγ.If the multivariate normal distribution is nondegenerate, then it is absolutely continuous. Thus, one hasηγTx~N(μγTx,xTDγDγTx).LetG(x)=Pr⁡{ηγTx≥z}. If xTDγDγTx>0, by standardization, G(x) has the following equivalent representation: (17)Pr⁡ηγTx≥z=Pr⁡ηγTx-μγTxxTDγDγTx≥z-μγTxxTDγDγTx=ΦμγTx-zxTDγDγTx.Taking all the cases ofG(x) into account, the following analytical expression is obtained: (18)Gx=1,ifxTDγDγTx=0,μγTx≥z,0,ifxTDγDγTx=0,μγTx<z,ΦμγTx-zxTDγDγTx,ifxTDγDγTx≠0.For any givenγ∈Γ, if xTDγDγTx≠0, then one has (19)Gx≥α⟺ΦμγTx-zxTDγDγTx≥α⟺μγTx-zxTDγDγTx≥Φ-1α⟺μγTx≥Φ-1αxTDγDγTx+z.From the above analysis, ifE[η1,γ],E[η2,γ],…,E[ηn,γ] are comonotonic or mutual independent, then model (12) can be turned into the following credibilistic programming problem: (20)max∑i=1nxiEμis.t.Crγ∣μγTx≥Φ-1αxTDγDγTx+z≥βz≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.


## 3.3. The Convexity of Feasible Region

Section3.2 has dealt with the probability risk. In this subsection, we will address the equivalent representation of credibilistic constraint and discuss the convexity of feasible region.LetC={x∣Cr{μTx≥Φ-1(α)xTDDTx+z}≥β}. In real-life portfolio selection problems, small confidence levels are meaningless. Thus, in this paper, we consider the case where α,β are in the interval [0.5,1].Theorem 5.
Assume thatη~N(μ,Σ), where Σ is deterministic matrix, μi, i=1,2,…,n, are mutual independent fuzzy variables, and μTx is continuous. In the case of α,β≥0.5, the following results hold.(i)
Cr{μTx≥Φ-1(α)xTDDTx+z}≥β⇔Φ-1(α)xTDDTx+z-(μTx)sup(β)≤0.(ii)
C={x∣Cr{μTx≥Φ-1(α)xTDDTx+z}≥β} is a convex set.Proof.
We first prove the necessity of assertion (i). By the definition of the optimistic value of fuzzy variableμTx, one has (21)μTxsupβ≥Φ-1αxTDDTx+z,which implies (22)Φ-1αxTDDTx+z-μTxsupβ≤0.The necessity of assertion (i) holds true.
We next prove the sufficiency of assertion (i). According to the definition of optimistic value, we know the credibilistic constraintCr{μTx≥(μTx)sup(β)}≥β holds. Since (23)μTxsupβ≥Φ-1αxTDDTx+zis equivalent to (24)Φ-1αxTDDTx+z-μTxsupβ≤0,we have (25)CrμTx≥Φ-1αxTDDTx+z≥βdue to the continuity of μTx. The sufficiency is proved.
In what follows, we prove assertion (ii).
According to the property of optimistic value, the equality(μTx)sup(β)=∑i=1nxiμi,sup(β) holds. According to assertion (i), the set C is equivalent to(26)x∣Φ-1αxTDDTx+z-∑i=1nxiμi,supβ≤0.
It is easy to know thatΦ-1(α)≥0 for α≥0.5, xTDDTx is a convex function with respect to x, and ∑i=1nxiμi,sup(β) is a linear function with respect to xi. According to the properties of convex function, {x∣Φ-1(α)xTDDTx+z-∑i=1nxiμi,sup(β)≤0} is a convex set. The proof of assertion (ii) is complete.Based on the above theoretical analysis, the following theorem gives the deterministic equivalent model of model (20).Theorem 6.
Assume thatη~N(μ,Σ), where Σ is a deterministic matrix, μi, i=1,2,…,n, are mutual independent fuzzy variables, and μTx is continuous. Then, in the case of α,β≥0.5, EV-ERV model (20) is equivalent to the following programming model: (27)max∑i=1nxiEμis.t.Φ-1αxTDDTx+z-∑i=1nxiμi,supβ≤0z≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.On the basis of Theorems5 and 6, we arrive at the following result.Theorem 7.
Assume thatη~N(μ,Σ), where Σ is a deterministic matrix, μi, i=1,2,…,n, are mutual independent fuzzy variables, and μTx is continuous. Then, in the case of α,β≥0.5, model (27) is a convex programming model.Proof.
In the case ofα,β≥0.5, by Theorem 5, C={x∣Cr{μTx≥Φ-1(α)xTDDTx+z}≥β} is a convex set. Then the feasible region is convex. Due to the linearity of the objective function and the convexity of the feasible region, model (27) is a convex programming model.


## 4. Equivalent Deterministic Convex Programming Models

In this section, we specify the possibility distributions of fuzzy parametersμi, i=1,2,…,n, so that the analytical expressions of E[μi] and μi,sup(β) are available. We consider the cases that fuzzy parameters μi, i=1,2,…,n, are characterized by trapezoidal, triangular, and normal fuzzy variables, respectively.Case 1 (trapezoidal fuzzy variables).
Letη~N(μ,Σ), μi=(ri(1),ri(2),ri(3),ri(4)), i=1,2,…,n, be mutual independent trapezoidal fuzzy variables. According to [34], the expected return of μTx is (28)EμTx=14∑i=1nri1xi+∑i=1nri2xi+∑i=1nri3xi+∑i=1nri4xi.
The optimistic value ofμTx is (29)∑i=1nxiμi,supβ=21-β∑i=1nri2xi+2β-1∑i=1nri1xi,where β≥0.5.
As a consequence, for any given parametersα,β≥0.5, by Theorem 6, model (27) is equivalent to the following deterministic convex programming model: (30)max14∑i=1nri1xi+∑i=1nri2xi+∑i=1nri3xi+∑i=1nri4xis.t.Φ-1αxTDDTx+z-21-β∑i=1nri2xi-2β-1∑i=1nri1xi≤0z≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.
Triangular fuzzy variable is a special case of trapezoidal fuzzy variable. The following remark is about the result related to triangular possibility distributions.Remark 8.
Letη~N(μ,Σ), μi=(ri(1),ri(2),ri(3)), i=1,2,…,n, be mutual independent triangular fuzzy variables. According to model (30), for any given parameters α,β≥0.5, the equivalent deterministic convex programming model of model (27) is (31)max14∑i=1nri1xi+2∑i=1nri2xi+∑i=1nri3xis.t.Φ-1αxTDDTx+z-21-β∑i=1nri2xi-2β-1∑i=1nri1xi≤0z≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.Case 2 (normal fuzzy variables).
Letη~N(μ,Σ), μi=n(mi,σi), i=1,2,…,n, be mutual independent normal fuzzy variables. According to [34], the expected return of μTx is E[μTx]=∑i=1nmixi. The optimistic value of μTx is (32)∑i=1nxiμi,supβ=∑i=1nmixi+-2ln⁡21-β∑i=1nxi2σi,where β≥0.5.
Therefore, according to Theorem6, for any given parameters α,β≥0.5, the equivalent deterministic convex programming model of model (27) reads (33)max∑i=1nmixis.t.Φ-1αxTDDTx+z-∑i=1nmixi--2ln⁡21-β∑i=1nxi2σi≤0z≥κ∑i=1nxi=1xi≥0,i=1,2,…,n.


## 5. Numerical Experiments and Comparison Study

In this section, we will conduct some numerical experiments to demonstrate the feasibility and effectiveness of the developed equilibrium optimization method. We first give some descriptions about our portfolio selection problem in the next subsection.


### 5.1. Problem Description

Assume that there are 20 potential risky assets for an investor. In this portfolio selection problem, the return rates have twofold uncertainty and are represented by random fuzzy variables. For the sake of presentation, we suppose that the prescribed confidence levelsα,β take their values from the interval [0.5,1].Letη=(η1,η2,…,η20)~N(μ,Σ) and μi=(ri(1),ri(2),ri(3),ri(4)), i=1,2,…,20, be mutual independent trapezoidal fuzzy variables. The possibility distributions of μi, i=1,2,…,20, are provided in Table 1, and the covariance matrix Σ=(σij)20×20 is expressed as 10-2(Σ1Σ2), where the matrices Σ1 and Σ2 are defined as follows: (34)Σ1=0.61980.11550.1096-0.06850.00380.12220.0049-0.07700.0319-0.03680.11550.5989-0.09370.07570.05310.1086-0.00410.05980.10670.03510.1096-0.09370.68240.02150.0389-0.0539-0.00040.0912-0.0550-0.0010-0.06850.07570.02150.64810.07180.01300.1047-0.0651-0.0494-0.06730.00380.05310.03890.07180.6583-0.1011-0.0604-0.0524-0.00970.04410.12220.1086-0.05390.0130-0.10110.50620.0014-0.0576-0.0785-0.01820.0049-0.0041-0.00040.1047-0.06040.00140.5376-0.0959-0.0997-0.0202-0.07700.05980.0912-0.0651-0.0524-0.0576-0.09590.75620.0317-0.08620.03190.1067-0.0550-0.0494-0.0097-0.0785-0.09970.03170.5243-0.0096-0.03680.0351-0.0010-0.06730.0441-0.0182-0.0202-0.0862-0.00960.81990.0536-0.05660.0058-0.0055-0.09260.1329-0.1258-0.0200-0.11450.12640.0596-0.07670.0744-0.0118-0.0230-0.0668-0.1022-0.17420.0360-0.06470.0506-0.1416-0.15210.05060.0718-0.0597-0.0279-0.05990.17100.00260.02780.0399-0.11900.1091-0.0933-0.0491-0.00450.0752-0.09810.02300.1505-0.03530.0332-0.0742-0.0080-0.03790.07430.06300.05470.11860.0356-0.01050.0371-0.02140.0032-0.1460-0.0696-0.05140.0008-0.06550.0969-0.0331-0.14520.01230.0205-0.09580.05670.0483-0.0491-0.0479-0.0039-0.0670-0.0077-0.0698-0.0421-0.05140.10340.03450.02410.0215-0.0289-0.06190.0332-0.0900-0.0145-0.0699-0.1032-0.0703-0.06120.03260.08250.01930.01450.04660.11610.1231-0.07210.02140.04320.0527,Σ2=0.05360.05960.05060.02780.15050.03560.0969-0.0039-0.02890.0825-0.0566-0.0767-0.14160.0399-0.0353-0.0105-0.0331-0.0670-0.06190.01930.00580.0744-0.1521-0.11900.03320.0371-0.1452-0.00770.03320.0145-0.0055-0.01180.05060.1091-0.0742-0.02140.0123-0.0698-0.09000.0466-0.0926-0.02300.0718-0.0933-0.00800.00320.0205-0.0421-0.01450.11610.1329-0.0668-0.0597-0.0491-0.0379-0.1460-0.0958-0.0514-0.06990.1231-0.1258-0.1022-0.0279-0.00450.0743-0.06960.05670.1034-0.1032-0.0721-0.0200-0.1742-0.05990.07520.0630-0.05140.04830.0345-0.07030.0214-0.11450.03600.1710-0.09810.05470.0008-0.04910.0241-0.06120.04320.1264-0.06470.00260.02300.1186-0.0655-0.04790.02150.03260.05270.4256-0.0703-0.0141-0.0424-0.11220.09070.01710.09200.04710.0074-0.07030.6835-0.06550.0812-0.0195-0.08300.04980.0173-0.04540.0768-0.0141-0.06550.58840.04050.06400.0214-0.10600.0359-0.0390-0.0281-0.04240.08120.04050.72220.02760.0692-0.0715-0.01770.12220.0913-0.1122-0.01950.06400.02760.6265-0.00790.0838-0.0921-0.0681-0.03010.0907-0.08300.02140.0692-0.00790.7695-0.02830.0249-0.06010.05670.01710.0498-0.1060-0.07150.0838-0.02830.75370.02800.06670.05500.09200.01730.0359-0.0177-0.09210.02490.02800.7757-0.10080.09960.0471-0.0454-0.03900.1222-0.0681-0.06010.0667-0.10080.5766-0.01010.00740.0768-0.02810.0913-0.03010.05670.05500.0996-0.01010.5674.Table 1
The distributions of trapezoidal fuzzy parametersμi.
Optional risky assets
Possibility distributions
1
μ
1
=
(
0.005,0.036,0.038,0.048
)
2
μ
2
=
(
0.006,0.037,0.038,0.045
)
3
μ
3
=
(
0.007,0.039,0.040,0.047
)
4
μ
4
=
(
0.004,0.038,0.040,0.051
)
5
μ
5
=
(
0.005,0.037,0.039,0.045
)
6
μ
6
=
(
0.006,0.041,0.042,0.049
)
7
μ
7
=
(
0.004,0.036,0.039,0.048
)
8
μ
8
=
(
0.005,0.038,0.039,0.045
)
9
μ
9
=
(
0.006,0.038,0.040,0.046
)
10
μ
10
=
(
0.005,0.038,0.040,0.046
)
11
μ
11
=
(
0.004,0.036,0.039,0.047
)
12
μ
12
=
(
0.005,0.037,0.0385,0.0445
)
13
μ
13
=
(
0.006,0.039,0.0405,0.0465
)
14
μ
14
=
(
0.004,0.038,0.041,0.050
)
15
μ
15
=
(
0.004,0.037,0.0392,0.0448
)
16
μ
16
=
(
0.005,0.041,0.0425,0.0485
)
17
μ
17
=
(
0.004,0.036,0.040,0.047
)
18
μ
18
=
(
0.004,0.038,0.0385,0.0455
)
19
μ
19
=
(
0.005,0.038,0.0395,0.0465
)
20
μ
20
=
(
0.004,0.038,0.0395,0.0464
)In this case, model (30) becomes the following deterministic convex programming model:(35)max0.03175x1+0.0315x2+0.03325x3+0.03325x4+0.0315x5+0.0345x6+0.03175x7+0.03175x8+0.0325x9+0.03225x10+0.0315x11+0.03125x12+0.033x13+0.03325x14+0.03125x15+0.03425x16+0.03175x17+0.0315x18+0.03225x19+0.031975x20s.t.10-1Φ-1α∑i=120∑j=120σijxixj1/2+z-21-βV2x-2β-1V1x≤0z≥κ∑i=120xi=1xi≥0,i=1,2,…,20,where the coefficient of xi in the objective function is the expected value of trapezoidal fuzzy variable μi, σij in the ERV constraint is the element of matrix Σ at the ith row and the jth column, and the analytical expressions of functions V1(x) and V2(x) in the ERV constraint are defined as follows: (36)V1x=∑i=120ri1xi=0.005x1+0.006x2+0.007x3+0.004x4+0.005x5+0.006x6+0.004x7+0.005x8+0.006x9+0.005x10+0.004x11+0.005x12+0.006x13+0.004x14+0.004x15+0.005x16+0.004x17+0.004x18+0.005x19+0.004x20,V2x=∑i=120ri2xi=0.036x1+0.037x2+0.039x3+0.038x4+0.037x5+0.041x6+0.036x7+0.038x8+0.038x9+0.038x10+0.036x11+0.037x12+0.039x13+0.038x14+0.037x15+0.041x16+0.036x17+0.038x18+0.038x19+0.038x20.


### 5.2. Computational Results

Firstly, we set the confidence levelsα=0.8, β=0.8 and the prescribed value of ERV κ=0.006 and solve model (35) by LINGO software. After 31 iterations, we obtain the distributive investment ratios among 20 risky assets, which are reported in Table 2. Due to the convexity of model (35), the obtained investment scheme is the global optimal solution, and its corresponding optimal objective value is 0.03293.Table 2
Computational results withα=0.8, β=0.8, and κ=0.006.
Investment ratios
x
1
x
2
x
3
x
4
x
5
x
6
x
7
x
8
x
9
x
10
x
11
x
12
x
13
x
14
x
15
x
16
x
17
x
18
x
19
x
20
Values
0
0
0.06599
0.03981
0.03508
0.19412
0.07107
0.06188
0.07594
0.04273
0
0.05974
0.07214
0.02827
0
0.11136
0.04525
0.00917
0.08748
0Secondly, considering that investors may have different attitudes towards risk, model (35) is solved under various values of parameters α, β, and κ. The obtained computational results are provided in Table 3, from which we observe the effect of model parameters on the solution quality of the EV-ERV model. The computational results demonstrate that when the values of parameters α, β, and κ are changed, the investment ratios and objective values are also changed accordingly. If we fix two of the parameters α, β, and κ, the objective value will decrease with respect to the remaining parameter. For example, when α=0.8 and κ=0.006, the objective values are EV=0.03329,0.03302,0.03261, corresponding to β=0.75,0.79,0.82, respectively. In addition, the value of z is always equal to that of κ; that is, the ERV is attainable and equal to κ, which is evident from theoretical analysis.Table 3
The optimal solutions of equilibrium optimization model under various values of parameters.
Parameters
Investment ratios
Objective values
α
β
κ
x
1
x
2
x
3
x
4
x
5
EV
x
6
x
7
x
8
x
9
x
10
x
11
x
12
x
13
x
14
x
15
x
16
x
17
x
18
x
19
x
20
0.78
0.8
0.006
0
0
0.08332
0.04127
0.02539
0.03308
0.21540
0.05473
0.04693
0.07931
0.03955
0
0.04148
0.07848
0.04074
0
0.12020
0.05187
0.00402
0.07731
0
0.8
0.8
0.006
0
0
0.06599
0.03981
0.03508
0.03293
0.19412
0.07107
0.06188
0.07594
0.04273
0
0.05974
0.07214
0.02827
0
0.11136
0.04525
0.00917
0.08748
0
0.82
0.8
0.006
0
0.01039
0.04343
0.03564
0.04671
0.03271
0.16033
0.09309
0.08176
0.06673
0.04612
0.00108
0.08616
0.06753
0.00916
0
0.09753
0.03613
0.01668
0.10153
0
0.8
0.75
0.006
0
0
0.10802
0.04491
0.01061
0.03329
0.24738
0.02916
0.02402
0.08360
0.03457
0
0.01334
0.08674
0.05983
0
0.13391
0.06153
0
0.06240
0
0.8
0.79
0.006
0
0
0.07626
0.04095
0.02916
0.03302
0.20703
0.06120
0.05276
0.07790
0.04081
0
0.04857
0.07578
0.03599
0
0.11685
0.04931
0.06131
0.08131
0
0.8
0.82
0.006
0
0.01762
0.03796
0.03247
0.05033
0.03261
0.14456
0.09980
0.08373
0.06536
0.04170
0.01818
0.09175
0.06540
0.00562
0.00613
0.08971
0.03125
0.01620
0.10224
0
0.78
0.78
0.006
0
0
0.10076
0.04330
0.01526
0.03323
0.23748
0.03741
0.03131
0.08249
0.03619
0
0.02234
0.08461
0.05375
0
0.12957
0.05859
0
0.06693
0
0.78
0.78
0.008
0
0
0.07275
0.04044
0.03112
0.03299
0.20291
0.06426
0.05606
0.07705
0.04148
0
0.05253
0.07474
0.03300
0
0.11520
0.04764
0.00727
0.08355
0
0.78
0.78
0.009
0
0
0.05231
0.03833
0.04272
0.03282
0.17766
0.08354
0.07411
0.07299
0.04528
0
0.07457
0.06753
0.01768
0
0.10462
0.03948
0.01339
0.09579
0Thirdly, we demonstrate the relationship between ERV and EV in our equilibrium optimization model via numerical experiments. For this purpose, we set the confidence levelsα=0.7 and β=0.8. The computational results are plotted in Figure 1, from which we observe that all the points under the line satisfy the constraints in our problem, and are called feasible solutions. The points on the line are the optimal solutions for the given α and β. The area under the line is the effective coverage.Figure 1
Relationship between ERV and EV underα=0.7 and β=0.8.


### 5.3. Comparing with Stochastic Optimization Method

In this subsection, we compare the proposed equilibrium optimization method with classical stochastic optimization method. For the sake of comparison, we also employ the data provided in Section5.1, and replace the trapezoidal fuzzy variables in Table 1 by their expected values. As a result, the return rate of each asset is characterized by random variable. In this case, the random fuzzy return vector η degenerates into a random return vector η~N(μ,Σ), where μ,Σ are deterministic vector and matrix. That is, the return rates ηi, i=1,2,…,20, are considered as normal random variables, and their dependence are characterized by covariance matrix Σ.By calculation, the expected return rate vector is(37)μ=0.03175,0.0315,0.03325,0.03325,0.0315,0.0345,0.03175,0.03175,0.0325,0.03225,0.0315,0.03125,0.033,0.03325,0.03125,0.03425,0.03175,0.0315,0.03225,0.031975T.The objective function isE[ηTx]=∑i=1nxiμi, and the ERV constraint is expressed as(38)Φ-1αxTΣx+z-μTx≤0.As a consequence, our equilibrium portfolio optimization problem reduces to the following convex programming model:(39)max0.03175x1+0.0315x2+0.03325x3+0.03325x400000+0.0315x5+0.0345x6+0.03175x7+0.03175x800000+0.0325x9+0.03225x10+0.0315x11+0.03125x1200000+0.033x13+0.03325x14+0.03125x15+0.03425x1600000+0.03175x17+0.0315x1800000+0.03225x19+0.031975x20s.t.10-1Φ-1α∑i=120∑j=120σijxixj1/200000+z-0.03175x1+0.0315x200000+0.03325x3+0.03325x4+0.0315x5+0.0345x600000+0.03175x7+0.03175x8+0.0325x9+0.03225x1000000+0.0315x11+0.03125x12+0.033x1300000+0.03325x14+0.03125x15+0.03425x16+0.03175x1700000+0.0315x18+0.03225x19+0.031975x20≤000000z≥κ00000∑i=120xi=100000xi≥0,i=1,2,…,20.To solve model (39) by LINGO software, we set the confidence level α=0.8 and the prescribed value of ERV κ=0.006. The obtained optimal portfolio is collected in Table 4, with the optimal objective value 0.03398.Table 4
Computational results withα=0.8 and κ=0.006.
Investment ratios
x
1
x
2
x
3
x
4
x
5
x
6
x
7
x
8
x
9
x
10
x
11
x
12
x
13
x
14
x
15
x
16
x
17
x
18
x
19
x
20
Values
0
0
0.11689
0.04596
0
0.41539
0
0
0.01806
0
0
0
0.08538
0.07634
0
0.24198
0
0
0
0To identify the influences of model parameters on solution quality, the stochastic model (39) is also solved under various values of parameters α and κ. The optimal investment solutions are reported in Table 5, from which we find that the solution results are different under various values of α and κ. If we set parameter κ=0.006, then the objective function is decreasing with respect to confidence level α. Given the parameter α, the objective function is also decreasing with respect to parameter κ.Table 5
The optimal solutions of random model under various values of parameters.
Parameters
Investment ratios
Objective values
α
κ
x
1
x
2
x
3
x
4
x
5
EV
x
6
x
7
x
8
x
9
x
10
x
11
x
12
x
13
x
14
x
15
x
16
x
17
x
18
x
19
x
20
0.78
0.006
0
0
0.09906
0.03658
0
0.03408
0.46099
0
0
0
0
0
0
0.06700
0.06212
0
0.27424
0
0
0
0
0.8
0.006
0
0
0.11689
0.04596
0
0.03398
0.41539
0
0
0.01806
0
0
0
0.08538
0.07634
0
0.24198
0
0
0
0
0.82
0.006
0
0
0.12286
0.04955
0
0.03386
0.37902
0
0
0.04157
0
0
0
0.08867
0.08131
0
0.21614
0.01161
0
0.00927
0
0.78
0.008
0
0
0.11519
0.04488
0
0.03399
0.42154
0
0
0.01284
0
0
0
0.08459
0.07455
0
0.24641
0
0
0
0
0.78
0.009
0
0
0.12044
0.04822
0
0.03394
0.40256
0
0
0.02893
0
0
0
0.08701
0.08010
0
0.23274
0
0
0
0We now compare the computational results reported in Tables2 and 4 as well as their optimal objective values. When model parameters are set to the same values, our equilibrium optimization method and stochastic optimization method produce different investment schemes. The equilibrium method selects 15 risky assets, while stochastic method selects 7 risky assets. In addition, two optimization methods suggest different investment ratios for the same risky asset. For example, the equilibrium method provides investment ratio 0.06599 for the third risky asset, while the stochastic method provides investment ratio x3=0.11689. From the viewpoint of distributive investment, we conclude that our equilibrium method is superior to the classical stochastic method.We continue to compare the computational results reported in Tables3 and 5, from which we observe that the equilibrium method and stochastic method provide different combinations of risky assets and different optimal return rates. For example, when we set probability level α=0.82, credibility level β=0.8, and κ=0.006, the equilibrium method provides the following optimal investment plan:(40)x=0,0.01039,0.04343,0.03564,0.04671,0.16033,0.09309,0.08176,0.06673,0.04612,0.00108,0.08616,0.06753,0.00916,0,0.09753,0.03613,0.01668,0.10153,0T,whose objective value is 0.03271. The stochastic method only selects 9 risky assets i, i=3,4,6,9,13,14,16,17,19. The corresponding investment ratios are 0.12286, 0.04955, 0.37902, 0.04157, 0.08867, 0.08131, 0.21614, 0.01161, and 0.00927, respectively. The optimal return rate provided by the stochastic method is 0.03386.On the other hand, under the same values of model parameters, we observe that the optimal objective values in Table5 are larger than those provided in Table 3. For example, when α=0.78, β=0.8, and κ=0.006, the equilibrium method provides the optimal value 0.03308, while stochastic method provides the optimal value 0.03408. This is because the credibilistic constraints are relaxed when the random fuzzy return rates degenerate into random ones. Therefore, we may obtain higher expected return rate. Nevertheless, the superiority of our equilibrium method is reflected in terms of diversification.Finally, we want to point out that an investor may often encounter a hybrid uncertain environment in modern financial markets. In this situation, the investor cannot ignore the influence of fuzzy uncertainty on the solution quality. The computational results support our arguments. For example, when probability levelα=0.78 and κ=0.006, the stochastic model (39) suggests the optimal investment scheme x∗=(0,0,0.09906,0.03658,0,0.46099,0,0,0,0,0,0,0.0670,0.06212,0,0.27424,0,0,0,0)T. However, x∗ is not the feasible solution to equilibrium optimization model (35) with probability level α=0.78, credibility level β=0.78, and κ=0.006. In order to keep x∗ feasible to model (35), we have to reduce the value of parameter κ from 0.006 to -0.00717. By the definition of κ, the optimal solution corresponding to -0.00717 is meaningless from a practical viewpoint.As a consequence, we conclude from the computational results that our equilibrium optimization method is effective in modeling practical portfolio selection problem under hybrid uncertain environment, where randomness and fuzziness are the state of affairs.


## 5.1. Problem Description

Assume that there are 20 potential risky assets for an investor. In this portfolio selection problem, the return rates have twofold uncertainty and are represented by random fuzzy variables. For the sake of presentation, we suppose that the prescribed confidence levelsα,β take their values from the interval [0.5,1].Letη=(η1,η2,…,η20)~N(μ,Σ) and μi=(ri(1),ri(2),ri(3),ri(4)), i=1,2,…,20, be mutual independent trapezoidal fuzzy variables. The possibility distributions of μi, i=1,2,…,20, are provided in Table 1, and the covariance matrix Σ=(σij)20×20 is expressed as 10-2(Σ1Σ2), where the matrices Σ1 and Σ2 are defined as follows: (34)Σ1=0.61980.11550.1096-0.06850.00380.12220.0049-0.07700.0319-0.03680.11550.5989-0.09370.07570.05310.1086-0.00410.05980.10670.03510.1096-0.09370.68240.02150.0389-0.0539-0.00040.0912-0.0550-0.0010-0.06850.07570.02150.64810.07180.01300.1047-0.0651-0.0494-0.06730.00380.05310.03890.07180.6583-0.1011-0.0604-0.0524-0.00970.04410.12220.1086-0.05390.0130-0.10110.50620.0014-0.0576-0.0785-0.01820.0049-0.0041-0.00040.1047-0.06040.00140.5376-0.0959-0.0997-0.0202-0.07700.05980.0912-0.0651-0.0524-0.0576-0.09590.75620.0317-0.08620.03190.1067-0.0550-0.0494-0.0097-0.0785-0.09970.03170.5243-0.0096-0.03680.0351-0.0010-0.06730.0441-0.0182-0.0202-0.0862-0.00960.81990.0536-0.05660.0058-0.0055-0.09260.1329-0.1258-0.0200-0.11450.12640.0596-0.07670.0744-0.0118-0.0230-0.0668-0.1022-0.17420.0360-0.06470.0506-0.1416-0.15210.05060.0718-0.0597-0.0279-0.05990.17100.00260.02780.0399-0.11900.1091-0.0933-0.0491-0.00450.0752-0.09810.02300.1505-0.03530.0332-0.0742-0.0080-0.03790.07430.06300.05470.11860.0356-0.01050.0371-0.02140.0032-0.1460-0.0696-0.05140.0008-0.06550.0969-0.0331-0.14520.01230.0205-0.09580.05670.0483-0.0491-0.0479-0.0039-0.0670-0.0077-0.0698-0.0421-0.05140.10340.03450.02410.0215-0.0289-0.06190.0332-0.0900-0.0145-0.0699-0.1032-0.0703-0.06120.03260.08250.01930.01450.04660.11610.1231-0.07210.02140.04320.0527,Σ2=0.05360.05960.05060.02780.15050.03560.0969-0.0039-0.02890.0825-0.0566-0.0767-0.14160.0399-0.0353-0.0105-0.0331-0.0670-0.06190.01930.00580.0744-0.1521-0.11900.03320.0371-0.1452-0.00770.03320.0145-0.0055-0.01180.05060.1091-0.0742-0.02140.0123-0.0698-0.09000.0466-0.0926-0.02300.0718-0.0933-0.00800.00320.0205-0.0421-0.01450.11610.1329-0.0668-0.0597-0.0491-0.0379-0.1460-0.0958-0.0514-0.06990.1231-0.1258-0.1022-0.0279-0.00450.0743-0.06960.05670.1034-0.1032-0.0721-0.0200-0.1742-0.05990.07520.0630-0.05140.04830.0345-0.07030.0214-0.11450.03600.1710-0.09810.05470.0008-0.04910.0241-0.06120.04320.1264-0.06470.00260.02300.1186-0.0655-0.04790.02150.03260.05270.4256-0.0703-0.0141-0.0424-0.11220.09070.01710.09200.04710.0074-0.07030.6835-0.06550.0812-0.0195-0.08300.04980.0173-0.04540.0768-0.0141-0.06550.58840.04050.06400.0214-0.10600.0359-0.0390-0.0281-0.04240.08120.04050.72220.02760.0692-0.0715-0.01770.12220.0913-0.1122-0.01950.06400.02760.6265-0.00790.0838-0.0921-0.0681-0.03010.0907-0.08300.02140.0692-0.00790.7695-0.02830.0249-0.06010.05670.01710.0498-0.1060-0.07150.0838-0.02830.75370.02800.06670.05500.09200.01730.0359-0.0177-0.09210.02490.02800.7757-0.10080.09960.0471-0.0454-0.03900.1222-0.0681-0.06010.0667-0.10080.5766-0.01010.00740.0768-0.02810.0913-0.03010.05670.05500.0996-0.01010.5674.Table 1
The distributions of trapezoidal fuzzy parametersμi.
Optional risky assets
Possibility distributions
1
μ
1
=
(
0.005,0.036,0.038,0.048
)
2
μ
2
=
(
0.006,0.037,0.038,0.045
)
3
μ
3
=
(
0.007,0.039,0.040,0.047
)
4
μ
4
=
(
0.004,0.038,0.040,0.051
)
5
μ
5
=
(
0.005,0.037,0.039,0.045
)
6
μ
6
=
(
0.006,0.041,0.042,0.049
)
7
μ
7
=
(
0.004,0.036,0.039,0.048
)
8
μ
8
=
(
0.005,0.038,0.039,0.045
)
9
μ
9
=
(
0.006,0.038,0.040,0.046
)
10
μ
10
=
(
0.005,0.038,0.040,0.046
)
11
μ
11
=
(
0.004,0.036,0.039,0.047
)
12
μ
12
=
(
0.005,0.037,0.0385,0.0445
)
13
μ
13
=
(
0.006,0.039,0.0405,0.0465
)
14
μ
14
=
(
0.004,0.038,0.041,0.050
)
15
μ
15
=
(
0.004,0.037,0.0392,0.0448
)
16
μ
16
=
(
0.005,0.041,0.0425,0.0485
)
17
μ
17
=
(
0.004,0.036,0.040,0.047
)
18
μ
18
=
(
0.004,0.038,0.0385,0.0455
)
19
μ
19
=
(
0.005,0.038,0.0395,0.0465
)
20
μ
20
=
(
0.004,0.038,0.0395,0.0464
)In this case, model (30) becomes the following deterministic convex programming model:(35)max0.03175x1+0.0315x2+0.03325x3+0.03325x4+0.0315x5+0.0345x6+0.03175x7+0.03175x8+0.0325x9+0.03225x10+0.0315x11+0.03125x12+0.033x13+0.03325x14+0.03125x15+0.03425x16+0.03175x17+0.0315x18+0.03225x19+0.031975x20s.t.10-1Φ-1α∑i=120∑j=120σijxixj1/2+z-21-βV2x-2β-1V1x≤0z≥κ∑i=120xi=1xi≥0,i=1,2,…,20,where the coefficient of xi in the objective function is the expected value of trapezoidal fuzzy variable μi, σij in the ERV constraint is the element of matrix Σ at the ith row and the jth column, and the analytical expressions of functions V1(x) and V2(x) in the ERV constraint are defined as follows: (36)V1x=∑i=120ri1xi=0.005x1+0.006x2+0.007x3+0.004x4+0.005x5+0.006x6+0.004x7+0.005x8+0.006x9+0.005x10+0.004x11+0.005x12+0.006x13+0.004x14+0.004x15+0.005x16+0.004x17+0.004x18+0.005x19+0.004x20,V2x=∑i=120ri2xi=0.036x1+0.037x2+0.039x3+0.038x4+0.037x5+0.041x6+0.036x7+0.038x8+0.038x9+0.038x10+0.036x11+0.037x12+0.039x13+0.038x14+0.037x15+0.041x16+0.036x17+0.038x18+0.038x19+0.038x20.


## 5.2. Computational Results

Firstly, we set the confidence levelsα=0.8, β=0.8 and the prescribed value of ERV κ=0.006 and solve model (35) by LINGO software. After 31 iterations, we obtain the distributive investment ratios among 20 risky assets, which are reported in Table 2. Due to the convexity of model (35), the obtained investment scheme is the global optimal solution, and its corresponding optimal objective value is 0.03293.Table 2
Computational results withα=0.8, β=0.8, and κ=0.006.
Investment ratios
x
1
x
2
x
3
x
4
x
5
x
6
x
7
x
8
x
9
x
10
x
11
x
12
x
13
x
14
x
15
x
16
x
17
x
18
x
19
x
20
Values
0
0
0.06599
0.03981
0.03508
0.19412
0.07107
0.06188
0.07594
0.04273
0
0.05974
0.07214
0.02827
0
0.11136
0.04525
0.00917
0.08748
0Secondly, considering that investors may have different attitudes towards risk, model (35) is solved under various values of parameters α, β, and κ. The obtained computational results are provided in Table 3, from which we observe the effect of model parameters on the solution quality of the EV-ERV model. The computational results demonstrate that when the values of parameters α, β, and κ are changed, the investment ratios and objective values are also changed accordingly. If we fix two of the parameters α, β, and κ, the objective value will decrease with respect to the remaining parameter. For example, when α=0.8 and κ=0.006, the objective values are EV=0.03329,0.03302,0.03261, corresponding to β=0.75,0.79,0.82, respectively. In addition, the value of z is always equal to that of κ; that is, the ERV is attainable and equal to κ, which is evident from theoretical analysis.Table 3
The optimal solutions of equilibrium optimization model under various values of parameters.
Parameters
Investment ratios
Objective values
α
β
κ
x
1
x
2
x
3
x
4
x
5
EV
x
6
x
7
x
8
x
9
x
10
x
11
x
12
x
13
x
14
x
15
x
16
x
17
x
18
x
19
x
20
0.78
0.8
0.006
0
0
0.08332
0.04127
0.02539
0.03308
0.21540
0.05473
0.04693
0.07931
0.03955
0
0.04148
0.07848
0.04074
0
0.12020
0.05187
0.00402
0.07731
0
0.8
0.8
0.006
0
0
0.06599
0.03981
0.03508
0.03293
0.19412
0.07107
0.06188
0.07594
0.04273
0
0.05974
0.07214
0.02827
0
0.11136
0.04525
0.00917
0.08748
0
0.82
0.8
0.006
0
0.01039
0.04343
0.03564
0.04671
0.03271
0.16033
0.09309
0.08176
0.06673
0.04612
0.00108
0.08616
0.06753
0.00916
0
0.09753
0.03613
0.01668
0.10153
0
0.8
0.75
0.006
0
0
0.10802
0.04491
0.01061
0.03329
0.24738
0.02916
0.02402
0.08360
0.03457
0
0.01334
0.08674
0.05983
0
0.13391
0.06153
0
0.06240
0
0.8
0.79
0.006
0
0
0.07626
0.04095
0.02916
0.03302
0.20703
0.06120
0.05276
0.07790
0.04081
0
0.04857
0.07578
0.03599
0
0.11685
0.04931
0.06131
0.08131
0
0.8
0.82
0.006
0
0.01762
0.03796
0.03247
0.05033
0.03261
0.14456
0.09980
0.08373
0.06536
0.04170
0.01818
0.09175
0.06540
0.00562
0.00613
0.08971
0.03125
0.01620
0.10224
0
0.78
0.78
0.006
0
0
0.10076
0.04330
0.01526
0.03323
0.23748
0.03741
0.03131
0.08249
0.03619
0
0.02234
0.08461
0.05375
0
0.12957
0.05859
0
0.06693
0
0.78
0.78
0.008
0
0
0.07275
0.04044
0.03112
0.03299
0.20291
0.06426
0.05606
0.07705
0.04148
0
0.05253
0.07474
0.03300
0
0.11520
0.04764
0.00727
0.08355
0
0.78
0.78
0.009
0
0
0.05231
0.03833
0.04272
0.03282
0.17766
0.08354
0.07411
0.07299
0.04528
0
0.07457
0.06753
0.01768
0
0.10462
0.03948
0.01339
0.09579
0Thirdly, we demonstrate the relationship between ERV and EV in our equilibrium optimization model via numerical experiments. For this purpose, we set the confidence levelsα=0.7 and β=0.8. The computational results are plotted in Figure 1, from which we observe that all the points under the line satisfy the constraints in our problem, and are called feasible solutions. The points on the line are the optimal solutions for the given α and β. The area under the line is the effective coverage.Figure 1
Relationship between ERV and EV underα=0.7 and β=0.8.


## 5.3. Comparing with Stochastic Optimization Method

In this subsection, we compare the proposed equilibrium optimization method with classical stochastic optimization method. For the sake of comparison, we also employ the data provided in Section5.1, and replace the trapezoidal fuzzy variables in Table 1 by their expected values. As a result, the return rate of each asset is characterized by random variable. In this case, the random fuzzy return vector η degenerates into a random return vector η~N(μ,Σ), where μ,Σ are deterministic vector and matrix. That is, the return rates ηi, i=1,2,…,20, are considered as normal random variables, and their dependence are characterized by covariance matrix Σ.By calculation, the expected return rate vector is(37)μ=0.03175,0.0315,0.03325,0.03325,0.0315,0.0345,0.03175,0.03175,0.0325,0.03225,0.0315,0.03125,0.033,0.03325,0.03125,0.03425,0.03175,0.0315,0.03225,0.031975T.The objective function isE[ηTx]=∑i=1nxiμi, and the ERV constraint is expressed as(38)Φ-1αxTΣx+z-μTx≤0.As a consequence, our equilibrium portfolio optimization problem reduces to the following convex programming model:(39)max0.03175x1+0.0315x2+0.03325x3+0.03325x400000+0.0315x5+0.0345x6+0.03175x7+0.03175x800000+0.0325x9+0.03225x10+0.0315x11+0.03125x1200000+0.033x13+0.03325x14+0.03125x15+0.03425x1600000+0.03175x17+0.0315x1800000+0.03225x19+0.031975x20s.t.10-1Φ-1α∑i=120∑j=120σijxixj1/200000+z-0.03175x1+0.0315x200000+0.03325x3+0.03325x4+0.0315x5+0.0345x600000+0.03175x7+0.03175x8+0.0325x9+0.03225x1000000+0.0315x11+0.03125x12+0.033x1300000+0.03325x14+0.03125x15+0.03425x16+0.03175x1700000+0.0315x18+0.03225x19+0.031975x20≤000000z≥κ00000∑i=120xi=100000xi≥0,i=1,2,…,20.To solve model (39) by LINGO software, we set the confidence level α=0.8 and the prescribed value of ERV κ=0.006. The obtained optimal portfolio is collected in Table 4, with the optimal objective value 0.03398.Table 4
Computational results withα=0.8 and κ=0.006.
Investment ratios
x
1
x
2
x
3
x
4
x
5
x
6
x
7
x
8
x
9
x
10
x
11
x
12
x
13
x
14
x
15
x
16
x
17
x
18
x
19
x
20
Values
0
0
0.11689
0.04596
0
0.41539
0
0
0.01806
0
0
0
0.08538
0.07634
0
0.24198
0
0
0
0To identify the influences of model parameters on solution quality, the stochastic model (39) is also solved under various values of parameters α and κ. The optimal investment solutions are reported in Table 5, from which we find that the solution results are different under various values of α and κ. If we set parameter κ=0.006, then the objective function is decreasing with respect to confidence level α. Given the parameter α, the objective function is also decreasing with respect to parameter κ.Table 5
The optimal solutions of random model under various values of parameters.
Parameters
Investment ratios
Objective values
α
κ
x
1
x
2
x
3
x
4
x
5
EV
x
6
x
7
x
8
x
9
x
10
x
11
x
12
x
13
x
14
x
15
x
16
x
17
x
18
x
19
x
20
0.78
0.006
0
0
0.09906
0.03658
0
0.03408
0.46099
0
0
0
0
0
0
0.06700
0.06212
0
0.27424
0
0
0
0
0.8
0.006
0
0
0.11689
0.04596
0
0.03398
0.41539
0
0
0.01806
0
0
0
0.08538
0.07634
0
0.24198
0
0
0
0
0.82
0.006
0
0
0.12286
0.04955
0
0.03386
0.37902
0
0
0.04157
0
0
0
0.08867
0.08131
0
0.21614
0.01161
0
0.00927
0
0.78
0.008
0
0
0.11519
0.04488
0
0.03399
0.42154
0
0
0.01284
0
0
0
0.08459
0.07455
0
0.24641
0
0
0
0
0.78
0.009
0
0
0.12044
0.04822
0
0.03394
0.40256
0
0
0.02893
0
0
0
0.08701
0.08010
0
0.23274
0
0
0
0We now compare the computational results reported in Tables2 and 4 as well as their optimal objective values. When model parameters are set to the same values, our equilibrium optimization method and stochastic optimization method produce different investment schemes. The equilibrium method selects 15 risky assets, while stochastic method selects 7 risky assets. In addition, two optimization methods suggest different investment ratios for the same risky asset. For example, the equilibrium method provides investment ratio 0.06599 for the third risky asset, while the stochastic method provides investment ratio x3=0.11689. From the viewpoint of distributive investment, we conclude that our equilibrium method is superior to the classical stochastic method.We continue to compare the computational results reported in Tables3 and 5, from which we observe that the equilibrium method and stochastic method provide different combinations of risky assets and different optimal return rates. For example, when we set probability level α=0.82, credibility level β=0.8, and κ=0.006, the equilibrium method provides the following optimal investment plan:(40)x=0,0.01039,0.04343,0.03564,0.04671,0.16033,0.09309,0.08176,0.06673,0.04612,0.00108,0.08616,0.06753,0.00916,0,0.09753,0.03613,0.01668,0.10153,0T,whose objective value is 0.03271. The stochastic method only selects 9 risky assets i, i=3,4,6,9,13,14,16,17,19. The corresponding investment ratios are 0.12286, 0.04955, 0.37902, 0.04157, 0.08867, 0.08131, 0.21614, 0.01161, and 0.00927, respectively. The optimal return rate provided by the stochastic method is 0.03386.On the other hand, under the same values of model parameters, we observe that the optimal objective values in Table5 are larger than those provided in Table 3. For example, when α=0.78, β=0.8, and κ=0.006, the equilibrium method provides the optimal value 0.03308, while stochastic method provides the optimal value 0.03408. This is because the credibilistic constraints are relaxed when the random fuzzy return rates degenerate into random ones. Therefore, we may obtain higher expected return rate. Nevertheless, the superiority of our equilibrium method is reflected in terms of diversification.Finally, we want to point out that an investor may often encounter a hybrid uncertain environment in modern financial markets. In this situation, the investor cannot ignore the influence of fuzzy uncertainty on the solution quality. The computational results support our arguments. For example, when probability levelα=0.78 and κ=0.006, the stochastic model (39) suggests the optimal investment scheme x∗=(0,0,0.09906,0.03658,0,0.46099,0,0,0,0,0,0,0.0670,0.06212,0,0.27424,0,0,0,0)T. However, x∗ is not the feasible solution to equilibrium optimization model (35) with probability level α=0.78, credibility level β=0.78, and κ=0.006. In order to keep x∗ feasible to model (35), we have to reduce the value of parameter κ from 0.006 to -0.00717. By the definition of κ, the optimal solution corresponding to -0.00717 is meaningless from a practical viewpoint.As a consequence, we conclude from the computational results that our equilibrium optimization method is effective in modeling practical portfolio selection problem under hybrid uncertain environment, where randomness and fuzziness are the state of affairs.


## 6. Conclusions

On the basis of probability and credibility measures, a new risk index called the ERV of random fuzzy variable was introduced. Under the equilibrium risk criterion, a new EV-ERV portfolio optimization model was built for portfolio selection problems, where the return rates are characterized by both probability distributions and possibility distributions. In the case where both subjective consciousness and objective factors affect the current financial markets, the developed equilibrium framework provided a novel optimization method for depicting real-life portfolio selection problem.When the randomness of uncertain return rates follows normal distributions, the proposed equilibrium portfolio selection model was turned into an equivalent credibilistic portfolio optimization model. The convexity of the credibilistic portfolio optimization model was discussed in Theorem7, which facilitates finding the desired global optimal portfolio. Furthermore, when the fuzziness of uncertain return rates follows trapezoidal, triangular, and normal distributions, the credibilistic portfolio optimization model was turned into its equivalent deterministic convex programming models.We compared the proposed equilibrium optimization method with traditional stochastic optimization method via a portfolio selection problem. The computational results demonstrated that both optimization methods can provide diversified investment schemes. However, the obtained equilibrium optimal solutions are more superior in terms of diversification. That is, when the fuzziness of uncertain return rates is considered, the equilibrium optimal solution usually diversified the optimal solutions obtained by stochastic method. As a consequence, when the exact probability distributions of return rates are unavailable, the proposed equilibrium optimization method provided an effective way to model practical portfolio selection problems with hybrid uncertain return rates.



---

*Source: 9461021-2016-03-07.xml*