Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Multi-Granulation Neutrosophic Rough Sets on a Single Domain and Dual Domains with Applications
Next Article in Special Issue
Multiple-Attribute Decision-Making Method Using Similarity Measures of Hesitant Linguistic Neutrosophic Numbers Regarding Least Common Multiple Cardinality
Previous Article in Journal
Vegetation Greening for Winter Oblique Photography Using Cycle-Consistence Adversarial Networks
Previous Article in Special Issue
Selecting the Optimal Mine Ventilation System via a Decision Making Framework under Hesitant Linguistic Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Picture Hesitant Fuzzy Set and Its Application to Multiple Criteria Decision-Making

1
School of Transportation and Logistics, Southwest Jiaotong University, Chengdu 610031, China
2
National Lab of Railway Transportation, Southwest Jiaotong University, Chengdu 610031, China
*
Author to whom correspondence should be addressed.
Symmetry 2018, 10(7), 295; https://doi.org/10.3390/sym10070295
Submission received: 21 June 2018 / Revised: 17 July 2018 / Accepted: 18 July 2018 / Published: 20 July 2018
(This article belongs to the Special Issue Fuzzy Techniques for Decision Making 2018)

Abstract

:
To address the complex multiple criteria decision-making (MCDM) problems in practice, this article proposes the picture hesitant fuzzy set (PHFS) theory based on the picture fuzzy set and the hesitant fuzzy set. First, the concept of PHFS is put forward, and its operations are presented, simultaneously. Second, the generalized picture hesitant fuzzy weighted aggregation operators are developed, and some theorems and reduced operators of them are discussed. Third, the generalized picture hesitant fuzzy prioritized weighted aggregation operators are put forward to solve the MCDM problems that the related criteria are at different priorities. Fourth, two novel MCDM methods combined with the proposed operators are constructed to determine the best alternative in real life. Finally, two numerical examples and an application of web service selection are investigated to illustrate the effectiveness of the proposed methods. The sensitivity analysis shows that the different values of the parameter λ affect the ranking of alternatives, and the proposed operators are compared with several existing MCDM methods to illustrate their advantages.

1. Introduction

Multiple criteria decision-making (MCDM) problems occur in numerous practical fields [1,2,3]. For a specific purpose, several possible plans may be presented as the alternatives; then, decision makers assess the alternatives concerning the related criteria to determine the best one. Traditionally, crisp numbers are utilized to express the evaluation information. However, in real life, the data are inevitably incomplete and complex, and decision makers may be uncertain when evaluating the alternatives. To deal with the fuzziness of evaluation information, the fuzzy set (FS) [4] was proposed to improve the information form. During the past decades, many scholars devoted themselves to the study of the fuzzy MCDM problems [1]. Furthermore, in recent years, along with the complexity of the MCDM problems, how to improve the FS theory to deal with different specific situations has been a hot topic.
Although FS is a valid form to express the uncertain evaluation information, it cannot solve several complex situations in real life. For more effective expression of the evaluation information, many generalized forms of FS were proposed [5,6,7,8,9,10]. The purpose of this paper is to propose a new information form; the picture hesitant fuzzy set (PHFS) theory is put forward combined with the concepts of picture fuzzy set (PFS) [7] and hesitant fuzzy set (HFS) [8]. As a generalized form of FS, intuitionistic fuzzy set (IFS) [5], PFS, and HFS, PHFS can express the uncertainty and complexity of human opinions in practice; furthermore, the positive, neutral, negative, and refusal membership degrees are represented by several possible values that are given by decision makers.
In practice, the uncertain and complex evaluation information will be inevitably given by decision makers. For example, ten business managers discuss an investment project; five suggest agreement, two present disagreement, and the other business managers choose to abstain. Obviously, FS can only indicate the membership degree of evaluation information; thus, the opinions of the 10 business managers cannot be represented by FS. For overcoming the limitation of FS, Atanassov [5] put forward the non-membership function and developed the IFS. Then, the evaluation information in the aforementioned example can be expressed by IFS accurately. Later, the interval numbers were used to substitute the crisp numbers in IFS; then, the interval-valued intuitionistic fuzzy set (IVIFS) was developed [6]. To convey the indeterminate information of decision makers more effectively, Ye [9] and Liu and Yuan [11] extended the FS to triangular and trapezoidal intuitionistic fuzzy set, respectively. However, in some particular situations, it is not convincing to represent the evaluation information combined with IFS or IVIFS. For instance, there is a vote for a specific matter, the voting opinions of voters can be divided into four types, namely, vote for, abstain, vote against, and a refusal of the voting [12]. Therefore, Cuong [7,13] put forward the PFS, which is composed by the positive, neutral, negative, and refusal membership functions; thus, PFS can express the opinions of decision makers accurately in the example above. Subsequently, the correlation coefficient, distance measure, and cross-entropy measure of PFS were investigated in detail [14,15,16].
On the other hand, sometimes the accurate membership degree of evaluation information is difficult to be determined, which is also another shortcoming of FS. Therefore, the HFS was developed [17], in which the membership degrees are represented by several possible crisp numbers. Next, the interval numbers were introduced to extend the membership function of HFS and the interval-valued hesitant fuzzy set (IVHFS) theory was proposed [18]. According to the IFS and HFS, several potential membership and non-membership functions were expressed to put forward the dual hesitant fuzzy set (DHFS) [10]. Later, Farhadinia [19] constructed the dual interval-valued hesitant fuzzy set (DIVHFS) combined with DHFS. Nevertheless, HFS in the existing research cannot express all types of human opinions in the aforementioned example.
According to the evaluation information of the individual decision makers, the collective evaluation information of each alternative is obtained through the information fusion. Due to the important role of aggregation tools in MCDM problems, many scholars have investigated the aggregation operators of different fuzzy information. For example, Xu and Yager [20] developed the operations of intuitionistic fuzzy numbers (IFNs) and proposed the intuitionistic fuzzy geometric aggregation operators. Later, Xu [21] put forward the intuitionistic fuzzy weighted averaging aggregation operators to aggregate the IFNs. Next, several interval-valued intuitionistic fuzzy aggregation operators were constructed to deal with the MCDM [22,23,24]. With respect to the picture fuzzy (PF) evaluation information, Wei [25] defined the operations of picture fuzzy numbers (PFNs) according to the study of [21] and proposed the picture fuzzy weighted aggregation operators. In addition, several PF aggregation operators according to different operations were put forward [12,26]. Besides, a great time of hesitant fuzzy aggregation operators and their generalized forms were constructed [27], and several aggregation operators under dual hesitant fuzzy and dual interval-valued hesitant fuzzy environment were developed [28,29,30].
In some practical MCDM problems, the related criteria may be at different priority levels. For instance, a young couple wants to choose a toy for their child, the criteria of the toy they will consider are safety and price; obviously, the criteria safety has a higher priority than price. However, the aforementioned aggregation operators cannot fuse the aggregated arguments that are in different priority levels. In response to these situations, Yager [31] proposed the prioritized averaging (PA) operator. Inspired by Yager [31], Yu et al. [32,33] constructed the intuitionistic fuzzy prioritized fuzzy and interval-valued intuitionistic fuzzy prioritized fuzzy aggregation operators. Besides, the hesitant fuzzy prioritized aggregation operators were proposed to aggregate the evaluation information that is at different priorities [34]. Nevertheless, to our best knowledge, few researches have extended the PA operator to solve the MCDM problems under PF environment.
In summary, this paper defines the PHFS based on the PFS and HFS and develops the operations laws of picture hesitant fuzzy elements (PHFEs) according to the operations of IFNs [21]. Then, the generalized picture hesitant fuzzy aggregation operators and generalized picture hesitant fuzzy prioritized aggregation operators are put forward, and the properties and reduced operators of them are investigated. Furthermore, the proposed operators are utilized to solve diverse situations during MCDM processes under picture hesitant fuzzy (PHF) environment.
The rest of this paper is structured as follows. Definitions of the PFS, HFS, and PA operator are presented in Section 2. The concept of PHFS is defined, and the comparison method and operations of PHFEs are proposed in Section 3. Section 4 constructs the generalized picture hesitant fuzzy weighted averaging (GPHFWA), generalized picture hesitant fuzzy weighted geometric (GPHFWG), generalized picture hesitant fuzzy prioritized weighted averaging (GPHFPWA), and generalized picture hesitant fuzzy prioritized weighted geometric (GPHFPWG) operators. In Section 5, two MCDM methods are constructed according to the proposed operators. Section 6 applies the proposed methods into two numerical examples and an application of web service selection to show the effectiveness and advantages of the proposed methods. Finally, some conclusions are summarized in Section 7.

2. Preliminaries

To make this paper as self-contained as possible, we recall the definitions of the PFS, HFS, and PA operator, which will be utilized in the subsequent research.

2.1. PFS

Atanassov [5] applied the non-membership degree to extend FS; however, expressing the evaluation information depend on IFS is unreasonable in practice, at times. Therefore, Cuong [13] proposed the PFS theory based on FS and IFS, which can represent more information of decision makers, including yes, abstain, no, and refusal.
Definition 1.
Let X be a non-empty and finite set, a PFS P on X is defined by
P = { x , μ P ( x ) , η P ( x ) , v P ( x ) | x X } ,
where μ P ( x ) , η P ( x ) , and v P ( x ) are the positive, neutral, and negative membership functions that are belonging to [ 0 , 1 ] , respectively, and they meet the condition of 0 μ P ( x ) + η P ( x ) + v P ( x ) 1 . Furthermore, π P ( x ) = 1 μ P ( x ) η P ( x ) v P ( x ) is the refusal membership function.
Definition 2.
A picture fuzzy number (PFN) is represented by a = ( μ a , η a , v a ) , where μ a [ 0 , 1 ] , η a [ 0 , 1 ] , v a [ 0 , 1 ] , and μ a + η a + v a 1 [25].
Wei [25] proposed the operations of PFNs based on the operations of IFNs in [21].
Definition 3.
Let a 1 = ( μ 1 , η 1 , v 1 ) , a 2 = ( μ 2 , η 2 , v 2 ) , and a = ( μ , η , v ) be three PFNs, λ > 0 , and a c is the complementary set of a , then [25]
a c = ( v , η , μ ) ;  
a 1 a 2 = ( μ 1 + μ 2 μ 1 μ 2 , η 1 η 2 , v 1 v 2 ) ;  
a 1 a 2 = ( μ 1 μ 2 , η 1 + η 2 η 1 η 2 , v 1 + v 2 v 1 v 2 ) ;  
λ a = ( 1 ( 1 μ ) λ , η λ , v λ ) ;  
a λ = ( μ λ , 1 ( 1 η ) λ , 1 ( 1 v ) λ ) .

2.2. HFS

Due to the complexity of the evaluated object in practice, decision makers may have difficulty determining an accurate value of the membership level. To deal with this situation, Torra [8] developed the HFS theory in which the membership degree is expressed by several possible values.
Definition 4.
Let ( [ 0 , 1 ] ) be the set of all subsets of the unitary interval and X be a non- empty set. Let h A : X ( [ 0 , 1 ] ) , then an HFS A on X is defined by
A = { x , h A ( x ) | x X } .
Definition 5.
A hesitant fuzzy element (HFE) is a non-empty and finite subset of [ 0 , 1 ] [27].
Although a HFE can be given by any subset of [ 0 , 1 ] , in practice, HFS is commonly restricted to finite set in the MCDM problems [27]. Therefore, Bedregal at al. [35] proposed the typical hesitant fuzzy set (THFS), which is the finite and non-empty HFS. Later, Alcantud and Torra [36] defined the uniformly typical hesitant fuzzy set (UTHFS) that can simplify many theoretical and practical arguments, which is a generalized form of THFS. In this paper, the evaluation information of decision makers is expressed by UTHFS during the MCDM processes under hesitant fuzzy environment.
Definition 6.
Let H ( [ 0 , 1 ] ) be the set of all finite and non-empty subsets of [ 0 , 1 ] , and let X be a non- empty set. Then, a THFS A on X is defined by Equation (7), where h A : X H . Each h H is called a typical hesitant fuzzy element (THFE) [35].
Definition 7.
Let A be a THFS on X , if there is N such that the cardinality of the THFS l A ( x ) N for each x X . Then, the THFS A is an UTHFS. Each h H is called an uniformly typical hesitant fuzzy element (UTHFE) [36].
To aggregate the hesitant fuzzy evaluation information, Xia and Xu [27] investigated the operations of HFEs, which is also valid for fusing UTHFEs.
Definition 8.
Let h , h 1 , and h 2 be three UTHFEs, λ > 0 , then
h 1 h 2 = γ 1 h 1 , γ 2 h 2 { γ 1 + γ 2 γ 1 γ 2 } ;  
h 1 h 2 = γ 1 h 1 , γ 2 h 2 { γ 1 γ 2 } ;  
λ h = γ h { 1 ( 1 γ ) λ } ;  
h λ = γ h { γ λ } .

