Abstract

In order to determine the membership of an element to a set owing to ambiguity between a few different values, the hesitant fuzzy set (HFS) has been proposed and widely diffused to deal with vagueness and uncertainty involved in the process of multiple criteria group decision making (MCGDM) problems. In this paper, we develop novel definitions of score function and distance measure for HFSs. Some examples are given to illustrate that the proposed definitions are more reasonable than the traditional ones. Furthermore, our study extends the MULTIMOORA (Multiple Objective Optimization on the basis of Ratio Analysis plus Full Multiplicative Form) method with HFSs. The proposed method thus provides the means for multiple criteria decision making (MCDM) regarding uncertain assessments. Utilization of hesitant fuzzy power aggregation operators also enables facilitating the process of MCGDM. A numerical example of software selection demonstrates the possibilities of application of the proposed method.

1. Introduction

The fuzzy set (FS) theory, which is a generalization of classical set theory introduced by Zadeh [1], has drawn the attention of many researchers. Several classical extensions such as the interval-valued fuzzy set (IVFS) [2], intuitionistic fuzzy set (IFS) [36], interval-valued intuitionistic fuzzy set (IVIFS) [7], linguistic fuzzy set (IFS) [8, 9], type-2 fuzzy set (T-2FS) [1013], type- fuzzy set [14], and fuzzy multiset [15] have been developed. In general, many scholars have noticed that the difficulty in establishing the degree of membership of an element in a set does not arise from a margin of error (as in IFS or IVFSs) or a specified possibility distribution of the possible values (as in T-2FSs) but instead originates from our hesitation between a few different values [16]. In fact, it is usually difficult to reach a final agreement during the decision making process when people are hesitant and irresolute for one thing or another. For example, two decision makers (DMs) may differ on the membership degree of an element to a set to be assigned, one perhaps wants to assign 0.5 but the other 0.6 and they cannot persuade each other. As such, there is a set of possible values when describing the membership degree of an element. In such cases, the concept of a hesitant fuzzy set (HFS) naturally introduced by Torra and Narukawa [17] and Torra [18] was aimed at dealing with the inherent hesitation in humankind when making a decision. The membership degree of an element to a HFS, called the hesitant fuzzy element (HFE), can be represented by a collection of possible values between 0 and 1 [19]. In this sense, it can provide a relative accurate representation of people’s hesitancy when they are describing their preferences over objects, whereas the fuzzy set or its classical extensions cannot handle this situation well. Recently, a great number of other extensions of the HFSs have been developed such as dual hesitant fuzzy sets (DHFSs) [20], interval-valued hesitant fuzzy sets (IVHFS) [21], generalized hesitant fuzzy sets (GHFSs) [22, 23], and hesitant fuzzy linguistic term sets (HFLTSs) [8, 24]. It is worth noticing that the HFS theory has been applied to many practical problems, especially in the area of decision making [9, 2128].

The main task of multiple criteria group decision making (MCGDM) is to find the most desirable alternative(s) among a set of feasible alternatives based on the preferences provided by a group of decision makers or experts [27, 28]. In recent years, MCGDM with HFSs has attracted a great deal of attention by experts in decision making area [9, 23, 2628]. Xia and Xu [25] proposed several aggregation operators for hesitant fuzzy information and discussed the relationship between IFS and HFS, based on which they developed some operations and aggregation operators for HFEs. They also illustrated the correlations among the aggregation operators and gave their application in solving decision making problems. Zhang [16] developed a series of hesitant fuzzy power aggregation operators for hesitant fuzzy information, demonstrated several useful properties of the operators and discussed the relationships between them. The new aggregation operators are utilized to develop techniques for MCGDM with hesitant fuzzy information. To differently process and effectively aggregate hesitant fuzzy information and capture their interrelationship Zhou [19] proposed the hesitant fuzzy reducible weighted Bonferroni mean (HFRWBM) and studied some of its prominent characteristics. They further investigated its generalized form, that is, the generalized hesitant fuzzy reducible weighted Bonferroni mean (GHFRWBM). They finally applied the proposed aggregation operators to multicriteria aggregation. Yu et al. [26] investigated aggregation methods for prioritized hesitant fuzzy elements and their application on personnel evaluation. The generalized hesitant fuzzy prioritized weighted average (GHFPWA) and generalized hesitant fuzzy prioritized weighted geometric (GHFPWG) operators are proposed and some desirable properties of the methods are discussed. A procedure and algorithm for MCGDM was provided using GHFPWA or GHFPWG and applied to a representative personnel evaluation problem that involves a prioritization relationship over the evaluation index. Xu and Zhang [28] developed a novel approach based on TOPSIS (Technique for Order Preference by Similarity to Ideal Solution) and the maximizing deviation method for solving multiple criteria decision making (MCDM) problems, in which the evaluation information provided by the decision maker is expressed in HFEs and the information about criterion weights is incomplete.

