Elsevier

Fuzzy Sets and Systems

Volume 257, 16 December 2014, Pages 169-183
Fuzzy Sets and Systems

Solution of linear differential equations with fuzzy boundary values

https://doi.org/10.1016/j.fss.2013.08.008Get rights and content

Abstract

We investigate linear differential equations with boundary values expressed by fuzzy numbers. In contrast to most approaches, which search for a fuzzy-valued function as the solution, we search for a fuzzy set of real functions as the solution. We define a real function as an element of the solution set if it satisfies the differential equation and its boundary values are in intervals determined by the corresponding fuzzy numbers. The membership degree of the real function is defined as the lowest value among membership degrees of its boundary values in the corresponding fuzzy sets. To find the fuzzy solution, we use a method based on the properties of linear transformations. We show that the fuzzy problem has a unique solution if the corresponding crisp problem has a unique solution. We prove that if the boundary values are triangular fuzzy numbers, then the value of the solution at a given time is also a triangular fuzzy number. The defined solution is the same as one of the solutions obtained by Zadeh's extension principle. For a second-order differential equation with constant coefficients, the solution is expressed in analytical form. Examples are given to describe the proposed approach and to compare it to a method that uses the generalized Hukuhara derivative, which demonstrates the advantages of our method.

Introduction

Approaches to fuzzy boundary value problems (FBVPs) and other fuzzy problems can be of three types. The first approach assumes that even if only the boundary values are fuzzy, the solution is a fuzzy function and consequently the derivatives in the differential equation must be considered as fuzzy derivatives. These can be either Hukuhara or generalized derivatives. Bede demonstrated that a large class of BVPs have no solution if the Hukuhara derivative is used [4]. To overcome this difficulty, the concept of a generalized derivative was developed [5], [12] and fuzzy differential equations were investigated using this concept [7], [11], [13], [29]. Khastan and Nieto found solutions for a large enough class of boundary value problems using the generalized derivative [20]. However, their examples and examples presented here show that these solutions are difficult to interpret, because the four different problems obtained using the generalized first and second derivatives often do not reflect the nature of the problem. Liu reduced these four problems to two when the right-hand-side function is monotonic [23].

Khastan et al. demonstrated that a class of first-order linear differential equations subject to periodic boundary conditions can be solved by switching between two types of generalized derivative over the time domain [21]. Rodríguez-López improved this result for equations whose coefficients may change their sign finitely many times [28]. Using novel generalizations of the Hukuhara difference for fuzzy sets, Bede and Stefanini investigated new generalized differentiability concepts for fuzzy-valued functions [8].

To avoid difficulties with fuzzy derivatives, some researchers propose that an equivalent fuzzy integral equation should be solved instead of a fuzzy differential equation [1], [14], [27].

The second approach is based on Zadeh's extension principle. For a fuzzy initial value problem, the associated crisp problem is solved and in the solution the initial fuzzy value is substituted instead of the real constant. In the final solution, arithmetic operations are considered to be operations on fuzzy numbers [9], [10]. To find an exact solution to the periodic boundary value problem for a first-order linear fuzzy differential equation with impulses, Nieto et al. used the crisp solution [26].

In the third approach, the fuzzy problem is transformed to a crisp problem. There are two ways to realize this approach. The first, suggested by Hüllermeier [18], uses the concept of differential inclusion. In this way, by taking an α-cut of the initial value and the solution, the given differential equation is converted to a differential inclusion and the solution is accepted as the α-cut of the fuzzy solution. Li et al. used this approach in a recent study [22]. In the paper the concept of big solutions is introduced and some existence and uniqueness theorems are obtained. Misukoshi et al. [24] have proved that, under certain conditions, this approach is equivalent to the second one for the initial value problem. The second way is offered by Gasilov et al. [15]. In this way the fuzzy problem is considered to be a set of crisp problems.

In this study, we investigate a differential equation with fuzzy boundary values. We interpret the problem as a set of crisp problems. For linear equations, we propose a method based on the properties of linear transformations. We show that, if the solution of the corresponding crisp problem exists and is unique, the fuzzy problem also has unique solution. Moreover, we prove that if the boundary values are triangular fuzzy numbers, then the value of the solution is a triangular fuzzy number at each time. We explain the proposed method on examples. We find analytical expression for solution of second-order linear differential equation with constant coefficients. We demonstrate with an example the advantages of the proposed method over the method which uses the generalized Hukuhara derivative.

