# A Parametric Sharpe Ratio Optimization Approach for Fuzzy Portfolio Selection Problem

Authors: Ying Liu; Ya-Nan Li
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

When facing to make a portfolio decision, investors may care more about every portfolio’s performance on a return and risk trade-off. In this paper, a new low partial moment measurement that only punishes the loss risk is defined for selection variables based on L-S integral. Furthermore, a new performance measure for portfolio evaluation is proposed to generalize the Sharpe ratio in the fuzzy context. With the optimal performance criterion, a new parametric Sharpe ratio portfolio optimization model is developed wherein uncertain returns are presented as parametric interval-valued fuzzy variables. To make the proposed model easy to solve, we transform the fractional programming into an equivalent form and solve it with domain decomposition method (DDM). Finally, we apply the proposed performance measure into a portfolio selection problem, compare the computational results in different cases, and analyze the influence of different parameters on the optimal portfolio.


---

## Body

## 1. Introduction

Markowitz’s seminal paper on portfolio selection [1] has a major impact on financial research field and introduced a theory popularly referred to asmodern portfolio theory. The core of the theory is telling us the truth that investors should consider the risk and return together and determine the allocations of funds among investment alternatives on the basis of the trade-off between them [2]. In the past three decades, numerous researchers had contributed to the development of modern portfolio theory and proposed many optimization models and selection approaches for portfolio selection problem based on the trade-off analysis between reward and risk; see [3–7].The Sharpe ratio [8, 9], also calledreward-to-risk ratio, is often used to measure the risk-adjusted return of an investment asset or a portfolio. From the ratio, investors can distinguish how well the return of a portfolio compensates them for the risk taken. The expression of conventional Sharpe ratio was defined as the mean of excess return divided by its variance; therein the excess return is the differential return between a portfolio and the benchmark.A critique to this measure is that variance treats the upside and downside equally as risks and penalizes volatility with respect to expected value in both directions. However, positive excess returns (i.e., windfall profits) should not be regarded as something to avoid by investors. The occurrence of the low partial moment (LPM) [10, 11] is the development of downside risk measure in the portfolio research, and it also represents the attitude of decision maker from seeking to risk neutral to risk aversion. After that, many scholars have modified the Sharpe ratio with the LPM and presented a series of risk-adjusted performance measures. Sortino and van der Meer [12] presented the Sortino ratio as a variation of Sharpe ratio that adjusted the expected return for the risk of falling short of the risk-free return. Keating and Shadwick [13] proposed the Omega function as a natural performance measure to analyze returns distributions with the risk of the downside, lower partial moment, and gain-loss. Kaplan and Knowles [14] developed a Kappa ratio as a generalized downside risk-adjusted performance measurement based on the concept of LPM with parameter n. The three reward-to-risk ratios are all based on the partial moments, and the difference of them is Omega; Sortino and Kappa ratio use LPMs of orders 1, 2, and n, respectively. Value-at-Risk (VaR) and conditional Value-at-Risk (CVaR) are two downside risk measurements which directly gauge the investors’ losses. As an alternative risk measure, VaR [15] was incorporated into the ratio of the performance evaluation as the VaR ratio. Martin et al. [16] presented the CVaR ratio in which CVaR was substituted for variance. Chow and Lai [17] introduced the conditional Sharpe ratios (CSR) that measured lower partial risk-adjusted excess returns of an asset with respect to return distribution on the benchmark based on the conditional stochastic dominance. The authors summarized the monotonic transformations of Sharpe ratio with eleven downside risks in [18].It is worth noting that, in real life, the future asset return sometimes needs to be given by the expert’s judgement and estimation for no sufficient historical data or the noisy original data. Therefore, many scholars tried to model the uncertain parameter with fuzzy uncertainty [19–22], and a number of researchers have developed some downside risk measures in fuzzy decision problem [23, 24]. Since the ability to measure the risk-adjusted performance of a portfolio, the Sharpe ratio was introduced to assess the portfolio in fuzzy circumstance. Hasuike and Katagiri [25] considered a random fuzzy portfolio selection problem with general distributions and transformed it into the deterministic standard mathematical programming problem by introducing Sharpe ratio and fuzzy goals. Raei and Jahromi [26] proposed a model based on multiple criteria decision making and then by implementation of fuzzy analytical network analysis found appropriate weights of the asset. To measure the efficiency of the proposed model, the Sharpe ratio was computed for performance evaluation of the optimal portfolio. Based on [27], Nguyen et al. [28] modeled stock returns by the LR triangular fuzzy random variables and introduced a portfolio performance ratio named “reward-to-uncertainty” ratio to measure a portfolio performance.As the authors pointed in [29, 30], a common feature of most related studies is to assume the specific distributions of uncertain parameters known. However, in the financial market, the fully distributional information about the uncertain return affected by many factors is not always available. For this reason, it is reasonable to consider a more effective and flexible way to describe the uncertainty on return. Based on fuzzy possibility theory [31], Y. Liu and Y.-K. Liu [32] proposed the concept of parametric interval-valued fuzzy variable, which can handle the case of uncertainty contained in fuzzy possibility distribution. Specifically, the distribution of parametric interval-valued fuzzy variable is a variable possibility distribution. When the exact distribution of uncertain return is unavailable, the depictive approach of variable possibility distribution that depends on the choice of the parameter’s value is valid.Therefore, in the present paper, we consider a portfolio optimization problem in an incomplete market where there are multiple risk assets and one risk-free asset. We use a parametric optimization method to deal with portfolio performance evaluation; the uncertain return rate is characterized by variable possibility distribution and presented as a parametric interval-valued fuzzy variable. Based on L-S integral [33], we propose a new risk measure which meets the nature requirement only punishing the real risk which brings loss. We, therefore, modify the Sharpe ratio and formulate a parametric Sharpe ratio (PSR) to assess a portfolio in the fuzzy context. Subsequently, we proposed a fuzzy Sharpe ratio portfolio selection model with the viewpoint of portfolio performance evaluation. We also analyze the sensibility of parameters in the model and give a numerical example to illustrate the efficiency of the proposed model.The rest of this paper is organized as follows. For a better understanding, in Section2, we will briefly introduce the necessary knowledge about parametric interval-valued fuzzy variable and its selection variable. Then, based on L-S integral, the semimoment of selection variable is presented. According to the new risk measure, we refine the Sharpe ratio and a fuzzy parametric Sharpe ratio portfolio selection model will be developed in Section 3. In Section 4, we will derive the formulas of semimoment of excess return under the assumption that uncertain returns are parametric interval-valued trapezoidal/triangular fuzzy variables. After that, we analyze the proposed model, and the equivalent parametric forms of the model will be presented in Section 5. In Section 6, we will design a domain decomposition method to solve the model. For the sake of illustration, in Section 7, the application of the model will be discussed by a numerical example, and then we study the results with parameter analysis. Finally, some concluding remarks will be given in Section 8.


## 2. The Semimoment of Selection Variable

