A Novel Hybrid Fuzzy Grey TOPSIS Method: Supplier Evaluation of a Collaborative Manufacturing Enterprise

: Recently, there is of signiﬁcant interest in developing multi-criteria decision making (MCDM) techniques with large applications for real-life problems. Making a reasonable and accurate decision on MCDM problems can help develop enterprises better. The existing MCDM methods, such as the grey comprehensive evaluation (GCE) method and the technique for order preference by similarity to an ideal solution (TOPSIS), have their one-sidedness and shortcomings. They neither consider the di ﬀ erence of shape and the distance of the evaluation sequence of alternatives simultaneously nor deal with the interaction that universally exists among criteria. Furthermore, some enterprises cannot consult the best professional expert, which leads to inappropriate decisions. These reasons motivate us to contribute a novel hybrid MCDM technique called the grey fuzzy TOPSIS (FGT). It applies fuzzy measures and fuzzy integral to express and integrate the interaction among criteria, respectively. Fuzzy numbers are employed to help the experts to make more reasonable and accurate evaluations. The GCE method and the TOPSIS are combined to improve their one-sidedness. A case study of supplier evaluation of a collaborative manufacturing enterprise veriﬁes the e ﬀ ectiveness of the hybrid method. The evaluation result of di ﬀ erent methods shows that the proposed approach overcomes the shortcomings of GCE and TOPSIS. The proposed hybrid decision-making model provides a more accurate and reliable method for evaluating the fuzzy system MCDM problems with interaction criteria.


Introduction
Multi-criteria decision making (MCDM) problems are utterly fatal and can be frequently observed academically in many domains, such as in manufacturing [1,2], finance [3], logistics [4], and the supply chain [5]. The rationality and correctness of the decision are quite crucial for a company. The development of an effective and efficient MCDM method to solve the optimization problems correctly is significant. The GCE method depends on the grey correlation degree, which can reflect the closeness degree between the evaluation sequences of objects [41]. It does not take into account the number and distribution of samples and is easy to implement. However, most studies focus on the construction of the grey correlation coefficients. If the shapes of evaluation sequences of some objects are very similar, the GCE method may not be able to evaluate them properly. For example, assume that the evaluation sequences U1, U2, and U3 are shown in Figure 2. The correlation degree between U2 and U1 is the same as that between U3 and U1. The GCE method depends on the grey correlation degree, which can reflect the closeness degree between the evaluation sequences of objects [41]. It does not take into account the number and distribution of samples and is easy to implement. However, most studies focus on the construction of the grey correlation coefficients. If the shapes of evaluation sequences of some objects are very similar, the GCE method may not be able to evaluate them properly. For example, assume that the evaluation The GCE method depends on the grey correlation degree, which can reflect the closeness degree between the evaluation sequences of objects [41]. It does not take into account the number and distribution of samples and is easy to implement. However, most studies focus on the construction of the grey correlation coefficients. If the shapes of evaluation sequences of some objects are very similar, the GCE method may not be able to evaluate them properly. For example, assume that the evaluation sequences U1, U2, and U3 are shown in Figure 2. The correlation degree between U2 and U1 is the same as that between U3 and U1. Other similar traditional methods, such as TOPSIS and GCE, assume that criteria are independent of each other. In this regard, to a great extent, it leads to limiting them when addressing some practical problems with interactive criteria. By considering these features of the TOPSIS and GCE approaches, this study develops a novel hybrid MCDM method called the fuzzy grey TOPSIS (FGT). By comparing the exiting methods, this research makes the following contributions.
First, the proposed method combines the GCE and TOPSIS, which improves the one-sidedness of the GCE and TOPSIS methods. Decision-makers can make a decision according to their preference through adjusting a parameter that reflects the weights of the shape and the distance factors.
Second, fuzzy measures and fuzzy integral are employed to express and integrate the interaction between the criteria that cannot be properly expressed through the existing methods.
Third, fuzzy numbers are used to help the experts make a more reasonable and accurate evaluation of some qualitative criteria.
The remainder of this paper is organized as follows: Section 2 is a literature review. Section 3 introduces the preliminaries. Section 4 presents the novel hybrid MCDM method that combines the fuzzy integral, the grey correlation, and TOPSIS. Section 5 utilizes an example of supplier evaluation of a collaborative manufacturing enterprise to illustrate the procedure, the feasibility, and effectiveness of the proposed method. Finally, conclusions are reached in Section 6.

