Abstract
FUZZY INFORMATION AND ENGINEERING 2019, VOL. 11, NO. 1, 1–11 https://doi.org/10.1080/16168658.2021.1886809 Majid Amirfakhrian and S. R. Fayyaz Behrooz Department of Mathematics, Central Tehran Branch, Islamic Azad University, Tehran, Iran ABSTRACT ARTICLE HISTORY Received 24 December 2016 In this article, we try to solve a full fuzzy algebraic equation with a Revised 17 April 2017 fuzzy variable and fuzzy coefficients. This can be done by a numerical Accepted 6 April 2018 iterative process. We offer an algorithm to produce a sequence that may converge to a root of such an equation. KEYWORDS Approximation; algebraic fuzzy equation; nonlinear system 1. Introduction In many purposes of a system, we may need to solve a fuzzy equation. Zadeh introduced some properties of fuzzy equations in [1] and some other types of such equations were discussed by some researchers [2, 3] and they tried to solve them analytically. Solving an algebraic equation of degree at most 3 has no analytical method even in crisp form. So it is necessary to find the numerical solution of an equation. In [4], fuzzy nonlinear equa- tions with fuzzy unknown were considered. A research effort has been done on algebraic fuzzy equations by using neural network [5]. Two numerical iterative methods have been presented in [6, 7] to find the roots of an algebraic fuzzy equation of degree n, one with fuzzy coefficients and crisp unknown and the other one with crisp coefficients and fuzzy unknown. In the present essay, we introduce a numerical method to solve a full fuzzy alge- braic equation, such that all numbers are fuzzy, both the unknown variable and coefficients. The organisation of the present paper is as follows. In Section 2, we represent a full fuzzy algebraic equation and we discuss solving a nonlinear system of equations. In Section 3, an algorithm is introduced to numerically solving a full fuzzy algebraic fuzzy equation. Some examples are given in Section 4, and Section 5 contains the conclusion of this work. 2. Preliminaries The presentation of a fuzzy number f by a pair of functions as (f, f ) is called the parametric form of f, such that over [0, 1], we have: f is a left continuous and monotonically increasing function, f is a left continuous and monotonically decreasing function, and f ≤ f [8–12]. We consider the set of all fuzzy numbers by F. CONTACT Majid Amirfakhrian amirfakhrian@iauctb.ac.ir © 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. 2 M. AMIRFAKHRIAN AND S. R. FAYYAZ BEHROOZ Suppose f = (f, f ) and g ˜ = (g, g) arebothbelongingto F. Then some operations are defined as follows: (kf, xf ), k > 0, • kf = (kf, xf ), k < 0; • f + g ˜ = (f + g, f + g); • f − g ˜ = (f − g, f − g). • If w ˜ = f.g ˜, then w = min{f.g, f.g, f.g, f.g} w = max{f.g, f.g, f.g, f.g} ˜ ˜ A fuzzy number f is called non-negative (non-positive) if f ≥ 0 (f ≤ 0) on [0, 1]. Also f ∈ F is called positive (negative), if f > 0 (f < 0) on [0, 1]. It is obvious that for non-negative fuzzy ˜ ˜ ˜ numbers f and g ˜ we have f g ˜ = (f.g, f.g), and for non-positive fuzzy numbers f and g ˜ we have f g ˜ = (f.g, f.g). Definition 2.1: r−cut of a fuzzy number v ˜ is the interval [v(r), v(r)], and we use the notation [v ˜] for it. Definition 2.2: P (x ˜) is a fuzzy polynomial of degree at most n with fuzzy coefficients, if there are some fuzzy numbers a ˜ , ... , a ˜ ,suchthat 0 n P (x ˜) = a ˜ x ˜.(1) n j j=0 For a positive integer n, an algebraic equation with full fuzzy feature of degree n,can be defined as ˜ ˜ P (x ˜) = b,(2) where a ˜ , ... a ˜ , x ˜ ∈ F, a ˜ = 0and b ∈ F. 0 n n Definition 2.3: A‘m-degree polynomial form’ fuzzy number v ˜ is defined by a pair v ˜ = (p , q ),suchthat p and q are two polynomials of degree at most m [13–15]. m m m m We use PF for the set of all m-degree polynomial form fuzzy numbers [13, 14, 16]. Let s d T T F : R −→ R be a function such that maps Z = (y , ... , y ) to F (Z) = (f (Z), ... , f (Z)), 1 s 1 d where d ≥ s. To solve the system of equations F(Z) = 0, (3) ∂F (Z) by Gauss–Newton method [17], we consider A(Z) as the Jacobian of F: A(Z) = [ ], ∂y (0) j = 1, ... , s and i = 1, ... , d. Considering Z as an initial vector, we try to improve this (k) (k) (k) guess, so the system A(Z )H =−F(Z ) should be solved in each iteration and then the (k+1) (k) (k) approximated solution will be improved by considering Z = Z + H .Thissystemof FUZZY INFORMATION AND ENGINEERING 3 equation can be solved by a least square method as the following: (k) T (k) (k) (k) T (k) A(Z ) A(Z )H =−A(Z ) F(Z ).(4) Some conditions on convergence and uniqueness of the solution of (3) are proposed [18–20]. Let b = (b, b) and a ˜ = (a , a ) for j = 0, ... , n are known fuzzy numbers. We try j j to solve the following full fuzzy algebraic equation from degree n: a ˜ x ˜ + ··· + a ˜ x ˜ + a ˜ = b,(5) n 1 0 in which we have a ˜ = 0. We consider x ˜ = (x, x) as the unknown fuzzy root of Equation (5). 3. Solving the Algebraic Fuzzy Equation Let the unknown x ˜ be a fuzzy number belonging to PF . Also let we have a ˜ ∈ PF and j m b ∈ PF ,inwhich l ≤ nk + m. We discuss about Equation (5), by four cases: • First we consider the case that the unknown fuzzy number x ˜ be non-negative and for all j = 0, ... , n,wehave a ≥ 0. By these assumptions, Equation (5) changes to the two following equations: a x = b,(6) j=0 a x = b.(7) j=0 Now by the assumptions on fuzzy numbers, we consider that m m i i a = a r , a = a r , j ij j ij i=0 i=0 l l i i b = b r , b = b r , i=0 i=0 k k i i x = α r , x = β r . i i i=0 i=0 By these assumptions, Equation (6) changes to ⎧ ⎫ ⎛ ⎞ n m k l ⎨ ⎬ i i j i ⎝ ⎠ a r α r ) = b r.(8) ij i ⎩ ⎭ j=0 i=0 i=0 i=0 ˜ ˜ Since a ∈ PF , x ∈ PF ,wehave j m P (x ˜) = a ˜ x ˜ ∈ PF.(9) n j nk+m j=0 4 M. AMIRFAKHRIAN AND S. R. FAYYAZ BEHROOZ By rewriting (8) by an ordering on the power of r and by considering (9), for i = 0, ... , nk + m, the coefficient of r is a function L and we have nk+m l i i L (α , ... , α )r = b r . (10) i 0 k i=0 i=0 After taking the corresponding coefficients in both sides equal, for i = 0, ... , nk + m,we have an equation: b , i = 0, ... , l, L (α , ... , α ) = (11) i 0 k 0, i = l + 1, ... , nk + m. Similarly, by rewriting (7) by an ordering on the power of r and by considering (9), for i = 0, ... , nk + m, the coefficient of r is a function U and we have nk+m l i i U (β , ... , β )r = b r , (12) i 0 k i i=0 i=0 and by taking the corresponding coefficients in both sides equal, for i = 0, ... , nk + m, we have an equation: b , i = 0, ... , l, U (β , ... , β ) = (13) i 0 k 0, i = l + 1, ... , nk + m. ˜ ˜ • Let x be non-negative and for some j, a be negative. Thus Equation (5) changes to the two following equations: j j a x + a x = b, (14) j j a ≥0 a <0 j j j j a x + a x = b. (15) j j a ≥0 a <0 j j • In the case that a ˜ ’s are all non-negative and x ˜ is negative, Equation (5) changes to the following equations: j j a x + a x = b, (16) j even j odd j j a x + a x = b. (17) j even j odd ˜ ˜ • If x be negative and for some j, a be negative, then Equation (5) changes to the following equations: j j j j a x + a x + a x + a x = b (18) j j j j a ≥0 a ≥0 a <0 a <0 j j j j j even j even j odd j odd and j j j j a x + a x + a x + a x = b. (19) j j j j a ≥0 a ≥0 a <0 a <0 j j j j j even j even j odd j odd FUZZY INFORMATION AND ENGINEERING 5 In the last three cases, rewriting the equations by an ordering on the power of r leads us ˆ ˆ ˆ to introduce two functions L and U for i = 0, ... , nk + m, which are the coefficients of r (L i i i respect to (14), (16) or (18) and U respect to (15), (17) or (19)). For i = 0, ... , nk + m both ˆ ˆ functions L and U are functions of α ’s and β ’s. Thus we have i i j j b , i = 0, ... , l, L (α , ... , α , β , ... , β ) = (20) i 0 k 0 k 0, i = l + 1, ... , nk + m, and b , i = 0, ... , l, U (α , ... , α , β , ... , β ) = (21) i 0 0 k k 0, i = l + 1, ... , nk + m. Let us introduce two functions L and U as follows: i i L − b , i = 0, ... , l, L = (22) L , i = l + 1, ... , nk + m, and U − b , i = 0, ... , l, U = (23) U , i = l + 1, ... , nk + m. Therefore we have L = 0, i = 0, 1, ... , nk + m, (24) U = 0, i = 0, 1, ... , nk + m. Equation (24) is a nonlinear system of equations with d = 2(nk + m + 1) equations and s = 2k + 2 unknowns, so it can be solved by an iterative Gauss–Newton method. Suppose that t = nk + m Defining ⎛ ⎞ ∂L ∂L ∂L ∂L 0 0 0 0 ... ... ⎜ ⎟ ∂α ∂α ∂β ∂β 0 k 0 k ⎜ ⎟ . . . . . . ⎜ ⎟ . . . . . . . . ⎜ . . . . ⎟ ⎜ ⎟ ⎜ ⎟ ∂L ∂L ∂L ∂L t t t t ⎜ ⎟ ... ... ⎜ ⎟ ∂α ∂α ∂β ∂β 0 0 ⎜ k k ⎟ A = (25) ⎜ ⎟ ⎜ ⎟ ∂U ∂U ∂U ∂U 0 0 0 0 ⎜ ⎟ ... ... ⎜ ⎟ ∂α ∂α ∂β ∂β 0 0 k k ⎜ ⎟ ⎜ . . . . ⎟ . . . . . . . . ⎜ ⎟ . . . . . . ⎜ ⎟ ⎝ ⎠ ∂U ∂U ∂U ∂U t t t t ... ... ∂α ∂α ∂β ∂β 0 k 0 k and ⎛ ⎞ ⎛ ⎞ L α 0 0 ⎜ ⎟ ⎜ ⎟ ⎛ ⎞ . . ⎜ . ⎟ ⎜ . ⎟ . . ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ L α ⎜ t ⎟ ⎜ ⎟ ⎜ ⎟ B =− , Z = , H = , (26) ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ U ⎟ ⎜ β ⎟ ⎝ ⎠ 0 0 ⎜ ⎟ ⎜ ⎟ . . ⎜ ⎟ ⎜ ⎟ . . ⎝ ⎠ ⎝ ⎠ 2k+2 . . U β k 6 M. AMIRFAKHRIAN AND S. R. FAYYAZ BEHROOZ we have AH = B. (27) A and B areusedinstead of A(Z) and B(Z), respectively. By solving this system with an initial (0) (k) (0) vector Z , a sequence {Z } will be obtained as follows: we consider an initial vector Z . (k) (k) (k) (k) (k) For k = 1, 2, ..., we compute A and B and then we solve A H = B by a least square method: T T (k) (k) (k) (k) (k) A A H = A B , (28) and the obtaining vector will be improved by (k+1) (k) (k) Z = Z + H . (29) (k) (k) (k) (k) (k) Again A and B areusedinstead of A(Z ) and B(Z ), respectively. If {Z } converges ∗ ∗ to a vector Z ; then x is a solution of (5) where n n ∗ ∗ i ∗ ∗ i x = α r , x = β r . (30) i i i=0 i=0 Therefore the algorithm will be as follows: Algorithm 3.1: (0) (1) k = 0 and specify an initial vector Z . (2) Do the following steps until convergence condition is yield . (k) (k) (3) Compute A and B . T T (k) (k) (k) (k) (k) (k) (4) Solve linear equations A A H = A B and Compute H . (k+1) (k) (k) (5) Z = Z + H . (6) k = k + 1, and go to step 3. ∗ ∗ 1 Lemma 3.1: If the fuzzy Equation (5) has a root x ˜ then [x ˜ ] is a root of the equation 1 j 1 1 [a ˜ ] [x ˜ ] = [b] . j=1 Corollary 3.2: The algebraic fuzzy equation of degree n (5) has at most n fuzzy roots with distinct cores. Theorem 3.3: If the equation is crisp then the sequence is given by the following recurrence equation: a x + ... + a x + a − b n 1 ν 0 x = x − . (31) ν+1 ν n−1 na x + ... + 2a x + a n−1 ν 2 ν 1 Proof: For a crisp function algebraic equation P (x) = a x + ··· + a x + a = b,wehave n n 1 0 x = α and x = β , where α = β . For this equation, we have 0 0 0 0 L = a α + ··· + a α + a − b, 0 n 1 0 0 (32) U = a β + ··· + a β + a − b. 0 n 1 0 0 0 FUZZY INFORMATION AND ENGINEERING 7 The system (27) for this equation is as follows: ⎛ ⎞ dL ⎜ ⎟ dα h L 0 1 0 ⎜ ⎟ =− ⎝ ⎠ h U dU 2 0 dβ (ν) (ν) P (x ) n ν Using α = β = x,wehave L = U ,and h = h =− , therefore in the iteration 0 0 0 0 1 2 P (x ) n ν P (x ) n ν ν + 1, we have x = x − . ν+1 ν P (x ) By using the following lemma, we find that Algorithm 3.1 always generate a sequence. Lemma 3.4 ([21]): The system (28) always has a solution, for any k. 4. Numerical Examples In this section, we present some examples of full fuzzy algebraic equations. All examples are solved by commercial software Mathematica 10. Example 4.1: Consider the following equation: (r,2 − r)x ˜ + (2r,3 − r)x ˜ + (r,2 − r) = b, such that 3 2 3 2 b = (r + 4r + 4r, −r + 9r − 28r + 29). For this equation, we have n = 2, l = 3and m = 1. The least acceptable value of k is 1. By considering k = 1 and the initial fuzzy number x ˜ = (0, 0), after 6 iterations and 12 digits ∗ ∗ 1 2 ˜ ˜ accuracy, we obtain x = (r + 1, 3 − r). In this case, [x ] is the solution of x + 2x + 1 = 9. Example 4.2: For equation: (r,2 − r)x ˜ + (2r,3 − r)x ˜ + (r,2 − r) = b, with 5 3 3 2 b = (r + 2r + r, −4r + 22r − 43r + 29), we have n = 2, l = 5and m = 1. The least acceptable value of k is 2. By considering k = 2 and the initial fuzzy number x ˜ = (1, 1), after 6 iterations and 11 digits accuracy, we obtain ∗ 2 ∗ 1 2 x ˜ = (r ,3 − 2r). In this case, [x ˜ ] is the solution of x + 2x + 1 = 4. Example 4.3: Consider the following equation: ˜ ˜ (r,2 − r)x + (r − 3, −2r)x + (r,2 − r) = b, such that 3 2 3 2 b = (r + 3r + 12r − 12, −r + 8r − 37r + 34). For this equation, we have n = 2, l = 3and m = 1 and one of the coefficients is non- positive. The least acceptable value of k is 1. By considering k = 1 and the initial fuzzy 8 M. AMIRFAKHRIAN AND S. R. FAYYAZ BEHROOZ ˜ ˜ number x = (1, 1), after 7 iterations and 8 digits accuracy, we obtain x = (r + 2, 4 − r).In ∗ 1 2 this case, [x ˜ ] is the solution of x − 2x + 1 = 4. Example 4.4: Consider the following equation: (r + 1, 4 − 2r)x ˜ + (r,2 − r)x ˜ + (1, 1) = b, such that 3 2 3 2 b = (9r + 8r + 6r − 7, −2r + 17r − 64r + 65). For this equation, we have n = 2, l = 3and m = 1. The least acceptable value of k is 1. By considering k = 1 and the initial fuzzy number x ˜ = (−1, −1), after 7 iterations and 8 digits ∗ ∗ 1 2 accuracy, we obtain x ˜ = (r − 4, −3r). In this case, [x ˜ ] is the solution of 2x + x + 1 = 16. Here the root of equation is non-positive. Example 4.5: For the following equation: (r + 2, 5 − 2r)x ˜ + (2r − 6, −2r − 2)x ˜ + (2, 2) = b, with 3 2 3 2 b = (r + 4r + 2r + 2, −2r + 15r − 38r + 34), we have n = 2, l = 3and m = 1 and one of the coefficients is negative. The least accept- able value of k is 1. By considering k = 1 and the initial fuzzy number x ˜ = (−1, −1),after 8 ∗ ∗ 1 iterations and 8 digits accuracy, we obtain x ˜ = (r − 2, −r). In this case, [x ˜ ] is the solution of 3x − 4x + 2 = 9. Here the root of equation is non-positive. Example 4.6: Suppose that we want to solve the following equation: 4 3 (r,2 − r)x ˜ + (r,3 − 2r)x ˜ + (r,2 − r) 5 4 5 4 3 2 = (r + r + r, −r + 12r − 55r + 122r − 133r + 58). For this equation, we have n = 4, l = 5and m = 1. The least acceptable value of k is 1. By considering k = 1 and the initial fuzzy number x ˜ = (1, 1), after 10 iterations we obtain ∗ ∗ 1 4 3 x ˜ = (r,2 − r). In this case, [x ˜ ] is the solution of x + x + 1 = 3. Example 4.7: Suppose that a car is moving in a straight line with the velocity v ˜.Ifthiscar breaks, after a time t it will stop. The distance satisfies in the following relation: ˜ ˜ a ˜t + v ˜t = d, (33) in which a ˜ is the constant negative acceleration and t is the stopping time, and all variables m m are linguistic (Figure 1). If acceleration is about −6 and the velocity is about 30 , then by solving the equation, it will be clear that the car will stop after about 5 s (Figures 2 and 3). FUZZY INFORMATION AND ENGINEERING 9 Figure 1. Thetimeofbreaking. Figure 2. The coefficients of quadratic equation: (a) Velocity, (b) Acceleration and (c) Distance. Figure 3. Stopping time. Example 4.8: Consider the following fuzzy linear differential equation of order n: d x ˜ dx ˜ a ˜ + ··· + a ˜ + a ˜ x ˜ = 0, n 1 0 dt dt where t is the independent crisp variable. In order to find the solution of this ODE one must solve the following fuzzy algebraic equation of degree n: a ˜ w ˜ + ··· + a ˜ w ˜ + a ˜ = 0. n 1 0 Figure 4. The coefficients of quadratic equation: (a) a ˜ ,(b) a ˜ and (c) a ˜ . 2 1 0 10 M. AMIRFAKHRIAN AND S. R. FAYYAZ BEHROOZ Figure 5. Root of equation. For example by considering a ˜ , a ˜ and a ˜ as shown in Figure 4, the solution of quadratic 2 1 0 equation is the fuzzy number which is shown in Figure 5 and the solution of related differential equation is x ˜(t) = exp((1 + 0.35r, 1.6 − 0.25r)t) (Figures 4 and 5). 5. Conclusion In this work, we presented a method to find a numerical solution of a full fuzzy alge- braic equation. The method can be used for any kind of fuzzy equations, with non-positive or non-negative roots and with non-positive or non-negative coefficients. Some possible future research directions may consist of the more general case than the four cases which considered in this paper, and using the proposed for multiple roots of an equation. Disclosure statement No potential conflict of interest was reported by the authors. Notes on contributors Majid Amirfakhrian was born and raised in Tehran, Iran. He went to the University of Guilan, Rasht, Iran for his post–secondary education. He received his Master of Science in Applied Mathematics from the Sharif University of Technology and completed his Ph.D. through the Azad University, Science and Research branch. He joined the Azad University, Central Tehran Branch (IAUCTB) in 2002. He was promoted to a full professorship of Applied Mathematics at IAUCTB in 2015. His main research inter- ests lie in Numerical Analysis and Approximation, Fuzzy Approximation, Data Sciences, Multivariate Approximation, Image Processing, and Numerical Partial Differential Equations. He has side interests in Optimization and Statistics. S. R. Fayyaz Behrooz was born on November 22, 1978, and grown up in Tabriz, Iran. She received an M.S. degree from Tabriz University in 2004 and a B.S. degree from Islamic Azad University, Tabriz Branch in 2002. She is currently a Ph.D. student at the Islamic Azad University, Central Tehran Branch, Iran. She is a faculty member at Islamic Azad University, Osku Branch since 2008. Her main research interest includes Numerical and Computational Methods in Science and Fuzzy problems. ORCID Majid Amirfakhrian http://orcid.org/0000-0001-9934-6412 FUZZY INFORMATION AND ENGINEERING 11 References [1] Zadeh LA. Toward a generalized theory of uncertainty (GTU) an outline. Inform Sci. 2005;172:1–40. [2] Buckley JJ, Qu Y. Solving linear and quadratic fuzzy equations. Fuzzy Sets Syst. 1990;38:43–59. [3] Zhao R, Govind R. Solutions of algebraic equations involving generalized fuzzy numbers. Inform Sci. 1991;56:199–243. [4] Abbasbandy S, Asady B. Newton’s method for solving fuzzy nonlinear equations. Appl Math Comput. 2004;159:349–356. [5] Abbasbandy S, Otadi M. Numerical solution of fuzzy polynomials by a fuzzy neural network. Appl Math Comput. 2006;181:1084–1089. [6] Amirfakhrian M. Numerical solution of algebraic fuzzy equations with crisp variable by Gauss–Newton method. Appl Math Model. 2008;32:1859–1868. [7] Amirfakhrian M. An iterative Gauss–Newton method to solve an algebraic fuzzy equation with crisp coefficients. J Intel Fuzzy Syst. 2011;22:207–216. [8] Friedman M, Ma M, Kandel A. Fuzzy linear systems. Fuzzy Sets Syst. 1998;96:201–209. [9] Ma M, Friedman M, Kandel A. A new fuzzy arithmetic. Fuzzy Sets Syst. 1999;108:83–90. [10] Ma M, Kandel A, Friedman M. A new approach for defuzzification. Fuzzy Sets Syst. 2000;111:351–356. [11] Stefanini L, Sorini L, Guerra ML. Parametric representation of fuzzy numbers and application to fuzzy calculus. Fuzzy Sets Syst. 2006;157:2423–2455. [12] Goestschel R, Voxman W. Elementary calculus. Fuzzy Sets Syst. 1986;18:31–43. [13] Abbasbandy S, Amirfakhrian M. The nearest approximation of a fuzzy quantity in parametric form. Appl Math Comput. 2006;172:624–632. [14] Abbasbandy S, Amirfakhrian M. The nearest trapezoidal form of a generalized left right fuzzy number. Int J Approx Reason. 2006;43:166–178. [15] Amirfakhrian M. Some approximation methods in fuzzy logic: some method for finding approx- imation of fuzzy functions and fuzzy numbers. Germany: Lambert Academic Publishing; 2012. [16] Amirfakhrian M. Numerical solution of a system of polynomial parametric form fuzzy linear equations. Chapter 24, in Book ‘Ferroelectrics’, INTECH Publisher; 2010. [17] Deuflhard P. Newton methods for nonlinear problems. Berlin: Springer-Verlag Berlin Heidelberg; [18] Chen J, Li W. Convergence of Gauss–Newton’s method and uniqueness of the solution. Appl Math Comput. 2005;170:686–705. [19] Chen J, Li W. Local convergence results of Gauss–Newton’s like method in weak conditions. J Math Anal Appl. 2006;324:1381–1394. [20] Li C, Zhang WH, Jin XQ. Convergence and uniqueness properties of Gauss–Newton’s method. Comput Math Appl. 2004;47(6–7):1057–1067. [21] Gill PE, Murray W, Wright MH. Numerical linear algebra and optimization. Redwood City (CA): Addison–Wesley Publishing Company; 1991.
Journal
Fuzzy Information and Engineering
– Taylor & Francis
Published: Jan 2, 2019
Keywords: Approximation; algebraic fuzzy equation; nonlinear system