TOPSIS Methods for Probabilistic Hesitant Fuzzy MAGDM and Application to Performance Evaluation of Public Charging Service Quality

. The performance evaluation of public charging service quality is frequently viewed as the multiple attribute group decision-making (MAGDM) issue. In this paper, an extended TOPSIS model is established to provide new means to solve the performance evaluation of public charging service quality. The TOPSIS method integrated with FUCOM method in probabilistic hesitant fuzzy circumstance is applied to rank the optional alternatives and a numerical example for performance evaluation of public charging service quality is used to test the newly proposed method’s practicability in comparison with other methods. The results display that the approach is uncomplicated, valid and simple to compute. The main results of this paper: (1) a novel PHF-TOPSIS method is proposed; (2) the extended TOPSIS method is developed in the probabilistic hesitant fuzzy environment; (3) the FUCOM method is used to obtain the attribute weight; (4) the normalization process of the original data has adapted the latest method to verify the precision; (5) The built models and methods are useful for other selection issues and evaluation issues.


Introduction
Many management decision-making problems in the real world, such as logistics park location, supplier selection, medical service evaluation, fault diagnosis, etc., can be considered from the perspective of MADM (Garg et al., 2018;Akram et al., 2021;Waseem et al., 2019;Lu et al., 2021;Wei et al., 2022).As an important branch of management science and modern decision science, MADM theory and methods have been widely used in many practical decision-making problems (Yang and Pang, 2019;Xu and Zhang, 2019;Zavadskas et al., 2013;Ning et al., 2022).In the actual MAGDM process, due to the complexity and uncertainty of objective things, the limitations of human cognition and the ambiguity of thinking, it is difficult to use quantitative and accurate information to describe decision objects (Wang et al., 2022;Zhang H. et al., 2022;Liu et al., 2019;Li et al., 2021).In 1965, Zadeh (1965) was the first to define a novel fuzzy sets (FSs) to cope with information in the fuzzy new domain (Garg and Kumar, 2018;Garg, 2018b;Zhang and Xu, 2015;Su et al., 2022;Jiang et al., 2022;Lei et al., 2022).To extend the FSs, the intuitionistic fuzzy sets (IFSs) (Atanassov, 1989) were also defined.Subsequently, FSs and its related extension knowledges were predominantly exploited in decision analysis domains (Yu et al., 2017;Wan and Li, 2014;Zhang D. et al., 2022;Zhang et al., 2022a).Su et al. (2011) proposed the interactive method for dynamic IF-MAGDM.Arya and Yadav (2018) defined the intuitionistic fuzzy super-efficiency slack-based measure.Tian et al. (2017) studied the partial derivative and complete differential of binary IF-mathematical functions.Garg (2018a) proposed the improved cosine similarity measure for IFSs.Tan (2011) constructed the Choquet integral-based TOPSIS method for IF-MADM.Zhao et al. (2017) defined the Interactive intuitionistic fuzzy algorithms for multilevel programming problems.Li (2011) built the GOWA operator to MADM using IFSs.Buyukozkan et al. (2018) selected the transportation schemes with integrated intuitionistic fuzzy Choquet integral method.Joshi et al. (2018) defined the Jensen-alpha-Norm dissimilarity measure for IFSs.De and Sana (2018) defined the The (p, q, r, l) method for random demand with Bonferroni mean under IFSs.Li et al. (2018) defined the time-preference and VIKOR-based dynamic method for IF-MADM.Niroomand (2018) defined the multiobjective based direct solution method for linear programming along with intuitionistic fuzzy parameters.Zhao et al. (2021) perfected TODIM for IF-MAGDM on the strength of cumulative prospect theory.Yu S. et al. (2017) defined the derivatives and differentials for multiplicative IFSs.Yu (2012) defined the generalized prioritized geometric operators under IFSs.Xiao et al. (2020) built the intuitionistic fuzzy taxonomy method.Wu and Zhang (2011) built the IF-MADM based on weighted entropy.Verma and Sharma (2014) defined the measure of inaccuracy IF-MADM.Iakovidis and Papageorgiou (2011) defined the cognitive maps for medical decision making under IFSs.
Then the hesitant fuzzy element (HFE) proposed by Xia and Xu (2011) is to solve the problem of determining the element's membership to a set on account of the uncertainty between different numbers and then prove the intuitionistic fuzzy set and hesitant fuzzy set.With the proposition of the HFE, the idea of correspondent operators to aggregate hesitant fuzzy information was obtained.Not long after this, Xu and Xia (2011) raised the idea of the score function, deviation function and the comparison rule, and set the basis on the calculation.Xu and Cai (2010) provided the aggregating operators to integrate the hesitant fuzzy information.Nevertheless, HFE can be regarded as a particular equivalent form whose occurring probabilities of the possible value are equal.The probabilistic hesitant fuzzy set and the corresponding score function, deviation function and its comparison law were proposed by Xu and Zhou (2017).Moreover, the probabilistic hesitant fuzzy weighted averaging geometric operators were introduced by Xu and Zhou (2017) to process PHFE information.Then the improved PHFS was introduced by Zhang et al. (2017) to give more space for hesitation, the integrations of the improved PHES can be calculated by the improved operators.Farhadinia and Xu (2021) gave the comparison techniques of PHFEs.Krishankumar et al. (2021a) built the COPRAS approach to PHFSs.Krishankumar et al. (2021b) proposed to extend a well-known VIKOR method to the PHFS context.Lin et al. (2021) put forward a novel probabilistic hesitant fuzzy MULTIMOORA method.Liu et al. (2021) defined the DEA cross-efficiency with probabilistic hesitant fuzzy preference relations.Yang and Xu (2021) defined the measure of probabilistic hesitant I-fuzzy sets and decision making for strategy choice.Song and Chen (2021) extended the COPRAS method to solve MADM problems under probabilistic hesitant fuzzy environment.Liu and Guan (2021) devised a new PHFE comparison method and then defined the comprehensive characteristic distance measure based on four characteristics.
Technique for order performance by similarity to ideal solution (TOPSIS) was originally developed by Hwang and Yoon (1981) for the sake of addressing a MADM problem.Lai et al. (1994) expanded TOPSIS to deal with a diverse objective decision making problem for Bow River Valley water quality management.Chen (2000) proposed the TOPSIS approach for group decision-making within the fuzzy environment.Wang and Elhag (2006) employed the fuzzy TOPSIS approach on the basis of alpha level sets for bridge risk assessment.Taleizadeh et al. (2009) designed a novel method which combined Pareto, TOPSIS and genetic algorithm to solve the multi-product multi-constraint inventory control systems with random fuzzy replenishments.Zhang et al. (2022b) defined the TOPSIS method for spherical fuzzy MAGDM based on cumulative prospect theory.Wei (2010) developed the TOPSIS method to cope with 2-tuple linguistic MAGDM with incomplete weight information.Nilashi et al. (2019) used two MADM techniques, Decision Making Trial and Evaluation Laboratory (DEMATEL) and Fuzzy TOPSIS, to reveal the interrelationships among the factors and to find the relative importance of these factors in the decision making model.In this paper, we extend the TOPSIS method to probabilistic hesitant fuzzy (PHF) environment based on the FUCOM method to deal with the flexible and complicated decision-making circumstance.The following is the innovation of this paper: (1) a novel PHF-TOPSIS method is proposed; (2) the extended TOPSIS method is developed in the probabilistic hesitant fuzzy environment; (3) the FUCOM method is used to obtain the attribute weight; (4) the normalization process of the original data has adapted the latest method to verify the precision.
The whole thread of the article is as follows: Section 2 gives a simple introduction of the PHF information, Section 3 structures the model of TOPSIS and Section 4 illustrates an example for performance evaluation of public charging service quality to prove the practicability of this new method.Section 5 gives a sensitivity analysis and comparison analysis with other existing models.