Literature Review
Academically, many studies have been conducted on MCDM problems, and substantial  Other similar traditional methods, such as TOPSIS and GCE, assume that criteria are independent of each other. In this regard, to a great extent, it leads to limiting them when addressing some practical problems with interactive criteria. By considering these features of the TOPSIS and GCE approaches, this study develops a novel hybrid MCDM method called the fuzzy grey TOPSIS (FGT). By comparing the exiting methods, this research makes the following contributions.
First, the proposed method combines the GCE and TOPSIS, which improves the one-sidedness of the GCE and TOPSIS methods. Decision-makers can make a decision according to their preference through adjusting a parameter that reflects the weights of the shape and the distance factors.
Second, fuzzy measures and fuzzy integral are employed to express and integrate the interaction between the criteria that cannot be properly expressed through the existing methods.
Third, fuzzy numbers are used to help the experts make a more reasonable and accurate evaluation of some qualitative criteria.
The remainder of this paper is organized as follows: Section 2 is a literature review. Section 3 introduces the preliminaries. Section 4 presents the novel hybrid MCDM method that combines the fuzzy integral, the grey correlation, and TOPSIS. Section 5 utilizes an example of supplier evaluation of a collaborative manufacturing enterprise to illustrate the procedure, the feasibility, and effectiveness of the proposed method. Finally, conclusions are reached in Section 6.

Literature Review
Academically, many studies have been conducted on MCDM problems, and substantial progress has been made in this field. For example, Bouzarour-Amokrane et al. formulated a consensus bipolar method to solve the collaborative group MCDM problems considering the impact of human behavior, e.g., individualism, fear, caution [42]. The analysis results of the application example about real size wind farm implantation problem showed that this proposed method can be a useful decision-making tool to solve this problem [42]. Tchangani et al. proposed a bipolar aggregation method that combines weighted cardinal fuzzy measures, which can effectively overcome difficulties that dissuade the use of Choquet integral in practices to solve the fuzzy nominal classification in the MCDM problem [43]. Zhang et al. formulated a hybrid optimization approach combining the best worst method, grey relational analysis, and visekriterijumsko kompromisno rangiranje (VIKOR) to solve the MCDM problem in rail transit [44]. A multi-cell thin-walled aluminum energy-absorbing structure was applied to verify that this integrated method is valid and practical [44]. The results proved that this method provides an accurate and effective tool for the structural decision-making problem in rail transit [44]. Mousavi-Nasab and Sotoudeh-Anvari presented a new multi-criteria decision-making approach for the sustainable material selection problem [45]. Ten examples were applied to verify this proposed method [45].
TOPSIS and GCE are two of the several important methods in dealing with MCDM problems. Many researchers have a contribution to hybridization of MCDM techniques with some other method Appl. Sci. 2019, 9, 3770 4 of 25 to make the best use of their advantages and bypass their disadvantages, as shown in Tables 1 and 2. In addition, the fuzzy integral and fuzzy set are used to address the interaction among criteria which have been applied to many domains, as exposed in Table 3.

The λ-Fuzzy-Measure and the Choquet Integral
In the FGT method, the weights of all subsets of a criterion set are expressed through the λ-fuzzy-measure.
. , x n } be a finite set. If the fuzzy density of x i is g i , then g λ can be calculated as follows: Choquet integral [79]: Without loss of generality, assume that the evaluation values on n criteria where g λ X j = g λ x 1 , x 2 , . . . , x j represents the fuzzy measure of the corresponding criterion set.

The Correlation Coefficients
Assume that the evaluation matrix is where the i-th row represents the evaluation sequence of the i-th object. {C * (.)} = {C * (1), C * (2), . . . ,C * (n)} is selected as the reference sequence where C * (l) represents the best of the l-th criterion. The correlation coefficient vector of the i-th alternative sequence related to the reference sequence is denoted as ξ i * = ξ i * (1), ξ i * (2), . . . , ξ i * (n) . We then have [80]: where the resolution ratio ρ ∈ [0, 1] is usually assigned as 0.5.
The absolute correlation degree can be found by The bigger the degree value is, the better the object is.