2.3. The PA Operator

Aggregation operator plays a crucial role in the process of information fusion. Sometimes, the criteria have different priorities according to their important degree; thus, Yager [31] constructed the PA operator to address these situations.
Definition 9.
Let C = { C 1 , C 2 , , C n } be a set of criteria, which are divided into several priority levels, i.e., the priority of C p is higher than C q when p < q . The C j ( x ) [ 0 , 1 ] is the evaluation value of the alternative x concerning the criteria C j . Thus, the PA operator is expressed by
P A ( C 1 ( x ) , C 2 ( x ) , , C n ( x ) ) = j = 1 n w j C j ( x ) .
where w j = T j / j = 1 n T j , T j = k = 1 j 1 C k ( x ) , and T 1 = 1 .

3. PHFS

According to the PFS and UTHFS, we can define the PHFS that is composed by four membership functions, namely, positive, neutral, negative, and refusal membership functions. The four membership degrees are denoted by several values belonging to [ 0 , 1 ] , respectively, which can convey the hesitancy of decision makers.
Definition 10.
Let X be a non-empty and finite set, a PHFS N on X is defined by
N = { x , μ ˜ ( x ) , η ˜ ( x ) , v ˜ ( x ) | x X } ,
where μ ˜ ( x ) = { α | α μ ˜ ( x ) } , η ˜ ( x ) = { β | β η ˜ ( x ) } , and v ˜ ( x ) = { γ | γ v ˜ ( x ) } are three sets of several values in [ 0 , 1 ] , representing the potential positive, neutral, and negative membership degrees. The degrees above satisfy the condition of 0 α + + β + + γ + 1 , where α + = α μ ˜ ( x ) max { α } , β + = β η ˜ ( x ) max { β } , and γ + = γ v ˜ ( x ) max { γ } . For convenience, we call n ˜ = { μ ˜ ( x ) , η ˜ ( x ) , v ˜ ( x ) } is a PHFE, denoted by n ˜ = { μ ˜ , η ˜ , v ˜ } .
During the process of applying the PHFEs to the practical MCDM problems, it is necessary to rank the PHFEs; thus, we develop the score and accuracy functions of PHFEs.
Definition 11.
Let n ˜ = { μ ˜ , η ˜ , v ˜ } be a PHFE, the numbers of values in μ ˜ , η ˜ , v ˜ are l , p , q , respectively. Thus, the score function is defined as
s ( n ˜ ) = ( 1 + 1 l i = 1 l α i 1 p i = 1 p β i 1 q i = 1 q γ i ) / 2 , s ( n ˜ ) [ 0 , 1 ] .
the accuracy function is expressed as
h ( n ˜ ) = 1 l i = 1 l α i + 1 p i = 1 p β i + 1 q i = 1 q γ i , h ( n ˜ ) [ 0 , 1 ] .
Based on the score and accuracy values of PHFEs, we can determine the order relations between two PHFEs as in the following.
Definition 12.
Let n ˜ 1 = { μ ˜ 1 , η ˜ 1 , v ˜ 1 } and n ˜ 2 = { μ ˜ 2 , η ˜ 2 , v ˜ 2 } be two PHFEs, then
(1) 
If s ( n ˜ 1 ) > s ( n ˜ 2 ) , then n ˜ 1 > n ˜ 2 ;
(2) 
If s ( n ˜ 1 ) = s ( n ˜ 2 ) , then
a.
If h ( n ˜ 1 ) > h ( n ˜ 2 ) , then n ˜ 1 > n ˜ 2 ;
b.
If h ( n ˜ 1 ) = h ( n ˜ 2 ) , then n ˜ 1 = n ˜ 2 ;
c.
If h ( n ˜ 1 ) < h ( n ˜ 2 ) , then n ˜ 1 < n ˜ 2
For example, let n ˜ 1 = { { 0.3 , 0.4 } , { 0.2 } , { 0.2 , 0.3 } } ± and n ˜ 2 = { { 0.3 } , { 0.2 , 0.3 } , { 0.1 , 0.2 } } be two PHFEs, according to the Definition 11, we have s ( n ˜ 1 ) = s ( n ˜ 2 ) = 0.45 , h ( n ˜ 1 ) = 0.4 , and h ( n ˜ 2 ) = 0.35 , then n ˜ 1 > n ˜ 2 .
Inspired by the operational laws of PFNs and UTHFEs, i.e., the Definition 3 and 8, we propose the operational laws of PHFEs as follows.
Definition 13.
Let n ˜ = { μ ˜ , η ˜ , v ˜ } , n ˜ 1 = { μ ˜ 1 , η ˜ 1 , v ˜ 1 } , and n ˜ 2 = { μ ˜ 2 , η ˜ 2 , v ˜ 2 } be three PHFEs, λ > 0 , and n ˜ c is the complementary set of n ˜ , and the operations of PHFEs are represented as
n ˜ c = α μ ˜ , β η ˜ , γ v ˜ { { γ } , { β } , { α } } ;  
n ˜ 1 n ˜ 2 = { μ ˜ 1 μ ˜ 2 , η ˜ 1 η ˜ 2 , v ˜ 1 v ˜ 2 } = α 1 μ ˜ 1 , β 1 η ˜ 1 , γ 1 v ˜ 1 , α 2 μ ˜ 2 , β 2 η ˜ 2 , γ 2 v ˜ 2 { { α 1 + α 2 α 1 α 2 } , { β 1 β 2 } , { γ 1 γ 2 } } ;  
n ˜ 1 n ˜ 2 = { μ ˜ 1 μ ˜ 2 , η ˜ 1 η ˜ 2 , v ˜ 1 v ˜ 2 } = α 1 μ ˜ 1 , β 1 η ˜ 1 , γ 1 v ˜ 1 , α 2 μ ˜ 2 , β 2 η ˜ 2 , γ 2 v ˜ 2 { { α 1 α 2 } , { β 1 + β 2 β 1 β 2 } , { γ 1 + γ 2 γ 1 γ 2 } } ;  
λ n ˜ = α μ ˜ , β η ˜ , γ v ˜ { { 1 ( 1 α ) λ } , { β λ } , { γ λ } } ;  
n ˜ λ = α μ ˜ , β η ˜ , γ v ˜ { { α λ } , { 1 ( 1 β ) λ } , { 1 ( 1 γ ) λ } } .
For example, let n ˜ 1 = { { 0.3 , 0.4 } , { 0.2 } , { 0.2 , 0.3 } } and n ˜ 2 = { { 0.3 } , { 0.2 , 0.3 } , { 0.1 , 0.2 } } be two PHFEs, λ = 2 , then
(1)
n ˜ 1 c = { { 0.2 , 0.3 } , { 0.2 } , { 0.3 , 0.4 } } , n ˜ 2 c = { { 0.1 , 0.2 } , { 0.2 , 0.3 } , { 0.3 } } ;
(2)
n ˜ 1 n ˜ 2 = { { 0.51 , 0.58 } , { 0.04 , 0.06 } , { 0.02 , 0.03 , 0.04 , 0.06 } } ;
(3)
n ˜ 1 n ˜ 2 = { { 0.09 , 0.12 } , { 0.36 , 0.44 } , { 0.28 , 0.36 , 0.37 , 0.44 } } ;
(4)
λ n ˜ 1 = { { 0.51 , 0.64 } , { 0.04 } , { 0.04 , 0.09 } } , λ n ˜ 2 = { { 0.51 } , { 0.04 , 0.09 } , { 0.01 , 0.04 } } ;
(5)
n ˜ 1 λ = { { 0.09 , 0.16 } , { 0.36 } , { 0.36 , 0.51 } } , n ˜ 2 λ = { { 0.09 } , { 0.36 , 0.51 } , { 0.19 , 0.36 } } .
Obviously, the following theorem can be obtained based on the Definition 13.
Theorem 1.
Let n ˜ = { μ ˜ , η ˜ , v ˜ } , n ˜ 1 = { μ ˜ 1 , η ˜ 1 , v ˜ 1 } , and n ˜ 2 = { μ ˜ 2 , η ˜ 2 , v ˜ 2 } be three PHFEs, λ , λ 1 , λ 2 > 0 , then
(1) 
n ˜ 1 n ˜ 2 = n ˜ 2 n ˜ 1 ;
(2) 
n ˜ 1 n ˜ 2 = n ˜ 2 n ˜ 1 ;
(3) 
λ ( n ˜ 1 n ˜ 2 ) = λ n ˜ 1 λ n ˜ 2 ;
(4) 
( n ˜ 1 n ˜ 2 ) λ = n ˜ 1 λ n ˜ 2 λ ;
(5) 
λ 1 n ˜ λ 2 n ˜ = ( λ 1 + λ 2 ) n ˜ ;
(6) 
n ˜ λ 1 n ˜ λ 2 = n ˜ ( λ 1 + λ 2 ) ;
(7) 
( n ˜ λ 1 ) λ 2 = n ˜ λ 1 λ 2 .

4. Generalized Picture Hesitant Fuzzy Aggregation Operators

Combined with the concept and operations of PHFS, the GPHFWA, GPHFWG, GPHFPWA, and GPHFPWG operators are developed. Then, several properties of them are discussed, and some other aggregation operators under PHF environment that reduced by the proposed operators are presented.

4.1. The GPHFWA Operator

Definition 14.
Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, the GPHFWA operator is a mapping Ω n Ω as
G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( w 1 n ˜ 1 λ w 2 n ˜ 2 λ w n n ˜ n λ ) 1 / λ = j = 1 n ( w j n ˜ j λ ) 1 / λ .
where w = ( w 1 , w 2 , , w n ) is the weight vector of PHFEs n ˜ j , and satisfies the conditions of w j > 0 and j = 1 n w j = 1 .
Based on the Definition 13, we can obtain the theorems as follows.
Theorem 2.
Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, then their aggregated value by using the GPHFWA operator is also a PHFE, and
G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = α 1 μ ˜ 1 , α 2 μ ˜ 2 , , α n μ ˜ n , β 1 η ˜ 1 , β 2 η ˜ 2 , , β n η ˜ n , γ 1 v ˜ 1 , γ 2 v ˜ 2 , , γ n v ˜ n { { ( 1 j = 1 n ( 1 α j λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 β j ) λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 γ j ) λ ) w j ) 1 / λ } } .
Proof. 
See Appendix A. □
Theorem 3.
(Idempotency) Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, if all the PHFEs are equal, i.e., n ˜ j = n ˜ = { μ ˜ , η ˜ , v ˜ } , μ ˜ = α , η ˜ = β , v ˜ = γ , then
G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = n ˜ = { μ ˜ , η ˜ , v ˜ } .
Proof. 
See Appendix B. □
Theorem 4.
(Boundedness) Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, if n ˜ = { { α } , { β + } , { γ + } } and n ˜ + = { { α + } , { β } , { γ } } , where α = α j μ ˜ j min { α j } , β = β j η ˜ j min { β j } , γ = γ j v ˜ j min { γ j } , α + = α j μ ˜ j max { α j } , β + = β j η ˜ j max { β j } , and γ = γ j v ˜ j min { γ j } , thus
n ˜ G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) n ˜ + .
Proof. 
See Appendix C. □
Theorem 5.
(Monotonicity) Let n ˜ j ( j = 1 , 2 , , n ) and n ˜ j * ( j = 1 , 2 , , n ) be two collections of PHFEs, if n ˜ j n ˜ j * , then
G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) G P H F W A λ ( n ˜ 1 * , n ˜ 2 * , , n ˜ n * ) .
Proof. 
Theorem 5 can be obtained by the Theorem 4. □
Under some specific situations, we can obtain the reduced operators of the GPHFWA operator.
Case 1.
If λ = 1 , then the GPHFWA operator is reduced to the picture hesitant fuzzy weighted averaging (PHFWA) operator
P H F W A ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( w 1 n ˜ 1 w 2 n ˜ 2 w n n ˜ n ) = j = 1 n ( w j n ˜ j ) .
Case 2.
If λ = 1 and w = ( 1 / n , 1 / n , , 1 / n ) , then the GPHFWA operator is reduced to the picture hesitant fuzzy arithmetic averaging (PHFAA) operator
P H F A A ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( 1 n n ˜ 1 1 n n ˜ 2 1 n n ˜ n ) .

4.2. The GPHFWG Operator