Let(Γ,P(Γ),P~os) be a fuzzy possibility space [31]. A map ξ~=(ξ~1,ξ~2,…,ξ~n):Γ↦Rn is called an n-ary type 2 (T2) fuzzy vector. When n=1, the map ξ~ is called a T2 fuzzy variable.The secondary possibility distribution functionμ~ξ(x) of ξ is defined as(1)μ~ξ~x=P~osγ∈Γ∣ξ~γ=x,x∈Rn,and type 2 possibility distribution function μξ~(x,u) of ξ is defined as(2)μξ~x,u=Posμ~ξ~x=u,x,u∈Rn×Jx, where Jx⊂[0,1] is the support of μ~ξ~(x).As a special case of T2 fuzzy variable, if for anyx∈R, u∈Jx⊆[0,1], the T2 possibility distribution function μξ~(x,u)=1, then ξ~ is an interval T2 fuzzy variable [32].Furthermore, if the secondary possibility distributionμ~ξ~(x) is a single subinterval of [0,1], which means if for any x∈R, μ~ξ~(x) is a subinterval [μξL(x;θl),μξU(x;θr)] of [0,1] with parameters θl,θr∈[0,1], then ξ~ is called a parametric interval-valued fuzzy variable [32].Ifθl=θr=0, then the reduced secondary possibility distribution is called the principle possibility distribution of ξ~, and the fuzzy variable associated with the principle possibility distribution is denoted by ξp [34].Letr1<r2≤r3<r4 be real numbers. Then a map ξ~ is called a parametric interval-valued trapezoidal fuzzy variable if its secondary possibility distribution μ~ξ~(r) is the subinterval(3)r-r1r2-r1-θlmin⁡r-r1r2-r1,r2-rr2-r1,r-r1r2-r1+θrmin⁡r-r1r2-r1,r2-rr2-r1 of [0,1] for r∈[r1,r2], the interval [1,1] for r∈[r2,r3], and the subinterval(4)r4-rr4-r3-θlmin⁡r4-rr4-r3,r-r3r4-r3,r4-rr4-r3+θrmin⁡r4-rr4-r3,r-r3r4-r3 of [0,1] for x∈[r3,r4], where θl,θr∈[0,1] are two parameters characterizing the degree of uncertainty that ξ~ takes on the value r. We denote the parametric interval-valued trapezoidal fuzzy variable ξ~ with the above distribution by [r1,r2,r3,r4;θl,θr].Letr1<r2<r3 be real numbers. Then a map ξ~ is called a parametric interval-valued triangular fuzzy variable if its secondary possibility distribution μ~ξ~(r) is the parametric subinterval(5)r-r1r2-r1-θlmin⁡r-r1r2-r1,r2-rr2-r1,r-r1r2-r1+θrmin⁡r-r1r2-r1,r2-rr2-r1 of [0,1] for r∈[r1,r2] and the parametric subinterval(6)r3-rr3-r2-θlmin⁡r3-rr3-r2,r-r2r3-r2,r3-rr3-r2+θrmin⁡r3-rr3-r2,r-r2r3-r2of [0,1] for r∈[r2,r3], where θ=(θl,θr), and θl,θr∈[0,1] are two parameters characterizing the degree of uncertainty that ξ~ takes on the value r. For simplicity, we denote the parametric interval-valued triangular fuzzy variable ξ~ with the above distribution by [r1,r2,r3;θl,θr].A fuzzy variableξ is called λ selection variable of parametric interval-valued fuzzy variable ξ~ [32], if ξ has the following parametric possibility distribution function:(7)μξr;θ=1-λμξLr;θl+λμξUr;θr,λ∈0,1,where θ=(θl,θr).Ifξ is a selection variable with finite expected value m, the definition of (ξ-m)- is(8)ξ-m-=m-ξ,ifξ≤m0,ifξ>m.Based on L-S integral [33], we define the semimoment of ξ as follows.Definition 1.
Ifξ is a selection variable with finite expected value m, then the integral(9)SMξ=∫-∞,+∞ξ-m-2dCrξ≤r=∫-∞,mm-ξ2dCrξ≤r,is called the semimoment of ξ, where Cr{ξ≤r} is the credibility distribution of ξ.The credibility distribution ofξ is defined in [35] and computed by(10)Crξ≤r=121+supt≤r⁡μξt;θ-supt>r⁡μξt;θ.


## 3. A Parametric Sharpe Ratio Portfolio Optimization Model

Assume an investor wants to make an allocation decision in a financial market, where there aren risk assets/securities and one risk-free asset. The investor wants a portfolio of the n risk assets plus the risk-free asset based on the portfolio performance ratio.The return of the risk-free asset is denoted byrf, and the investor can deposit his/her money or invest in treasury bond at this rate. Usually, rf is as a benchmark to the returns of risky assets. That means the investor demands a premium over the risk-free return by investing in risky assets, and the risk-free return is often less than risk return. If a risk-free asset performs better than the risky assets, it is advisable for investor to put all money on it.For various economic and noneconomic factors, the returns of risky assets are uncertain and their fixed possibility distributions are difficultly obtained in many situations. Thus, in this paper, we consider the case that the uncertain return of theith risky asset is presented as a parametric interval-valued fuzzy variable ξ~i(i=1,2,…,n), and the uncertainty return vector is ξ~=(ξ~1,…,ξ~n)T.Letxi be the percentage of wealth invested in the ith risky asset, and x=(x1,…,xn)T is the weight vector of a portfolio. Note that x does not have to sum to 1, and the remaining part 1-xT1 can be invested in the risk-free asset, where 1 is an n-dimension column vector and every component is 1. If short-selling does not exist in this asset market, the investment proportion x≥0. This is a common restriction on portfolio weights. In practice, for investment policy or legal reasons, or sometimes just because of particular asset classes, the investor cannot sell assets short.For a portfoliox, the total return rp is equal to the sum of the returns of two kinds of assets; that is,(11)rp=xTξ~+1-xT1rf.Using the risk-free returnrf as the benchmark, the excess return is(12)rp-rf=xTξ~+1-xT1rf-rf=xTξ~-rfxT1.Becauseξ~ is a parametric interval-valued fuzzy return vector, the excess return rate of a portfolio is a linear combination of parametric interval-valued fuzzy variables. According to [32], xTξ~-rfxT1 is still a parametric interval-valued fuzzy variable. The λ selection of the excess return rate is denoted as xTξ-rfxT1. Then the portfolio’s excess expected return rate is E[xTξ-rfxT1].To make a best trade-off between risk and return, the core problem is to find an appropriate risk measure according to the risk attitudes of different investors. In this paper, we use semimoment to measure the risk, and the semimoment of a portfolio’s excess return rate isSM[xTξ-rfxT1].In general, the population Sharpe ratio (SR) formula is demonstrated by(13)SR=excess  expected  returnunit  risk  deviation.With the expression above, we can formally define a new performance measure, that is, the fuzzy parametric Sharpe ratio (PSR) of excess return ratexTξ-rfxT1, as follows:(14)PSRx,ξ;λ=ExTξ-rfxT1SMxTξ-rfxT1.According to the different risk attitudes, there exist two optimization directions. We can choose maximizing the excess expected return per unit of risk or minimizing the risk per unit of return. In the present paper, we take the viewpoint of the investor who is willing to minimize the risk per unit of return. Thus, minimizing fuzzy parametric Sharpe ratio of excess return ratexTξ-rfxT1 is expressed as(15)minx⁡PSRx,ξ;λ=SMxTξ-rfxT1ExTξ-rfxT1.For a portfoliox, xTξ~ is a parametric interval-valued fuzzy variable and rfxT1 is a constant. If xTξ is the λ selection of xTξ~, by the definitions of λ selection and semimoment, we have (16)ExTξ-rfxT1=ExTξ-rfxT1,SMxTξ-rfxT1=SMxTξ.Thus, under the assumption of no short-selling, a fuzzy parametric Sharpe ratio portfolio selection model can be built as follows:(17)minSMxTξExTξ-rfxT1subject  toxT1≤1xi≥0,i=1,2,…,n.Model (17) is a fractional programming; the complexity of solution mainly focuses on the computation of the objective function. In the following section, we will discuss how to deal with this computational issue.


## 4. The Analytic Expression of Semimoment

