Interval-valued Pythagorean fuzzy multi-criteria decision-making method based on the set pair analysis theory and Choquet integral

This paper proposes a novel fuzzy multi-criteria decision-making method based on an improved score function of connection numbers and Choquet integral under interval-valued Pythagorean fuzzy environment. To do so, we first introduce a method to convert interval-valued Pythagorean fuzzy numbers into connection numbers based on the set pair analysis theory. Then an improved score function of connection numbers is proposed to make the ranking order of connection numbers more in line with reality in multi-criteria decision-making process. In addition, some properties of the proposed score function of connection numbers and some examples have been given to illustrate the advantages of conversion method proposed in the paper. Then, considering interactions among different criteria, we propose a fuzzy multi-criteria decision-making approach based on set pair analysis and Choquet integral under interval-valued Pythagorean fuzzy environment. Finally, a case of online learning satisfaction survey and a brief comparative analysis with other existing approaches are studied to show that the proposed method is simple,convenient and easy to implement. Comparing with previous studies, the method in this paper, from a new perspective, effectively deals with multi-criteria decision-making problems that the alternatives cannot be reasonably ranked in the decision-making process under interval-valued Pythagorean fuzzy environment.


Introduction
An important issue in multi-criteria decision-making (MCDM) is to obtain a reasonable ranking order of all alternatives. Due to the complexity of reality, fuzzy and uncertain information is naturally involved in MCDM process. For this reason, the theory of fuzzy set [1], intuitionistic fuzzy set (IFS) [2], interval-valued intuitionistic fuzzy set (IVIFS) [3] and their applications have been put forward one after another with the development of research [4][5][6][7][8][9][10] Yager [11] proposed such an example in real life: a decision maker may express his satisfaction with an alternative on a criterion is 0.6, but his dissatisfaction is 0.5. Because 0.6 + 0.5 > 1, the above special case cannot be modeled by the theory of IFS or IVIFS , which requires the sum of membership degree and the non-membership degree less than or equal to one [12]. Therefore, a concept of Pythagorean fuzzy set (PFS) is introduced by Yager, of which the square sum of membership degree and non-membership degree is less than or equal to one [11]. As extensions of PFS, Smarandache [13] introduced the refined Pythagorean fuzzy sets, Ünver [14] defined Spherical Fuzzy Sets and Zhang [15] proposed a concept of interval-valued Pythagorean fuzzy set (IVPFS). As powerful tools to deal with vagueness and uncertainty involved in MCDM problems, theories and applications of these sets have recently been extensively studied in the literature. For instance, Ejegwa [16] solved career placement problems under the Pythagorean fuzzy environment. Saeed et al. [17] showed the properties, set-theoretic operations and axiomatic results for the refined Pythagorean Fuzzy Sets. For more details, please refer to [18][19][20][21][22][23][24][25][26].
It is noteworthy that the ranking technique of fuzzy numbers is crucial in the fuzzy MCDM process [15]. That is, if the ranking technique is not appropriate then no matter what fuzzy MCDM method you use, the results are unreliable or even completely inconsistent with reality [27]. On the other hand, the interval-valued Pythagorean fuzzy numbers (IVPFNs) involves more uncertainties than other fuzzy numbers (e.g., intuitionistic fuzzy numbers (IFNs), intervalvalued intuitionistic fuzzy numbers (IVIFNs), Pythagorean fuzzy numbers (PFNs),etc) which are usually able to adapt to higher degrees of uncertainty [28]. In order to make the solutions to MCDM problems more reliable, it is necessary to develop a ranking method which not only ranks IVPFNs intuitively but also loses useful information as little as possible [27]. Ever since IVPFNs' appearance, many studies have focused on the ranking problems under interval-valued Pythagorean fuzzy environment. For instance, Zhang [15] proposed a ranking method based on the closeness index of PFNs and IVPFNs and presented a Pythagorean fuzzy hierarchical qualitative flexible multiple criteria approach (QUAL-IFLEX) to solve the fuzzy MCDM problems. Although the method is relatively simple, it relies too much on the definition of distance for PFNs and IVPFNs. That is, different distances will get different ranking results which will bring some inconveniences in fuzzy MCDM process. Moreover, it is noteworthy that the score function and the accuracy function are important tools for ranking PFNs and IVPFNs [15]. Zhang [15] introduced the score function and accuracy function of IVPFNs, which generalized the definition for PFNs in [29]. However, these definitions will lead to a certain loss of information, because they failed to consider the influence of hesitation of IVPFNs in fuzzy MCDM process under the interval-valued Pythagorean fuzzy environment. Therefore, a novel accuracy function of IVPFNs [30], an improved accuracy function [12] and an improved score function of IVPFNs [28] were proposed by Garg considering the effect of hesitation interval index of IVPFNs. These above ranking methods have been widely used in the field of interval-valued Pythagorean fuzzy MCDM.
However, by browsing the literature, we find that there is one type of IVPFNs, the elements in which cannot be reasonably ranked in MCDM process using existing methods. This type of IVPFNs satisfies the following two properties: first, for each IVPFN, the lower limit of its membership degree is equal to the lower limit of its non-membership degree, and the upper limit of its membership degree is equal to the upper limit of its non-membership degree; second, the square sum of lower limit and upper limit of the membership degree of the one set is equal to the square sum of lower limit and upper limit of the membership degree of the other one. According to the existing methods, it can be concluded that these sets are equivalent even if they are completely different ones (see Example 4 for details). In addition, we find the value of improved accuracy function for some IVPFNs may exceed one even if it is not the largest IVPFN (see Example 5 for details). Obviously, all these results are not in line with reality. In view of the above analysis, it is necessary to propose a new ranking approach for IVPFNs from a new perspective to obtain a reasonable order between them.
The set pair analysis (SPA) theory is a new framework combining the certainty and uncertainty into a unified way [31]. The connection number (CN) is a principal mathematical tool of SPA [32]. It uses the degree of 'identity', 'discrepancy', and 'contrary' to indicate the certainty, hesitancy, and uncertainty of a system, respectively [33]. Since the SPA theory was proposed, researchers have done a lot of in-depth studies on its theory and applications under IFS and IVIFS environment [34][35][36][37]. For example, based on the SPA theory, Garg and Kumar proposed some similarity measures of IFSs [38] and some series of distance measures for IFSs [39]. In [40], they introduced a TOPSIS IVIFS MADM method in decision-making process using the SPA theory. And to rank different IVIFNs, Garg and Kumar proposed a new possibility measure of IVIFS based on the CNs of SPA [41]. Kumar and Chen [42] proposed a multi-attribute decision making method based on SPA under the interval-valued intuitionistic fuzzy environment and introduced a score function of connection numbers, and so on. Since IVPFS is the generalization of IVIFS, it can be inferred that the SPA theory can also be a useful tool to deal with uncertainty in MCDM process under interval-valued Pythagorean fuzzy environment. Unfortunately, we have not yet found any application of the SPA theory under the interval-valued Pythagorean fuzzy environment, let alone the research on ranking techniques and multi-attribute decision making methods under the interval-valued Pythagorean fuzzy environment. Moreover, the proposed score function of CNs [42] has some shortcomings which are unable to get the reasonable sorting of alternatives in some MCDM processes under the interval-valued Pythagorean fuzzy environment (for details, see Example 2).
Motivated by above analysis, this paper first develops a novel ranking method for IVPFNs based on the SPA theory. That is, in order to get a reasonable order of IVPFNs in MCDM process, a technique to convert IVPFNs into CNs based on SPA is introduced at first by taking the hesitation interval index and Pythagorean property of IVPFNs into consideration properly. Then we propose an improved score function of CNs which can make the ranking order of CNs more in line with reality. The properties of the proposed score function of CNs and some examples are also given to illustrate the advantage of our proposed ranking method. Next, considering interactions among different criteria in the decision-making process, we propose a fuzzy MCDM approach under interval-valued Pythagorean fuzzy environment based on SPA and Choquet integral which is used to aggregate the evaluation information of criteria for each alternative.
Therefore, the innovations of this paper are summarized as follows : (i) A novel technique of converting IVPFNs into CNs based on SPA is proposed to rank IVPFNs in MCDM process from a new perspective for the first time in the literature. In addition, the idea of transformation fully takes into account the influence of hesitation interval index and Pythagorean property on information uncertainty under the interval-valued Pythagorean fuzzy environment. (ii) An improved score function of CNs is presented which can make the ranking order of CNs more in line with reality, and overcome the shortcomings of existing formulas; (iii) The aggregation of alternative evaluation information using SPA and Choquet integral considers uncertainty and interactions among criteria in MCDM under interval-valued Pythagorean fuzzy environment simultaneously.
This paper is organized as follows. The basic concepts of PFS, IVPFS, the SPA theory, the score function and accuracy function of CNs as well as fuzzy measure and Choquet integral are reviewed in the next section. The technique to convert IVPFNs into CNs and an improved score function of CNs are proposed and some examples are given to verify the advantage of the proposed converting method from IVPFNS into CNs in the following section. A novel MCDM approach based on SPA and Choquet integral under intervalvalued Pythagorean fuzzy environment is introduced in the next section. An example of online learning satisfaction survey and a brief discussion and a comparative analysis with other existing methods are studied to illustrate the simplicity and viability of the proposed fuzzy MCDM approach in the following section. Conclusion of this paper is given in the last part.

Preliminaries
In this section, we first review the basic concepts of PFS, IVPFS, the SPA theory, fuzzy measure and Choquet integral.

Interval-valued Pythagorean fuzzy set
Definition 1 [11]. Let X be a universe of discourse, then a Pythagorean fuzzy set (PFS) P in X can be denoted as The t(x) is the membership degree and m(x) is the non-membership degree of x, respectively, and π p (x) = 1 − (t(x)) 2 − (m(x)) 2 is said to be the hesitancy degree of x to P. For convenience, t p , m p is called a Pythagorean fuzzy number (PFN) by Zhang and Xu [29] with 0 ≤ t p , m p ≤ 1 and 0 ≤ t 2 p + m 2 p ≤ 1. Definition 2 [15]. Let X be a universe of discourse. An interval-valued Pythagorean fuzzy set (IVPFS) P in X is defined as For convenience, an IVPFS can be expressed by an interval-valued Pythagorean fuzzy number be an IVPFN, a score function T of P can be shown as and an accuracy function of IVPFN P can be shown as Suppose there are two IVPFNs Definition 4 [28]. An improved score function R of an It is clear that are two IVPFNs, then an order between P 1 and P 2 can be denoted as follows: Definition 5 [12].
be an IVPFN. An improved accuracy function C of P is denoted by be two IVPFNs, then the following order holds:

Set pair analysis and connection numbers
Definition 6 [31]. Let L(S 1 , S 2 ) is a set pair consisting of two sets S 1 and S 2 for a given problem G . Suppose that there are total N features, I identical features , C contrary features , and D = N − I − C discrepancy features in the problem G, respectively. The connection number (CN) μ of the set pair L(S 1 , S 2 ) for the problem G is defined as where a = I N , b = D N and c = C N show 'identity', 'discrepancy', and 'contrary' degree between S 1 and S 2 , respectively.

Fuzzy measure and Choquet integral
In the MCDM process, the criteria usually have interactions with each other [43]. To illustrate these interactions, a fuzzy measure (or a non-additive measure) of the criteria is proposed by Sugeno [44].
Definition 9 [44]. Let X is a universe of discourse. A fuzzy measure μ on X is a set function μ : £(X ) → [0, 1], satisfying the following conditions: where £(X ) is the power of universe X .

Ranking method for IVPFNs based on SPA
In this section, we deal with the issue of ranking IVPFNs from a new perspective. First, based on the SPA theory, we propose a method to convert IVPFNs into CNs by taking the hesitation interval index and Pythagorean property into consideration simultaneously. Then an improved score function of CNs is given in the following part to make the ranking order of IVPFNs more in line with reality in fuzzy MCDM process.

Conversion from IVPFNs into CNs
Through the literature, we have three findings as follows: (i) The membership degree and non-membership degree of IVPFN are very close to the identity degree and contrary degree of CN respectively; (ii) The hesitation interval index of IVPFN is an important influence factor for MCDM problems with intervalvalued Pythagorean fuzzy information. (iii) Since IVPFNs are the generalization of IVIFNs, the Pythagorean property should be fully considered to avoid the lack of uncertain information under intervalvalued Pythagorean fuzzy environment as much as possible.
In view of the above three findings, we introduce the following definition of conversion from IVPFNs into CNs.

Definition 12 For an IVPFN
where and Then the C N μ p = a p + b p i + c p j of P given by Eqs. (12)(13)(14) in Definition 12 is a rational CN.
Proof For an IVPFN P = [a, b], [c, d] , to prove that the CN μ p = a p + b p i + c p j given by Eqs. (12)(13)(14) is a rational CN, we need to prove that 0 ≤ a p ≤ 1, 0 ≤ b p ≤ 1 and 0 ≤ c p ≤ 1. Since a p ≥ 0, c p ≥ 0, then a p + c p ≥ 0. In view of b p = 1−a p −c p , we just need to prove a p +c p ≤ 1.
In the following, a simple example is given to show how to calculate the CN of an IVPFN.

Improved score function of CNs
Generally speaking, the methods for ranking CNs can be divided into two categories. One is to compare the identity degree and the discrepancy degree between CNs as mentioned before in Definition 7. However, by this method, only a partial order of IVPFNs can be obtained rather than a total order of IVPFNs. The other is to compare the values of score function of CNs to determine the ranking order of CNs as mentioned in Definition 8. However, as shown below, the score function in Definition 8 is unable to rank CNs correctly in some cases.
Therefore, to get a reasonable ordering of CNs, an improved score function of CNs is presented as follows.

Definition 13
Let μ = a + bi + cj be a CN denoted by Eq. (5). An improved score function S F of μ is defined as According to Definition 13, we can obtain the following properties of the proposed score function S F (μ): Obviously, from above, we can obtain that the value of score function S F (μ) of μ is between -1 and 1 for any CN μ = a + bi + cj. In addition, the larger the value of S F (μ), the more forward position of CNs will take in the ranking order. That is, Definition 14 For two CNs μ 1 = a 1 + b 1 i + c 1 j and μ 2 = a 2 +b 2 i +c 2 j, the following order between μ 1 and μ 2 holds: Example 3 If we apply the score function as given in Eq. (15) to the Example 2, we get S F (μ 1 ) = 0.2, S F (μ 2 ) = 0.125 and S F (μ 3 ) = 0.096. Since S F(μ 1 ) > S F(μ 2 ) > S F(μ 3 ), we obtain that the ranking order is μ 1 μ 2 μ 3 . It is clear that the result is consistent with the result by Definition 7.
In the following, a few examples are presented to illustrate the advantage of proposed approach in this section.   ) in Definition 4, we obtain R(P 3 ) = R(P 4 ) = R(P 5 ) = R(P 6 ) = 0. Therefore, we can get P 3 ∼ P 4 ∼ P 5 ∼ P 6 . Obviously, P 3 = P 4 = P 5 = P 6 .
In the following, we use the proposed approach in this paper to reorder the four IVPFNs. That is, Step 1. Convert the four IVPFNs into CNs by Eqs. (12)(13)(14), then get Step 2. Calculate the score value of CNs by Eq. (15), then get Step 3. Compare the above score values, then get Therefore, the order of CNs is μ 5 ≺ μ 6 ≺ μ 3 ≺ μ 4 . It is noticed that the order is also consistent with the result by Definition 7. Hence, based on above steps, the ranking order of the four IVPFNs is P 5 ≺ P 6 ≺ P 3 ≺ P 4 . and C(P 8 ) = 1.1631 > 1. Therefore, we obtain P 8 P 7 . However, it is obvious that P 8 ≺ P 7 , because P 7 = P max is the largest IVPFN as mentioned before. In the following, we use the proposed method in this section to reorder the two IVPFNs. Similar to the steps in Example 4,

Example 5 Let
Step 1. Convert P 7 and P 8 into CNs of them by Eqs. [12][13][14] and get Step 2. According to Eq. (15), calculate the score value of CNs, then get Step 3. Compare the above score values and get Therefore, we obtain the order P 7 P 8 . It is clear that it is a reasonable order more in line with reality.

In fact, let P = {P
, by those previous methods [12,28,30], we get P i ∼ P j even if P i = P j for any i = j. However, by the approach proposed in this section, a reasonable order between them can be obtained as shown in Example 4 above. In fact, the proposed approach in this section is suitable for all IVPFNs ranking problems in interval-valued Pythagorean fuzzy MCDM process. For example, by Eq. (15), the score values of CNs μ p 1 and μ p 2 in Example 1 are S F (μ p 1 ) = −0.1247 and S F (μ p 2 ) = −0.1712. Recall the score values of CNs obtained above in then we obtain the ranking order P 2 ≺ P 1 ≺ P 5 ≺ P 6 ≺ P 3 ≺ P 4 ≺ P 8 ≺ P 7 . That is, this proposed ranking approach is robust and it is more suitable for the ranking problem of IVPFNs than those previous methods. At the same time, as a ranking technique from a new perspective, the proposed approach is simple, effective and easy to implement.

Interval-valued Pythagorean fuzzy multicriteria decision making method
In this section, considering interactions among different criteria, we introduce a fuzzy MCDM method under interval-valued Pythagorean fuzzy environment by taking the advantages of proposed ranking technique above and Choquet integral.
Assume that, there are m alternatives, denoted by A = {A 1 , A 2 , ..., A m }, and n criteria, indicated as C = {C 1 , C 2 , · · · , C n }. Let P i j = [a i j , b i j ], [c i j , d i j ] ∈ P be the estimation for alternative i on criteria j. Here [a i j , b i j ] shows the degree of satisfaction, while [c i j , d i j ] is the degree of dissatisfaction of alternative i with respect to criteria j, where i = 1, 2, ..., m, j = 1, 2, ..., n. In the following, to find the most desirable alternative in interval-valued Pythagorean fuzzy multi-criteria decision process, the specific steps of the proposed fuzzy MCDM method are given as follows: Step 1. Construct and normalize the interval-valued Pythagorean fuzzy decision matrix D( P i j ) m×n : where of which, B represents the benefit type criteria subset, and T is the cost type criteria subset. Step 2. Convert the above matrix into a CNs matrix: where μ i j = a i j + b i j i + c i j j is the CN of P i j . In addition, a i j , b i j , c i j are calculated by Eqs. (12)(13)(14) in Definition 12, respectively, for i = 1, 2, · · · , m, j = 1, 2, · · · , n. Step 3. Calculate the value of score function S F (μ i j ) of μ i j (i = 1, 2, · · · , m; j = 1, 2, · · · , n) using Eq. (15) to obtain the score function matrix as follows: Step 4. Aggregate the above score values on each criterion to each alternative by Choquet integral. The specific implementation process is as follows: 1) Calculate the λ-fuzzy measure of each criterion set.
In general, to obtain the λ-fuzzy measure of each criterion set, we should first get the weight, or the importance, of each criterion [43]. At present, there are three typical methods to determine the weight of criteria, i.e., the subjective method, the objective method and hybrid method [47]. No matter which method is used, the weight is obtained by a numerical value from 0 to 1 [46]. In this paper,we denote the weight of criteria j by g λ (C j ), j = 1, · · · , n. According to the Eq. (9), we get Then we can obtain the parameter λ by resolving Eq. (17) if we have get all g λ (C j )( j = 1, · · · , n) in some way. Based on the parameter λ, the fuzzy measure of each criteria set can be derived by Eq. (7). 2) Calculate the Choquet integral of the score values of CNs for each alternative. According to the above fuzzy measures, we can aggregate the score values of CNs for each alternative with respect to each criterion as follows: for any i ∈ {1, 2, · · · , m}. Here the subscript ((1), (2), · · · , (n)) is a permutation on X such that S F (μ i(1) ) ≤ S F (μ i (2) ) ≤ · · · ≤ S F (μ i(n) ) and Step 5. According to the aggregated values above, obtain a reasonable ranking order of all alternatives to find the most desirable one.
The complete flow chart of the proposed method is given in Fig. 1.

Discussion and comparative study
To illustrate the convenience and feasibility of our proposed approach, an example about the online learning satisfaction survey is given at first. Then the relevant discussion and comparative analysis between the proposed approach and other existing approaches [12,15,28,30,43,49] are provided to validate the performance of the proposed approach in the paper.

Illustrative example
During the COVID-19 pandemic, a large number of Chinese universities implement large-scale online learning, which is a new model that can break the constraints of time and space for students. In the post-epidemic era, how to improve the effect of online learning model has become a very important topic. The online learning satisfaction survey is an important way that can help educators understand the real needs of students. The data in this section come from an online learning satisfaction survey conducted during the COVID-19 pandemic on a questionnaire survey platform [48]. By survey, it is found that the involved students can be divided into four categories according to the nature of their schools: research-oriented university ( A 1 ), application-oriented university (A 2 ), higher vocational college (A 3 ), and other types (A 4 ). In this paper, we study the ranking order of online learning satisfaction for the four types of schools based on SPA and Choquet integral. We select four criteria for research, that is, the service of learning resource (C 1 ), the communication between teachers and students (C 2 ), the design of the course (C 3 ) and the communication between students (C 4 ). By analysis, it is found that the four criteria are not independent(or say there are interactions among them). For example, the more extensive the service of learning resource, and the more interesting the teacher's design of course, the more active the communication between teachers and students, and vice versa. In addition, the communication between teachers and students can also improve the communication between students, and vice versa. Therefore, the decision-making steps are described as below: (Step 1:) Through analysis, it is found that all four criteria are benefit-type. Each type of school A i (i = 1, 2, 3, 4) is given four evaluation values by their students which represent the degree of satisfaction on each criterion C j ( j = 1, 2, 3, 4). The values are provided in the form of IVPFNs. Thus, the normalized Pythagorean fuzzy decision matrix D( P i j ) 4×4 can be obtained by Eq. (16): where the interval [0.5103, 0.7001] means the certain degree that alternative A 1 (research-oriented universities) satisfies C 1 (the service of learning resource), and [0.3565, 0.4103] is the uncertain degree that the alternative A 1 dissatisfies C 1 . The other values in D( P i j ) 4×4 have similar meanings. (Step 2:) Convert the above matrix into the CNs matrix U (μ i j ) 4×4 by Eqs. (12)(13)(14), that is Step 3:) Convert the above matrix into the score function matrix S F (μ i j ) 4×4 by Eq. (15): be calculated by the method in [47]. This is an objective method by training a lot of data from the questionnaire survey platform [48]. The fuzzy measure of each criterion is obtained as follows [47]: Then, by Eq. (17), we get the parameter λ = −0.64. And by Eq. (7), the fuzzy measure of each criteria subset can be obtained in Table 1. (Step 5:) Calculate the overall satisfaction degree of all types schools based on Choquet integral by Eq. (18) as follows: That is, the students from A 2 (application-oriented university) have the highest degree of satisfaction towards Table 1 The fuzzy measures of criteria online learning. This result is consistent with the result of questionnaire survey.

Discussion and comparison analysis
To illustrate the advantages of proposed approach in this paper, we conduct a brief discussion and a comparative analysis with some of the existing approaches [12,15,28,30,43,49] under interval-valued Pythagorean fuzzy environment. These approaches corresponding to them are performed on the considered data and the rankings of alternatives which are summarized in Table 2. As shown in Table 2, it is easily observed that: (i) The approach proposed in this paper, from a new perspective, that is, SPA, deals with fuzzy information in the interval-valued Pythagorean fuzzy multi-criteria decision-making problems. To our knowledge, Zhang [15] proposed a closeness index for IVPFNs and presented the closeness index-based ranking method (QUALIFLEX MCDM approach) for ranking IVPFNs.
In addition, Garg [12,28,30] proposed a novel accuracy function or an improved score function for IVPFNs for solving MCDM problems under interval-valued Pythagorean fuzzy environment. Khan et al. [43,49] introduced an extension of TOPSIS method and a fuzzy GRA method under interval-valued Pythagorean fuzzy environment. Essentially, they are still ranking IVPFNs based on score function and accuracy function of IVPFNs [15] before integration process by the Choquet integral. However, the main idea of the approach proposed in this paper is to convert IVPFNs into CNs using SPA theory before integrating the values of alternatives by Choquet integral. Comparing with the existing approaches, this idea of transformation is first given under the general interval-valued Pythagorean fuzzy environment in the literature. According to the decision-making processes and results in Table 2, the proposed approach in this paper not only can effectively solve the MCDM problems under interval-valued Pythagorean fuzzy environment, but also it is simple, convenient and easy to implement. Meanwhile, the result of our proposed approach is consistent with most existing approaches. (ii) The approach proposed in this paper can deal with the interval-valued Pythagorean fuzzy multi-criteria decision-making problems not only under the general interval-valued Pythagorean fuzzy environment, but also in some special cases which cannot be solved by existing methods [15,28,43,49]. For example, suppose it is still the online learning satisfaction survey question in the above section and the fuzzy measure of all four criteria is also given by Eq. 19. Suppose that the evaluation values of an alternative A 0 on four criteria are P 3 , P 4 , P 5 , P 6 in Example 4, respectively. Since the values of score function for the four IVPFNs are T (P 3 ) = T (P 4 ) = T (P 5 ) = T (P 6 ) = 0 and the accuracy functions of them are V (P 3 ) = V (P 4 ) = V (P 5 ) = V (P 6 ) = 0.58, and the closeness index P(P 3 ) = P(P 4 ) = P(P 5 ) = P(P 6 ) = 0.5 (for details,see [15]), we can get P 3 ∼ P 4 ∼ P 5 ∼ P 6 . In this case, the Choquet integral cannot be performed, because IVPFNs cannot be ranked reasonably. Therefore, the approaches [15,28,43,49] are invalid. However, according to the proposed approach in this paper, we get S F (μ p 3 ) = 0.0882, S F (μ p 4 ) = 0.0892, S F (μ p 5 ) = 0.0782, S F (μ p 6 ) = 0.0824, then P 5 ≺ P 6 ≺ P 3 ≺ P 4 as men- tioned in Example 4. By simple calculation, we get the Choquet integral of the alternative A 0 is C g λ (A 0 ) = 0.0782 * (g λ (C 1 , C 2 , C 3 , C 4 ) − g λ (C 1 , C 2 , C 4 )) + 0.0824 * (g λ (C 1 , C 2 , C 4 ) − g λ (C 1 , C 2 )) + 0.0882 * (g λ (C 1 , C 2 ) − g λ (C 2 )) + 0.0892 * g λ (C 2 ) = 0.0782 * (1 − 0.7) + 0.0824 * (0.7 − 0.63) + 0.0882 * (0.63 − 0.25) + 0.0892 * 0.25 = 0.0850. Therefore, the proposed approach in this paper is robust and can obtain the overall satisfaction evaluation from a new perspective. (iii) The approach proposed in this paper takes the interactions among criteria into account using Choquet integral in MCDM under interval-valued Pythagorean fuzzy environment. While the existing approaches [12,28,30] can only deal with the MCDM problem in which the relationships among criteria are assumed to be independent in advance. At the same time, from the result of comparison which is shown in Table 2, we find that the result of our proposed approach is consistent with the method made by Zhang [15] and Khan [43] which takes into account the interaction between the criteria in MCDM under interval-valued Pythagorean fuzzy environment.

Conclusions
IVPFS is a useful tool for dealing with uncertain information in MCDM process, while the SPA theory has the advantage of combining the certainty and uncertainty into a unified way. In view of these advantages, we propose a method to convert IVPFNs into CNs and introduce an improved score function of CNs for the first time. Furthermore, taking into account the uncertainty in decision-making and interactions between criteria simultaneously, we propose a fuzzy MCDM approach based on SPA and Choquet integral to get the overall degree of satisfaction for each alternative. The key contribution of this study is (1) a novel technique of converting IVPFNs into CNs is proposed considering the uncertainty and Pythagorean property of IVPFNs; (2) an improved score function of CNs is developed to overcome the shortcomings of previous formulas; (3) an effective ranking method for IVPFNs is introduced from a new perspective for the first time; (4) an interval-valued pythagorean fuzzy multi-critria decision-making method based on SPA and Choquet intergral is proposed by considering the interactions among criteria.
From examples shown above, we can see that the proposed approach successfully overcomes the drawbacks presented in [12,15,28,30,42,43,49]. The proposed approach provides us with a very useful way to deal with MCDM problems in the IVPFSs context. At the same time, examples show that the proposed approach is simple and easy to implement in interval-valued Pythagorean fuzzy MCDM process. However, the tendency of decision-makers towards fuzzy indexes (the membership, non-membership or hesitation) of IVPFSs has a significant impact on the outcome of the decision-making in many fields. Here the tendency refers to the decision-makers' attitudes, habits, knowledge or experience, etc. In this paper, we failed to take into account the tendency of decision makers towards different fuzzy indexes of IVPFSs in the process of converting IVPFNs into CNs.In the future, on one hand, we will consider the tendency of decision makers in the process of converting IVPFNs into CNs, and on the other hand, we will also apply the proposed approach to more fields, such as pattern recognition, the security of industrial control, the big data of education, etc.