Based on the MOORA (Multi-Objective Optimization by Ratio Analysis), the MULTIMOORA (Multi-Objective Optimization by Ratio Analysis plus the Full Multiplicative From) method was initially proposed by Brauers and Zavadskas [29, 30]. In general, the MULTIMOORA is regarded as a supplement for MOORA, which consists of the Ratio System and the Reference Point, with the Full Multiplicative Form. Brauers and Zavadskas [30] applied the MULTIMOORA method in MCDM process. Brauers and Zavadskas [31, 32] discussed the issues of information aggregation and data normalization in MULTIMOORA. The MOORA method has been employed for material selection [33, 34]. Brauers and Zavadskas [35] extended MULTIMOORA with triangular fuzzy numbers and Baležentis et al. [36] further proposed fuzzy MULTIMOORA for group decision making. Baležentis and Zeng [37] updated the MULTIMOORA method with T-2FSs, namely, generalized interval-valued trapezoidal fuzzy numbers. Balezentiene et al. [38] offered a MCDM framework for prioritization of energy crops based on fuzzy MULTIMOORA method which enables to tackle uncertain information. Streimikiene and Balezentis [39] applied the MULTIMOORA method for prioritization of the climate change mitigation strategies. Given selection of sustainable energy sources involves some conflicting criteria, Streimikiene et al. [40] employed multicriteria decision methods MULTIMOORA and TOPSIS to facilitate the analysis. In this paper, we propose an improved MULTIMOORA method with HFEs. By ameliorating the definitions of score function and distance measure for HFEs, the proposed method can address the hesitant fuzzy MCGDM problems and therefore reach more accurate decision results. In addition, the comparisons with other methods show the reasonability and efficiency of our algorithms.

The paper proceeds as follows. Section 2 discusses preliminaries for some basic notions, definitions, and properties of PA operator and HFSs. Section 3 presents novel definitions of the score function and the distance measure for HFEs. Section 4 introduces the crisp MULTIMOORA. Section 5 treats the MULTIMOORA method updated with HFEs and puts forward two algorithms for addressing MCGDM problems. Finally, Section 6 presents an application of the proposed methods in group decision making.

2. Preliminaries

This section is devoted to introduce some basic notions, definitions, and properties of PA operator and HFSs, which will facilitate our further analysis.

2.1. Power Average Operator

Information aggregation is a process that fuses data from multiple resources by using a proper aggregation technique. As one of the well-known aggregation techniques, the power-average (PA) operator, which is a nonlinear weighted-average aggregation tools whose weighting vectors depend on the input arguments, was initially introduced by Yager [41] and has been receiving extensive attention from researchers and practitioners over the past decades. This subsection aims to briefly introduce the notion of PA operator.

Definition 1 (see [41]). Let be a collection of crisp numbers; then, a nonlinear weighted-average aggregation tool, which is called PA operator, can be defined as follows: where and is the support for from , which satisfies the following properties: (1);(2);(3), if .

