The Solution of Fully Fuzzy Quadratic Equation Based on Optimization

0 downloads 0 Views 2MB Size Report
May 17, 2014 - To this purpose we find and as optimum values which construct the best spreads. ... [12] obtained the real valued roots of fuzzy polynomials using ... numerical and neural net solutions to solve fuzzy equations. ... The set of all these fuzzy numbers is denoted by F. ...... Submit your manuscripts at.
Hindawi Publishing Corporation ξ€ e Scientific World Journal Volume 2014, Article ID 156203, 6 pages http://dx.doi.org/10.1155/2014/156203

Research Article The Solution of Fully Fuzzy Quadratic Equation Based on Optimization Theory T. Allahviranloo and L. Gerami Moazam Department of Mathematics, Islamic Azad University, Science and Research Branch, Tehran, Iran Correspondence should be addressed to T. Allahviranloo; [email protected] Received 16 April 2014; Accepted 17 May 2014; Published 9 June 2014 Academic Editor: Mohsen Vaez-ghasemi Copyright Β© 2014 T. Allahviranloo and L. Gerami Moazam. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Firstly in this paper we introduce a new concept of the 2nd power of a fuzzy number. It is exponent to production (EP) method that Μƒ 𝐢. Μƒ Μƒ = 𝐷, Μƒ where 𝐹(𝑋) Μƒ =𝐴 ̃𝑋 Μƒ2 + 𝐡̃𝑋+ provides an analytical and approximate solution for fully fuzzy quadratic equation (FFQE) : 𝐹(𝑋) To use the mentioned EP method, at first the 1-cut solution of FFQE as a real root is obtained and then unknown manipulated unsymmetrical spreads are allocated to the core point. To this purpose we find πœ† and πœ‡ as optimum values which construct the best spreads. Finally to illustrate easy application and rich behavior of EP method, several examples are given.

1. Introduction

2. Basic Concepts

The problem of finding the roots of equations like the quadratic equation has many applications in applied sciences like finance [1, 2], economy [3–6], and mechanics [7]. Sevastjanov and Dymova [8] proposed a new method for solving interval and fuzzy linear equations. In [9–11] Buckly discussed solving fuzzy equations. Abbasbandy and Otadi in [12] obtained the real valued roots of fuzzy polynomials using fuzzy neural networks. In [13–16] the authors have introduced numerical and neural net solutions to solve fuzzy equations. In the current paper we propose a new method to solve Μƒ+𝐢 Μƒ = 𝐷 Μƒ that the complication of arithmetics ̃𝑋 Μƒ2 + 𝐡̃𝑋 𝐴 Μƒ 𝐡, Μƒ and 𝑋 Μƒ being positive or negative. does not depend on 𝐴, This method solves some problems that have no analytical solution and this is an advantage of EP method because, as we know, numerical methods need initial guess to continue and the new method provides this neediness. The rest of the paper is set out as follows. In the second section some related basic definitions of fuzzy mathematics for the analysis are recalled. In Section 3, a new method, called EP, for solving fully fuzzy quadratic equation is presented. In Section 4, this method is used for an analytical approximate solution of FFQE. In Section 5, the conclusions are drawn.