In this section, we focus on the computational issue of Sharpe ratio of excess return rate. Under the assumption that the uncertain returns are some special parametric interval-valued fuzzy variables, we present the determinate equivalent parametric forms of semimoment and expected value.Suppose a fuzzy return vector isξ~=(ξ~1,…,ξ~n)T, where ξ~i=[r1i,r2i,r3i,r4i;θli,θri] are parametric interval-valued trapezoidal fuzzy variables and x=(x1,…,xn)T is the portfolio’s weight vector, xi∈[0,1]. Then the coefficient matrix of the fuzzy return vector is presented as (18)C=r11r21r31r41r12r22r32r42⋯⋯⋯⋯r1nr2nr3nr4nn×4,and denote cj=(rj1,…,rjn)T, j=1,2,3,4. Obviously, we have C=(c1,c2,c3,c4).By [32], the expected value of fuzzy return xTξ is(19)ExTξ=14∑i=1nxir1i+r2i+r3i+r4i+18λθr-1-λθl∑i=1nxir1i-r2i-r3i+r4i=14xTc1+c2+c3+c4+18λθr-1-λθlxTc1-c2-c3+c4,where parameters θl=max1≤j≤n⁡θlj and θr=min1≤j≤n⁡θrj.DenotingE[xTξ]=m, we can draw the following conclusions of semimoment.Theorem 2.
Let the uncertain return of asseti be parametric interval-valued trapezoidal fuzzy variable ξ~i=[r1i,r2i,r3i,r4i;θli,θri] and xi real numbers for i=1,2,…,n. If the principle possibility distributions of ξ~i are mutually independent, then the semimoment of fuzzy return sum xTξ is as follows.(i)
Ifm<xTc2, then(20)SMxTξ=12rTQr-1-λθr+1-λθl2xTc1+c2-c3-c4-λθr-1-λθlxTc1-c2-c3+c433072xTc2-c1,wherer=CTx=(∑i=1nxir1i,∑i=1nxir2i,∑i=1nxir3i,∑i=1nxir4i)T, and the elements of the symmetric matrix Q include(21)Q11=1128λθr-1-λθl3-7128λθr-1-λθl2+796λθr-1-λθl+1396,Q12=-132λθr-1-λθl3+332λθr-1-λθl2+112λθr-1-λθl+124,Q13=-132λθr-1-λθl3+532λθr-1-λθl2-116λθr-1-λθl+14,Q14=132λθr-1-λθl3-132λθr-1-λθl2-516λθr-1-λθl-14,Q22=1128λθr-1-λθl3+1128λθr-1-λθl2+196λθr-1-λθl+196,Q23=132λθr-1-λθl3+732λθr-1-λθl2+316λθr-1-λθl,Q24=-132λθr-1-λθl3-332λθr-1-λθl2-116λθr-1-λθl,Q33=1128λθr-1-λθl3-3128λθr-1-λθl2+132,Q34=-132λθr-1-λθl3-132λθr-1-λθl2+18λθr-1-λθl+18,Q44=1128λθr-1-λθl3+5128λθr-1-λθl2+116λθr-1-λθl+132.(ii)
IfxTc2≤m<xTc3, then(22)SMxTξ=12rTPr,wherer=CTx=(∑i=1nxir1i,∑i=1nxir2i,∑i=1nxir3i,∑i=1nxir4i)T, and the elements of the symmetric matrix P include(23)P11=-364λθr-1-λθl2+116λθr-1-λθl+748,P12=316λθr-1-λθl2-112,P13=P24=116λθr-1-λθl2-14,P14=-116λθr-1-λθl2-14λθr-1-λθl-14,P22=-364λθr-1-λθl2-116λθr-1-λθl-748,P23=-116λθr-1-λθl2+34λθr-1-λθl-14,P33=164λθr-1-λθl2-116λθr-1-λθl+132,P34=-116λθr-1-λθl2+14,P44=164λθr-1-λθl2+116λθr-1-λθl+132.(iii)
Ifm≥xTc3, then(24)SMxTξ=12rTPr-1-λθr+1-λθl2xTc1+c2-3c3+c4-λθr-1-λθlxTc1-c2-c3+c433072xTc4-c3,wherer=CTx=(∑i=1nxir1i,∑i=1nxir2i,∑i=1nxir3i,∑i=1nxir4i)T, and the elements of the symmetric matrix P are the same as case (ii), with the parameters θl=max1≤i≤n⁡θli and θr=min1≤i≤n⁡θri.
Moreover, the semimomentSM[xTξ] in case (ii) is a parametric quadratic convex function with respect to decision vector x∈Rn.On the other side, if fuzzy returnsξ~i=[r1i,r2i,r3i;θli,θri] are parametric interval-valued triangular fuzzy variables and x=(x1,…,xn)T is the portfolio’s weight vector, xi∈[0,1], then the coefficient matrix of the fuzzy return vector is presented as (25)C=r11r21r31r12r22r32⋯⋯⋯r1nr2nr3nn×3,and denote cj=(rj1,…,rjn)T, j=1,2,3. Obviously, we have C=(c1,c2,c3).By [32], the expected value of λ selection variable xTξ is(26)ExTξ=∑i=1nxir1i+2r2i+r3i4+λθr-1-λθl∑i=1nxir1i-2r2i+r3i8=xTc1+2c2+c34+λθr-1-λθlxTc1-2c2+c38,where the parameters θl=max1≤j≤n⁡θlj and θr=min1≤j≤n⁡θrj. Denoting E[xTξ]=m, we have the following corollary.Corollary 3.
Let the uncertain return of asseti be parametric interval-valued triangular fuzzy variable ξ~i=[r1i,r2i,r3i;θli,θri] and xi real numbers for i=1,2,…,n. If the principle possibility distributions of ξ~i are mutually independent, then the semimoment of fuzzy return sum xTξ is as follows. (i)
Ifm<xTc2, then(27)SMxTξ=12rTRr-1-λθr+1-λθl2xTc1-c3-λθr+1-λθlxTc1-2c2+c333072xTc2-c1,wherer=CTx=(∑i=1nxir1i,∑i=1nxir2i,∑i=1nxir3i)T, and the elements of the symmetric matrix R include (28)
R
11
=
1
128
λ
θ
r
-
1
-
λ
θ
l
3
-
7
128
λ
θ
r
-
1
-
λ
θ
l
2
+
7
96
λ
θ
r
-
1
-
λ
θ
l
+
13
96
,
R
12
=
-
1
16
λ
θ
r
-
1
-
λ
θ
l
3
+
1
4
λ
θ
r
-
1
-
λ
θ
l
2
+
1
48
λ
θ
r
-
1
-
λ
θ
l
-
7
24
,
R
13
=
1
32
λ
θ
r
-
1
-
λ
θ
l
3
-
1
32
λ
θ
r
-
1
-
λ
θ
l
2
-
5
16
λ
θ
r
-
1
-
λ
θ
l
-
1
4
,
R
22
=
1
32
λ
θ
r
-
1
-
λ
θ
l
3
-
1
32
λ
θ
r
-
1
-
λ
θ
l
2
-
1
48
λ
θ
r
-
1
-
λ
θ
l
+
1
24
,
R
23
=
-
1
16
λ
θ
r
-
1
-
λ
θ
l
3
-
1
8
λ
θ
r
-
1
-
λ
θ
l
2
+
1
16
λ
θ
r
-
1
-
λ
θ
l
+
1
8
,
R
33
=
1
128
λ
θ
r
-
1
-
λ
θ
l
3
+
5
128
λ
θ
r
-
1
-
λ
θ
l
2
+
1
16
λ
θ
r
-
1
-
λ
θ
l
+
1
32
.(ii)
Ifm≥xTc2, then(29)SMxTξ=12rTSr-1-λθr+1-λθl2+λθr-1-λθl3xTc1-2c2+c333072xTc3-c2,wherer=CTx=(∑i=1nxir1i,∑i=1nxir2i,∑i=1nxir3i)T, and the elements of the symmetric matrix S include(30)S11=-364λθr-1-λθl2+116λθr-1-λθl+748,S12=14λθr-1-λθl2-13,S13=-116λθr-1-λθl2-14λθr-1-λθl-14,S22=2-116λθr-1-λθl2+112,S23=0,S33=164λθr-1-λθl2+116λθr-1-λθl+116,with the parameters θl=max1≤i≤n⁡θli and θr=min1≤i≤n⁡θri.


## 5. The Equivalent Parametric Programming

In this section, we firstly illustrate how the fractional formulation of fuzzy parametric Sharpe ratio portfolio selection model can be transformed into an equivalent nonlinear form and then derive its equivalent parametric submodels.


### 5.1. The Transformation of Fuzzy Sharpe Ratio Model

In mathematical programming, if fractional programming meets the four conditions(A1)
the feasible regionB is bounded and nonempty;(A2)
for everyx∈B, g(x)>0;(A3)
for everyx∈B, f(x)>0;(A4)
g(x) is a continuous and positively homogeneous risk measure, then fractional programming(31)mingxfxsubject  tox∈Bcan be transformed into easily solving programming(32)mingy,tsubject  tofy,t≥1Ay-bt=0y≥0t>0,where t=1/f(x)>0 and y=tx (the Charnes-Cooper transformation [36]).Subsequently, we will check the conditions of model (17).(A1)
x is the weight vector of a portfolio, and xi≥0,i=1,2,…,n, so B is bounded and not negative.(A2)
For anyx, obviously, g(x)=SM[xTξ]>0.(A3)
The aim of investment is to gain the high return, so the return of risk asset should be higher than that of risk-free asset; namely, for anyx, E[xTξ]-rfxT1>0.(A4)
By the definition of semimoment,g(x)=SM[xTξ] is continuous and positive homogeneous.So far, we verify that the model meets four terms, and that means model (17) can be translated to its equivalent form:(33)minSMyTξsubject  toEyTξ-rfyT1-1≥0yT1-t≤0yi≥0,i=1,2,…,nt>0.According to [37], we can derive the relationship of optimal solutions between model (17) and model (33) easily. The decision variable of model (17) is x∈Rn, while the decision variables of model (33) are y∈Rn and t. If (y∗,t∗)∈Rn+1 is a globally optimal solution for model (33), then x∗=y∗/t∗ is a globally optimal solution for model (17).


### 5.2. The Equivalent Parametric Submodels