Generally, the support (i.e., Sup) measure is recognized as a similarity index. The more the similarity, the closer the two crisp numbers are and the more they support each other. Following Definition 1, the nonlinear weights of the input arguments are , . If the inherent information of is closer to the aggregated information, that is, the support of is considerably larger than others, then should be assigned more weight. Conversely, if the support of is considerably smaller than others, then should be assigned less weight. With the aid of the PA operator, the input arguments are allowed to be aggregated to support and reinforce each other. This operator allows for aggregations in which a subset of data clustered around a common value can combine in a nonlinear fashion to a ct in concert in determining the aggregated value. As such, the influence of the input arguments with large deviation will be reduced.

Two useful types of support function are proposed by Yager [41], and they both lie in the unit interval.

Theorem 2 (see [41]). Let be a collection of crisp numbers, and they will be aggregated by the PA operator. The support functions of the PA operator can be defined as (1);(2),where is the support for from and , .

To facilitate further application, this work assumes that is the calculation form of the support function when the input arguments take the form of crisp numbers.

2.2. Hesitant Fuzzy Sets

In this subsection, we review the basic concepts of hesitant fuzzy sets (HFSs) and hesitant fuzzy elements (HFEs).

Definition 3 (see [18]). A HFS on is defined in terms of a function when applied to , which returns a finite subset of ; that is, where is a set of some different values in , representing the possible membership degrees of the element to .

For convenience, we call the hesitant fuzzy element (HFE) and the set of all HFEs.

Given three HFEs , , and , Torra and Narukawa [17, 42] defined the following HFE operations.

Definition 4 (see [17, 42]). Let , , and be three HFEs defined on ; then, we have(1); (2); (3).

To facilitate the aggregation of hesitant fuzzy information, Xia and Xu [25] presented the following new operations on HFEs , , and .

Definition 5 (see [25]). Let , , and be three HFEs defined on ; then, we have(1), ;(2), ;(3); (4).

In order to compare HFEs, the comparison rules for HFEs defined by Xia and Xu [25] have been widely applied.

Definition 6 (see [25]). Let be a HFE defined on ; then, is called the score function of , where is the number of elements in . For any two HFEs, and , if , then ; if , then .

Let and be two HFEs. In most cases, . To make a comparison of these two HFEs, Xu and Xia [25] extended the shorter HFE until the length of both HFEs was the same, which in fact simply extended the shorter HFE by appending the same value repeatedly. In real application, the appended values reflect the risk preferences of DMs. Optimistic DMs expect desirable outcomes and would append the maximum value, and those with pessimistic anticipate unfavorable results and would append the minimum value [25]. In this paper, we assume that all of the DMs are pessimistic.

Another important concept for HFEs is the distance measure of HFEs and it has been applied to areas such as pattern recognition, cluster analysis, approximate reasoning, image processing, and decision making [8, 9, 2128, 43, 44]. We first introduce the axiom for distance measure.

Definition 7 (see [27]). Let and be two HFEs defined on ; then, the distance measure between and is defined as , which satisfies the following properties:(1);(2) if and only if ;(3).

A variety of distance measures of HFEs have been proposed based on Definition 7; two most commonly used of them, which are called hesitant normalized Hamming distance and hesitant normalized Euclidean distance, are introduced as follows.

Definition 8 (see [16]). Let and be two HFEs defined on ; then, is defined as the hesitant normalized Hamming distance between and , where and are the ith largest values in and and .

Definition 9 (see [16]). Let and be two HFEs defined on ; then, is defined as the hesitant normalized Euclidean distance between and , where and are the th largest values in and and .

2.3. Hesitant Fuzzy Power Average Operators

Generally, the power aggregation operators have been widely diffused to situations in which the input arguments are exact numerical values. Recently, Zhang [16] has extended the power average operator to accommodate hesitant fuzzy information as input.

Definition 10 (see [16]). Let be a collection of HFEs. If then is called the hesitant fuzzy power average (HFPA) operator, where and is the support for from , which satisfies the following three properties:(1); (2); (3), if , where is a distance measure between any two HFEs.