The basic definitions are given as follows. Definition 1 (see [17–20]). A fuzzy number is a function 𝑒̃ : R β†’ [0, 1] which satisfies the following: (1) 𝑒̃ is upper semicontinuous on R; (2) 𝑒̃ is normal; that is, βˆƒπ‘₯0 ∈ R with 𝑒̃(π‘₯0 ) = 1; (3) 𝑒̃ is convex fuzzy set; (4) {π‘₯ ∈ R | 𝑒̃(π‘₯) > 0} is compact, where 𝐴 denotes the closure of 𝐴. Note: in this paper we consider fuzzy numbers which have a unique π‘₯0 ∈ R with 𝑒̃(π‘₯0 ) = 1 [18]. The set of all these fuzzy numbers is denoted by F. Obviously, R βŠ† F. For 0 < π‘Ÿ ≀ 1, we define π‘Ÿ-cut of 𝑒]0 = fuzzy number 𝑒̃ as [Μƒ 𝑒]π‘Ÿ = {π‘₯ ∈ R : 𝑒̃(π‘₯) β‰₯ π‘Ÿ} and [Μƒ {π‘₯ ∈ R : 𝑒̃(π‘₯) β‰₯ 0}. In [8] from (4)–(10) it follows that [Μƒ 𝑒]π‘Ÿ is a bounded closed interval for each π‘Ÿ ∈ [0, 1]. We denote the π‘Ÿ-cut of fuzzy number 𝑒̃ as [Μƒ 𝑒]π‘Ÿ = [𝑒(π‘Ÿ), 𝑒(π‘Ÿ)]. Definition 2 (see [18]). A fuzzy number 𝑒̃ is positive (negative) if 𝑒̃(π‘₯) = 0 for all π‘₯ < 0 (π‘₯ > 0).

2

The Scientific World Journal

Definition 3 (see [21, 22]). A fuzzy number 𝑒̃ in parametric form is a pair (𝑒, 𝑒) of functions 𝑒(π‘Ÿ) and 𝑒(π‘Ÿ), 0 ≀ π‘Ÿ ≀ 1, which satisfy the following requirements: (1) 𝑒(π‘Ÿ) is a bounded nondecreasing left continuous function in [0, 1]; (2) 𝑒(π‘Ÿ) is a bounded nonincreasing left continuous function in [0, 1]; (3) 𝑒(π‘Ÿ) ≀ 𝑒(π‘Ÿ), 0 ≀ π‘Ÿ ≀ 1. Definition 4 (see [21]). For arbitrary 𝑒̃ = (𝑒(π‘Ÿ), 𝑒(π‘Ÿ)) and ΜƒV = (V(π‘Ÿ), V(π‘Ÿ)), 0 ≀ π‘Ÿ ≀ 1, and scalar π‘˜, we define addition, subtraction, and scalar product by π‘˜ and multiplication is, respectively, as follows. Addition: 𝑒 + V(π‘Ÿ) = 𝑒(π‘Ÿ)+V(π‘Ÿ); 𝑒 + V(π‘Ÿ) = 𝑒(π‘Ÿ)+V(π‘Ÿ). Subtraction: 𝑒 βˆ’ V(π‘Ÿ) = 𝑒(π‘Ÿ) βˆ’ V(π‘Ÿ); 𝑒 βˆ’ V(π‘Ÿ) = 𝑒(π‘Ÿ) βˆ’ V(π‘Ÿ). Scalar product: (π‘˜π‘’ (π‘Ÿ) , π‘˜π‘’ (π‘Ÿ)) , π‘˜Μƒ 𝑒={ (π‘˜π‘’ (π‘Ÿ) , π‘˜π‘’ (π‘Ÿ)) ,

π‘˜β‰₯0 π‘˜ < 0.

(1)

Multiplication: 𝑒V (π‘Ÿ) = min {𝑒 (π‘Ÿ) V (π‘Ÿ) , 𝑒 (π‘Ÿ) V (π‘Ÿ) , 𝑒 (π‘Ÿ) V (π‘Ÿ) , 𝑒 (π‘Ÿ) V (π‘Ÿ)} 𝑒V (π‘Ÿ) = max {𝑒 (π‘Ÿ) V (π‘Ÿ) , 𝑒 (π‘Ÿ) V (π‘Ÿ) , 𝑒 (π‘Ÿ) V (π‘Ÿ) , 𝑒 (π‘Ÿ) V (π‘Ÿ)} . (2) For two important cases multiplication of two fuzzy numbers is defined by the following terms. If 𝑒̃ β‰₯ 0 and ΜƒV β‰₯ 0, then 𝑒V(π‘Ÿ) = 𝑒(π‘Ÿ)V(π‘Ÿ) and 𝑒V(π‘Ÿ) = 𝑒(π‘Ÿ)V(π‘Ÿ). If 𝑒̃ ≀ 0 and ΜƒV ≀ 0, then 𝑒V(π‘Ÿ) = 𝑒(π‘Ÿ)V(π‘Ÿ) and 𝑒V(π‘Ÿ) = 𝑒(π‘Ÿ)V(π‘Ÿ). If 𝑒̃ β‰₯ 0 and ΜƒV ≀ 0, then 𝑒V(π‘Ÿ) = 𝑒(π‘Ÿ)V(π‘Ÿ) and 𝑒V(π‘Ÿ) = 𝑒(π‘Ÿ)V(π‘Ÿ). If 𝑒̃ ≀ 0 and ΜƒV β‰₯ 0, then 𝑒V(π‘Ÿ) = 𝑒(π‘Ÿ)V(π‘Ÿ) and 𝑒V(π‘Ÿ) = 𝑒(π‘Ÿ)V(π‘Ÿ). Arithmetics of π‘Ÿ-cuts is similar to arithmetics of the parametric form recalled previously [23, 24]. Definition 5 (see [21]). Two fuzzy numbers 𝑒̃ and ΜƒV are said to be equal, if and only if 𝑒(π‘Ÿ) = V(π‘Ÿ) and 𝑒(π‘Ÿ) = V(π‘Ÿ), for each π‘Ÿ ∈ [0, 1]. A crisp number 𝛼 in parametric form is 𝑒(π‘Ÿ) = 𝑒(π‘Ÿ) = 𝛼, 0 ≀ π‘Ÿ ≀ 1. A triangular fuzzy number is popular and represented by 𝑒̃ = (π‘š, 𝛼, 𝛽), where 𝛼 > 0 and 𝛽 > 0, which has the parametric form as follows: 𝑒 (π‘Ÿ) = π‘š βˆ’ 𝛼 + π‘Ÿπ›Ό,

𝑒 (π‘Ÿ) = π‘š + 𝛽 βˆ’ π‘Ÿπ›½.

Μƒ, ΜƒV ⨁ 𝑀 Μƒ) = 𝐷(Μƒ Μƒ ∈ F; (1) 𝐷(Μƒ 𝑒⨁𝑀 𝑒, ΜƒV), for all 𝑒̃, ΜƒV, 𝑀 (2) 𝐷(π‘˜ ⨀ 𝑒̃, π‘˜ ⨀ ΜƒV) = |π‘˜|𝐷(Μƒ 𝑒, ΜƒV), for all π‘˜ ∈ R and 𝑒̃, ΜƒV ∈ F; Μƒ ⨁ 𝑒̃) ≀ 𝐷(Μƒ Μƒ) + 𝐷(ΜƒV, 𝑒̃), for all 𝑒̃, ΜƒV, 𝑀 Μƒ, (3) 𝐷(Μƒ 𝑒 ⨁ ΜƒV, 𝑀 𝑒, 𝑀 𝑒̃ ∈ F. Therefore (F, 𝐷) is a complete metric space.

3. Exponent to Production Method Μƒ = (π‘Ž(π‘Ÿ), π‘Ž(π‘Ÿ)), 𝐡̃ = (𝑏(π‘Ÿ), 𝑏(π‘Ÿ)), 𝐢 Μƒ = (𝑐(π‘Ÿ), 𝑐(π‘Ÿ)), 𝐷 Μƒ = Let 𝐴 Μƒ (𝑑(π‘Ÿ), 𝑑(π‘Ÿ)), 𝑋 = (π‘₯(π‘Ÿ), π‘₯(π‘Ÿ)) and Μƒ = 𝐷, Μƒ 𝐹 (𝑋)

(4)

Μƒ = (𝐹 (π‘₯, π‘₯, π‘Ÿ) , 𝐹 (π‘₯, π‘₯, π‘Ÿ)) , 𝐹 (𝑋)

(5)

Μƒ =𝐴 ̃𝑋 Μƒ2 + 𝐡̃𝑋 Μƒ + 𝐢. Μƒ In this method, we convert the where 𝐹(𝑋) 2nd exponent of a fuzzy number to a product of two fuzzy numbers in parametric form. By this conversion we obtain an analytical and approximate solution for a FFQE. In EP method, at first, we find 𝛼, as a real root of crisp 1-cut equation Μƒ = (π›Όβˆ’π‘“1 (𝛼1 (π‘Ÿ), 𝛼2 (π‘Ÿ)), 𝛼+𝑓2 (𝛼1 (π‘Ÿ), 𝛼2 (π‘Ÿ))) and then we get 𝑋 Μƒ2 β‰… (𝛼 βˆ’ as solution of (4) and apply the approximation 𝑋 𝛼1 (π‘Ÿ), 𝛼 + 𝛼1 (π‘Ÿ))(𝛼 βˆ’ 𝛼2 (π‘Ÿ), 𝛼 + 𝛼2 (π‘Ÿ)), in which 𝛼1 (π‘Ÿ) β©Ύ 0, 𝛼2 (π‘Ÿ) β©Ύ 0, 𝛼1σΈ€  (π‘Ÿ) < 0, and 𝛼2σΈ€  (π‘Ÿ) < 0. To find 𝑓𝑖 (𝛼1 (π‘Ÿ), 𝛼2 (π‘Ÿ)), for 𝑖 = 1, 2, we consider two cases: (1) FFQE has analytical solution; (2) FFQE does not have analytical solution. 3.1. Case (1). In this case we have π‘₯(0) β©½ π‘₯(1) = π‘₯(1) β©½ π‘₯(0). Therefore we construct conditions that provide a good approximation. These conditions are as follows: (1) 𝑓𝑖 (𝛼1 (1), 𝛼2 (1)) = 0, for 𝑖 = 1, 2; (2) 𝛼 βˆ’ 𝑓1 (𝛼1 (0), 𝛼2 (0)) = π‘₯(0) and 𝛼 + 𝑓2 (𝛼1 (0), 𝛼2 (0)) = π‘₯(0). Μƒ = (𝛼 βˆ’ To find 𝑓𝑖 (𝛼1 (π‘Ÿ), 𝛼2 (π‘Ÿ)), 𝑖 = 1, 2, at first we substitute 𝑋 2 Μƒ 𝑓1 (𝛼1 (π‘Ÿ), 𝛼2 (π‘Ÿ)), 𝛼 + 𝑓2 (𝛼1 (π‘Ÿ), 𝛼2 (π‘Ÿ))) and 𝑋 β‰… (𝛼 βˆ’ 𝛼1 (π‘Ÿ), 𝛼 + 𝛼1 (π‘Ÿ))(𝛼 βˆ’ 𝛼2 (π‘Ÿ), 𝛼 + 𝛼2 (π‘Ÿ)), in parametric form of (4), and then we obtain (π‘Ž (π‘Ÿ) , π‘Ž (π‘Ÿ)) (𝛼 βˆ’ 𝛼1 (π‘Ÿ) , 𝛼 + 𝛼1 (π‘Ÿ)) (𝛼 βˆ’ 𝛼2 (π‘Ÿ) , 𝛼 + 𝛼2 (π‘Ÿ)) + (𝑏 (π‘Ÿ) , 𝑏 (π‘Ÿ)) Γ— (𝛼 βˆ’ 𝑓1 (𝛼1 (π‘Ÿ) , 𝛼2 (π‘Ÿ)) , 𝛼 + 𝑓2 (𝛼1 (π‘Ÿ) , 𝛼2 (π‘Ÿ))) + (𝑐 (π‘Ÿ) , 𝑐 (π‘Ÿ)) = (𝑑 (π‘Ÿ) , 𝑑 (π‘Ÿ)) .

(3)

Definition 6 (see [25]). Let 𝐷 : F Γ— F β†’ R βˆͺ {0} and let 𝐷(Μƒ 𝑒, ΜƒV) = supπ‘Ÿβˆˆ[0,1] max{|𝑒(π‘Ÿ) βˆ’ V(π‘Ÿ)|, |𝑒(π‘Ÿ) βˆ’ V(π‘Ÿ)|} be the Hausdorff distance between fuzzy numbers, where [Μƒ 𝑒]π‘Ÿ = [𝑒(π‘Ÿ), 𝑒(π‘Ÿ)] and [ΜƒV]π‘Ÿ = [V(π‘Ÿ), V(π‘Ÿ)]. The following properties are well known:

(6) Set π‘Ÿ = 1. Using condition (1) and π‘Ž(1) = π‘Ž(1), 𝑏(1) = 𝑏(1), 𝑐(1) = 𝑐(1), and 𝑑(1) = 𝑑(1), we obtain 𝛼 (𝛼1 (1) + 𝛼2 (1)) = 0,

(7)

The Scientific World Journal

3 Μƒ < 0 and 𝐡̃ < 0, then (4) if 𝐴

for each 𝛼 ∈ R; that means 𝛼1 (1) + 𝛼2 (1) = 0.

(8)

πœ’ (π‘Ÿ) = βˆ’

π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) βˆ’ π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) . π‘Ž (π‘Ÿ) π‘ˆπ‘ˆ+ (πœ†, π‘Ÿ) + π‘Ž (π‘Ÿ) 𝐿𝐿+ (πœ‡, π‘Ÿ)