The Entropy Weight Method
With the premise that there are m objects and n criteria for evaluation, which form the evaluation matrix Z = z ij ; i = 1, 2, . . . , m; j = 1, 2, . . . , n , the weights of all criteria can be calculated through the entropy weight method [82], as shown in Equation (7). (7)

Background
Assume that there are m objects to be evaluated. The objects set is A = {a 1 , a 2 , . . . , a m }. There are p experts involved in the evaluation and decision. The experts set is denoted as To evaluate these objects, the experts should construct a reasonable and effective criterion system according to the feature of the considered problem and the purpose and preference of the decision-maker. As it is widely acknowledged, a real-life system usually contains some unknown information, i.e., it is a grey system. Therefore, some qualitative criteria need to be employed to evaluate the system. In addition, we also need a few quantitative criteria for some measurable aspects.
Assume that among n criteria, some qualitative and some quantitative criteria are considered. The criterion set is X = {x 1 , x 2 , . . . , x n }. For the convenience of expression, assume that x 1 , x 2 , . . . , x i are quantitative criteria and x i+1 , . . . , x n are qualitative.

The Procedure of Fuzzy Grey TOPSIS
Step 1: Find the evaluation data. Step1.1: Obtain evaluation data of qualitative criteria.

Obtain the Fuzzy Linguistic Value and Transform it into Fuzzy Numbers
To make the experts' evaluation more accurate, a set of fuzzy linguistic values is set up and noted as Superb, Good, Normal, Bad, Terrible ({S, G, N, B, T} for short). The evaluation data of qualitative criteria are given by experts in the form of fuzzy linguistic values that correspond to fuzzy numbers. Mapping rules are shown in Table 4. Table 4. Mapping rules of linguistic variables and fuzzy numbers.
The linguistic variable evaluation matrixes are transformed as fuzzy number matrixes, as shown below.
where J i represents the fuzzy number matrix of the i-th object j k ij represents the fuzzy number corresponding to the k-th expert's evaluation of the i-th object on the j-th criterion.
The fuzzy number evaluation values of k experts need to be integrated into one fuzzy number according to Equation (6). The fuzzy evaluation matrix J can be given as where j ij represents the fuzzy number of the i-th object on the j-th criterion.