In the previous subsection, fuzzy Sharpe ratio portfolio selection model (17) has been transformed into equivalent nonfractional programming (33). In the following, we will discuss the equivalent parametric submodels of model (33) under the assumption that uncertain returns are parametric interval-valued trapezoidal fuzzy variables.Based on Theorem2 and (19), if E[xTξ]<xTc2, we have (34)ExTξ-rfxT1<xTc2-rfxT1.Let E[xTξ]-rfxT1=1/t, by the discussion in Section 5.1; then yTc2-rfyT1>1.Thus, ifE[xTξ]<xTc2, that is, yTc2-rfyT1>1, we can obtain the equivalent parametric programming of model (33) as follows:(35)min12yTCQCTy-1-λθr+1-λθl2yTc1+c2-c3-c4-λθr-1-λθlyTc1-c2-c3+c433072yTc2-c1subject  toyTc1+c2+c3+c44+λθr-1-λθlyTc1-c2-c3+c48-rfyT1-1≥0yTc2>1+rfyT1yT1-t≤0yi≥0,i=1,2,…,nt>0.Analogously, ifxTc2≤E[xTξ]<xTc3, that is, yTc2≤1+rfyT1<yTc3, then model (33) can be transformed into the following equivalent parametric programming:(36)min12yTCPCTysubject  toyTc1+c2+c3+c44+λθr-1-λθlyTc1-c2-c3+c48-rfyT1-1≥0yTc2≤1+rfyT1<yTc3yT1-t≤0yi≥0,i=1,2,…,nt>0.IfE[xTξ]≥xTc3, that is, yTc3≤1+rfyT1, then we can solve model (33) by solving the following equivalent parametric programming:(37)min12yTCPCTy-1-λθr+1-λθl2xTc1+c2-3c3+c4-λθr-1-λθlxTc1-c2-c3+c433072xTc4-c3subject  toyTc1+c2+c3+c44+λθr-1-λθlyTc1-c2-c3+c48-rfyT1-1≥0yTc3≤1+rfyT1yT1-t≤0yi≥0,i=1,2,…,nt>0.Models (35)–(37) are three equivalent parametric submodels of model (33) in the different partitions of feasible region. For any given parameter values, based on the previous discussion, submodels (35) and (37) are determinate nonlinear programming. Meanwhile, submodel (36) with given parameters is a quadratic convex programming. They all can be solved by standard software.


## 5.1. The Transformation of Fuzzy Sharpe Ratio Model

In mathematical programming, if fractional programming meets the four conditions(A1)
the feasible regionB is bounded and nonempty;(A2)
for everyx∈B, g(x)>0;(A3)
for everyx∈B, f(x)>0;(A4)
g(x) is a continuous and positively homogeneous risk measure, then fractional programming(31)mingxfxsubject  tox∈Bcan be transformed into easily solving programming(32)mingy,tsubject  tofy,t≥1Ay-bt=0y≥0t>0,where t=1/f(x)>0 and y=tx (the Charnes-Cooper transformation [36]).Subsequently, we will check the conditions of model (17).(A1)
x is the weight vector of a portfolio, and xi≥0,i=1,2,…,n, so B is bounded and not negative.(A2)
For anyx, obviously, g(x)=SM[xTξ]>0.(A3)
The aim of investment is to gain the high return, so the return of risk asset should be higher than that of risk-free asset; namely, for anyx, E[xTξ]-rfxT1>0.(A4)
By the definition of semimoment,g(x)=SM[xTξ] is continuous and positive homogeneous.So far, we verify that the model meets four terms, and that means model (17) can be translated to its equivalent form:(33)minSMyTξsubject  toEyTξ-rfyT1-1≥0yT1-t≤0yi≥0,i=1,2,…,nt>0.According to [37], we can derive the relationship of optimal solutions between model (17) and model (33) easily. The decision variable of model (17) is x∈Rn, while the decision variables of model (33) are y∈Rn and t. If (y∗,t∗)∈Rn+1 is a globally optimal solution for model (33), then x∗=y∗/t∗ is a globally optimal solution for model (17).


## 5.2. The Equivalent Parametric Submodels

In the previous subsection, fuzzy Sharpe ratio portfolio selection model (17) has been transformed into equivalent nonfractional programming (33). In the following, we will discuss the equivalent parametric submodels of model (33) under the assumption that uncertain returns are parametric interval-valued trapezoidal fuzzy variables.Based on Theorem2 and (19), if E[xTξ]<xTc2, we have (34)ExTξ-rfxT1<xTc2-rfxT1.Let E[xTξ]-rfxT1=1/t, by the discussion in Section 5.1; then yTc2-rfyT1>1.Thus, ifE[xTξ]<xTc2, that is, yTc2-rfyT1>1, we can obtain the equivalent parametric programming of model (33) as follows:(35)min12yTCQCTy-1-λθr+1-λθl2yTc1+c2-c3-c4-λθr-1-λθlyTc1-c2-c3+c433072yTc2-c1subject  toyTc1+c2+c3+c44+λθr-1-λθlyTc1-c2-c3+c48-rfyT1-1≥0yTc2>1+rfyT1yT1-t≤0yi≥0,i=1,2,…,nt>0.Analogously, ifxTc2≤E[xTξ]<xTc3, that is, yTc2≤1+rfyT1<yTc3, then model (33) can be transformed into the following equivalent parametric programming:(36)min12yTCPCTysubject  toyTc1+c2+c3+c44+λθr-1-λθlyTc1-c2-c3+c48-rfyT1-1≥0yTc2≤1+rfyT1<yTc3yT1-t≤0yi≥0,i=1,2,…,nt>0.IfE[xTξ]≥xTc3, that is, yTc3≤1+rfyT1, then we can solve model (33) by solving the following equivalent parametric programming:(37)min12yTCPCTy-1-λθr+1-λθl2xTc1+c2-3c3+c4-λθr-1-λθlxTc1-c2-c3+c433072xTc4-c3subject  toyTc1+c2+c3+c44+λθr-1-λθlyTc1-c2-c3+c48-rfyT1-1≥0yTc3≤1+rfyT1yT1-t≤0yi≥0,i=1,2,…,nt>0.Models (35)–(37) are three equivalent parametric submodels of model (33) in the different partitions of feasible region. For any given parameter values, based on the previous discussion, submodels (35) and (37) are determinate nonlinear programming. Meanwhile, submodel (36) with given parameters is a quadratic convex programming. They all can be solved by standard software.


## 6. The Domain Decomposition Method

The domain decomposition method (DDM) is a decomposition-based computation technique. It searches every partition of the feasible region to find a global optimum of the optimization problem. To begin with, the feasible region should be reasonably split into multiple subregions using the different cases existing in the original problem or different values of parameters. In each individual subregion, the subproblem corresponding to the original problem is not identical, and the differences between subproblems lie in the objective function or constraints. Next, to find the optimal solution of the original problem, it is necessary to design or select an appropriate algorithm and solve the subproblems according to the characteristics of the subproblems. Finally, we obtain the global optimum of original problem by comparing the obtained local optimal values of all subproblems. To be specific, the solution process mainly comprises the following three steps.Step 1.
Solve three equivalent parametric submodels (35)–(37), and find the local optimal solutions.Step 2.
Compute the objective values of three submodels (35)–(37) corresponding to these local optimal solutions, respectively.Step 3.
Compare the obtained optimal values of submodels (35)–(37), and choose the minimum as the global optimal value of original model (17).As mentioned in Section5.2, the feasible region of model (33) is decomposed into three disjoint subregions. Models (35)–(37) are the equivalent parametric subproblems of model (33) in the different subregions, and for any given parameters θ=(θl,θr) and λ they just determinate nonlinear programming which can be solved by classic algorithms. Lingo is a comprehensive tool designed to build and solve linear, nonlinear (convex or nonconvex), quadratic, and integer optimization models. However, the ability to solve programming is also constrained by the number of variables. Therefore, it is necessary to develop a more targeted efficient algorithm to the proposed models. In this paper, we employ the DDM combined Lingo software to solve model (33). In the next section, we will demonstrate the validity and practicability of this method via numerical experiments.


## 7. Numerical Example

Consider an asset market which has twenty risk assets and a risk-free asset, and an investor wants to decide on an appropriate mix of assets as his/her investment portfolio. The risk-free asset is usually a treasury bond with the fixed return raterf. The short-selling is assumed to be prohibited for any assets in this market. Due to the impacts of the economic and political factors, the returns of these assets are uncertain, and their exact possibility distributions cannot be obtained. Thus, in this paper, these twenty risk assets’ returns are presented by parametric interval-valued trapezoidal fuzzy variables ξ~i=[r1i,r2i,r3i,r4i;θli,θri],i=1,…,20 with variable possibility distributions, and the uncertain returns of twenty candidate assets are given in Table 1.Table 1
The parametric interval-valued trapezoidal fuzzy returns of twenty assets.
Asseti Parametric interval-valued fuzzy return 1 [0.9946, 0.9967, 1.0012, 1.0061;θl1,θr1] 2 [0.9911, 1.0002, 1.0019, 1.0052;θl2,θr2] 3 [0.9966, 1.0073, 1.0081, 1.0109;θl3,θr3] 4 [0.9983, 1.0090, 1.0122, 1.0135;θl4,θr4] 5 [1.0059, 1.0122, 1.0262, 1.0279;θl5,θr5] 6 [1.0099, 1.0159, 1.0248, 1.0429;θl6,θr6] 7 [1.0209, 1.0225, 1.0416, 1.0553;θl7,θr7] 8 [1.0291, 1.0299, 1.0468, 1.0678;θl8,θr8] 9 [1.0201, 1.0504, 1.0651, 1.0709;θl9,θr9] 10 [1.0305, 1.0399, 1.0469, 1.1093;θl10,θr10] 11 [1.0304, 1.0468, 1.0712, 1.0908;θl11,θr11] 12 [1.0365, 1.0529, 1.0698, 1.1096;θl12,θr12] 13 [1.0214, 1.0699, 1.0778, 1.1004;θl13,θr13] 14 [1.0391, 1.0699, 1.0969, 1.1096;θl14,θr14] 15 [1.0422, 1.0766, 1.0977, 1.1178;θl15,θr15] 16 [1.0373, 1.0914, 1.0972, 1.1171;θl16,θr16] 17 [1.0604, 1.0786, 1.1183, 1.1282;θl17,θr17] 18 [1.0460, 1.0952, 1.1098, 1.1361;θl18,θr18] 19 [1.0415, 1.0885, 1.1155, 1.1475;θl19,θr19] 20 [1.0326, 1.0896, 1.1330, 1.1392;θl20,θr20]