(16)

By this conclusion we decide to get Μƒ 0,

(9)

Μƒ > 0 and 𝐡̃ < 0, then (1) if 𝐴

πœ†, πœ‡ ∈ R.

πœ’ (π‘Ÿ) = βˆ’

Using condition (2) we have 𝛼 βˆ’ πœ† (𝛼1 (0) + 𝛼2 (0)) = π‘₯ (0) , 𝛼 + πœ‡ (𝛼1 (0) + 𝛼2 (0)) = π‘₯ (0) .

(10)

πœ’ (π‘Ÿ) = βˆ’

+ (𝑏 (π‘Ÿ) , 𝑏 (π‘Ÿ)) (𝛼 βˆ’ πœ† (𝛼1 + 𝛼2 ) (π‘Ÿ) , 𝛼 + πœ‡ (𝛼1 + 𝛼2 ) (π‘Ÿ))

πœ’ (π‘Ÿ) =

πœ’ (π‘Ÿ) = (11)

(12)

as a solution of (4). Notice that always we have real spreads, which means 𝛼1 (π‘Ÿ) + 𝛼2 (π‘Ÿ) ∈ R, because 𝛼1 (π‘Ÿ) and 𝛼2 (π‘Ÿ) are the roots of 𝑍2 βˆ’ πœ’(π‘Ÿ)𝑍 + ](π‘Ÿ) = 0. Using the proposed method we obtain the following set of expressions for πœ’(π‘Ÿ).