Defuzzify the Fuzzy Numbers
The fuzzy numbers need to be defuzzified into numerical values to take part in the calculation procedure later. There are many defuzzification methods, each of which has its advantages and limits. Three methods are adopted in this paper. In the case that the operation is too simple, and its effectiveness cannot be verified [83], the mean value of the results of the three methods is taken. Fuzzy numbers j ik = (a ik , b ik , c ik ) can be defuzzicated through the three methods stated as follows: (1) Distance measure method [81]: If the optimal fuzzy evaluation value is defined as j ik * = (1, 1, 1) and the worst fuzzy evaluation value is defined as j ik − = (0, 0, 0), then (2) Central value method [84]: (3) Gravity method [84]: The numerical value j ik can be found by taking the mean value of the three results, i.e., The fuzzy evaluation matrix J can be transformed into a numerical evaluation matrix J * .
The evaluation data of a few numerical criteria j k1 , j k2 , . . . , j ki k = 1, 2, . . . , m are obtained through statistics and measurement methods. By seaming the numerical evaluation matrix of qualitative and quantitative criteria, the evaluation matrix J can be found, as shown below, where J kl represents the evaluation value of the k-th object on the l-th criterion, k = 1, 2, . . . , m, l = 1, 2, . . . , n.
To eliminate the influence of dimension and order of magnitude, the evaluation matrix J needs to be standardized according to Equation (14): where j k is the mean and S k is the standard deviation of the evaluation values of the k-th criterion.
In that way, the standardized evaluation matrix C can be found.
Theorem 1. The standardization transformation satisfies the following two axioms: (1) Isotonicity: If j xk < j yk , then c xk < c yk . If j xk > j yk , then c xk > c yk .
Proof. If j xk < j yk , then j xk − j k < j yk − j k . Thus (2) Difference remaining: Step 3: Obtain weights of criteria.
To aggregate the evaluation data, some series of weights that are reasonable and accurate as much as possible is certainly needed. The weights are considered in two aspects. On the one hand, the objective weights are calculated through the entropy weight method such that the influence of data can be considered. On the other hand, the subjective weights are given by the experts to express the decision preference. Finally, the weighted average value of objective and subjective weights is taken.
Step 3.1: Calculate objective weights of criteria.
Step 3.2: Obtain subjective weights of criteria.
where the parameter θ ∈ (0, 1) is determined by the decision maker according to the degree of trust in objective data and subjective judgment of experts. We set it as 0.5 in this paper.
It is generally believed that many interactions exist among criteria. If we integrate the evaluation matrix by weights only, the interactions are ignored, which will lead to an inadequate decision. Thus, the λ-fuzzy-measures are employed here to express the interactions.
Step 5: Calculate the grey fuzzy correlation degree.
To describe the difference in shapes of all objects, their positive and the negative correlation coefficients need to be calculated. First of all, the positive and the negative ideal solution schemes C + = c + 1 , c + 2 , . . . , c + n and C − = c − 1 , c − 2 , . . . , c − n should be selected. With the premise that J + represents the set of some criteria which are the bigger the better, and J − represents the set of some criteria which are the smaller the better, the vector C + and C − can be determined according to Equations (19) and (20).
To take the interactions among criteria into account, the Choquet fuzzy integral is used to integrate the correlation coefficients of every object. The positive grey fuzzy integral vector R + = R + 1 , R + 2 , . . . , R + m and the negative grey fuzzy integral vector R − = R − 1 , R − 2 , . . . , R − m can be calculated as follows.
The positive grey fuzzy Choquet integral of the i-th object is defined as where ω 12...k represents the weight of the set {x 1 , x 2 , . . . , x k } (k ≤ n) after rearrangement.
The negative grey fuzzy Choquet integral of the i-th object is defined as Step 6: Calculate fuzzy distance to ideal solution.
To embody the difference the in distance of all objects, the distance between all objects to the positive and negative ideal solution schemes should be obtained. The Choquet fuzzy integral should also be employed to amend the Euclidean distance such that the interactions can be expressed in the distance.
The positive fuzzy ideal solution distance of the i-th object is defined as The positive ideal solution fuzzy distance of the i-th object is defined as Step 7: Get grey fuzzy TOPSIS.
The positive grey fuzzy Choquet integral and the negative ideal solution distance depend on the close degree to the positive ideal solution of each scheme on the shape and distance, respectively. Their combination can represent the comprehensive close degree to the positive ideal solution of each scheme. Analogously, the combination of negative grey fuzzy Choquet integral and the positive ideal solution distance can characterize the close degree to the negative ideal solution of each scheme. To express the effects of R + , R − , D + , and D − , they need to be normalized before combination.
where N i represent R + (R − ; D + ; D − ), respectively. For the conciseness, the normalized values of them are still denoted as R + , R − , D + , D − , respectively. The positive comprehensive proximity degree vector S + = S + 1 , S + 2 , . . . , S + m and the negative comprehensive proximity degree vector S − = S − 1 , S − 2 , . . . , S − m can be obtained by combining R + , R − , D + , and D − with a proper weight parameter p according to Equation (28).
The m objects can be evaluated reasonably according to the comprehensive evaluation vector CS = [CS 1 , CS 2 , . . . , CS m ]. The object whose comprehensive evaluation value is big is better than the one whose comprehensive evaluation value is small.
The flow chart of the proposed method is given in Figure 3.
The objects can be evaluated reasonably according to the comprehensive evaluation vector = [ , , … , ]. The object whose comprehensive evaluation value is big is better than the one whose comprehensive evaluation value is small. The flow chart of the proposed method is given in Figure 3.:

An Illustrative Example
To verify the feasibility and effectiveness of the proposed method, an illustrative example of supplier evaluation of a collaborative manufacturing enterprise is given in this paper. The result was compared with GCE and TOPSIS. This example has six alternatives = { , , , , , } and 12 experts = { , , … , }.
According to the feature of the supplier evaluation of the collaborative manufacturing enterprise and the purpose and preference of the decision-maker, a reasonable and effective criterion system was constructed by experts, as shown in Figure 4.
Step 5: Get grey fuzzy correlation degree Step 4: Identificate λ-fuzzy-measures Step 3: Get weights of criteria Step 1: Get evaluation matrix Step 1.1：Get the evaluation data of qualitative criteria 1. Get the fuzzy linguistic value and transform it into fuzzy numbers