The remainder of the paper is organized as follows. Section 2 provides preliminary information. In Section 3 we describe the fuzzy boundary value problem and the solution concept. The methodology is proposed in Section 4. We apply this methodology to solve examples in Section 5. In Section 6 we compare our approach with a method that uses a generalized Hukuhara derivative. Section 7 concludes and identifies further research directions.

Section snippets

Basic concepts of fuzzy set theory

The notion of a fuzzy set is an extension of the classical notion of a set. In classical set theory, an element either belongs or does not belong to a given set. By contrast, in fuzzy set theory, an element has a degree of membership, which is a real number from [0,1], in a given fuzzy set. In fuzzy set theory, classical sets are usually called crisp sets.

A fuzzy set A˜ can be defined as a pair of the universal set U and the membership function μ:U[0,1]. If the universal set U is fixed, a

Fuzzy boundary value problem

In some applications, the behavior of an object is determined by physics laws and is crisp. However, if the boundary values are obtained from measurements, for example, these values can be uncertain and often it is more suitable to model them using fuzzy numbers. This gives rise to BVPs with crisp dynamics but fuzzy boundary values. In this paper, we consider a BVP involving a crisp differential equation but with fuzzy boundary values:{x+a1(t)x+a2(t)x=f(t)x(0)=A˜x(T)=B˜.

Let U (the universal

Solution methodology

Lemma 1

X˜(t), the value of the solution to (8) at a given time t, is well defined, that is, it is a fuzzy number.

Proof

According to (9), an α-cut of X˜ is expressed asXα={x()|x+a1(t)x+a2(t)x=0;x(0)=a;x(T)=b;a[aα̲,aα¯];b[bα̲,bα¯]}. Let linear independent solutions of the crisp differential equation (8), x1(t) and x2(t), be known. Then we can constitute the vector-function w(t) using (2). According to (4) and (10), we haveXα={x()|x=wu;u=(a,b);a[aα̲,aα¯];b[bα̲,bα¯]}.

We consider a fixed time t and let v

Examples

In this section, we solve two examples to demonstrate the proposed method.

Example 1

Solve the FBVP:{x4x+4x=12t2x(0)=(2,3,4)x(1)=(1,2,2.5).

Solution

We represent the solution as x˜(t)=xcr(t)+x˜un(t).

  • 1.

    We solve the crisp non-homogeneous problem{x4x+4x=12t2x(0)=3x(1)=2 and find the crisp solution xcr(t)=12(t+1)2+3.5(1t)e2t+4te2(t1) (dashed line in Fig. 1).

  • 2.

    To determine the uncertainty, we consider the fuzzy homogeneous problem:{x4x+4x=0x(0)=(1,0,1)x(1)=(1,0,0.5).

    x1(t)=e2t and x2(t)=te2t are linear

Comparison

In this section, we compare the solution obtained by our method with the solution obtained by a method that uses the generalized Hukuhara (gH) derivative. For this purpose we use a numerical example.

Example 3

Solve the FBVP:{x=16xx(0)=(1,0,0.5)x(2)=(0.5,0,0.5).

Solution

We first solve the problem using our method. We note that the problem is homogeneous and the boundary values are fuzzy numbers with vertices at 0. Therefore, the solution can be calculated according to (14).

x1(t)=cos4t and x2(t)=sin4t are

Conclusion

We searched for an FBVP solution as a fuzzy set of real functions. Examples demonstrated that this approach is more effective in some cases than other approaches. We explained the proposed method for a second-order linear differential equation. However, the approach can be generalized for nth-order m-point linear FBVPs.

In future we plan to apply the approach to linear BVPs with a fuzzy forcing function and to non-linear equations.

Acknowledgements

We are grateful to the reviewers for comments and suggestions that helped to improve the quality of this paper. We also thank Associate Professor Dr. Süleyman Tosun (Ankara University, Turkey) for helpful comments.

References (30)

Cited by (48)

  • Similarity solutions of fractional parabolic boundary value problems with uncertainty

    2021, Communications in Nonlinear Science and Numerical Simulation
  • On interactive fuzzy boundary value problems

    2019, Fuzzy Sets and Systems
    Citation Excerpt :

    Their approach may not produce fuzzy solution on the entire domain, that is, the solution obtained at a given instant t may not be a fuzzy number. Moreover, Gasilov et al. argued that the solutions obtained using Khastan and Nieto's method are difficult to interpret because the solutions of the four different problems may not reflect the nature of the studied phenomenon [4]. Other approaches to FBVPs are based on the theory of differential inclusion [5] and modification/adaptation of classical numerical methods for the fuzzy case, such as the undetermined coefficients method [6], the finite difference method [7], and the finite element method [8].

View all citing articles on Scopus
View full text