π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) βˆ’ π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) ; π‘Ž (π‘Ÿ) π‘ˆπΏ+ (πœ‡, π‘Ÿ) + π‘Ž (π‘Ÿ) πΏπ‘ˆ+ (πœ†, π‘Ÿ)

(19)

π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) βˆ’ π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) ; π‘Ž (π‘Ÿ) π‘ˆπΏβˆ’ (πœ†, π‘Ÿ) + π‘Ž (π‘Ÿ) πΏπ‘ˆβˆ’ (πœ‡, π‘Ÿ)

(20)

where 𝐺 (π‘₯, π‘₯, π‘Ÿ) = 𝐹 (π‘₯, π‘₯, π‘Ÿ) βˆ’ 𝑑 (π‘Ÿ) , 𝐺 (π‘₯, π‘₯, π‘Ÿ) = 𝐹 (π‘₯, π‘₯, π‘Ÿ) βˆ’ 𝑑 (π‘Ÿ) , π‘ˆπ‘ˆβ€  (πœ‰, π‘Ÿ) = π›Όπ‘Ž (π‘Ÿ) β€ πœ‰π‘ (π‘Ÿ) , 𝐿𝐿† (πœ‰, π‘Ÿ) = π›Όπ‘Ž (π‘Ÿ) β€ πœ‰π‘ (π‘Ÿ) ,

(21)

π‘ˆπΏβ€  (πœ‰, π‘Ÿ) = π›Όπ‘Ž (π‘Ÿ) β€ πœ‰π‘ (π‘Ÿ) , πΏπ‘ˆβ€  (πœ‰, π‘Ÿ) = π›Όπ‘Ž (π‘Ÿ) β€ πœ‰π‘ (π‘Ÿ) , in which † ∈ {+, βˆ’} and πœ‰ ∈ {πœ†, πœ‡}.

Μƒ>0: Case (1). 𝑋 Μƒ > 0 and 𝐡̃ < 0, then (1) if 𝐴 π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) βˆ’ π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) ; π‘Ž (π‘Ÿ) π‘ˆπ‘ˆβˆ’ (πœ†, π‘Ÿ) + π‘Ž (π‘Ÿ) πΏπΏβˆ’ (πœ‡, π‘Ÿ)

(13)

π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) βˆ’ π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) ; π‘Ž (π‘Ÿ) π‘ˆπ‘ˆ+ (πœ‡, π‘Ÿ) + π‘Ž (π‘Ÿ) 𝐿𝐿+ (πœ†, π‘Ÿ)

(14)

π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) βˆ’ π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) ; π‘Ž (π‘Ÿ) π‘ˆπ‘ˆβˆ’ (πœ‡, π‘Ÿ) + π‘Ž (π‘Ÿ) πΏπΏβˆ’ (πœ†, π‘Ÿ)