Defuzzify the fuzzy numbers
Step1.2： Obtain the evaluation data of quantitative criteria Step 2: Standardize evaluation matrix Step 3.1: Get the objective weights of criteria Step 3.2: Get the subjective weights of criteria Step 3.3: Integrate objective and subjective weights Step 5.1: Get correlation coefficients Step 5.2: Get grey fuzzy Choquet integral Step 6: Get fuzzy distance to ideal solution Step 7: Get grey fuzzy TOPSIS

An Illustrative Example
To verify the feasibility and effectiveness of the proposed method, an illustrative example of supplier evaluation of a collaborative manufacturing enterprise is given in this paper. The result was compared with GCE and TOPSIS. This example has six alternatives A = {A 1 , A 2 , A 3 , A 4 , A 5 , A 6 } and 12 experts E = {E 1 , E 2 , . . . , E 12 }.
According to the feature of the supplier evaluation of the collaborative manufacturing enterprise and the purpose and preference of the decision-maker, a reasonable and effective criterion system was constructed by experts, as shown in Figure 4.

Conducting the Proposed Method
The detailed process of the proposed method to solve this example is illustrated below. Each level is addressed one by one. We illustrate the quality level in detail only for economy of space.
Step 1: Get the evaluation matrix.
The fuzzy linguistic value evaluation matrices {̃ ,̃ ,̃ ,̃ ,̃ ,̃ } are shown as follows, in which the mapping rules of linguistic variables and fuzzy number are given in Table 4.

B B B T T B B B B T B B
J N N N G N N N G N N G G  N G G N G S G G G G N N N N G N N N G N B N N  N N N B N B B B N N N B  N N N B N N N

S S G S S S G G S G S S J
Financial position x 61 Inventory control level x 62 Ratio of professional personnel x 63 Staff training level x 64 Employee discipline level x 65 Corporate culture x 66 Delivery period x 71 Ability of delivering on time x 72 Proficiency of technics x 32 Ability of planning and control x 33 Investment in research x 41 Research capacity x 42 Technique ripeness level x 43 Attitude towards customers x 51 Cycle of problem feedback x 52 After sales system level x 53 Supply flexibility x 73 Information sharing level x 74 Risk sharing ability x 75 Cooperation stability x 76 Communication cost x 24 Advanced equipment level x 31 Quality certification level x 11 Pass-rate of production x 12 Production experience x 13 Quality control ability x 14 Technical personnel level x 15 Detailed level of quote Cooperation level

Conducting the Proposed Method
The detailed process of the proposed method to solve this example is illustrated below. Each level is addressed one by one. We illustrate the quality level in detail only for economy of space.
Step 1: Get the evaluation matrix.
The fuzzy linguistic value evaluation matrices j 1 , j 2 , j 3 , j 4 , j 5 , j 6 are shown as follows, in which the mapping rules of linguistic variables and fuzzy number are given in Table 4. N N N G N N N G N N N N N G N N N G N B N N  N N N B N B B B N N N B  N N N B N N N N N G G N N N G N G N N  N  G N G N N N N N G N N  N N G N N N G N N N N where each matrix represents the evaluation data of an alternative, and each row in the matrix represents the evaluation data of one qualitative criterion given by 12 experts expect for the second row of each matrix that is the evaluation data of the quantitative criterion and has just one value. In fact, we can give the evaluation of qualitative criteria and quantitative criteria, respectively, but for the convenience of writing, we give them in one matrix.
We transformed the fuzzy linguistic values into fuzzy numbers according to (6). Then these fuzzy numbers were aggregated according to Equation (6) and defuzzified according to Equations (8)- (13). Then the evaluation matrix J was found as follows: where each row represents the evaluation data of an alternative, and each column represents the evaluation data of a criterion.
The 3-dimensional histogram of J is portrayed in Figure 5. Step 2: Standardize evaluation matrix.
The evaluation matrix can be standardized as according to Equation (14), as shown as below. Step 2: Standardize evaluation matrix.
The evaluation matrix J can be standardized as C according to Equation (14), as shown as below.
The 3-dimensional histograms of C is depicted in Figure 6. Step 2: Standardize evaluation matrix.
The evaluation matrix can be standardized as according to Equation (14), as shown as below.
The 3-dimensional histograms of is depicted in Figure 6. Step 3: Obtain weights of criteria.
Step 3.1: Calculated the objective weights of criteria. The objective weight vector which was calculated by Equation (15)  Step 3.3: Integrate objective and subjective weights. Step 3: Obtain weights of criteria.
Step 3.1: Calculated the objective weights of criteria. The objective weight vector which was calculated by Equation (15)  Step 3.2: Obtain the subjective weights of criteria. The subjective weight vector which was given by experts was ω 2 = [0.1990, 0.6040, 0.2080, 0.2990, 0.3100].
Step 3.3: Integrate objective and subjective weights.
The comprehensive weights vector can be obtained by combining the objective weights vector and the subjective weights vector according to Equation (16). It was ω 0 = [0.2103, 0.4069, 0.2119, 0.2404, 0.2405].
Step 4: Identify λ-fuzzy-measures. The parameter λ is identified according to Equation (17), which was −0.5262. The weights of all subsets of the power set of the criterion set can be calculated according to Equation (18), as shown in Table 5. Step 5: Calculate the grey fuzzy correlation degree.
Step  (21) and (22), the positive correlation matrix ξ + and the negative correlation matrix ξ − were found as follows.
Step 5.2: Calculated grey fuzzy Choquet integral. The positive and negative grey fuzzy integral vectors can be obtained by applying Equations (23) and (24)  Step 7: Calculate the comprehensive proximity degree. Normalize R + , R − , D + , D − according to Equation (27). Set the parameter p as 0.5, which means the importance of shape and distance is the same. The positive and negative comprehensive proximity degree vectors can be calculated by applying Equation (28)  Intuitively, Figure 6 shows that enterprise 1 performed very well on all criteria; while enterprise 2 performed very badly on all criteria. Enterprise 4 performed averagely on criteria 1, 2, 4, and 5 but performed excellently on criterion 3. Enterprises 5, 6, and 3 all performed averagely and enterprise 5 performed slightly better than enterprise 6, and enterprise 6 performed slightly better than enterprise 3 on most criteria. The sequence coincided with our intuitionistic judgment.
The bar chart of R + , R − , D + , and D − is depicted in Figure 7. It shows that the values of R + , R − , D + , and D − were not coincident. If we concentrate on different aspects, we can obtain different evaluation results.
Intuitively, Figure 6 shows that enterprise 1 performed very well on all criteria; while enterprise 2 performed very badly on all criteria. Enterprise 4 performed averagely on criteria 1, 2, 4, and 5 but performed excellently on criterion 3. Enterprises 5, 6, and 3 all performed averagely and enterprise 5 performed slightly better than enterprise 6, and enterprise 6 performed slightly better than enterprise 3 on most criteria. The sequence coincided with our intuitionistic judgment.
The bar chart of , , , and is depicted in Figure 7. It shows that the values of , , , and were not coincident. If we concentrate on different aspects, we can obtain different evaluation results.

Comparisons and Sensitivity Analysis
To compare the effects of the proposed method with the existing methods, the result of GCE and TOPSIS is given below.
By standardizing , the weights of GCE can be obtained. The absolute correlation degree (ACD) can be calculated according to the weights and the positive correlation matrix. The ideal solution distance and can be found by the standardized evaluation matrix. The relative closeness degree (RCD) can be obtained through and . To carry out a sensitivity analysis of parameter , a series of experiments was conducted and the value of was between 0 and 1, and the interval was 0.2. The result is shown in Table 6 and Figure 8.