The support measure is essentially a similarity measure. It can be used to measure the proximity of a hesitant fuzzy preference value by one DM to another hesitant fuzzy preference value provided by a different DM. Without loss of generality, the support function for HFEs being used is in this study.

Theorem 11 (see [16]). Let be a collection of HFEs. The aggregated value using the operator is also a HFE, and

Definition 12 (see [16]). Let be a collection of HFEs. A generalized hesitant fuzzy power average (GHFPA) operator is a mapping such that where satisfies conditions (1) to (3) in Definition 7.

Theorem 13 (see [16]). Let be a collection of HFEs. The aggregated value using the HFPA operator is also a HFE, and

Definition 14 (see [16]). Let be a collection of HFEs. A weighted generalized hesitant fuzzy power average (WGHFPA) operator is a mapping such that where with the conditions that , for , and .

Theorem 15 (see [16]). Let be a collection of HFEs. The aggregated value using the operator is also a HFE, and

3. Novel Definitions of the Score Function and the Distance Measure for HFEs

In this section, we propose novel definitions of the score function and the distance measure for HFEs.

3.1. Score Function for HFEs

In fact, the score function advocated in Definition 6 can be interpreted as the simple arithmetic average of the HFEs if we assume that all the HFEs are of same weights. For a HFS , where the HFEs are listed in the set , it is obvious that four out of the five DMs believe the membership degree of to is greater than 0.82, and only one of them recognizes the membership degree as 0.1. In most cases, we can infer that the membership degree of to is probably greater than 0.82 rather than around 0.1. Therefore, the score function in Definition 6 may be unreasonable in some cases. In what follows, we define a novel definition of the score function of HFSs based on the PA operator.

Definition 16. Let be a HFE defined on ; then, we define as the score function of , where is the power average of .

The proposed score function can be simply interpreted as the power average of the elements of a HFE. An important characteristic of this score function is that it can reduce the influence caused by the HFE(s) of unduly large (or small) deviation. In another words, if a HFE or some of the HFSs is/are far from most of the others, then it/they will be assigned less weight when aggregated by the PA operator. From this point of view, the score function defined in this study will be more reasonable. In addition, the comparison rules of HFEs are same to that of Definition 6.

Example 17. Suppose that there are two HFEs and . According to Definition 6, the scores for and are, respectively, and . However, if we use the score function in Definition 16 to calculate their scores, we have , , , and for , and , , and for , based on which we obtain and .

From Example 17, we can know that the weights of each element in are 0.25 and the weights of each value in are 0.33 according to Definition 6. However, following Definition 16, we obtain that the weights of each element in are 0.2413, 0.2571, 0.2547, and 0.2469, respectively, and besides, the weights of each value in are 0.3167, 0.3537, and 0.3296, respectively. Obviously, the elements of unduly large deviation in , that is, 0.8 and 0.2, are assigned lower weights if we utilize the proposed score function, which means that the score function defined in this study can alleviate the influence of unduly large (or small) deviations on the aggregation results.

3.2. Distance Measure for HFEs

The hesitant normalized Hamming distance and the hesitant normalized Euclidean distance are widely applied in decision making. However, for two HFEs and , if , the hesitant normalized Hamming distance in Definition 8 and hesitant normalized Euclidean distance in Definition 9 will be somewhat unreasonable because some elements of the shorter HFE are appended subjectively. Hence, this work utilizes a new distance measure defined in the following.

Definition 18. Let and be two HFEs defined on ; then, is defined as distance measure between and , where and are, respectively, the power average of and .

The proof of the fact that the proposed distance measure satisfies the properties in Definition 7 is omitted here because it is trivial. Noteworthy, the proposed distance measure is more objective than those in Definitions 8 and 9, and it completely incorporates the inherent information of any two HFEs and , which can be illustrated in Example 19.

Example 19. Suppose that there are two HFEs and . According to Definitions 7 and 8, the distances between and are, respectively, and . However, if we use the distance measure in Definition 18 to calculate the distance between and , we have .