Μƒ > 0. Let Proof. Without loss of generality suppose that 𝑋 Μƒ = (𝛼 βˆ’ 𝛼1 (π‘Ÿ) , 𝛼 + 𝛼1 (π‘Ÿ)) (𝛼 βˆ’ 𝛼2 (π‘Ÿ) , 𝛼 + 𝛼2 (π‘Ÿ)) , ̃𝑍 π‘Œ

Μƒ < 0 and 𝐡̃ > 0, then (3) if 𝐴 πœ’ (π‘Ÿ) = βˆ’

3.2. Case (2). In this case we do not have analytical solution, and we do not have π‘₯(0) and π‘₯(0) in which π‘₯(0) β©½ π‘₯(1) = π‘₯(1) β©½ π‘₯(0); therefore we propose that πœ† = πœ‡ = 0.5 because of Lemma 7. Μƒ =𝐷 Μƒ and we do not Lemma 7. If 𝛼 is real core point of 𝐹(𝑋) have analytical solution, then πœ† = πœ‡ = 0.5 is the best choice in EP method.

Μƒ > 0 and 𝐡̃ > 0, then (2) if 𝐴 πœ’ (π‘Ÿ) =

(18)

Μƒ < 0 and 𝐡̃ < 0, then (4) if 𝐴

+ (𝑐 (π‘Ÿ) , 𝑐 (π‘Ÿ)) = (𝑑 (π‘Ÿ) , 𝑑 (π‘Ÿ)) . Now we can find πœ†, πœ‡, and πœ’(0). To construct solution of the new method, in the above parametric form, let πœ’(π‘Ÿ) = 𝛼1 (π‘Ÿ) + 𝛼2 (π‘Ÿ) and ](π‘Ÿ) = 𝛼1 (π‘Ÿ) Γ— 𝛼2 (π‘Ÿ); then, by solving a 2 Γ— 2 system via πœ’(π‘Ÿ) and ](π‘Ÿ), we obtain πœ’(π‘Ÿ) and we set

π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) βˆ’ π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) ; π‘Ž (π‘Ÿ) π‘ˆπΏβˆ’ (πœ‡, π‘Ÿ) + π‘Ž (π‘Ÿ) πΏπ‘ˆβˆ’ (πœ†, π‘Ÿ)

Μƒ < 0 and 𝐡̃ > 0, then (3) if 𝐴

(π‘Ž (π‘Ÿ) , π‘Ž (π‘Ÿ)) (𝛼 βˆ’ 𝛼1 (π‘Ÿ) , 𝛼 + 𝛼1 (π‘Ÿ)) (𝛼 βˆ’ 𝛼2 (π‘Ÿ) , 𝛼 + 𝛼2 (π‘Ÿ))

πœ’ (π‘Ÿ) =

(17)

Μƒ > 0 and 𝐡̃ > 0, then (2) if 𝐴

Up to now we have two equations and three unknown πœ†, πœ‡, and πœ’(0) = 𝛼1 (0) + 𝛼2 (0). The third equation comes from the equality below, for π‘Ÿ = 0,

Μƒ = (𝛼 βˆ’ πœ†πœ’ (π‘Ÿ) , 𝛼 + πœ‡πœ’ (π‘Ÿ)) 𝑋

π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) βˆ’ π‘Ž (π‘Ÿ) 𝐺 (𝛼, 𝛼, π‘Ÿ) ; π‘Ž (π‘Ÿ) π‘ˆπΏ+ (πœ†, π‘Ÿ) + π‘Ž (π‘Ÿ) πΏπ‘ˆ+ (πœ‡, π‘Ÿ)

Μƒ = (𝛼 βˆ’ πœ† (𝛼1 (π‘Ÿ) + 𝛼2 (π‘Ÿ)) , 𝛼 + πœ‡ (𝛼1 (π‘Ÿ) + 𝛼2 (π‘Ÿ))) . 𝑋 (22) (15)

Μƒ Μƒ2 β‰… π‘Œ Μƒ 𝑍. In EP method we use the approximation 𝑋

4

The Scientific World Journal

To find the optimum parameters πœ† and πœ‡, we must solve the minimization problem as follows: Min

Μƒ Μƒ 𝑍) Μƒ2 , π‘Œ 𝐷 (𝑋

st.

πœ†, πœ‡ ∈ R,

πœ† > 0,

(23)

Μƒ < 0 and 𝐡̃ > 0, then (3) if 𝐴 𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) β‰₯ βˆ’ [π‘ˆπ‘ˆβˆ’ (0.5, 0) + πΏπΏβˆ’ (0.5, 0)] πœ’ (0) ;

(30)

Μƒ < 0 and 𝐡̃ < 0, then (4) if 𝐴

πœ‡ > 0,

Μƒ = sup Μƒ2 , π‘Œ Μƒ 𝑍) where 𝐷(𝑋 π‘Ÿβˆˆ[0,1] max{|𝛼(2πœ† βˆ’ 1)πœ’(π‘Ÿ) + ](π‘Ÿ) βˆ’ 2 2 πœ† πœ’ (π‘Ÿ)|, |𝛼(1 βˆ’ 2πœ‡)πœ’(π‘Ÿ) + ](π‘Ÿ) βˆ’ πœ‡2 πœ’2 (π‘Ÿ)|}. To delimitate maximum error for any 𝛼 ∈ R, we choose πœ† = πœ‡ = 0.5. This completes the proof. Lemma 8. Necessary condition for existence of EP solution with πœ† = πœ‡ = 0.5 is π‘Ž (0) 𝐺 (𝛼, 𝛼, 0) βˆ’ π‘Ž (0) 𝐺 (𝛼, 𝛼, 0) β‰₯ 0,