Similarly, the GPHFWG operator can be defined as in the following.
Definition 15.
Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, the GPHFWG operator is a mapping Ω n Ω as
G P H F W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = 1 λ ( λ n ˜ 1 w 1 λ n ˜ 2 w 2 λ n ˜ n w n ) = 1 λ j = 1 n ( λ n ˜ j w j ) .
where w = ( w 1 , w 2 , , w n ) is the weight vector of PHFEs n ˜ j , and satisfies the conditions of w j > 0 and j = 1 n w j = 1 .
According to the operational laws of PHFEs, the theorem can be obtained as follows.
Theorem 6.
Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, then their aggregated value by using the GPHFWG operator is also a PHFE, and
G P H F W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = α 1 μ ˜ 1 , α 2 μ ˜ 2 , , α n μ ˜ n , β 1 η ˜ 1 , β 2 η ˜ 2 , , β n η ˜ n , γ 1 v ˜ 1 , γ 2 v ˜ 2 , , γ n v ˜ n { { 1 ( 1 j = 1 n ( 1 ( 1 α j ) λ ) w j ) 1 / λ } , { ( 1 j = 1 n ( 1 β j λ ) w j ) 1 / λ } , { ( 1 j = 1 n ( 1 γ j λ ) w j ) 1 / λ } .
It can be proven by the same process as Theorem 3–5 that the GPHFWG operator also has several properties.
Theorem 7.
(Idempotency) Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, if all the PHFEs are equal, i.e., n ˜ j = n ˜ = { μ ˜ , η ˜ , v ˜ } , μ ˜ = α , η ˜ = β , v ˜ = γ , then
G P H F W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = n ˜ = { μ ˜ , η ˜ , v ˜ } .
Theorem 8.
(Boundedness) Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, if n ˜ = { { α } , { β + } , { γ + } } and n ˜ + = { { α + } , { β } , { γ } } , where α = α j μ ˜ j min { α j } , β = β j η ˜ j min { β j } , γ = γ j v ˜ j min { γ j } , α + = α j μ ˜ j max { α j } , β + = β j η ˜ j max { β j } , and γ = γ j v ˜ j min { γ j } , thus
n ˜ G P H F W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) n ˜ + .
Theorem 9.
(Monotonicity) Let n ˜ j ( j = 1 , 2 , , n ) and n ˜ j * ( j = 1 , 2 , , n ) be two collections of PHFEs, if n ˜ j n ˜ j * , then
G P H F W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) G P H F W G λ ( n ˜ 1 * , n ˜ 2 * , , n ˜ n * ) .
Several reduced operators of the GPHFWG operator are presented as:
Case 3.
If λ = 1 , then the GPHFWG operator is reduced to the picture hesitant fuzzy weighted geometric (PHFWG) operator
P H F W G ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( n ˜ 1 w 1 n ˜ 2 w 2 n ˜ n w n ) = j = 1 n ( n ˜ j w j ) .
Case 4.
If λ = 1 and w = ( 1 / n , 1 / n , , 1 / n ) , then the GPHFWG operator is reduced to the picture hesitant fuzzy geometric averaging (PHFGA) operator
P H F G A ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( n ˜ 1 n ˜ 2 n ˜ n ) 1 / n .

4.3. The GPHFPWA Operator

In real life, the criteria sometimes have different priority levels. For example, safety has a higher priority than price when a couple chooses a toy for their child. Obviously, the GPHFWA and GPHFWG operators cannot deal with this situation; then, the GPHFPWA and GPHFPWG operators are developed according to the PA operator proposed by Yager [31].
Definition 16.
Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, the GPHFPWA operator is a mapping Ω n Ω as
G P H F P W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( T 1 j = 1 n T j n ˜ 1 λ T 2 j = 1 n T j n ˜ 2 λ T n j = 1 n T j n ˜ n λ ) 1 / λ .
where T j = k = 1 j 1 s ( n ˜ k ) ( j = 2 , , n ) , T 1 = 1 , and s ( n ˜ k ) is the score value of PHFE n ˜ k .
Similarly, the following theorem can be put forward.
Theorem 10.
Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, then their aggregated value by using the GPHFPWA operator is also a PHFE, and
G P H F P W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = α 1 μ ˜ 1 , α 2 μ ˜ 2 , , α n μ ˜ n , β 1 η ˜ 1 , β 2 η ˜ 2 , , β n η ˜ n , γ 1 v ˜ 1 , γ 2 v ˜ 2 , , γ n v ˜ n { { ( 1 j = 1 n ( 1 α j λ ) T j j = 1 n T j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 β j ) λ ) T j j = 1 n T j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 γ j ) λ ) T j j = 1 n T j ) 1 / λ } } .
The GPHFPWA operator also has the properties as follows.
Theorem 12.
(Idempotency) Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, if all the PHFEs are equal, i.e., n ˜ j = n ˜ = { μ ˜ , η ˜ , v ˜ } , μ ˜ = α , η ˜ = β , v ˜ = γ , then
G P H F P W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = n ˜ = { μ ˜ , η ˜ , v ˜ } .
Theorem 13.
(Boundedness) Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, n ˜ = { { α } , { β + } , { γ + } } and n ˜ + = { { α + } , { β } , { γ } } , where α = α j μ ˜ j min { α j } , β = β j η ˜ j min { β j } , γ = γ j v ˜ j min { γ j } , α + = α j μ ˜ j max { α j } , β + = β j η ˜ j max { β j } , and γ = γ j v ˜ j min { γ j } , thus
n ˜ G P H F P W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) n ˜ + .
Theorem 14.
(Monotonicity) Let n ˜ j ( j = 1 , 2 , , n ) and n ˜ j * ( j = 1 , 2 , , n ) be two collections of PHFEs, if n ˜ j n ˜ j * , then
G P H F P W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) G P H F P W A λ ( n ˜ 1 * , n ˜ 2 * , , n ˜ n * ) .
Then, the reduced operators of the GPHFPWA operator can be obtained.
Case 5.
If λ = 1 , then the GPHFPWA operator is reduced to the picture hesitant fuzzy prioritized weighted averaging (PHFPWA) operator
P H F P W A ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( T 1 j = 1 n T j n ˜ 1 T 2 j = 1 n T j n ˜ 2 T n j = 1 n T j n ˜ n ) .
Case 6.
If λ = 1 and the criteria are at the same priority, then the GPHFPWA operator is reduced to the PHFWA operator
P H F W A ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( w 1 n ˜ 1 w 2 n ˜ 2 w n n ˜ n ) = j = 1 n ( w j n ˜ j ) .
Case 7.
If λ = 1 , w = ( 1 / n , 1 / n , , 1 / n ) , and the criteria are at the same priority, then the GPHFPWA operator is reduced to the PHFAA operator
P H F A A ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( 1 n n ˜ 1 1 n n ˜ 2 1 n n ˜ n ) .

4.4. The GPHFPWG Operator

Similarly, the GPHFPWG operator is constructed as below.
Definition 17.
Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, the GPHFPWG operator is a mapping Ω n Ω as
G P H F P W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = 1 λ ( ( λ n ˜ 1 ) T 1 j = 1 n T j ( λ n ˜ 2 ) T 2 j = 1 n T j ( λ n ˜ n ) T n j = 1 n T j ) .
where T j = k = 1 j 1 s ( n ˜ k ) ( j = 2 , , n ) , T 1 = 1 , and s ( n ˜ k ) is the score value of PHFE n ˜ k .
Combined with the operations of PHFEs, the following theorems are obtained.
Theorem 15.
Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, then their aggregated value by using the GPHFPWG operator is also a PHFE, and
G P H F P W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = α 1 μ ˜ 1 , α 2 μ ˜ 2 , , α n μ ˜ n , β 1 η ˜ 1 , β 2 η ˜ 2 , , β n η ˜ n , γ 1 v ˜ 1 , γ 2 v ˜ 2 , , γ n v ˜ n { { 1 ( 1 j = 1 n ( 1 ( 1 α j ) λ ) T j j = 1 n T j ) 1 / λ } , { ( 1 j = 1 n ( 1 β j λ ) T j j = 1 n T j ) 1 / λ } , { ( 1 j = 1 n ( 1 γ j λ ) T j j = 1 n T j ) 1 / λ } .
Theorem 16.
(Idempotency) Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, if all the PHFEs are equal, i.e., n ˜ j = n ˜ = { μ ˜ , η ˜ , v ˜ } , μ ˜ = α , η ˜ = β , v ˜ = γ , then
G P H F P W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = n ˜ = { μ ˜ , η ˜ , v ˜ } .
Theorem 17.
(Boundedness) Let n ˜ j ( j = 1 , 2 , , n ) be a collection of PHFEs, if n ˜ = { { α } , { β + } , { γ + } } and n ˜ + = { { α + } , { β } , { γ } } , where α = α j μ ˜ j min { α j } , β = β j η ˜ j min { β j } , γ = γ j v ˜ j min { γ j } , α + = α j μ ˜ j max { α j } , β + = β j η ˜ j max { β j } , and γ = γ j v ˜ j min { γ j } , thus
n ˜ G P H F P W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) n ˜ + .
Theorem 18.
(Monotonicity) Let n ˜ j ( j = 1 , 2 , , n ) and n ˜ j * ( j = 1 , 2 , , n ) be two collections of PHFEs, if n ˜ j n ˜ j * , then
G P H F P W G λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) G P H F P W G λ ( n ˜ 1 * , n ˜ 2 * , , n ˜ n * ) .
Several reduced operators of the GPHFPWG operator are presented as below:
Case 8.
If λ = 1 , then the GPHFPWG operator is reduced to the picture hesitant fuzzy prioritized weighted geometric (PHFPWG) operator
P H F P W G ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( ( n ˜ 1 ) T 1 j = 1 n T j ( n ˜ 2 ) T 2 j = 1 n T j ( n ˜ n ) T n j = 1 n T j ) .
Case 9.
If λ = 1 and the criteria are at the same priority, then the GPHFPWG operator is reduced to the PHFWG operator
P H F W G ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( n ˜ 1 w 1 n ˜ 2 w 2 n ˜ n w n ) = j = 1 n ( n ˜ j w j ) .
Case 10.
If λ = 1 , w = ( 1 / n , 1 / n , , 1 / n ) and the criteria are at the same priority, then the GPHFPWG operator is reduced to the PHFGA operator
P H F G A ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( n ˜ 1 n ˜ 2 n ˜ n ) 1 / n .

5. MCDM Methods under PHF Environment