Using Definitions 8 and 9, the distance between and can only be calculated by extending the until the length of and was the same, which in fact simply extended the by appended the minimum value 0.1. Therefore, these two distances and are obtained based on the two HFEs and , which is subjective to some extent. Whereas the new distance does not need to append values to the shorter HFE until the length of the HFEs are the same, which makes the calculation more objective. In addition, the PA of and can reduce the influence of unduly large (or small) deviations on the aggregation results, so and are capable of representing the people’s hesitancy more reasonable and more accurate.

4. The Crisp MULTIMOORA Method

This section is devoted to briefly introduce the MULTIMOORA (MOORA plus full multiplicative form) method. As already mentioned earlier, Multi-Objective Optimization by Ratio Analysis (MOORA) method was first introduced in [29], and it was then extended to MULTIMOORA because the extension made it more robust.

Let be a matrix where its elements denote th alternative of jth objective as

In general, the criteria comprise benefit type and cost type. If is of benefit type, then the greater is the better it will be; if is of cost type, then the smaller is the better it will be. MOORA method consists of two parts: the Ratio System and the Reference Point approach. The MULTIMOORA method includes internal normalization and treats originally all the objectives equally important. In principle, all stakeholders interested in the issue only could give more importance to an objective. Therefore they could either multiply the dimensionless number representing the response on an objective with a significance coefficient or they could decide beforehand to split an objective into different subobjectives [2937, 45].

The Ratio System of MOORA. Ratio system defines data normalization by comparing alternative of an objective to all values of the objective: where denotes th alternative of th objective. These numbers usually belong to the interval ; that is, . The normalized matrix is denoted as . These indicators are added (if desirable value of indicator is maximum) or subtracted (if desirable value is minimum). Thus, the summarizing index of each alternative is derived in this way: where denotes number of objectives to be maximized. Then, every ratio is given the rank: the higher the index , the higher the rank.

The Reference Point of MOORA. Reference point approach is based on the Ratio System. The Maximal Objective Reference Point (vector) is found according to ratios found by employing (18). The th coordinate of the reference point can be described as in case of maximization. Every coordinate of this vector represents maximum or minimum of certain objective (indicator). Then, every element of normalized response matrix is recalculated and final rank is given according to deviation from the reference point and the Min–Max Metric of Tchebycheff:

The Full Multiplicative Form and MULTIMOORA. Brauers and Zavadskas [29] proposed MOORA to be updated by the Full Multiplicative Form method embodying maximization as well as minimization of purely multiplicative utility function. Overall utility of the ith alternative can be expressed as dimensionless number: where , denotes the product of objectives of the th alternative to be maximized with being the number of objectives to be maximized and where denotes the product of objectives of the th alternative to be minimized with being the number of objectives (indicators) to be minimized. Thus MULTIMOORA summarizes MOORA (i.e., Ratio System and Reference point) and the Full Multiplicative Form. The Dominance theory was proposed to summarize three ranks provided by respective parts of MULTIMOORA into a single one [35].

5. MULTIMOORA Based upon HFEs

The MULTIMOORA method is aimed at helping DMs analyze the problem as part of the decision-making process. In this section, we investigate an MCGDM problem, where the criterion values take the form of HFEs, the weight information on criteria is unknown. First, some notations are used to denote the indices, sets, and variables in an MCGDM problem in the following.(i): Total number of alternatives.(ii): Total number of criteria.(iii): Total number of DMs involved in the decision process.(iv): Index of alternatives.(v): Index of criteria.(vi): Index of DMs.(vii): The th alternative.(viii): A set of alternatives.(ix): The th criterion.(x): A set of criteria, which are considered to be independent.(xi): The th DM.(xii): A set of DMs.(xiii): Weight of the th DM, where , and .(xiv): Significant coefficient vector of criteria (which indicates the important degree of each criterion); suppose that , , and .(xv): A collection of benefit criteria (i.e., large values of s indicate great preference).(xvi): A collection of cost criteria (i.e., small values of s indicate great preference) such that and .(xvii): Evaluation on alternative concerning criterion that is given by decision maker and is a HFE.(xviii): Vector of hesitant fuzzy decision matrices with respect to all DMs.