### 7.1. Computational Result under Parametric Interval-Valued Fuzzy Returns

Firstly, when the exact possibility distributions of asset returns are unavailable, we conduct the experiments under variable possibility distribution to find the optimal portfolio selection policy.In our experiments, we set the values of parametersθli and θri, i=1,2,…,20 as (38)θl1,θl2,…,θl20=0.8,0.2,0.3,0.9,0.7,0.5,0.4,0.6,0.8,0.1,0.6,0.3,0.5,0.3,0.7,0.4,0.2,0.8,0.5,0.2,θr1,θr2,…,θr20=0.3,0.8,0.2,0.6,0.3,0.5,0.3,0.9,0.7,0.5,0.4,0.6,0.7,0.4,0.9,0.8,0.5,0.2,0.8,0.2. Since θl=max1≤i≤nθli and θr=min1≤i≤nθri, we have θl=0.9, θr=0.2.Assume the investor selects the values ofλ as 0, 0.25, 0.5, 0.75, and 1. To test the effect of the risk-free return rate on the portfolio decision, we adjust the value of rf and do some experiments, and the returns of the risk-free asset are set as follows rf= 0.33%, 0.93%, 1.93%, 4.38%, 8.48%, 9.25%, 9.40%, 9.65%, and 9.70%. The obtained optimal allocation proportions to twenty risk assets are reported in Tables 2–6. Note that, according to x∗=y∗/t∗, we have the allocation proportions xi easily by reducing yi with the same ratio t. For simplicity, we give directly the results of xi.Table 2
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=0 (%).
r f x 1 x 3 x 4 x 8 x 16 x 17 x 18 x 19 0.0033 37.08982 51.43058 11.47960 0.000000 0.000000 0.000000 0.000000 0.000000 0.0093 0.000000 53.50372 37.74321 8.753070 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 47.97375 36.58445 0.000000 2.875511 12.56628 0.000000 0.000000 0.0438 0.000000 27.74888 26.18283 0.000000 25.06405 21.00424 0.000000 0.000000 0.0848 0.000000 0.000000 22.61346 0.000000 62.15479 35.32249 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 11.84791 13.55289 74.59920 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 0.000000 3.144038 86.73522 10.12074 0.0965 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.145002 97.85500 0.0970 — — — — — — — —Table 3
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=0.25 (%).
r f x 1 x 3 x 4 x 8 x 16 x 17 x 18 x 19 0.0033 38.99274 60.41775 0.000000 0.589515 0.000000 0.000000 0.000000 0.000000 0.0093 5.941417 61.16051 22.99126 9.906809 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 50.47323 34.22256 0.000000 3.955556 11.34866 0.000000 0.000000 0.0438 0.000000 53.58670 0.000000 0.000000 26.61964 19.79366 0.000000 0.000000 0.0848 0.000000 2.979972 0.000000 0.000000 65.03688 31.98315 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 17.58623 12.23930 70.17447 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 3.478543 6.432623 90.08883 0.000000 0.0965 0.000000 0.000000 0.000000 0.000000 0.000000 28.36709 7.005321 64.62759 0.0970 0.000000 0.000000 0.000000 0.000000 0.000000 5.780347 0.000000 94.21965Table 4
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=0.5 (%).
rf x 1 x 3 x 4 x 8 x 16 x17 x18 x19 0.0033 38.02403 61.84081 0.000000 0.135164 0.000000 0.000000 0.000000 0.000000 0.0093 3.779654 85.59803 0.000000 10.62232 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 84.02355 0.000000 0.000000 4.756157 11.22030 0.000000 0.000000 0.0438 0.000000 53.91651 0.000000 0.000000 28.45683 17.62665 0.000000 0.000000 0.0848 0.000000 3.533319 0.000000 0.000000 68.11919 28.34749 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 23.96391 10.77935 65.25674 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 9.371306 5.083677 85.54502 0.000000 0.0965 0.000000 0.000000 0.000000 0.000000 0.000000 69.14263 0.000000 30.85737 0.0970 0.000000 0.000000 0.000000 0.000000 0.000000 37.53457 0.000000 62.46543Table 5
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=0.75 (%).
rf x 1 x 3 x 4 x 8 x 16 x17 x18 x19 0.0033 37.73585 62.26415 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0093 2.467786 87.52520 0.000000 10.00701 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 54.51039 0.000000 42.69392 2.795687 0.000000 0.000000 0.000000 0.0438 0.000000 3.033993 0.000000 73.24945 23.71656 0.000000 0.000000 0.000000 0.0848 0.000000 0.000000 0.000000 6.626968 73.65104 19.72199 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 31.09411 9.147131 59.75876 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 15.95938 3.575565 80.46506 0.000000 0.0965 0.000000 0.000000 0.000000 0.000000 0.000000 88.80842 0.000000 11.19158 0.0970 0.000000 0.000000 0.000000 0.000000 0.000000 61.21745 0.000000 38.78255Table 6
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=1 (%).
rf x 1 x 3 x 4 x 8 x 16 x17 x18 x19 0.0033 37.73585 62.26415 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0093 0.960099 89.74003 0.000000 9.299870 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 46.90265 0.000000 53.09735 0.000000 0.000000 0.000000 0.000000 0.0438 0.000000 0.000000 0.000000 77.39837 22.60163 0.000000 0.000000 0.000000 0.0848 0.000000 0.000000 0.000000 0.000000 48.43750 51.56250 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 39.11838 7.310251 53.57137 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 23.37353 1.878349 74.74812 0.000000 0.0965 0.000000 0.000000 0.000000 0.153667 0.000000 99.84633 0.000000 0.000000 0.0970 0.000000 0.000000 0.000000 0.000000 0.000000 79.55936 0.000000 20.44064Next, we illustrate the solving process of DDM in detail with the caserf=0.0033 and λ=0. To begin with, we solve the equivalent parametric subproblems (35)–(37) with the fixed parameters values, respectively. By Lingo software, the objective value of subproblem (35) is 0.7501525×10-6 with solution (x1,x3,x4)=(37.08982%,51.43058%,11.47960%). In the meantime, after 21 solver iterations, the obtained objective value of subproblem (36) is 0.2788316×10-5 with solution (x1,x3,x8)=(39.85010%,59.15827%,0.9916381%). Moreover, the objective value and solution of subproblem (37) in this case are 0.1276898×10-4 and (x3,x10)=(81.86165%,17.00404%). Obviously, 0.7501525×10-6 is less than the other two objective values. According to DDM, the optimal solution of model (33) is the results of subproblem (35); that is, x∗=(x1,x3,x4)=(37.08982%,51.43058%,11.47960%). Under the different values of rf, how to find the optimal portfolio through comparing three submodels determined by DDM are shown in Table 7.Table 7
The obtained optimal portfolio results by comparing three submodels whenθl=0.9, θr=0.2, and λ=0.
rf 0.0033 0.0093 0.0193 0.0438 0.0848 0.0925 0.0940 0.0965 0.0970 Model (35) ∗ ∗ ∗ ∗ ∗ ∗ ∗ — — Model (36) √ √ √ √ √ √ √ ∗ — Model (37) √ √ √ √ — — — — — Notations:∗: optimal; √: feasible but nonoptimal; —: unfeasible.


### 7.2. Computational Result under Fuzzy Returns