We utilize the proposed operators to deal the different MCDM problems under PHF environment in this section. Let A = { A 1 , A 2 , , A m } be a collection of alternatives and C = { C 1 , C 2 , , C n } be a set of criteria; decision maker evaluates the m alternatives concerning the n criteria by using the PHFEs. Thus, suppose that N = ( n ˜ i j ) ( i = 1 , 2 , , m ; j = 1 , 2 , , n ) is the PHF evaluation matrix, and n ˜ i j = { μ ˜ i j , η ˜ i j , v ˜ i j } is the evaluation information when the alternative A i is evaluated concerning the criteria C j . In general, the criteria can be divided into two types in practice, namely, the cost criteria and benefit criteria; therefore, the evaluation information concerning the cost criteria should be transformed into the evaluation information concerning the benefit criteria to obtain the standardized PHF evaluation matrix N ¯ = ( n ¯ i j ) as
n ¯ i j = { n ˜ i j , for   the   benefit   criteria ; ( n ˜ i j ) c , for   the   cost   criteria .
According to the aforementioned assumptions, when the criteria of a specific MCDM problem are in same priority level, and let w = ( w 1 , w 2 , , w n ) be the weight vector of the criteria. We can construct a novel approach, i.e., Algorithm 1 to solve it based on the GPHFWA or the GPHFWG operator. The flow diagram of the Algorithm 1 is presented in Figure 1, and the ranking result can be obtained by the following steps.
Algorithm 1. MCDM method based on the GPHFWA or the GPHFWG operator.
1: Normalize the PHF evaluation matrix N to obtain the standardized PHF evaluation matrix N ¯ combined with Equation (51).
2: Utilize the GPHFWA operator
G P H F W A λ ( n ¯ i 1 , n ¯ i 2 , , n ¯ i n ) = ( w 1 n ¯ i 1 λ w 2 n ¯ i 2 λ w n n ¯ i n λ ) 1 / λ = j = 1 n ( w j n ¯ i j λ ) 1 / λ = α i 1 μ ˜ i 1 , α i 2 μ ˜ i 2 , , α i n μ ˜ i n , β i 1 η ˜ i 1 , β i 2 η ˜ i 2 , , β i n η ˜ i n , γ i 1 v ˜ i 1 , γ i 2 v ˜ i 2 , , γ i n v ˜ i n { { ( 1 j = 1 n ( 1 α ¯ i j λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 β ¯ i j ) λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 γ ¯ i j ) λ ) w j ) 1 / λ } }
or the GPHFWG operator
G P H F W G λ ( n ¯ i 1 , n ¯ i 2 , , n ¯ i n ) = 1 λ ( λ n ¯ i 1 w 1 λ n ¯ i 2 w 2 λ n ¯ i n w n ) = 1 λ j = 1 n ( λ n ¯ i j w j ) = α i 1 μ ˜ i 1 , α i 2 μ ˜ i 2 , , α i n μ ˜ i n , β i 1 η ˜ i 1 , β i 2 η ˜ i 2 , , β i n η ˜ i n , γ i 1 v ˜ i 1 , γ i 2 v ˜ i 2 , , γ i n v ˜ i n { { 1 ( 1 j = 1 n ( 1 ( 1 α ¯ i j ) λ ) w j ) 1 / λ } , { ( 1 j = 1 n ( 1 β ¯ i j λ ) w j ) 1 / λ } , { ( 1 j = 1 n ( 1 γ ¯ i j λ ) w j ) 1 / λ }
to aggregate the standardized PHF evaluation matrix N ¯ to obtain the collective evaluation information of each alternative, i.e., n ˜ i = { μ ˜ i , η ˜ i , v ˜ i } .
3: Compute the score and accuracy values of each alternative using Equation (14) and (15).
4: Based on the comparison method of PHFEs, rank the alternatives.
When the criteria are in different priorities, we can solve the MCDM problem combined with the Algorithm 2 based on the GPHFPWA or the GPHFPWG operator. The flow diagram of Algorithm 2 is presented in Figure 2, and the ranking result can be obtained by the following steps.
Algorithm 2. MCDM method based on the GPHFPWA or the GPHFPWG operator.
1: Normalize the PHF evaluation matrix N to obtain the standardized PHF evaluation matrix combined with Equation (51).
2: Compute the values of T i j using the equations as
T i j = k = 1 j 1 s ( n ˜ i k ) , T i 1 = 1 .
3: Utilize the GPHFPWA operator
G P H F P W A λ ( n ¯ i 1 , n ¯ i 2 , , n ¯ i n ) = ( T i 1 j = 1 n T i j n ¯ i 1 λ T i 2 j = 1 n T i j n ¯ i 2 λ T i n j = 1 n T i j n ¯ i n λ ) 1 / λ = α i 1 μ ˜ i 1 , α i 2 μ ˜ i 2 , , α i n μ ˜ i n , β i 1 η ˜ i 1 , β i 2 η ˜ i 2 , , β i n η ˜ i n , γ i 1 v ˜ i 1 , γ i 2 v ˜ i 2 , , γ i n v ˜ i n { { ( 1 j = 1 n ( 1 α ¯ i j λ ) T i j j = 1 n T i j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 β ¯ i j ) λ ) T i j j = 1 n T i j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 γ ¯ i j ) λ ) T i j j = 1 n T i j ) 1 / λ } }
or the GPHFPWG operator
G P H F P W G λ ( n ¯ i 1 , n ¯ i 2 , , n ¯ i n ) = 1 λ ( ( λ n ¯ i 1 ) T i 1 j = 1 n T i j ( λ n ¯ i 2 ) T i 2 j = 1 n T i j ( λ n ¯ i n ) T i n j = 1 n T i j ) = α i 1 μ ˜ i 1 , α i 2 μ ˜ i 2 , , α i n μ ˜ i n , β i 1 η ˜ i 1 , β i 2 η ˜ i 2 , , β i n η ˜ i n , γ i 1 v ˜ i 1 , γ i 2 v ˜ i 2 , , γ i n v ˜ i n { { 1 ( 1 j = 1 n ( 1 ( 1 α ¯ i j ) λ ) T i j j = 1 n T i j ) 1 / λ } , { ( 1 j = 1 n ( 1 β ¯ i j λ ) T i j j = 1 n T i j ) 1 / λ } , { ( 1 j = 1 n ( 1 γ ¯ i j λ ) T i j j = 1 n T i j ) 1 / λ }
to aggregate the standardized PHF evaluation matrix N ¯ to obtain the collective evaluation information of each alternative, i.e., n ˜ i = { μ ˜ i , η ˜ i , v ˜ i } .
4: Compute the score and accuracy values of each alternative using Equation (14) and (15).
5: Based on the comparison method of PHFEs, rank the alternatives.

6. Numerical Examples

We adopt two numerical examples of MCDM problems from the study of [25] and [34] and an application of web service selection [37] to show the feasibility and advantages of the proposed methods.

6.1. Implementation

Example 1.
Suppose that an organization wants to construct the enterprise resource planning (ERP) system [25]. After investigating the existing vendors of ERP systems on the market, five potential ERP systems are primary determined to be chosen from, i.e., Ai (i = 1, 2, 3, 4, 5). Decision makers utilize the PHFEs to evaluate the five alternatives with respect to four criteria, namely, function and technology (C1), strategic fitness (C2), ability of vendor (C3), and reputation of vendor (C4), and the weight vector of the criteria is w = (0.2, 0.1, 0.3, 0.4). Subsequently, the PHF evaluation matrix N = ( n ˜ i j ) is obtained as shown in Table 1.
Then, we can determine the ranking of the five potential ERP systems using the Algorithm 1, which are presented as below.
  • Step 1: Because of all the criteria are the benefit type, the standardized PHF evaluation matrix N ¯ is as same as the PHF evaluation matrix N .
  • Step 2: Use the GPHFWA ( λ = 1 ) operator to aggregate the standardized PHF evaluation matrix N ¯ , and the collective evaluation information of each alternative is obtained as
    n ˜ 1 = { { 0.3636 , 0.3877 , 0.4113 , 0.4336 } , { 0.3862 , 0.4048 , 0.4136 , 0.4335 } , { 0.0444 , 0.0482 , 0.0502 , 0.0544 } } ; n ˜ 2 = { { 0.4061 , 0.4401 , 0.4684 , 0.5023 , 0.5308 , 0.5545 } , { 0.2563 , 0.2612 , 0.2612 , 0.2659 , 0.2662 , 0.2709 , 0.2709 , 0.2761 } , { 0.0891 , 0.0942 , 0.0964 , 0.1019 } } ; n ˜ 3 = { { 0.4014 , 0.4789 , 0.5180 , 0.5230 , 0.5804 , 0.6159 } , { 0.1627 , 0.1677 , 0.1725 , 0.1763 , 0.1779 , 0.1870 , 0.1884 , 0.1943 , 0.1999 , 0.2042 , 0.2061 , 0.2166 } , { 0.0448 , 0.0551 } } ; n ˜ 4 = { { 0.3549 , 0.3738 , 0.4077 , 0.4251 } , { 0.3834 , 0.3989 , 0.4018 , 0.4181 , 0.4248 , 0.4420 , 0.4452 , 0.5632 } , { 0.0576 } } ; n ˜ 5 = { { 0.3468 , 0.4038 , 0.4756 } , { 0.3321 , 0.3444 , 0.3516 , 0.3647 } , { 0.0612 , 0.0642 } } .
  • Step 3: Compute the score values of each alternative combined with Equation (14):
    s ( n ˜ 1 ) = 0.4701 , s ( n ˜ 2 ) = 0.5611 , s ( n ˜ 3 ) = 0.6409 , s ( n ˜ 4 ) = 0.4553 , s ( n ˜ 5 ) = 0.4989 .
  • Step 4: According to the score values, the ranking result of the five ERP systems is determined as A 3 A 2 A 5 A 1 A 4 .
If the GPHFWG operator is utilized in the steps above to complete the information fusion, the ranking procedures are presented as follows.
  • Step 1′: See Step 1.
  • Step 2′: Use the GPHFWG ( λ = 1 ) operator to aggregate the standardized PHF evaluation matrix N ¯ , and the collective evaluation information of each alternative is obtained as
    n ˜ 1 = { { 0.2307 , 0.2343 , 0.2405 , 0.2443 } , { 0.5365 , 0.6663 , 0.5380 , 0.6673 } , { 0.0600 , 0.0797 , 0.0661 , 0.0856 } } ; n ˜ 2 = { { 0.2017 , 0.2101 , 0.2151 , 0.2212 , 0.2304 , 0.2358 } , { 0.4526 , 0.4550 , 0.4670 , 0.4693 , 0.4914 , 0.4936 , 0.5047 , 0.5070 } , { 0.0901 , 0.0993 , 0.0999 , 0.1090 } } ; n ˜ 3 = { { 0.1986 , 0.2057 , 0.2081 , 0.2342 , 0.2427 , 0.2454 } , { 0.3334 , 0.3363 , 0.3602 , 0.3630 , 0.3766 , 0.3792 , 0.4016 , 0.4042 , 0.4830 , 0.4852 , 0.5038 , 0.5059 } , { 0.0481 , 0.0634 } } ; n ˜ 4 = { { 0.1024 , 0.1037 , 0.1044 , 0.1058 } , { 0.5949 , 0.5977 , 0.6709 , 0.6732 , 0.7038 , 0.7058 , 0.7594 , 0.7611 } , { 0.0591 } } ; n ˜ 5 = { { 0.1613 , 0.1660 , 0.1696 } , { 0.5850 , 0.6372 , 0.6503 , 0.6943 } , { 0.0716 , 0.0805 } } .
  • Step 3′: Compute the score values of each alternative combined with Equation (14):
    s ( n ˜ 1 ) = 0.2813 , s ( n ˜ 2 ) = 0.3197 , s ( n ˜ 3 ) = 0.3778 , s ( n ˜ 4 ) = 0.1808 , s ( n ˜ 5 ) = 0.2240 .
  • Step 4′: According to the score values, the ranking result of the five ERP systems is determined as A 3 A 2 A 1 A 5 A 4 .
Example 2.
Suppose a university wants to introduce excellent foreign professors to improve the level of teaching and scientific research [34]. There are five foreign professors who are selected by the University’s human resources department. Based on the priority level, the criteria of investigation is successively morality (C1), research ability (C2), teaching capacity (C3), and educational experience (C4); a priority relationship C 1 C 2 C 3 C 4 exists between the criteria. Then, the PHF evaluation matrix N = ( n ˜ i j ) is presented in Table 2.
Subsequently, we can determine the ranking of the five foreign professors using the Algorithm 2, which are presented as follows.
  • Step 1: Because of all the criteria are the benefit type, the standardized PHF evaluation matrix N ¯ is as same as the PHF evaluation matrix N .
  • Step 2: Compute the values of T i j using the Equation (54)
    T i j = [ 1.000 0.6667 0.4783 0.3157 1.000 0.7600 0.5225 0.4311 1.000 0.7475 0.3526 0.1992 1.000 0.7100 0.3905 0.3417 1.000 0.7000 0.4620 0.3719 ] .
  • Step 3: Use the GPHFPWA ( λ = 1 ) operator to aggregate the standardized PHF evaluation matrix N ¯ , and the collective evaluation information of each alternative is obtained as
    n ˜ 1 = { { 0.5003 , 0.5177 , 0.5360 , 0.5382 , 0.5522 , 0.6230 , 0.6361 , 0.6516 , 0.6516 } , { 0.0495 , 0.0521 , 0.0562 , 0.0592 } , { 0.1176 , 0.1345 , 0.1558 , 0.1783 } } ; n ˜ 2 = { { 0.6290 , 0.6396 , 0.6446 , 0.6548 , 0.6723 , 0.6816 , 0.6861 , 0.6950 } , { 0.0460 , 0.0559 , 0.0594 , 0.0722 } , { 0.1084 , 0.1175 , 0.1186 , 0.1287 , 0.1316 , 0.1427 , 0.1441 , 0.1562 } } ; n ˜ 3 = { { 0.5335 , 0.5533 , 0.5537 , 0.5727 , 0.5996 , 0.6169 } , { 0.0494 , 0.0550 , 0.0617 , 0.0686 } , { 0.1692 , 0.1732 , 0.1857 , 0.1902 , 0.2018 , 0.2066 , 0.2216 , 0.2269 } } ; n ˜ 4 = { { 0.5593 , 0.5820 , 0.5978 , 0.6185 , 0.6425 , 0.6609 } , { 0.0921 , 0.1015 , 0.1055 , 0.1162 } , { 0.0947 , 0.1044 , 0.1159 , 0.1277 } } ; n ˜ 5 = { { 0.5888 , 0.6181 , 0.6429 , 0.6683 } , { 0.0605 , 0.0670 , 0.0796 , 0.0881 } , { 0.1670 , 0.1742 , 0.1871 , 0.1951 } } .
  • Step 4: Compute the score values of each alternative combined with Equation (14)
    s ( n ˜ 1 ) = 0.6888 , s ( n ˜ 2 ) = 0.7368 , s ( n ˜ 3 ) = 0.6580 , s ( n ˜ 4 ) = 0.6978 , s ( n ˜ 5 ) = 0.6874 .
  • Step 5: According to the score values, the ranking result of the five foreign professors is determined as A 2 A 4 A 1 A 5 A 3 .
If the GPHFPWG operator is utilized in the steps above to complete the information fusion, the ranking procedures are presented as follows.
  • Step 1′: See Step 1.
  • Step 2′: See Step 2.
  • Step 3′: Use the GPHFPWG ( λ = 1 ) operator to aggregate standardized the PHF evaluation matrix N ¯ , and the collective evaluation information of each alternative is obtained as
    n ˜ 1 = { { 0.4753 , 0.4963 , 0.5142 , 0.5204 , 0.5434 , 0.5966 , 0.6231 , 0.6455 , 0.6455 } , { 0.0527 , 0.0597 , 0.0609 , 0.0678 } , { 0.1203 , 0.1402 , 0.1614 , 0.1804 } } ; n ˜ 2 = { { 0.6049 , 0.6124 , 0.6267 , 0.6345 , 0.6376 , 0.6456 , 0.6606 , 0.6689 } , { 0.0668 , 0.0739 , 0.0809 , 0.0878 } , { 0.1176 , 0.1230 , 0.1350 , 0.1403 , 0.1462 , 0.1515 , 0.1630 , 0.1682 } } ; n ˜ 3 = { { 0.4297 , 0.4424 , 0.4902 , 0.5047 , 0.5788 , 0.5959 } , { 0.0526 , 0.0571 , 0.0703 , 0.0748 } , { 0.2020 , 0.2110 , 0.2216 , 0.2304 , 0.2410 , 0.2495 , 0.2596 , 0.2680 } } ; n ˜ 4 = { { 0.5026 , 0.5363 , 0.5416 , 0.5769 , 0.5779 , 0.6155 } , { 0.1119 , 0.1178 , 0.1264 , 0.1322 } , { 0.0994 , 0.1236 , 0.1297 , 0.1531 } } ; n ˜ 5 = { { 0.5596 , 0.5733 , 0.6141 , 0.6292 } , { 0.0640 , 0.0801 , 0.0838 , 0.0995 } , { 0.1791 , 0.1975 , 0.1985 , 0.2164 } } .
  • Step 4′: Compute the score values of each alternative combined with Equation (14)
    s ( n ˜ 1 ) = 0.6757 , s ( n ˜ 2 ) = 0.7080 , s ( n ˜ 3 ) = 0.6040 , s ( n ˜ 4 ) = 0.6550 , s ( n ˜ 5 ) = 0.6572 .
  • Step 5′: According to the score values, the ranking result of the five foreign professors is obtained as A 2 A 1 A 5 A 4 A 3 .

6.2. Sensitivity Analysis

To explore the impact of the parameter λ on the ranking results, different possible values of λ are used in the algorithms of two aforementioned numerical examples, such as 0.001, 0.5, 1, 2, 3, 5, 10, 20, and 50. Then, combined with the proposed methods, the different rankings of alternatives are presented in Table 3, Table 4, Table 5 and Table 6. From Table 3 and Table 4, we can find that the best potential ERP system in Example 1 is always A 3 using both the GPHFWA operator and GPHFWG operator; however, some differences exist between the ranking results concerning different values of λ . Table 5 and Table 6 show that when we utilize the GPHFPWA operator to complete the information fusion, the best foreign professor is A 2 for 0.001 λ 10 , but the best alternative is A 4 for 20 λ 50 . In addition, when the GPHFPWG operator is used in Algorithm 2, the best foreign professor is A 2 for 0.001 λ 3 , but the best alternative is A 1 for 5 λ 50 . On the other hand, the score values of all the alternatives vary with different values of λ ; the reason is that the aggregation processes of the proposed operators have changed. For instance, when λ = 2 , the GPHFWA operator can be reduced to the picture hesitant fuzzy weighted quadratic averaging (PHFWQA) operator as
P H F W Q A ( n ¯ i 1 , n ¯ i 2 , , n ¯ i n ) = ( w 1 n ¯ i 1 2 w 2 n ¯ i 2 2 w n n ¯ i n 2 ) 1 / 2 .
when λ = 3 , the GPHFWA operator can be reduced to the picture hesitant fuzzy weighted cubic averaging (PHFWCA) operator as
P H F W C A ( n ¯ i 1 , n ¯ i 2 , , n ¯ i n ) = ( w 1 n ¯ i 1 3 w 2 n ¯ i 2 3 w n n ¯ i n 3 ) 1 / 3 .
Besides, the following results can be obtained from Table 3, Table 4, Table 5 and Table 6:
(1)
In Example 1, the score values of each alternative obtained by the GPHFWA operator are bigger than those obtained by the GPHFWG operator, and the difference between them increases along with the increasing of λ . It means that the GPHFWA operator is more suitable to aggregate the PHFEs of optimistic decision makers, while the GPHFWG operator can reflect the opinion of pessimistic decision makers. Furthermore, the level of optimism and pessimism are greater with the bigger value of λ .
(2)
In Example 2, the score values of each alternative obtained by the GPHFPWA and GPHFPWG operators are relatively stable when the different values of λ are used; the parameter λ cannot reflect the attitude of decision makers. In addition, the best alternative varies when the value of λ is relatively high, while the best alternative is always the same in Example 1. It means that the rankings obtained by the GPHFPWA and GPHFPWG operators are more affected by the parameter λ than those obtained by the GPHFWA and GPHFWG operators.
The aforementioned sensitivity analysis results show that the value of λ plays a very important role in MCDM problems, especially when the value of λ is relatively high. The value of λ can be determined based on the personal preference of decision makers to obtain different ranking results; thus, the proposed methods are highly flexible to deal with different situations in practice.

6.3. Comparative Analysis

To prove the feasibility of the proposed MCDM methods, the rankings of Example 1 in this paper are compared with the rankings obtained by the existing MCDM methods as presented in Table 7; including the PFWA and PFWG operators [25], and the picture fuzzy cross-entropy method [16]. Similarly, a comparison of Example 2 between the GPHFPWA and GPHFPWG operators and the HFPWA and HFPWG operators [34] is presented in Table 8.
Table 7 shows that the best alternative of Example 2 obtained by the MCDM methods based on the GPHFWA and GPHFWG operators is always A 3 , which is consistent with the existing methods; the results can demonstrate the feasibility of the proposed method. Compared with the PFS that is used in the study of [25] and [16], PHFS proposed in this paper can convey the human opinions more effectively, including yes, abstain, no, and refusal. For instance, the evaluation information of the alternative A 1 concerning the criteria C 1 that are given by decision maker is expressed as a PFN (0.53,0.33,0.09) [16,25]. In practice, decision maker may feel doubtful to determine an exact value of each membership level. Obviously, PFS cannot deal with this situation; however, we can use PHFS to represent the evaluation information as a PHFE {{0.43,0.53}, {0.33}, {0.06,0.09}} as shown in Table 1. Consequently, the proposed method can solve the MCDM problems when decision makers feel difficulty to determine the accurate value of each membership level. On the other hand, when the numbers of the criteria are relatively large, the aggregation process of the proposed operators will be more complicated than the existing methods and the data size will be relatively large; it is the limitation of the proposed method. Table 8 shows that the best alternative of Example 2 obtained by the HFPWA operator is A 5 , but the result of other MCDM methods is A 2 . The main reason of the difference is that the MCDM methods combined with the HFPWA and HFPWG operators ignore some complex evaluation information of decision makers in practice. UTHFS allows the decision makers to give several values of positive membership level, for instance, the evaluation information of the alternative A 1 concerning the criteria C 1 that are given by decision maker is expressed as a UTHFE (0.4,0.5,0.7) [34]. Nevertheless, in some particular situations, it is not convincing to express the evaluation information that only considers the positive membership level of decision makers; many scholars have focused on this problem and made some improvements to UTHFS [10,19]. Thus, we can overcome the limitation of UTHFS combined with the proposed method. It is worth noting that the GPHFPWA and GPHFPWG operators also have the same disadvantage as the GPHFWA and GPHFWG operators.
According to the aforementioned comparison results, we can summarize the advantages and disadvantages of the different MCDM methods (see Table 9), as well as their respective fields of application (see Table 10). In addition, the benefits of the aggregation process by using the proposed operators are presented as in the following.
(1) The Expansion of the Evaluation Information
The GPHFWA, GPHFWG, GPHFPWA, and GPHFPWG operators can solve the MCDM problems under PHF environment. PHFS proposed in this paper can express the different human opinions in real life and allow the decision makers to give several possible values of the different membership levels; thus, it can simultaneously depict the uncertainty and hesitancy of decision makers’ evaluation information, which cannot be achieved by PFS and UTHFS. Therefore, when decision makers are not fully aware of the evaluation target and feel doubtful about each membership level, it is reasonable to deal with these MCDM problems combined with the proposed methods. Furthermore, as a generalized form of FS, IFS, PFS, and UTHFS, we can transform the proposed methods into the existing MCDM methods if necessary.
(2) The Flexibility of Information Aggregation with Different Values of λ
Recall the sensitivity analysis in Section 6.2, the proposed operators can be reduced to other specific PHF aggregation operators by varying the value of λ ; thus, the proposed methods are highly flexible to deal with different situations. Furthermore, the parameter λ can also be regarded as a measure of the optimism and pessimism level of decision makers in the information fusion of the GPHFWA and GPHFWG operators; and the value of λ can be determined by decision makers according to their preferences in practice.
(3) The Simplicity of Dealing with Different Types of Criteria
During the MCDM process, the weight values of criteria play an important role and will affect the final ranking results. The criteria can be divided into two categories: one is in the same priority, the other is in different priorities. On the one hand, when the criteria have the same priority level, we can utilize the proposed method based on the GPHFWA and GPHFWG operators combined with the weight vector of criteria to solve the MCDM problem. On the other hand, when the criteria have different priority levels, the GPHFPWA and GPHFPWG operators can be introduced to determine the ranking of alternatives. In practice, we can use different aggregation operators in this paper to deal with different situations.

6.4. Application of Web Service Selection

To investigate the applications of the proposed methods in a more realistic scenario, we use the proposed methods to solve the Quality of Service (QoS) based web service selection problem [37]. According to the study of [37], the evaluation information of QoS is measured by a crisp number scale of 1–9, and the related criteria are availability ( C 1 ), throughput ( C 2 ), successability ( C 3 ), reliability ( C 4 ), compliance ( C 5 ), best practices ( C 6 ), documentation ( C 7 ), latency ( C 8 ), and response time ( C 9 ). Due to the criteria latency and response time are the cost type criteria, the closer the evaluation values concerning these two criteria are to 1, the better the alternative.
Suppose there are 20 web services to be evaluated concerning the aforementioned nine criteria, i.e., W S i ( i = 1 , 2 , , 20 ) ; the evaluation information of each web service is presented in Table 11. As each evaluation value in [37] is expressed by an exact crisp number, the PHFS can be reduced to the PFS to represent the evaluation information of each web service. Based on the relationship between the linguistic variables and IFNs [38], we develop the transformation relationship between the linguistic variables and PFNs as presented in Table 12. Then, the evaluation information in Table 11 can be transformed into a PF evaluation matrix A = ( a i j ) ( i = 1 , 2 , , 20 ;   j = 1 , 2 , , 9 ) , and the ranking of the 20 web services can be obtained by the Algorithm 1 in this paper. Subsequently, the ranking result will be compared with the rankings determined by AHP, TOPSIS, COPRAS, VIKOR, and SAW methods in [37]. It is worth noting that, in order to compare different MCDM methods, more effectively we suppose each criteria is considered equally important, i.e., w j = 1 / 9 ( j = 1 , 2 , , 9 ) . Then, the ranking of the 20 web services can be determined by the following steps.
  • Step 1: According to the Definition 3, normalize the PF evaluation matrix A = ( a i j ) to the standardized PF evaluation matrix A ¯ = ( a ¯ i j ) as
    a ¯ i j = { a i j , for   the   benefit   criteria ; ( a i j ) c , for   the   cost   criteria .
  • Step 2: Utilize the GPFWA ( λ = 1 ) operator
    G P F W A λ = 1 ( a ¯ i 1 , a ¯ i 2 , , a ¯ i 9 ) = a i = ( 1 j = 1 9 ( 1 μ ¯ i j ) w j , j = 1 9 ( η ¯ i j ) w j , j = 1 9 ( v ¯ i j ) w j )
    to aggregated the PF evaluation matrix A ¯ = ( a ¯ i j ) , and the collective PFNs of each web service are obtained.
  • Step 3: Compute the score values of each web service using the equation
    s ( a i ) = ( 1 + μ i η i v i ) / 2 .
Then, the ranking of the 20 web services can be determined; the lager the score value, the better the web service. The related data of the ranking are presented in Table 13.
To verify the accuracy of the ranking obtained by the proposed method, we use AHP, TOPSIS, COPRAS, VIKOR, and SAW methods to solve the web service selection problem combined with the evaluation information in Table 11. Subsequently, the ranking results of different MCDM methods are presented in Table 14. The Spearman’s rank correlation coefficient is a powerful tool for measuring the similarity between two MCDM methods [39]. Then, we can calculate the Spearman’s rank correlation coefficients between the proposed method and the other five MCDM methods as shown in Table 15. Table 15 shows that the Spearman’s rank correlation coefficients between the proposed method and AHP and TOPSIS are 0.9722 and 0.9549, respectively, which demonstrate that the proposed method is highly correlated with these two methods. AHP and TOPSIS methods have been approved to be the most suitable two methods to solve web service selection problems [39]; thus, the comparison results above illustrate the feasibility of the proposed method.
From the information aggregation of the proposed method, we can find that the calculating procedure of the proposed method is more complicated than AHP and TOPSIS methods. In addition, TOPSIS method does not require the transformation of the evaluation information concerning cost and benefit type criteria. However, when decision makers are not sure if it is 3 or 4 about the evaluation information of the web service W S 1 concerning the criteria C 1 , AHP and TOPSIS methods cannot deal with this situation in practice; we can use PHFS to express the evaluation information above, i.e., {{0.25,0.35},{0.05},{0.55,0.65}}. On the other hand, when the criteria are in different priorities, the GPHFPWA and GPHFPWG operators can be used to aggregation the evaluation information. Thus, the AHP, TOPSIS, and proposed methods have their own advantages and disadvantages; in real life, decision makers can determine to utilize which MCDM methods to solve problems according to the actual situations.

7. Conclusions

Combined with the picture fuzzy set and uniformly typical hesitant fuzzy set, this paper develops the picture hesitant fuzzy set, in which the positive, neutral, negative, and refusal membership degrees are expressed by several possible values. Then, the operations and comparison method of picture hesitant fuzzy elements are developed. To solve the multiple criteria decision-making problems under picture hesitant fuzzy environment, the generalized picture hesitant fuzzy weighted averaging and generalized picture hesitant fuzzy weighted geometric operators are put forward to aggregate the picture hesitant elements given by decision maker. Furthermore, considering the different priorities between the related criteria in practice, the generalized picture hesitant fuzzy prioritized weighted averaging and generalized picture hesitant fuzzy prioritized weighted geometric operators are proposed. Meanwhile, some desirable properties and the reduced operators of them are investigated in detail. Finally, two kinds of multiple criteria decision-making methods combined with the proposed operators are constructed to solve the multiple criteria decision-making problems in different situations. Subsequently, two numerical examples and an application of web service selection are provided to indicate the applications and advantages of the proposed methods.
In future research, we will investigate other operations of picture hesitant fuzzy elements and develop different aggregation operators to aggregate picture hesitant fuzzy elements. In addition, we will propose the consensus model to improve the proposed methods; then, the non-consensus evaluation information of decision makers will be revised to obtain a more accurate ranking result.

Author Contributions

R.W. put forward the picture hesitant fuzzy set, and explored the operational laws and comparison method of the picture hesitant fuzzy set. R.W. and Y.L. developed the aggregation operators under picture hesitant fuzzy environment. Then, R.W. wrote the original manuscript, and Y.L. improved the writing.

Funding

This study was supported by the National Natural Science Foundation of China (no. 71371156) and the Doctoral Innovation Fund Program of Southwest Jiaotong University (D-CX201727).

Acknowledgments

Thank the editor and reviewers for the positive comments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof. 
a. For n = 1 , according to Theorem 1, since
G P H F W A λ ( n ˜ 1 ) = ( w 1 n ˜ 1 λ ) 1 / λ = ( n ˜ 1 λ ) 1 / λ = n ˜ 1 .
Obviously, Equation (22) holds for n = 1 .
b. For n = 2 , since
n ˜ 1 = α 1 μ ˜ 1 , β 1 η ˜ 1 , γ 1 v ˜ 1 { { α 1 λ } , { 1 ( 1 β 1 ) λ } , { 1 ( 1 γ 1 ) λ } } , n ˜ 2 = α 2 μ ˜ 2 , β 2 η ˜ 2 , γ 2 v ˜ 2 { { α 2 λ } , { 1 ( 1 β 2 ) λ } , { 1 ( 1 γ 2 ) λ } } .
we have
w 1 n ˜ 1 λ = α 1 μ ˜ 1 , β 1 η ˜ 1 , γ 1 v ˜ 1 { { 1 ( 1 α 1 λ ) w 1 } , { ( 1 ( 1 β 1 ) λ ) w 1 } , { ( 1 ( 1 γ 1 ) λ ) w 1 } } , w 2 n ˜ 2 λ = α 2 μ ˜ 2 , β 2 η ˜ 2 , γ 2 v ˜ 2 { { 1 ( 1 α 2 λ ) w 2 } , { ( 1 ( 1 β 2 ) λ ) w 2 } , { ( 1 ( 1 γ 2 ) λ ) w 2 } } .
then,
w 1 n ˜ 1 λ w 2 n ˜ 2 λ = α 1 μ ˜ 1 , β 1 η ˜ 1 , γ 1 v ˜ 1 , α 2 μ ˜ 2 , β 2 η ˜ 2 , γ 2 v ˜ 2 { { 1 ( 1 α 1 λ ) w 1 + 1 ( 1 α 2 λ ) w 2 ( 1 ( 1 α 1 λ ) w 1 ) ( 1 ( 1 α 2 λ ) w 2 ) } , { ( 1 ( 1 β 1 ) λ ) w 1 ( 1 ( 1 β 2 ) λ ) w 2 } , { ( 1 ( 1 γ 1 ) λ ) w 1 ( 1 ( 1 γ 2 ) λ ) w 2 } } . w 1 n ˜ 1 λ w 2 n ˜ 2 λ = α 1 μ ˜ 1 , β 1 η ˜ 1 , γ 1 v ˜ 1 , α 2 μ ˜ 2 , β 2 η ˜ 2 , γ 2 v ˜ 2 { { 1 ( 1 α 1 λ ) w 1 ( 1 α 2 λ ) w 2 } , { ( 1 ( 1 β 1 ) λ ) w 1 ( 1 ( 1 β 2 ) λ ) w 2 } , { ( 1 ( 1 γ 1 ) λ ) w 1 ( 1 ( 1 γ 2 ) λ ) w 2 } } .
and
G P H F W A λ ( n ˜ 1 , n ˜ 2 ) = ( w 1 n ˜ 1 λ w 2 n ˜ 2 λ ) 1 / λ = α 1 μ ˜ 1 , β 1 η ˜ 1 , γ 1 v ˜ 1 , α 2 μ ˜ 2 , β 2 η ˜ 2 , γ 2 v ˜ 2 { { ( 1 ( 1 α 1 λ ) w 1 ( 1 α 2 λ ) w 2 ) 1 / λ } , { 1 ( 1 ( 1 ( 1 β 1 ) λ ) w 1 ( 1 ( 1 β 2 ) λ ) w 2 ) 1 / λ } , { 1 ( 1 ( 1 ( 1 γ 1 ) λ ) w 1 ( 1 ( 1 γ 2 ) λ ) w 2 ) 1 / λ } } .
i.e., Equation (22) holds for n = 2 .
c. If Equation (22) holds for n = k , we have
G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ k ) = α 1 μ ˜ 1 , α 2 μ ˜ 2 , , α k μ ˜ k , β 1 η ˜ 1 , β 2 η ˜ 2 , , β k η ˜ k , γ 1 v ˜ 1 , γ 2 v ˜ 2 , , γ k v ˜ k { { ( 1 j = 1 k ( 1 α j λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 k ( 1 ( 1 β j ) λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 k ( 1 ( 1 γ j ) λ ) w j ) 1 / λ } } ,
when n = k + 1 , according to the operations of PHFEs, we have
w 1 n ˜ 1 λ w 2 n ˜ 2 λ w k n ˜ k λ w k + 1 n ˜ k + 1 λ = α 1 μ ˜ 1 , α 2 μ ˜ 2 , , α k μ ˜ k , β 1 η ˜ 1 , β 2 η ˜ 2 , , β k η ˜ k , γ 1 v ˜ 1 , γ 2 v ˜ 2 , , γ k v ˜ k { { 1 j = 1 k ( 1 α j λ ) w j } , { j = 1 k ( 1 ( 1 β j ) λ ) w j } , { j = 1 k ( 1 ( 1 γ j ) λ ) w j } } α k + 1 μ ˜ k + 1 , β k + 1 η ˜ k + 1 , γ k + 1 v ˜ k + 1 { { 1 ( 1 α k + 1 λ ) w k + 1 } , { ( 1 ( 1 β k + 1 ) λ ) w k + 1 } , { ( 1 ( 1 γ k + 1 ) λ ) w k + 1 } } = α 1 μ ˜ 1 , α 2 μ ˜ 2 , , α k + 1 μ ˜ k + 1 , β 1 η ˜ 1 , β 2 η ˜ 2 , , β k + 1 η ˜ k + 1 , γ 1 v ˜ 1 , γ 2 v ˜ 2 , , γ k + 1 v ˜ k + 1 { { 1 j = 1 k + 1 ( 1 α j λ ) w j } , { j = 1 k + 1 ( 1 ( 1 β j ) λ ) w j } , { j = 1 k + 1 ( 1 ( 1 γ j ) λ ) w j } } .
then,
G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ k ) = ( w 1 n ˜ 1 λ w 2 n ˜ 2 λ w k n ˜ k λ w k + 1 n ˜ k + 1 λ ) 1 / λ = α 1 μ ˜ 1 , α 2 μ ˜ 2 , , α k + 1 μ ˜ k + 1 , β 1 η ˜ 1 , β 2 η ˜ 2 , , β k + 1 η ˜ k + 1 , γ 1 v ˜ 1 , γ 2 v ˜ 2 , , γ k + 1 v ˜ k + 1 { { ( 1 j = 1 k + 1 ( 1 α j λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 k + 1 ( 1 ( 1 β j ) λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 k + 1 ( 1 ( 1 γ j ) λ ) w j ) 1 / λ } } .
i.e., Equation (22) holds for n = k + 1 ; we can demonstrate that Equation (22) holds for all values of n . □

Appendix B

Proof. 
According to Theorem 2, since n ˜ j = n ˜ = { μ ˜ , η ˜ , v ˜ } , we have
G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = ( w 1 n ˜ λ w 2 n ˜ λ w n n ˜ λ ) 1 / λ = α μ ˜ , β η ˜ , γ v ˜ { { ( 1 j = 1 n ( 1 α λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 β ) λ ) w j ) 1 / λ } , { 1 ( 1 j = 1 n ( 1 ( 1 γ ) λ ) w j ) 1 / λ } } = α μ ˜ , β η ˜ , γ v ˜ { { ( 1 ( 1 α λ ) ) 1 / λ } , { 1 ( 1 ( 1 ( 1 β ) λ ) ) 1 / λ } , { 1 ( 1 ( 1 ( 1 γ ) λ ) ) 1 / λ } } = α μ ˜ , β η ˜ , γ v ˜ { { ( α λ ) 1 / λ } , { 1 ( ( 1 β ) λ ) 1 / λ } , { 1 ( ( 1 γ ) λ ) 1 / λ } } = α μ ˜ , β η ˜ , γ v ˜ { { α } , { β } , { γ } } = n ˜ .

Appendix C

Proof. 
For λ ( 0 , ) , since α α j α + , then
α j λ ( α ) λ , 1 α j λ 1 ( α ) λ , ( 1 α j λ ) w j ( 1 ( α ) λ ) w j , j = 1 n ( 1 α j λ ) w j j = 1 n ( 1 ( α ) λ ) w j , 1 j = 1 n ( 1 α j λ ) w j 1 j = 1 n ( 1 ( α ) λ ) w j , ( 1 j = 1 n ( 1 α j λ ) w j ) 1 / λ ( 1 j = 1 n ( 1 ( α ) λ ) w j ) 1 / λ = α .
similarly, we have
( 1 j = 1 n ( 1 α j λ ) w j ) 1 / λ ( 1 j = 1 n ( 1 ( α + ) λ ) w j ) 1 / λ = α + .
As β β j β + , then
1 β j 1 β , ( 1 β j ) λ ( 1 β ) λ , 1 ( 1 β j ) λ 1 ( 1 β ) λ , ( 1 ( 1 β j ) λ ) w j ( 1 ( 1 β ) λ ) w j , j = 1 n ( 1 ( 1 β j ) λ ) w j j = 1 n ( 1 ( 1 β ) λ ) w j , 1 j = 1 n ( 1 ( 1 β j ) λ ) w j 1 j = 1 n ( 1 ( 1 β ) λ ) w j , ( 1 j = 1 n ( 1 ( 1 β j ) λ ) w j ) 1 / λ ( 1 j = 1 n ( 1 ( 1 β ) λ ) w j ) 1 / λ , 1 ( 1 j = 1 n ( 1 ( 1 β j ) λ ) w j ) 1 / λ 1 ( 1 j = 1 n ( 1 ( 1 β ) λ ) w j ) 1 / λ = β .
similarly, we have
1 ( 1 j = 1 n ( 1 ( 1 β j ) λ ) w j ) 1 / λ 1 ( 1 j = 1 n ( 1 ( 1 β + ) λ ) w j ) 1 / λ = β + .
and, as γ γ j γ + , we have
γ 1 ( 1 j = 1 n ( 1 ( 1 γ j ) λ ) w j ) 1 / λ γ + .
let G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) = n ˜ = { { α } , { β } , { γ } } , then
s ( n ˜ ) = 1 + 1 l i = 1 l α i 1 p i = 1 p β i 1 q i = 1 q γ i 2 1 + 1 l i = 1 l α 1 p i = 1 p β + 1 q i = 1 q γ + 2 = s ( n ˜ ) ,
and
s ( n ˜ ) = 1 + 1 l i = 1 l α i 1 p i = 1 p β i 1 q i = 1 q γ i 2 1 + 1 l + i = 1 l + α + 1 p + i = 1 p + β 1 q + i = 1 q + γ 2 = s ( n ˜ + ) .  
we obtain
n ˜ G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) n ˜ + , λ ( 0 , ) .
Similarly, we have
n ˜ G P H F W A λ ( n ˜ 1 , n ˜ 2 , , n ˜ n ) n ˜ + , λ ( , 0 ) .

References

  1. Mardani, A.; Jusoh, A.; Zavadskas, E.K. Fuzzy multiple criteria decision-making techniques and applications–two decades review from 1994 to 2014. Expert Syst. Appl. 2015, 42, 4126–4148. [Google Scholar] [CrossRef]
  2. Morente-Molinera, J.A.; Kou, G.; González-Crespo, R.; Corchado, J.M.; Herrera-Viedma, E. Solving multi-criteria group decision making problems under environments with a high number of alternatives using fuzzy ontologies and multi-granular linguistic modelling methods. Knowl. Syst. 2017, 137, 54–64. [Google Scholar] [CrossRef]
  3. Hernández, F.L.; Ory, E.G.D.; Aguilar, S.R.; González-Crespo, R. Residue properties for the arithmetical estimation of the image quantization table. Appl. Soft Comput. 2018, 67, 309–321. [Google Scholar] [CrossRef]
  4. Zadeh, L.A. Fuzzy sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef] [Green Version]
  5. Atanassov, K.T. Intuitionistic fuzzy sets. Fuzzy sets Syst. 1986, 20, 87–96. [Google Scholar] [CrossRef]
  6. Atanassov, K.T. Gargov, G. Interval valued intuitionistic fuzzy sets. Fuzzy Sets Syst. 1989, 31, 343–349. [Google Scholar] [CrossRef]
  7. Cường, B.C. Picture fuzzy sets. J. Comput. Sci.Cybern. 2015, 30, 409–420. [Google Scholar] [CrossRef]
  8. Torra, V. Hesitant fuzzy sets. Int. J. Intell. Syst. 2010, 25, 529–539. [Google Scholar] [CrossRef]
  9. Ye, J. Trapezoidal neutrosophic set and its application to multiple attribute decision-making. Neural Comput. Appl. 2015, 26, 1157–1166. [Google Scholar] [CrossRef]
  10. Zhu, B.; Xu, Z.S.; Xia, M.M. Dual hesitant fuzzy sets. J. Appl. Math. 2012, 2012, 2607–2645. [Google Scholar] [CrossRef]
  11. Liu, F.; Yuan, X.H. Fuzzy number intuitionistic fuzzy set. Fuzzy Syst. Math. 2007, 21, 88–91. [Google Scholar]
  12. Garg, H. Some picture fuzzy aggregation operators and their applications to multicriteria decision-making. Arab. J. Sci. Eng. 2017, 42, 1–16. [Google Scholar] [CrossRef]
  13. Cuong, B. Picture fuzzy sets-first results. Part 1. In Proceedings of the Third World Congress on Information and Communication WICT’2013, Hanoi, Vietnam, 15–18 December 2013; pp. 1–6. [Google Scholar]
  14. Singh, P. Correlation coefficients for picture fuzzy sets. J. Intell. Fuzzy Syst. 2015, 28, 1–12. [Google Scholar]
  15. Son, L.H. Generalized picture distance measure and applications to picture fuzzy clustering. Appl. Soft Comput. 2016, 46, 284–295. [Google Scholar] [CrossRef]
  16. Wei, G.W. Picture fuzzy cross-entropy for multiple attribute decision making problems. J. Bus. Econ. Manag. 2016, 17, 491–502. [Google Scholar] [CrossRef]
  17. Torra, V.; Narukawa, Y. On hesitant fuzzy sets and decision. In Proceedings of the IEEE International Conference on Fuzzy Systems, Jeju Island, Korea, 20–24 August 2009; pp. 1378–1382. [Google Scholar]
  18. Chen, N.; Xu, Z.S.; Xia, M.M. Interval-valued hesitant preference relations and their applications to group decision making. Knowl. Syst. 2013, 37, 528–540. [Google Scholar] [CrossRef]
  19. Farhadinia, B. Correlation for dual hesitant fuzzy sets and dual interval-valued hesitant fuzzy sets. Int. J. Intell. Syst. 2013, 29, 184–205. [Google Scholar] [CrossRef]
  20. Xu, Z.S.; Yager, R.R. Some geometric aggregation operators based on intuitionistic fuzzy sets. Int. J. Intell. Syst. 2006, 35, 417–433. [Google Scholar] [CrossRef]
  21. Xu, Z.S. Intuitionistic fuzzy aggregation operators. IEEE Trans. Fuzzy Syst. 2007, 15, 1179–1187. [Google Scholar]
  22. Xu, Z.S. Methods for aggregating interval-valued intuitionistic fuzzy information and their application to decision making. Control Decis. 2007, 22, 215–219. [Google Scholar]
  23. Xu, Z.S.; Chen, J. Approach to group decision making based on interval-valued intuitionistic judgment matrices. Syst. Eng. Theory Pract. 2007, 27, 126–133. [Google Scholar] [CrossRef]
  24. Xu, Z.S.; Chen, J. On geometric aggregation over interval-valued intuitionistic fuzzy information. In Proceedings of the International Conference on Fuzzy Systems and Knowledge Discovery, Haikou, China, 24–27 August 2007; pp. 466–471. [Google Scholar]
  25. Wei, G.W. Picture fuzzy aggregation operators and their application to multiple attribute decision making. J. Intell. Fuzzy Syst. 2017, 33, 713–724. [Google Scholar] [CrossRef]
  26. Wang, C.; Zhou, X.; Tu, H.; Tao, S. Some geometric aggregation operators based on picture fuzzy sets and their application in multiple attribute decision making. Ital. J. Pure Appl. Math. 2017, 37, 477–492. [Google Scholar]
  27. Xia, M.M.; Xu, Z.S. Hesitant fuzzy information aggregation in decision making. Int. J. Approx. Reason. 2011, 52, 395–407. [Google Scholar] [CrossRef] [Green Version]
  28. Yu, D.J.; Zhang, W.Y.; Huang, G. Dual hesitant fuzzy aggregation operators. Technol. Econ. Develop. Econ. 2016, 22, 194–209. [Google Scholar] [CrossRef]
  29. Wang, H.J.; Zhao, X.F.; Wei, G.W. Dual hesitant fuzzy aggregation operators in multiple attribute decision making. J. Intell. Fuzzy Syst. Appl. Eng. Technol. 2014, 26, 2281–2290. [Google Scholar]
  30. Zhang, W.K.; Li, X.; Ju, Y.B. Some aggregation operators based on Einstein operations under interval-valued dual hesitant fuzzy setting and their application. Math. Probl. Eng. 2014. [Google Scholar] [CrossRef]
  31. Yager, R.R. Prioritized aggregation operators. Int. J. Approx. Reason. 2008, 48, 263–274. [Google Scholar] [CrossRef]
  32. Yu, D.J. Intuitionistic fuzzy prioritized operators and their application in multi-criteria group decision making. Technol. Econ. Develop. Econ. 2013, 19, 1–21. [Google Scholar] [CrossRef]
  33. Yu, D.J.; Wu, Y.Y.; Lu, T. Interval-valued intuitionistic fuzzy prioritized operators and their application in group decision making. Knowl. Syst. 2012, 30, 57–66. [Google Scholar] [CrossRef]
  34. Wei, G.W. Hesitant fuzzy prioritized operators and their application to multiple attribute decision making. Knowl. Syst. 2012, 31, 176–182. [Google Scholar] [CrossRef]
  35. Bedregal, B.; Reiser, R.; Bustince, H.; Lopez-Molina, C.; Torra, V. Aggregation functions for typical hesitant fuzzy elements and the action of automorphisms. Infor. Sci. 2014, 255, 82–99. [Google Scholar] [CrossRef] [Green Version]
  36. Alcantud, J.C.R.; Torra, V. Decomposition theorems and extension principles for hesitant fuzzysets. Infor. Fusion. 2018, 41, 48–56. [Google Scholar] [CrossRef]
  37. Bagga, P.; Joshi, A.; Hans, R. QoS based web service selection and multi-criteria decision making methods. Int. J. Int. Multimed. Artif. Intell. 2017. [Google Scholar] [CrossRef]
  38. Zhao, H.; You, J.X.; Liu, H.C. Failure mode and effect analysis using MULTIMOORA method with continuous weighted entropy under interval-valued intuitionistic fuzzy environment. Soft Comput. 2017, 12, 5355–5367. [Google Scholar] [CrossRef]
  39. Kou, G.; Lu, Y.Q.; Peng, Y.; Shi, Y. Evaluation of classification algorithms using MCDM and rank correlation. Int. J. Inf. Technol. Decis. Mak. 2012, 11, 197–225. [Google Scholar] [CrossRef]