In this work, we focus on addressing an MCGDM problem aimed at ranking alternatives and selecting the most desirable alternative(s) from a finite set based on the set . The basic process for solving this MCGDM problem is illustrated in Figure 1. The algorithm involves the following steps.

Algorithm 20. Consider the following.
Step 1. Each DM provides his or her evaluations and constructs his or her own hesitant decision matrix with elements being responses of alternatives.
Step 2. Utilize the HFPA operator (9) to aggregate all of the individual hesitant fuzzy decision matrices into the collective hesitant fuzzy decision matrix .
Step 3. Because all of the criteria are expressed by the hesitant fuzzy elements, the normalization procedure for the collective hesitant fuzzy decision matrix is not necessary.
Step 4 (the ratio system). The hesitant fuzzy elements in are added up for the benefit criteria and subtracted for the cost criteria, which can be represented by and , respectively. Following Definition 16, the scores for and can be calculated and denoted and , and then we obtain the overall utility of the th alternative in terms of the Ratio System, which can be represented by
The alternatives are ranked by comparing their overall utility values. Specifically, alternatives with higher overall utility value receive higher ranks.
Step 5 (the reference point approach). The Maximal Objective Reference Vector (MORV) is obtained according to the following two principles: (1) if the criteria are of benefit type, that is, , then , ; (2) if the criteria are of cost type, that is, , then .
Thereafter, the novel distance measure proposed in Definition 18 can be utilized to identify the collective deviation from the MORV for each alternative: Then, the alternatives are ranked in ascending order of .
Step 6 (the full multiplicative form). The hesitant fuzzy elements in are, respectively, multiplied for the benefit criteria and the cost criteria, which can be represented by and , respectively. Following Definition 16, the scores for and can be calculated and denoted and , and then we obtain the utility of the th alternative in terms of the Full Multiplicative Form, which can be represented by where denotes the product of the criteria of the th alternative to the maximized with being the set of criteria to be maximized, and denotes the product of criteria of the th alternative to the minimized with being the set of criteria to be minimized. The alternatives are ranked in descending order of .
Step 7. Utilize the Dominance theory [31] to aggregate the three ranks provided by respective parts of MULTIMOORA.
Algorithm 20 is designed for situations in which the information regarding the weights of the DMs and the significant coefficients of the criteria is unknown. For the situations that the information regarding the weights of the DMs and the significant coefficients of the criteria has been obtained, we utilize the WGHFPA operator to aggregate the individual hesitant fuzzy decision matrices and construct Algorithm 21 to deal with these situations.

Algorithm 21. Consider the following.
Step 1. Each DM provides his or her evaluations and constructs his or her own hesitant decision matrix with elements being responses of alternatives.
Step 2. Utilize the WGHFPA operator (12)to aggregate all of the individual hesitant fuzzy decision matrices into the collective hesitant fuzzy decision matrix .
Step 3. Because all of the criteria are expressed by the hesitant fuzzy elements, the normalization procedure for the collective hesitant fuzzy decision matrix is not necessary.
Step 4 (the ratio system). The hesitant fuzzy elements in are added up for the benefit criteria and subtracted for the cost criteria, which can be represented by and , respectively. Following Definition 16, the scores for and can be calculated and denoted and , and then we obtain the overall utility of the th alternative in terms of the Ratio System, which can be represented by
The alternatives are ranked by comparing their overall utility values. Specifically, alternatives with higher overall utility value receive higher ranks.
Step 5 (the reference point approach). The Maximal Objective Reference Vector (MORV) is obtained according to the following two principles: (1) if the criteria are of benefit type, that is, , then ; (2) if the criteria are of cost type, that is, , then .
Thereafter, the novel distance measure proposed in Definition 18 can be utilized to identify the collective deviation from the MORV for each alternative: Then, the alternatives are ranked in ascending order of .
Step 6 (the full multiplicative form). The hesitant fuzzy elements in are, respectively, multiplied for the benefit criteria and the cost criteria, which can be represented by and , respectively. Following Definition 16, the scores for and can be calculated and denoted and , and then we obtain the utility of the th alternative in terms of the Full Multiplicative Form, which can be represented by where denotes the product of the criteria of the th alternative to the maximized with being the set of criteria to be maximized, and denotes the product of criteria of the th alternative to the minimized with being the set of criteria to be minimized. The alternatives are ranked in descending order of .
Step 7. Utilize the Dominance theory [31] to aggregate the three ranks provided by respective parts of MULTIMOORA.