Preliminaries
Definition 1 (Xu and Zhou, 2017).Assume q is a fixed set, and probabilistic hesitant fuzzy sets on E, which range from 0 to 1, and the probabilistic hesitant fuzzy element (PHFE) is described as follows: V q = v q t a q g a q t a q , g a q , (1) where t a q ∈ R, 0 t a q 1, a = 1, 2, . . ., #t, and #t represents the total number of elements, t a q shows the degree of membership, while g a q is the probability of the membership degree, #t a=1 t a q = 1.
The first and significant step is the normalization process and we adapt the normalization approach proposed by Li et al. (2020) to break the limitation when processing multiplication of the sets which include different probabilities.
Definition 2 (Li et al., 2020) Step 1. Define the first element.If Definition 3 (Xu and Zhou, 2017).Calculate the score function = by Eq. ( 2): where ti shows the i-th largest elements of normalized PHFE, and ḡi is the probability of occurrence of the corresponding element.

PHF-TOPSIS Method for MAGDM
The MAGDM decision matrix is V c = [v c rx (g rx )] s×y , and the optional alternatives are defined as Q r = {Q 1 , Q 2 , . . ., Q s } and the attribute is shown as I x = {I 1 , I 2 , . . ., I y } and the decision makers (DMs) are defined as c = {c 1 , c 2 , . . ., c l }, while the weighting vector between the DMs is defined as u c = {0.3,0.3, 0.4} and the weighting vector among the criterions is j x which is unknown, l c=1 u c = 1, With the above conventions, the operation of the PHF-TOPSIS is as follows: the whole operation flow chart is shown in Fig. 1.
Then using the introduction in Definition 2 to get normalized matrices.
Step 2. Acquire the criterion weights px using FUCOM method.
The full consistency (FUCOM) method proposed by Pamučar et al. (2018) is the latest model for weighting the coefficients of attributes.Compared to other methods, such as the best worst method (BWM) and the analytic hierarchy process (AHP) method, the FUCOM method can give more appropriate results when considering the relation between the criterions and the number of comparisons (only y − 1).The application of this new model for criterion weight has been used in numerous fields, such as assessment of alternative fuel vehicles for sustainable road transportation (Pamucar et al., 2021), safety evaluation of road sections (Simić et al., 2020), sustained academic quality assurance and ABET accreditation (Ahmad and Qahmash, 2020) and so on.Bozanic et al. (2020) built the MADM with Z-numbers based on FUCOM and MABAC.Durmić et al. (2020) combined FUCOM-Rough SAW model.Bozanić et al. (2021) built the FUCOM-Fuzzy RAFSI model for selecting the group of construction machines for enabling mobility with D-numbers.Pamucar and Ecer (2020) prioritized the weights of the evaluation criteria under fuzziness through the fuzzy full consistency method-FUCOM-F.Stević and Brković (2020) built novel integrated FUCOM-MARCOS model for evaluation of human resources in a transport company.Simić et al. (2020) built the CRITIC-fuzzy FUCOM-DEA-fuzzy MARCOS model for safety evaluation of road sections based on geometric parameters of road.Fazlollahtabar et al. (2019) defined the FUCOM method in group decision-making for selection of forklift in a warehouse.Durmić (2019) evaluated the criteria for sustainable supplier selection using FUCOM method.Baig et al. (2022) used FUCOM and FQFD for prioritizing the vulnerabilities and identifying those capabilities that can ensure protection against these vulnerabilities.
The specific process to get the weight is as follows: (i) Rank the criterions and get the set from I c x = {I c 1 , I c 2 , . . ., I c y } which is according to the relative importance of the criteria.Thus, the parameters rank is obtained by the values of the weight coefficient: where k denotes the order of the criterions.
(ii) The comparison priorities between the adjacent attributes ψ c x−1 x , x = 1, 2, . . ., y, which denotes the value of the N c y(x−1) relative to N c y(x) , then we get the set of the criterion comparative preference: , ψ c x−1 x .
(iii) Get the final weight coefficients of the attributes p c x = {p c 1 , p c 2 , . . ., p c y }, which should meet the conditions showing as follows: • The comparison priorities ψ c x−1 x calculated in (ii) are supposed to be equal to the ratio of the weight coefficient through Eq. ( 9): • The second condition is about the weight coefficients which should satisfy the following rule by using Eq. ( 10): Thus, the inequality constraints for this model are shown in Eq. ( 11): With the help of the MATLAB software, we get the final result of the weighting vector of the evaluation criterions for each DM.Then the integrated weight px is finally obtained by geometric means.
Step 3. Integrate the decision matrices by different DMS into one matrix vrx ( ĝrx ) = {t rx ( ĝrx )] s×y , using Eq. ( 12 Step 4. Figure out the score function of the integrated decision matrix by using Eq. ( 13) Step 5. Determine the v * x and v− x indexes by the following equation by using Eq. ( 14): Step 6. Calculate the positive and negative distance by using Eqs.( 15)-( 16) Step 7. Compute out the relative closeness to the ideal solution.
The relative closeness of alternative Q r with the probabilistic hesitant fuzzy positive ideal solution Q * is described by using Eq. ( 17): Step 8. Get the final ranking.
The final rank is determined by the value of relative closeness, and the sort rule is: the bigger the relative closeness is, the more appropriate the scheme is.
Step 2. Use FUCOM method to calculate the criterion weight.
• Calculate the criterion weight for DM 1.
(i) The DM 1 gives the ranking of significance of different attributes: (ii) Table 7 shows the priorities of each attribute which is range from 1 to 4, based on the comparison in the former step.According to the data in Table 7, we get the comparative priorities as follows.
Table 4 The standardized decision matrix by the first DM.

Table 7
The priorities of criteria of DM 1.

Criteria
I 1 1 (iii) A finite model for criterion weight coefficient meeting the condition which is introduced in the above: The weight can be calculated by the software Lingo, and the result is p i x = {0.531,0.197, 0.140, 0.133}, and the result of χ 1 is 0.00.
• Calculate the criterion weight for DM 2 (see Table 8).
(i) The DM 1 gives the ranking of significance of different attributes (ii) According to the data in Table 8, we get the comparative priorities as follows: The weight can be calculated by the software LINGO, and the result is p 2 x = {0.160,0.106, 0.222, 0.511}, and the result of χ 2 is 0.00.
• Calculate the criterion weight for DM 3 (see Table 9).
(i) The DM 1 gives the ranking of significance of different attributes  (ii) According to data in Table 9, we get the comparative priorities as follows: The weight can be calculated by the software LINGO, and the result is p 3 x = {0.200,0.519, 0.124, 0.157}, and the result of χ 3 is 0.00.
The final criterion weight is obtained by the integration weight combined with experts' decision weight, and the result is px = {0.287,0.299, 0.158, 0.256}.
Step 3. Integrate the decision matrices by different DMS into one matrix (see Table 10).
Step 4. Figure out the score function of the integrated decision matrix (see Table 11).
Step 5. Determine the v * x and v− x indexes (see Tables 12-13).
Table 12 The positive index v * x .
Step 6. Calculate the positive and negative distance (see Table 14).
Step 7. Compute out the relative closeness to the ideal solution.
Then the result is H r = {0.396,1.000, 0.274, 0.105}, and we a get rank

Comparison and Discussion
In this section, TODIM (PHF-TODIM) method (Zhang et al., 2018), PHFWA operator (Xu and Zhou, 2017), PHFWG operator (Xu and Zhou, 2017) are utilized to compare with the PHF-FUCOM-TOPSIS method to test its feasibility and practicability.In order to compare the results more intuitively, we represent the result as a line chart in Fig. 2 and Table 15 where the original result is processed by the same manner in range 0 to 1.
From the above detailed analysis, it could be seen that these four given models have the same optimal choice Q 2 and the order of these four methods is the same.This verifies that the PHF-FUCOM-TOPSIS is reasonable & effective.These four given models have their given advantages: (1) PHFWA operator emphasizes group decision influences; (2) PHFWG operator emphasizes individual decision influences; (3) The PHF-TODIM method is an interactive multi-criteria decision-making method.The method is based on the value function of prospect theory, establishes the relative dominance function of a certain plan compared with other plans according to the psychological behaviour of decision makers, and selects the best plan according to the size of the dominance, so as to determine the optimal plan.At the moment, the TODIM method is continuously improved and Table 15 The sequence from different methods.

Method
The sequence The best alternative PHF-TODIM (Zhang et al., 2018) (Xu and Zhou, 2017) (Xu and Zhou, 2017) widely used in decision-making in various fields.(4) The "ideal solution" and "negative ideal solution" in the PHF-FUCOM-TOPSIS method are two basic concepts of the TOP-SIS method.The so-called ideal solution is an assumed optimal solution (scheme), and its various attribute values reach the best value among the alternative schemes; while the negative ideal solution is the assumed worst solution (scheme), and each of its attribute values achieve the worst value among the alternatives.The rule for sorting the schemes is to compare the alternatives with the ideal solution and the negative ideal solution.If one of the alternatives is closest to the ideal solution while far from the negative ideal solution, then this scheme is the best one among the alternatives.

Conclusions
In this study, we propose a new PHF-FUCOM-TOPSIS model for performance evaluation of public charging service quality and apply it in the probabilistic hesitant fuzzy environ-ment.A novel extended TOPSIS model integrated with FUCOM method was proposed to evaluate green selection supplier.Finally, we apply this method in a numerical study for performance evaluation of public charging service quality and compare the results with other methods to test its validity.The specific contributions of it are as follows: (1) It integrates classical TOPSIS method and FUCOM method in the probabilistic hesitant fuzzy environment including more information to make the decision-making process more reasonable.(2) It extends the FUCOM method to calculate criterion weight in the probabilistic hesitant fuzzy environment.
In the future, we firmly believe that PHF-FUCOM-TOPSIS method will be applied in a larger number of fields.Meanwhile, we should consider the attributes of the actual situation when solving the performance evaluation of public charging service quality and apply this new model in more fields.

Table 1
Decision matrix c 1 given by the first DM.

Table 2
Decision matrix c 2 given by the second DM.

Table 3
Decision matrix c 3 given by the third DM.

Table 9
The priorities of criteria of DM 3.

Table 11
The score of the integrated decision matrix.