Comparisons and Sensitivity Analysis
To compare the effects of the proposed method with the existing methods, the result of GCE and TOPSIS is given below.
By standardizing w 0 , the weights of GCE can be obtained. The absolute correlation degree (ACD) can be calculated according to the weights and the positive correlation matrix. The ideal solution distance D + TOPSIS and D − TOPSIS can be found by the standardized evaluation matrix. The relative closeness degree (RCD) can be obtained through D + TOPSIS and D − TOPSIS . To carry out a sensitivity analysis of parameter p, a series of experiments was conducted and the value of p was between 0 and 1, and the interval was 0.2. The result is shown in Table 6 and Figure 8. Table 6 and Figure 8 show that the sequence sorted according to GCE was always 1 5 4 6 3 2 while the sequence sorted due to TOPSIS was always 1 4 6 5 3 2. The result means that enterprise 4 is bad in terms of shape but good in terms of distance since it performed well on criterion 3 but badly on criteria 1 and 2. Enterprise 5 is bad in terms of distance, but good in terms of shape since it performed slightly better than average on most criteria, which is similar to enterprise 1.

Result Analysis
The sequence sorted according to FGT is slightly related to parameter p. When p = 0, the sequence was 1 4 5 6 2 3. when p > 0, the sequence was 1 4 5 6 3 2. Enterprise 2 has slight advantages on distance, since it performed barely satisfactorily on criteria 4 and 5. However, its performance was instable. When the factor of shape is taken into consideration, enterprise 3 will be better than enterprise 2. Figure 8 clearly shows that the comprehensive evaluation value of enterprise 1 (CS 1 ) obviously decreased, CS 4 slightly decreased, CS 2 , CS 3 , and CS 5 obviously increased, and CS 6 slightly increased with the increase of parameter p, i.e., the increase consideration of shape.
The result shows that the proposed method can reflect both the difference of shape and distance of objects, which overcomes the one-sidedness of the GCE method and TOPSIS method.   Table 6 and Figure 8 show that the sequence sorted according to GCE was always 1 ≻ 5 ≻ 4 ≻ 6 ≻ 3 ≻ 2 while the sequence sorted due to TOPSIS was always 1 ≻ 4 ≻ 6 ≻ 5 ≻ 3 ≻ 2. The result means that enterprise 4 is bad in terms of shape but good in terms of distance since it performed well on criterion 3 but badly on criteria 1 and 2. Enterprise 5 is bad in terms of distance, but good in terms of shape since it performed slightly better than average on most criteria, which is similar to enterprise 1.
Enterprise 2 has slight advantages on distance, since it performed barely satisfactorily on criteria 4 and 5. However, its performance was instable. When the factor of shape is taken into consideration, enterprise 3 will be better than enterprise 2. Figure 8 clearly shows that the comprehensive evaluation value of enterprise 1 ( ) obviously decreased, slightly decreased, , , and obviously increased, and slightly increased with the increase of parameter , i.e., the increase consideration of shape.
The result shows that the proposed method can reflect both the difference of shape and distance of objects, which overcomes the one-sidedness of the GCE method and TOPSIS method.

Conclusions
Generally speaking, the MCDM problems are essential in modern society for businesses and private individuals. The existing MCDM methods, such as GCE and TOPSIS, have some

Conclusions
Generally speaking, the MCDM problems are essential in modern society for businesses and private individuals. The existing MCDM methods, such as GCE and TOPSIS, have some disadvantages which cannot take the interaction of criteria into consideration. In the view of this, a novel MCDM method that combines the GCE, TOPSIS, and fuzzy integral, called the fuzzy grey TOPSIS, was proposed to improve the deficiency mentioned above. The effectiveness of the proposed method was verified by the illustrative example of supplier evaluation of a collaborative manufacturing enterprise as well as some sensitivities.
There are some recommendations for future studies. One idea is to focus on other types of fuzzy numbers and correlation coefficients and even redefine them. We will compare them by conducting more experiments and refine the proposed method [85][86][87]. As such, hybridizing the proposed MCDM method with recent advances in heuristics and metaheuristics is another good continuation of this work [88][89][90][91][92][93].

Conflicts of Interest:
The authors declare no conflict of interest.

j ik
The evaluation value of the i-th object on the k-th qualitative criterion in the form of fuzzy number, i ∈ A, k ∈ x k x k is a qualitative criterion . For the convenience, j a i x k is presented in a simplified form as j ik .
The numerical results of defuzzification of j ik through methods 1, 2, and 3.
Numerical evaluation matrix of qualitative criteria.
Evaluation matrix.