For the sake of comparison, in this subsection, we solve model (33) with fuzzy returns. The fuzzy returns are characterized by the principle possibility distributions of interval-valued fuzzy variables, which are obtained by setting θlj=θrj=0 in parametric possibility distributions. Then, the returns of twenty securities are trapezoidal fuzzy variables ξjp=(r1j,r2j,r3j,r4j), j=1,2,…,20, and their expected values and semimoments are given in Table 8.Table 8
The values ofE[ξjp] and SM[ξjp] of twenty fuzzy returns.
Asset Fuzzy return(r1j,r2j,r3j,r4j) E [ ξ j p ] SM [ ξ j p ] 1 ξ 1 p = ( 0.9946,0.9967,1.0012,1.0061 ) 0.999650 8.18375 × 1 0 - 6 2 ξ 2 p = ( 0.9911,1.0002,1.0019,1.0052 ) 0.999600 1.12477 × 1 0 - 5 3 ξ 3 p = ( 0.9966,1.0073,1.0081,1.0109 ) 1.005725 1.18349 × 1 0 - 5 4 ξ 4 p = ( 0.9983,1.0090,1.0122,1.0135 ) 1.008250 1.53438 × 1 0 - 5 5 ξ 5 p = ( 1.0059,1.0122,1.0262,1.0279 ) 1.018050 4.21537 × 1 0 - 5 6 ξ 6 p = ( 1.0099,1.0159,1.0248,1.0429 ) 1.023375 5.63628 × 1 0 - 5 7 ξ 7 p = ( 1.0209,1.0225,1.0416,1.0553 ) 1.035075 8.95520 × 1 0 - 5 8 ξ 8 p = ( 1.0291,1.0299,1.0468,1.0678 ) 1.043400 9.66317 × 1 0 - 5 9 ξ 9 p = ( 1.0201,1.0504,1.0651,1.0709 ) 1.051625 1.72324 × 1 0 - 4 10 ξ 10 p = ( 1.0305,1.0399,1.0469,1.1093 ) 1.056650 2.33733 × 1 0 - 4 11 ξ 11 p = ( 1.0304,1.0468,1.0712,1.0908 ) 1.059800 2.35927 × 1 0 - 4 12 ξ 12 p = ( 1.0365,1.0529,1.0698,1.1096 ) 1.067200 2.64332 × 1 0 - 4 13 ξ 13 p = ( 1.0214,1.0699,1.0778,1.1004 ) 1.067375 3.33943 × 1 0 - 4 14 ξ 14 p = ( 1.0391,1.0699,1.0969,1.1096 ) 1.078875 3.36597 × 1 0 - 4 15 ξ 15 p = ( 1.0422,1.0766,1.0977,1.1178 ) 1.083575 3.41522 × 1 0 - 4 16 ξ 16 p = ( 1.0373,1.0914,1.0972,1.1171 ) 1.085750 3.50375 × 1 0 - 4 17 ξ 17 p = ( 1.0604,1.0786,1.1183,1.1282 ) 1.096375 3.74934 × 1 0 - 4 18 ξ 18 p = ( 1.0460,1.0952,1.1098,1.1361 ) 1.096775 4.43425 × 1 0 - 4 19 ξ 19 p = ( 1.0415,1.0885,1.1155,1.1475 ) 1.098250 6.44823 × 1 0 - 2 20 ξ 20 p = ( 1.0326,1.0896,1.1330,1.1392 ) 1.098600 8.38500 × 1 0 - 2Employing Lingo software, we solve model (33) and obtain the optimal allocation policy in Table 9. The results imply that the investor should put the money to different assets under different risk-free return levels.Table 9
The allocation proportions of model (33) with θl=0, θr=0 (%).
rf x 1 x 3 x 8 x 16 x 17 x 18 x 19 0.0033 37.73585 62.26415 0.000000 0.000000 0.000000 0.000000 0.000000 0.0093 2.077768 88.09815 9.824086 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 46.90265 53.09735 0.000000 0.000000 0.000000 0.000000 0.0438 0.000000 0.000000 77.39837 22.60163 0.000000 0.000000 0.000000 0.0848 0.000000 0.000000 6.882793 74.62438 18.49283 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 33.18711 8.668010 58.14488 0.000000 0.0940 0.000000 0.000000 0.000000 17.89324 3.132873 78.97389 0.000000 0.0965 0.000000 0.000000 0.000000 0.000000 93.33333 0.000000 6.666667 0.0970 0.000000 0.000000 0.000000 0.000000 66.66667 0.000000 33.33333


### 7.3. Parameter Analysis and Results Comparison

In this subsection, we will analyze the parameters in the fuzzy PSR portfolio selection model. In order to identify the parameter’s influence on the optimal portfolio selection, we solved the model with different values of parameters.Case 1 (the influence of parameterθ).
The parameterθ=(θl,θr) characterizes the degree of uncertainty of a parametric interval-valued fuzzy return. For example, regarding the parametric interval-valued fuzzy variable [1,2,3,4;θl,θr], the changeable range of parametric possibility distributions with different parameters θ=(θl,θr) is shown in Figure 1. From the figure, we know that the value of parameter θ is bigger; the uncertainty degree of fuzzy return is higher.Figure 1
The changeable range of parametric possibility distribution of parametric interval-valued fuzzy variable[1,2,3,4;θl,θr].
(a)
θ l = 0.8, θr=0.6 (b)
θ l = 0.3, θr=0.2We have computed the model with two sets of values:θl=0.9,θr=0.2, and θl=θr=0 in Tables 2 and 9. These results show that the optimal allocation policies are different with the variation of parameter θ. For example, if rf=0.0033, when θl=0.9, θr=0.2, and λ=0, the optimal allocation policy is selecting three assets with numbers 1, 3, and 4. In the meantime, if rf=0.0033, when θl=θr=0, the investor should choose two assets with numbers 1 and 3. The allocation proportions are distinctly different.Case 2 (the influence of parameterλ).
The parameterλ is a real number in [0,1], and it determines the location of variable possibility distribution. The λ selection is a variable representative of a parametric interval-valued fuzzy variable. Given the parameter θ, the variation of λ can make the possibility distribution run over the entire support of a parametric interval-valued fuzzy variable. Figure 2 shows the parametric possibility distributions of five selection variables of parametric interval-valued fuzzy variable [0,2,3,5;0.9,0.2] with parameter λ=0,0.25,0.5,0.75,1.Figure 2
The possibility distributions of[0,2,3,5;0.9,0.2] with five different values of λ.Assume parametersθl=0.9, θr=0.2; Tables 2–6 give the results of asset allocation with different λ. And we know that the selected assets of optimal portfolio vary in the candidate assets, and allocation proportions are different. From the results of the optimal value, we observe that the value of parameter λ is bigger; the portfolio risk is smaller. So for the proposed model we can reduce the portfolio risk by choosing the bigger value of parameter λ. Thus we draw the conclusion that the model is sensitive to parameter, and the portfolio risk is negative correlation with parameter λ.


### 7.4. Discussion

In this subsection, we will carry out a discussion about the difference between the variable possibility distribution and fixed possibility distribution. From the computational results in Sections7.1 and 7.2, we obtain the following findings.Firstly, when the fixed possibility distributions of uncertain returns are available, the portfolio selection model is a fuzzy optimization model. In this case, the optimal allocation policies are listed in Table9. From the computational results, we can observe that when the return level of the risk-free asset rf increases, the optimal portfolio proportions are different. For example, when rf is 0.0093, the investor should allocate his/her money to three assets 1, 3, and 8; when rf is 0.0848, the investor should select assets 8, 16, and 17. It is worth noting that if the decision maker cannot obtain the precise distribution information before making decisions, this method will be invalid.Secondly, when the fixed possibility distributions cannot be obtained during the portfolio selection process, we can characterize the uncertain returns with parametric interval-valued fuzzy distributions. In this situation, the portfolio selection model is a fuzzy parametric optimization model, in which two types of parameters are embedded. From the computational results in Tables2–6, we can see that the selected assets and allocation proportions change with the variation of parameter λ in [0,1]. Specifically, when the return level of the risk-free asset rf is 0.0193, if λ=0 and 0.25, the decision maker should allocate his/her money to four assets with numbers 3, 4, 16, and 17, while if λ=0.5 and 0.75, the money is allocated to three assets; and if λ=1, the selected assets are only two ones with numbers 3 and 8. Obviously, the optimal portfolio policy is closely related to parameter λ. Furthermore, comparing the variable possibility distribution method with the fixed possibility distribution method, we know the two methods can provide different investment while the return level rf varies from 0.0033 to 0.0970. Even though they provide diversified investment for the same assets, the allocation proportions to them are different.Finally, in other methods of interval fuzzy optimization, the fixed lower/upper bounds or midpoint of the distributional interval is often chosen as representative of the family of possibility distributions. In our parametric interval-valued method, the possibility distribution ofλ selection variable runs over the entire fuzzy support as the parameter λ varies in the unit interval [0,1]. Namely, the variable possibility distribution of selection variable can be any distribution in the fuzzy support. Thus the parametric optimization method based on the variable possibility distribution is flexible and comprehensive for decision makers. Particularly, if the investor cannot obtain the exact distributional information, the method of variable possibility distribution provides an effectual way for decision makers to make their informed portfolio selection policies.


## 7.1. Computational Result under Parametric Interval-Valued Fuzzy Returns

Firstly, when the exact possibility distributions of asset returns are unavailable, we conduct the experiments under variable possibility distribution to find the optimal portfolio selection policy.In our experiments, we set the values of parametersθli and θri, i=1,2,…,20 as (38)θl1,θl2,…,θl20=0.8,0.2,0.3,0.9,0.7,0.5,0.4,0.6,0.8,0.1,0.6,0.3,0.5,0.3,0.7,0.4,0.2,0.8,0.5,0.2,θr1,θr2,…,θr20=0.3,0.8,0.2,0.6,0.3,0.5,0.3,0.9,0.7,0.5,0.4,0.6,0.7,0.4,0.9,0.8,0.5,0.2,0.8,0.2. Since θl=max1≤i≤nθli and θr=min1≤i≤nθri, we have θl=0.9, θr=0.2.Assume the investor selects the values ofλ as 0, 0.25, 0.5, 0.75, and 1. To test the effect of the risk-free return rate on the portfolio decision, we adjust the value of rf and do some experiments, and the returns of the risk-free asset are set as follows rf= 0.33%, 0.93%, 1.93%, 4.38%, 8.48%, 9.25%, 9.40%, 9.65%, and 9.70%. The obtained optimal allocation proportions to twenty risk assets are reported in Tables 2–6. Note that, according to x∗=y∗/t∗, we have the allocation proportions xi easily by reducing yi with the same ratio t. For simplicity, we give directly the results of xi.Table 2
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=0 (%).
r f x 1 x 3 x 4 x 8 x 16 x 17 x 18 x 19 0.0033 37.08982 51.43058 11.47960 0.000000 0.000000 0.000000 0.000000 0.000000 0.0093 0.000000 53.50372 37.74321 8.753070 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 47.97375 36.58445 0.000000 2.875511 12.56628 0.000000 0.000000 0.0438 0.000000 27.74888 26.18283 0.000000 25.06405 21.00424 0.000000 0.000000 0.0848 0.000000 0.000000 22.61346 0.000000 62.15479 35.32249 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 11.84791 13.55289 74.59920 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 0.000000 3.144038 86.73522 10.12074 0.0965 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.145002 97.85500 0.0970 — — — — — — — —Table 3
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=0.25 (%).
r f x 1 x 3 x 4 x 8 x 16 x 17 x 18 x 19 0.0033 38.99274 60.41775 0.000000 0.589515 0.000000 0.000000 0.000000 0.000000 0.0093 5.941417 61.16051 22.99126 9.906809 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 50.47323 34.22256 0.000000 3.955556 11.34866 0.000000 0.000000 0.0438 0.000000 53.58670 0.000000 0.000000 26.61964 19.79366 0.000000 0.000000 0.0848 0.000000 2.979972 0.000000 0.000000 65.03688 31.98315 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 17.58623 12.23930 70.17447 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 3.478543 6.432623 90.08883 0.000000 0.0965 0.000000 0.000000 0.000000 0.000000 0.000000 28.36709 7.005321 64.62759 0.0970 0.000000 0.000000 0.000000 0.000000 0.000000 5.780347 0.000000 94.21965Table 4
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=0.5 (%).
rf x 1 x 3 x 4 x 8 x 16 x17 x18 x19 0.0033 38.02403 61.84081 0.000000 0.135164 0.000000 0.000000 0.000000 0.000000 0.0093 3.779654 85.59803 0.000000 10.62232 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 84.02355 0.000000 0.000000 4.756157 11.22030 0.000000 0.000000 0.0438 0.000000 53.91651 0.000000 0.000000 28.45683 17.62665 0.000000 0.000000 0.0848 0.000000 3.533319 0.000000 0.000000 68.11919 28.34749 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 23.96391 10.77935 65.25674 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 9.371306 5.083677 85.54502 0.000000 0.0965 0.000000 0.000000 0.000000 0.000000 0.000000 69.14263 0.000000 30.85737 0.0970 0.000000 0.000000 0.000000 0.000000 0.000000 37.53457 0.000000 62.46543Table 5
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=0.75 (%).
rf x 1 x 3 x 4 x 8 x 16 x17 x18 x19 0.0033 37.73585 62.26415 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0093 2.467786 87.52520 0.000000 10.00701 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 54.51039 0.000000 42.69392 2.795687 0.000000 0.000000 0.000000 0.0438 0.000000 3.033993 0.000000 73.24945 23.71656 0.000000 0.000000 0.000000 0.0848 0.000000 0.000000 0.000000 6.626968 73.65104 19.72199 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 31.09411 9.147131 59.75876 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 15.95938 3.575565 80.46506 0.000000 0.0965 0.000000 0.000000 0.000000 0.000000 0.000000 88.80842 0.000000 11.19158 0.0970 0.000000 0.000000 0.000000 0.000000 0.000000 61.21745 0.000000 38.78255Table 6
The allocation proportions of model (33) with θl=0.9, θr=0.2, and λ=1 (%).
rf x 1 x 3 x 4 x 8 x 16 x17 x18 x19 0.0033 37.73585 62.26415 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0093 0.960099 89.74003 0.000000 9.299870 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 46.90265 0.000000 53.09735 0.000000 0.000000 0.000000 0.000000 0.0438 0.000000 0.000000 0.000000 77.39837 22.60163 0.000000 0.000000 0.000000 0.0848 0.000000 0.000000 0.000000 0.000000 48.43750 51.56250 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 0.000000 39.11838 7.310251 53.57137 0.000000 0.0940 0.000000 0.000000 0.000000 0.000000 23.37353 1.878349 74.74812 0.000000 0.0965 0.000000 0.000000 0.000000 0.153667 0.000000 99.84633 0.000000 0.000000 0.0970 0.000000 0.000000 0.000000 0.000000 0.000000 79.55936 0.000000 20.44064Next, we illustrate the solving process of DDM in detail with the caserf=0.0033 and λ=0. To begin with, we solve the equivalent parametric subproblems (35)–(37) with the fixed parameters values, respectively. By Lingo software, the objective value of subproblem (35) is 0.7501525×10-6 with solution (x1,x3,x4)=(37.08982%,51.43058%,11.47960%). In the meantime, after 21 solver iterations, the obtained objective value of subproblem (36) is 0.2788316×10-5 with solution (x1,x3,x8)=(39.85010%,59.15827%,0.9916381%). Moreover, the objective value and solution of subproblem (37) in this case are 0.1276898×10-4 and (x3,x10)=(81.86165%,17.00404%). Obviously, 0.7501525×10-6 is less than the other two objective values. According to DDM, the optimal solution of model (33) is the results of subproblem (35); that is, x∗=(x1,x3,x4)=(37.08982%,51.43058%,11.47960%). Under the different values of rf, how to find the optimal portfolio through comparing three submodels determined by DDM are shown in Table 7.Table 7
The obtained optimal portfolio results by comparing three submodels whenθl=0.9, θr=0.2, and λ=0.
rf 0.0033 0.0093 0.0193 0.0438 0.0848 0.0925 0.0940 0.0965 0.0970 Model (35) ∗ ∗ ∗ ∗ ∗ ∗ ∗ — — Model (36) √ √ √ √ √ √ √ ∗ — Model (37) √ √ √ √ — — — — — Notations:∗: optimal; √: feasible but nonoptimal; —: unfeasible.


## 7.2. Computational Result under Fuzzy Returns

