After definition of the discrete grey stochastic variable and its expected value, the expected probability degree is defined. For multicriteria decisionmaking problems, in which the criteria weights are incompletely certain and the criteria values of alternatives are in the form of grey stochastic variables, a grey stochastic multicriteria decisionmaking approach is proposed. In this method, the evaluation value of each alternative under each criterion can be transformed to comprise the expected probability degree judgment matrix, based on which, a nonlinear programming model can be enacted. In the end, the genetic algorithm is used to solve the model to attain the criteria weights, and the ranking of alternatives can be produced consequently. The feasibility and validity of this approach are illustrated by an example.
Grey stochastic variable ; Expected probability degree ; Alternative similarity scale ; Multicriteria decisionmaking ; Genetic algorithm
Multicriteria decisionmaking theories and methods have become one of the most active subjects in many disciplines, such as decisionmaking science, system engineering, management and logistics, etc. Since the rapid development of society and economics, the fuzziness of the decisionmaking environment has been realized by more and more people, in addition to the complexity and uncertainty existing commonly in decisionmaking problems. Hence, information should be presented by fuzziness, randomness and uncertainty in a real decisionmaking process. So far, there has been considerable research unto multicriteria decisionmaking problems, which are in the form of the above three types of uncertain decisionmaking information, and, meanwhile, some studies have also focused on problems with multiple kinds of uncertainty. The multicriteria decisionmaking problems, in which criteria values have randomness and fuzziness simultaneously are discussed in [1] , [2] and [3] . Problems with criterion values taking the form of grey and fuzziness in the meantime are discussed in [4] , [5] and [6] . However, there has been relatively little research conducted on multicriteria decisionmaking problems with criteria value in the form of grey and randomness at the same time. In [7] , the complementary problems of grey theory and random theory are studied, showing that some methods for random problems contain the thoughts of grey methods and concepts, while grey problems can be also realized and solved from the randomness perspective. The objective function can be established by minimization of the comprehensive weighed distances of alternatives, which are to the positive alternative and the negative alternative. To solve that function, a method is provided to deal with some multicriteria decisionmaking problems, in which the criterion value is in the form of a random variable. However, it only takes into consideration the determining of the criteria value in [8] . As for the grey risk decisionmaking problems with uncertain criteria weights, the twobasepoint method is proposed by applying the thought of a positive and negative ideal point to the grey stochastic domain in [9] . Considering the dynamic hybrid multiattribute decision making problems under risk, in which the weights are uncertain and criteria values take the form of precise numbers and interval numbers, as well as languagetype fuzzy numbers at the same time, an approach, based on grey matrix relative degree, is provided to solve them in [10] . The main thought of this method is transforming the risk judgment matrix to a riskfree judgment matrix, and then ranking the alternatives using the grey interval relative analysis method. According to the characteristics of multicriteria decision making problems under risk, a reasonable solution is given to these problems by introducing the concept of probability preference in [11] . Even though the abovementioned research points out the problems existing in grey stochastic, there are some kinds of preliminary studies on these issues. Since there is little research on grey stochastic multicriteria decisionmaking problems, and these kinds of problems can be more objective in describing some decisionmaking situations in reality, they are worthy of more attention. Thus, in this paper, a method is given to meet the demand of practical decisionmaking.
The rest of this paper is organized as follows. Section 2 gives the definition of a grey number, grey stochastic variable and expected probability degree. Section 3 proposes the grey stochastic multicriteria decisionmaking method, based on expected probability degree, and illustrates the procedures in detail. Section 4 applies the proposed method to a practical example to explain its rationality and effectiveness. Section 5 is the conclusion.
The grey number can be defined as the number with a general range, but the exact value of this number cannot be known. In application, the grey number is an uncertain number which takes the value in a scope or a particular number set. It can be denoted as .
An interval grey number can be defined as the grey number with the lower limit, , and the upper limit, . The interval grey number can be denoted as .
Assume interval grey numbers, . The interval grey numbers’ operational rules can be defined as follows [12] :
If the stochastic variables are countable values in the form of interval grey numbers, and the corresponding probability with regard to each value can be attained, then, this kind of variable can be defined as the discrete grey stochastic variable, which is called a grey stochastic variable, for short, in this paper.
The grey stochastic variable is denoted as , and its th value can be presented by . Table 1 shows the probability distribution of .
ζ(⊗)  ⊗1  ⊗2  …  ⊗i  …  ⊗n 
P  P1  P2  …  Pi  …  Pn 
In Table 1 , is a grey stochastic variable. is the th possible value that would be taken by is the probability with respect to , while is the number of values that a grey stochastic variable can have.
The probability density function can be denoted as

For better understanding of the above definitions, an example is given. Assume that an investment alternative, , may gain an annual revenue of 270–280 million RMB with the probability of 0.4, 290–300 million with 0.2 possibility and 310–340 million with 0.4. This information can be presented by the grey stochastic variable, , in Table 2 .
ζ(⊗)  [2.7, 2.8]  [2.9, 3.0]  [3.1, 3.4] 
P  0.4  0.2  0.4 
Meanwhile, the probability density function can be denoted as follows:

Assume that is a grey stochastic variable, then, the expected value of the grey stochastic variable is defined by , if the value of the formula can be attained. And the expected value can be denoted as , which satisfies:

According to the grey stochastic variable’s operational rules (1) and (4) , it can be concluded that is an interval grey number.
The concept of the expected value of the grey stochastic variable can be also illustrated with the abovementioned example.

Let be the length of the grey number, , and be the length of the grey number, , then,

is called the expected probability degree of against . Therefore, the relationship between and can be determined as follows:
Assume that and are grey stochastic variables, and their probability density functions can be denoted as and , respectively. The expected probability degree of grey stochastic variable, , against , can be denoted as :

where and are the values that the variables and may take separately. , and . are the possibility degrees of the interval grey number, , which is more than or equal to .
The expected probability degree is the measurement of the average probability degree in which the grey stochastic variable, , is better than [11] . The properties of expected probability degree can be listed as follows:
The proof procedure for these properties will be shown as follows.
As ,

and .
Furthermore, considering any value that or may take, there always exist and , so, one can reach the conclusion that .
Due to , as well as , we have:

Since (see Ref. [9] ), then:

According to the property (2) , if , then:

Therefore, .
Assume that is a discrete set of grey stochastic variables, then the expected probability degree judgment matrix can be defined as a matrix that consists of all the expected probability degrees in which each grey stochastic variable is better than other variables separately. It can be denoted as , where: .
The expected probability degree judgment matrix is a complementary judgment matrix.
According to the definition for the expected probability degree, it is easy to reach the conclusion that in the expected probability degree judgment matrix. Furthermore, there always exists the equality , so, the judgment matrix of the expected probability degree is a complementary judgment matrix.
There is a grey stochastic multicriteria decisionmaking problem. Assume that is a discrete alternative set of possible alternatives, and is a set of criteria. is the vector of criteria weights, and the weights are subject to the constraints as follows:

The criterion value of alternative , with respect to criterion , is denoted as , and is a grey stochastic variable. The incompletely certain information of the criteria weights can be represented by symbol and it can be also shown in the form of linear equalities and inequalities, as follows [14] :
The order of the alternatives is needed to be listed under the above conditions.
For the mentioned grey stochastic multicriteria decisionmaking problems, the solving procedure can be summarized as follows:
Step 1. Establish the expected probability degree judgment matrix.
Calculate the expected probability degree of every alternative, with respect to each criterion, and form the expected probability degree judgment matrix. As far as criterion is concerned, the expected probability degree judgment matrix for the alternative set can be denoted as :

where .
According to Theorem 1 , is a complementary judgment matrix.
Step 2. Form the comprehensive judgment matrix of expected probability degree.
After gaining the expected probability degree judgment matrix, with respect to each criterion, the comprehensive judgment matrix of expected probability degree can be denoted by , where .
It can be proved that the comprehensive judgment matrix of the expected probability degree is also a complementary judgment matrix, and the process proof can be shown as follows.
According to property (3) of the expected probability degree, there always exists , so does the expression . Besides,

thus, matrix is a complementary judgment matrix.
Step 3. Calculate the criteria weights.
According to the sorting vector, was introduced in [15] for solving the ranking problem of the complementary judgment matrix. can be calculated by the method as follows:

( 1) 
where .
If the criteria weights are already known, the order of the alternative set can be attained by comparing . However, the situation of criteria weights, which are incompletely certain, is common in real multicriteria decisionmaking. Consequently, by drawing lessons from Refs. [16] , [17] and [18] , establishing an optimization model based on the closeness degree of is a major way of obtaining the criteria weights in this paper.
For Formula (1) , it can be transformed into Formula (2) , as follows:

( 2) 
where .
Therefore, the comprehensive sorting value of alternative can be represented by . Furthermore, the comprehensive sorting vector of the alternative set can be denoted by . As for given , obviously, if the value of is bigger, the corresponding alternative, , would be better. Consequently, the multiobjective decisionmaking model can be built as follows:

( 3) 
Solve the programming model as follows:

( 4) 
Since model (4) is a linear programming model, it can be solved to obtain the optimal weight vector , with respect to alternative , using the simplex method. Then, the positive comprehensive sorting value of alternative can be attained in terms of the , namely, .
Then, we solve the programming model as follows:

( 5) 
Since model (5) is also a linear programming model, it can be solved to obtain the optimal weight vector, , with respect to alternative , using the simplex method. Then, the negative comprehensive sorting value of alternative can be attained in terms of , namely, .
Vectors and are called the positive and negative ideal points, respectively, in this multicriteria decisionmaking problem.
Consequently, the closeness degree function for the alternative is defined by the following expression, as:

( 6) 
For any weight vector, , which belongs to , if the value that can have is bigger, then is closer to , which means the alternatives are closer to the optimal status from the whole view, and vice versa. Therefore, the programming model can be built as follows:

( 7) 
As model (7) is a nonlinear programming model, and the objective function is complicated, it is difficult to solve it by common ways. In this paper, the genetic algorithm [19] is introduced to deal with the mentioned model, consequently, the optimal criteria weight vector can be calculated and denoted as .
Step 4. Rank the alternatives.
According to and Formula (2) , the order of the alternatives can be obtained.
An investment bank is planning to invest in three listed companies, which are denoted as and , accordingly. There are three criteria taken into account, namely, annual product income, , social benefit, , and environmental pollution degree, . Among these three criteria, and belong to the benefit type of criterion, while is a cost criterion. All three companies would have three possible values, which are in the form of a grey interval number under each criterion, and the corresponding probabilities are known. The vector of criteria weights are denoted with . The ranking of the alternatives need to be given under the abovementioned conditions. The data for the alternatives are shown in Table 3 . In Table 3 , “ ” is the abbreviation of “alternative”, “ ” is “possible value”, “ ” is “Probability”, and “ ” is “Criterion”.
C1  0.4  [2.7, 2.7]  [2.5, 2.5]  [3.1, 3.1] 
0.2  [3.0, 3.0]  [2.1, 2.1]  [3.5, 3.5]  
0.4  [2.8, 2.8]  [2.7, 2.7]  [2.9, 2.9]  
C2  0.4  [3.5, 4.0]  [3.5, 3.9]  [3.3, 3.5] 
0.2  [3.9, 4.4]  [4.4, 4.5]  [2.6, 3.1]  
0.4  [3.3, 3.8]  [3.8, 4.1]  [3.2, 3.7]  
C3  0.4  [0.25, 0.4]  [0.4, 0.6]  [0.4, 0.6] 
0.2  [0.1, 0.25]  [0.6, 0.75]  [0.25, 0.4]  
0.4  [0.4, 0.6]  [0.25, 0.4]  [0.6, 0.75] 
The steps to solve the above problem can be displayed as follows:
Step 1. Calculate the expected probability degree of the alternative, with regard to each criterion, and establish the expected probability degree judgment matrix.
In Table 2 , and are benefit criteria, but is a cost criterion. Consequently, should be transformed to the benefit criterion for the unification of evaluation. Table 4 shows the results.
C1  0.4  [2.7, 2.7]  [2.5, 2.5]  [3.1, 3.1] 
0.2  [3.0, 3.0]  [2.1, 2.1]  [3.5, 3.5]  
0.4  [2.8, 2.8]  [2.7, 2.7]  [2.9, 2.9]  
C2  0.4  [3.5, 4.0]  [3.5, 3.9]  [3.3, 3.5] 
0.2  [3.9, 4.4]  [4.4, 4.5]  [2.6, 3.1]  
0.4  [3.3, 3.8]  [3.8, 4.1]  [3.2, 3.7]  
C3  0.4  [−0.4, −0.25]  [−0.6, −0.4]  [−0.6, −0.4] 
0.2  [−0.25, −0.1]  [−0.75, −0.6]  [−0.4, −0.25]  
0.4  [−0.6, −0.4]  [−0.4, −0.25]  [−0.75, −0.6] 
The expected probability degree judgment matrix, with regard to each criterion, can be attained conveniently in Matlab 7.0, and the matrices can be denoted as , and . represents the judgment matrix with regard to with regard to , and with regard to .

Step 2. Form the comprehensive judgment matrix of the expected probability degree.
Since the criteria weights are incompletely certain, the comprehensive judgment matrix of the expected probability degree can be denoted by matrix as given in Box I :

Step 3. Calculate the criteria weights.
According to Formula (2) , the comprehensive sorting value of alternative can be shown as follows:

Consequently, the programming model for alternative , which maximizes the value of , can be established as follows:

Using the simplex method to solve the above models, the results can be shown as follows:

And the programming model for alternative , which minimizes the value of , can be established as follows:

Consequently, .
Therefore, the nonlinear programming model can be built as follows:

The optimal criteria weights can be easily obtained using the genetic algorithm toolbox in Matlab 7.0, and they are .
Step 4: Rank the alternatives.
According to and Formula (2) , we calculate the sorting vector, which is . Thus, the order of alternatives is , which means should be the priority for the investment bank, while is the second choice, and is the last one. If we use the approach proposed in [20] , the same ranking order of the alternatives is obtained. Thus, the proposed approach is valid.
An approach is proposed in this paper for grey stochastic multicriteria decisionmaking problems with incompletely uncertain criteria weights. Firstly, the definition of the expected probability degree of a grey stochastic variable is given by introducing the concept of probability preference, and its properties are discussed in this paper. Then, the evaluations of the alternatives, with respect to the criteria, can be transformed into the expected probability degree judgment matrices. The comprehensive judgment matrix was consequently formed. As the criteria weights are incompletely certain, an optimal programming model, which is based on the closeness degree of the sorting vector, is built, and solved by the genetic algorithm to obtain optimal criteria weights. Therefore, the order of the alternative is listed. This approach is verified by an example.
The authors thank the editors and anonymous reviewers for their helpful comments and suggestions. This work was supported by the National Natural Science Foundation of China (Nos. 71271218 and 71221061 ) and the Layout Foundation of Humanities and Social Sciences of PR.C Ministry of Education (No. 11YJCZH227 ).
Published on 06/10/16
Licence: Other
Are you one of the authors of this document?