Figure 1. Flow diagram of the Algorithm 1.
Figure 1. Flow diagram of the Algorithm 1.
Symmetry 10 00295 g001
Figure 2. Flow diagram of the Algorithm 2.
Figure 2. Flow diagram of the Algorithm 2.
Symmetry 10 00295 g002
Table 1. PHF evaluation matrix of Example 1.
Table 1. PHF evaluation matrix of Example 1.
Alternatives C 1 C 2 C 3 C 4
A 1 {{0.43,0.53},{0.33},
{0.06,0.09}}
{{0.76,0.89},
{0.05,0.08},{0.03}}
{{0.42},{0.35},
{0.12,0.18}}
{{0.08},{0.75,0.89},
{0.02}}
A 2 {{0.53,0.65,0.73},
{0.10,0.12},{0.08}}
{{0.13},{0.53,0.64},
{0.12,0.21}}
{{0.03},{0.77,0.82},
{0.10,0.13}}
{{0.58,0.73},{0.15},
{0.08}}
A 3 {{0.72,0.86,0.91},
{0.03},{0.02}}
{{0.07},{0.05,0.09},
{0.05}}
{{0.04},{0.65,0.72,0.85},{0.05,0.10}}{{0.45,0.68},{0.18,0.26},{0.06}}
A 4 {{0.77,0.85},{0.09},
{0.05}}
{{0.65,0.74},{0.10,0.16},{0.10}}{{0.02},{0.78,0.89},{0.05}}{{0.08},{0.65,0.84},
{0.06}}
A 5 {{0.70,0.81,0.90},
{0.05},{0.02}}
{{0.68},{0.08},
{0.13,0.21}}
{{0.05},{0.77,0.87},{0.06}}{{0.13},{0.65,0.75},
{0.09}}
Table 2. PHF evaluation matrix in Example 2.
Table 2. PHF evaluation matrix in Example 2.
Alternatives C 1 C 2 C 3 C 4
A 1 {{0.40,0.50,0.70},{0.05},{0.10,0.20}}{{0.65},{0.05,0.08},
{0.15}}
{{0.40,0.50,0.60},{0.03},{0.10,0.20}}{{0.55},{0.10,0.15},
{0.15}}
A 2 {{0.65,0.75},{0.02,0.04},{0.15}}{{0.60},{0.05,0.10},
{0.10,0.20}}
{{0.75,0.80},{0.06},
{0.05,0.08}}
{{0.40,0.50},{0.20},
{0.15,0.25}}
A 3 {{0.70},{0.06,0.10},
{0.10,0.15}}
{{0.20,0.30,0.50},{0.04},{0.30,0.40}}{{0.50},{0.03,0.06},
{0.30,0.35}}
{{0.50,0.70},{0.10},
{0.10}}
A 4 {{0.50,0.60,0.70},{0.08},{0.10}}{{0.40,0.50},{0.20},
{0.10,0.20}}
{{0.85},{0.03,0.07},
{0.05}}
{{0.45},{0.10,0.20},
{0.15,0.30}}
A 5 {{0.65},{0.05,0.10},
{0.15,0.20}}
{{0.50,0.70},{0.08},
{0.20}}
{{0.70,0.80},{0.04},
{0.10}}
{{0.35},{0.10,0.20},{0.30,0.40}}
Table 3. Sensitivity analysis results obtained by the GPHFWA operator.
Table 3. Sensitivity analysis results obtained by the GPHFWA operator.
Values s ( n ˜ 1 ) s ( n ˜ 2 ) s ( n ˜ 3 ) s ( n ˜ 4 ) s ( n ˜ 5 ) Ranking
λ = 0.001 0.41170.49880.57870.33880.4058 A 3 A 2 A 1 A 5 A 4
λ = 0.5 0.44120.53280.61240.39770.4526 A 3 A 2 A 5 A 1 A 4
λ = 1 0.47010.56110.64090.45530.4989 A 3 A 2 A 5 A 1 A 4
λ = 2 0.51970.60010.68120.54260.5743 A 3 A 2 A 5 A 4 A 1
λ = 3 0.55780.62420.70740.59810.6255 A 3 A 5 A 2 A 4 A 1
λ = 5 0.61330.65220.74060.66130.6852 A 3 A 5 A 4 A 2 A 1
λ = 10 0.69650.68310.78370.72740.7481 A 3 A 5 A 4 A 1 A 2
λ = 20 0.76790.70620.81970.77220.7916 A 3 A 5 A 4 A 1 A 2
λ = 50 0.82430.72740.85270.80630.8280 A 3 A 5 A 1 A 4 A 2
Table 4. Sensitivity analysis results obtained by the GPHFWG operator.
Table 4. Sensitivity analysis results obtained by the GPHFWG operator.
Values s ( n ˜ 1 ) s ( n ˜ 2 ) s ( n ˜ 3 ) s ( n ˜ 4 ) s ( n ˜ 5 ) Ranking
λ = 0.001 0.32270.37910.44860.21420.2720 A 3 A 2 A 1 A 5 A 4
λ = 0.5 0.30190.34920.41200.19600.2457 A 3 A 2 A 1 A 5 A 4
λ = 1 0.28130.31970.37780.18080.2240 A 3 A 2 A 1 A 5 A 4
λ = 2 0.24510.26900.32370.15870.1926 A 3 A 2 A 1 A 5 A 4
λ = 3 0.21680.23140.28580.14390.1716 A 3 A 2 A 1 A 5 A 4
λ = 5 0.17800.18290.23890.12490.1441 A 3 A 2 A 1 A 5 A 4
λ = 10 0.13010.12650.18780.09910.1093 A 3 A 1 A 2 A 5 A 4
λ = 20 0.12520.11580.19230.12160.1095 A 3 A 1 A 4 A 2 A 5
λ = 50 0.14580.13770.16920.10410.1351 A 3 A 1 A 2 A 5 A 4
Table 5. Sensitivity analysis results obtained by the GPHFPWA operator.
Table 5. Sensitivity analysis results obtained by the GPHFPWA operator.
Values s ( n ˜ 1 ) s ( n ˜ 2 ) s ( n ˜ 3 ) s ( n ˜ 4 ) s ( n ˜ 5 ) Ranking
λ = 0.001 0.68660.73370.64770.69110.6829 A 2 A 4 A 1 A 5 A 3
λ = 0.5 0.68770.73520.65290.69440.6852 A 2 A 4 A 1 A 5 A 3
λ = 1 0.68880.73680.65800.69780.6874 A 2 A 4 A 1 A 5 A 3
λ = 2 0.69130.73990.66800.70540.6920 A 2 A 4 A 5 A 1 A 3
λ = 3 0.69400.74320.67690.71350.6964 A 2 A 4 A 5 A 1 A 3
λ = 5 0.69940.74950.69150.72980.7045 A 2 A 4 A 5 A 1 A 3
λ = 10 0.71080.76340.71510.76280.7203 A 2 A 4 A 5 A 3 A 1
λ = 20 0.72430.78310.73690.79760.7407 A 4 A 2 A 5 A 3 A 1
λ = 50 0.74080.80920.75700.83320.7693 A 4 A 2 A 5 A 3 A 1
Table 6. Sensitivity analysis results obtained by the GPHFPWG operator.
Table 6. Sensitivity analysis results obtained by the GPHFPWG operator.
Values s ( n ˜ 1 ) s ( n ˜ 2 ) s ( n ˜ 3 ) s ( n ˜ 4 ) s ( n ˜ 5 ) Ranking
λ = 0.001 0.68140.72300.62420.67330.6706 A 2 A 1 A 4 A 5 A 3
λ = 0.5 0.67890.71630.61470.66470.6646 A 2 A 1 A 4 A 5 A 3
λ = 1 0.67570.70800.60400.65500.6572 A 2 A 1 A 5 A 4 A 3
λ = 2 0.66800.68850.58150.63440.6393 A 2 A 1 A 5 A 4 A 3
λ = 3 0.65970.66920.56110.61550.6194 A 2 A 1 A 5 A 4 A 3
λ = 5 0.64420.63780.53060.58730.5815 A 1 A 2 A 4 A 5 A 3
λ = 10 0.61900.59310.49210.55240.5213 A 1 A 2 A 4 A 5 A 3
λ = 20 0.67270.57740.51220.54510.5002 A 1 A 2 A 4 A 3 A 5
λ = 50 0.78250.73550.67320.72140.6872 A 1 A 2 A 4 A 5 A 3
Table 7. Comparison result of Example 1.
Table 7. Comparison result of Example 1.
MCDM methodRanking
The GPHFWA operator ( λ = 1 ) A 3 A 2 A 5 A 1 A 4
The GPHFWG operator ( λ = 1 ) A 3 A 2 A 1 A 5 A 4
The PFWA operator A 3 A 2 A 1 A 5 A 4
The PFWG operator A 3 A 1 A 2 A 5 A 4
The picture fuzzy cross-entropy method A 3 A 1 A 2 A 5 A 4
Table 8. Comparison result of Example 2.
Table 8. Comparison result of Example 2.
MCDM methodRanking
The GPHFPWA operator ( λ = 1 ) A 2 A 4 A 1 A 5 A 3
The GPHFPWG operator ( λ = 1 ) A 2 A 1 A 5 A 4 A 3
The HFPWA operator A 5 A 2 A 1 A 4 A 3
The HFPWG operator A 2 A 5 A 1 A 4 A 3
Table 9. Comparison of each MCDM methods.
Table 9. Comparison of each MCDM methods.
MethodsAdvantagesDisadvantages
The GPHFWA/GPHFWG operator• The human opinions including yes, abstain, no, and refusal can be expressed, and each membership functions can be represented by several possible values.
• The operators can be reduced to other forms by varying the value of λ .
• The PHFS can be transformed into its special cases, i.e., FS, IFS, PFS, and UTHFS.
• The calculating process is complex when the numbers of criteria are relatively large.
• The size of data is relatively large.
The PFWA/PFWG operator• The human opinions including yes, abstain, no, and refusal can be expressed.
• The PFS can be transformed into IFS and FS.
• It cannot express the evaluation information when decision makers have difficulty determining an accurate value of each membership level.
Picture fuzzy cross-entropy• The human opinions including yes, abstain, no, and refusal can be expressed.
• The PFS can be transformed into IFS and FS.
• The ranking is obtained without aggregating the evaluation information; it can avoid the loss of information.
• The step of normalizing the evaluation information can be omitted.
• It cannot express the evaluation information when decision makers have difficulty determining an accurate value of each membership level.
• It cannot solve the multiple criteria group decision-making problems.
The GPHFPWA/GPHFPWG operator• The human opinions including yes, abstain, no, and refusal can be expressed, and each membership functions can be represented by several possible values.
• The operators can be reduced to other forms by varying the value of the λ .
• The PHFS can be transformed into its special cases, i.e., FS, IFS, PFS, and UTHFS.
• It can solve the MCDM problem that the criteria are in different priorities.
• The calculating process is complex when the numbers of criteria are relatively large.
• The size of data is relatively large.
The HFPWA/HFPWG operator• The positive membership function can be expressed by several possible values.
• It can solve the MCDM problem that the criteria are in different priorities.
• It cannot express the human opinions including abstain, no, and refusal.
Table 10. MCDM application fields of each MCDM method.
Table 10. MCDM application fields of each MCDM method.
MethodsMCDM Application Fields
The GPHFWA/GPHFWG operator• The evaluation information of decision makers is diverse.
• Decision makers feel doubtful to determine the accurate value of each membership level.
• The numbers of the criteria are relatively small.
The PFWA/PFWG operator• The evaluation information of decision makers is diverse.
Picture fuzzy cross entropy• The evaluation information of decision maker is diverse.
• The alternatives are evaluated by an individual decision maker.
The GPHFPWA/GPHFPWG operator• The evaluation information of decision makers is diverse.
• Decision makers feel doubtful to determine the accurate value of each membership level.
• The numbers of the criteria are relatively small.
• The criteria are in different priorities.
The HFPWA/HFPWG operator• Decision makers feel doubtful to determine the accurate value of positive membership level.
• The criteria are in different priorities.
Table 11. Evaluation information of each web service.
Table 11. Evaluation information of each web service.
Alternatives C 1 C 2 C 3 C 4 C 5 C 6 C 7 C 8 C 9
W S 1 342567628
W S 2 452346834
W S 3 356832245
W S 4 445562678
W S 5 562477883
W S 6 443768686
W S 7 445787743
W S 8 676657863
W S 9 532265247
W S 10 887622343
W S 11 584457458
W S 12 455568766
W S 13 647645445
W S 14 465454674
W S 15 756276552
W S 16 364723824
W S 17 448445248
W S 18 435354466
W S 19 544667672
W S 20 623546555
Table 12. Transformation between linguistic variables and PFNs.
Table 12. Transformation between linguistic variables and PFNs.
Crisp numbersLinguistic variablesPFNs
1Extremely low (EL)(0.05,0.05,0.85)
2Very low (VL)(0.15,0.05,0.75)
3Low (L)(0.25,0.05,0.65)
4Medium low (ML)(0.35,0.05,0.55)
5Medium (M)(0.45,0.05,0.45)
6Medium high (MH)(0.55,0.05,0.35)
7High (H)(0.65,0.05,0.25)
8Very high (VH)(0.75,0.05,0.15)
9Extremely high (EH)(0.85,0.05,0.05)
Table 13. Ranking results obtained by the proposed method.
Table 13. Ranking results obtained by the proposed method.
AlternativesCollective Evaluation InformationScore ValuesRanking
W S 1 (0.4677,0.0500,0.4210)0.498312
W S 2 (0.4833,0.0500,0.4068)0.51339
W S 3 (0.4313,0.0500,0.4579)0.461715
W S 4 (0.3776,0.0500,0.5192)0.404218
W S 5 (0.5244,0.0500,0.3647)0.55485
W S 6 (0.4736,0.0500,0.4159)0.503811
W S 7 (0.5817,0.0500,0.3119)0.60992
W S 8 (0.5871,0.0500,0.3079)0.61461
W S 9 (0.3485,0.0500,0.5475)0.375520
W S 10 (0.5447,0.0500,0.3414)0.57673
W S 11 (0.4683,0.0500,0.4222)0.498013
W S 12 (0.5045,0.0500,0.3879)0.53338
W S 13 (0.4828,0.0500,0.4145)0.509110
W S 14 (0.4371,0.0500,0.4609)0.463114
W S 15 (0.5426,0.0500,0.3498)0.57144
W S 16 (0.5191,0.0500,0.3668)0.55116
W S 17 (0.4154,0.0500,0.4744)0.445516
W S 18 (0.3535,0.0500,0.5459)0.378819
W S 19 (0.5186,0.0500,0.3737)0.54747
W S 20 (0.4179,0.0500,0.4798)0.444017
Table 14. Ranking of web services of different MCDM methods.
Table 14. Ranking of web services of different MCDM methods.
AlternativesThe Proposed MethodAHPTOPSISCOPRASVIKORSAW
W S 1 121314121514
W S 2 9101161212
W S 3 151615111719
W S 4 181818201918
W S 5 568999
W S 6 111112161413
W S 7 222132
W S 8 111511
W S 9 202020192020
W S 10 3444810
W S 11 13111015138
W S 12 8671046
W S 13 1086825
W S 14 1413131367
W S 15 433373
W S 16 68921111
W S 17 161716171816
W S 18 191919181017
W S 19 745754
W S 20 171517141615
Table 15. Spearman’s rank correlation coefficients between the proposed method and the other MCDM methods.
Table 15. Spearman’s rank correlation coefficients between the proposed method and the other MCDM methods.
Existing MCDM MethodsSpearman’s Rank Correlation Coefficients
AHP0.9722
TOPSIS0.9549
COPRAS0.9023
VIKOR0.7429
SAW0.8165

Share and Cite

MDPI and ACS Style

Wang, R.; Li, Y. Picture Hesitant Fuzzy Set and Its Application to Multiple Criteria Decision-Making. Symmetry 2018, 10, 295. https://doi.org/10.3390/sym10070295

AMA Style

Wang R, Li Y. Picture Hesitant Fuzzy Set and Its Application to Multiple Criteria Decision-Making. Symmetry. 2018; 10(7):295. https://doi.org/10.3390/sym10070295

Chicago/Turabian Style

Wang, Rui, and Yanlai Li. 2018. "Picture Hesitant Fuzzy Set and Its Application to Multiple Criteria Decision-Making" Symmetry 10, no. 7: 295. https://doi.org/10.3390/sym10070295

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop