Abstract
FUZZY INFORMATION AND ENGINEERING 2019, VOL. 11, NO. 1, 105–120 https://doi.org/10.1080/16168658.2021.1886815 ORIGINAL ARTICLE A New Method for Optimal Solution of Intuitionistic Fuzzy Transportation Problems via Generalized Trapezoidal Intuitionistic Fuzzy Numbers a,b b c Darunee Hunwisai , Poom Kumam and Wiyada Kumam Department of Applied Mathematics, Faculty of Science and Technology, Valaya Alongkorn Rajabhat University under the Royal Patronage, Pathumthani, Thailand; KMUTT-Fixed Point Theory and Applications Research Group (KMUTT-FPTA), Theoretical and Computational Science Center (TaCS), Science Laboratory Building, Faculty of Science, King Mongkut’s University of Technology Thonburi (KMUTT), Bangkok, Thailand; Department of Mathematics and Computer Science, Faculty of Science and Technology, Rajamangala University of Technology Thanyaburi, Pathumthani, Thailand ABSTRACT ARTICLE HISTORY Received 23 June 2016 In this paper, we introduce the new method for solving the intu- Revised 18 January 2018 itionistic fuzzy transportation problem (IFTP), by using north-west Accepted 22 January 2019 corner method and modified distribution method to find the optimal solution for IFTP. KEYWORDS Intuitionistic fuzzy set; matrix game; linear programming; intuitionistic fuzzy transportation problem 1. Introduction In 1956, Zadeh [1] firstly defined the concept of fuzzy set theory. The concept of an intuition- istic fuzzy set was proposed by Atanassov in 1986 [2]. This concept referred to the reflection of the relation among ‘1 minus the degree of membership’, ‘the degree of non-membership’ and ‘the degree of hesitation’. The intuitionistic fuzzy set was rasterised by the degree of membership and the degree of non-membership. The intuitionistic fuzzy set had more abundant and flexible than the fuzzy set with uncertain information. Many researchers have also used fuzzy and intuitionistic fuzzy set for solving real world optimisation problems such as transportation problem. The transportation problem is a special kind of optimisation problem. Transportation problem is interested in finding the least total transportation cost of goods in order to satisfy demand at destinations using available supplies at the sources. In usual, transporta- tion problems are solved with the hypothesis that values of supplies and demands and the transportation costs are specified in a precise way. In the real world, in many cases, the decision-maker has no crisp information about the coefficients belonging to the transporta- tion problem. In this situation, the corresponding elements defining the problem can be formulated by mean of fuzzy set, and the fuzzy transportation problem appears in a natural way. In 1941, Hitchcock [3] originally developed the basic transportation problem. Dantzig [4] applied linear programming to solving the transportation problem. Several authors CONTACT Wiyada Kumam wiyada.kum@rmutt.ac.th © 2021 The Author(s). Published by Taylor & Francis Group on behalf of the Fuzzy Information and Engineering Branch of the Operations Research Society, Guangdong Province Operations Research Society of China. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/ licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 106 D. HUNWISAI ET AL. have carried out an examination about fuzzy transportation problem [5–9]. Moreover, sev- eral authors have used intuitionistic fuzzy set theory for solving transportation problems. Hussain and Kumar [10] investigate the transportation problem with the aid of triangular intuitionistic fuzzy numbers (TIFN). Pramila and Uthra [11] presented optimal solution of an IFTP. Antony et al. [12] studied method for solving the transportation problem by using TIFN. Singh and Yadav [13] discussed new approach for solving IFTP of type-2 where the supply, demand are fixed crisp numbers and the cost is TIFN. In this paper, we using a linear ranking function for generalised trapezoidal intuitionis- tic fuzzy numbers (GTrIFNs) to find the IBFS and optimal solution of GTrIFNs based on the allocation of demands and availabilities are real numbers and cots are GTrIFNs. This paper is organised as follows. Section 2 gives the concept of mathematics preliminaries. Section 3 presents ranking of GTrIFN. Section 4 describes a mathematics formulation for IFTP. Section 5 details some numerical example. In the final section, the paper is concluded in Section 6. 2. Mathematical Preliminaries In this section, we give some basic definitions and concepts of cut sets of trapezoidal intuitionistic fuzzy number (TrIFN). 2.1. Some Definitions of TrIFNs Definition 2.1: [1]: Let X be an arbitrary nonempty set of the universe. A fuzzy set A in X is a function with domain X and values in [0, 1]. If A is a fuzzy set and x ∈ X, then the function value μ (x) is called the membership function of x in A. A fuzzy set can be written as order pair, given by {x, μ (x)|x ∈ X} where 0 ≤ μ (x) ≤ 1. A A Definition 2.2: [2]: Let X be an arbitrary nonempty set of the universe. If there are two mapping on the set X: μ (x) : X → [0, 1] and ν (x) : X → [0, 1] with the condition 0 ≤ μ (x) + ν (x) ≤ 1. The μ and ν are called determining and intu- A A A A itionistic fuzzy set A on the universal set X, denote by {x; μ (x), ν (x)|x ∈X} we called μ A A A and ν are membership function and nonmembership function of A, respectively. μ (x) A A and ν (x) are called the membership degree and nonmembership degree of an element x belonging to A ⊆ X, respectively. IF(X) is called the set of the intuitionistic fuzzy set on the universal set X. Definition 2.3: An intuitionistic fuzzy number (IFN) A is (i) subset of the real line. (ii) convex for the membership function μ (x), that is, μ (αx + (1 − α)x ) ≥ min(μ (x ), A A 1 2 A 1 μ (x )) for all x , x ∈ R, α ∈ [0, 1] A 2 1 2 (iii) concave for the non-membership function ν (x), ν (αx + (1 − α)x ) ≤ max(ν (x ), A A 1 2 A 1 ν (x )) for all x , x ∈ R, α ∈ [0, 1] A 2 1 2 FUZZY INFORMATION AND ENGINEERING 107 Figure 1. ATrIFN A = (l, c, d, r); t , z . A A (iv) normal, that is, μ (x ) = 1, ν (x ) = 0 for some x ∈ R. A 0 A 0 0 Definition 2.4: ATrIFN A = (l, c, d, r); t , z is called GTrIFN, is shown Figure 1 if its mem- A A bership and nonmembership functions are defined as follows: 0if x < l t (x − l)/(c − l) if l ≤ x < c μ (x) = t if c ≤ x ≤ d (1) A A t (r − x)/(r − d) if d < x ≤ r 0if x > r And ⎪1if x < l [c − x + z (x − l)]/(c − l) if l ≤ x < c z if c ≤ x ≤ d ν (x) = (2) z (r − x) ⎪ A x − d + if d < x ≤ r ⎪ r − d 1if x > r respectively, where l ≤ c ≤ d ≤ r, the values t and z sub/sub are maximum mem- A A bership degree and minimum nonmembership degree of A, respectively, such that they satisfy the following condition: t ∈ [0, 1], z ∈ [0, 1] and t + z ∈ [0, 1]. A A A A Let π (x) = 1 − μ (x) − ν (x) (3) A A A π (x) is called the hesitancy degree of an element x ∈ A. It is the degree of indeterminacy membership of the element x to A. 108 D. HUNWISAI ET AL. From Definition 2.4, it is obvious that μ (x) + ν (x) = 1 for any x ∈ R if t = 1and A A A z = 0. Hence, the TrIFN A = (l, c, d, r); t , z degenerates to A = (l, c, d, r);1,0, which is a A A A trapezoidal fuzzy number (TrFN) [14]. Therefore, the concept of the TrIFN is generalisation of that of the TrFN. From A = (l, c, d, r); t , z if c = d = p then A = (l, p, r); t , z that is A = (l, p, r); t , z is a A A A A A A TIFN, which is particular case of TrIFN. Likewise to algebraic operations of TIFN and TrIFN are defined as follows. Definition 2.5: Let A = (l , c , d , r ); t , z and B = (l , c , d , r ); t , z be two GTrIFNs with 1 1 1 1 A A 2 2 2 2 B B t = t , z = z and γ = 0 be any real number. Then, the algebraic operations of GTrIFNs A B A B are defined as follows: A ⊕ B = (l + l , c + c , d + d , r + 2 ); t ∧ t , z ∨ z 1 2 1 2 1 2 1 2 A B A B A
B = (l − r , c − d , d − c , r − l ); t ∧ t , z ∨ z 1 2 1 2 1 2 1 2 A B A B (l l , c c , d d , r r ); t ∧ t , z ∨ z if A > 0, B > 0 ⎨ 1 2 1 2 1 2 1 2 A B A B A ⊗ B = (l r , c d , d c , r l ); t ∧ t , z ∨ z if A 0, B 0 1 2 1 2 1 2 1 2 A B A B (r r , d d , c c , l l ); t ∧ t , z ∨ z if A < 0, B < 0 1 2 1 2 1 2 1 2 A B A B (l /r , c /d , d /c , r /l ); t ∧ t , z ∨ z if A > 0, B > 0 ⎨ 1 2 1 2 1 2 1 2 A B A B A B = (r /r , d /d , c /c , l /l ); t ∧ t , z ∨ z if A 0, B 0 1 2 1 2 1 2 1 2 A B A B (r /l , d /c , c /d , l /r ); t ∧ t , z ∨ z if A < 0, B < 0 1 2 1 2 1 2 1 2 A B A B (γ l , γ c , γ d , γ r ); t , z if γ> 0 1 1 1 1 A A γ A = (γ r , γ d , γ c , γ l ); t , z if γ< 0 1 1 1 1 A A ∗−1 A = (1/r ,1/d ,1/c ,1/l ); t , z if A = 0 1 1 1 1 A A where the symbols and is the minimum operator and ∨ is the maximum operator. 2.2. Cut Sets of TrIFN Definition 2.6: [15]: A (α, λ)− cut set of A = (l, c, d, r); t , z is a crisp subset of R,which is A A defined as follows: ∗λ A ={x|μ (x) ≥ α, ν (x) ≤ λ} A A where 0 ≤ α ≤ t , z ≤ λ ≤ 1and 0 ≤ α + λ ≤ 1. A A Definition 2.7: [15]: The α− cut set and λ-cut set of A = (l, c, d, r); t , z are a crisp subset A A of R, which is defined as follows: A ={x|μ (x) ≥ α} α A and ∗λ A ={x|ν (x) ≤ λ} respectively. FUZZY INFORMATION AND ENGINEERING 109 Using the membership function of A = (l, c, d, r); t , z and Definition 2.7 such that A = A A α ∗λ {x|μ (x) ≥ α} and A ={x|ν (x) ≤ λ} are closed interval and calculated as follows: A A (t − α)l + αc (t − α)r + αd A A A = [L (α), R (α)] = , (10) α A A t t A A and (1 − λ)c + (λ − z )l (1 − λ)d + (λ − z )r A A A = [L (λ), R (λ)] = , (11) A A 1 − z 1 − z A A respectively. 3. Ranking of TrIFN This section briefly reviews the ambiguities and the accuracy function of a GTrIFN. Definition 3.1: Let A be an arbitrary IFN. The score function for the IFN A for membership and non-membership functions are denoted by M(μ ) and M(ν ),respectively. M(μ ) and A A A M(ν ) are defined by M(μ ) = [L (α) + R (α)]h(α)d(α) (12) A A A and M(ν ) = [L (λ) + R (λ)]g(λ)d(λ) (13) A A where h(α) and g(λ) satisfy the following conditions: (i) h(α)and g(λ) are monotonic increasing of α ∈ [0, t ] and monotonic decreasing of λ ∈ [z ,1]. (ii) h(α) ∈ [0, 1] and g(λ) ∈ [0, 1]. (iii) h(0) = 0and g(1) = 0. Let A be an arbitrary IFN. The ambiguities for IFN A for membership and nonmembership functions are denote by V (µ )and V (ν ), respectively. respectively. V(μ ) and V(ν ) are A A A A defined by V(μ ) = [L (α) + R (α)]h(α)d(α) (14) A A A and V(ν ) = [L (λ) + R (λ)]g(λ)d(λ) (15) A A Next, we find score, accuracy and ambiguities function of a GTrIFN. Let a GTrIFN A = (l, c, d, r); t , z the score function of a GTrIFN A for membership and A A non-membership functions can be written as follows: from Equations (10), (12) and h(α) = 110 D. HUNWISAI ET AL. α,weget l + 2c + 2d + r M(μ ) = t (16) Similarly, from Equations (11), (13) and g(λ) = λ,wehave l + 2c + 2d + r M(ν ) = (1 − z ) (17) A A The accuracy function of a GTrIFN A is denoted by M(μ ) + M(ν ) (l + 2c + 2d + r)t + (l + 2c + 2d + r)(1 − z ) A A A (A) = = (18) 2 12 from Equations (10), (14) and h(α) = α,weget r − l + 2d − 2c V(μ ) = t (19) Similarly, from Equations (11), (15) and g(λ) = λ,weget r − l + 2d − 2c V(ν ) = (1 − z ) . (20) A A The accuracy function of a GTrIFN A is denoted by 2 2 V(μ ) + V(ν ) (r − l + 2d − 2c)t + (r − l + 2d − 2c)(1 − z ) A A ∇(A) = = . (21) 2 12 Example 3.1: Let A = (155, 165, 175, 180); 0.7, 0.2and B = (130, 146, 150, 165); 0.6, 0.3 be two GTrIFNs then, (155 + 2(165) + 2(175) + 180)(0.7) M(μ ) = = 82.892 (155 + 2(165) + 2(175) + 180)(1 − 0.2) M(ν ) = = 108.267 82.892 + 108.267 ∴ (A) = = 95.580 (180 − 155 + 2(175) − 2(165))(0.7) V(μ ) = = 3.675 (180 − 155 + 2(175) − 2(165))(1 − 0.2) V(ν ) = = 4.8 3.675 + 4.8 ∴ ∇(A) = = 4.328 (130 + 2(146) + 2(150) + 165)(0.6) M(μ ) = = 8.87 6 FUZZY INFORMATION AND ENGINEERING 111 (130 + 2(146) + 2(150) + 165)(1 − 0.3) M(ν ) = = 72.438 8.87 + 72.438 ∴ (B) = = 40.65 (165 − 130 + 2(150) − 2(146)))(0.6) V(μ ) = = 2.58 (165 − 130 + 2(150) − 2(146)))(1 − 0.3) V(ν ) = = 3.512 2.58 + 3.512 ∴ ∇(B) = = 3.046 Theorem 3.1: Let A = (a , a , a , a ); t , z and B = (b , b , b , b ); t , z be GTrIFNs with 1 2 3 4 A A 1 2 3 4 B B t = t and z = z . The accuracy function : GIF(R) → R is a linear function. A B A B Proof: LetA ={(a , a , a , a ); t , z } and B ={(b , b , b , b ); t , z } then γ ≥ 0, β ≥ 0, 1 2 3 4 A A 1 2 3 4 B B we have Δ(γ + β ) A B = Δ{(γ a , γ a , γ a , γ a ); t , z }+ {(βb , βb , βb , βb ); t , z } 1 2 3 4 A A 1 2 3 4 B B = Δ{{(γ a + βb , γ a + βb , γ a + βb , γ a + βb ); t ∧ t , z ∨ z }} 1 1 2 2 3 3 4 4 A B A B = {{(γ a + βb ) + 2(γ a + βb ) + 2(γ a + βb ) 1 1 2 2 3 3 + (γ a + βb ))(t ∧ t ) }} + {{(γ a + βb ) + 2(γ a + βb ) 4 4 A B 1 1 2 2 + 2(γ a + βb ) + (γ a + βb ))(1 − (z ∨ z ) }} 3 3 4 4 A B = {{(γ a + 2γ a + 2γ a + γ a ) + (βb + 2βb + 2βb + βb )))(t ∧ t ) }} 1 2 3 4 1 2 3 4 A B + {{(γ a + 2γ a + 2γ a + γ a ) + (βb + 2βb + 2βb + βb )))(1 − (z ∨ z ) }} 1 2 3 4 1 2 3 4 A B 1 1 2 2 = γ (a + 2a + 2a + a )(t ) + β (b + 2b + 2b + b )(t ) 1 2 3 4 A 1 2 3 4 B 12 12 1 1 2 2 + γ (a + 2a + 2a + a )(1 − z ) + β (b + 2b + 2b + b )(1 − z ) 1 2 3 4 A 1 2 3 4 B 12 12 2 2 = γ (a + 2a + 2a + a )(t ) + (a + 2a + 2a + a )(1 − z ) 1 2 3 4 A 1 2 3 4 A 2 2 + β (b + 2b + 2b + b )(t ) + (b + 2b + 2b + b )(1 − z ) 1 2 3 4 B 1 2 3 4 B = γΔ(A) + βΔ(B). In the same way, if γ< 0,β< 0 we can prove (γ A + B) = γ(A) + (B). Therefore, is a linear function. 112 D. HUNWISAI ET AL. Theorem 3.2: LetA = (a , a , a , a ); t , z and B = (b , b , b , b ); t , z be GTrIFNs with 1 2 3 4 A A 1 2 3 4 B B t = t and z = z . The ambiguities function : GIF(R) → R is a linear function. A B A B (The rest of the proof is similar to proof of Theorem 3.1). Definition 3.2: Let A = (a , a , a , a ); t , z and B = (b , b , b , b ); t , z be GTrIFNs. The 1 2 3 4 A A 1 2 3 4 B B ranking order of AandB is stipulated as follows: (i) if A >B, then A > B (ii) if A <B, then A < B (iii) if A = B, then (iiia) if ∇(A) =∇(B), then A = B (iiib) if ∇(A)> ∇(B), then A < B (iiic) if ∇(A) ∇(B), then A B 4. Mathematical Formulation for IFTP This section, first introduces the mathematical formulation of the IFTP. Later, we find IBFS by NWCM and we use MODIM for finding optimal solution. The mathematical formulation of the IFTP is of the following form: m n (IFTP:1) Minimize = c x i,j i,j i=1 j=1 subject to x ≤ a , i ∈{1, 2, ... , m} ij i j=1 x ≥ b , j ∈{1, 2, ... , n} ij j i=1 x ≥ 0 for all i and j ij th th where c be GTrIFN cost of transportation one unit of the goods from i source to the j ij th th destination. x be the quantity transportation from i source to the j destination, is shown ij Table 1. th Here, a be the total availability of the goods at i source. th b be the total demand of the goods at j destination. m n x c be total intuitionistic fuzzy transportation cost. ij ij i=1 j=1 m n If a = b then IFTP is said to be balanced. i j i=1 j=1 m n If a = b then IFTP is said to be unbalanced (Table 1). i j i=1 j=1 From IFTP:1 can be written as the following linear programming problem (LPP): Minimize (LPP): Minimize (X) = C (X) subject to AX = b X ≥ 0, FUZZY INFORMATION AND ENGINEERING 113 Table 1. The intuitionistic fuzzy transportation table. 12 ... N a 1 c c ... c a 11 12 1n 1 2 c c ... c a 21 22 2n 2 ... ... ... ... ... ... m n b b b ... b a = b j 1 2 n i j i=1 j=1 where A be an m × n matrix, X be an n − vector, b be an m − vector, and c = (c , c , ... , c , ... , c , ... , c ) . 11 12 1n m1 mn Theorem 4.1: Let the intuitionistic fuzzy linear programming problem (IFLPP) be given as Minimize (X) = C (X) subject to AX = b (22) X ≥ 0, where A = (a ) = (A , A , ... , A ), c = (c , c , ... , c ) and b = (b ) , c , j = 1, 2, ... , n ij m×n 1 2 n 1 2 n i m×1 j are GTrIFNs. If for BFS X ,all T −1 = (c ) B A , j B j j=1 then X is optimal solution, where jare given by T −1 −1 j = (c ) B A and B A = ξ . j B j j j=1 Proof: We need to prove (X ) ≤ (Z).Let c = (c , c , ... , c ), B = (A , A , ... , A ), B B 1 2 m 1 2 m X = (x , x , ... , x ), (X ) = C (X ), where x (i = 1, 2, ... , m), is some x (j = 1, 2, ... , m),. B 1 2 m B B i j Let Z = (z , z , ... , z , ... , z ), any other feasible solution with z (i = 1, 2, ... , m, ... ,, n) 1 2 m n i some x (j = 1, 2, ... , m).Since B is basis, we have j j j A = ξ A + ξ A + ··· + ξ A , j ∈{1, 2, ... , n} (23) j 1 2 m m 1 2 Also, Z is a feasible solution. This refer to z A + z A +···+ z A = b (24) 1 1 2 2 n n from Equations (23)and (24), we get 1 2 n 1 2 n (z ξ + z ξ + ··· + z ξ )A + ··· + (z ξ + z ξ + ··· + z ξ )A = b 1 2 n 1 1 2 n m 1 1 1 m m m Since XB is a solution, that is x A + x A +···+ x A = b (26) 1 1 2 2 n n 114 D. HUNWISAI ET AL. Then Equations (25) and (26), together imply that x = z ξ , x (i = 1, 2, ... , m) i i i j=1 Since (c ) ≥ 0and is linear, therefore, j j (z)) = (c z ⊕ c z ⊕ ... c z ) 1 1 2 2 n n ≥( z ⊕ z ⊕ ... z ) 1 1 2 2 n n ⎛ ⎞ T j ⎝ ⎠ = (c) ξ z j=1 n m = c ξ z i j j=1 i=1 m n = (c ) z ξ i j i=1 j=1 = (c )x i i i=1 = c x i i i=1 = ( (X )) This implies that ( (X )) ≤ ( (Z)) and therefore (X ) ≤ (Z).So,X is optimal solu- B B B tion. The dual of the IFTP:1 can be written as m n Maximize = a u ⊕ b v (D) i i i j i=1 j=1 subject to u ⊕ v ≤ c , i ∈{1, 2, ... , m}; j ∈{1, 2, ... , n} i j ij That is Maximize = b Z (D) Subject to A Z ≤ c Z ≥ 0, where Z = (u , u , ... , v , v ... , v ) . 1 2 1 2 n 4.1. Algorithm to Find an Initial Basic Feasible Solution (IBFS) of IFTP In this section, we use intuitionistic fuzzy NWCM to compute IBFS of IFTP. Step 1: Set up the formulated intuitionistic fuzzy linear programming problem into the tab- ular form know as intuitionistic fuzzy transportation table (IFTT). An we approximate cost by GTrIFNs. FUZZY INFORMATION AND ENGINEERING 115 Step 2: Examine that the IFTP is balanced or unbalanced, if unbalanced, make it balanced. Step 3: Choose the north-west corner cell (NWCC) of the IFTT. Let it be the cell(i, j).Find x = min(a , b ). ij i j case (i) If a = min(a , b ), then allocate x = a in the (i, j)th cell of m × n IFTT. Delete the ith i i j ij i row to obtain a new IFTT of order (m − 1) × n. Replace b by b − a in obtained IFTT. j j i Go to step 4. case (ii) If b = min(a , b ), then allocate x = b in the (i, j)th cell of m × n IFTT. Delete the j i j ij j jth column to obtain a new allocate IFTT of order (m) × (n − 1).Replace a by a − b in i i j obtained IFTT. Go to step 4. case (iii) If a = b , then either follow case(i)or case(ii) but not both together. Go to step 4. i j Step 4: Calculate the penalties for the reduced IFTT obtain in step 3. Repeat step 3 until the IFTT is reduced to 1 × 1. Step 5: Allocate all x in the (i, j)th cell of the given IFTT. ij Step 6: The obtained IBFS and initial intuitionistic fuzzy transportation cost are x and ij m n x c respectively. ij ij i=1 j=1 4.2. Modified Distribution Method for Finding Optimal Solution In this section, we use generalised intuitionistic modified distribution method (GIMODIM) to find the optimal solution for IFTP. Algorithm of GIMODIM is illustrated as follows: Step 1: Find IBFS by propose IFNWCM. Step 2: Compute IF dual variables u and v for all row and column, respectively, sat- i j isfying (c ) = (u ⊕ v ) for all occupied cell. To start with. take any v or u as ij i j j i (−1, 0, 0, 1; 1, 0). Step 3: For unoccupied cells, find opportunity e = e , where = u ⊕ v . Step 4: ij ij ij ij i j Consider valued of (e ). ij case (i) IBFS is the intuitionistic fuzzy optimal solution, if (e ) ≥ 0 for all unoccupied cells. ij case (ii) IBFS is not the intuitionistic fuzzy optimal solution, for at least one (e )< 0. ij Go to step 5. Step 5: Choose the unoccupied cell for the most negative value of (e ). ij Step 6: We construct the closed loop below. At first, start the closed loop with choose the unoccupied cell and move vertically and hor- izontally with corner cells occupied and come back to choose the unoccupied cell to complete the loop. Use sign ‘+’and‘−’ at the corners of the closed loop, by assigning the ‘+’ sign to the selected unoccupied cell first. Step 7: Look for the least allocation value from the cells which have ‘−’sign. After that, allocate this value to the choose empty cell and subtract it to the other occupied cell having ‘−’ sign and add it to the other occupied cells having ‘+’ sign. Step 8: Allocation in Step 7 will result an improved basic feasible solution (BFS). Step 9: Test the optimality condition for improved BFS. The process is complete when (e ) ≥ 0 for all the unoccupied cell. ij 5. Numerical Example Next, we present some examples to illustrate our result. 116 D. HUNWISAI ET AL. Table 2. Data of the Example 5.1: the fuzzy transportation table. Source Phangnga Phuket Krabi Availability Si (3, 5, 7, 14); 0.6, 0.3 (2,4,8,13);0.7,0.2 (3,5,9,15);0.5,0.3 35 Yanok (2, 5, 8, 10); 0.8, 0.2 (3,6,9,12);0.5,0.4 (4,7,10,16);0.6,0.3 40 PhiPhi (3, 6, 8, 13); 0.8, 0.1 (4,8,10,15);0.6,0.2 (5,9,13,15);0.7,0.3 50 Demand 45 55 25 125 Table 3. The ﬁrst iteration choose the north-west corner cell. Source Phangnga Phuket Krabi Availability Si (3,5,7,14);0.6,0.3 35 (2,4,8,13);0.7,0.2 (3,5,9,15);0.5,0.3 35 Yanok (2,5,8,10);0.8,0.2 (3,6,9,12);0.5,0.4 (4,7,10,16);0.6,0.3 40 PhiPhi (3,6,8,13);0.8,0.1 (4,8,10,15);0.6,0.2 (5,9,13,15);0.7,0.3 50 Demand 45 55 25 125 Table 4. Last iteration by the north-west corner method. Source Phangnga Phuket Krabi availability Si (3,5,7,14);0.6,0.3 35 (2,4,8,13);0.7,0.2 (3,5,9,15);0.5,0.3 35 Yanok (2,5,8,10);0.8,0.2 10 (3,6,9,12);0.5,0.4 30 (4,7,10,16);0.6,0.3 40 PhiPhi (3,6,8,13);0.8,0.1 (4,8,10,15);0.6,0.2 25 (5,9,13,15);0.7,0.3 25 50 Demand 45 55 25 125 Example 5.1: Packing company a bird’s nest concession for nesting island three islands include Si, Yanok, and Phi Phi Island. Every week, the Bird’s Nest is transported to the three plants, which is located on the banks include Phangnga, Phuket and Krabi. Each island can collect nest up to 35, 40 and 50 kg, respectively. While, Phangnga, Phuket and Krabi were able to get a nest, cleaning and packing 45, 55 and 25 kg, respectively, shown in Table 2. For transportation costs from island to plant are as follows: (unit: 10 Baht per one kilogram of bird’s nest). From table 4, we will find out the minimum cost of total fuzzy transportation. 3 3 Since a = b = 125, the FTP is balanced. i j i=1 j=1 Finding IBFS of IFTP by IFNWCM. Now, transfer this allocation to the FTT. The first allocation is shown in Table 3 and the final allocation is shown in Table 4. Therefore, IBFS is x = 35, x = 10, x = 30, x = 25, x = 25, and total intuitionistic 11 21 22 32 33 fuzzy transportation cost is 35 (3, 5, 7, 14); 0.6, 0.3 ⊕ 10 (2, 5, 8, 10); 0.8, 0.2 ⊕ 30 (3, 6, 9, 12); 0.5, 0.4 ⊕ 25(4, 8, 10, 15); 0.6, 0.2⊕ 25((5.9.13.15) : 0.7.0.3) = ((440.830.1170.1700) : 0.5, 0.4) Now, we apply GIMODIM to compute the optimal solution. Algorithm of modified distribu- tion method as shown in Section 4.2. Firstly, we compute intuitionistic fuzzy dual variables u and v for each row and i j column, respectively, satisfying u ⊕ v = c for each occupied cell. Therefore, let v = i j ij 1 (−1, 0, 0, 1);1,0. FUZZY INFORMATION AND ENGINEERING 117 Table 5. Construction of loop. For each occupied cell,u ⊕ v = c we have i j ij c = u ⊕ v ; u = (2, 5, 7, 15); 0.6, 0.3 11 1 1 1 c = u ⊕ v ; u = (1, 5, 8, 11); 0.8, 0.2 21 2 1 2 c = u ⊕ v ; v = (−8, −2, 4, 11); 0.5, 0.4 22 2 2 2 c = u ⊕ v ; u = (−7, 4, 12, 23); 0.5, 0.4 32 3 2 3 c = u ⊕ v ; v = (−18, −3, 9, 22); 0.5, 0.4 . 33 3 3 3 Hence, we obtain e = c
(u ⊕ v ); = (−24, −7, 5, 19); 0.5, 0.3 12 12 1 2 e = c
(u ⊕ v ); = (−34, −11, 7, 31); 0.5, 0.4 13 13 1 3 e = c
(u ⊕ v ); = (−29, −10, 8, 33); 0.5, 0.4 . 23 23 2 3 e = c
(u ⊕ v ); = (−21, −6, 4, 21); 0.5, 0.4 31 31 3 1 From above, we found that the value of e is most negative, so IBFS is not intuitionistic fuzzy optimal. In Table 5, construct of loop. We use sign ‘+’in (1, 3)th cell, (2, 1)thcell and (3, 2)th cell. And use sign ‘−’in (1, 1)th cell, (2, 2)th cell and (3, 3)th cell. Check e again, if e ≥ 0 for all unoccupied cells, then the solution is intuitionistic ij ij fuzzy optimal solution. If e < 0, go to Step 5. ij Next, improved Basic Feasible Solution. Let v = (−1, 0, 0, 1);1,0. For each occupied cell, u ⊕ v = c , we compute i j ij c = u ⊕ v ; u = (2, 5, 7, 15); 0.6, 0.3 11 1 1 1 c = u ⊕ v ; v = (−13, −3, 3, 11); 0.6, 0.3 . 12 1 2 2 c = u ⊕ v ; v = (−12, −2, 4, 13); 0.5, 0.3 . 13 1 3 3 c = u ⊕ v ; u = (1, 5, 8, 11); 0.8, 0.2 21 2 1 2 c = u ⊕ v ; u = (−7, 5, 13, 28); 0.6, 0.3 32 3 2 3 118 D. HUNWISAI ET AL. Table 6. Improved basic feasible solution. Source Phangnga Phuket Krabi Availability Si (3,5,7,14);0.6,0.3 5 (2,4,8,13);0.7,0. 5 (3,5,9,15);0.5,0.3 25 35 Y anok (2,5,8,10);0.8,0.2 40 (3,6,9,12);0.5,0.4(4,7,10,16);0.6,0.340 PhiPhi (3,6,8,13);0.8,0.1(4,8,10,15);0.6,0.2 50 (5,9,13,15);0.7,0.350 Demand 45 55 25 125 Hence, we observe that e = c
(u ⊕ v ); = (−10, −3, 3, 10); 0.8, 0.2 21 21 2 1 e = c
(u ⊕ v ); = (−19, −5, 7, 24); 0.5, 0.4 22 22 2 2 e = c
(u ⊕ v ); = (−35, −8, 8, 35); 0.6, 0.3 32 32 3 2 e = c
(u ⊕ v ); = (−36, −8, 10, 34); 0.5, 0.3 . 33 33 3 3 From above, we found that the value of e ≥ 0 for all unoccupied cells, so optimal ij solution is x = 5, x = 5, x = 25, x = 40, x = 50 11 12 13 21 32 shown in Table 6, and the minimum transportation intuitionistic fuzzy cost is = (15, 25, 35, 70); 0.6, 0.3 ⊕ (10, 20, 40, 65); 0.7, 0.2 ⊕ (75, 125, 225, 375); 0.5, 0.3 ⊕ (80, 200, 320, 400); 0.8, 0.2 ⊕ (200, 400, 500, 750); 0.5, 0.3 = (380, 770, 1120, 1660); 0.5, 0.3 The minimum transportation intuitionistic fuzzy cost can be interpreted as follows: the min- imum transportation intutionistic fuzzy costs stay in the ranges [380, 1660], when (α, λ) = (0.5, 0.3). That is, the degree of acceptance of the transportation cost for the decision making increases if the cost increases from 380 to 770. The degree of acceptance of the transportation cost for the decision making is stationary when the costs are in the range 770–1120, while it decreases if the cost in increases from1120 to 1660. The transportation cost is totally acceptable if transportation cost stays in the ranges [770, 1120]. The degree of non-acceptance of the transportation cost for the decision making decreases if the cost increases from 380 to 770. The degree of un-acceptance of the transportation cost for the decision making is stationary when the costs are in the range 770–1120, while it increases if the cost increases from 1120 to 1660. 6. Conclusion In this paper, we are defined a new concept of linear ranking function for GTrIFNs. This new method is proposed to find the IBFS and the optimal solution of GTrIFNs based on the both demands and availabilities are real numbers. In addition, the cost is always GTrIFNs under the condition of the linear transportation problem. The advantages of this method can be used to solve for all kinds of IFTP, whether triangular fuzzy number, TrFN, TIFN, TrIFN or GTrIFN which this method is obtained solution is always optimal. Moreover, this method can use both the maximum and minimum values of an objective function. FUZZY INFORMATION AND ENGINEERING 119 However, this method has a limit for the linear multi-objective transportation problem and including other (nonlinear) shapes for membership functions, such as exponential membership function and hyperbolic membership function etc. Acknowledgements The authors would like to thank the referees for their esteemed comments and suggestions. Wiyada Kumam was financially supported by the Rajamangala University of Technology Thanyaburi (RMUTTT) [grant number NSF62D0604]. Disclosure statement No potential conflict of interest was reported by the author(s). Funding Darunee Hunwisai was financially supported by the Valaya Alongkorn Rajabhat University under the Royal Patronage (VRU) and King Mongkut’s University of Technology Thonburi (KMUTT). Notes on contributors DaruneeHunwisai was born in Bangkok, Thailand. She received a B.Ed. (Mathematics) degree from the Phranakhon Rajabhat University, Bangkok, Thailand, in 2000, the M.Ed. (Mathematics) degree from Phranakhon Rajabhat University, Thailand, in 2006 and the Ph.D. (Applied Mathematics) degree from King Mongkut’s University of Technology Thonburi, Thailand, in 2018. Currently, She is working at the Department of Mathematics and Statistics, Faculty of Science and Technology, Valaya Alongkorn Rajabhat University under the Royal Patronage. Her research interests are in the field of fuzzy fixed point, fuzzy mathematical models and optimization. Poom Kumam received the Ph.D. degree in mathematics from Naresuan University, Thailand. He is currently a Full Professor with the Department of Mathematics, King Mongkut’s University of Technol- ogy Thonburi (KMUTT). He is also the Head of the Fixed Point Theory and Applications Research Group, KMUTT, and also with the Theoretical and Computational Science Center (TaCS-Center), KMUTT. He is also the Director of the Computational and Applied Science for Smart Innovation Cluster (CLAS- SIC Research Cluster), KMUTT. He has successfully advised five master’s, and 38 Ph.D. graduates. His research targeted fixed point theory, variational analysis, random operator theory, optimization the- ory, and approximation theory. Also, fractional differential equations, differential game, entropy and quantum operators, fuzzy soft set, mathematical modeling for fluid dynamics and areas of inter- est inverse problems, dynamic games in economics, traffic network equilibria, bandwidth allocation problem, wireless sensor networks, image restoration, signal and image processing, game theory, and cryptology. He has provided and developed many mathematical tools in his fields productively over the past years. He has over 800 scientific articles and projects either presented or published. More- over, he is editorial board journals more than 50 journals and also he delivers many invited talks on different international conferences every year all around the world. Wiyada Kumam received the Ph.D. degree in applied mathematics from the King Mongkut’s Uni- versity of Technology Thonburi (KMUTT). She is currently an Associate Professor at the Program in Applied Statistics, Department of Mathematics and Computer Science, Faculty of Science and Tech- nology, Rajamangala University of Technology Thanyaburi (RMUTT). Her research interests include fuzzy optimization, fuzzy regression, fuzzy nonlinear mappings, leastsquares method, optimization problems, and image processing. 120 D. HUNWISAI ET AL. ORCID Poom Kumam http://orcid.org/0000-0002-5463-4581 Wiyada Kumam http://orcid.org/0000-0001-8773-4821 References [1] Zadeh LA. Fuzzy sets. Inf Control. 1965;8:338–353. [2] Atanassov KT. Intuitionistic fuzzy sets. Fuzzy Sets Syst. 1986;20(1):87–96. [3] Hitchcock FL. The distribution of a product from several sources to numerous localities. J Math Phys. 1941;20:224–230. [4] Dantzig GB. Application of the simplex method to a transportation problem, activity analysis of production and allocation. In: Koopmans TC, editor. New York: Wiley; 1951. p. 359–373. [5] Shanmugasundari M, Ganesan K. A novel approach for the fuzzy optimal solution of fuzzy transportion problem. Int J Eng Res Appl. 2013;3(1):1416–1421. [6] Liu P, Yang L, Wang L, et al. A solid transportation problem with type-2 fuzzy variables. Appl Soft Comput. 2014;24:543–558. [7] Giri PK, Maiti MK, Maiti M. Entropy based solid transportation problems with discounted unit costs under fuzzy random environment. OPSEARCH. 2014;51:479–532. [8] Basirzadeh H. An approach for solving fuzzy transportation problem. Appl Math Sci. 2011; 5(32):1549–1566. [9] Dubey D, Chandra S, Mehra A. Fuzzy linear programming under interval uncertainty based on IFS representation. Fuzzy Sets Syst. 2012;188(1):68–87. [10] Hussain RJ, Kumar PS. Algorithmic approach for solving intuitionistic fuzzy transportation prob- lem. Appl Math Sci. 2012;6(80):3981–3989. [11] Pramila K, Uthra G. Optimal solution of an intuitionistic fuzzy transportation problem. Ann Pure Appl Math. 2014;8(2):67–73. [12] Antony RJP, Savarimuthu SJ, Pathinathan T. Method for solving the transportation problem using triangular intuitionistic fuzzy number. Int J Comput Algorithm. 2014;3:590–605. [13] Singh SK, Yadav SP. A new approach for solving intuitionistic fuzzy transportation problem of type-2. Ann Oper Res. 2014: 1–15. [14] Dubois D, Prade H. Fuzzy set and systems theory and application. New York: Academic Press; [15] Nan JX, Li CF, Zhang MJ. A lexicographic method for matrix games wign payoffs of triangular intuitionistic fuzzy numbers. Int J Comput Intell Syst. 2010;3:280–289. [16] Li DF. Decision and game theory management with intuitionistic fuzzy sets. New York; 2014.
Journal
Fuzzy Information and Engineering
– Taylor & Francis
Published: Jan 2, 2019
Keywords: Intuitionistic fuzzy set; matrix game; linear programming; intuitionistic fuzzy transportation problem