For the sake of comparison, in this subsection, we solve model (33) with fuzzy returns. The fuzzy returns are characterized by the principle possibility distributions of interval-valued fuzzy variables, which are obtained by setting θlj=θrj=0 in parametric possibility distributions. Then, the returns of twenty securities are trapezoidal fuzzy variables ξjp=(r1j,r2j,r3j,r4j), j=1,2,…,20, and their expected values and semimoments are given in Table 8.Table 8
The values ofE[ξjp] and SM[ξjp] of twenty fuzzy returns.
Asset Fuzzy return(r1j,r2j,r3j,r4j) E [ ξ j p ] SM [ ξ j p ] 1 ξ 1 p = ( 0.9946,0.9967,1.0012,1.0061 ) 0.999650 8.18375 × 1 0 - 6 2 ξ 2 p = ( 0.9911,1.0002,1.0019,1.0052 ) 0.999600 1.12477 × 1 0 - 5 3 ξ 3 p = ( 0.9966,1.0073,1.0081,1.0109 ) 1.005725 1.18349 × 1 0 - 5 4 ξ 4 p = ( 0.9983,1.0090,1.0122,1.0135 ) 1.008250 1.53438 × 1 0 - 5 5 ξ 5 p = ( 1.0059,1.0122,1.0262,1.0279 ) 1.018050 4.21537 × 1 0 - 5 6 ξ 6 p = ( 1.0099,1.0159,1.0248,1.0429 ) 1.023375 5.63628 × 1 0 - 5 7 ξ 7 p = ( 1.0209,1.0225,1.0416,1.0553 ) 1.035075 8.95520 × 1 0 - 5 8 ξ 8 p = ( 1.0291,1.0299,1.0468,1.0678 ) 1.043400 9.66317 × 1 0 - 5 9 ξ 9 p = ( 1.0201,1.0504,1.0651,1.0709 ) 1.051625 1.72324 × 1 0 - 4 10 ξ 10 p = ( 1.0305,1.0399,1.0469,1.1093 ) 1.056650 2.33733 × 1 0 - 4 11 ξ 11 p = ( 1.0304,1.0468,1.0712,1.0908 ) 1.059800 2.35927 × 1 0 - 4 12 ξ 12 p = ( 1.0365,1.0529,1.0698,1.1096 ) 1.067200 2.64332 × 1 0 - 4 13 ξ 13 p = ( 1.0214,1.0699,1.0778,1.1004 ) 1.067375 3.33943 × 1 0 - 4 14 ξ 14 p = ( 1.0391,1.0699,1.0969,1.1096 ) 1.078875 3.36597 × 1 0 - 4 15 ξ 15 p = ( 1.0422,1.0766,1.0977,1.1178 ) 1.083575 3.41522 × 1 0 - 4 16 ξ 16 p = ( 1.0373,1.0914,1.0972,1.1171 ) 1.085750 3.50375 × 1 0 - 4 17 ξ 17 p = ( 1.0604,1.0786,1.1183,1.1282 ) 1.096375 3.74934 × 1 0 - 4 18 ξ 18 p = ( 1.0460,1.0952,1.1098,1.1361 ) 1.096775 4.43425 × 1 0 - 4 19 ξ 19 p = ( 1.0415,1.0885,1.1155,1.1475 ) 1.098250 6.44823 × 1 0 - 2 20 ξ 20 p = ( 1.0326,1.0896,1.1330,1.1392 ) 1.098600 8.38500 × 1 0 - 2Employing Lingo software, we solve model (33) and obtain the optimal allocation policy in Table 9. The results imply that the investor should put the money to different assets under different risk-free return levels.Table 9
The allocation proportions of model (33) with θl=0, θr=0 (%).
rf x 1 x 3 x 8 x 16 x 17 x 18 x 19 0.0033 37.73585 62.26415 0.000000 0.000000 0.000000 0.000000 0.000000 0.0093 2.077768 88.09815 9.824086 0.000000 0.000000 0.000000 0.000000 0.0193 0.000000 46.90265 53.09735 0.000000 0.000000 0.000000 0.000000 0.0438 0.000000 0.000000 77.39837 22.60163 0.000000 0.000000 0.000000 0.0848 0.000000 0.000000 6.882793 74.62438 18.49283 0.000000 0.000000 0.0925 0.000000 0.000000 0.000000 33.18711 8.668010 58.14488 0.000000 0.0940 0.000000 0.000000 0.000000 17.89324 3.132873 78.97389 0.000000 0.0965 0.000000 0.000000 0.000000 0.000000 93.33333 0.000000 6.666667 0.0970 0.000000 0.000000 0.000000 0.000000 66.66667 0.000000 33.33333


## 7.3. Parameter Analysis and Results Comparison

In this subsection, we will analyze the parameters in the fuzzy PSR portfolio selection model. In order to identify the parameter’s influence on the optimal portfolio selection, we solved the model with different values of parameters.Case 1 (the influence of parameterθ).
The parameterθ=(θl,θr) characterizes the degree of uncertainty of a parametric interval-valued fuzzy return. For example, regarding the parametric interval-valued fuzzy variable [1,2,3,4;θl,θr], the changeable range of parametric possibility distributions with different parameters θ=(θl,θr) is shown in Figure 1. From the figure, we know that the value of parameter θ is bigger; the uncertainty degree of fuzzy return is higher.Figure 1
The changeable range of parametric possibility distribution of parametric interval-valued fuzzy variable[1,2,3,4;θl,θr].
(a)
θ l = 0.8, θr=0.6 (b)
θ l = 0.3, θr=0.2We have computed the model with two sets of values:θl=0.9,θr=0.2, and θl=θr=0 in Tables 2 and 9. These results show that the optimal allocation policies are different with the variation of parameter θ. For example, if rf=0.0033, when θl=0.9, θr=0.2, and λ=0, the optimal allocation policy is selecting three assets with numbers 1, 3, and 4. In the meantime, if rf=0.0033, when θl=θr=0, the investor should choose two assets with numbers 1 and 3. The allocation proportions are distinctly different.Case 2 (the influence of parameterλ).
The parameterλ is a real number in [0,1], and it determines the location of variable possibility distribution. The λ selection is a variable representative of a parametric interval-valued fuzzy variable. Given the parameter θ, the variation of λ can make the possibility distribution run over the entire support of a parametric interval-valued fuzzy variable. Figure 2 shows the parametric possibility distributions of five selection variables of parametric interval-valued fuzzy variable [0,2,3,5;0.9,0.2] with parameter λ=0,0.25,0.5,0.75,1.Figure 2
The possibility distributions of[0,2,3,5;0.9,0.2] with five different values of λ.Assume parametersθl=0.9, θr=0.2; Tables 2–6 give the results of asset allocation with different λ. And we know that the selected assets of optimal portfolio vary in the candidate assets, and allocation proportions are different. From the results of the optimal value, we observe that the value of parameter λ is bigger; the portfolio risk is smaller. So for the proposed model we can reduce the portfolio risk by choosing the bigger value of parameter λ. Thus we draw the conclusion that the model is sensitive to parameter, and the portfolio risk is negative correlation with parameter λ.


## 7.4. Discussion

In this subsection, we will carry out a discussion about the difference between the variable possibility distribution and fixed possibility distribution. From the computational results in Sections7.1 and 7.2, we obtain the following findings.Firstly, when the fixed possibility distributions of uncertain returns are available, the portfolio selection model is a fuzzy optimization model. In this case, the optimal allocation policies are listed in Table9. From the computational results, we can observe that when the return level of the risk-free asset rf increases, the optimal portfolio proportions are different. For example, when rf is 0.0093, the investor should allocate his/her money to three assets 1, 3, and 8; when rf is 0.0848, the investor should select assets 8, 16, and 17. It is worth noting that if the decision maker cannot obtain the precise distribution information before making decisions, this method will be invalid.Secondly, when the fixed possibility distributions cannot be obtained during the portfolio selection process, we can characterize the uncertain returns with parametric interval-valued fuzzy distributions. In this situation, the portfolio selection model is a fuzzy parametric optimization model, in which two types of parameters are embedded. From the computational results in Tables2–6, we can see that the selected assets and allocation proportions change with the variation of parameter λ in [0,1]. Specifically, when the return level of the risk-free asset rf is 0.0193, if λ=0 and 0.25, the decision maker should allocate his/her money to four assets with numbers 3, 4, 16, and 17, while if λ=0.5 and 0.75, the money is allocated to three assets; and if λ=1, the selected assets are only two ones with numbers 3 and 8. Obviously, the optimal portfolio policy is closely related to parameter λ. Furthermore, comparing the variable possibility distribution method with the fixed possibility distribution method, we know the two methods can provide different investment while the return level rf varies from 0.0033 to 0.0970. Even though they provide diversified investment for the same assets, the allocation proportions to them are different.Finally, in other methods of interval fuzzy optimization, the fixed lower/upper bounds or midpoint of the distributional interval is often chosen as representative of the family of possibility distributions. In our parametric interval-valued method, the possibility distribution ofλ selection variable runs over the entire fuzzy support as the parameter λ varies in the unit interval [0,1]. Namely, the variable possibility distribution of selection variable can be any distribution in the fuzzy support. Thus the parametric optimization method based on the variable possibility distribution is flexible and comprehensive for decision makers. Particularly, if the investor cannot obtain the exact distributional information, the method of variable possibility distribution provides an effectual way for decision makers to make their informed portfolio selection policies.


## 8. Conclusions

In portfolio management, an important role of performance measure is to help investors evaluate their portfolios and make a portfolio selection decision. This paper presents a theoretically sound portfolio performance to optimize the investment policy, and the main conclusions are summarized as follows:(i)
A new downside risk measure: the semimoment of the selection variable is presented to characterize the loss risk. Using the proposed measure, we modify the conventional Sharpe ratio and define a parametric Sharpe ratio for fuzzy portfolio performance evaluation.(ii)
A novel fuzzy PSR portfolio selection model is built wherein the uncertain returns are characterized by interval-valued possibility distributions. When the exact distributions of uncertain returns are difficult to be determined in advance, the proposed approach could be employed to characterize the uncertainty in practice.(iii)
By performing an equivalent transformation, the proposed model is transformed into nonfractional programming. Furthermore, according to structure characteristics, a domain decomposition method is introduced and the proposed model is solved by computing three equivalent parametric submodels.(iv)
We conduct the experiments and compare the computational results under variable possibility distribution with fixed possibility distribution and analyze the influence of different parameters on the optimal portfolio.Future academical innovation directions might address the following two topics. On the one hand, the current risk measure used in this paper is semimoment; an extension to other downside risk measure is possible. On the other hand, the present model finds the optimal portfolio decision via parametric optimization method; the case dealing with distributionally robust optimization method based on the uncertain distribution set is a future avenue for our research.



---

*Source: 6279859-2017-01-31.xml*