6. Numerical Example

In this section, the use of the proposed algorithm is illustrated by the following example.

Wang and Lee [46] considered a software selection problem in which the alternatives are the software packages to be selected and the criteria are the criteria under consideration. Let us reconsider this problem. Suppose that the manager of a computer center at a university wants to select a new information system to improve work productivity. Four alternatives remain on the candidate list after preliminary screening. Three decision makers are chosen from the computer center and form a decision-making committee with the weight vector . The computer center must make a decision according to the following four criteria: (1) the costs of the hardware/software investment ; (2) the contribution to the performance of the organization ; (3) the effort to transfer from the current system ; and (4) the reliability of the outsourcing software developer . Among the considered criteria, is of cost type, and are of benefit type. The significant coefficient vector of the criteria is .

In this context, the information regarding the weights of the DMs and the significant coefficients of the criteria has been obtained. However, it is common to encounter with situations in which the information regarding the weights of the DMs and the significant coefficients of the criteria is unknown; we first apply Algorithm 20 to support the DMs when making a decision by conducting the following steps.

Step 1. The decision makers evaluate the software packages with respect to the criteria and construct three hesitant fuzzy decision matrices (see Tables 1, 2, and 3), where is a HFE that denotes all of the possible values for the alternative regarding the criterion .

Step 2. Calculate the supports for each HFE from the other HFEs and they are shown in Table 4.
Utilize the HFPA operator to aggregate all of the individual hesitant fuzzy decision matrices into the collective hesitant fuzzy decision matrix (see Table 5).

Step 3 (the ratio system). Calculate the scores for and and we have
Then, we obtain the overall utility of the th alternative in terms of the Ratio System, which can be represented by
The alternatives are ranked by comparing their overall utility values, and we obtain (“” means “more important than”).

Step 4 (the reference point approach). The MORV can be obtained from the collective hesitant fuzzy matirx, which is as follows:
Thereafter, the proposed distance measure is utilized to identify the collective deviation from the MORV for each alternative: Then, the alternatives are ranked in descending order of and we have .

Step 5 (the full multiplicative form). The product of the benefit criteria of the th alternative and the product of cost criteria of the ith alternative are calculated and Then, The alternatives are ranked in descending order of and we have .

Step 6. According to ranks obtained by the ratio system, the reference point and the Full Multiplicative Form (see Table 6), the Dominance theory is employed to summarize these three ranks provided by respective parts of MULTIMOORA.

By conducting the multiple criteria evaluation process using the MULTIMOORA method, alternative is identified as the best choice for the DMs.

Back to the concrete context of the software selection problem, if the information regarding the weights of the DMs and the significant coefficients of the criteria has been obtained, then we can apply Algorithm 21 to support the DMs when making a decision by conducting the following steps.

Step 1. It is the same as Step 1 in Algorithm 20 and hence is omitted here.

Step 2. Calculate the supports for each HFE from the other HFEs and they are shown in Table 7.
Utilize the WGHFPA operator (let ) to aggregate all of the individual hesitant fuzzy decision matrices into the collective hesitant fuzzy decision matrix (see Table 8).

Step 3 (the ratio system). Calculate the scores for and and we have
Then, we obtain the overall utility of the ith alternative in terms of the Ratio System, which can be represented by
The alternatives are ranked by comparing their overall utility values, and we obtain (“” means “more important than”).

