Intuitionistic Fuzzy Three-Way Decision Model Based on the Three-Way Granular Computing Method

: Three-way decisions, as a general model for uncertain information processing and decisions, mainly utilize the threshold generated by the decision cost matrix to determine the decision category of the object. However, the determination of the threshold is usually accompanied by varying degrees of subjectivity. In addition, the potential symmetrical relationship between the advantages and disadvantages of the decision cost is also a problem worthy of attention. In this study, we propose a novel intuitionistic fuzzy three-way decision (IFTWD) model based on a three-way granular computing method. First, we present the calculation methods for the possibility of membership state and non-membership state, as well as prove the related properties. Furthermore, we investigate the object information granules, i.e., the ﬁne-grained, medium-grained, and coarse-grained objects, by combining the state probability distribution and probability distribution. Then, for decision and evaluation issues, we deﬁne the superiority-compatibility relation and inferiority-compatibility relation for IFTWD model construction. In addition, we use the superiority degree and inferiority degree instead of the original thresholds and design a new method for evaluating decision cost. Finally, we focus on the algorithm research of the proposed model and present an empirical study of agricultural ecological investment in Hubei Province to demonstrate the e ﬀ ectiveness of our


Introduction
Pawlak proposed the rough set theory [1], which rapidly attracted the attention of scholars and was successfully applied to machine learning [2,3], formal concept analysis [4][5][6], knowledge engineering [7,8], granular computing [9,10], etc. The rough set mainly utilizes the equivalence relation to express the classification ability of domain knowledge. In detail, the domain of discourse is divided into three disjoint subsets according to the upper and lower approximations, i.e., the positive domain, the negative domain, and the boundary domain [11]. However, such strict equivalence relations will easily cause excessive noise and redundant information in the boundary domain. From the perspective of improving the fault tolerance of the boundary, scholars have introduced various improved rough set models [12][13][14][15][16].
In particular, based on the decision-theoretic rough sets (DTRS) model [13], Yao proposed the three-way decision theory [17,18], which provided a reasonable semantic interpretation for the positive, boundary, and negative domains of the rough set. In simple terms, corresponding to these three domains, different strategies are adopted, i.e., acceptance, non-commitment, and rejection [19]. As a more general decision model and information-processing mode, three-way decision theory has received extensive attention from scholars since it was proposed. To date, the research on three-way decisions can be categorized into the following two aspects, namely, modelling and application. In terms of model construction, the traditional method is to study three-way decisions by constructing different (1) The concepts of membership state possibility and non-membership state possibility are defined, and their related properties are proved. Then, according to the relationship between the state probability distribution and the probability distribution, we give an approximate division rule for the object information granularity. (2) By combining object information granules and dominance relation, we propose superiority-compatibility relation and inferiority-compatibility relation. Furthermore, we deduced the IFTWD model on the basis of the superiority-compatibility class and inferiority-compatibility class. (3) To objectively evaluate the decision-making results and analyze the potential symmetrical relationship between decision costs, we define the concepts of advantage cost and disadvantage cost, and create a secondary decision strategy for boundary domain objects.
As we stated above, the main purpose of this investigation is to construct an intuitionistic fuzzy three-way decision model through the three-way granular computing method. The rest of the structure of this paper is as follows. In Section 2, we simply review some basic knowledge of three-way decisions and intuitionistic fuzzy sets. In Section 3, we propose a method for calculating the object information granularity and then define the superiority-compatibility relation to construct the IFTWD model. Finally, a new strategy for analysing decision costs is given. In Section 4, the validity of the model proposed in this paper is verified by an example. In Section 5, the paper ends with conclusions.

Preliminaries
This section briefly reviews some basic concepts [1,13,17,18,41,54,55]. Throughout this paper, we assume that the domain of discourse U and attribute set A are both non-empty finite sets.

Pawlak Rough Set
Pawlak believes that knowledge can be divided by equivalence relations to form an easy-to-express formal language that helps people understand and calculate. Definition 1. [1] Let U be the universe of discourse; R is an equivalent relationship defined on U, and apr = (U, R) is an approximate space. The division of the universe U under the equivalence relation R can be described as U/R = [x] R x ∈ U , where[x] R is an equivalence class containing x. For any X ⊆ U, its upper and lower approximate sets can be expressed as: It can be found that both the upper and lower approximations are crisp sets, and the elements in the upper approximation overlap with the equivalent classes, while the elements in the lower approximation are completely contained in the equivalent classes. Therefore, the discourse domain can be divided into three disjointed parts by upper and lower approximations, namely, the positive domain POS(X), the negative domain NEG(X), and the boundary domain BND(X) are as follows: POS(X) = apr R (X), NEG(X) = U − apr R (X), BND(X) = apr R (X) − apr R (X). (2)

Three-Way Decision
Due to the strict equivalence relations, the positive and negative domains obtained by the Pawlak method are correct and meet the requirements, i.e., there is no wrong acceptance or rejection. Furthermore, in practice, the BND region often has too much redundant information. To overcome its shortcomings, Yao proposed three-way decision theories combined with DTRS and Bayesian minimum risk theory, which gave rough sets a new semantic interpretation. Definition 2. [13,17,18] Let Ω = {X, ¬X} be a set representing two complementary states; i.e., element x belongs to object X or does not belong to object X. Γ = {a P , a B , a N } represents three different decisions, where a P is the acceptance decision, a B is the delay decision, and a N is the rejection decision. The decision losses caused by actions taken in different states are shown in Table 1. Table 1. Cost function matrix.

Action
Cost Function In Table 1, λ PP , λ BP , and λ NP represent the cost losses of adopting a P , a B , and a N decisions when in state X. Similarly, λ PN , λ BN , and λ NN denote the decision loss under state ¬X. Therefore, the expected loss of object x under different actions is R(a i [x])(i = P, B, N) , which can be expressed as follows: where Pr(X [x] R ) and Pr(¬X [x] R ) represent the probabilities that the equivalence class belongs to X and ¬X.
According to the Bayesian minimum risk decision process, we can obtain the following decision rules: By a reasonable hypothesis, λ PP ≤ λ BP < λ NP and λ NN ≤ λ BN < λ PN ; if the BND domain exists, (4) can be simplified as: Here, α and β represent: Through the above analysis, we can find that the three-way decision rules are closely related to the (α, β) thresholds. In the uncertain environment, eliminating the constraints of DTRS and constructing decision evaluation methods from other perspectives is a topic worthy of further research.

Intuitionistic Fuzzy Set and Dominance Relation
Definition 3. [54] Let U be the universe and and ν S (x) are the membership and non-membership degrees of object x, respectively, and η S ( is called the hesitancy degree or intuitionistic fuzzy index. It is easy to find that when η S (x) = 0 , the intuitionistic fuzzy set degenerates into a fuzzy set. For convenience, S = < x, µ S (x), ν S (x) > x ∈ U is abbreviated as S = (µ S (x), ν S (x)), and all intuitionistic fuzzy sets on U are denoted as IF(U).

Definition 4.
[54] Let U be the universe; x ∈ U are two intuitionistic fuzzy sets on U. Then, the basic operation is as follows: Then, (Ps, Ng) represents an intuitionistic fuzzy possibility measure on U.

Definition 6.
[55] Let quaternion Ψ = (U, A, V, f ) represent an information system, where U = {x 1 , x 2 , · · · , x n } is the universe, A = {B 1 , B 2 , · · · , B m } is the attribute set, V is the value domain of attribute B j , and f : U × A → V is an information mapping function. The superiority and inferiority relations are expressed as follows: Superiority relation: For the convenience of description, if V B j ∈ IF(U), we use Ψ = (U, A, V, f ) to represent an intuitionistic fuzzy information system, abbreviated as Ψ.

IFTWD Derived from the Three-Way Granular Computing Method
The TAO model is a generalized framework for three-way granular computing. To enrich its connotation, we first construct object information granules by fusing state probability functions of uncertain information features to solve the trisecting problem. Second, for the problems of acting and outcome, the concept of a superiority-compatible relation is given, and the decision results are evaluated by the superiority degree and inferiority degree. Finally, the secondary decision is made on the objects in the boundary domain, and the IFTWD model decision rules and algorithms are derived.

State Possibility and Information Granularity
Reference [41] considers the composition of intuitionistic fuzzy probability measured from a macro perspective, but from a micro perspective, if the overlap of the membership degree of elements is not large or the conflict of non-membership degree is not strong, then the useful information of used elements cannot be integrated and utilized. If attribute B ∈ A has multiple states under object x i ∈ U, the possibility in a certain state is simply referred to as the state possibility. Definition 7. Let the triplet Θ = U, A, Sq be a state possibility space, where U is the universe, A is the attribute set, and Sp(Sq + (x i ), Sq − (x i )) is the state possibility function established between U and A. Specifically, Sq + and Sq − indicate the membership state possibility and non-membership state possibility, respectively. For any x i ∈ U, B, C ∈ A, satisfying the following conditions: Then, Sq is called the state possibility function. On the basis of (M1)-(M3), if Sq satisfies: Then, Sq is called the weakened non-membership state possibility function. Similarly, if Sq satisfies: Then, Sq is called the enhanced membership state possibility function.
(M1) indicates that the state probability of any object is between 0 and 1. In addition, since probability only describes the difficulty of an event, (M2) gives the extreme value condition of state possibility, and (M3) is the monotonic condition of state possibility. Moreover, it is easy to determine that the membership range of the state possibility is between [µ(x i ), µ(x i ) + η(x i )] and that the non-membership state possibility is between [ν(x i ), ν(x i ) + η(x i )], so (M4) and (M5) can be regarded as the weakening and strengthening conditions of the state possibility.
The following will take max = µ B (x i ) as an example to prove. Since i.e., they have a negative correlation.
; it is easy to determine that κ B and ν B (x i ) change in the same direction, and the logarithmic function with e as the base monotonically increases in the domain, which leads to a positive correlation between ν B (x i ) and ln(ε B + 1).
In summary, is the state possibility function. Remark. Equation (7) considers the effects of the degree of hesitation converted into the membership state possibility and the non-membership state possibility. On the one hand, the degree of membership and non-membership represent deterministic information to a certain extent. On the other hand, due to the existence of unknown information, the degree of hesitation has the arbitrary possibility of turning into membership and non-membership. Therefore, the following points need to be considered in actual calculations: (1) Satisfy the weak and strong conditions of (M4) and (M5). Specifically, there is always a possibility that (1,0) satisfies the enhanced membership state and (0,1) satisfies the weakened non-membership state. (2) The following two state values need special consideration to make the calculation results meaningful. Case 1: when the state value is (1,0), we only need to consider the membership state possibility. Case 2: when the state value is (0,1), we only need to consider the non-membership state possibility. (3) When the state value is (0, 0), the conversion of hesitation degree is completely uncertain. Relative to the determined membership and non-membership degrees, the intuitionistic fuzzy index is caused by the existence of a variety of unknown information. For the sake of analysis, it is assumed that these unknown factors are independent. According to the central limit theorem, these independent unknown factors obey the Gaussian distribution (hypothetical standard Gaussian distribution), and then the transformation degree of the intuitionistic fuzzy index can be expressed as follows: Proof. It is straightforward to prove Lemma 1 by use of Definition 7 and Theorem 1.
Proof. When µ B (x i ) = 0, the possibility of a membership state is completely determined by the hesitation degree, and the amount of membership information contained in the hesitation degree is uncertain. According to Equation (7), Example 1. Reference [56] calculated the risk levels of different audit objects in the intuitionistic fuzzy security audit evaluation system, where the universe is U = {x 1 , x 2 , x 3 , x 4 } and the attribute set is A = B, C, D, E, F, G , as shown in Table 2. According to Theorem 1, we can obtain the membership state possibility and non-membership state possibility of each object under attribute set A, as shown in Figure 1. Example 1 demonstrates that the higher the possibility of the membership state of the object x i under the attribute B ∈ A is, the lower the possibility of the corresponding non-membership state. In addition, the membership state probability and the non-membership state possibility are given based on the maximum transformation condition of hesitation degree, which can only be quantitatively analysed from the degree of difficulty of a certain state. Therefore, the similarity or consistency between the probability distribution and the state probability distribution needs further analysis. as shown in Table 2.
According to Theorem 1, we can obtain the membership state possibility and non-membership state possibility of each object under attribute set A, as shown in Figure 1. Example 1 demonstrates that the higher the possibility of the membership state of the object i x under the attribute B A ∈  is, the lower the possibility of the corresponding non-membership state. In addition, the membership state probability and the non-membership state possibility are given based on the maximum transformation condition of hesitation degree, which can only be quantitatively analysed from the degree of difficulty of a certain state. Therefore, the similarity or consistency between the probability distribution and the state probability distribution needs further analysis.  Table 3. Table 3. Relationship between state probability distribution and probability distribution. Definition 8. Let quaternion Λ = U, A, Sp, P be a state probability-probability space, where U is the universe of discourse, A is the attribute set, Sp is the state probability function built between U and A, and P is the probability conversion function. ∀x i ∈ U, ∀ B ∈ A , if the state value (µ B (x i ), ν B (x i )) of the object x i with respect to the attribute B is regarded as a two-dimensional random variable about the membership function and the non-membership function, then the relationship between the probability distribution and the state probability distribution is shown in Table 3. Table 3. Relationship between state probability distribution and probability distribution.

Comparison and Classification of Values Probability Distributions
If the value is (0,0), let P k (0, 0) = 0. Definition 9. Let ψ = (U, A, V, f ) be an intuitionistic fuzzy information system. For any x i ∈ U, B j ∈ A, there is a state value (µ B (x ij ), ν B (x ij )) corresponding to f, and the matrix form is M = (µ B (x ij ), ν B (x ij )) n×m = (x ij ) n×m . The expression ability of object x i in the complete set A, i.e., the information granularity G, can be expressed as follows: If there is consistently more than one Sp η (x ij ), the mean value is taken as the final result. In Equation (9), θ η ensures that the information uncertainty when the state probability is converted into a probability value remains unchanged, and G(x ij ) provides the degree of fitting between the state probability and probability distributions after adding uncertain information, i.e., the smaller the information granularity G of an object is, the stronger its expressive power. Since probability is generally understood as the difficulty of an event, probability is understood as the frequency of an event. In the intuitionistic fuzzy system ψ, a high probability of a certain state means that it has a high probability of occurrence. Conversely, a low probability of a certain state implies that it has a low probability of occurrence. Therefore, the assumptions without loss of generality are P big ( Combining Table 2 and Equation (9), we can obtain the approximate division of object information granularity, as shown in Table 4. Table 4. Approximate division of information granularity based on state probability.
"O" and "×" indicate that the conditions are satisfied and not satisfied, respectively. The approximate division rules for the granularity of object information obtained from Table 4 are as follows: (1). If C 1 and C 2 , then x i ∈ G F . (2). If C 1 and ¬C 2 , or ¬C 1 and C 2 , then x i ∈ G M . (3). If ¬C 1 and ¬C 2 , then x i ∈ G C .
It is worth noting that the conditions C 1 and C 2 are just a trend judgement. In practice, we also need to perform quantitative analysis and calculation. In addition, to increase the noise tolerance of the C 1 and C 2 conditions, the Euclidean distance of the state probability and probability value is used as the elasticity judgement interval, which is expressed as follows: Sq + max (x i ) represents the maximum membership state possibility of object x i under attribute set A. Due to the existence of uncertain information, the difference between the Euclidean distance of Sq + max (x i ), Sq η (x i ) and P k (x i ) is used as the lower bound of the judgment condition P big (x i ) → Sq big (x i ) . Similarly, the sum of Euclidean distance between the Sq − min (x i ), Sq η (x i ), and P k (x i ) is used as the upper bound as the judgment condition Sq small (x i ) → P small (x i ) . Thus, the approximate division conditions C 1 and C 2 in Table 4 can be adjusted as: Example 2. According to the relationship between the state probability distribution and the probability distribution as well as the division conditions, the information granularity of each object in Table 3 can be calculated as shown in Table 5. Table 5. Information granularity of each domain object in Table 2. The granular classification algorithm of the universe object is given below. First, the membership state possibility and non-membership state possibility of the universe object x i under attribute set A are calculated. Then, according to the relationship between the state probability distribution and the probability distribution, the information granularity value of object x i is obtained. Finally, according to the C * 1 and C * 2 judgement conditions, the information granularity of the universe object is approximately classified.
Algorithm 1 Granularity classification of objects in the universe of discourse Object x i is Medium-grained: x i ∈ G M // satisfies either C * 1 or C *

Superiority-Compatibility Relation and Decision Evaluation
According to the view of three-way granular computing [51][52][53], decision and evaluation are another focus after trisecting work. Therefore, in this section, the concept of the superiority-compatibility relation and the calculation method of the superiority degree are given. Accordingly, a new decision cost evaluation strategy is proposed, and the secondary decision is made for the BND domain objects. Definition 10. Let Ψ = (U, A, V, f ) be an intuitionistic fuzzy information system; then, R ≤ δ is called the superiority-compatibility relation if and only if the following condition is satisfied: Analogously, inferiority-compatibility R ≥ ε is as follows: Equations (12) and (13) use the superiority-compatibility relation and inferiority-compatibility relation to replace the traditional equivalence relation, and ||[x i ] R ≤ and [x i ] R ≥ represent the number of elements in the maximum superiority-compatibility class and the maximum inferiority-compatibility class of object x i , respectively. Moreover, R ≤ is defined as an object with a membership degree that is not greater than δ= min{ and R ≥ is defined as an object with a membership degree that is not less than ε= max{ Thus, the superiority-compatibility class and inferiority-compatibility class are as follows: Theorem 2.
represents the largest superiority-compatibility classes.
represents the largest inferiority-compatibility classes.
Proof. Since the methods are similar, only (1) is demonstrated below.
From Equations (12)- (14), it is easy to find that the superiority-compatibility relation introduces the dominance relation in the information granularity of the object and integrates the compatibility relation of the attribute set. Therefore, ∀X ⊆ U, the superiority lower approximation R(X) and inferiority upper approximation R(X) can be obtained as follows: It can be seen that the lower approximation R(X) is a further conditional restriction on fine-grained objects and that the upper approximation R(X) is a stronger judgement condition on coarse-grained objects. Theorem 3. ∀X, Y ⊆ U, we have the following conclusions: Proof. It is straightforward to prove Theorem 1 by use of Definition 10.
Proof. Suppose x i ∈ POS, according to formula (14), we can obtain [ Based on further analysis, for any X ⊆ U, the lower approximate set expression ability is strong and has full advantages, while the upper approximate set expression ability is weak and completely at a disadvantage. Therefore, the POS(X) domain, BND(X) domain and NEG(X) domain can be obtained as follows: Corresponding to Equation (16), the division rules of the POS domain, NEG domain, and BND domain can be expressed as follows: Example 3. Combined with the granularity information obtained in Example 2, the calculation of the superiority-compatibility class, the inferiority-compatibility class, the upper approximation set, the lower approximation set, and the classification result of each object is shown in Table 6.
From Example 3, on the one hand, from the perspective of object information granularity, the partition results obtained from R ≤ δ and R ≥ ε are related to the original granularity of the object to a certain extent. For example, x 1 ∈ G F , and the final decision result is x 1 ∈ POS. On the other hand, due to the characteristics of R ≤ δ and R ≥ ε , there is generally only one object in the POS domain and NEG domain, which is quite different from the classical three-way decision division rules. In this way, although the correctness of the POS domain and NEG domain is guaranteed, there will be too many delayed decision objects in the BND domain, so it is particularly important to make circular decisions on the BND domain objects. Moreover, in the three-way granular computing method, apart from the trisecting and decision, evaluation also plays an important role. To solve the above problems and evaluate the decision results reasonably and objectively, the following will analyse the decision result through the advantage cost and disadvantage cost of the object.

Definition 11.
Let G ς (x i )(ς ∈ {F, M, C}) be the granularity information, and [x i ] R ≤ δ is the superiority-compatibility class; then, the advantage cost h(x i ) of object x i is as follows: Correspondingly, the disadvantage cost of object x i is as follows: where |·| represents the cardinality of the set. It can be known that the coarser the granularity of the object x i is, the fewer disadvantaged elements contained in the corresponding superiority-compatibility class, which results in a smaller advantage cost. For the convenience of calculation, the number of , and the superiority degree can be obtained as follows: Similarly, the inferiority degree of each object in the universe is as follows: Equations (20) and (21) can be used to calculate the difference of superiority degree and inferiority degree to judge the advantages and disadvantages of object x i . There are two main cases: Proof. It can be directly proven by Theorem 4 and Definition 11. Definition 12. Let U = {x 1 , x 2 , · · · , x n } be the discourse domain, and the superiority degree matrix and the inferiority degree matrix on U are H(x) and L(x), respectively. Then, the secondary decision rules for objects in the BND domain are as follows: The secondary decision in the BND domain is a cyclic process. Each P3 and N3 rule can obtain the object with the highest superiority degree and the highest inferiority degree, respectively. Therefore, the BND domain set becomes U − R − R − P3 − N3 after one cycle, and the above process is repeated until the BND domain is empty. In addition, when there is only one object left in the BND domain, Based on Algorithm 1, the superiority-compatibility relation and inferiority-compatibility relation will be used to calculate the superior-compatibility class, inferior-compatibility class, and their upper and lower approximate set of the object. Then, the POS domain, NEG domain and BND domain are obtained. Finally, the decision results are evaluated according to the superiority degree and inferiority degree of the object, and a secondary decision is made for the BND domain.
then Calculate the dominance cost and disadvantages cost of each object do ), i j 12: end for 13: for Get the advantages and disadvantages degree of each object do 14: 17: x i ∈ POS // having the maximum superiority degree 18: if 19:  (7) x 7 , Southwest of Hubei. In addition, the evaluation indicators can be divided into 3 categories: (1) Q 1 , ecological benefit; (2) Q 2 , economic benefit; and (3) Q 3 , social benefit. The costs of decision assessments for different agroecological regions are also different. Hence, the investment decisions of the agroecological regions of Hubei Province are consistent with three-way decisions [58]. The committee evaluates the performance of agroecological regions in the years 2004-2006 [57], and the weight ω = (ω(t 1 ), ω(t 2 ), ω(t 3 )) T is (0.2, 0.3, 0.5), where t 1 , t 2 , and t 3 represent "2004 year", "2005 year" and "2006 year", respectively. In the following, we will make decisions on the above problems through the IFTWD model proposed in this paper.
Step 1 Use Equations (7)- (11) to divide the grain size of different agroecological regions.
Step 2 According to Equations (16) and (17), calculate the superiority-compatibility class and inferiority-compatibility class of each agroecological region from t 1 to t 3 , and make three-way decision divisions, as shown in Tables 7-9.  [ Table 9. Division results of the agroecological region in t 3 . [ Step 3 Based on Equations (18)- (21) and combined with the weights taken over the years, calculate the superiority degree and inferiority degree of each agroecological region.
Step 4 Make a secondary decision on the BND domain by Equation (22), as shown in Tables 10-12. Table 10. Agroecological region secondary decision in t 1 . Table 11. Agroecological region secondary decision in t 2  Tables 7-9, it can be seen that the number of superiority-compatibility class elements is the largest in the POS domain, while the inferiority-compatibility class elements are the least in the NEG domain. According to Theorem 5, from t 1 to t 3 , the element with the most intersections in the superiority-compatibility class is x 3 , and the element with the most intersections in the inferiority-compatibility class is x 1 . The results are consistent with x 1 ∈ POS and x 3 ∈ NEG. Moreover, the relationship between the superiority degree of the object and inferiority degree obtained according to R ≤ δ and R ≥ ε is shown in Figure 2, where the upper semi-axis of x represents the superiority degree and the lower semi-axis represents the inferiority degree. Taking t 1 as an example, element x 1 in the POS domain is completely in an advantageous position, and element x 3 in the NEG domain is completely in an inferior position. The superiority degree and inferiority degree of other elements are not balanced in t 1 − t 3 , which requires further analysis. From Tables 10-12 show the secondary decision evaluation of the BND domain. 12,1068 in the comprehensive evaluation results of the committee, the total number o groecological regions within 1 3 t t is compared, as shown in Figure 3.  To obtain the comprehensive evaluation results of the committee, the total number of decisions of different agroecological regions within t 1 − t 3 is compared, as shown in Figure 3.

Initialization First Cycle Second Cycle Final Result
In Figure 3, object x 1 has the highest acceptance rate, x 7 has a rejection rate of 70%, and the remaining delayed decision parts are also converted to reject decisions after the secondary decision.
Similarly, object x 4 is transformed into an accept decision, and objects x 3 and x 5 are transformed into a reject decision after the second decision. Combined with Tables 10-12, the final decision result is  shown in Table 13. re 2. Advantages and disadvantages of each agroecological region between gure 3. Acceptance rate, rejection rate and delay rate in the secondary decisio re 3, object 1 x has the highest acceptance rate, 7 x has a rejection rate of 70 elayed decision parts are also converted to reject decisions after the seconda ject 4 x is transformed into an accept decision, and objects 3 x and 5 x are t decision after the second decision. Combined with Tables 10-12, the final dec  Table 13.  The final decision result of our model is that x 1 , x 2 , and x 4 regions can be invested, and x 3 , x 5 , x 6 , and x 7 regions cannot be directly invested. If we can only invest in one agroecological region in the POS domain, we can choose x 1 according to the relationship between superiority degree and inferiority degree, which is consistent with the results given in reference [57]. Although reference [58] did not further analyse the BND domain, based on the results of the POS domain and NEG domain, we believe that x 1 and x 4 are suitable for investment immediately and that x 7 is not suitable for investment. The only difference is in the decision division of the object x 2 . There are two possible reasons for our analysis: one is that reference [58] uses hypothetical conditional probability, which may cause our final results to be different, and the other is the effect of the state probability, as shown in Figure 4.
State possibility and advantage cost we believe that 1 and 4 are suitable for investment immediately and that 7 is not suitable for investment. The only difference is in the decision division of the object 2 x . There are two possible reasons for our analysis: one is that reference [58] uses hypothetical conditional probability, which may cause our final results to be different, and the other is the effect of the state probability, as shown in Figure 4. State possibility and advantage cost   Figure 4 shows the effect of state probability on advantage cost and information granularity by taking 1 Q as an example. We can find that the state possibility affects the formation of information granularity and indirectly affects the calculation of advantage cost. In addition, Figure 2 also shows that the change from 1 t to 3 t of object 2 x is unstable. All of these factors may cause inconsistency in the final decision of 2 x .  Figure 4 shows the effect of state probability on advantage cost and information granularity by taking Q 1 as an example. We can find that the state possibility affects the formation of information granularity and indirectly affects the calculation of advantage cost. In addition, Figure 2 also shows that the change from t 1 to t 3 of object x 2 is unstable. All of these factors may cause inconsistency in the final decision of x 2 .

Conclusions
In this paper, we propose an IFTWD model based on the granular computing method. First, we construct the object information granularity by the relationship between the state probability distribution and the probability distribution. Second, to solve the problem of action and evaluation, a superiority-compatibility relation and inferiority-compatibility relation are proposed, and the decisions are divided. Finally, the calculation methods of superiority degree and inferiority degree are given and used for the secondary decision. This study extends the TAO model by using information granulation, decision stratification, and evaluation subdivision to analyse the uncertainty problem and construct a corresponding IFTWD model, which provides new ideas for the study of three-way granular computing and the three-way decision model under uncertain conditions. The next step will focus on the flexible properties and operator construction of three-way granular computing on intuitionistic fuzzy multiscale problems and the corresponding three-way decision model. Author Contributions: X.X. initiated the research and wrote the paper, J.S. supervised the research work and provided helpful suggestions, and W.P. participated in some of these search work. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.