Μƒ > 0, if 𝐴

(24)

π‘Ž (0) 𝐺 (𝛼, 𝛼, 0) βˆ’ π‘Ž (0) 𝐺 (𝛼, 𝛼, 0) ≀ 0,

Μƒ < 0. if 𝐴

(25)

Proof. Since πœ’(π‘Ÿ) is sum of 𝛼1 (π‘Ÿ) and 𝛼2 (π‘Ÿ) and we want 𝛼1 (π‘Ÿ) and 𝛼2 (π‘Ÿ) to be nonnegative, then necessary condition for existence of EP solution with πœ† = πœ‡ = 0.5 is πœ’(π‘Ÿ) β‰₯ 0 and since [𝛼 βˆ’ 0.5πœ’(π‘Ÿ), 𝛼 + 0.5πœ’(π‘Ÿ)] βŠ† [𝛼 βˆ’ 0.5πœ’(0), 𝛼 + 0.5πœ’(0)], for all 0 ≀ π‘Ÿ ≀ 1, it is sufficient to have πœ’(0) β‰₯ 0. Considering denominators of (13)–(20) we find that (24) and (25) hold if πœ’(0) β‰₯ 0, and this completes the proof. Lemma 9. The EP method with πœ† = πœ‡ = 0.5 does not have solution, if π‘Ž (0) 𝐺 (𝛼, 𝛼, 0) βˆ’ π‘Ž (0) 𝐺 (𝛼, 𝛼, 0) < 0,

Μƒ > 0, when 𝐴 (26)

𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) β‰₯ βˆ’ [π‘ˆπ‘ˆ+ (0.5, 0) + 𝐿𝐿+ (0.5, 0)] πœ’ (0) ,

(31)

Μƒ < 0; with 𝑋 Μƒ > 0 and 𝐡̃ < 0, then (1) if 𝐴 𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) β‰₯ βˆ’ [πΏπ‘ˆ+ (0.5, 0) + π‘ˆπΏ+ (0.5, 0)] πœ’ (0) ;

(32)

Μƒ > 0 and 𝐡̃ > 0, then (2) if 𝐴 𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) β‰₯ βˆ’ [πΏπ‘ˆβˆ’ (0.5, 0) + π‘ˆπΏβˆ’ (0.5, 0)] πœ’ (0) ;

(33)

Μƒ < 0 and 𝐡̃ > 0, then (3) if 𝐴 𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) β‰₯ [πΏπ‘ˆ+ (0.5, 0) + π‘ˆπΏ+ (0.5, 0)] πœ’ (0) ; (34) Μƒ < 0 and 𝐡̃ < 0, then (4) if 𝐴

Μƒ < 0. when 𝐴 (27)

𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) β‰₯ [πΏπ‘ˆβˆ’ (0.5, 0) + π‘ˆπΏβˆ’ (0.5, 0)] πœ’ (0) . (35)

Lemma 10. Suppose πœ† = πœ‡ = 0.5 in EP method and πœ’(0) β‰₯ 0, then sufficient condition for existence of solution is ](0) β‰₯ 0.

Μƒ> Proof. We consider only one case to discuss. We consider 𝐴 Μƒ Μƒ 0, 𝐡 < 0, and 𝑋 > 0 and solve (11) via πœ’(π‘Ÿ) and ](π‘Ÿ) in π‘Ÿ = 0; we obtain

if π‘Ž (0) 𝐺 (𝛼, 𝛼, 0) βˆ’ π‘Ž (0) 𝐺 (𝛼, 𝛼, 0) > 0, Proof. This lemma is conclusion of Lemma 8.

Proof. We know 𝛼1 (π‘Ÿ) and 𝛼2 (π‘Ÿ) are nonnegative if πœ’(π‘Ÿ) β‰₯ 0 and ](π‘Ÿ) β‰₯ 0. Because of construction of EP method for πœ† = πœ‡ = 0.5 and Lemma 8, it is obvious that we must have ](0) β‰₯ 0. Lemma 11. Suppose πœ’(0) β‰₯ 0 and πœ† = πœ‡ = 0.5 in EP method, Μƒ > 0 as follows: then we have solution with 𝑋 Μƒ > 0 and 𝐡̃ < 0, then (1) if 𝐴 𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) β‰₯ [π‘ˆπ‘ˆβˆ’ (0.5, 0) + πΏπΏβˆ’ (0.5, 0)] πœ’ (0) ; (28) Μƒ > 0 and 𝐡̃ > 0, then (2) if 𝐴 𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) β‰₯ [π‘ˆπ‘ˆ+ (0.5, 0) + 𝐿𝐿+ (0.5, 0)] πœ’ (0) ; (29)

] (0) = (𝐺 (𝛼, 𝛼, 0) βˆ’ 𝐺 (𝛼, 𝛼, 0) βˆ’ (π‘ˆπ‘ˆβˆ’ (0.5, 0) + πΏπΏβˆ’ (0.5, 0)) πœ’ (0) )

(36)

βˆ’1