Step 4 (the reference point approach). The MORV can be obtained from the collective hesitant fuzzy matrix, which is as follows:
Thereafter, the proposed distance measure is utilized to identify the collective deviation from the MORV for each alternative: Then, the alternatives are ranked in descending order of and we have .

Step 5 (the full multiplicative form). The product of the benefit criteria of the th alternative and the product of cost criteria of the th alternative are calculated and Then, The alternatives are ranked in descending order of and we have .

Step 6. According to ranks obtained by the ratio system, the reference point and the Full Multiplicative Form (see Table 9), the Dominance theory is employed to summarize these three ranks provided by respective parts of MULTIMOORA.

By conducting the multiple criteria evaluation process using the MULTIMOORA method, alternative is identified as the best choice for the DMs.

When varies, we can obtain different results (see Table 10). The DMs can choose the values of according to their preferences.

In [16], the application of Zhang’s method in this example shows that alternative is identified as the best choice for the DMs when , which is different from the results obtained through the two algorithms proposed in this study. Particularly, the aggregation operator utilized in Algorithm 21 is the same to that used in Approach 1 in [16]. If , by comparing the final ranking of the alternatives obtained by our method and Zhang’s method, we find out that the rankings obtained by our approach, say , is totally different from that of Zhang, say . With the change of the attitudes of the DMs, we can clearly observe that the most desirable alternative will be varied from one to another; this result is therefore more reasonable than that in [16] because the attitudes of the DMs can be effectively reflected while using our method. We provide below a summary of the desirable advantages of our approach over the method of [16].(i)This study refines the definitions of the score function and distance measure for HFEs and establishes their novel formulas based on PA operator. The proposed score function and distance measure are remarkably different from those proposed in [16] and are investigated highly effective and reasonable because they can alleviate the influence caused by the HFE(s) that is/are of unduly large (or small) deviation. Therefore, the proposed novel definitions will enhance the accuracy of the improved MULTIMOORA method with HFEs and finally lead to precise decision results.(ii)We propose two algorithms for addressing MCGDM problems to respectively accommodate situations in which the information regarding the weights of the DMs and the significant coefficients of the criteria is unknown or known. That is, our approach is also capable of dealing with situations where the weight and significant coefficient information of the alternatives is unknown or completely known. Furthermore, the improved MULTIMOORA framework proposed in this paper can be flexibly extended to other MCGDM processes, in which the criteria take the forms of high-order fuzzy sets, random variables, fuzzy random variables, and so forth, for obtaining satisfactory decision results. As such, the proposed approaches have a wide applied future.(iii)We utilize the HFPA and WGHFPA operators in this study to fuse all of the individual hesitant fuzzy decision matrices into the collective hesitant fuzzy decision matrix. Our approach can make the ranking of the alternatives convertible according to the risk attitudes of DMs because the values of the parameters represent the optimistic, pessimistic or neural levels of DMs.

7. Conclusions

In this paper, we have defined the novel score function and distance measure for HFEs based on PA operator in order to alleviate the influence of unduly large (or small) deviation during the aggregation process. We further extend the MULTIMOORA method by employing HFEs and propose two algorithms for addressing MCGDM issues. Similar to the crisp MULTIMOORA method, the proposed MULTIMOORA method based on HFEs also consists of the three parts, namely the Ratio System, the Reference Point, and the Full Multiplicative Form, representing different approaches of data aggregation and ranking of alternatives. Therefore, the proposed method could not only provide the means for MCDM related to uncertain assessments, but also facilitate MCGDM with the aid of utilization of information aggregation operators.

A numerical example of software selection demonstrates the possibilities of application of the proposed method. Comparison with other methods further illustrates the reasonability and efficiency of our algorithms. Future studies should focus on development of the MULTIMOORA method in terms of interval-valued hesitant fuzzy sets.

Conflict of Interests

The author declares that they have no competing interests.

Acknowledgment

This study was financially supported by the one-hundred-talents scheme of Southwest Jiaotong University of China for Yong Scholars.