Γ— [π‘Ž (0) βˆ’ π‘Ž (0)] . This is obvious that ](0) β‰₯ 0, if the numerator is nonnegative and this completes the proof.

4. Numerical Examples In the next examples we use round numbers with approximation less than 10βˆ’4 . Μƒ = (4/1/1), 𝐡̃ = (2/1/1), 𝐢 Μƒ = (1/1/1), and Example 1. Let 𝐴 Μƒ 𝐷 = (3/2/2) [9]. Μƒ β‰₯ 0. We will look for a solution where 𝑋

The Scientific World Journal

5

Μƒ = 𝐷 Μƒ becomes in parametric form as Equation 𝐹(𝑋) follows: (3 + π‘Ÿ, 5 βˆ’ π‘Ÿ) (π‘₯2 (π‘Ÿ) , π‘₯2 (π‘Ÿ)) + (1 + π‘Ÿ, 3 βˆ’ π‘Ÿ) (π‘₯ (π‘Ÿ) , π‘₯ (π‘Ÿ)) + (π‘Ÿ, 2 βˆ’ π‘Ÿ) = (1 + 2π‘Ÿ, 5 βˆ’ 2π‘Ÿ) . (37) The real roots of 1-cut and 0-cut equations are 𝛼 = π‘₯(1) = π‘₯(1) = 0.5, π‘₯(0) = 0.4343, and π‘₯(0) = 0.5307. Therefore we have analytical solution and EP solution by (24) and (29). By (10) and (14), we find πœ† = 0.7067, πœ‡ = 0.3296, and πœ’ (π‘Ÿ) =

βˆ’π‘Ÿ2

2 βˆ’ 2π‘Ÿ . (38) + 2π‘Ÿ + 15 + πœ‡ (9 βˆ’ π‘Ÿ2 ) βˆ’ πœ† (π‘Ÿ2 βˆ’ 4π‘Ÿ βˆ’ 5)

Notice that the numerical methods needed π‘₯(0) and π‘₯(0) to obtain initial guess and often these values achieve analytical solution, but in Example 3 these values achieve EP method because in this example we do not have analytical solution.

5. Conclusion In this paper we introduced a new method to solve a fully fuzzy quadratic equation. To this purpose we found the optimum spreads to decrease maximum error. One of the advantages of this method is that complications do not depend on the sign of the coefficients and variable. It is possible that these equations do not have any analytical solution, but the proposed method gives us an approximate analytical solution.

Μƒ = (0.5βˆ’πœ†πœ’(π‘Ÿ), 0.5+πœ‡πœ’(π‘Ÿ)), Hausdorff In this example, with 𝑋 Μƒ 𝐷) Μƒ = 0.0228. metric is 𝐷(𝐹(𝑋),

Conflict of Interests

Μƒ = (4/2/2), 𝐡̃ = (2/2/2), and 𝐷 Μƒ = Example 2. Letting 𝐴 (1/0.5/0.5), we have

The authors declare that there is no conflict of interests regarding the publication of this paper.

(2 + 2π‘Ÿ, 6 βˆ’ 2π‘Ÿ) (π‘₯2 (π‘Ÿ) , π‘₯2 (π‘Ÿ)) + (2π‘Ÿ, 4 βˆ’ 2π‘Ÿ) (π‘₯ (π‘Ÿ) , π‘₯ (π‘Ÿ)) = (0.5 + 0.5π‘Ÿ, 1.5 βˆ’ 0.5π‘Ÿ) . (39) The real roots of 1-cut and 0-cut equations are 𝛼 = π‘₯(1) = π‘₯(1) = (βˆ’1 + √5)/4 = 0.3090, π‘₯(0) = 0.5, and π‘₯(0) = 0.2676. Therefore this example does not have analytical solution. We look for EP solution. By (26) we find that this example does not have EP solution with πœ† = πœ‡ = 0.5 too. Μƒ < 0, 𝐡̃ > 0, and Now we consider an example with 𝐴 Μƒ > 0. 𝑋 Μƒ = (βˆ’2/1/1), 𝐡̃ = (3/1/2), and 𝐷 Μƒ = Example 3. Letting 𝐴 (1/1/2), we have (βˆ’3 + π‘Ÿ, βˆ’1 βˆ’ π‘Ÿ) (π‘₯2 (π‘Ÿ) , π‘₯2 (π‘Ÿ)) + (2 + π‘Ÿ, 5 βˆ’ 2π‘Ÿ) (π‘₯ (π‘Ÿ) , π‘₯ (π‘Ÿ)) = (π‘Ÿ, 3 βˆ’ 2π‘Ÿ) . (40) The real roots of 1-cut and 0-cut equations are 𝛼 = π‘₯(1) = π‘₯(1) = 0.5, 1, π‘₯(0) = 0.7838, 1.2527, and π‘₯(0) = 0.7229, 0.914. Therefore, this example does not have analytical solution. We look for EP solution. By (25), (30), and (15), we find that this example has EP solution with πœ† = πœ‡ = 0.5 as follows: Μƒ = (0.5 βˆ’ 0.5πœ’(π‘Ÿ), 0.5 + 0.5πœ’(π‘Ÿ)) that 𝑋 πœ’ (π‘Ÿ) =

βˆ’π‘Ÿ2 + 6π‘Ÿ βˆ’ 5 , βˆ’π‘Ÿ2 + 4π‘Ÿ βˆ’ 23

(41)

and 0-cut of EP solution is π‘₯(0) = 0.3913 and π‘₯(0) = 0.6087. In this example by using Hausdorff metric we have Μƒ 𝐷) Μƒ = 0.3290. 𝐷(𝐹(𝑋),

References [1] J. J. Buckley, β€œThe fuzzy mathematics of finance,” Fuzzy Sets and Systems, vol. 21, no. 3, pp. 257–273, 1987. [2] M. Li Calzi, β€œTowards a general setting for the fuzzy mathematics of finance,” Fuzzy Sets and Systems, vol. 35, no. 3, pp. 265–280, 1990. [3] J. J. Buckley, β€œSolving fuzzy equations in economics and finance,” Fuzzy Sets and Systems, vol. 48, no. 3, pp. 289–296, 1992. [4] J. J. Buckley and Y. Qu, β€œSolving systems of linear fuzzy equations,” Fuzzy Sets and Systems, vol. 43, no. 1, pp. 33–43, 1991. [5] Q.-X. Li and S.-F. Liu, β€œThe foundation of the grey matrix and the grey input-output analysis,” Applied Mathematical Modelling, vol. 32, no. 3, pp. 267–291, 2008. [6] C. C. Wu and N. B. Chang, β€œGrey input-output analysis and its application for environmental cost allocation,” European Journal of Operational Research, vol. 145, no. 1, pp. 175–201, 2003. [7] S.-H. Chen and X.-W. Yang, β€œInterval finite element method for beam structures,” Finite Elements in Analysis and Design, vol. 34, no. 1, pp. 75–88, 2000. [8] P. Sevastjanov and L. Dymova, β€œA new method for solving interval and fuzzy equations: linear case,” Information Sciences, vol. 179, no. 7, pp. 925–937, 2009. [9] J. J. Buckley and Y. Qu, β€œSolving linear and quadratic fuzzy equations,” Fuzzy Sets and Systems, vol. 38, no. 1, pp. 43–59, 1990. [10] J. J. Buckley and Y. Qu, β€œOn using 𝛼-cuts to evaluate fuzzy equations,” Fuzzy Sets and Systems, vol. 38, no. 3, pp. 309–312, 1990. [11] J. J. Buckley and Y. Qu, β€œSolving fuzzy equations: a new solution concept,” Fuzzy Sets and Systems, vol. 39, no. 3, pp. 291–301, 1991. [12] S. Abbasbandy and M. Otadi, β€œNumerical solution of fuzzy polynomials by fuzzy neural network,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1084–1089, 2006. [13] S. Abbasbandy and B. Asady, β€œNewton's method for solving fuzzy nonlinear equations,” Applied Mathematics and Computation, vol. 159, no. 2, pp. 349–356, 2004.

6 [14] T. Allahviranloo, M. Otadi, and M. Mosleh, β€œIterative method for fuzzy equations,” Soft Computing, vol. 12, no. 10, pp. 935–939, 2008. [15] J. J. Buckley and E. Eslami, β€œNeural net solutions to fuzzy problems: the quadratic equation,” Fuzzy Sets and Systems, vol. 86, no. 3, pp. 289–298, 1997. [16] J. J. Buckley, E. Eslami, and Y. Hayashi, β€œSolving fuzzy equations using neural nets,” Fuzzy Sets and Systems, vol. 86, no. 3, pp. 271– 278, 1997. [17] D. Dubois and H. Prade, β€œOperations on fuzzy numbers,” Journal of Systems Sience, vol. 9, pp. 613–626, 1978. [18] D. Dubois and H. Prade, Theory and Application, Fuzzy Sets and Systems, Academic Press, 1980. [19] L. A. Zadeh, β€œFuzzy sets,” Information and Control, vol. 8, no. 3, pp. 338–353, 1965. [20] H. J. Zimmermann, Fuzzy Set Theory and Its Applications, Kluwer Academic Press, Dordrecht, The Netherlands, 1991. [21] W. Cong-Xin and M. Ming, β€œEmbedding problem of fuzzy number space: part I,” Fuzzy Sets and Systems, vol. 44, no. 1, pp. 33–38, 1991. [22] R. Goetschel Jr. and W. Voxman, β€œElementary fuzzy calculus,” Fuzzy Sets and Systems, vol. 18, no. 1, pp. 31–43, 1986. [23] L. Jaulin, M. Kieffir, O. Didrit, and E. Walter, Applied Interval Analysis, Springer, London, UK, 2001. [24] R. E. Moore, Interval Analysis, Prentice-Hall, Englewood Cliffs, NJ, USA, 1966. [25] P. Diamond and P. Kloeden, Metric Space of Fuzzy Sets, WordScientific, Singapore, 1994.

The Scientific World Journal

Advances in

Operations Research Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Advances in

Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Applied Mathematics

Algebra

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Probability and Statistics Volume 2014

The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Differential Equations Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Submit your manuscripts at http://www.hindawi.com International Journal of

Advances in

Combinatorics Hindawi Publishing Corporation http://www.hindawi.com

Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of Mathematics and Mathematical Sciences

Mathematical Problems in Engineering

Journal of

Mathematics Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Discrete Mathematics

Journal of

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Discrete Dynamics in Nature and Society

Journal of

Function Spaces Hindawi Publishing Corporation http://www.hindawi.com

Abstract and Applied Analysis

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Journal of

Stochastic Analysis

Optimization

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Suggest Documents